Science.gov

Sample records for 2d compositional modeling

  1. MODELING THE TRANSVERSE THERMAL CONDUCTIVITY OF 2-D SICF/SIC COMPOSITES MADE WITH WOVEN FABRIC

    SciTech Connect

    Youngblood, Gerald E; Senor, David J; Jones, Russell H

    2004-06-01

    The hierarchical two-layer (H2L) model describes the effective transverse thermal conductivity (Keff) of a 2D-SiCf/SiC composite plate made from stacked and infiltrated woven fabric layers in terms of constituent properties and microstructural and architectural variables. The H2L model includes the effects of fiber-matrix interfacial conductance, high fiber packing fractions within individual tows and the non-uniform nature of 2D fabric/matrix layers that usually include a significant amount of interlayer porosity. Previously, H2L model Keff-predictions were compared to measured values for two versions of 2D Hi-Nicalon/PyC/ICVI-SiC composite, one with a “thin” (0.11m) and the other with a “thick” (1.04m) pyrocarbon (PyC) fiber coating, and for a 2D Tyranno SA/”thin” PyC/FCVI-SIC composite. In this study, H2L model Keff-predictions were compared to measured values for a 2D-SiCf/SiC composite made using the ICVI-process with Hi-Nicalon type S fabric and a “thin” PyC fiber coating. The values of Keff determined for the latter composite were significantly greater than the Keff-values determined for the composites made with either the Hi-Nicalon or the Tyranno SA fabrics. Differences in Keff-values were expected for the different fiber types, but major differences also were due to observed microstructural and architectural variations between the composite systems, and as predicted by the H2L model.

  2. Modeling the Elastic Modulus of 2D Woven CVI SiC Composites

    NASA Technical Reports Server (NTRS)

    Morscher, Gregory N.

    2006-01-01

    The use of fiber, interphase, CVI SiC minicomposites as structural elements for 2D-woven SiC fiber reinforced chemically vapor infiltrated (CVI) SiC matrix composites is demonstrated to be a viable approach to model the elastic modulus of these composite systems when tensile loaded in an orthogonal direction. The 0deg (loading direction) and 90deg (perpendicular to loading direction) oriented minicomposites as well as the open porosity and excess SiC associated with CVI SiC composites were all modeled as parallel elements using simple Rule of Mixtures techniques. Excellent agreement for a variety of 2D woven Hi-Nicalon(TradeMark) fiber-reinforced and Sylramic-iBN reinforced CVI SiC matrix composites that differed in numbers of plies, constituent content, thickness, density, and number of woven tows in either direction (i.e, balanced weaves versus unbalanced weaves) was achieved. It was found that elastic modulus was not only dependent on constituent content, but also the degree to which 90deg minicomposites carried load. This depended on the degree of interaction between 90deg and 0deg minicomposites which was quantified to some extent by composite density. The relationships developed here for elastic modulus only necessitated the knowledge of the fractional contents of fiber, interphase and CVI SiC as well as the tow size and shape. It was concluded that such relationships are fairly robust for orthogonally loaded 2D woven CVI SiC composite system and can be implemented by ceramic matrix composite component modelers and designers for modeling the local stiffness in simple or complex parts fabricated with variable constituent contents.

  3. 2D and 3D numerical models on compositionally buoyant diapirs in the mantle wedge

    NASA Astrophysics Data System (ADS)

    Hasenclever, Jörg; Morgan, Jason Phipps; Hort, Matthias; Rüpke, Lars H.

    2011-11-01

    We present 2D and 3D numerical model calculations that focus on the physics of compositionally buoyant diapirs rising within a mantle wedge corner flow. Compositional buoyancy is assumed to arise from slab dehydration during which water-rich volatiles enter the mantle wedge and form a wet, less dense boundary layer on top of the slab. Slab dehydration is prescribed to occur in the 80-180 km deep slab interval, and the water transport is treated as a diffusion-like process. In this study, the mantle's rheology is modeled as being isoviscous for the benefit of easier-to-interpret feedbacks between water migration and buoyant viscous flow of the mantle. We use a simple subduction geometry that does not change during the numerical calculation. In a large set of 2D calculations we have identified that five different flow regimes can form, in which the position, number, and formation time of the diapirs vary as a function of four parameters: subduction angle, subduction rate, water diffusivity (mobility), and mantle viscosity. Using the same numerical method and numerical resolution we also conducted a suite of 3D calculations for 16 selected parameter combinations. Comparing the 2D and 3D results for the same model parameters reveals that the 2D models can only give limited insights into the inherently 3D problem of mantle wedge diapirism. While often correctly predicting the position and onset time of the first diapir(s), the 2D models fail to capture the dynamics of diapir ascent as well as the formation of secondary diapirs that result from boundary layer perturbations caused by previous diapirs. Of greatest importance for physically correct results is the numerical resolution in the region where diapirs nucleate, which must be high enough to accurately capture the growth of the thin wet boundary layer on top of the slab and, subsequently, the formation, morphology, and ascent of diapirs. Here 2D models can be very useful to quantify the required resolution, which we

  4. Modeling the Transverse Thermal Conductivity of 2-D SiCf/SiC Composites Made with Woven Fabric

    SciTech Connect

    Youngblood, Gerald E.; Senor, David J.; Jones, Russell H.

    2004-06-30

    The hierarchical two-layer (H2L) model was developed to describe the effective transverse thermal conductivity, Keff, of a 2D-SiCf/SiC composite made from stacked and infiltrated woven fabric layers in terms of constituent properties and microstructural and architectural variables. The H2L model includes the expected effects of fiber-matrix interfacial conductance as well as the effects of high fiber packing fractions within individual tows and the non-uniform nature of 2D-fabric layers that usually include a significant amount of interlayer porosity. Previously, H2L model predictions were compared to measured values of Keff for two versions of DuPont 2D-Hi NicalonÔ/PyC/ICVI-SiC composite, one with a “thin” (0.110 μm) and the other with a “thick” (1.040 μm) pyrocarbon (PyC) fiber coating, and for a 2D-TyrannoÔ SA/”thin” PyC/FCVI-SIC composite made by ORNL. In this study, H2L model predictions are compared to measured Keff-values for a 2D-SiCf/SiC composite made by GE Power Systems (formerly DuPont Lanxide) using the ICVI-process with Hi-NicalonÔ type S fabric. The values of Keff determined for the composite made with the Hi-NicalonÔ type S fabric were significantly greater than Keff-values determined for the composites made with either the Hi-NicalonÔor the TyrannoÔ SA fabrics. Differences in Keff-values were expected for using different fiber types, but major differences also were due to observed microstructural variations between the systems, and as predicted by the H2L model.

  5. Investigations on void morphology in CFRP composite materials and ultrasonic scattering attenuation based on a 2D random void model

    NASA Astrophysics Data System (ADS)

    Lin, L.; Ding, S. S.; Chen, J.; Liang, X. Y.; Li, X. M.

    2012-05-01

    A 2D random void model (RVM) is proposed to describe voids morphology in Carbon Fiber Reinforced Plastic (CFRP) composite materials and used to investigate Ultrasonic Scattering Attenuation Coefficient (USAC). Void morphology simulations from RVM present good matches to micrographic observations. The fluctuations of USAC due to the randomness of void morphology and their dependence on the frequency have been discussed, which are significantly helpful to clarify ultrasonic scattering attenuation mechanism from voids in nature.

  6. Multi-scale Model of Residual Strength of 2D Plain Weave C/SiC Composites in Oxidation Atmosphere

    NASA Astrophysics Data System (ADS)

    Chen, Xihui; Sun, Zhigang; Sun, Jianfen; Song, Yingdong

    2016-06-01

    Multi-scale models play an important role in capturing the nonlinear response of woven carbon fiber reinforced ceramic matrix composites. In plain weave carbon fiber/silicon carbon (C/SiC) composites, the carbon fibers and interphases will be oxidized at elevated temperature and the strength of the composite will be degraded when oxygen enters micro-cracks formed in the as-produced parts due to the mismatch in thermal properties between constituents. As a result of the oxidation on fiber surface, fiber shows a notch-like morphology. In this paper, the change rule of fiber notch depth is fitted by circular function. And a multi-scale model based upon the change rule of fiber notch depth is developed to simulate the residual strength and post-oxidation stress-strain curves of the composite. The multi-scale model is able to accurately predict the residual strength and post-oxidation stress-strain curves of the composite. Besides, the simulated residual strength and post-oxidation stress-strain curves of 2D plain weave C/SiC composites in oxidation atmosphere show good agreements with experimental results. Furthermore, the oxidation time and temperature of the composite are investigated to show their influences upon the residual strength and post-oxidation stress-strain curves of plain weave C/SiC composites.

  7. A hierarchical lattice spring model to simulate the mechanics of 2-D materials-based composites

    NASA Astrophysics Data System (ADS)

    Brely, Lucas; Bosia, Federico; Pugno, Nicola

    2015-07-01

    In the field of engineering materials, strength and toughness are typically two mutually exclusive properties. Structural biological materials such as bone, tendon or dentin have resolved this conflict and show unprecedented damage tolerance, toughness and strength levels. The common feature of these materials is their hierarchical heterogeneous structure, which contributes to increased energy dissipation before failure occurring at different scale levels. These structural properties are the key to exceptional bioinspired material mechanical properties, in particular for nanocomposites. Here, we develop a numerical model in order to simulate the mechanisms involved in damage progression and energy dissipation at different size scales in nano- and macro-composites, which depend both on the heterogeneity of the material and on the type of hierarchical structure. Both these aspects have been incorporated into a 2-dimensional model based on a Lattice Spring Model, accounting for geometrical nonlinearities and including statistically-based fracture phenomena. The model has been validated by comparing numerical results to continuum and fracture mechanics results as well as finite elements simulations, and then employed to study how structural aspects impact on hierarchical composite material properties. Results obtained with the numerical code highlight the dependence of stress distributions on matrix properties and reinforcement dispersion, geometry and properties, and how failure of sacrificial elements is directly involved in the damage tolerance of the material. Thanks to the rapidly developing field of nanocomposite manufacture, it is already possible to artificially create materials with multi-scale hierarchical reinforcements. The developed code could be a valuable support in the design and optimization of these advanced materials, drawing inspiration and going beyond biological materials with exceptional mechanical properties.

  8. 2D photochemical modeling of Saturn's stratosphere. Part II: Feedback between composition and temperature

    NASA Astrophysics Data System (ADS)

    Hue, V.; Greathouse, T. K.; Cavalié, T.; Dobrijevic, M.; Hersant, F.

    2016-03-01

    Saturn's axial tilt of 26.7° produces seasons in a similar way as on Earth. Both the stratospheric temperature and composition are affected by this latitudinally varying insolation along Saturn's orbital path. The atmospheric thermal structure is controlled and regulated by the amount of hydrocarbons in the stratosphere, which act as absorbers and coolants from the UV to the far-IR spectral range, and this structure has an influence on the amount of hydrocarbons. We study here the feedback between the chemical composition and the thermal structure by coupling a latitudinal and seasonal photochemical model with a radiative seasonal model. Our results show that the seasonal temperature peak in the higher stratosphere, associated with the seasonal increase of insolation, is shifted earlier than the maximum insolation peak. This shift is increased with increasing latitudes and is caused by the low amount of stratospheric coolants in the spring season. At 80° in both hemispheres, the temperature peak at 10-2 mbar is seen to occur half a season (3-4 Earth years) earlier than was previously predicted by radiative seasonal models that assumed spatially and temporally uniform distribution of coolants. This shift progressively decreases with increasing pressure, up to around the 0.5 mbar pressure level where it vanishes. On the opposite, the thermal field has a small feedback on the abundance distributions. Accounting for that feedback modifies the predicted equator-to-pole temperature gradient. The meridional gradients of temperature at the mbar pressure levels are better reproduced when this feedback is accounted for. At lower pressure levels, Saturn's stratospheric thermal structure seems to depart from pure radiative seasonal equilibrium as previously suggested by Guerlet et al. (2014). Although the agreement with the absolute value of the stratospheric temperature observed by Cassini is moderate, it is a mandatory step toward a fully coupled GCM-photochemical model.

  9. Numerical Modeling of Oxidized 2D C/SiC Composites in Air Environments Below 900 °C: Microstructure and Elastic Properties

    NASA Astrophysics Data System (ADS)

    Sun, Zhigang; Chen, Xihui; Shao, Hongyan; Song, Yingdong

    2016-08-01

    A numerical model is presented for simulation of the oxidation-affected behaviors of two dimensional carbon fiber-reinforced silcon carbide matrix composite (2D C/SiC) exposed to air oxidizing environments below 900 °C, which incorporates the modeling of oxidized microstructure and computing of degraded elastic properties. This model is based upon the analysis of the representative volume cell (RVC) of the composite. The multi-scale model of 2D C/SiC composites is concerned in the present study. Analysis results of such a composite can provide a guideline for the real 2D C/SiC composite. The micro-structure during oxidation process is firstly modeled in the RVC. The elastic moduli of oxidized composite under non-stress oxidation environment is computed by finite element analysis. The elastic properties of 2D-C/SiC composites in air oxidizing environment are evaluated and validated in comparison to experimental data. The oxidation time, temperature and fiber volume fractions of C/SiC composite are investigated to show their influences upon the elastic properties of 2D C/SiC composites.

  10. Numerical Modeling of Oxidized 2D C/SiC Composites in Air Environments Below 900 °C: Microstructure and Elastic Properties

    NASA Astrophysics Data System (ADS)

    Sun, Zhigang; Chen, Xihui; Shao, Hongyan; Song, Yingdong

    2016-04-01

    A numerical model is presented for simulation of the oxidation-affected behaviors of two dimensional carbon fiber-reinforced silcon carbide matrix composite (2D C/SiC) exposed to air oxidizing environments below 900 °C, which incorporates the modeling of oxidized microstructure and computing of degraded elastic properties. This model is based upon the analysis of the representative volume cell (RVC) of the composite. The multi-scale model of 2D C/SiC composites is concerned in the present study. Analysis results of such a composite can provide a guideline for the real 2D C/SiC composite. The micro-structure during oxidation process is firstly modeled in the RVC. The elastic moduli of oxidized composite under non-stress oxidation environment is computed by finite element analysis. The elastic properties of 2D-C/SiC composites in air oxidizing environment are evaluated and validated in comparison to experimental data. The oxidation time, temperature and fiber volume fractions of C/SiC composite are investigated to show their influences upon the elastic properties of 2D C/SiC composites.

  11. Modeling the Effect of Multiple Matrix Cracking Modes on Cyclic Hysteresis Loops of 2D Woven Ceramic-Matrix Composites

    NASA Astrophysics Data System (ADS)

    Longbiao, Li

    2016-08-01

    In this paper, the effect of multiple matrix cracking modes on cyclic loading/unloading hysteresis loops of 2D woven ceramic-matrix composites (CMCs) has been investigated. The interface slip between fibers and the matrix existed in matrix cracking mode 3 and mode 5, in which matrix cracking and interface debonding occurred in longitudinal yarns, are considered as the major reason for hysteresis loops of 2D woven CMCs. The effects of fiber volume content, peak stress, matrix crack spacing, interface properties, matrix cracking mode proportion and interface wear on interface slip and hysteresis loops have been analyzed. The cyclic loading/unloading hysteresis loops of 2D woven SiC/SiC composite corresponding to different peak stresses have been predicted using the present analysis. It was found that the damage parameter, i.e., the proportion of matrix cracking mode 3 in the entire cracking modes of the composite, increases with increasing peak stress.

  12. Modeling the Effect of Multiple Matrix Cracking Modes on Cyclic Hysteresis Loops of 2D Woven Ceramic-Matrix Composites

    NASA Astrophysics Data System (ADS)

    Longbiao, Li

    2016-02-01

    In this paper, the effect of multiple matrix cracking modes on cyclic loading/unloading hysteresis loops of 2D woven ceramic-matrix composites (CMCs) has been investigated. The interface slip between fibers and the matrix existed in matrix cracking mode 3 and mode 5, in which matrix cracking and interface debonding occurred in longitudinal yarns, are considered as the major reason for hysteresis loops of 2D woven CMCs. The effects of fiber volume content, peak stress, matrix crack spacing, interface properties, matrix cracking mode proportion and interface wear on interface slip and hysteresis loops have been analyzed. The cyclic loading/unloading hysteresis loops of 2D woven SiC/SiC composite corresponding to different peak stresses have been predicted using the present analysis. It was found that the damage parameter, i.e., the proportion of matrix cracking mode 3 in the entire cracking modes of the composite, increases with increasing peak stress.

  13. 2D photochemical modeling of Saturn's stratosphere. Part I: Seasonal variation of atmospheric composition without meridional transport

    NASA Astrophysics Data System (ADS)

    Hue, V.; Cavalié, T.; Dobrijevic, M.; Hersant, F.; Greathouse, T. K.

    2015-09-01

    Saturn's axial tilt of 26.7° produces seasons in a similar way as on Earth. Both the stratospheric temperature and composition are affected by this latitudinally varying insolation along Saturn's orbital path. A new time-dependent 2D photochemical model is presented to study the seasonal evolution of Saturn's stratospheric composition. This study focuses on the impact of the seasonally variable thermal field on the main stratospheric C2-hydrocarbon chemistry (C2H2 and C2H6) using a realistic radiative climate model. Meridional mixing and advective processes are implemented in the model but turned off in the present study for the sake of simplicity. The results are compared to a simple study case where a latitudinally and temporally steady thermal field is assumed. Our simulations suggest that, when the seasonally variable thermal field is accounted for, the downward diffusion of the seasonally produced hydrocarbons is faster due to the seasonal compression of the atmospheric column during winter. This effect increases with increasing latitudes which experience the most important thermal changes in the course of the seasons. The seasonal variability of C2H2 and C2H6 therefore persists at higher-pressure levels with a seasonally-variable thermal field. Cassini limb-observations of C2H2 and C2H6 (Guerlet, S. et al. [2009]. Icarus 203, 214-232) are reasonably well-reproduced from the equator to 40° in both hemispheres in the 0.1-1 mbar pressure range. At lower pressure levels, the models only fit the Cassini observations in the northern hemisphere, from the equator to 40°N. Beyond 40° in both hemispheres, deviations from the pure photochemical predictions, mostly in the southern hemisphere, suggest the presence of large-scale stratospheric dynamics.

  14. A model for calculating the errors of 2D bulk analysis relative to the true 3D bulk composition of an object, with application to chondrules

    NASA Astrophysics Data System (ADS)

    Hezel, Dominik C.

    2007-09-01

    Certain problems in Geosciences require knowledge of the chemical bulk composition of objects, such as, for example, minerals or lithic clasts. This 3D bulk chemical composition (bcc) is often difficult to obtain, but if the object is prepared as a thin or thick polished section a 2D bcc can be easily determined using, for example, an electron microprobe. The 2D bcc contains an error relative to the true 3D bcc that is unknown. Here I present a computer program that calculates this error, which is represented as the standard deviation of the 2D bcc relative to the real 3D bcc. A requirement for such calculations is an approximate structure of the 3D object. In petrological applications, the known fabrics of rocks facilitate modeling. The size of the standard deviation depends on (1) the modal abundance of the phases, (2) the element concentration differences between phases and (3) the distribution of the phases, i.e. the homogeneity/heterogeneity of the object considered. A newly introduced parameter " τ" is used as a measure of this homogeneity/heterogeneity. Accessory phases, which do not necessarily appear in 2D thin sections, are a second source of error, in particular if they contain high concentrations of specific elements. An abundance of only 1 vol% of an accessory phase may raise the 3D bcc of an element by up to a factor of ˜8. The code can be queried as to whether broad beam, point, line or area analysis technique is best for obtaining 2D bcc. No general conclusion can be deduced, as the error rates of these techniques depend on the specific structure of the object considered. As an example chondrules—rapidly solidified melt droplets of chondritic meteorites—are used. It is demonstrated that 2D bcc may be used to reveal trends in the chemistry of 3D objects.

  15. Brittle damage models in DYNA2D

    SciTech Connect

    Faux, D.R.

    1997-09-01

    DYNA2D is an explicit Lagrangian finite element code used to model dynamic events where stress wave interactions influence the overall response of the system. DYNA2D is often used to model penetration problems involving ductile-to-ductile impacts; however, with the advent of the use of ceramics in the armor-anti-armor community and the need to model damage to laser optics components, good brittle damage models are now needed in DYNA2D. This report will detail the implementation of four brittle damage models in DYNA2D, three scalar damage models and one tensor damage model. These new brittle damage models are then used to predict experimental results from three distinctly different glass damage problems.

  16. Matrix models of 2d gravity

    SciTech Connect

    Ginsparg, P.

    1991-01-01

    These are introductory lectures for a general audience that give an overview of the subject of matrix models and their application to random surfaces, 2d gravity, and string theory. They are intentionally 1.5 years out of date.

  17. Matrix models of 2d gravity

    SciTech Connect

    Ginsparg, P.

    1991-12-31

    These are introductory lectures for a general audience that give an overview of the subject of matrix models and their application to random surfaces, 2d gravity, and string theory. They are intentionally 1.5 years out of date.

  18. Comparison of 2D Finite Element Modeling Assumptions with Results From 3D Analysis for Composite Skin-Stiffener Debonding

    NASA Technical Reports Server (NTRS)

    Krueger, Ronald; Paris, Isbelle L.; OBrien, T. Kevin; Minguet, Pierre J.

    2004-01-01

    The influence of two-dimensional finite element modeling assumptions on the debonding prediction for skin-stiffener specimens was investigated. Geometrically nonlinear finite element analyses using two-dimensional plane-stress and plane-strain elements as well as three different generalized plane strain type approaches were performed. The computed skin and flange strains, transverse tensile stresses and energy release rates were compared to results obtained from three-dimensional simulations. The study showed that for strains and energy release rate computations the generalized plane strain assumptions yielded results closest to the full three-dimensional analysis. For computed transverse tensile stresses the plane stress assumption gave the best agreement. Based on this study it is recommended that results from plane stress and plane strain models be used as upper and lower bounds. The results from generalized plane strain models fall between the results obtained from plane stress and plane strain models. Two-dimensional models may also be used to qualitatively evaluate the stress distribution in a ply and the variation of energy release rates and mixed mode ratios with delamination length. For more accurate predictions, however, a three-dimensional analysis is required.

  19. Evaluation of 2D ceramic matrix composites in aeroconvective environments

    NASA Technical Reports Server (NTRS)

    Riccitiello, Salvatore R.; Love, Wendell L.; Balter-Peterson, Aliza

    1992-01-01

    An evaluation is conducted of a novel ceramic-matrix composite (CMC) material system for use in the aeroconvective-heating environments encountered by the nose caps and wing leading edges of such aerospace vehicles as the Space Shuttle, during orbit-insertion and reentry from LEO. These CMCs are composed of an SiC matrix that is reinforced with Nicalon, Nextel, or carbon refractory fibers in a 2D architecture. The test program conducted for the 2D CMCs gave attention to their subsurface oxidation.

  20. WFR-2D: an analytical model for PWAS-generated 2D ultrasonic guided wave propagation

    NASA Astrophysics Data System (ADS)

    Shen, Yanfeng; Giurgiutiu, Victor

    2014-03-01

    This paper presents WaveFormRevealer 2-D (WFR-2D), an analytical predictive tool for the simulation of 2-D ultrasonic guided wave propagation and interaction with damage. The design of structural health monitoring (SHM) systems and self-aware smart structures requires the exploration of a wide range of parameters to achieve best detection and quantification of certain types of damage. Such need for parameter exploration on sensor dimension, location, guided wave characteristics (mode type, frequency, wavelength, etc.) can be best satisfied with analytical models which are fast and efficient. The analytical model was constructed based on the exact 2-D Lamb wave solution using Bessel and Hankel functions. Damage effects were inserted in the model by considering the damage as a secondary wave source with complex-valued directivity scattering coefficients containing both amplitude and phase information from wave-damage interaction. The analytical procedure was coded with MATLAB, and a predictive simulation tool called WaveFormRevealer 2-D was developed. The wave-damage interaction coefficients (WDICs) were extracted from harmonic analysis of local finite element model (FEM) with artificial non-reflective boundaries (NRB). The WFR-2D analytical simulation results were compared and verified with full scale multiphysics finite element models and experiments with scanning laser vibrometer. First, Lamb wave propagation in a pristine aluminum plate was simulated with WFR-2D, compared with finite element results, and verified by experiments. Then, an inhomogeneity was machined into the plate to represent damage. Analytical modeling was carried out, and verified by finite element simulation and experiments. This paper finishes with conclusions and suggestions for future work.

  1. Homogenization models for 2-D grid structures

    NASA Technical Reports Server (NTRS)

    Banks, H. T.; Cioranescu, D.; Rebnord, D. A.

    1992-01-01

    In the past several years, we have pursued efforts related to the development of accurate models for the dynamics of flexible structures made of composite materials. Rather than viewing periodicity and sparseness as obstacles to be overcome, we exploit them to our advantage. We consider a variational problem on a domain that has large, periodically distributed holes. Using homogenization techniques we show that the solution to this problem is in some topology 'close' to the solution of a similar problem that holds on a much simpler domain. We study the behavior of the solution of the variational problem as the holes increase in number, but decrease in size in such a way that the total amount of material remains constant. The result is an equation that is in general more complex, but with a domain that is simply connected rather than perforated. We study the limit of the solution as the amount of material goes to zero. This second limit will, in most cases, retrieve much of the simplicity that was lost in the first limit without sacrificing the simplicity of the domain. Finally, we show that these results can be applied to the case of a vibrating Love-Kirchhoff plate with Kelvin-Voigt damping. We rely heavily on earlier results of (Du), (CS) for the static, undamped Love-Kirchhoff equation. Our efforts here result in a modification of those results to include both time dependence and Kelvin-Voigt damping.

  2. Integrated Design for Manufacturing of Braided Preforms for Advanced Composites Part I: 2D Braiding

    NASA Astrophysics Data System (ADS)

    Gao, Yan Tao; Ko, Frank K.; Hu, Hong

    2013-12-01

    This paper presents a 2D braiding design system for advanced textile structural composites was based on dynamic models. A software package to assist in the design of braided preform manufacturing has been developed. The package allows design parameters (machine speeds, fiber volume fraction, tightness factor, etc.) to be easily obtained and the relationships between said parameters to be demonstrated graphically. The fabirc geometry model (FGM) method was adopted to evaluate the mechanical properties of the composites. Experimental evidence demonstrates the success of the use of dynamic models in the design software for the manufacture of braided fabric preforms.

  3. 2D model of the Nucleus

    NASA Astrophysics Data System (ADS)

    Lach, Theodore M.

    2003-10-01

    The CBM (model) of the nucleus has resulted in the prediction of two new quarks, an "up" quark of mass 237.31 MeV/c2 and a "dn" quark of mass 42.392 MeV/c2. These two new predicted quarks helped to determine that the masses of the quarks and leptons are all related by a geometric progression relationship. The mass of each quark or lepton is just the "geometric mean" of two related elementary particles, either in the same generation or in the same family. This numerology predicts the following masses for the electron family: 0.511000 (electron), 7.74 (predicted), 117.3, 1778.4 (tau), 26950.1 MeV. The geometric ratio of this progression is 15.154 (e to the power e). The mass of the tau in this theory agrees very well with accepted values. This theory suggests that all the "dn like" quarks have a mass of just 10X multiples of 4.24 MeV (the mass of the "d" quark). The first 3 "up like" quark masses are 38, 237.31 and 1500 MeV. This theory also predicts a new heavy generation with a lepton mass of 27 GeV, a "dn like" quark of 42.4 GeV, and an "up like" quark of 65 GeV. Significant evidence already exists for the existence of these new quarks, and lepton. Ref. Masses of the Sub-Nuclear Particles, nucl-th/ 0008026, @ http://xxx.lanl.gov. Infinite Energy, Vol 5, issue 30.

  4. Technical Review of the UNET2D Hydraulic Model

    SciTech Connect

    Perkins, William A.; Richmond, Marshall C.

    2009-05-18

    The Kansas City District of the US Army Corps of Engineers is engaged in a broad range of river management projects that require knowledge of spatially-varied hydraulic conditions such as velocities and water surface elevations. This information is needed to design new structures, improve existing operations, and assess aquatic habitat. Two-dimensional (2D) depth-averaged numerical hydraulic models are a common tool that can be used to provide velocity and depth information. Kansas City District is currently using a specific 2D model, UNET2D, that has been developed to meet the needs of their river engineering applications. This report documents a tech- nical review of UNET2D.

  5. Influence of 2D Finite Element Modeling Assumptions on Debonding Prediction for Composite Skin-stiffener Specimens Subjected to Tension and Bending

    NASA Technical Reports Server (NTRS)

    Krueger, Ronald; Minguet, Pierre J.; Bushnell, Dennis M. (Technical Monitor)

    2002-01-01

    The influence of two-dimensional finite element modeling assumptions on the debonding prediction for skin-stiffener specimens was investigated. Geometrically nonlinear finite element analyses using two-dimensional plane-stress and plane strain elements as well as three different generalized plane strain type approaches were performed. The computed deflections, skin and flange strains, transverse tensile stresses and energy release rates were compared to results obtained from three-dimensional simulations. The study showed that for strains and energy release rate computations the generalized plane strain assumptions yielded results closest to the full three-dimensional analysis. For computed transverse tensile stresses the plane stress assumption gave the best agreement. Based on this study it is recommended that results from plane stress and plane strain models be used as upper and lower bounds. The results from generalized plane strain models fall between the results obtained from plane stress and plane strain models. Two-dimensional models may also be used to qualitatively evaluate the stress distribution in a ply and the variation of energy release rates and mixed mode ratios with lamination length. For more accurate predictions, however, a three-dimensional analysis is required.

  6. An Intercomparison of 2-D Models Within a Common Framework

    NASA Technical Reports Server (NTRS)

    Weisenstein, Debra K.; Ko, Malcolm K. W.; Scott, Courtney J.; Jackman, Charles H.; Fleming, Eric L.; Considine, David B.; Kinnison, Douglas E.; Connell, Peter S.; Rotman, Douglas A.; Bhartia, P. K. (Technical Monitor)

    2002-01-01

    A model intercomparison among the Atmospheric and Environmental Research (AER) 2-D model, the Goddard Space Flight Center (GSFC) 2-D model, and the Lawrence Livermore National Laboratory 2-D model allows us to separate differences due to model transport from those due to the model's chemical formulation. This is accomplished by constructing two hybrid models incorporating the transport parameters of the GSFC and LLNL models within the AER model framework. By comparing the results from the native models (AER and e.g. GSFC) with those from the hybrid model (e.g. AER chemistry with GSFC transport), differences due to chemistry and transport can be identified. For the analysis, we examined an inert tracer whose emission pattern is based on emission from a High Speed Civil Transport (HSCT) fleet; distributions of trace species in the 2015 atmosphere; and the response of stratospheric ozone to an HSCT fleet. Differences in NO(y) in the upper stratosphere are found between models with identical transport, implying different model representations of atmospheric chemical processes. The response of O3 concentration to HSCT aircraft emissions differs in the models from both transport-dominated differences in the HSCT-induced perturbations of H2O and NO(y) as well as from differences in the model represent at ions of O3 chemical processes. The model formulations of cold polar processes are found to be the most significant factor in creating large differences in the calculated ozone perturbations

  7. 2-D Magnetohydrodynamic Modeling of A Pulsed Plasma Thruster

    NASA Technical Reports Server (NTRS)

    Thio, Y. C. Francis; Cassibry, J. T.; Wu, S. T.; Rodgers, Stephen L. (Technical Monitor)

    2002-01-01

    Experiments are being performed on the NASA Marshall Space Flight Center (MSFC) MK-1 pulsed plasma thruster. Data produced from the experiments provide an opportunity to further understand the plasma dynamics in these thrusters via detailed computational modeling. The detailed and accurate understanding of the plasma dynamics in these devices holds the key towards extending their capabilities in a number of applications, including their applications as high power (greater than 1 MW) thrusters, and their use for producing high-velocity, uniform plasma jets for experimental purposes. For this study, the 2-D MHD modeling code, MACH2, is used to provide detailed interpretation of the experimental data. At the same time, a 0-D physics model of the plasma initial phase is developed to guide our 2-D modeling studies.

  8. 2 1/2 -D compressible reconnection model

    NASA Astrophysics Data System (ADS)

    Skender, M.; Vršnak, B.

    The exact solution of the jump conditions on the RD/SMS discontinuity system in a two-and-half-dimensional (2 1/2 -D) symmetrical reconnection model enables one to analyse the outflowing jet characteristics in dependence on the inflow velocity, and to follow changes in transition to the two-dimensional model. Implications arising from the exact solution and its relevance for solar flares are discussed.

  9. 2D face database diversification based on 3D face modeling

    NASA Astrophysics Data System (ADS)

    Wang, Qun; Li, Jiang; Asari, Vijayan K.; Karim, Mohammad A.

    2011-05-01

    Pose and illumination are identified as major problems in 2D face recognition (FR). It has been theoretically proven that the more diversified instances in the training phase, the more accurate and adaptable the FR system appears to be. Based on this common awareness, researchers have developed a large number of photographic face databases to meet the demand for data training purposes. In this paper, we propose a novel scheme for 2D face database diversification based on 3D face modeling and computer graphics techniques, which supplies augmented variances of pose and illumination. Based on the existing samples from identical individuals of the database, a synthesized 3D face model is employed to create composited 2D scenarios with extra light and pose variations. The new model is based on a 3D Morphable Model (3DMM) and genetic type of optimization algorithm. The experimental results show that the complemented instances obviously increase diversification of the existing database.

  10. Mechanical properties of 2D and 3D braided textile composites

    NASA Technical Reports Server (NTRS)

    Norman, Timothy L.

    1991-01-01

    The purpose of this research was to determine the mechanical properties of 2D and 3D braided textile composite materials. Specifically, those designed for tension or shear loading were tested under static loading to failure to investigate the effects of braiding. The overall goal of the work was to provide a structural designer with an idea of how textile composites perform under typical loading conditions. From test results for unnotched tension, it was determined that the 2D is stronger, stiffer, and has higher elongation to failure than the 3D. It was also found that the polyetherether ketone (PEEK) resin system was stronger, stiffer, and had higher elongation at failure than the resin transfer molding (RTM) epoxy. Open hole tension tests showed that PEEK resin is more notch sensitive than RTM epoxy. Of greater significance, it was found that the 3D is less notch sensitive than the 2D. Unnotched compression tests indicated, as did the tension tests, that the 2D is stronger, stiffer, and has higher elongation at failure than the RTM epoxy. The most encouraging results were from compression after impact. The 3D braided composite showed a compression after impact failure stress equal to 92 percent of the unimpacted specimen. The 2D braided composite failed at about 67 percent of the unimpacted specimen. Higher damage tolerance is observed in textiles over conventional composite materials. This is observed in the results, especially in the 3D braided materials.

  11. Modelling RF sources using 2-D PIC codes

    SciTech Connect

    Eppley, K.R.

    1993-03-01

    In recent years, many types of RF sources have been successfully modelled using 2-D PIC codes. Both cross field devices (magnetrons, cross field amplifiers, etc.) and pencil beam devices (klystrons, gyrotrons, TWT`S, lasertrons, etc.) have been simulated. All these devices involve the interaction of an electron beam with an RF circuit. For many applications, the RF structure may be approximated by an equivalent circuit, which appears in the simulation as a boundary condition on the electric field (``port approximation``). The drive term for the circuit is calculated from the energy transfer between beam and field in the drift space. For some applications it may be necessary to model the actual geometry of the structure, although this is more expensive. One problem not entirely solved is how to accurately model in 2-D the coupling to an external waveguide. Frequently this is approximated by a radial transmission line, but this sometimes yields incorrect results. We also discuss issues in modelling the cathode and injecting the beam into the PIC simulation.

  12. Modelling RF sources using 2-D PIC codes

    SciTech Connect

    Eppley, K.R.

    1993-03-01

    In recent years, many types of RF sources have been successfully modelled using 2-D PIC codes. Both cross field devices (magnetrons, cross field amplifiers, etc.) and pencil beam devices (klystrons, gyrotrons, TWT'S, lasertrons, etc.) have been simulated. All these devices involve the interaction of an electron beam with an RF circuit. For many applications, the RF structure may be approximated by an equivalent circuit, which appears in the simulation as a boundary condition on the electric field ( port approximation''). The drive term for the circuit is calculated from the energy transfer between beam and field in the drift space. For some applications it may be necessary to model the actual geometry of the structure, although this is more expensive. One problem not entirely solved is how to accurately model in 2-D the coupling to an external waveguide. Frequently this is approximated by a radial transmission line, but this sometimes yields incorrect results. We also discuss issues in modelling the cathode and injecting the beam into the PIC simulation.

  13. Unitary matrix models and 2D quantum gravity

    SciTech Connect

    Dalley, S. . Joseph Henry Labs.); Johnson, C.V.; Morris, T.R. . Dept. of Physics); Watterstam, A. )

    1992-09-21

    In this paper the KdV and modified KdV integrable hierarchies are shown to be different descriptions of the same 2D gravitational system - open-closed string theory. Non-perturbative solutions of the multicritical unitary matrix models map to non-singular solutions of the renormalization group equation for the string susceptibility, [P, Q] = Q. The authors also demonstrate that the large-N solutions of unitary matrix integrals in external fields, studied by Gross and Newman, equal the non-singular pure closed-string solutions of [[bar P], Q] = Q.

  14. Brane brick models and 2 d (0 , 2) triality

    NASA Astrophysics Data System (ADS)

    Franco, Sebastián; Lee, Sangmin; Seong, Rak-Kyeong

    2016-05-01

    We provide a brane realization of 2 d (0 , 2) Gadde-Gukov-Putrov triality in terms of brane brick models. These are Type IIA brane configurations that are T-dual to D1-branes over singular toric Calabi-Yau 4-folds. Triality translates into a local transformation of brane brick models, whose simplest representative is a cube move. We present explicit examples and construct their triality networks. We also argue that the classical mesonic moduli space of brane brick model theories, which corresponds to the probed Calabi-Yau 4-fold, is invariant under triality. Finally, we discuss triality in terms of phase boundaries, which play a central role in connecting Calabi-Yau 4-folds to brane brick models.

  15. 2D numerical modelling of meandering channel formation

    NASA Astrophysics Data System (ADS)

    XIAO, Y.; ZHOU, G.; YANG, F. S.

    2016-03-01

    A 2D depth-averaged model for hydrodynamic sediment transport and river morphological adjustment was established. The sediment transport submodel takes into account the influence of non-uniform sediment with bed surface armoring and considers the impact of secondary flow in the direction of bed-load transport and transverse slope of the river bed. The bank erosion submodel incorporates a simple simulation method for updating bank geometry during either degradational or aggradational bed evolution. Comparison of the results obtained by the extended model with experimental and field data, and numerical predictions validate that the proposed model can simulate grain sorting in river bends and duplicate the characteristics of meandering river and its development. The results illustrate that by using its control factors, the improved numerical model can be applied to simulate channel evolution under different scenarios and improve understanding of patterning processes.

  16. Experimental validation of 2D profile photoresist shrinkage model

    NASA Astrophysics Data System (ADS)

    Bunday, Benjamin; Cordes, Aaron; Self, Andy; Ferry, Lorena; Danilevsky, Alex

    2011-03-01

    For many years, lithographic resolution has been the main obstacle in allowing the pace of transistor densification to meet Moore's Law. For the 32 nm node and beyond, new lithography techniques will be used, including immersion ArF (iArF) lithography and extreme ultraviolet lithography (EUVL). As in the past, these techniques will use new types of photoresists with the capability to print smaller feature widths and pitches. These smaller feature sizes will also require the use of thinner layers of photoresists, such as under 100 nm. In previous papers, we focused on ArF and iArF photoresist shrinkage. We evaluated the magnitude of shrinkage for both R&D and mature resists as a function of chemical formulation, lithographic sensitivity, scanning electron microscope (SEM) beam condition, and feature size. Shrinkage results were determined by the well accepted methodology described in SEMATECH's CD-SEM Unified Specification. In other associated works, we first developed a 1-D model for resist shrinkage for the bottom linewidth and then a 2-D profile model that accounted for shrinkage of all aspects of a trapezoidal profile along a given linescan. A fundamental understanding of the phenomenology of the shrinkage trends was achieved, including how the shrinkage behaves differently for different sized and shaped features. In the 1-D case, calibration of the parameters to describe the photoresist material and the electron beam was all that was required to fit the models to real shrinkage data, as long as the photoresist was thick enough that the beam could not penetrate the entire layer of resist. The later 2-D model included improvements for solving the CD shrinkage in thin photoresists, which is now of great interest for upcoming realistic lithographic processing to explore the change in resist profile with electron dose and to predict the influence of initial resist profile on shrinkage characteristics. The 2-D model also included shrinkage due to both the primary

  17. A 2D channel-clogging biofilm model.

    PubMed

    Winstanley, H F; Chapwanya, M; Fowler, A C; O'Brien, S B G

    2015-09-01

    We present a model of biofilm growth in a long channel where the biomass is assumed to have the rheology of a viscous polymer solution. We examine the competition between growth and erosion-like surface detachment due to the flow. A particular focus of our investigation is the effect of the biofilm growth on the fluid flow in the pores, and the issue of whether biomass can grow sufficiently to shut off fluid flow through the pores, thus clogging the pore space. Net biofilm growth is coupled along the pore length via flow rate and nutrient transport in the pore flow. Our 2D model extends existing results on stability of 1D steady state biofilm thicknesses to show that, in the case of flows driven by a fixed pressure drop, full clogging of the pore can indeed happen in certain cases dependent on the functional form of the detachment term. PMID:25240390

  18. Mass loss in 2D rotating stellar models

    SciTech Connect

    Lovekin, Caterine; Deupree, Bob

    2010-10-05

    Radiatively driven mass loss is an important factor in the evolution of massive stars . The mass loss rates depend on a number of stellar parameters, including the effective temperature and luminosity. Massive stars are also often rapidly rotating, which affects their structure and evolution. In sufficiently rapidly rotating stars, both the effective temperature and radius vary significantly as a function of latitude, and hence mass loss rates can vary appreciably between the poles and the equator. In this work, we discuss the addition of mass loss to a 2D stellar evolution code (ROTORC) and compare evolution sequences with and without mass loss. Preliminary results indicate that a full 2D calculation of mass loss using the local effective temperature and luminosity can significantly affect the distribution of mass loss in rotating main sequence stars. More mass is lost from the pole than predicted by 1D models, while less mass is lost at the equator. This change in the distribution of mass loss will affect the angular momentum loss, the surface temperature and luminosity, and even the interior structure of the star. After a single mass loss event, these effects are small, but can be expected to accumulate over the course of the main sequence evolution.

  19. 2D Quantum Transport Modeling in Nanoscale MOSFETs

    NASA Technical Reports Server (NTRS)

    Svizhenko, Alexei; Anantram, M. P.; Govindan, T. R.; Biegel, Bryan

    2001-01-01

    With the onset of quantum confinement in the inversion layer in nanoscale MOSFETs, behavior of the resonant level inevitably determines all device characteristics. While most classical device simulators take quantization into account in some simplified manner, the important details of electrostatics are missing. Our work addresses this shortcoming and provides: (a) a framework to quantitatively explore device physics issues such as the source-drain and gate leakage currents, DIBL, and threshold voltage shift due to quantization, and b) a means of benchmarking quantum corrections to semiclassical models (such as density- gradient and quantum-corrected MEDICI). We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions, oxide tunneling and phase-breaking scattering are treated on equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. Quantum simulations are focused on MIT 25, 50 and 90 nm "well- tempered" MOSFETs and compared to classical and quantum corrected models. The important feature of quantum model is smaller slope of Id-Vg curve and consequently higher threshold voltage. These results are quantitatively consistent with I D Schroedinger-Poisson calculations. The effect of gate length on gate-oxide leakage and sub-threshold current has been studied. The shorter gate length device has an order of magnitude smaller current at zero gate bias than the longer gate length device without a significant trade-off in on-current. This should be a device design consideration.

  20. Cascading rainfall uncertainties into 2D inundation impact models

    NASA Astrophysics Data System (ADS)

    Souvignet, Maxime; de Almeida, Gustavo; Champion, Adrian; Garcia Pintado, Javier; Neal, Jeff; Freer, Jim; Cloke, Hannah; Odoni, Nick; Coxon, Gemma; Bates, Paul; Mason, David

    2013-04-01

    Existing precipitation products show differences in their spatial and temporal distribution and several studies have presented how these differences influence the ability to predict hydrological responses. However, an atmospheric-hydrologic-hydraulic uncertainty cascade is seldom explored and how, importantly, input uncertainties propagate through this cascade is still poorly understood. Such a project requires a combination of modelling capabilities, runoff generation predictions based on those rainfall forecasts, and hydraulic flood wave propagation based on the runoff predictions. Accounting for uncertainty in each component is important in decision making for issuing flood warnings, monitoring or planning. We suggest a better understanding of uncertainties in inundation impact modelling must consider these differences in rainfall products. This will improve our understanding of the input uncertainties on our predictive capability. In this paper, we propose to address this issue by i) exploring the effects of errors in rainfall on inundation predictive capacity within an uncertainty framework, i.e. testing inundation uncertainty against different comparable meteorological conditions (i.e. using different rainfall products). Our method cascades rainfall uncertainties into a lumped hydrologic model (FUSE) within the GLUE uncertainty framework. The resultant prediction uncertainties in discharge provide uncertain boundary conditions, which are cascaded into a simplified shallow water 2D hydraulic model (LISFLOOD-FP). Rainfall data captured by three different measurement techniques - rain gauges, gridded data and numerical weather predictions (NWP) models are used to assess the combined input data and model parameter uncertainty. The study is performed in the Severn catchment over the period between June and July 2007, where a series of rainfall events causing record floods in the study area). Changes in flood area extent are compared and the uncertainty envelope is

  1. Duality Between Spin Networks and the 2D Ising Model

    NASA Astrophysics Data System (ADS)

    Bonzom, Valentin; Costantino, Francesco; Livine, Etera R.

    2016-06-01

    The goal of this paper is to exhibit a deep relation between the partition function of the Ising model on a planar trivalent graph and the generating series of the spin network evaluations on the same graph. We provide respectively a fermionic and a bosonic Gaussian integral formulation for each of these functions and we show that they are the inverse of each other (up to some explicit constants) by exhibiting a supersymmetry relating the two formulations. We investigate three aspects and applications of this duality. First, we propose higher order supersymmetric theories that couple the geometry of the spin networks to the Ising model and for which supersymmetric localization still holds. Secondly, after interpreting the generating function of spin network evaluations as the projection of a coherent state of loop quantum gravity onto the flat connection state, we find the probability distribution induced by that coherent state on the edge spins and study its stationary phase approximation. It is found that the stationary points correspond to the critical values of the couplings of the 2D Ising model, at least for isoradial graphs. Third, we analyze the mapping of the correlations of the Ising model to spin network observables, and describe the phase transition on those observables on the hexagonal lattice. This opens the door to many new possibilities, especially for the study of the coarse-graining and continuum limit of spin networks in the context of quantum gravity.

  2. Effects of Agent's Repulsion in 2d Flocking Models

    NASA Astrophysics Data System (ADS)

    Moussa, Najem; Tarras, Iliass; Mazroui, M'hammed; Boughaleb, Yahya

    In nature many animal groups, such as fish schools or bird flocks, clearly display structural order and appear to move as a single coherent entity. In order to understand the complex behavior of these systems, many models have been proposed and tested so far. This paper deals with an extension of the Vicsek model, by including a second zone of repulsion, where each agent attempts to maintain a minimum distance from the others. The consideration of this zone in our study seems to play an important role during the travel of agents in the two-dimensional (2D) flocking models. Our numerical investigations show that depending on the basic ingredients such as repulsion radius (R1), effect of density of agents (ρ) and noise (η), our nonequilibrium system can undergo a kinetic phase transition from no transport to finite net transport. For different values of ρ, kinetic phase diagrams in the plane (η ,R1) are found. Implications of these findings are discussed.

  3. 2-D Model for Normal and Sickle Cell Blood Microcirculation

    NASA Astrophysics Data System (ADS)

    Tekleab, Yonatan; Harris, Wesley

    2011-11-01

    Sickle cell disease (SCD) is a genetic disorder that alters the red blood cell (RBC) structure and function such that hemoglobin (Hb) cannot effectively bind and release oxygen. Previous computational models have been designed to study the microcirculation for insight into blood disorders such as SCD. Our novel 2-D computational model represents a fast, time efficient method developed to analyze flow dynamics, O2 diffusion, and cell deformation in the microcirculation. The model uses a finite difference, Crank-Nicholson scheme to compute the flow and O2 concentration, and the level set computational method to advect the RBC membrane on a staggered grid. Several sets of initial and boundary conditions were tested. Simulation data indicate a few parameters to be significant in the perturbation of the blood flow and O2 concentration profiles. Specifically, the Hill coefficient, arterial O2 partial pressure, O2 partial pressure at 50% Hb saturation, and cell membrane stiffness are significant factors. Results were found to be consistent with those of Le Floch [2010] and Secomb [2006].

  4. Ab initio modeling of 2D layered organohalide lead perovskites.

    PubMed

    Fraccarollo, Alberto; Cantatore, Valentina; Boschetto, Gabriele; Marchese, Leonardo; Cossi, Maurizio

    2016-04-28

    A number of 2D layered perovskites A2PbI4 and BPbI4, with A and B mono- and divalent ammonium and imidazolium cations, have been modeled with different theoretical methods. The periodic structures have been optimized (both in monoclinic and in triclinic systems, corresponding to eclipsed and staggered arrangements of the inorganic layers) at the DFT level, with hybrid functionals, Gaussian-type orbitals and dispersion energy corrections. With the same methods, the various contributions to the solid stabilization energy have been discussed, separating electrostatic and dispersion energies, organic-organic intralayer interactions and H-bonding effects, when applicable. Then the electronic band gaps have been computed with plane waves, at the DFT level with scalar and full relativistic potentials, and including the correlation energy through the GW approximation. Spin orbit coupling and GW effects have been combined in an additive scheme, validated by comparing the computed gap with well known experimental and theoretical results for a model system. Finally, various contributions to the computed band gaps have been discussed on some of the studied systems, by varying some geometrical parameters and by substituting one cation in another's place. PMID:27131557

  5. 2-D Inhomogeneous Modeling of the Solar CO Bands

    NASA Astrophysics Data System (ADS)

    Ayres, T. R.

    1996-05-01

    The recent discovery of off-limb emissions in the mid-IR ( ~ 5 mu m) vibration-rotation bands of solar carbon monoxide (CO) has sparked new interest in the formation of the molecular lines, and their ability to diagnose thermal conditions at high altitudes. The off-limb extensions of the strong CO lines indicate the penetration of cool material (T ~ 3500 K) several hundred kilometers into the otherwise hot (T ~ 6000 K) chromosphere. The origin of the cool gas, and its role in the thermal energy balance, remain controversial. The interpretation of the CO observations must rely heavily upon numerical modeling, in particular highly-inhomogeneous thermal structures arrayed in a 2-D scheme that can properly treat the geometry of the grazing rays at the solar limb. The radiation transport, itself, is especially simple for the CO off-limb emissions, because the fundamental bands form quite close to LTE (high collision rates; low spontaneous decay rates) and the background continuum is purely thermal as well (f--f transitions in H(-) and H). Thus, the geometrical aspects of the problem can be treated in considerably more detail than would be practical for typical NLTE scattering lines. I describe the recent modeling efforts, and the diagnostic potential of the CO bands for future observational studies of inhomogeneous surface structure on the Sun, and on other stars of late spectral type. This work was supported by NSF grant AST-9218063 to the University of Colorado.

  6. Ab initio modeling of 2D layered organohalide lead perovskites

    NASA Astrophysics Data System (ADS)

    Fraccarollo, Alberto; Cantatore, Valentina; Boschetto, Gabriele; Marchese, Leonardo; Cossi, Maurizio

    2016-04-01

    A number of 2D layered perovskites A2PbI4 and BPbI4, with A and B mono- and divalent ammonium and imidazolium cations, have been modeled with different theoretical methods. The periodic structures have been optimized (both in monoclinic and in triclinic systems, corresponding to eclipsed and staggered arrangements of the inorganic layers) at the DFT level, with hybrid functionals, Gaussian-type orbitals and dispersion energy corrections. With the same methods, the various contributions to the solid stabilization energy have been discussed, separating electrostatic and dispersion energies, organic-organic intralayer interactions and H-bonding effects, when applicable. Then the electronic band gaps have been computed with plane waves, at the DFT level with scalar and full relativistic potentials, and including the correlation energy through the GW approximation. Spin orbit coupling and GW effects have been combined in an additive scheme, validated by comparing the computed gap with well known experimental and theoretical results for a model system. Finally, various contributions to the computed band gaps have been discussed on some of the studied systems, by varying some geometrical parameters and by substituting one cation in another's place.

  7. 2D modeling of electromagnetic waves in cold plasmas

    SciTech Connect

    Crombé, K.; Van Eester, D.; Koch, R.; Kyrytsya, V.

    2014-02-12

    The consequences of sheath (rectified) electric fields, resulting from the different mobility of electrons and ions as a response to radio frequency (RF) fields, are a concern for RF antenna design as it can cause damage to antenna parts, limiters and other in-vessel components. As a first step to a more complete description, the usual cold plasma dielectric description has been adopted, and the density profile was assumed to be known as input. Ultimately, the relevant equations describing the wave-particle interaction both on the fast and slow timescale will need to be tackled but prior to doing so was felt as a necessity to get a feeling of the wave dynamics involved. Maxwell's equations are solved for a cold plasma in a 2D antenna box with strongly varying density profiles crossing also lower hybrid and ion-ion hybrid resonance layers. Numerical modelling quickly becomes demanding on computer power, since a fine grid spacing is required to capture the small wavelengths effects of strongly evanescent modes.

  8. Numerical modeling of seismogram envelopes in 2-D random media

    NASA Astrophysics Data System (ADS)

    Fehler, Michael

    2002-11-01

    Several portions of seismograms recorded from regional earthquakes cannot be easily explained as resulting from waves propagating along deterministic paths within the Earth. For example, seismic coda, which is the tail portion of the seismogram of an earthquake recorded at distances of less than 100 km, is considered as resulting from waves that are multiply scattered from random heterogeneities in the Earth's lithosphere. At greater distances, observations that the duration of the initial arriving wave packet is much longer than the source-time duration is explained as being due to multiple forward scattering along the path between the source and the receiver. To investigate these phenomena, we use a finite difference method to numerically simulate 2-D scalar-waves that propagate through random media characterized by a von Karman autocorrelation function. Such media are considered to be appropriate models for the random component of the structure of the Earth's lithosphere. We investigate the characteristics of the resulting wavefields and compare them with those of observed seismograms.

  9. A 2D electrohydrodynamic model for electrorotation of fluid drops.

    PubMed

    Feng, James Q

    2002-02-01

    A theoretical analysis of spontaneous electrorotation of deformable fluid drops in a DC electric field is presented with a 2D electrohydrodynamic model. The fluids in the system are assumed to be leaky dielectric and Newtonian. If the rotating flow is dominant over the cellular convection type of electrohydrodynamic flow, closed-form solutions for drops of small deformations can be obtained. Because the governing equations are in general nonlinear even when drop deformations are ignored, the general solution for even undeformed drop takes a form of infinite series and can only be evaluated by numerical means. Both closed-form solutions for special cases and numerical solutions for more general cases are obtained here to describe steady-state field variables and first-order drop deformations. In a DC electric field of strength beyond the threshold value, spontaneous electrorotation of a drop is shown to occur when charge relaxation in the surrounding fluid is faster than the fluid inside the drop. With increasing the strength of the applied electric field from the threshold for onset of electrorotation, the axis of drop contraction deviates from from that of the applied electric field in the direction of the rotating flow with an angle increasing with the field strength. PMID:16290391

  10. VAM2D: Variably saturated analysis model in two dimensions

    SciTech Connect

    Huyakorn, P.S.; Kool, J.B.; Wu, Y.S. )

    1991-10-01

    This report documents a two-dimensional finite element model, VAM2D, developed to simulate water flow and solute transport in variably saturated porous media. Both flow and transport simulation can be handled concurrently or sequentially. The formulation of the governing equations and the numerical procedures used in the code are presented. The flow equation is approximated using the Galerkin finite element method. Nonlinear soil moisture characteristics and atmospheric boundary conditions (e.g., infiltration, evaporation and seepage face), are treated using Picard and Newton-Raphson iterations. Hysteresis effects and anisotropy in the unsaturated hydraulic conductivity can be taken into account if needed. The contaminant transport simulation can account for advection, hydrodynamic dispersion, linear equilibrium sorption, and first-order degradation. Transport of a single component or a multi-component decay chain can be handled. The transport equation is approximated using an upstream weighted residual method. Several test problems are presented to verify the code and demonstrate its utility. These problems range from simple one-dimensional to complex two-dimensional and axisymmetric problems. This document has been produced as a user's manual. It contains detailed information on the code structure along with instructions for input data preparation and sample input and printed output for selected test problems. Also included are instructions for job set up and restarting procedures. 44 refs., 54 figs., 24 tabs.

  11. A 2-D modeling contribution to river training design

    NASA Astrophysics Data System (ADS)

    Anselmo, V.; Coccato, M.; Frank, E.; Guiot, E.

    2003-04-01

    In the last ten years, two major floods (1994 and 2000) occurred in North-western Italy and a few questions arose about the hydraulic behavior of the streams as well about the choice and design of protection works. The River Po Authority is oriented to assign "design flows" in selected cross sections of the main rivers, as a design constraint to land management and river training in the upstream areas. Since the region has been fully developed in the last century and somewhere it is overcrowded, space for spreading flood flows is strongly reduced, while large partially developed areas are prone to flooding and residents ask for being protected. A first question regards the contribution to flood peak reduction of the still existing flood prone undeveloped areas beside the main channels, and a second question is about the best way to improve such a behavior. A 2-D unsteady model (Sobek, originated by Delft Hydraulics) was applied to a 25 km reach of the upper River Po. The effects of major floods was investigated, proving that the reduction of the peak flow is minor mainly because of the rather high slope (0.0015) and of the flood volume (500·106 m3). Aiming to enhance the role of the flooded areas, a few types of river training schemes were checked, with particular attention to the so called "Po system". Depth and extension of compartments are the main variables. Results are interesting, but must be evaluated in front of the cost-benefit analysis. The investigation is being extended to more steep stream reaches (up to 0.01), which are representative of the main upper Po tributaries.

  12. Predicting Fracture Using 2D Finite Element Modeling

    PubMed Central

    MacNeil, J.A.M.; Adachi, J.D; Goltzman, D; Josse, R.G; Kovacs, C.S; Prior, J.C; Olszynski, W; Davison, K.S.; Kaiser, S.M

    2013-01-01

    A decrease in bone density at the hip or spine has been shown to increase the risk of fracture. A limitation of the bone mineral density (BMD) measurement is that it provides only a measure of a bone samples average density when projected onto a 2D surface. Effectively, what determines bone fracture is whether an applied load exceeds ultimate strength, with both bone tissue material properties (can be approximated through bone density), and geometry playing a role. The goal of this project was to use bone geometry and BMD obtained from radiographs and DXA measurements respectively to estimate fracture risk, using a two-dimensional finite element model (FEM) of the sagittal plane of lumbar vertebrae. The Canadian Multicenter Osteoporosis Study (CaMos) data was used for this study. There were 4194 men and women over the age of 50 years, with 786 having fractures. Each subject had BMD testing and radiographs of their lumbar vertebrae. A single two dimensional FEM of the first to fourth lumbar vertebra was automatically generated for each subject. Bone tissue stiffness was assigned based on the BMD of the individual vertebrae, and adjusted for patient age. Axial compression boundary conditions were applied with a force proportional to body mass. The resulting overall strain from the applied force was found. Men and women were analyzed separately. At baseline, the sensitivity of BMD to predict fragility fractures in women and men was 3.77 % and 0.86 %, while the sensitivity of FEM to predict fragility fractures for women and men was 10.8 % and 11.3 %. The FEM ROC curve demonstrated better performance compared to BMD. The relative risk of being considered at high fracture risk using FEM at baseline, was a better predictor of 5 year incident fragility fracture risk compared to BMD. PMID:21959170

  13. A 2D simulation model for urban flood management

    NASA Astrophysics Data System (ADS)

    Price, Roland; van der Wielen, Jonathan; Velickov, Slavco; Galvao, Diogo

    2014-05-01

    The European Floods Directive, which came into force on 26 November 2007, requires member states to assess all their water courses and coast lines for risk of flooding, to map flood extents and assets and humans at risk, and to take adequate and coordinated measures to reduce the flood risk in consultation with the public. Flood Risk Management Plans are to be in place by 2015. There are a number of reasons for the promotion of this Directive, not least because there has been much urban and other infrastructural development in flood plains, which puts many at risk of flooding along with vital societal assets. In addition there is growing awareness that the changing climate appears to be inducing more frequent extremes of rainfall with a consequent increases in the frequency of flooding. Thirdly, the growing urban populations in Europe, and especially in the developing countries, means that more people are being put at risk from a greater frequency of urban flooding in particular. There are urgent needs therefore to assess flood risk accurately and consistently, to reduce this risk where it is important to do so or where the benefit is greater than the damage cost, to improve flood forecasting and warning, to provide where necessary (and possible) flood insurance cover, and to involve all stakeholders in decision making affecting flood protection and flood risk management plans. Key data for assessing risk are water levels achieved or forecasted during a flood. Such levels should of course be monitored, but they also need to be predicted, whether for design or simulation. A 2D simulation model (PriceXD) solving the shallow water wave equations is presented specifically for determining flood risk, assessing flood defense schemes and generating flood forecasts and warnings. The simulation model is required to have a number of important properties: -Solve the full shallow water wave equations using a range of possible solutions; -Automatically adjust the time step and

  14. Tension Strength, Failure Prediction and Damage Mechanisms in 2D Triaxial Braided Composites with Notch

    NASA Technical Reports Server (NTRS)

    Norman, Timothy L.; Anglin, Colin

    1995-01-01

    The unnotched and notched (open hole) tensile strength and failure mechanisms of two-dimensional (2D) triaxial braided composites were examined. The effect of notch size and notch position were investigated. Damage initiation and propagation in notched and unnotched coupons were also examined. Theory developed to predict the normal stress distribution near an open hole and failure for tape laminated composites was evaluated for its applicability to 2D triaxial braided textile composite materials. Four different fiber architectures were considered; braid angle, yarn and braider size, percentage of longitudinal yarns and braider angle varied. Tape laminates equivalent to textile composites were also constructed for comparison. Unnotched tape equivalents were stronger than braided textiles but exhibited greater notch sensitivity. Notched textiles and tape equivalents have roughly the same strength at large notch sizes. Two common damage mechanisms were found: braider yarn cracking and near notch longitudinal yarn splitting. Cracking was found to initiate in braider yarns in unnotched and notched coupons, and propagate in the direction of the braider yarns until failure. Damage initiation stress decreased with increasing braid angle. No significant differences in prediction of near notch strain between textile and tape equivalents could be detected for small braid angle, but the correlations were weak for textiles with large braid angle. Notch strength could not be predicted using existing anisotropic theory for braided textiles due to their insensitivity to notch.

  15. Experimental Investigation of the Strain Rate Dependent Behaviour of 2D Biaxially and Triaxially Reinforced Braided Composites

    NASA Astrophysics Data System (ADS)

    Böhm, R.; Hornig, A.; Luft, J.; Becker, M.; Koch, I.; Grüber, B.; Hufenbach, W.

    2014-04-01

    The performance of 2D biaxially and triaxially reinforced braided carbon fibre composites under dynamic loading is evaluated in the presented study. The accurate manufacturing of tensile specimen made of braided sleeves is explained particularly with regard to efficiency and reproducibility. In order to determine reliable strain rate dependent properties, the high-speed testing procedure is discussed. Using five materials, the parameter identification is described and relevant material data is provided. The measured stiffnesses and strengths are used to predict the non-linear stress-strain behaviour with an earlier proposed phenomenological damage model for textile composites. The gained orthotropic property-profile provides the input parameters for a numerical analysis of braided composite components using the calibrated model.

  16. A 2D simulation model for urban flood management

    NASA Astrophysics Data System (ADS)

    Price, Roland; van der Wielen, Jonathan; Velickov, Slavco; Galvao, Diogo

    2014-05-01

    The European Floods Directive, which came into force on 26 November 2007, requires member states to assess all their water courses and coast lines for risk of flooding, to map flood extents and assets and humans at risk, and to take adequate and coordinated measures to reduce the flood risk in consultation with the public. Flood Risk Management Plans are to be in place by 2015. There are a number of reasons for the promotion of this Directive, not least because there has been much urban and other infrastructural development in flood plains, which puts many at risk of flooding along with vital societal assets. In addition there is growing awareness that the changing climate appears to be inducing more frequent extremes of rainfall with a consequent increases in the frequency of flooding. Thirdly, the growing urban populations in Europe, and especially in the developing countries, means that more people are being put at risk from a greater frequency of urban flooding in particular. There are urgent needs therefore to assess flood risk accurately and consistently, to reduce this risk where it is important to do so or where the benefit is greater than the damage cost, to improve flood forecasting and warning, to provide where necessary (and possible) flood insurance cover, and to involve all stakeholders in decision making affecting flood protection and flood risk management plans. Key data for assessing risk are water levels achieved or forecasted during a flood. Such levels should of course be monitored, but they also need to be predicted, whether for design or simulation. A 2D simulation model (PriceXD) solving the shallow water wave equations is presented specifically for determining flood risk, assessing flood defense schemes and generating flood forecasts and warnings. The simulation model is required to have a number of important properties: -Solve the full shallow water wave equations using a range of possible solutions; -Automatically adjust the time step and

  17. Study of the mechanical behavior of a 2-D carbon-carbon composite

    NASA Technical Reports Server (NTRS)

    Avery, W. B.; Herakovich, C. T.

    1987-01-01

    The out-of-plane fracture of a 2-D carbon-carbon composite was observed and characterized to gain an understanding of the factors influencing the stress distribution in such a laminate. Finite element analyses of a two-ply carbon-carbon composite under in-plane, out-of-plane, and thermal loading were performed. Under in-plane loading all components of stress were strong functions of geometry. Additionally, large thermal stresses were predicted. Out-of-plane tensile tests revealed that failure was interlaminar, and that cracks propagated along the fiber-matrix interface. An elasticity solution was utilized to analyze an orthotropic fiber in an isotropic matrix under uniform thermal load. The analysis reveals that the stress distributions in a transversely orthotropic fiber are radically different than those predicted assuming the fiber to be transversely isotropic.

  18. The Implementation of C-ID, R2D2 Model on Learning Reading Comprehension

    ERIC Educational Resources Information Center

    Rayanto, Yudi Hari; Rusmawan, Putu Ngurah

    2016-01-01

    The purposes of this research are to find out, (1) whether C-ID, R2D2 model is effective to be implemented on learning Reading comprehension, (2) college students' activity during the implementation of C-ID, R2D2 model on learning Reading comprehension, and 3) college students' learning achievement during the implementation of C-ID, R2D2 model on…

  19. Image restoration using 2D autoregressive texture model and structure curve construction

    NASA Astrophysics Data System (ADS)

    Voronin, V. V.; Marchuk, V. I.; Petrosov, S. P.; Svirin, I.; Agaian, S.; Egiazarian, K.

    2015-05-01

    In this paper an image inpainting approach based on the construction of a composite curve for the restoration of the edges of objects in an image using the concepts of parametric and geometric continuity is presented. It is shown that this approach allows to restore the curved edges and provide more flexibility for curve design in damaged image by interpolating the boundaries of objects by cubic splines. After edge restoration stage, a texture restoration using 2D autoregressive texture model is carried out. The image intensity is locally modeled by a first spatial autoregressive model with support in a strongly causal prediction region on the plane. Model parameters are estimated by Yule-Walker method. Several examples considered in this paper show the effectiveness of the proposed approach for large objects removal as well as recovery of small regions on several test images.

  20. An Incompressible 2D Didactic Model with Singularity and Explicit Solutions of the 2D Boussinesq Equations

    NASA Astrophysics Data System (ADS)

    Chae, Dongho; Constantin, Peter; Wu, Jiahong

    2014-09-01

    We give an example of a well posed, finite energy, 2D incompressible active scalar equation with the same scaling as the surface quasi-geostrophic equation and prove that it can produce finite time singularities. In spite of its simplicity, this seems to be the first such example. Further, we construct explicit solutions of the 2D Boussinesq equations whose gradients grow exponentially in time for all time. In addition, we introduce a variant of the 2D Boussinesq equations which is perhaps a more faithful companion of the 3D axisymmetric Euler equations than the usual 2D Boussinesq equations.

  1. 2-D model of the streamer zone of a leader

    NASA Astrophysics Data System (ADS)

    Milikh, G. M.; Likhanskii, A. V.; Shneider, M. N.; Raina, A.; George, A.

    2016-02-01

    Formation of the streamer zone of a leader is an outstanding problem in the physics of electric discharges which is relevant to laboratory leaders, as well as to the leaders formed by lightning. Despite substantial progress in the theoretical understanding of this complicated phenomenon, significant puzzles, such as the low propagation velocity of a leader compared to the fast streamers, remain. The objective of this paper is to present 2-D plasma simulations of the formation and propagation of the streamer zone of a leader. In these simulations we will generate a group of streamers that propagate in a discharge gap while interacting with each other. It is shown that interaction between the streamers significantly reduces their propagation velocity. This explains why the leader, which consists of many streamers, is much slower than a single streamer formed in the same discharge gap. It is shown that the mean velocity suppression of the group of streamers is determined by the inter-streamer distance. The critical value of the packing factor of the streamers at which the interactions between them can be neglected, and thus the discussed process can be treated as caused by a single streamer, is obtained.

  2. 2D Transition-Metal-Dichalcogenide-Nanosheet-Based Composites for Photocatalytic and Electrocatalytic Hydrogen Evolution Reactions.

    PubMed

    Lu, Qipeng; Yu, Yifu; Ma, Qinglang; Chen, Bo; Zhang, Hua

    2016-03-01

    Hydrogen (H2) is one of the most important clean and renewable energy sources for future energy sustainability. Nowadays, photocatalytic and electrocatalytic hydrogen evolution reactions (HERs) from water splitting are considered as two of the most efficient methods to convert sustainable energy to the clean energy carrier, H2. Catalysts based on transition metal dichalcogenides (TMDs) are recognized as greatly promising substitutes for noble-metal-based catalysts for HER. The photocatalytic and electrocatalytic activities of TMD nanosheets for the HER can be further improved after hybridization with many kinds of nanomaterials, such as metals, oxides, sulfides, and carbon materials, through different methods including the in situ reduction method, the hot-injection method, the heating-up method, the hydro(solvo)thermal method, chemical vapor deposition (CVD), and thermal annealing. Here, recent progress in photocatalytic and electrocatalytic HERs using 2D TMD-based composites as catalysts is discussed. PMID:26676800

  3. 2D quantum double models from a 3D perspective

    NASA Astrophysics Data System (ADS)

    Bernabé Ferreira, Miguel Jorge; Padmanabhan, Pramod; Teotonio-Sobrinho, Paulo

    2014-09-01

    In this paper we look at three dimensional (3D) lattice models that are generalizations of the state sum model used to define the Kuperberg invariant of 3-manifolds. The partition function is a scalar constructed as a tensor network where the building blocks are tensors given by the structure constants of an involutory Hopf algebra A. These models are very general and are hard to solve in its entire parameter space. One can obtain familiar models, such as ordinary gauge theories, by letting A be the group algebra {C}(G) of a discrete group G and staying on a certain region of the parameter space. We consider the transfer matrix of the model and show that quantum double Hamiltonians are derived from a particular choice of the parameters. Such a construction naturally leads to the star and plaquette operators of the quantum double Hamiltonians, of which the toric code is a special case when A={C}({{{Z}}_{2}}). This formulation is convenient to study ground states of these generalized quantum double models where they can naturally be interpreted as tensor network states. For a surface Σ, the ground state degeneracy is determined by the Kuperberg 3-manifold invariant of \\Sigma \\times {{S}^{1}}. It is also possible to obtain extra models by simply enlarging the allowed parameter space but keeping the solubility of the model. While some of these extra models have appeared before in the literature, our 3D perspective allows for an uniform description of them.

  4. A fully coupled 2D model of equiaxed eutectic solidification

    SciTech Connect

    Charbon, Ch.; LeSar, R.

    1995-12-31

    We propose a model of equiaxed eutectic solidification that couples the macroscopic level of heat diffusion with the microscopic level of nucleation and growth of the eutectic grains. The heat equation with the source term corresponding to the latent heat release due to solidification is calculated numerically by means of an implicit finite difference method. In the time stepping scheme, the evolution of solid fraction is deduced from a stochastic model of nucleation and growth which uses the local temperature (interpolated from the FDM mesh) to determine the local grain density and the local growth rate. The solid-liquid interface of each grain is tracked by using a subdivision of each grain perimeter in a large number of sectors. The state of each sector (i.e. whether it is still in contact with the liquid or already captured by an other grain) and the increase of radius of each grain during one time step allows one to compute the increase of solid fraction. As for deterministic models, the results of the model are the evolution of temperature and of solid fraction at any point of the sample. Moreover the model provides a complete picture of the microstructure, thus not limiting the microstructural information to the average grain density but allowing one to compute any stereological value of interest. We apply the model to the solidification of gray cast iron.

  5. Improvement of a 2D numerical model of lava flows

    NASA Astrophysics Data System (ADS)

    Ishimine, Y.

    2013-12-01

    I propose an improved procedure that reduces an improper dependence of lava flow directions on the orientation of Digital Elevation Model (DEM) in two-dimensional simulations based on Ishihara et al. (in Lava Flows and Domes, Fink, JH eds., 1990). The numerical model for lava flow simulations proposed by Ishihara et al. (1990) is based on two-dimensional shallow water model combined with a constitutive equation for a Bingham fluid. It is simple but useful because it properly reproduces distributions of actual lava flows. Thus, it has been regarded as one of pioneer work of numerical simulations of lava flows and it is still now widely used in practical hazard prediction map for civil defense officials in Japan. However, the model include an improper dependence of lava flow directions on the orientation of DEM because the model separately assigns the condition for the lava flow to stop due to yield stress for each of two orthogonal axes of rectangular calculating grid based on DEM. This procedure brings a diamond-shaped distribution as shown in Fig. 1 when calculating a lava flow supplied from a point source on a virtual flat plane although the distribution should be circle-shaped. To improve the drawback, I proposed a modified procedure that uses the absolute value of yield stress derived from both components of two orthogonal directions of the slope steepness to assign the condition for lava flows to stop. This brings a better result as shown in Fig. 2. Fig. 1. (a) Contour plots calculated with the original model of Ishihara et al. (1990). (b) Contour plots calculated with a proposed model.

  6. Stratosphere chemistry in a 2-D model with residual circulation

    NASA Technical Reports Server (NTRS)

    Guthrie, Paul D.; Jackman, Charles H.

    1990-01-01

    The objective of this research was to examine the effects of chemical perturbations on the stratosphere using models which can incorporate fully interactive radiative, chemical, and dynamical responses, in the context of a zonally averaged model. Model runs for the unperturbed, chlorine-perturbed and simultaneously chlorine-and CO2-perturbed cases were completed using the JPL-87 chemical kinetics data. The base case was analyzed and submitted for publication. The perturbed cases show substantial sensitivity of the predicted column ozone depletion to the perturbations affecting lower stratosphere temperature, but less to far dynamical perturbations. The column ozone distribution changed substantially when the kinetics data was changed. This implies a greater-than-expected uncertainty in predicted latitude distributions of ozone depletion, due to uncertainty about the accuracy and completeness of the chemical kinetics data set.

  7. Development of CCHE2D embankment break model

    Technology Transfer Automated Retrieval System (TEKTRAN)

    Earthen embankment breach often results in detrimental impact on downstream residents and infrastructure, especially those located in the flooding zone. Embankment failures are most commonly caused by overtopping or internal erosion. This study is to develop a practical numerical model for simulat...

  8. Phase Structure of the Random Zq Models in 2D

    NASA Astrophysics Data System (ADS)

    Sasamoto, T.; Nishimori, H.

    We discuss the phase diagram of the random Zq models in two dimensions. It is argued that, when q is large enough, there exist three phases in the phase diagram with two axes being the temperature and the strength of randomness. Our conlusions are derived based on the application of the duality arguments for random systems, which have been formulated recently by Maillard et al.

  9. Simulation of subgrid orographic precipitation with an embedded 2-D cloud-resolving model

    NASA Astrophysics Data System (ADS)

    Jung, Joon-Hee; Arakawa, Akio

    2016-03-01

    By explicitly resolving cloud-scale processes with embedded two-dimensional (2-D) cloud-resolving models (CRMs), superparameterized global atmospheric models have successfully simulated various atmospheric events over a wide range of time scales. Up to now, however, such models have not included the effects of topography on the CRM grid scale. We have used both 3-D and 2-D CRMs to simulate the effects of topography with prescribed "large-scale" winds. The 3-D CRM is used as a benchmark. The results show that the mean precipitation can be simulated reasonably well by using a 2-D representation of topography as long as the statistics of the topography such as the mean and standard deviation are closely represented. It is also shown that the use of a set of two perpendicular 2-D grids can significantly reduce the error due to a 2-D representation of topography.

  10. Analytic Differentiation of Barlat's 2D Criteria for Inverse Modeling

    SciTech Connect

    Endelt, Benny; Nielsen, Karl Brian; Danckert, Joachim

    2005-08-05

    The demand for alternative identification schemes for identification of constitutive parameters is getting more pronounced as the complexity of the constitutive equations increases, i.e. the number of parameters subject to identification. A general framework for inverse identification of constitutive parameters associated with sheet metal forming is proposed in the article. The inverse problem is solved, through minimization of the least square error between an experimental punch force sampled from a deep drawing and a predicted punch force produced from a coherent finite element model.

  11. Conservation laws and LETKF with 2D Shallow Water Model

    NASA Astrophysics Data System (ADS)

    Zeng, Yuefei; Janjic, Tijana

    2016-04-01

    Numerous approaches have been proposed to maintain physical conservation laws in the numerical weather prediction models. However, to achieve a reliable prediction, adequate initial conditions are also necessary, which are produced by a data assimilation algorithm. If an ensemble Kalman filters (EnKF) is used for this purpose, it has been shown that it could yield unphysical analysis ensemble that for example violates principles of mass conservation and positivity preservation (e.g. Janjic et al 2014) . In this presentation, we discuss the selection of conservation criteria for the analysis step, and start with testing the conservation of mass, energy and enstrophy. The simple experiments deal with nonlinear shallow water equations and simulated observations that are assimilated with LETKF (Localized Ensemble Transform Kalman Filter, Hunt et al. 2007). The model is discretized in a specific way to conserve mass, angular momentum, energy and enstrophy. The effects of the data assimilation on the conserved quantities (of mass, energy and enstrophy) depend on observation covarage, localization radius, observed variable and observation operator. Having in mind that Arakawa (1966) and Arakawa and Lamb (1977) showed that the conservation of both kinetic energy and enstrophy by momentum advection schemes in the case of nondivergent flow prevents systematic and unrealistic energy cascade towards high wave numbers, a cause of excessive numerical noise and possible eventual nonlinear instability, we test the effects on prediction depending on the type of errors in the initial condition. The performance with respect to nonlinear energy cascade is assessed as well.

  12. Google Earth as a tool in 2-D hydrodynamic modeling

    NASA Astrophysics Data System (ADS)

    Chien, Nguyen Quang; Keat Tan, Soon

    2011-01-01

    A method for coupling virtual globes with geophysical hydrodynamic models is presented. Virtual globes such as Google TM Earth can be used as a visualization tool to help users create and enter input data. The authors discuss techniques for representing linear and areal geographical objects with KML (Keyhole Markup Language) files generated using computer codes (scripts). Although virtual globes offer very limited tools for data input, some data of categorical or vector type can be entered by users, and then transformed into inputs for the hydrodynamic program by using appropriate scripts. An application with the AnuGA hydrodynamic model was used as an illustration of the method. Firstly, users draw polygons on the Google Earth screen. These features are then saved in a KML file which is read using a script file written in the Lua programming language. After the hydrodynamic simulation has been performed, another script file is used to convert the resulting output text file to a KML file for visualization, where the depths of inundation are represented by the color of discrete point icons. The visualization of a wind speed vector field was also included as a supplementary example.

  13. Bond Order Correlations in the 2D Hubbard Model

    NASA Astrophysics Data System (ADS)

    Moore, Conrad; Abu Asal, Sameer; Yang, Shuxiang; Moreno, Juana; Jarrell, Mark

    We use the dynamical cluster approximation to study the bond correlations in the Hubbard model with next nearest neighbor (nnn) hopping to explore the region of the phase diagram where the Fermi liquid phase is separated from the pseudogap phase by the Lifshitz line at zero temperature. We implement the Hirsch-Fye cluster solver that has the advantage of providing direct access to the computation of the bond operators via the decoupling field. In the pseudogap phase, the parallel bond order susceptibility is shown to persist at zero temperature while it vanishes for the Fermi liquid phase which allows the shape of the Lifshitz line to be mapped as a function of filling and nnn hopping. Our cluster solver implements NVIDIA's CUDA language to accelerate the linear algebra of the Quantum Monte Carlo to help alleviate the sign problem by allowing for more Monte Carlo updates to be performed in a reasonable amount of computation time. Work supported by the NSF EPSCoR Cooperative Agreement No. EPS-1003897 with additional support from the Louisiana Board of Regents.

  14. A 2D model to design MHD induction pumps

    NASA Astrophysics Data System (ADS)

    Stieglitz, R.; Zeininger, J.

    2006-09-01

    Technical liquid metal systems accompanied by a thermal transfer of energy such as reactor systems, metallurgical processes, metal refinement, casting, etc., require a forced convection of the fluid. The increased temperatures and more often the environmental conditions as, e.g., in a nuclear environment, pumping principles are required, in which rotating parts are absent. Additionally, in many applications a controlled atmosphere is indispensable, in order to ensure the structural integrity of the duct walls. An interesting option to overcome the sealing problem of a mechanical pump towards the surrounding is offered by induction systems. Although their efficiency compared to that of turbo machines is quite low, they have several advantages, which are attractive to the specific requirements in liquid metal applications such as: - low maintenance costs due to the absence of sealings, bearings and moving parts; - low degradation rate of the structural material; - simple replacement of the inductor without cut of the piping system; - fine regulation of flow rate by different inductor connections; - change of pump characteristics without change of the mechanical set-up. Within the article, general design requirements of electromagnetic pumps (EMP) are elaborated. The design of two annular linear induction pumps operating with sodium and lead-bismuth are presented and the calculated pump characteristics and experimentally obtained data are compared. In this context, physical effects leading to deviations between the model and the real data are addressed. Finally, the main results are summarized. Tables 4, Figs 4, Refs 12.

  15. Stress analysis of a rectangular implant in laminated composites using 2-D and 3-D finite elements

    NASA Technical Reports Server (NTRS)

    Chow, Wai T.; Graves, Michael J.

    1992-01-01

    An analysis method using the FEM based on the Hellinger-Reissner variation principle has been developed to determine the 3-D stresses and displacements near a rectangular implant inside a laminated composite material. Three-dimensional elements are employed in regions where the interlaminar stress is considered to be significant; 2-D elements are used in other areas. Uniaxially loaded graphite-epoxy laminates have been analyzed; the implant was modeled as four plies of 3501/6 epoxy located in the middle of the laminate. It is shown that the interlaminar stresses are an order of magnitude lower than the stress representing the applied far-field load. The stress concentration factors of both the interlaminar and in-plane stresses depend on the stacking sequence of the laminate.

  16. GEO2D - Two-Dimensional Computer Model of a Ground Source Heat Pump System

    DOE Data Explorer

    James Menart

    2013-06-07

    This file contains a zipped file that contains many files required to run GEO2D. GEO2D is a computer code for simulating ground source heat pump (GSHP) systems in two-dimensions. GEO2D performs a detailed finite difference simulation of the heat transfer occurring within the working fluid, the tube wall, the grout, and the ground. Both horizontal and vertical wells can be simulated with this program, but it should be noted that the vertical wall is modeled as a single tube. This program also models the heat pump in conjunction with the heat transfer occurring. GEO2D simulates the heat pump and ground loop as a system. Many results are produced by GEO2D as a function of time and position, such as heat transfer rates, temperatures and heat pump performance. On top of this information from an economic comparison between the geothermal system simulated and a comparable air heat pump systems or a comparable gas, oil or propane heating systems with a vapor compression air conditioner. The version of GEO2D in the attached file has been coupled to the DOE heating and cooling load software called ENERGYPLUS. This is a great convenience for the user because heating and cooling loads are an input to GEO2D. GEO2D is a user friendly program that uses a graphical user interface for inputs and outputs. These make entering data simple and they produce many plotted results that are easy to understand. In order to run GEO2D access to MATLAB is required. If this program is not available on your computer you can download the program MCRInstaller.exe, the 64 bit version, from the MATLAB website or from this geothermal depository. This is a free download which will enable you to run GEO2D..

  17. Comparison of 1D and 2D modelling with soil erosion model SMODERP

    NASA Astrophysics Data System (ADS)

    Kavka, Petr; Weyskrabova, Lenka; Zajicek, Jan

    2013-04-01

    The contribution presents a comparison of a runoff simulated by profile method (1D) and spatially distributed method (2D). Simulation model SMODERP is used for calculation and prediction of soil erosion and surface runoff from agricultural land. SMODERP is physically based model that includes the processes of infiltration (Phillips equation), surface runoff (kinematic wave based equation), surface retention, surface roughness and vegetation impact on runoff. 1D model was developed in past, new 2D model was developed in last two years. The model is being developed at the Department of Irrigation, Drainage and Landscape Engineering, Civil Engineering Faculty, CTU in Prague. 2D model was developed as a tool for widespread GIS software ArcGIS. The physical relations were implemented through Python script. This script uses ArcGIS system tools for raster and vectors treatment of the inputs. Flow direction is calculated by Steepest Descent algorithm in the preliminary version of 2D model. More advanced multiple flow algorithm is planned in the next version. Spatially distributed models enable to estimate not only surface runoff but also flow in the rills. Surface runoff is described in the model by kinematic wave equation. Equation uses Manning roughness coefficient for surface runoff. Parameters for five different soil textures were calibrated on the set of forty measurements performed on the laboratory rainfall simulator. For modelling of the rills a specific sub model was created. This sub model uses Manning formula for flow estimation. Numerical stability of the model is solved by Courant criterion. Spatial scale is fixed. Time step is dynamically changed depending on how flow is generated and developed. SMODERP is meant to be used not only for the research purposes, but mainly for the engineering practice. We also present how the input data can be obtained based on available resources (soil maps and data, land use, terrain models, field research, etc.) and how can

  18. 2D-photochemical modeling of Saturn’s stratosphere: hydrocarbon and water distributions

    NASA Astrophysics Data System (ADS)

    Hue, Vincent; Cavalié, Thibault; Hersant, Franck; Dobrijevic, Michel; Greathouse, Thomas; Lellouch, Emmanuel; Hartogh, Paul; Cassidy, Timothy; Spiga, Aymeric; Guerlet, Sandrine; Sylvestre, Melody

    2014-11-01

    Saturn’s axial tilt of 27° produces seasons in a similar way as on Earth. The seasonal forcing over Saturn’s 30 years period influences the production/loss of the major atmospheric absorbers and coolants through photochemistry, and influences therefore Saturn’s stratospheric temperatures. We have developed a 2D time-dependent photochemical model of Saturn’s atmosphere [Hue et al., in prep.], coupled to a radiative-climate model [Greathouse et al., 2008] to study seasonal effects on its atmospheric composition. Cassini spacecraft has revealed that the distribution of hydrocarbons in Saturn’s stratosphere [Guerlet et al., 2009] differs from pure photochemical predictions, i.e. without meridional transport [Moses et al., 2005]. Differences between the observed distribution of hydrocarbons and 2D-photochemical predictions are likely to be an indicator of dynamical forcing.Disentangling the origin of water in the stratosphere of this planet has been a long-term issue. Due to Saturn’s cold tropopause trap, which acts as a transport barrier, the water vapor observed by the Infrared Space Observatory (ISO) [Feuchtgruber et al., 1997] has an external origin. Three external sources have been identified: (i) permanent flux from interplanetary dust particles, (ii) local sources form planetary environments (rings, satellites), (iii) large cometary impacts, similar to Shoemaker-Levy 9 on Jupiter. Previous observations of Saturn with Herschel’s Hsso program [Hartogh et al., 2009] led to the detection of a water torus around Saturn [Hartogh et al., 2011], fed by Enceladus’ geysers. A substantial fraction of this torus is predicted to be a local source of water for Saturn’s and its satellites, as it will spread in this system [Cassidy et al., 2010]. Using the new 2D-photochemical model, we test here the validity of Enceladus’ torus as the source of Saturn’s stratospheric water.References : Hue et al., in prep. Greathouse et al., 2008. AGU Fall Meeting

  19. Matrix Cracking in Four Different 2D SiC/SiC Composite Systems

    NASA Technical Reports Server (NTRS)

    Morscher, Gregory N.

    2003-01-01

    Silicon carbide fiber reinforced, silicon carbide matrix composites are some of the most advanced composite systems for high-temperature, high-stress applications in oxidizing environments. A basic area that needs to be understood for the purpose of material behavior modeling and optimization is the architectural, constituent, and mechanistic factors that contribute to non-linear stress-strain behavior. The mechanism that causes non-linear stress-strain in dense-matrix composites is the formation and propagation of bridged matrix cracks. In addition, the occurrence and propagation of matrix cracks controls the time-dependent strength-properties of these materials in oxidizing environments at elevated temperatures. A modal acoustic emission technique has been used to monitor and estimate the stress-dependent matrix cracking. Two different SiC matrix systems, chemical vapor infiltrated (CVI) and melt-infiltrated (MI), with two different SiC fiber reinforcement, Hi-Nicalon (trademark) and Sylramic (trademark) were compared. Even though the averages of the range where matrix cracking occurred for the composites varied by more than 0.1% in strain and almost 200 MPa in stress, the range or distribution for matrix cracking could be reduced to a narrow band of stress for CVI SiC and MI SiC composites if it were assumed that all matrix cracks emanate outside of the load-bearing fiber, interphase, CVI preform minicomposite. A simple relationship was determined to describe stress-dependent matrix cracking which can then be used to estimate the onset of large, bridged matrix cracks or for material behavior models.

  20. Spot size variation FCS in simulations of the 2D Ising model

    NASA Astrophysics Data System (ADS)

    Burns, Margaret C.; Nouri, Mariam; Veatch, Sarah L.

    2016-06-01

    Spot variation fluorescence correlation spectroscopy (svFCS) was developed to study the movement and organization of single molecules in plasma membranes. This experimental technique varies the size of an illumination area while measuring correlations in time using standard fluorescence correlation methods. Frequently, this data is interpreted using the assumption that correlation measurements reflect the dynamics of single molecule motions, and not motions of the average composition. Here, we explore how svFCS measurements report on the dynamics of components diffusing within simulations of a 2D Ising model with a conserved order parameter. Simulated correlation functions report on both the fast dynamics of single component mobility and the slower dynamics of the average composition. Over a range of simulation conditions, a conventional svFCS analysis suggests the presence of anomalous diffusion even though single molecule motions are nearly Brownian in these simulations. This misinterpretation is most significant when the surface density of the fluorescent label is elevated, therefore we suggest future measurements be made over a range of tracer densities. Some simulation conditions reproduce qualitative features of published svFCS experimental data. Overall, this work emphasizes the need to probe membranes using multiple complimentary experimental methodologies in order to draw conclusions regarding the nature of spatial and dynamical heterogeneity in these systems.

  1. Analysis of vegetation effect on waves using a vertical 2-D RANS model

    Technology Transfer Automated Retrieval System (TEKTRAN)

    A vertical two-dimensional (2-D) model has been applied in the simulation of wave propagation through vegetated water bodies. The model is based on an existing model SOLA-VOF which solves the Reynolds-Averaged Navier-Stokes (RANS) equations with the finite difference method on a staggered rectangula...

  2. Simulation of Cardiac Arrhythmias Using a 2D Heterogeneous Whole Heart Model.

    PubMed

    Balakrishnan, Minimol; Chakravarthy, V Srinivasa; Guhathakurta, Soma

    2015-01-01

    Simulation studies of cardiac arrhythmias at the whole heart level with electrocardiogram (ECG) gives an understanding of how the underlying cell and tissue level changes manifest as rhythm disturbances in the ECG. We present a 2D whole heart model (WHM2D) which can accommodate variations at the cellular level and can generate the ECG waveform. It is shown that, by varying cellular-level parameters like the gap junction conductance (GJC), excitability, action potential duration (APD) and frequency of oscillations of the auto-rhythmic cell in WHM2D a large variety of cardiac arrhythmias can be generated including sinus tachycardia, sinus bradycardia, sinus arrhythmia, sinus pause, junctional rhythm, Wolf Parkinson White syndrome and all types of AV conduction blocks. WHM2D includes key components of the electrical conduction system of the heart like the SA (Sino atrial) node cells, fast conducting intranodal pathways, slow conducting atriovenctricular (AV) node, bundle of His cells, Purkinje network, atrial, and ventricular myocardial cells. SA nodal cells, AV nodal cells, bundle of His cells, and Purkinje cells are represented by the Fitzhugh-Nagumo (FN) model which is a reduced model of the Hodgkin-Huxley neuron model. The atrial and ventricular myocardial cells are modeled by the Aliev-Panfilov (AP) two-variable model proposed for cardiac excitation. WHM2D can prove to be a valuable clinical tool for understanding cardiac arrhythmias. PMID:26733873

  3. Simulation of Cardiac Arrhythmias Using a 2D Heterogeneous Whole Heart Model

    PubMed Central

    Balakrishnan, Minimol; Chakravarthy, V. Srinivasa; Guhathakurta, Soma

    2015-01-01

    Simulation studies of cardiac arrhythmias at the whole heart level with electrocardiogram (ECG) gives an understanding of how the underlying cell and tissue level changes manifest as rhythm disturbances in the ECG. We present a 2D whole heart model (WHM2D) which can accommodate variations at the cellular level and can generate the ECG waveform. It is shown that, by varying cellular-level parameters like the gap junction conductance (GJC), excitability, action potential duration (APD) and frequency of oscillations of the auto-rhythmic cell in WHM2D a large variety of cardiac arrhythmias can be generated including sinus tachycardia, sinus bradycardia, sinus arrhythmia, sinus pause, junctional rhythm, Wolf Parkinson White syndrome and all types of AV conduction blocks. WHM2D includes key components of the electrical conduction system of the heart like the SA (Sino atrial) node cells, fast conducting intranodal pathways, slow conducting atriovenctricular (AV) node, bundle of His cells, Purkinje network, atrial, and ventricular myocardial cells. SA nodal cells, AV nodal cells, bundle of His cells, and Purkinje cells are represented by the Fitzhugh-Nagumo (FN) model which is a reduced model of the Hodgkin-Huxley neuron model. The atrial and ventricular myocardial cells are modeled by the Aliev-Panfilov (AP) two-variable model proposed for cardiac excitation. WHM2D can prove to be a valuable clinical tool for understanding cardiac arrhythmias. PMID:26733873

  4. Fast 2D flood modelling using GPU technology - recent applications and new developments

    NASA Astrophysics Data System (ADS)

    Crossley, Amanda; Lamb, Rob; Waller, Simon; Dunning, Paul

    2010-05-01

    In recent years there has been considerable interest amongst scientists and engineers in exploiting the potential of commodity graphics hardware for desktop parallel computing. The Graphics Processing Units (GPUs) that are used in PC graphics cards have now evolved into powerful parallel co-processors that can be used to accelerate the numerical codes used for floodplain inundation modelling. We report in this paper on experience over the past two years in developing and applying two dimensional (2D) flood inundation models using GPUs to achieve significant practical performance benefits. Starting with a solution scheme for the 2D diffusion wave approximation to the 2D Shallow Water Equations (SWEs), we have demonstrated the capability to reduce model run times in ‘real-world' applications using GPU hardware and programming techniques. We then present results from a GPU-based 2D finite volume SWE solver. A series of numerical test cases demonstrate that the model produces outputs that are accurate and consistent with reference results published elsewhere. In comparisons conducted for a real world test case, the GPU-based SWE model was over 100 times faster than the CPU version. We conclude with some discussion of practical experience in using the GPU technology for flood mapping applications, and for research projects investigating use of Monte Carlo simulation methods for the analysis of uncertainty in 2D flood modelling.

  5. Dynamic Linkages Between the Transition Zone & Surface Plate Motions in 2D Models of Subduction

    NASA Astrophysics Data System (ADS)

    Arredondo, K.; Billen, M. I.

    2012-12-01

    Descending subducted slabs affect both plate tectonics at the surface and overall mantle flow (e.g. Conrad and Lithgow-Bertelloni, 2002). For time-dependent numerical models, the potential evolution of these slabs, ranging from immediate penetration into the lower mantle to prior buckling and stagnation, are affected by parameters such as the plate age, the viscosity jump into the lower mantle, the presence of phase transitions, trench motion and the chosen governing equation approximation (e.g. Billen and Hirth, 2007). Similarly, the overall deviatoric stress within the slab, especially where modified by the phase transitions, may explain the uneven distribution of deep earthquakes with depth (e.g. Bina, 1997). Better understanding of these processes may arise from a more realistic 2-D model that is fully-dynamic, with an overriding plate, freely-moving trench, compositionally-layered slab and seven major phase transitions, in addition to using the compressible (TALA) form of the governing equations. Though the thermodynamic parameters of certain phase transitions may be uncertain, this study aims to test the latest data and encourage further mineralogical research. We will present fully-dynamic models, which explore the importance of the phase transitions, especially those that have been previously excluded such as the wadsleyite to ringwoodite and the pyroxene and garnet phase transitions. These phase transitions, coupled with the modeled compositionally distinct crust, harzburgite, and pyrolite lithosphere layers, may produce new large-scale dynamic behavior not seen in past numerical models, as well as stress variations within the slab related to deep slab seismicity. Feedback from the compositionally complex slab to the dynamic trench may provide further insight on the mechanics of slab stagnation and behavior in the upper and lower mantle. Billen, M. I., and G. Hirth, Rheologic controls on slab dynamics, Geochemistry, Geophysics and Geosystems, 8 (Q08012

  6. A multispeed Discrete Boltzmann Model for transcritical 2D shallow water flows

    NASA Astrophysics Data System (ADS)

    La Rocca, Michele; Montessori, Andrea; Prestininzi, Pietro; Succi, Sauro

    2015-03-01

    In this work a Discrete Boltzmann Model for the solution of transcritical 2D shallow water flows is presented and validated. In order to provide the model with transcritical capabilities, a particular multispeed velocity set has been employed for the discretization of the Boltzmann equation. It is shown that this particular set naturally yields a simple and closed procedure to determine higher order equilibrium distribution functions needed to simulate transcritical flow. The model is validated through several classical benchmarks and is proven to correctly and accurately simulate both 1D and 2D transitions between the two flow regimes.

  7. A 2-D dynamical model of mesospheric temperature inversions in winter

    SciTech Connect

    Hauchecorne, A.; Maillard, A. )

    1990-11-01

    A 2-D stratospheric and mesospheric dynamical model including drag and diffusion due to gravity wave breaking is used to simulate winter mesospheric temperature inversions similar to those observed by Rayleigh lidar. It is shown that adiabatic heating associated to descending velocities in the mesosphere is the main mechanism involved in the formation of such inversions. Sensitivity tests are performed with the model and confirm this assumption. It is also explained why other previous similar studies with 2-D models did not show mesospheric inversion layers.

  8. Merging of RVR meander with CONCEPTS: Simplified 2D model for long-term meander evolution

    Technology Transfer Automated Retrieval System (TEKTRAN)

    RVR Meander is a simplified two-dimensional (2D) hydrodynamic and migration model (Abad and Garcia, 2006) while CONCEPTS (CONservational Channel Evolution and Pollutant Transport System) is a one-dimensional (1D) hydrodynamic and morphodynamic model (Langendoen and Alonso, 2008; Langendoen and Simon...

  9. Introducing the R2D2 Model: Online Learning for the Diverse Learners of This World

    ERIC Educational Resources Information Center

    Bonk, Curtis J.; Zhang, Ke

    2006-01-01

    The R2D2 method--read, reflect, display, and do--is a new model for designing and delivering distance education, and in particular, online learning. Such a model is especially important to address the diverse preferences of online learners of varied generations and varied Internet familiarity. Four quadrants can be utilized separately or as part…

  10. Evaluation of 2D shallow-water model for spillway flow with a complex geometry

    Technology Transfer Automated Retrieval System (TEKTRAN)

    Although the two-dimensional (2D) shallow water model is formulated based on several assumptions such as hydrostatic pressure distribution and vertical velocity is negligible, as a simple alternative to the complex 3D model, it has been used to compute water flows in which these assumptions may be ...

  11. Knight shift and spin relaxation in the single band 2D Hubbard model

    NASA Astrophysics Data System (ADS)

    Leblanc, James; Chen, Xi; Gull, Emanuel

    We study in detail the roles of spin and charge fluctuations in the single band 2D Hubbard model. Using dynamical mean field theory and cluster extensions such as the dynamical cluster approximation (DCA), we compute the full two particle susceptibilities in the spin and charge representations. By performing analytic continuations we obtain the temperature and doping dependence of the spin-lattice relaxation (T1- 1) and knight shift in the 2D Hubbard model relevant to NMR results on doped cuprates and connect these to RPA results in weak coupling limits.

  12. Error propagation for velocity and shear stress prediction using 2D models for environmental management

    NASA Astrophysics Data System (ADS)

    Pasternack, Gregory B.; Gilbert, Andrew T.; Wheaton, Joseph M.; Buckland, Evan M.

    2006-08-01

    SummaryResource managers, scientists, government regulators, and stakeholders are considering sophisticated numerical models for managing complex environmental problems. In this study, observations from a river-rehabilitation experiment involving gravel augmentation and spawning habitat enhancement were used to assess sources and magnitudes of error in depth, velocity, and shear velocity predictions made at the 1-m scale with a commercial two-dimensional (depth-averaged) model. Error in 2D model depth prediction averaged 21%. This error was attributable to topographic survey resolution, which at 1 pt per 1.14 m 2, was inadequate to resolve small humps and depressions influencing point measurements. Error in 2D model velocity prediction averaged 29%. More than half of this error was attributable to depth prediction error. Despite depth and velocity error, 56% of tested 2D model predictions of shear velocity were within the 95% confidence limit of the best field-based estimation method. Ninety percent of the error in shear velocity prediction was explained by velocity prediction error. Multiple field-based estimates of shear velocity differed by up to 160%, so the lower error for the 2D model's predictions suggests such models are at least as accurate as field measurement. 2D models enable detailed, spatially distributed estimates compared to the small number measurable in a field campaign of comparable cost. They also can be used for design evaluation. Although such numerical models are limited to channel types adhering to model assumptions and yield predictions only accurate to ˜20-30%, they can provide a useful tool for river-rehabilitation design and assessment, including spatially diverse habitat heterogeneity as well as for pre- and post-project appraisal.

  13. 2D condensation model for the inner Solar Nebula: an enstatite-rich environment

    NASA Astrophysics Data System (ADS)

    Pignatale, F. C.; Liffman, Kurt; Maddison, Sarah T.; Brooks, Geoffrey

    2016-04-01

    Infrared observations provide the dust composition in the protoplanetary discs surface layers, but cannot probe the dust chemistry in the mid-plane, where planet formation occurs. Meteorites show that dynamics was important in determining the dust distribution in the Solar Nebula and needs to be considered if we are to understand the global chemistry in discs. 1D radial condensation sequences can only simulate one disc layer at a time and cannot describe the global chemistry or the complexity of meteorites. To address these limitations, we compute for the first time the 2D distribution of condensates in the inner Solar Nebula using a thermodynamic equilibrium model, and derive time-scales for vertical settling and radial migration of dust. We find two enstatite-rich zones within 1 AU from the young Sun: a band ˜0.1 AU thick in the upper optically-thin layer of the disc interior to 0.8 AU, and in the optically-thick disc mid-plane out to ˜0.4 AU. The two enstatite-rich zones support recent evidence that Mercury and enstatite chondrites (ECs) shared a bulk material with similar composition. Our results are also consistent with infrared observation of protoplanetary disc which show emission of enstatite-rich dust in the inner surface of discs. The resulting chemistry and dynamics suggests that the formation of the bulk material of ECs occurred in the inner surface layer of the disc, within 0.4 AU. We also propose a simple alternative scenario in which gas fractionation and vertical settling of the condensates lead to an enstatite-chondritic bulk material.

  14. Validation of DYSTOOL for unsteady aerodynamic modeling of 2D airfoils

    NASA Astrophysics Data System (ADS)

    González, A.; Gomez-Iradi, S.; Munduate, X.

    2014-06-01

    From the point of view of wind turbine modeling, an important group of tools is based on blade element momentum (BEM) theory using 2D aerodynamic calculations on the blade elements. Due to the importance of this sectional computation of the blades, the National Renewable Wind Energy Center of Spain (CENER) developed DYSTOOL, an aerodynamic code for 2D airfoil modeling based on the Beddoes-Leishman model. The main focus here is related to the model parameters, whose values depend on the airfoil or the operating conditions. In this work, the values of the parameters are adjusted using available experimental or CFD data. The present document is mainly related to the validation of the results of DYSTOOL for 2D airfoils. The results of the computations have been compared with unsteady experimental data of the S809 and NACA0015 profiles. Some of the cases have also been modeled using the CFD code WMB (Wind Multi Block), within the framework of a collaboration with ACCIONA Windpower. The validation has been performed using pitch oscillations with different reduced frequencies, Reynolds numbers, amplitudes and mean angles of attack. The results have shown a good agreement using the methodology of adjustment for the value of the parameters. DYSTOOL have demonstrated to be a promising tool for 2D airfoil unsteady aerodynamic modeling.

  15. 2D-photochemical model for forbidden oxygen line emission for comet 1P/Halley

    NASA Astrophysics Data System (ADS)

    Cessateur, G.; De Keyser, J.; Maggiolo, R.; Rubin, M.; Gronoff, G.; Gibbons, A.; Jehin, E.; Dhooghe, F.; Gunell, H.; Vaeck, N.; Loreau, J.

    2016-08-01

    We present here a 2D-model of photochemistry for computing the production and loss mechanisms of the O(1S) and O(1D) states, which are responsible for the emission lines at 577.7 nm, 630 nm, and 636.4 nm, in case of the comet 1P/Halley. The presence of O2 within cometary atmospheres, measured by the in-situ ROSETTA and GIOTTO missions, necessitates a revision of the usual photochemical models. Indeed, the photodissociation of molecular oxygen also leads to a significant production of oxygen in excited electronic states. In order to correctly model the solar UV flux absorption, we consider here a 2D configuration. While the green to red-doublet ratio is not affected by the solar UV flux absorption, estimates of the red-doublet and green lines emissions are, however, overestimated by a factor of two in the 1D model compared to the 2D model. Considering a spherical symmetry, emission maps can be deduced from the 2D model in order to be directly compared to ground and/or in-situ observations.

  16. Oriented Gaussian mixture models for nonrigid 2D/3D coronary artery registration.

    PubMed

    Baka, N; Metz, C T; Schultz, C J; van Geuns, R-J; Niessen, W J; van Walsum, T

    2014-05-01

    2D/3D registration of patient vasculature from preinterventional computed tomography angiography (CTA) to interventional X-ray angiography is of interest to improve guidance in percutaneous coronary interventions. In this paper we present a novel feature based 2D/3D registration framework, that is based on probabilistic point correspondences, and show its usefulness on aligning 3D coronary artery centerlines derived from CTA images with their 2D projection derived from interventional X-ray angiography. The registration framework is an extension of the Gaussian mixture model (GMM) based point-set registration to the 2D/3D setting, with a modified distance metric. We also propose a way to incorporate orientation in the registration, and show its added value for artery registration on patient datasets as well as in simulation experiments. The oriented GMM registration achieved a median accuracy of 1.06 mm, with a convergence rate of 81% for nonrigid vessel centerline registration on 12 patient datasets, using a statistical shape model. The method thereby outperformed the iterative closest point algorithm, the GMM registration without orientation, and two recently published methods on 2D/3D coronary artery registration. PMID:24770908

  17. The simulation of 3D mass models in 2D digital mammography and breast tomosynthesis

    SciTech Connect

    Shaheen, Eman De Keyzer, Frederik; Bosmans, Hilde; Ongeval, Chantal Van; Dance, David R.; Young, Kenneth C.

    2014-08-15

    Purpose: This work proposes a new method of building 3D breast mass models with different morphological shapes and describes the validation of the realism of their appearance after simulation into 2D digital mammograms and breast tomosynthesis images. Methods: Twenty-five contrast enhanced MRI breast lesions were collected and each mass was manually segmented in the three orthogonal views: sagittal, coronal, and transversal. The segmented models were combined, resampled to have isotropic voxel sizes, triangularly meshed, and scaled to different sizes. These masses were referred to as nonspiculated masses and were then used as nuclei onto which spicules were grown with an iterative branching algorithm forming a total of 30 spiculated masses. These 55 mass models were projected into 2D projection images to obtain mammograms after image processing and into tomographic sequences of projection images, which were then reconstructed to form 3D tomosynthesis datasets. The realism of the appearance of these mass models was assessed by five radiologists via receiver operating characteristic (ROC) analysis when compared to 54 real masses. All lesions were also given a breast imaging reporting and data system (BIRADS) score. The data sets of 2D mammography and tomosynthesis were read separately. The Kendall's coefficient of concordance was used for the interrater observer agreement assessment for the BIRADS scores per modality. Further paired analysis, using the Wilcoxon signed rank test, of the BIRADS assessment between 2D and tomosynthesis was separately performed for the real masses and for the simulated masses. Results: The area under the ROC curves, averaged over all observers, was 0.54 (95% confidence interval [0.50, 0.66]) for the 2D study, and 0.67 (95% confidence interval [0.55, 0.79]) for the tomosynthesis study. According to the BIRADS scores, the nonspiculated and the spiculated masses varied in their degrees of malignancy from normal (BIRADS 1) to highly

  18. Impact of high speed civil transports on stratospheric ozone: A 2-D model investigation

    SciTech Connect

    Kinnison, D.E.; Connell, P.S.

    1996-12-01

    This study investigates the effect on stratospheric ozone from a fleet of proposed High Speed Civil Transports (HSCTs). The new LLNL 2-D operator-split chemical-radiative-transport model of the troposphere and stratosphere is used for this HSCT investigation. This model is integrated in a diurnal manner, using an implicit numerical solver. Therefore, rate coefficients are not modified by any sort of diurnal average factor. This model also does not make any assumptions on lumping of chemical species into families. Comparisons to previous model-derived HSCT assessment of ozone change are made, both to the previous LLNL 2-D model and to other models from the international assessment modeling community. The sensitivity to the NO{sub x} emission index and sulfate surface area density is also explored.

  19. An enhanced CCRTM (E-CCRTM) damage imaging technique using a 2D areal scan for composite plates

    NASA Astrophysics Data System (ADS)

    He, Jiaze; Yuan, Fuh-Gwo

    2016-04-01

    A two-dimensional (2-D) non-contact areal scan system was developed to image and quantify impact damage in a composite plate using an enhanced zero-lag cross-correlation reverse-time migration (E-CCRTM) technique. The system comprises a single piezoelectric actuator mounted on the composite plate and a laser Doppler vibrometer (LDV) for scanning a region to capture the scattered wavefield in the vicinity of the PZT. The proposed damage imaging technique takes into account the amplitude, phase, geometric spreading, and all of the frequency content of the Lamb waves propagating in the plate; thus, the reflectivity coefficients of the delamination can be calculated and potentially related to damage severity. Comparisons are made in terms of damage imaging quality between 2-D areal scans and linear scans as well as between the proposed and existing imaging conditions. The experimental results show that the 2-D E-CCRTM performs robustly when imaging and quantifying impact damage in large-scale composites using a single PZT actuator with a nearby areal scan using LDV.

  20. Extension and application of the Preissmann slot model to 2D transient mixed flows

    NASA Astrophysics Data System (ADS)

    Maranzoni, Andrea; Dazzi, Susanna; Aureli, Francesca; Mignosa, Paolo

    2015-08-01

    This paper presents an extension of the Preissmann slot concept for the modeling of highly transient two-dimensional (2D) mixed flows. The classic conservative formulation of the 2D shallow water equations for free surface flows is adapted by assuming that two fictitious vertical slots, aligned along the two Cartesian plane directions and normally intersecting, are added on the ceiling of each integration element. Accordingly, transitions between free surface and pressurized flow can be handled in a natural and straightforward way by using the same set of governing equations. The opportunity of coupling free surface and pressurized flows is actually useful not only in one-dimensional (1D) problems concerning sewer systems but also for modeling 2D flooding phenomena in which the pressurization of bridges, culverts, or other crossing hydraulic structures can be expected. Numerical simulations are performed by using a shock-capturing MUSCL-Hancock finite volume scheme combined with the FORCE (First-Order Centred) solver for the evaluation of the numerical fluxes. The validation of the mathematical model is accomplished on the basis of both exact solutions of 1D discontinuous initial value problems and reference radial solutions of idealized test cases with cylindrical symmetry. Furthermore, the capability of the model to deal with practical field-scale applications is assessed by simulating the transit of a bore under an arch bridge. Numerical results show that the proposed model is suitable for the prediction of highly transient 2D mixed flows.

  1. Comparison of 3-D finite element model of ashlar masonry with 2-D numerical models of ashlar masonry

    NASA Astrophysics Data System (ADS)

    Beran, Pavel

    2016-06-01

    3-D state of stress in heterogeneous ashlar masonry can be also computed by several suitable chosen 2-D numerical models of ashlar masonry. The results obtained from 2-D numerical models well correspond to the results obtained from 3-D numerical model. The character of thermal stress is the same. While using 2-D models the computational time is reduced more than hundredfold and therefore this method could be used for computation of thermal stresses during long time periods with 10 000 of steps.

  2. 2D-Raman-THz spectroscopy: A sensitive test of polarizable water models

    NASA Astrophysics Data System (ADS)

    Hamm, Peter

    2014-11-01

    In a recent paper, the experimental 2D-Raman-THz response of liquid water at ambient conditions has been presented [J. Savolainen, S. Ahmed, and P. Hamm, Proc. Natl. Acad. Sci. U. S. A. 110, 20402 (2013)]. Here, all-atom molecular dynamics simulations are performed with the goal to reproduce the experimental results. To that end, the molecular response functions are calculated in a first step, and are then convoluted with the laser pulses in order to enable a direct comparison with the experimental results. The molecular dynamics simulation are performed with several different water models: TIP4P/2005, SWM4-NDP, and TL4P. As polarizability is essential to describe the 2D-Raman-THz response, the TIP4P/2005 water molecules are amended with either an isotropic or a anisotropic polarizability a posteriori after the molecular dynamics simulation. In contrast, SWM4-NDP and TL4P are intrinsically polarizable, and hence the 2D-Raman-THz response can be calculated in a self-consistent way, using the same force field as during the molecular dynamics simulation. It is found that the 2D-Raman-THz response depends extremely sensitively on details of the water model, and in particular on details of the description of polarizability. Despite the limited time resolution of the experiment, it could easily distinguish between various water models. Albeit not perfect, the overall best agreement with the experimental data is obtained for the TL4P water model.

  3. 2D-Raman-THz spectroscopy: A sensitive test of polarizable water models

    SciTech Connect

    Hamm, Peter

    2014-11-14

    In a recent paper, the experimental 2D-Raman-THz response of liquid water at ambient conditions has been presented [J. Savolainen, S. Ahmed, and P. Hamm, Proc. Natl. Acad. Sci. U. S. A. 110, 20402 (2013)]. Here, all-atom molecular dynamics simulations are performed with the goal to reproduce the experimental results. To that end, the molecular response functions are calculated in a first step, and are then convoluted with the laser pulses in order to enable a direct comparison with the experimental results. The molecular dynamics simulation are performed with several different water models: TIP4P/2005, SWM4-NDP, and TL4P. As polarizability is essential to describe the 2D-Raman-THz response, the TIP4P/2005 water molecules are amended with either an isotropic or a anisotropic polarizability a posteriori after the molecular dynamics simulation. In contrast, SWM4-NDP and TL4P are intrinsically polarizable, and hence the 2D-Raman-THz response can be calculated in a self-consistent way, using the same force field as during the molecular dynamics simulation. It is found that the 2D-Raman-THz response depends extremely sensitively on details of the water model, and in particular on details of the description of polarizability. Despite the limited time resolution of the experiment, it could easily distinguish between various water models. Albeit not perfect, the overall best agreement with the experimental data is obtained for the TL4P water model.

  4. Analysis of Korean Students' International Mobility by 2-D Model: Driving Force Factor and Directional Factor

    ERIC Educational Resources Information Center

    Park, Elisa L.

    2009-01-01

    The purpose of this study is to understand the dynamics of Korean students' international mobility to study abroad by using the 2-D Model. The first D, "the driving force factor," explains how and what components of the dissatisfaction with domestic higher education perceived by Korean students drives students' outward mobility to seek foreign…

  5. Parallelized CCHE2D flow model with CUDA Fortran on Graphics Process Units

    Technology Transfer Automated Retrieval System (TEKTRAN)

    This paper presents the CCHE2D implicit flow model parallelized using CUDA Fortran programming technique on Graphics Processing Units (GPUs). A parallelized implicit Alternating Direction Implicit (ADI) solver using Parallel Cyclic Reduction (PCR) algorithm on GPU is developed and tested. This solve...

  6. Field lines and magnetic surfaces in a two-component slab/2D model of interplanetary magnetic fluctuations

    NASA Technical Reports Server (NTRS)

    Matthaeus, W. H.; Pontius, D. H., Jr.; Gray, P. C.; Bieber, J. W.

    1995-01-01

    A two-component model for the spectrum of interplanetary magnetic fluctuations was proposed on the basis of ISEE observations, and has found an intriguing level of application in other solar wind studies. The model fluctuations consist of a fraction of 'slab' fluctuations, varying only in the direction parallel to the locally uniform mean magnetic field B(0) and a complement of 2D (two-dimensional) fluctuations that vary in the directions transverse to B(0). We have developed an spectral method computational algorithm for computing the magnetic flux surfaces (flux tubes) associated with the composite model, based upon a precise analogy with equations for ideal transport of a passive scalar in planar two dimensional geometry. Visualization of various composite models will be presented, including the 80 percent 2D/ 20 percent slab model with delta B/B(0) approximately equals 1 and a minus 5/3 spectral law, that is thought to approximately represent a snapshot of solar wind turbulence. Characteristically, the visualizations show that flux tubes, even when defined as regular on some plane, shred and disperse rapidly as they are viewed along the parallel direction. This diffusive process, which generalizes the standard picture of field line random walk, will be discussed in detail. Evidently, the traditional picture that flux tubes randomize like strands of spaghetti with a uniform tangle along the axial direction is in need of modification.

  7. On craton thinning/destruction: Insight from 2D thermal-mechanical numerical modeling

    NASA Astrophysics Data System (ADS)

    Liao, J.

    2014-12-01

    Although most cratons maintain stable, some exceptions are present, such as the North China craton, North Atlantic craton, and Wyoming craton, which have experienced dramatic lithospheric deformation/thinning. Mechanisms triggering cratonic thinning remains enigmatic [Lee et al., 2011]. Using a 2D thermo-mechanical coupled numerical model [Gerya and Yuen, 2007], we investigate two possible mechanisms: (1) stratification of cratonic lithospheric mantle, and (2) rheological weakening due to hydration.Lithospheric mantle stratification is a common feature in cratonic areas which has been demonstrated by geophysical and geochemical studies [Thybo and Perchuc, 1997; Griffin et al., 2004; Romanowicz, 2009; Rychert and Shearer, 2009; Yuan and Romanowicz, 2010]. The influence of lithospheric mantle stratification during craton evolution remains poorly understood. A rheologically weak layer representing hydrated and/or metasomatized composition is implemented in the lithospheric mantle. Our results show that the weak mantle layer changes the dynamics of lithospheric extension by enhancing the deformation of the overlying mantle and crust and inhibiting deformation of the underlying mantle [Liao et al., 2013; Liao and Gerya, 2014]. Modeling results are compared with North China and North Atlantic cratons. Our work indicates that although the presence of a weak layer may not be sufficient to initiate craton deformation, it enhances deformation by lowering the required extensional plate boundary force. Rheological weakening due to hydration is a possible mechanism triggering/enhancing craton deformation, especially for cratons jaxtaposing with a subduction, since water can release from a subducting slab. We investigate the influence of wet mantle flow laws [Hirth and Kohlstedt, 2003], in which a water parameter (i.e. constant water content) is involved. Our results show that wet dislocation alone does not accelerate cratonic deformation significantly. However, if wet diffusion

  8. 2-D elasticity solution of layered composite beams with viscoelastic interlayers

    NASA Astrophysics Data System (ADS)

    Wu, Peng; Zhou, Ding; Liu, Weiqing

    2016-02-01

    This paper focuses on the mechanical properties of layered composite beams with viscoelastic interlayers. The exact two-dimensional elasticity theory is used to represent the deformation of each beam layer. The viscoelastic interlayer is described by the Maxwell-Wiechert model through the quasi-elastic approximation, which greatly simplifies the analytical process. The stress function with a series of undetermined coefficients depending on the time variable is derived for each beam layer. No matter how many layers the beam includes, the total solution can be obtained rapidly and efficiently by using the recursive matrix technique. The present method can give the exact stress and deformation distributions in the beam, which cannot be predicted by the approximate theories such as the one-dimensional Euler-Bernoulli theory. The convergence of the solution is numerically verified. A comparison study indicates that the present results are in agreement with those obtained from the finite element method; however, they have obvious differences from the results based on the Euler-Bernoulli theory for thick beams. Finally, the variations of stresses and displacements with respect to time in a five-layer beam are discussed in detail.

  9. Molecular Dynamics implementation of BN2D or 'Mercedes Benz' water model

    NASA Astrophysics Data System (ADS)

    Scukins, Arturs; Bardik, Vitaliy; Pavlov, Evgen; Nerukh, Dmitry

    2015-05-01

    Two-dimensional 'Mercedes Benz' (MB) or BN2D water model (Naim, 1971) is implemented in Molecular Dynamics. It is known that the MB model can capture abnormal properties of real water (high heat capacity, minima of pressure and isothermal compressibility, negative thermal expansion coefficient) (Silverstein et al., 1998). In this work formulas for calculating the thermodynamic, structural and dynamic properties in microcanonical (NVE) and isothermal-isobaric (NPT) ensembles for the model from Molecular Dynamics simulation are derived and verified against known Monte Carlo results. The convergence of the thermodynamic properties and the system's numerical stability are investigated. The results qualitatively reproduce the peculiarities of real water making the model a visually convenient tool that also requires less computational resources, thus allowing simulations of large (hydrodynamic scale) molecular systems. We provide the open source code written in C/C++ for the BN2D water model implementation using Molecular Dynamics.

  10. Justification for a 2D versus 3D fingertip finite element model during static contact simulations.

    PubMed

    Harih, Gregor; Tada, Mitsunori; Dolšak, Bojan

    2016-10-01

    The biomechanical response of a human hand during contact with various products has not been investigated in details yet. It has been shown that excessive contact pressure on the soft tissue can result in discomfort, pain and also cumulative traumatic disorders. This manuscript explores the benefits and limitations of a simplified two-dimensional vs. an anatomically correct three-dimensional finite element model of a human fingertip. Most authors still use 2D FE fingertip models due to their simplicity and reduced computational costs. However we show that an anatomically correct 3D FE fingertip model can provide additional insight into the biomechanical behaviour. The use of 2D fingertip FE models is justified when observing peak contact pressure values as well as displacement during the contact for the given studied cross-section. On the other hand, an anatomically correct 3D FE fingertip model provides a contact pressure distribution, which reflects the fingertip's anatomy. PMID:26856769

  11. The 2dF Galaxy Redshift Survey: voids and hierarchical scaling models

    NASA Astrophysics Data System (ADS)

    Croton, Darren J.; Colless, Matthew; Gaztañaga, Enrique; Baugh, Carlton M.; Norberg, Peder; Baldry, I. K.; Bland-Hawthorn, J.; Bridges, T.; Cannon, R.; Cole, S.; Collins, C.; Couch, W.; Dalton, G.; de Propris, R.; Driver, S. P.; Efstathiou, G.; Ellis, R. S.; Frenk, C. S.; Glazebrook, K.; Jackson, C.; Lahav, O.; Lewis, I.; Lumsden, S.; Maddox, S.; Madgwick, D.; Peacock, J. A.; Peterson, B. A.; Sutherland, W.; Taylor, K.

    2004-08-01

    We measure the redshift-space reduced void probability function (VPF) for 2dFGRS volume-limited galaxy samples covering the absolute magnitude range MbJ-5log10h=-18 to -22. Theoretically, the VPF connects the distribution of voids to the moments of galaxy clustering of all orders, and can be used to discriminate clustering models in the weakly non-linear regime. The reduced VPF measured from the 2dFGRS is in excellent agreement with the paradigm of hierarchical scaling of the galaxy clustering moments. The accuracy of our measurement is such that we can rule out, at a very high significance, popular models for galaxy clustering, including the lognormal distribution. We demonstrate that the negative binomial model gives a very good approximation to the 2dFGRS data over a wide range of scales, out to at least 20 h-1 Mpc. Conversely, the reduced VPF for dark matter in a Λ cold dark matter (ΛCDM) universe does appear to be lognormal on small scales but deviates significantly beyond ~4 h-1 Mpc. We find little dependence of the 2dFGRS reduced VPF on galaxy luminosity. Our results hold independently in both the North and South Galactic Pole survey regions.

  12. 1D/2D Carbon Nanomaterial-Polymer Dielectric Composites with High Permittivity for Power Energy Storage Applications.

    PubMed

    Dang, Zhi-Min; Zheng, Ming-Sheng; Zha, Jun-Wei

    2016-04-01

    With the development of flexible electronic devices and large-scale energy storage technologies, functional polymer-matrix nanocomposites with high permittivity (high-k) are attracting more attention due to their ease of processing, flexibility, and low cost. The percolation effect is often used to explain the high-k characteristic of polymer composites when the conducting functional fillers are dispersed into polymers, which gives the polymer composite excellent flexibility due to the very low loading of fillers. Carbon nanotubes (CNTs) and graphene nanosheets (GNs), as one-dimensional (1D) and two-dimensional (2D) carbon nanomaterials respectively, have great potential for realizing flexible high-k dielectric nanocomposites. They are becoming more attractive for many fields, owing to their unique and excellent advantages. The progress in dielectric fields by using 1D/2D carbon nanomaterials as functional fillers in polymer composites is introduced, and the methods and mechanisms for improving dielectric properties, breakdown strength and energy storage density of their dielectric nanocomposites are examined. Achieving a uniform dispersion state of carbon nanomaterials and preventing the development of conductive networks in their polymer composites are the two main issues that still need to be solved in dielectric fields for power energy storage. Recent findings, current problems, and future perspectives are summarized. PMID:26865507

  13. Electrical resistivity tomography applied to a complex lava dome: 2D and 3D models comparison

    NASA Astrophysics Data System (ADS)

    Portal, Angélie; Fargier, Yannick; Lénat, Jean-François; Labazuy, Philippe

    2015-04-01

    The study of volcanic domes growth (e.g. St. Helens, Unzen, Montserrat) shows that it is often characterized by a succession of extrusion phases, dome explosions and collapse events. Lava dome eruptive activity may last from days to decades. Therefore, their internal structure, at the end of the eruption, is complex and includes massive extrusions and lava lobes, talus and pyroclastic deposits as well as hydrothermal alteration. The electrical resistivity tomography (ERT) method, initially developed for environmental and engineering exploration, is now commonly used for volcano structure imaging. Because a large range of resistivity values is often observed in volcanic environments, the method is well suited to study the internal structure of volcanic edifices. We performed an ERT survey on an 11ka years old trachytic lava dome, the Puy de Dôme volcano (French Massif Central). The analysis of a recent high resolution DEM (LiDAR 0.5 m), as well as other geophysical data, strongly suggest that the Puy de Dôme is a composite dome. 11 ERT profiles have been carried out, both at the scale of the entire dome (base diameter of ~2 km and height of 400 m) on the one hand, and at a smaller scale on the summit part on the other hand. Each profile is composed of 64 electrodes. Three different electrode spacing have been used depending on the study area (35 m for the entire dome, 10 m and 5 m for its summit part). Some profiles were performed with half-length roll-along acquisitions, in order to keep a good trade-off between depth of investigation and resolution. Both Wenner-alpha and Wenner-Schlumberger protocols were used. 2-D models of the electrical resistivity distribution were computed using RES2DINV software. In order to constrain inversion models interpretation, the depth of investigation (DOI) method was applied to those results. It aims to compute a sensitivity index on inversion results, illustrating how the data influence the model and constraining models

  14. Simplified 2D Bidomain Model of Whole Heart Electrical Activity and ECG Generation

    NASA Astrophysics Data System (ADS)

    Sovilj, Siniša; Magjarević, Ratko; Abed, Amr Al; Lovell, Nigel H.; Dokos, Socrates

    2014-06-01

    The aim of this study was the development of a geometrically simple and highly computationally-efficient two dimensional (2D) biophysical model of whole heart electrical activity, incorporating spontaneous activation of the sinoatrial node (SAN), the specialized conduction system, and realistic surface ECG morphology computed on the torso. The FitzHugh-Nagumo (FHN) equations were incorporated into a bidomain finite element model of cardiac electrical activity, which was comprised of a simplified geometry of the whole heart with the blood cavities, the lungs and the torso as an extracellular volume conductor. To model the ECG, we placed four electrodes on the surface of the torso to simulate three Einthoven leads VI, VII and VIII from the standard 12-lead system. The 2D model was able to reconstruct ECG morphology on the torso from action potentials generated at various regions of the heart, including the sinoatrial node, atria, atrioventricular node, His bundle, bundle branches, Purkinje fibers, and ventricles. Our 2D cardiac model offers a good compromise between computational load and model complexity, and can be used as a first step towards three dimensional (3D) ECG models with more complex, precise and accurate geometry of anatomical structures, to investigate the effect of various cardiac electrophysiological parameters on ECG morphology.

  15. CLFE2D: A generalized plane strain finite element program laminated composites subject to mechanical and hygrothermal loading

    NASA Technical Reports Server (NTRS)

    Buczek, M. B.; Gregory, M. A.; Herakovich, C. T.

    1983-01-01

    CLFE2D is a two dimensional generalized plane strain finite element code, using a linear, four node, general quadrilateral, isoparametric element. The program is developed to calculate the displacements, strains, stresses, and strain energy densities in a finite width composite laminate. CLFE2D offers any combination of the following load types: nodal displacements, nodal forces, uniform normal strain, or hygrothermal. The program allows the user to input one set of three dimensional orthotropic material properties. The user can then specify the angle of material principal orientation for each element in the mesh. Output includes displacements, stresses, strains and strain densities at points selected by the user. An option is also available to plot the underformed and deformed finite element meshes.

  16. Approximate analytic solutions to 3D unconfined groundwater flow within regional 2D models

    NASA Astrophysics Data System (ADS)

    Luther, K.; Haitjema, H. M.

    2000-04-01

    We present methods for finding approximate analytic solutions to three-dimensional (3D) unconfined steady state groundwater flow near partially penetrating and horizontal wells, and for combining those solutions with regional two-dimensional (2D) models. The 3D solutions use distributed singularities (analytic elements) to enforce boundary conditions on the phreatic surface and seepage faces at vertical wells, and to maintain fixed-head boundary conditions, obtained from the 2D model, at the perimeter of the 3D model. The approximate 3D solutions are analytic (continuous and differentiable) everywhere, including on the phreatic surface itself. While continuity of flow is satisfied exactly in the infinite 3D flow domain, water balance errors can occur across the phreatic surface.

  17. TRENT2D WG: a smart web infrastructure for debris-flow modelling and hazard assessment

    NASA Astrophysics Data System (ADS)

    Zorzi, Nadia; Rosatti, Giorgio; Zugliani, Daniel; Rizzi, Alessandro; Piffer, Stefano

    2016-04-01

    Mountain regions are naturally exposed to geomorphic flows, which involve large amounts of sediments and induce significant morphological modifications. The physical complexity of this class of phenomena represents a challenging issue for modelling, leading to elaborate theoretical frameworks and sophisticated numerical techniques. In general, geomorphic-flows models proved to be valid tools in hazard assessment and management. However, model complexity seems to represent one of the main obstacles to the diffusion of advanced modelling tools between practitioners and stakeholders, although the UE Flood Directive (2007/60/EC) requires risk management and assessment to be based on "best practices and best available technologies". Furthermore, several cutting-edge models are not particularly user-friendly and multiple stand-alone software are needed to pre- and post-process modelling data. For all these reasons, users often resort to quicker and rougher approaches, leading possibly to unreliable results. Therefore, some effort seems to be necessary to overcome these drawbacks, with the purpose of supporting and encouraging a widespread diffusion of the most reliable, although sophisticated, modelling tools. With this aim, this work presents TRENT2D WG, a new smart modelling solution for the state-of-the-art model TRENT2D (Armanini et al., 2009, Rosatti and Begnudelli, 2013), which simulates debris flows and hyperconcentrated flows adopting a two-phase description over a mobile bed. TRENT2D WG is a web infrastructure joining advantages offered by the software-delivering model SaaS (Software as a Service) and by WebGIS technology and hosting a complete and user-friendly working environment for modelling. In order to develop TRENT2D WG, the model TRENT2D was converted into a service and exposed on a cloud server, transferring computational burdens from the user hardware to a high-performing server and reducing computational time. Then, the system was equipped with an

  18. A Convective Vorticity Vector Associated With Tropical Convection: A 2D Cloud-Resolving Modeling Study

    NASA Technical Reports Server (NTRS)

    Gao, Shou-Ting; Ping, Fan; Li, Xiao-Fan; Tao, Wei-Kuo

    2004-01-01

    Although dry/moist potential vorticity is a useful physical quantity for meteorological analysis, it cannot be applied to the analysis of 2D simulations. A convective vorticity vector (CVV) is introduced in this study to analyze 2D cloud-resolving simulation data associated with 2D tropical convection. The cloud model is forced by the vertical velocity, zonal wind, horizontal advection, and sea surface temperature obtained from the TOGA COARE, and is integrated for a selected 10-day period. The CVV has zonal and vertical components in the 2D x-z frame. Analysis of zonally-averaged and mass-integrated quantities shows that the correlation coefficient between the vertical component of the CVV and the sum of the cloud hydrometeor mixing ratios is 0.81, whereas the correlation coefficient between the zonal component and the sum of the mixing ratios is only 0.18. This indicates that the vertical component of the CVV is closely associated with tropical convection. The tendency equation for the vertical component of the CVV is derived and the zonally-averaged and mass-integrated tendency budgets are analyzed. The tendency of the vertical component of the CVV is determined by the interaction between the vorticity and the zonal gradient of cloud heating. The results demonstrate that the vertical component of the CVV is a cloud-linked parameter and can be used to study tropical convection.

  19. Parameterising root system growth models using 2D neutron radiography images

    NASA Astrophysics Data System (ADS)

    Schnepf, Andrea; Felderer, Bernd; Vontobel, Peter; Leitner, Daniel

    2013-04-01

    Root architecture is a key factor for plant acquisition of water and nutrients from soil. In particular in view of a second green revolution where the below ground parts of agricultural crops are important, it is essential to characterise and quantify root architecture and its effect on plant resource acquisition. Mathematical models can help to understand the processes occurring in the soil-plant system, they can be used to quantify the effect of root and rhizosphere traits on resource acquisition and the response to environmental conditions. In order to do so, root architectural models are coupled with a model of water and solute transport in soil. However, dynamic root architectural models are difficult to parameterise. Novel imaging techniques such as x-ray computed tomography, neutron radiography and magnetic resonance imaging enable the in situ visualisation of plant root systems. Therefore, these images facilitate the parameterisation of dynamic root architecture models. These imaging techniques are capable of producing 3D or 2D images. Moreover, 2D images are also available in the form of hand drawings or from images of standard cameras. While full 3D imaging tools are still limited in resolutions, 2D techniques are a more accurate and less expensive option for observing roots in their environment. However, analysis of 2D images has additional difficulties compared to the 3D case, because of overlapping roots. We present a novel algorithm for the parameterisation of root system growth models based on 2D images of root system. The algorithm analyses dynamic image data. These are a series of 2D images of the root system at different points in time. Image data has already been adjusted for missing links and artefacts and segmentation was performed by applying a matched filter response. From this time series of binary 2D images, we parameterise the dynamic root architecture model in the following way: First, a morphological skeleton is derived from the binary

  20. Validation of a 2-D semi-coupled numerical model for fluid-structure-seabed interaction

    NASA Astrophysics Data System (ADS)

    Ye, Jianhong; Jeng, Dongsheng; Wang, Ren; Zhu, Changqi

    2013-10-01

    A 2-D semi-coupled model PORO-WSSI 2D (also be referred as FSSI-CAS 2D) for the Fluid-Structure-Seabed Interaction (FSSI) has been developed by employing RANS equations for wave motion in fluid domain, VARANS equations for porous flow in porous structures; and taking the dynamic Biot's equations (known as "u - p" approximation) for soil as the governing equations. The finite difference two-step projection method and the forward time difference method are adopted to solve the RANS, VARANS equations; and the finite element method is adopted to solve the "u - p" approximation. A data exchange port is developed to couple the RANS, VARANS equations and the dynamic Biot's equations together. The analytical solution proposed by Hsu and Jeng (1994) and some experiments conducted in wave flume or geotechnical centrifuge in which various waves involved are used to validate the developed semi-coupled numerical model. The sandy bed involved in these experiments is poro-elastic or poro-elastoplastic. The inclusion of the interaction between fluid, marine structures and poro-elastoplastic seabed foundation is a special point and highlight in this paper, which is essentially different with other previous coupled models The excellent agreement between the numerical results and the experiment data indicates that the developed coupled model is highly reliablefor the FSSI problem.

  1. 2D density model of the Chinese continental lithosphere along a NW-SE transect

    NASA Astrophysics Data System (ADS)

    Šimonová, Barbora; Bielik, Miroslav; Dérerová, Jana

    2015-06-01

    This paper presents a 2D density model along a transect from NW to SE China. The model was first constructed by the transformation of seismic velocity to density, revealed by previous deep seismic soundings (DSS) investigations in China. Then, the 2D density model was updated using the GM-SYS software by fitting the computed to the observed gravity data. Based on the density distribution of anomalous layers we divided the Chinese continental crust along the transect into three regions: north-western, central and south-eastern. The first one includes the Junggar Basin, Tianshan and Tarim Basin. The second part consists of the Qilian Orogen, the Qaidam Basin and the Songpan Ganzi Basin. The third region is represented by the Yangtze and the Cathaysia blocks. The low velocity body (vp =5.2 - 6.2 km/s) at the junction of the North-western and Central parts at a depth between 21 - 31 km, which was discovered out by DSS, was also confirmed by our 2D density modelling.

  2. A Neural-FEM tool for the 2-D magnetic hysteresis modeling

    NASA Astrophysics Data System (ADS)

    Cardelli, E.; Faba, A.; Laudani, A.; Lozito, G. M.; Riganti Fulginei, F.; Salvini, A.

    2016-04-01

    The aim of this work is to present a new tool for the analysis of magnetic field problems considering 2-D magnetic hysteresis. In particular, this tool makes use of the Finite Element Method to solve the magnetic field problem in real device, and fruitfully exploits a neural network (NN) for the modeling of 2-D magnetic hysteresis of materials. The NS has as input the magnetic inductions components B at the k-th simulation step and returns as output the corresponding values of the magnetic field H corresponding to the input pattern. It is trained by vector measurements performed on the magnetic material to be modeled. This input/output scheme is directly implemented in a FEM code employing the magnetic potential vector A formulation. Validations through measurements on a real device have been performed.

  3. Effects of preform architecture on modulus and strength of 2-D triaxially braided textile composites

    NASA Technical Reports Server (NTRS)

    Masters, John E.; Naik, Rajiv; Minguet, Pierre J.

    1995-01-01

    Laminates formed using braided fibrous preforms have been extensively investigated during the course of the past few years as alternatives to unidirectional prepreg tape systems. This paper focused on one aspect of that work. It defined the role of the fibrous preform architecture in controlling a laminate's mechanical properties. The presentation was divided into four sections as the outline listed above illustrates. The presentation began with a brief introduction which defined the objectives of the study and detailed the materials studied. This was followed by a review of empirical test results. The materials' moduli and strengths were measured in both tension and compression. Their shear moduli were also experimentally determined. The review of the empirical data comprised the bulk of the presentation. A comparison of the experimental data to results predicted analytically was then presented. The presentation concluded with a few summary remarks. The specimens studied in this investigation featured 2-D triaxially braided AS4 graphite fiber preforms impregnated with Shell 1895 epoxy resin.

  4. TMRPres2D: high quality visual representation of transmembrane protein models.

    PubMed

    Spyropoulos, Ioannis C; Liakopoulos, Theodore D; Bagos, Pantelis G; Hamodrakas, Stavros J

    2004-11-22

    The 'TransMembrane protein Re-Presentation in 2-Dimensions' (TMRPres2D) tool, automates the creation of uniform, two-dimensional, high analysis graphical images/models of alpha-helical or beta-barrel transmembrane proteins. Protein sequence data and structural information may be acquired from public protein knowledge bases, emanate from prediction algorithms, or even be defined by the user. Several important biological and physical sequence attributes can be embedded in the graphical representation. PMID:15201184

  5. Effects of Nesting on Compression-Loaded 2-D Woven Textile Composites

    NASA Technical Reports Server (NTRS)

    Adams, Daniel OHare; Breiling, Kurtis B.; Verhulst, Mark A.

    1995-01-01

    Layer nesting was investigated in five harness satin weave textile composite laminates under static compression loading. Two carbon/epoxy material systems, AS4/3501-6 and IM7/8551-7A were considered. Laminates were fabricated with three idealized nesting cases: stacked, split-span and diagonal. Similar compression strength reductions due to the effects of idealized nesting were identified for each material. The diagonal nesting geometry produced the largest reduction in static strength when compared to the compression strength of a conventional textile composite. All three nesting cases produced reductions in strength and ultimate strain due to the effects of idealized nesting. Finite element results showed consistent strength reduction trends for the idealized nesting cases, however the magnitudes of compressive strengths were overpredicted.

  6. A simple 2-D inundation model for incorporating flood damage in urban drainage planning

    NASA Astrophysics Data System (ADS)

    Pathirana, A.; Tsegaye, S.; Gersonius, B.; Vairavamoorthy, K.

    2008-11-01

    In this paper a new inundation model code is developed and coupled with Storm Water Management Model, SWMM, to relate spatial information associated with urban drainage systems as criteria for planning of storm water drainage networks. The prime objective is to achive a model code that is simple and fast enough to be consistently be used in planning stages of urban drainage projects. The formulation for the two-dimensional (2-D) surface flow model algorithms is based on the Navier Stokes equation in two dimensions. An Alternating Direction Implicit (ADI) finite difference numerical scheme is applied to solve the governing equations. This numerical scheme is used to express the partial differential equations with time steps split into two halves. The model algorithm is written using C++ computer programming language. This 2-D surface flow model is then coupled with SWMM for simulation of both pipe flow component and surcharge induced inundation in urban areas. In addition, a damage calculation block is integrated within the inundation model code. The coupled model is shown to be capable of dealing with various flow conditions, as well as being able to simulate wetting and drying processes that will occur as the flood flows over an urban area. It has been applied under idealized and semi-hypothetical cases to determine detailed inundation zones, depths and velocities due to surcharged water on overland surface.

  7. 2D Cross Sectional Analysis and Associated Electrochemistry of Composite Electrodes Containing Dispersed Agglomerates of Nanocrystalline Magnetite, Fe₃O₄.

    PubMed

    Bock, David C; Kirshenbaum, Kevin C; Wang, Jiajun; Zhang, Wei; Wang, Feng; Wang, Jun; Marschilok, Amy C; Takeuchi, Kenneth J; Takeuchi, Esther S

    2015-06-24

    When electroactive nanomaterials are fully incorporated into an electrode structure, characterization of the crystallite sizes, agglomerate sizes, and dispersion of the electroactive materials can lend insight into the complex electrochemistry associated with composite electrodes. In this study, composite magnetite electrodes were sectioned using ultramicrotome techniques, which facilitated the direct observation of crystallites and agglomerates of magnetite (Fe3O4) as well as their dispersal patterns in large representative sections of electrode, via 2D cross sectional analysis by Transmission Electron Microscopy (TEM). Further, the electrochemistry of these electrodes were recorded, and Transmission X-ray Microscopy (TXM) was used to determine the distribution of oxidation states of the reduced magnetite. Unexpectedly, while two crystallite sizes of magnetite were employed in the production of the composite electrodes, the magnetite agglomerate sizes and degrees of dispersion in the two composite electrodes were similar to each other. This observation illustrates the necessity for careful characterization of composite electrodes, in order to understand the effects of crystallite size, agglomerate size, and level of dispersion on electrochemistry. PMID:26024206

  8. Comparison between 2D and 3D Numerical Modelling of a hot forging simulative test

    SciTech Connect

    Croin, M.; Ghiotti, A.; Bruschi, S.

    2007-04-07

    The paper presents the comparative analysis between 2D and 3D modelling of a simulative experiment, performed in laboratory environment, in which operating conditions approximate hot forging of a turbine aerofoil section. The plane strain deformation was chosen as an ideal case to analyze the process because of the thickness variations in the final section and the consequent distributions of contact pressure and sliding velocity at the interface that are closed to the conditions of the real industrial process. In order to compare the performances of 2D and 3D approaches, two different analyses were performed and compared with the experiments in terms of loads and temperatures peaks at the interface between the dies and the workpiece.

  9. Nested 1D-2D approach for urban surface flood modeling

    NASA Astrophysics Data System (ADS)

    Murla, Damian; Willems, Patrick

    2015-04-01

    Floods in urban areas as a consequence of sewer capacity exceedance receive increased attention because of trends in urbanization (increased population density and impermeability of the surface) and climate change. Despite the strong recent developments in numerical modeling of water systems, urban surface flood modeling is still a major challenge. Whereas very advanced and accurate flood modeling systems are in place and operation by many river authorities in support of flood management along rivers, this is not yet the case in urban water management. Reasons include the small scale of the urban inundation processes, the need to have very high resolution topographical information available, and the huge computational demands. Urban drainage related inundation modeling requires a 1D full hydrodynamic model of the sewer network to be coupled with a 2D surface flood model. To reduce the computational times, 0D (flood cones), 1D/quasi-2D surface flood modeling approaches have been developed and applied in some case studies. In this research, a nested 1D/2D hydraulic model has been developed for an urban catchment at the city of Gent (Belgium), linking the underground sewer (minor system) with the overland surface (major system). For the overland surface flood modelling, comparison was made of 0D, 1D/quasi-2D and full 2D approaches. The approaches are advanced by considering nested 1D-2D approaches, including infiltration in the green city areas, and allowing the effects of surface storm water storage to be simulated. An optimal nested combination of three different mesh resolutions was identified; based on a compromise between precision and simulation time for further real-time flood forecasting, warning and control applications. Main streets as mesh zones together with buildings as void regions constitute one of these mesh resolution (3.75m2 - 15m2); they have been included since they channel most of the flood water from the manholes and they improve the accuracy of

  10. Mechanical Modelling of Pultrusion Process: 2D and 3D Numerical Approaches

    NASA Astrophysics Data System (ADS)

    Baran, Ismet; Hattel, Jesper H.; Akkerman, Remko; Tutum, Cem C.

    2015-02-01

    The process induced variations such as residual stresses and distortions are a critical issue in pultrusion, since they affect the structural behavior as well as the mechanical properties and geometrical precision of the final product. In order to capture and investigate these variations, a mechanical analysis should be performed. In the present work, the two dimensional (2D) quasi-static plane strain mechanical model for the pultrusion of a thick square profile developed by the authors is further improved using generalized plane strain elements. In addition to that, a more advanced 3D thermo-chemical-mechanical analysis is carried out using 3D quadratic elements which is a novel application for the numerical modelling of the pultrusion process. It is found that the 2D mechanical models give relatively reasonable and accurate stress and displacement evolutions in the transverse direction as compared to the 3D model. Moreover, the generalized plane strain model predicts the longitudinal process induced stresses more similar to the ones calculated in the 3D model as compared with the plane strain model.

  11. Brane brick models, toric Calabi-Yau 4-folds and 2d (0,2) quivers

    NASA Astrophysics Data System (ADS)

    Franco, Sebastián; Lee, Sangmin; Seong, Rak-Kyeong

    2016-02-01

    We introduce brane brick models, a novel type of Type IIA brane configurations consisting of D4-branes ending on an NS5-brane. Brane brick models are T-dual to D1-branes over singular toric Calabi-Yau 4-folds. They fully encode the infinite class of 2 d (generically) {N}=(0,2) gauge theories on the worldvolume of the D1-branes and streamline their connection to the probed geometries. For this purpose, we also introduce new combinatorial procedures for deriving the Calabi-Yau associated to a given gauge theory and vice versa.

  12. Canonical vs. micro-canonical sampling methods in a 2D Ising model

    SciTech Connect

    Kepner, J.

    1990-12-01

    Canonical and micro-canonical Monte Carlo algorithms were implemented on a 2D Ising model. Expressions for the internal energy, U, inverse temperature, Z, and specific heat, C, are given. These quantities were calculated over a range of temperature, lattice sizes, and time steps. Both algorithms accurately simulate the Ising model. To obtain greater than three decimal accuracy from the micro-canonical method requires that the more complicated expression for Z be used. The overall difference between the algorithms is small. The physics of the problem under study should be the deciding factor in determining which algorithm to use. 13 refs., 6 figs., 2 tabs.

  13. Complex zeros of the 2 d Ising model on dynamical random lattices

    NASA Astrophysics Data System (ADS)

    Ambjørn, J.; Anagnostopoulos, K. N.; Magnea, U.

    1998-04-01

    We study the zeros in the complex plane of the partition function for the Ising model coupled to 2 d quantum gravity for complex magnetic field and for complex temperature. We compute the zeros by using the exact solution coming from a two matrix model and by Monte Carlo simulations of Ising spins on dynamical triangulations. We present evidence that the zeros form simple one-dimensional patterns in the complex plane, and that the critical behaviour of the system is governed by the scaling of the distribution of singularities near the critical point.

  14. 2D-3D Registration of CT Vertebra Volume to Fluoroscopy Projection: A Calibration Model Assessment

    NASA Astrophysics Data System (ADS)

    Bifulco, P.; Cesarelli, M.; Allen, R.; Romano, M.; Fratini, A.; Pasquariello, G.

    2009-12-01

    This study extends a previous research concerning intervertebral motion registration by means of 2D dynamic fluoroscopy to obtain a more comprehensive 3D description of vertebral kinematics. The problem of estimating the 3D rigid pose of a CT volume of a vertebra from its 2D X-ray fluoroscopy projection is addressed. 2D-3D registration is obtained maximising a measure of similarity between Digitally Reconstructed Radiographs (obtained from the CT volume) and real fluoroscopic projection. X-ray energy correction was performed. To assess the method a calibration model was realised a sheep dry vertebra was rigidly fixed to a frame of reference including metallic markers. Accurate measurement of 3D orientation was obtained via single-camera calibration of the markers and held as true 3D vertebra position; then, vertebra 3D pose was estimated and results compared. Error analysis revealed accuracy of the order of 0.1 degree for the rotation angles of about 1 mm for displacements parallel to the fluoroscopic plane, and of order of 10 mm for the orthogonal displacement.

  15. Effect of open hole on tensile failure properties of 2D triaxial braided textile composites and tape equivalents

    NASA Technical Reports Server (NTRS)

    Norman, Timothy L.; Anglin, Colin; Gaskin, David; Patrick, Mike

    1995-01-01

    The unnotched and notched (open hole) tensile strength and failure mechanisms of two-dimensional (2D) triaxial braided composites were examined. The effect of notch size and notch position were investigated. Damage initiation and propagation in notched and unnotched coupons were also examined. Theory developed to predict the normal stress distribution near an open hole and failure for tape laminated composites was evaluated for its applicability to triaxial braided textile composite materials. Four fiber architectures were considered with different combinations of braid angle, longitudinal and braider yam size, and percentage of longitudinal yarns. Tape laminates equivalent to textile composites were also constructed for comparison. Unnotched tape equivalents were stronger than braided textiles but exhibited greater notch sensitivity. Notched textiles and tape equivalents have roughly the same strength at large notch sizes. Two common damage mechanisms were found: braider yams cracking and near notch longitudinal yarn splitting. Cracking was found to initiate in braider yarns in unnotched and notched coupons, and propagate in the direction of the braider yarns until failure. Longitudinal yarn splitting occurred in three of four architectures that were longitudinally fiber dominated. Damage initiation stress decreased with increasing braid angle. No significant differences in prediction of near notch stress between measured and predicted stress were weak for textiles with large braid angle. Notch strength could not be predicted using existing anisotropic theory for braided textiles due to their insensitivity to notch.

  16. A 2-D semi-analytical model of double-gate tunnel field-effect transistor

    NASA Astrophysics Data System (ADS)

    Huifang, Xu; Yuehua, Dai; Ning, Li; Jianbin, Xu

    2015-05-01

    A 2-D semi-analytical model of double gate (DG) tunneling field-effect transistor (TFET) is proposed. By aid of introducing two rectangular sources located in the gate dielectric layer and the channel, the 2-D Poisson equation is solved by using a semi-analytical method combined with an eigenfunction expansion method. The expression of the surface potential is obtained, which is a special function for the infinite series expressions. The influence of the mobile charges on the potential profile is taken into account in the proposed model. On the basis of the potential profile, the shortest tunneling length and the average electrical field can be derived, and the drain current is then constructed by using Kane's model. In particular, the changes of the tunneling parameters Ak and Bk influenced by the drain—source voltage are also incorporated in the predicted model. The proposed model shows a good agreement with TCAD simulation results under different drain—source voltages, silicon film thicknesses, gate dielectric layer thicknesses, and gate dielectric layer constants. Therefore, it is useful to optimize the DG TFET and this provides a physical insight for circuit level design. Project supported by the National Natural Science Foundation of China (No. 61376106) and the Graduate Innovation Fund of Anhui University.

  17. Gender and ethnicity specific generic elastic models from a single 2D image for novel 2D pose face synthesis and recognition.

    PubMed

    Heo, Jingu; Savvides, Marios

    2012-12-01

    In this paper, we propose a novel method for generating a realistic 3D human face from a single 2D face image for the purpose of synthesizing new 2D face images at arbitrary poses using gender and ethnicity specific models. We employ the Generic Elastic Model (GEM) approach, which elastically deforms a generic 3D depth-map based on the sparse observations of an input face image in order to estimate the depth of the face image. Particularly, we show that Gender and Ethnicity specific GEMs (GE-GEMs) can approximate the 3D shape of the input face image more accurately, achieving a better generalization of 3D face modeling and reconstruction compared to the original GEM approach. We qualitatively validate our method using publicly available databases by showing each reconstructed 3D shape generated from a single image and new synthesized poses of the same person at arbitrary angles. For quantitative comparisons, we compare our synthesized results against 3D scanned data and also perform face recognition using synthesized images generated from a single enrollment frontal image. We obtain promising results for handling pose and expression changes based on the proposed method. PMID:22201062

  18. Microstructure and mechanical properties of 2D SiC-reinforced carbon matrix composites

    SciTech Connect

    Bourrat, X.; Evans, A.G.

    1995-10-01

    Fiber-reinforced composites consisting of combinations of SiC and C have been explored as materials having thermostructural characteristics suitable for elevated temperature applications. There are two opposing strategies. One uses C fibers and a matrix of SiC, usually produced by CVI. The other uses SiC fibers and a matrix consisting of C; those are designated CERACARB{trademark} and were studied in this paper. Among different grades of CERACARB{trademark} examined in this work, one has a high tensile strength and good notch resistance. This material has a relatively low interface sliding resistance, {tau} (10--15 MPa) and satisfies global load sharing (GLS) criteria for the ultimate tensile strength (UTS). The present designations for this material is LT. The second material, designated HT, has lower tensile strength and greater notch sensitivity. This material has a relatively large interface sliding resistance ({tau} {approx} 90MPa). This UTS loss is related to the addition of oxidation inhibitors (B{sub 4}C) as filler in the matrix. This work was focused on the structural investigation of this effect and the addition of a pyrocarbon interface to recover the right mechanical behavior (material LTC).

  19. An Integrative Model of Excitation Driven Fluid Flow in a 2D Uterine Channel

    NASA Astrophysics Data System (ADS)

    Maggio, Charles; Fauci, Lisa; Chrispell, John

    2009-11-01

    We present a model of intra-uterine fluid flow in a sagittal cross-section of the uterus by inducing peristalsis in a 2D channel. This is an integrative multiscale computational model that takes as input fluid viscosity, passive tissue properties of the uterine channel and a prescribed wave of membrane depolarization. This voltage pulse is coupled to a model of calcium dynamics inside a uterine smooth muscle cell, which in turn drives a kinetic model of myosin phosphorylation governing contractile muscle forces. Using the immersed boundary method, these muscle forces are communicated to a fluid domain to simulate the contractions which occur in a human uterus. An analysis of the effects of model parameters on the flow properties and emergent geometry of the peristaltic channel will be presented.

  20. Global regularity for the 2D Oldroyd-B model in the corotational case

    NASA Astrophysics Data System (ADS)

    Ye, Zhuan; Xu, Xiaojing

    2016-09-01

    This paper is dedicated to the Oldroyd-B model with fractional dissipation $(-\\Delta)^{\\alpha}\\tau$ for any $\\alpha>0$. We establish the global smooth solutions to the Oldroyd-B model in the corotational case with arbitrarily small fractional powers of the Laplacian in two spatial dimensions. The methods described here are quite different from the tedious iterative approach used in recent paper \\cite{XY}. Moreover, in the Appendix we provide some a priori estimates to the Oldroyd-B model in the critical case which may be useful and of interest for future improvement. Finally, the global regularity to to the Oldroyd-B model in the corotational case with $-\\Delta u$ replaced by $(-\\Delta)^{\\gamma}u$ for $\\gamma>1$ are also collected in the Appendix. Therefore our result is more closer to the resolution of the well-known global regularity issue on the critical 2D Oldroyd-B model.

  1. A solidification constitutive model for NIKE2D and NIKE3D

    SciTech Connect

    Raboin, P.J.

    1994-03-17

    This memo updates the current status of a solidification material model development which has been underway for more than a year. Significant modeling goals such as predicting cut-off stresses, thermo-elasto-plasticity, strain rate dependent plasticity and dynamic recovery have been completed. The model is called SOLMAT for solidification material model, and while developed for NIKE2D, it has already been implemented in NIKE3D and NIT03D by B. Maker. This memo details the future development strategy of SOLMAT including liquid and solid constitutive improvements, coupling of deviatoric and dilatational deformation and a plan to switch between constitutive theories. It explains some of the difficulties associated solidification modeling and proposes two experiments to measure properties for using SOLMAT. Due to the sensitive nature of these plans in relation to programmatic and CRADA concerns, this memo should be treated as confidential document.

  2. A velocity-dependent anomalous radial transport model for (2-D, 2-V) kinetic transport codes

    NASA Astrophysics Data System (ADS)

    Bodi, Kowsik; Krasheninnikov, Sergei; Cohen, Ron; Rognlien, Tom

    2008-11-01

    Plasma turbulence constitutes a significant part of radial plasma transport in magnetically confined plasmas. This turbulent transport is modeled in the form of anomalous convection and diffusion coefficients in fluid transport codes. There is a need to model the same in continuum kinetic edge codes [such as the (2-D, 2-V) transport version of TEMPEST, NEO, and the code being developed by the Edge Simulation Laboratory] with non-Maxwellian distributions. We present an anomalous transport model with velocity-dependent convection and diffusion coefficients leading to a diagonal transport matrix similar to that used in contemporary fluid transport models (e.g., UEDGE). Also presented are results of simulations corresponding to radial transport due to long-wavelength ExB turbulence using a velocity-independent diffusion coefficient. A BGK collision model is used to enable comparison with fluid transport codes.

  3. Exact solution of an anisotropic 2D random walk model with strong memory correlations

    NASA Astrophysics Data System (ADS)

    Cressoni, J. C.; Viswanathan, G. M.; da Silva, M. A. A.

    2013-12-01

    Over the last decade, there has been progress in understanding one-dimensional non-Markovian processes via analytic, sometimes exact, solutions. The extension of these ideas and methods to two and higher dimensions is challenging. We report the first exactly solvable two-dimensional (2D) non-Markovian random walk model belonging to the family of the elephant random walk model. In contrast to Lévy walks or fractional Brownian motion, such models incorporate memory effects by keeping an explicit history of the random walk trajectory. We study a memory driven 2D random walk with correlated memory and stops, i.e. pauses in motion. The model has an inherent anisotropy with consequences for its diffusive properties, thereby mixing the dominant regime along one dimension with a subdiffusive walk along a perpendicular dimension. The anomalous diffusion regimes are fully characterized by an exact determination of the Hurst exponent. We discuss the remarkably rich phase diagram, as well as several possible combinations of the independent walks in both directions. The relationship between the exponents of the first and second moments is also unveiled.

  4. Evaluation of Hydrus-2D model for solute distribution in subsurface drip

    NASA Astrophysics Data System (ADS)

    Souza, Claudinei; Bizari, Douglas; Grecco, Katarina

    2015-04-01

    The competition for water use between agriculture, industry and population has become intense over the years, requiring a rational use of this resource for food production. The subsurface drip irrigation can help producers with the optimization of operating parameters such as frequency and duration of irrigation, flow, spacing and depth of the dripper installation. This information can be obtained by numerical simulations using mathematical models, thus the aim of this study was to evaluate the HYDRUS-2D model from experimental data to predict the size of the wet bulbs generated by emitters of different application rates (1.0 and 1.6 L h-1). The results showed that horizontal displacement (bulb diameter) remained the largest in all the bulbs, observed both in experimental trials and estimated by the model and the correlation between them was high, above 0.90 to below 16% error. We conclude that the HYDRUS-2D model can be used to estimate the dimensions of the wet bulb getting new information on the sizing of the irrigation system.

  5. Molecular-dynamics of a 2D Model of the Shape Memory Effect

    NASA Astrophysics Data System (ADS)

    Kastner, Oliver

    2006-08-01

    This work investigates the thermodynamic properties of a qualitative atomistic model for austenite martensite transitions. The model, still in 2D, employs Lennard-Jones potentials for the determination of the atomic interactions. By use of two atom species it is possible to identify three stable lattice structures in 2D, interpreted as austenite and two variants of martensite. The model is described in the first part of the work [6] in detail. The present work studies the thermodynamic properties of the model concerning a small, 2-dimensional test assembly consisting of 41 atoms. The phase stability is investigated by exploitation of the condition of minimal free energy. The free energy is calculated from the thermal equation of state, which is measured in numerical tensile tests. In the second part of this work a chain of eleven 41-atom assemblies is investigated. The chain is interpreted as an idealized larger body, where the individual crystallites represent crystallographic layers allowing for the creation of micro structure. By use of tensile tests at various temperature conditions we sketch how such chain may exhibit quasi-plasticity, pseudo-elasticity and the shape memory effect.

  6. Momentum Transport: 2D and 3D Cloud Resolving Model Simulations

    NASA Technical Reports Server (NTRS)

    Tao, Wei-Kuo

    2001-01-01

    The major objective of this study is to investigate the momentum budgets associated with several convective systems that developed during the TOGA COARE IOP (west Pacific warm pool region) and GATE (east Atlantic region). The tool for this study is the improved Goddard Cumulas Ensemble (GCE) model which includes a 3-class ice-phase microphysical scheme, explicit cloud radiative interactive processes and air-sea interactive surface processes. The model domain contains 256 x 256 grid points (with 2 km resolution) in the horizontal and 38 grid points (to a depth of 22 km) in the vertical. The 2D domain has 1024 grid points. The simulations were performed over a 7-day time period (December 19-26, 1992, for TOGA COARE and September 1-7, 1994 for GATE). Cyclic literal boundary conditions are required for this type of long-term integration. Two well organized squall systems (TOGA, COARE February 22, 1993, and GATE September 12, 1994) were also simulated using the 3D GCE model. Only 9 h simulations were required to cover the life time of the squall systems. the lateral boundary conditions were open for these two squall systems simulations. the following will be examined: (1) the momentum budgets in the convective and stratiform regions, (2) the relationship between momentum transport and cloud organization (i.e., well organized squall lines versus less organized convective), (3) the differences and similarities in momentum transport between 2D and 3D simulated convective systems, and (4) the differences and similarities in momentum budgets between cloud systems simulated with open and cyclic lateral boundary conditions. Preliminary results indicate that there are only small differences between 2D and 3D simulated momentum budgets. Major differences occur, however, between momentum budgets associated with squall systems simulated using different lateral boundary conditions.

  7. A quasi 2D semianalytical model for the potential profile in hetero and homojunction tunnel FETs

    NASA Astrophysics Data System (ADS)

    Villani, F.; Gnani, E.; Gnudi, A.; Reggiani, S.; Baccarani, G.

    2015-11-01

    A quasi 2D semianalytical model for the potential profile in hetero and homojunction tunnel FETs is developed and compared with full-quantum simulation results. It will be shown that the pure analytical solution perfectly matches results at high VDS. However, a coupling with the numerical solution of the 1D Poisson equation in the radial direction is necessary at low VDS, in order to properly account for the charge density in equilibrium with the drain contact. With such an approach we are able to correctly predict the potential profile for both the linear and saturation regimes.

  8. Vector chiral phases in the frustrated 2D XY model and quantum spin chains.

    PubMed

    Schenck, H; Pokrovsky, V L; Nattermann, T

    2014-04-18

    The phase diagram of the frustrated 2D classical and 1D quantum XY models is calculated analytically. Four transitions are found: the vortex unbinding transitions triggered by strong fluctuations occur above and below the chiral transition temperature. Vortex interaction is short range on small and logarithmic on large scales. The chiral transition, though belonging to the Ising universality class by symmetry, has different critical exponents due to nonlocal interaction. In a narrow region close to the Lifshitz point a reentrant phase transition between paramagnetic and quasiferromagnetic phase appears. Applications to antiferromagnetic quantum spin chains and multiferroics are discussed. PMID:24785067

  9. Effect of River Training Project on Hydrodynamics Flow Circumstances by 2D Finite Element Numerical Model

    NASA Astrophysics Data System (ADS)

    Zou, B.; Li, D. F.; Hu, H. J.; Zhang, H. W.; Lou, L. H.; Chen, M.; Lv, Z. Y.

    Based on the verified two dimensional(2D) finite element model for river flow simulation, the effect of estuary training levees on the water flow and sediment movement in the Yellow River estuary is analyzed. For disclosing the effect of setting the two training levees on the flow and sediment motion, the calculation and analysis for the two projects, (one is no levees, the other is setting up two no levees) are given. The results show that when setting up two training levees, water flow is bound by levees and the water flows become more concentrated. As a result, velocity increases in the main channel, sediment carrying capacity of water flow increases correspondingly.

  10. Kosaki-Longo index and classification of charges in 2D quantum spin models

    NASA Astrophysics Data System (ADS)

    Naaijkens, Pieter

    2013-08-01

    We consider charge superselection sectors of two-dimensional quantum spin models corresponding to cone localisable charges, and prove that the number of equivalence classes of such charges is bounded by the Kosaki-Longo index of an inclusion of certain observable algebras. To demonstrate the power of this result we apply the theory to the toric code on a 2D infinite lattice. For this model we can compute the index of this inclusion, and conclude that there are four distinct irreducible charges in this model, in accordance with the analysis of the toric code model on compact surfaces. We also give a sufficient criterion for the non-degeneracy of the charge sectors, in the sense that Verlinde's matrix S is invertible.

  11. Synthesis of 2D/2D Structured Mesoporous Co3O4 Nanosheet/N-Doped Reduced Graphene Oxide Composites as a Highly Stable Negative Electrode for Lithium Battery Applications.

    PubMed

    Sennu, Palanichamy; Kim, Hyo Sang; An, Jae Youn; Aravindan, Vanchiappan; Lee, Yun-Sung

    2015-08-01

    Mesoporous Co3O4 nanosheets (Co3 O4 -NS) and nitrogen-doped reduced graphene oxide (N-rGO) are synthesized by a facile hydrothermal approach, and the N-rGO/Co3O4 -NS composite is formulated through an infiltration procedure. Eventually, the obtained composites are subjected to various characterization techniques, such as XRD, Raman spectroscopy, surface area analysis, X-ray photoelectron spectroscopy (XPS), and TEM. The lithium-storage properties of N-rGO/Co3O4 -NS composites are evaluated in a half-cell assembly to ascertain their suitability as a negative electrode for lithium-ion battery applications. The 2D/2D nanostructured mesoporous N-rGO/Co3O4 -NS composite delivered a reversible capacity of about 1305 and 1501 mAh g(-1) at a current density of 80 mA g(-1) for the 1st and 50th cycles, respectively. Furthermore, excellent cyclability, rate capability, and capacity retention characteristics are noted for the N-rGO/Co3O4 -NS composite. This improved performance is mainly related to the existence of mesoporosity and a sheet-like 2D hierarchical morphology, which translates into extra space for lithium storage and a reduced electron pathway. Also, the presence of N-rGO and carbon shells in Co3O4 -NS should not be excluded from such exceptional performance, which serves as a reliable conductive channel for electrons and act as synergistically to accommodate volume expansion upon redox reactions. Ex-situ TEM, impedance spectroscopy, and XPS, are also conducted to corroborate the significance of the 2D morphology towards sustained lithium storage. PMID:26033848

  12. Accelerating numerical modeling of wave propagation through 2-D anisotropic materials using OpenCL.

    PubMed

    Molero, Miguel; Iturrarán-Viveros, Ursula

    2013-03-01

    We present an implementation of the numerical modeling of elastic waves propagation, in 2D anisotropic materials, using the new parallel computing devices (PCDs). Our study is aimed both to model laboratory experiments and explore the capabilities of the emerging PCDs by discussing performance issues. In the experiments a sample plate of an anisotropic material placed inside a water tank is rotated and, for every angle of rotation it is subjected to an ultrasonic wave (produced by a large source transducer) that propagates in the water and through the material producing some reflection and transmission signals that are recording by a "point-like" receiver. This experiment is numerically modeled by running a finite difference code covering a set of angles θ∈[-50°, 50°], and recorded the signals for the transmission and reflection results. Transversely anisotropic and weakly orthorhombic materials are considered. We accelerated the computation using an open-source toolkit called PyOpenCL, which lets one to easily access the OpenCL parallel computation API's from the high-level programming environment of Python. A speedup factor over 19 using the GPU is obtained when compared with the execution of the same program in parallel using a CPU multi-core (in this case we use the 4-cores that has the CPU). The performance for different graphic cards and operating systems is included together with the full 2-D finite difference code with PyOpenCL. PMID:23290584

  13. Self-Organization in 2D Traffic Flow Model with Jam-Avoiding Drive

    NASA Astrophysics Data System (ADS)

    Nagatani, Takashi

    1995-04-01

    A stochastic cellular automaton (CA) model is presented to investigate the traffic jam by self-organization in the two-dimensional (2D) traffic flow. The CA model is the extended version of the 2D asymmetric exclusion model to take into account jam-avoiding drive. Each site contains either a car moving to the up, a car moving to the right, or is empty. A up car can shift right with probability p ja if it is blocked ahead by other cars. It is shown that the three phases (the low-density phase, the intermediate-density phase and the high-density phase) appear in the traffic flow. The intermediate-density phase is characterized by the right moving of up cars. The jamming transition to the high-density jamming phase occurs with higher density of cars than that without jam-avoiding drive. The jamming transition point p 2c increases with the shifting probability p ja. In the deterministic limit of p ja=1, it is found that a new jamming transition occurs from the low-density synchronized-shifting phase to the high-density moving phase with increasing density of cars. In the synchronized-shifting phase, all up cars do not move to the up but shift to the right by synchronizing with the move of right cars. We show that the jam-avoiding drive has an important effect on the dynamical jamming transition.

  14. Stochastic dynamics of phase singularities under ventricular fibrillation in 2D Beeler-Reuter model

    NASA Astrophysics Data System (ADS)

    Suzuki, Akio; Konno, Hidetoshi

    2011-09-01

    The dynamics of ventricular fibrillation (VF) has been studied extensively, and the initiation mechanism of VF has been elucidated to some extent. However, the stochastic dynamical nature of sustained VF remains unclear so far due to the complexity of high dimensional chaos in a heterogeneous system. In this paper, various statistical mechanical properties of sustained VF are studied numerically in 2D Beeler-Reuter-Drouhard-Roberge (BRDR) model with normal and modified ionic current conductance. The nature of sustained VF is analyzed by measuring various fluctuations of spatial phase singularity (PS) such as velocity, lifetime, the rates of birth and death. It is found that the probability density function (pdf) for lifetime of PSs is independent of system size. It is also found that the hyper-Gamma distribution serves as a universal pdf for the counting number of PSs for various system sizes and various parameters of our model tissue under VF. Further, it is demonstrated that the nonlinear Langevin equation associated with a hyper-Gamma process can mimic the pdf and temporal variation of the number of PSs in the 2D BRDR model.

  15. Adaptive finite element modeling of direct current resistivity in 2-D generally anisotropic structures

    NASA Astrophysics Data System (ADS)

    Yan, Bo; Li, Yuguo; Liu, Ying

    2016-07-01

    In this paper, we present an adaptive finite element (FE) algorithm for direct current (DC) resistivity modeling in 2-D generally anisotropic conductivity structures. Our algorithm is implemented on an unstructured triangular mesh that readily accommodates complex structures such as topography and dipping layers and so on. We implement a self-adaptive, goal-oriented grid refinement algorithm in which the finite element analysis is performed on a sequence of refined grids. The grid refinement process is guided by an a posteriori error estimator. The problem is formulated in terms of total potentials where mixed boundary conditions are incorporated. This type of boundary condition is superior to the Dirichlet type of conditions and improves numerical accuracy considerably according to model calculations. We have verified the adaptive finite element algorithm using a two-layered earth with azimuthal anisotropy. The FE algorithm with incorporation of mixed boundary conditions achieves high accuracy. The relative error between the numerical and analytical solutions is less than 1% except in the vicinity of the current source location, where the relative error is up to 2.4%. A 2-D anisotropic model is used to demonstrate the effects of anisotropy upon the apparent resistivity in DC soundings.

  16. Inhibitory effects of phytochemicals on metabolic capabilities of CYP2D6*1 and CYP2D6*10 using cell-based models in vitro

    PubMed Central

    Qu, Qiang; Qu, Jian; Han, Lu; Zhan, Min; Wu, Lan-xiang; Zhang, Yi-wen; Zhang, Wei; Zhou, Hong-hao

    2014-01-01

    Aim: Herbal products have been widely used, and the safety of herb-drug interactions has aroused intensive concerns. This study aimed to investigate the effects of phytochemicals on the catalytic activities of human CYP2D6*1 and CYP2D6*10 in vitro. Methods: HepG2 cells were stably transfected with CYP2D6*1 and CYP2D6*10 expression vectors. The metabolic kinetics of the enzymes was studied using HPLC and fluorimetry. Results: HepG2-CYP2D6*1 and HepG2-CYP2D6*10 cell lines were successfully constructed. Among the 63 phytochemicals screened, 6 compounds, including coptisine sulfate, bilobalide, schizandrin B, luteolin, schizandrin A and puerarin, at 100 μmol/L inhibited CYP2D6*1- and CYP2D6*10-mediated O-demethylation of a coumarin compound AMMC by more than 50%. Furthermore, the inhibition by these compounds was dose-dependent. Eadie-Hofstee plots demonstrated that these compounds competitively inhibited CYP2D6*1 and CYP2D6*10. However, their Ki values for CYP2D6*1 and CYP2D6*10 were very close, suggesting that genotype-dependent herb-drug inhibition was similar between the two variants. Conclusion: Six phytochemicals inhibit CYP2D6*1 and CYP2D6*10-mediated catalytic activities in a dose-dependent manner in vitro. Thus herbal products containing these phytochemicals may inhibit the in vivo metabolism of co-administered drugs whose primary route of elimination is CYP2D6. PMID:24786236

  17. Modeling and Control of 2-D Grasping of an Object with Arbitrary Shape under Rolling Contact

    NASA Astrophysics Data System (ADS)

    Arimoto, Suguru; Yoshida, Morio; Sekimoto, Masahiro; Tahara, Kenji

    Modeling, control, and stabilization of dynamics of two-dimensional object grasping by using a pair of multi-joint robot fingers are investigated under rolling contact constraints and an arbitrary geometry of the object and fingertips. First, a fundamental testbed problem of modeling and control of rolling motion between 2-D rigid bodies with an arbitrary shape is treated under the assumption that the two contour curves coincide at the contact point and share the same tangent. The rolling constraint induces the Euler equation of motion that is parameterized by a common arclength parameter and constrained onto the kernel space orthogonally complemented to the image space spanned from the constraint gradient. By extending the analysis to the problem of stable grasp of a 2-D object with an arbitrary shape by a pair of robot fingers, the Euler-Lagrange equation of motion of the overall fingers/object system parametrized by arclength parameters is derived, together with a couple of first-order differential equations that express evolutions of contact points in terms of the second fundamental form. It is shown that 2-D rolling constraints are integrable in the sense of Frobonius even if their Pfaffian forms are characterized by arclength parameters. A control signal called “blind grasping” is introduced and shown to be effective in stabilization of grasping without using the details of the object shape and parameters or external sensing. An extension of the Dirichlet-Lagrange stability theorem to a class of systems with DOF-redundancy under constraints is suggested by using a Morse-Bott-Lyapunov function.

  18. Uncertainties in modelling Mt. Pinatubo eruption with 2-D AER model and CCM SOCOL

    NASA Astrophysics Data System (ADS)

    Kenzelmann, P.; Weisenstein, D.; Peter, T.; Luo, B. P.; Rozanov, E.; Fueglistaler, S.; Thomason, L. W.

    2009-04-01

    Large volcanic eruptions may introduce a strong forcing on climate. They challenge the skills of climate models. In addition to the short time attenuation of solar light by ashes the formation of stratospheric sulphate aerosols, due to volcanic sulphur dioxide injection into the lower stratosphere, may lead to a significant enhancement of the global albedo. The sulphate aerosols have a residence time of about 2 years. As a consequence of the enhanced sulphate aerosol concentration both the stratospheric chemistry and dynamics are strongly affected. Due to absorption of longwave and near infrared radiation the temperature in the lower stratosphere increases. So far chemistry climate models overestimate this warming [Eyring et al. 2006]. We present an extensive validation of extinction measurements and model runs of the eruption of Mt. Pinatubo in 1991. Even if Mt. Pinatubo eruption has been the best quantified volcanic eruption of this magnitude, the measurements show considerable uncertainties. For instance the total amount of sulphur emitted to the stratosphere ranges from 5-12 Mt sulphur [e.g. Guo et al. 2004, McCormick, 1992]. The largest uncertainties are in the specification of the main aerosol cloud. SAGE II, for instance, could not measure the peak of the aerosol extinction for about 1.5 years, because optical termination was reached. The gap-filling of the SAGE II [Thomason and Peter, 2006] using lidar measurements underestimates the total extinctions in the tropics for the first half year after the eruption by 30% compared to AVHRR [Rusell et. al 1992]. The same applies to the optical dataset described by Stenchikov et al. [1998]. We compare these extinction data derived from measurements with extinctions derived from AER 2D aerosol model calculations [Weisenstein et al., 2007]. Full microphysical calculations with injections of 14, 17, 20 and 26 Mt SO2 in the lower stratosphere were performed. The optical aerosol properties derived from SAGE II

  19. The concept models and implementations of multiport neural net associative memory for 2D patterns

    NASA Astrophysics Data System (ADS)

    Krasilenko, Vladimir G.; Nikolskyy, Aleksandr I.; Yatskovskaya, Rimma A.; Yatskovsky, Victor I.

    2011-04-01

    The paper considers neural net models and training and recognizing algorithms with base neurobiologic operations: p-step autoequivalence and non-equivalenc The Modified equivalently models (MEMs) of multiport neural net associative memory (MNNAM) are offered with double adaptive - equivalently weighing (DAEW) for recognition of 2D-patterns (images). It is shown, the computing process in MNNAM under using the proposed MEMs, is reduced to two-step and multi-step algorithms and step-by-step matrix-matrix (tensor-tensor) procedures. The given results of computer simulations confirmed the perspective of such models. Besides the result was received when MNNAM capacity on base of MEMs exceeded the amount of neurons.

  20. Longtime Well-posedness for the 2D Groma-Balogh Model

    NASA Astrophysics Data System (ADS)

    Wan, Renhui; Chen, Jiecheng

    2016-07-01

    In this paper, we consider the cauchy problem for the 2D Groma-Balogh model (Acta Mater 47:3647-3654, 1999). From the works Cannone et al. (Arch Ration Mech Anal 196:71-96, 2010) and El Hajj (Ann Inst Henri Poincaré Anal Nonlinéaire 27:21-35, 2010), one can see global well-posedness for this model is an open question. However, we can prove longtime well-posedness. In particular, we show that this model admits a unique solution with the lifespan T^star satisfying T^star log ^2(1+T^star )≳ ɛ ^{-2} if the initial data is of size ɛ . To achieve this, we first establish some new decay estimates concerning the operator e^{-{R}_{12}^2t} . Then, we prove the longtime well-posedness by utilizing the weak dissipation to deal with the nonlinear terms.

  1. Well-posedness and generalized plane waves simulations of a 2D mode conversion model

    NASA Astrophysics Data System (ADS)

    Imbert-Gérard, Lise-Marie

    2015-12-01

    Certain types of electro-magnetic waves propagating in a plasma can undergo a mode conversion process. In magnetic confinement fusion, this phenomenon is very useful to heat the plasma, since it permits to transfer the heat at or near the plasma center. This work focuses on a mathematical model of wave propagation around the mode conversion region, from both theoretical and numerical points of view. It aims at developing, for a well-posed equation, specific basis functions to study a wave mode conversion process. These basis functions, called generalized plane waves, are intrinsically based on variable coefficients. As such, they are particularly adapted to the mode conversion problem. The design of generalized plane waves for the proposed model is described in detail. Their implementation within a discontinuous Galerkin method then provides numerical simulations of the process. These first 2D simulations for this model agree with qualitative aspects studied in previous works.

  2. Interpretation of gravity data using 2-D continuous wavelet transformation and 3-D inverse modeling

    NASA Astrophysics Data System (ADS)

    Roshandel Kahoo, Amin; Nejati Kalateh, Ali; Salajegheh, Farshad

    2015-10-01

    Recently the continuous wavelet transform has been proposed for interpretation of potential field anomalies. In this paper, we introduced a 2D wavelet based method that uses a new mother wavelet for determination of the location and the depth to the top and base of gravity anomaly. The new wavelet is the first horizontal derivatives of gravity anomaly of a buried cube with unit dimensions. The effectiveness of the proposed method is compared with Li and Oldenburg inversion algorithm and is demonstrated with synthetics and real gravity data. The real gravity data is taken over the Mobrun massive sulfide ore body in Noranda, Quebec, Canada. The obtained results of the 2D wavelet based algorithm and Li and Oldenburg inversion on the Mobrun ore body had desired similarities to the drill-hole depth information. In all of the inversion algorithms the model non-uniqueness is the challenging problem. Proposed method is based on a simple theory and there is no model non-uniqueness on it.

  3. Modeling floods in a dense urban area using 2D shallow water equations

    NASA Astrophysics Data System (ADS)

    Mignot, E.; Paquier, A.; Haider, S.

    2006-07-01

    SummaryA code solving the 2D shallow water equations by an explicit second-order scheme is used to simulate the severe October 1988 flood in the Richelieu urban locality of the French city of Nîmes. A reference calculation using a detailed description of the street network and of the cross-sections of the streets, considering impervious residence blocks and neglecting the flow interaction with the sewer network provides a mean peak water elevation 0.13 m lower than the measured flood marks with a standard deviation between the measured and computed water depths of 0.53 m. Sensitivity analysis of various topographical and numerical parameters shows that globally, the results keep the same level of accuracy, which reflects both the stability of the calculation method and the smoothening of results. However, the local flow modifications due to change of parameter values can drastically modify the local water depths, especially when the local flow regime is modified. Furthermore, the flow distribution to the downstream parts of the city can be altered depending on the set of parameters used. Finally, a second event, the 2002 flood, was simulated with the calibrated model providing results similar to 1988 flood calculation. Thus, the article shows that, after calibration, a 2D model can be used to help planning mitigation measures in a dense urban area.

  4. Modeling and 2-D discrete simulation of dislocation dynamics for plastic deformation of metal

    NASA Astrophysics Data System (ADS)

    Liu, Juan; Cui, Zhenshan; Ou, Hengan; Ruan, Liqun

    2013-05-01

    Two methods are employed in this paper to investigate the dislocation evolution during plastic deformation of metal. One method is dislocation dynamic simulation of two-dimensional discrete dislocation dynamics (2D-DDD), and the other is dislocation dynamics modeling by means of nonlinear analysis. As screw dislocation is prone to disappear by cross-slip, only edge dislocation is taken into account in simulation. First, an approach of 2D-DDD is used to graphically simulate and exhibit the collective motion of a large number of discrete dislocations. In the beginning, initial grains are generated in the simulation cells according to the mechanism of grain growth and the initial dislocation is randomly distributed in grains and relaxed under the internal stress. During the simulation process, the externally imposed stress, the long range stress contribution of all dislocations and the short range stress caused by the grain boundaries are calculated. Under the action of these forces, dislocations begin to glide, climb, multiply, annihilate and react with each other. Besides, thermal activation process is included. Through the simulation, the distribution of dislocation and the stress-strain curves can be obtained. On the other hand, based on the classic dislocation theory, the variation of the dislocation density with time is described by nonlinear differential equations. Finite difference method (FDM) is used to solve the built differential equations. The dislocation evolution at a constant strain rate is taken as an example to verify the rationality of the model.

  5. 2D and 3D shape based segmentation using deformable models.

    PubMed

    El-Baz, Ayman; Yuksel, Seniha E; Shi, Hongjian; Farag, Aly A; El-Ghar, Mohamed A; Eldiasty, Tarek; Ghoneim, Mohamed A

    2005-01-01

    A novel shape based segmentation approach is proposed by modifying the external energy component of a deformable model. The proposed external energy component depends not only on the gray level of the images but also on the shape information which is obtained from the signed distance maps of objects in a given data set. The gray level distribution and the signed distance map of the points inside and outside the object of interest are accurately estimated by modelling the empirical density function with a linear combination of discrete Gaussians (LCDG) with positive and negative components. Experimental results on the segmentation of the kidneys from low-contrast DCE-MRI and on the segmentation of the ventricles from brain MRI's show how the approach is accurate in segmenting 2-D and 3-D data sets. The 2D results for the kidney segmentation have been validated by a radiologist and the 3D results of the ventricle segmentation have been validated with a geometrical phantom. PMID:16686036

  6. A 2D Electromechanical Model of Human Atrial Tissue Using the Discrete Element Method

    PubMed Central

    Brocklehurst, Paul; Adeniran, Ismail; Yang, Dongmin; Sheng, Yong; Zhang, Henggui; Ye, Jianqiao

    2015-01-01

    Cardiac tissue is a syncytium of coupled cells with pronounced intrinsic discrete nature. Previous models of cardiac electromechanics often ignore such discrete properties and treat cardiac tissue as a continuous medium, which has fundamental limitations. In the present study, we introduce a 2D electromechanical model for human atrial tissue based on the discrete element method (DEM). In the model, single-cell dynamics are governed by strongly coupling the electrophysiological model of Courtemanche et al. to the myofilament model of Rice et al. with two-way feedbacks. Each cell is treated as a viscoelastic body, which is physically represented by a clump of nine particles. Cell aggregations are arranged so that the anisotropic nature of cardiac tissue due to fibre orientations can be modelled. Each cell is electrically coupled to neighbouring cells, allowing excitation waves to propagate through the tissue. Cell-to-cell mechanical interactions are modelled using a linear contact bond model in DEM. By coupling cardiac electrophysiology with mechanics via the intracellular Ca2+ concentration, the DEM model successfully simulates the conduction of cardiac electrical waves and the tissue's corresponding mechanical contractions. The developed DEM model is numerically stable and provides a powerful method for studying the electromechanical coupling problem in the heart. PMID:26583141

  7. Locally adaptive 2D-3D registration using vascular structure model for liver catheterization.

    PubMed

    Kim, Jihye; Lee, Jeongjin; Chung, Jin Wook; Shin, Yeong-Gil

    2016-03-01

    Two-dimensional-three-dimensional (2D-3D) registration between intra-operative 2D digital subtraction angiography (DSA) and pre-operative 3D computed tomography angiography (CTA) can be used for roadmapping purposes. However, through the projection of 3D vessels, incorrect intersections and overlaps between vessels are produced because of the complex vascular structure, which makes it difficult to obtain the correct solution of 2D-3D registration. To overcome these problems, we propose a registration method that selects a suitable part of a 3D vascular structure for a given DSA image and finds the optimized solution to the partial 3D structure. The proposed algorithm can reduce the registration errors because it restricts the range of the 3D vascular structure for the registration by using only the relevant 3D vessels with the given DSA. To search for the appropriate 3D partial structure, we first construct a tree model of the 3D vascular structure and divide it into several subtrees in accordance with the connectivity. Then, the best matched subtree with the given DSA image is selected using the results from the coarse registration between each subtree and the vessels in the DSA image. Finally, a fine registration is conducted to minimize the difference between the selected subtree and the vessels of the DSA image. In experimental results obtained using 10 clinical datasets, the average distance errors in the case of the proposed method were 2.34±1.94mm. The proposed algorithm converges faster and produces more correct results than the conventional method in evaluations on patient datasets. PMID:26824922

  8. 2-D Modeling of Nanoscale MOSFETs: Non-Equilibrium Green's Function Approach

    NASA Technical Reports Server (NTRS)

    Svizhenko, Alexei; Anantram, M. P.; Govindan, T. R.; Biegel, Bryan

    2001-01-01

    We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions and oxide tunneling are treated on an equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. Electron-electron interaction is treated within Hartree approximation by solving NEGF and Poisson equations self-consistently. For the calculations presented here, parallelization is performed by distributing the solution of NEGF equations to various processors, energy wise. We present simulation of the "benchmark" MIT 25nm and 90nm MOSFETs and compare our results to those from the drift-diffusion simulator and the quantum-corrected results available. In the 25nm MOSFET, the channel length is less than ten times the electron wavelength, and the electron scattering time is comparable to its transit time. Our main results are: (1) Simulated drain subthreshold current characteristics are shown, where the potential profiles are calculated self-consistently by the corresponding simulation methods. The current predicted by our quantum simulation has smaller subthreshold slope of the Vg dependence which results in higher threshold voltage. (2) When gate oxide thickness is less than 2 nm, gate oxide leakage is a primary factor which determines off-current of a MOSFET (3) Using our 2-D NEGF simulator, we found several ways to drastically decrease oxide leakage current without compromising drive current. (4) Quantum mechanically calculated electron density is much smaller than the background doping density in the poly silicon gate region near oxide interface. This creates an additional effective gate voltage. Different ways to. include this effect approximately will be discussed.

  9. Be2D: A model to understand the distribution of meteoric 10Be in soilscapes

    NASA Astrophysics Data System (ADS)

    Campforts, Benjamin; Vanacker, Veerle; Vanderborght, Jan; Govers, Gerard

    2016-04-01

    Cosmogenic nuclides have revolutionised our understanding of earth surface process rates. They have become one of the standard tools to quantify soil production by weathering, soil redistribution and erosion. Especially Beryllium-10 has gained much attention due to its long half-live and propensity to be relatively conservative in the landscape. The latter makes 10Be an excellent tool to assess denudation rates over the last 1000 to 100 × 103 years, bridging the anthropogenic and geological time scale. Nevertheless, the mobility of meteoric 10Be in soil systems makes translation of meteoric 10Be inventories into erosion and deposition rates difficult. Here we present a coupled soil hillslope model, Be2D, that is applied to synthetic and real topography to address the following three research questions. (i) What is the influence of vertical meteoric Be10 mobility, caused by chemical mobility, clay translocation and bioturbation, on its lateral redistribution over the soilscape, (ii) How does vertical mobility influence erosion rates and soil residence times inferred from meteoric 10Be inventories and (iii) To what extent can a tracer with a half-life of 1.36 Myr be used to distinguish between natural and human-disturbed soil redistribution rates? The model architecture of Be2D is designed to answer these research questions. Be2D is a dynamic model including physical processes such as soil formation, physical weathering, clay migration, bioturbation, creep, overland flow and tillage erosion. Pathways of meteoric 10Be mobility are simulated using a two step approach which is updated each timestep. First, advective and diffusive mobility of meteoric 10Be is simulated within the soil profile and second, lateral redistribution because of lateral soil fluxes is calculated. The performance and functionality of the model is demonstrated through a number of synthetic and real model runs using existing datasets of meteoric 10Be from case-studies in southeastern US. Brute

  10. Estimating nitrogen losses in furrow irrigated soil amended by compost using HYDRUS-2D model

    NASA Astrophysics Data System (ADS)

    Iqbal, Shahid; Guber, Andrey; Zaman Khan, Haroon; ullah, Ehsan

    2014-05-01

    Furrow irrigation commonly results in high nitrogen (N) losses from soil profile via deep infiltration. Estimation of such losses and their reduction is not a trivial task because furrow irrigation creates highly nonuniform distribution of soil water that leads to preferential water and N fluxes in soil profile. Direct measurements of such fluxes are impractical. The objective of this study was to assess applicability of HYDRUS-2D model for estimating nitrogen balance in manure amended soil under furrow irrigation. Field experiments were conducted in a sandy loam soil amended by poultry manure compost (PMC) and pressmud compost (PrMC) fertilizers. The PMC and PrMC contained 2.5% and 0.9% N and were applied at 5 rates: 2, 4, 6, 8 and 10 ton/ha. Plots were irrigated starting from 26th day from planting using furrows with 1x1 ridge to furrow aspect ratio. Irrigation depths were 7.5 cm and time interval between irrigations varied from 8 to 15 days. Results of the field experiments showed that approximately the same corn yield was obtained with considerably higher N application rates using PMC than using PrMC as a fertilizer. HYDRUS-2D model was implemented to evaluate N fluxes in soil amended by PMC and PrMC fertilizers. Nitrogen exchange between two pools of organic N (compost and soil) and two pools of mineral N (soil NH4-N and soil NO3-N) was modeled using mineralization and nitrification reactions. Sources of mineral N losses from soil profile included denitrification, root N uptake and leaching with deep infiltration of water. HYDRUS-2D simulations showed that the observed increases in N root water uptake and corn yields associated with compost application could not be explained by the amount of N added to soil profile with the compost. Predicted N uptake by roots significantly underestimated the field data. Good agreement between simulated and field-estimated values of N root uptake was achieved when the rate of organic N mineralization was increased

  11. A guide to using material model No. 11 in NIKE2D: An internal variable, viscoplasticity model

    SciTech Connect

    Flower, E.C.; Nikkel, D.J. Jr.

    1990-10-30

    The need to accurately model the superplastic forming process which is highly rate and temperature dependent motivated the evaluation of Bammann's internal variable, viscoplasticity material model. The model is based upon the concepts of unified creep plasticity, but employs a yield surface for efficient implementation into large-scale numerical computer codes. It has proven elsewhere to be quite successful in describing large strain, thermal-mechanical behavior of crystalline materials. Features of the model enable it to simulate the apparent strain-rate behavior exhibited by many metals above one half the melt temperature. It is the efficient incorporation of features that make the model attractive for use in finite element modeling of metal deformation processes. Although this model was implemented into the Lawrence Livermore National Laboratory's NIKE2D finite element program in 1986, there have been no known reports of successful use by NIKE2D users. The purpose of this report is to provide the user the proper format to input model parameters, a procedure for determining appropriate values for material constants from experimental data, and supplemental information on the model relevant to the implementation in the NIKE2D finite element program. Detailed accounts of the theoretical aspects of the model can be found in the cited references. 4 refs., 8 figs.

  12. Singularities of the Partition Function for the Ising Model Coupled to 2D Quantum Gravity

    NASA Astrophysics Data System (ADS)

    Ambjørn, J.; Anagnostopoulos, K. N.; Magnea, U.

    We study the zeros in the complex plane of the partition function for the Ising model coupled to 2D quantum gravity for complex magnetic field and real temperature, and for complex temperature and real magnetic field, respectively. We compute the zeros by using the exact solution coming from a two-matrix model and by Monte-Carlo simulations of Ising spins on dynamical triangulations. We present evidence that the zeros form simple one-dimensional curves in the complex plane, and that the critical behaviour of the system is governed by the scaling of the distribution of the singularities near the critical point. Despite the small size of the systems studied, we can obtain a reasonable estimate of the (known) critical exponents.

  13. Comparative modeling of vertical and planar organic phototransistors with 2D drift-diffusion simulations

    NASA Astrophysics Data System (ADS)

    Bezzeccheri, E.; Colasanti, S.; Falco, A.; Liguori, R.; Rubino, A.; Lugli, P.

    2016-05-01

    Vertical Organic Transistors and Phototransistors have been proven to be promising technologies due to the advantages of reduced channel length and larger sensitive area with respect to planar devices. Nevertheless, a real improvement of their performance is subordinate to the quantitative description of their operation mechanisms. In this work, we present a comparative study on the modeling of vertical and planar Organic Phototransistor (OPT) structures. Computer-based simulations of the devices have been carried out with Synopsys Sentaurus TCAD in a 2D Drift-Diffusion framework. The photoactive semiconductor material has been modeled using the virtual semiconductor approach as the archetypal P3HT:PC61BM bulk heterojunction. It has been found that both simulated devices have comparable electrical and optical characteristics, accordingly to recent experimental reports on the subject.

  14. An investigation of DTNS2D for use as an incompressible turbulence modelling test-bed

    NASA Technical Reports Server (NTRS)

    Steffen, Christopher J., Jr.

    1992-01-01

    This paper documents an investigation of a two dimensional, incompressible Navier-Stokes solver for use as a test-bed for turbulence modelling. DTNS2D is the code under consideration for use at the Center for Modelling of Turbulence and Transition (CMOTT). This code was created by Gorski at the David Taylor Research Center and incorporates the pseudo compressibility method. Two laminar benchmark flows are used to measure the performance and implementation of the method. The classical solution of the Blasius boundary layer is used for validating the flat plate flow, while experimental data is incorporated in the validation of backward facing step flow. Velocity profiles, convergence histories, and reattachment lengths are used to quantify these calculations. The organization and adaptability of the code are also examined in light of the role as a numerical test-bed.

  15. Neutrino-electron Scattering in 2-D Models of Supernova Convection

    NASA Astrophysics Data System (ADS)

    DeNisco, K. R.; Swesty, F. D.; Myra, E. S.

    2005-12-01

    We present results from 2-D supernova simulations which include the effects of neutrino-electron scattering. The importance of neutrino-electron scattering in stellar collapse has been known for two decades. Yet it has often been neglected in multidimensional simulations due to the difficulty of implementing it consistently. The inclusion of this process is numerically challenging because of the extremely short scattering timescales involved. The stiffness resulting from this short timescale precludes an explicit numerical treatment of this phenomenon, such as those that have recently been utilized in some 2-D models. We describe our fully-implicit treatment of this process and present our initial results. This work was performed at the State University of New York at Stony Brook as part of the TeraScale Supernova Initiative, and is funded by SciDAC grant DE-FC02-01ER41185 from the U.S. Department of Energy, Office of Science High-Energy, Nuclear, and Advanced Scientific Computing Research Programs. We gratefully acknowledge support of the National Energy Research Scientific Computing Center (NERSC) for computational and consulting support.

  16. Implications of lack-of-ergodicity in 2D Potts model

    NASA Astrophysics Data System (ADS)

    Ota, Smita

    2015-03-01

    Microcanonical Monte Carlo simulation is used to study two dimensional (2D) q state Potts model. We consider a 2D square lattice having NxN spins with periodic boundary condition and simulated the system with N =15 and q =10. The demon energy distribution is found to be exponential for high system energy and large system size. For smaller system size and above the first order transition the demon energy distribution is found to deviate from exp(- βED) and has the form exp(- βED + γ ED2). Here β = 1/kBT and kB is the Boltzmann constant. It is found that γ is finite at higher temperatures. As the system energy is reduced γ becomes zero near the first order transition. It is found that during cooling γ changes sign from negative to positive and then to negative again near the 1st order transition. Therefore the demon energy distribution becomes exp(- βED) (or ergodic) at two values of system energy near the 1st order transition. Further cooling or at still lower temperatures the system shows lack of ergodicity. However, difference in heating cooling curves are apparent in E vs γ. The system energies for which γ is zero during cooling can represent the 'ergodic' states. This can be related to the two-level systems observed in glasses at low temperatures.

  17. GRAV2D: an interactive 2-1/2 dimensional gravity modeling program (user's guide and documentation for Rev. 1)

    SciTech Connect

    Nutter, C.

    1980-11-01

    GRAV2D is an interactive computer program used for modeling 2-1/2 dimensional gravity data. A forward algorithm is used to give the theoretical attraction of gravity intensity at a station due to a perturbing body given by the initial model. The resultant model can then be adjusted for a better fit by a combination of manual adjustment, one-dimensional automatic search, and Marquardt inversion. GRAV2D has an interactive data management system for data manipulation and display built around subroutines to do a forward problem, a one-dimensional direct search and an inversion. This is a user's guide and documentation for GRAV2D.

  18. Fabrication method and microstructural characteristics of coal-tar-pitch-based 2D carbon/carbon composites

    NASA Astrophysics Data System (ADS)

    Esmaeeli, Mohammad; Khosravi, Hamed; Mirhabibi, Alireza

    2015-02-01

    The lignin-cellulosic texture of wood was used to produce two-dimensional (2D) carbon/carbon (C/C) composites using coal tar pitch. Ash content tests were conducted to select two samples among the different kinds of woods present in Iran, including walnut, white poplar, cherry, willow, buttonwood, apricots, berry, and blue wood. Walnut and white poplar with ash contents of 1.994wt% and 0.351wt%, respectively, were selected. The behavior of these woods during pyrolysis was investigated by differential thermal analysis (DTA) and thermo gravimetric (TG) analysis. The bulk density and open porosity were measured after carbonization and densification. The microstructural characteristics of samples were investigated by scanning electron microscopy (SEM), X-ray diffraction (XRD), and Fourier-transform infrared (FT-IR) spectroscopy. The results indicate that the density of both the walnut and white poplar is increased, and the open porosity is decreased with the increasing number of carbonization cycles. The XRD patterns of the wood charcoal change gradually with increasing pyrolysis temperature, possibly as a result of the ultra-structural changes in the charcoal or the presence of carbonized coal tar pitch in the composite's body.

  19. 2D time-domain finite-difference modeling for viscoelastic seismic wave propagation

    NASA Astrophysics Data System (ADS)

    Fan, Na; Zhao, Lian-Feng; Xie, Xiao-Bi; Ge, Zengxi; Yao, Zhen-Xing

    2016-07-01

    Real Earth media are not perfectly elastic. Instead, they attenuate propagating mechanical waves. This anelastic phenomenon in wave propagation can be modeled by a viscoelastic mechanical model consisting of several standard linear solids. Using this viscoelastic model, we approximate a constant Q over a frequency band of interest. We use a four-element viscoelastic model with a tradeoff between accuracy and computational costs to incorporate Q into 2D time-domain first-order velocity-stress wave equations. To improve the computational efficiency, we limit the Q in the model to a list of discrete values between 2 and 1000. The related stress and strain relaxation times that characterize the viscoelastic model are pre-calculated and stored in a database for use by the finite-difference calculation. A viscoelastic finite-difference scheme that is second-order in time and fourth-order in space is developed based on the MacCormack algorithm. The new method is validated by comparing the numerical result with analytical solutions that are calculated using the generalized reflection/transmission coefficient method. The synthetic seismograms exhibit greater than 95 per cent consistency in a two-layer viscoelastic model. The dispersion generated from the simulation is consistent with the Kolsky-Futterman dispersion relationship.

  20. Transforming 2d Cadastral Data Into a Dynamic Smart 3d Model

    NASA Astrophysics Data System (ADS)

    Tsiliakou, E.; Labropoulos, T.; Dimopoulou, E.

    2013-08-01

    3D property registration has become an imperative need in order to optimally reflect all complex cases of the multilayer reality of property rights and restrictions, revealing their vertical component. This paper refers to the potentials and multiple applications of 3D cadastral systems and explores the current state-of-the art, especially the available software with which 3D visualization can be achieved. Within this context, the Hellenic Cadastre's current state is investigated, in particular its data modeling frame. Presenting the methodologies and specifications addressing the registration of 3D properties, the operating cadastral system's shortcomings and merits are pointed out. Nonetheless, current technological advances as well as the availability of sophisticated software packages (proprietary or open source) call for 3D modeling. In order to register and visualize the complex reality in 3D, Esri's CityEngine modeling software has been used, which is specialized in the generation of 3D urban environments, transforming 2D GIS Data into Smart 3D City Models. The application of the 3D model concerns the Campus of the National Technical University of Athens, in which a complex ownership status is established along with approved special zoning regulations. The 3D model was built using different parameters based on input data, derived from cadastral and urban planning datasets, as well as legal documents and architectural plans. The process resulted in a final 3D model, optimally describing the cadastral situation and built environment and proved to be a good practice example of 3D visualization.

  1. Dynamic Linkages Between the Transition Zone & Surface Plate Motion in 2D Models of Subduction

    NASA Astrophysics Data System (ADS)

    Arredondo, K.; Billen, M. I.

    2014-12-01

    Subduction zones exhibit a wide range of behavior, from slab stagnation at 660 km to direct penetration into the lower mantle. Due to uncertainties in the tectonic history of individual subduction zones, such as trench velocities, potential mechanisms for controlling slab behavior in the transition zone are explored using numerical models. Numerical simulations have utilized a range of assumptions to improve computational efficiency, such as ignoring latent heat, ignoring compositional effects or fixing the trench location: the net effect of these assumptions resulting modeled dynamics remains unclear. Additionally the eight major, composition-dependent, phase transitions for pyrolite, harzburgite and eclogite may be an important influence on subducting slab dynamics due to the additional forces that are dependent on depth and compositional layering within the slab (e.g., Ricard et al., 2005). With the goal of developing more complete, self-consistent, and less idealized simulations, we test the importance of various factors on slab behavior: the presence of shear, adiabatic and latent heating, compositional layering, composition-dependent phase transitions and explicit plate speeds versus dynamically evolving plate and trench velocities. Preliminary results indicate that individual components have a relatively minor effect, but produce large changes when combined together. The extent of slab folding and stagnation is overestimated by only modeling the 410 and 660 km phase transitions. Dynamic models with all seven composition-dependent phase transitions are very sensitive to the plate strength and weak zone viscosity, causing large changes in plate speed and slab detachment. Changes to the overriding plate buoyance and strength investigate the origin and influence of trench movement on slab deformation. These feedbacks and parameter-sensitive behavior indicate that the wide range of observed slab behavior may result from subtle differences in plate and plate

  2. The development and testing of a 2D laboratory seismic modelling system for heterogeneous structure investigations

    NASA Astrophysics Data System (ADS)

    Mo, Yike; Greenhalgh, Stewart A.; Robertsson, Johan O. A.; Karaman, Hakki

    2015-05-01

    Lateral velocity variations and low velocity near-surface layers can produce strong scattered and guided waves which interfere with reflections and lead to severe imaging problems in seismic exploration. In order to investigate these specific problems by laboratory seismic modelling, a simple 2D ultrasonic model facility has been recently assembled within the Wave Propagation Lab at ETH Zurich. The simulated geological structures are constructed from 2 mm thick metal and plastic sheets, cut and bonded together. The experiments entail the use of a piezoelectric source driven by a pulse amplifier at ultrasonic frequencies to generate Lamb waves in the plate, which are detected by piezoelectric receivers and recorded digitally on a National Instruments recording system, under LabVIEW software control. The 2D models employed were constructed in-house in full recognition of the similitude relations. The first heterogeneous model features a flat uniform low velocity near-surface layer and deeper dipping and flat interfaces separating different materials. The second model is comparable but also incorporates two rectangular shaped inserts, one of low velocity, the other of high velocity. The third model is identical to the second other than it has an irregular low velocity surface layer of variable thickness. Reflection as well as transmission experiments (crosshole & vertical seismic profiling) were performed on each model. The two dominant Lamb waves recorded are the fundamental symmetric mode (non-dispersive) and the fundamental antisymmetric (flexural) dispersive mode, the latter normally being absent when the source transducer is located on a model edge but dominant when it is on the flat planar surface of the plate. Experimental group and phase velocity dispersion curves were determined and plotted for both modes in a uniform aluminium plate. For the reflection seismic data, various processing techniques were applied, as far as pre-stack Kirchhoff migration. The

  3. Turbulence modeling for subsonic separated flows over 2-D airfoils and 3-D wings

    NASA Astrophysics Data System (ADS)

    Rosen, Aaron M.

    Accurate predictions of turbulent boundary layers and flow separation through computational fluid dynamics (CFD) are becoming more and more essential for the prediction of loads in the design of aerodynamic flight components. Standard eddy viscosity models used in many commercial codes today do not capture the nonequilibrium effects seen in a separated flow and thus do not generally make accurate separation predictions. Part of the reason for this is that under nonequilibrium conditions such as a strong adverse pressure gradient, the history effects of the flow play an important role in the growth and decay of turbulence. More recent turbulence models such as Olsen and Coakley's Lag model and Lillard's lagRST model seek to simulate these effects by lagging the turbulent variables when nonequilibrium effects become important. The purpose of the current research is to assess how these nonequilibrium turbulence models capture the separated regions on various 2-D airfoils and 3-D wings. Nonequilibrium models including the Lag model and the lagRST model are evaluated in comparison with three baseline models (Spalart-Allmaras, Wilcox's k-omega, and Menter's SST) using a modified version of the OVERFLOW code. Tuning the model coefficients of the Lag and lagRST models is also explored. Results show that the various lagRST formulations display an improvement in velocity profile predictions over the standard RANS models, but have trouble capturing the edge of the boundary layer. Experimental separation location measurements were not available, but several trends are noted which may be useful to tuning the model coefficients in the future.

  4. On the assimilation of flood extension images into 2D shallow-water models

    NASA Astrophysics Data System (ADS)

    Monnier, J.; Couderc, F.; Dartus, D.; Madec, R.; Vila, J.

    2012-12-01

    In river hydraulics, assimilation of water level measurements at gauging stations is well controlled, while assimilation of images (e.g. from satellite) is still delicate. In the present talk, we address the richness of satellite information to constraint a 2D shallow-water model, and present also related difficulties. A preliminary study done on Mosel river is presented in [LaMo] [HoLaMoPu]. On selected parts of the image, an 0th order model flow allows to obtain some reliable water levels with quantified uncertainties (C. Puech et al.). Next, variationnal sensitivities (based on a gradient computation and adjoint equations) reveal some difficulties that a model designer have to tackle (e.g. roughness parameters at open boundaries), and allow to better understand both the model and the flow. Next, a variational data assimilation algorithm (4D-var) shows that such data lead to a better calibration of the model (e.g. roughness coefficients) and potentially allows to identify the incoming and/or outgoing flow at open boundaries, [LaMo] [HoLaMoPu]. On the other side, the flood dynamic extension is difficult to represent accurately using a 2D SW model since the wet-dry front dynamics is difficult to compute. We compare some 2nd order finite volume solvers and obtain an accurate and stable scheme at wet-dry front. Then, we present some basic rules of compatibility between data and mesh resolution in order to be reliable enough to constraint the model with flood extension data, [CoMaMoViDa]. All the algorithms are implemented into DassFlow software (Fortran, MPI, adjoint) [Da]. [CoMaMoViDa] F. Couderc, R. Madec, J. Monnier, J.-P. Vila, D. Dartus. "Sensitivity analysis and variational data assimilation for geophysical shallow water flows". Submitted. [Da] DassFlow - Data Assimilation for Free Surface Flows. Open-source computational software http://www-gmm.insa-toulouse.fr/~monnier/DassFlow/ [HoLaMoPu] R. Hostache, X. Lai, J. Monnier, C. Puech. "Assimilation of spatial

  5. Coronary arteries motion modeling on 2D x-ray images

    NASA Astrophysics Data System (ADS)

    Gao, Yang; Sundar, Hari

    2012-02-01

    During interventional procedures, 3D imaging modalities like CT and MRI are not commonly used due to interference with the surgery and radiation exposure concerns. Therefore, real-time information is usually limited and building models of cardiac motion are difficult. In such case, vessel motion modeling based on 2-D angiography images become indispensable. Due to issues with existing vessel segmentation algorithms and the lack of contrast in occluded vessels, manual segmentation of certain branches is usually necessary. In addition, such occluded branches are the most important vessels during coronary interventions and obtaining motion models for these can greatly help in reducing the procedure time and radiation exposure. Segmenting different cardiac phases independently does not guarantee temporal consistency and is not efficient for occluded branches required manual segmentation. In this paper, we propose a coronary motion modeling system which extracts the coronary tree for every cardiac phase, maintaining the segmentation by tracking the coronary tree during the cardiac cycle. It is able to map every frame to the specific cardiac phase, thereby inferring the shape information of the coronary arteries using the model corresponding to its phase. Our experiments show that our motion modeling system can achieve promising results with real-time performance.

  6. 2-D modeling of laterally acoustically coupled thin film bulk acoustic wave resonator filters.

    PubMed

    Pensala, Tuomas; Meltaus, Johanna; Kokkonen, Kimmo; Ylilammi, Markku

    2010-11-01

    A 2-D model is developed for calculating lateral acoustical coupling between adjacent thin film BAW resonators forming an electrical N-port. The model is based on solution and superposition of lateral eigenmodes and eigenfrequencies in a structure consisting of adjacent regions with known plate wave dispersion properties. Mechanical and electrical response of the device are calculated as a superposition of eigenmodes according to voltage drive at one electrical port at a time while extracting current induced in the other ports, leading to a full Y-parameter description of the device. Exemplary cases are simulated to show the usefulness of the model in the study of the basic design rules of laterally coupled thin film BAW resonator filters. Model predictions are compared to an experimental 1.9-GHz band-pass filter based on aluminum nitride thin film technology and lateral acoustical coupling. Good agreement is obtained in prediction of passband behavior. The eigenmode-based model forms a useful tool for fast simulation of laterally coupled acoustic devices. It allows one to gain insight into basic device physics in a very intuitive fashion compared with more detailed but heavier finite element method. Shortcomings of this model and possible improvements are discussed. PMID:21041141

  7. The combined effect of attraction and orientation zones in 2D flocking models

    NASA Astrophysics Data System (ADS)

    Iliass, Tarras; Cambui, Dorilson

    2016-01-01

    In nature, many animal groups, such as fish schools or bird flocks, clearly display structural order and appear to move as a single coherent entity. In order to understand the complex motion of these systems, we study the Vicsek model of self-propelled particles (SPP) which is an important tool to investigate the behavior of collective motion of live organisms. This model reproduces the biological behavior patterns in the two-dimensional (2D) space. Within the framework of this model, the particles move with the same absolute velocity and interact locally in the zone of orientation by trying to align their direction with that of the neighbors. In this paper, we model the collective movement of SPP using an agent-based model which follows biologically motivated behavioral rules, by adding a second region called the attraction zone, where each particles move towards each other avoiding being isolated. Our main goal is to present a detailed numerical study on the effect of the zone of attraction on the kinetic phase transition of our system. In our study, the consideration of this zone seems to play an important role in the cohesion. Consequently, in the directional orientation, the zone that we added forms the compact particle group. In our simulation, we show clearly that the model proposed here can produce two collective behavior patterns: torus and dynamic parallel group. Implications of these findings are discussed.

  8. Models Ion Trajectories in 2D and 3D Electrostatic and Magnetic Fields

    SciTech Connect

    Dahl, David

    2000-02-21

    SIMION3D7.0REV is a C based ion optics simulation program that can model complex problems using Laplace equation solutions for potential fields. The program uses an ion optics workbench that can hold up to 200 2D and/or 3D electrostatic/magnetic potential arrays. Arrays can have up to 50,000,000 points. SIMION3D7.0''s 32 bit virtual Graphics User Interface provides a highly interactive advanced user environment. All potential arrays are visualized as 3D objects that the user can cut away to inspect ion trajectories and potential energy surfaces. User programs allow the user to customize the program for specific simulations. A geometry file option supports the definition of highly complex array geometry. Algorithm modifications have improved this version''s computational speed and accuracy.

  9. Robust autonomous model learning from 2D and 3D data sets.

    PubMed

    Langs, Georg; Donner, René; Peloschek, Philipp; Bischof, Horst

    2007-01-01

    In this paper we propose a weakly supervised learning algorithm for appearance models based on the minimum description length (MDL) principle. From a set of training images or volumes depicting examples of an anatomical structure, correspondences for a set of landmarks are established by group-wise registration. The approach does not require any annotation. In contrast to existing methods no assumptions about the topology of the data are made, and the topology can change throughout the data set. Instead of a continuous representation of the volumes or images, only sparse finite sets of interest points are used to represent the examples during optimization. This enables the algorithm to efficiently use distinctive points, and to handle texture variations robustly. In contrast to standard elasticity based deformation constraints the MDL criterion accounts for systematic deformations typical for training sets stemming from medical image data. Experimental results are reported for five different 2D and 3D data sets. PMID:18051152

  10. Calibration Of 2D Hydraulic Inundation Models In The Floodplain Region Of The Lower Tagus River

    NASA Astrophysics Data System (ADS)

    Pestanana, R.; Matias, M.; Canelas, R.; Araujo, A.; Roque, D.; Van Zeller, E.; Trigo-Teixeira, A.; Ferreira, R.; Oliveira, R.; Heleno, S.

    2013-12-01

    In terms of inundated area, the largest floods in Portugal occur in the Lower Tagus River. On average, the river overflows every 2.5 years, at times blocking roads and causing important agricultural damages. This paper focus on the calibration of 2D-horizontal flood simulation models for the floods of 2001 and 2006 on a 70-km stretch of the Lower Tagus River. Flood extent maps, derived from ERS SAR and ENVISAT ASAR imagery were compared with the flood extent maps obtained for each simulation, to calibrate roughness coefficients. The combination of the calibration results from the 2001 and 2006 floods provided a preliminary Manning coefficient map of the study area.

  11. Models Ion Trajectories in 2D and 3D Electrostatic and Magnetic Fields

    Energy Science and Technology Software Center (ESTSC)

    2000-02-21

    SIMION3D7.0REV is a C based ion optics simulation program that can model complex problems using Laplace equation solutions for potential fields. The program uses an ion optics workbench that can hold up to 200 2D and/or 3D electrostatic/magnetic potential arrays. Arrays can have up to 50,000,000 points. SIMION3D7.0''s 32 bit virtual Graphics User Interface provides a highly interactive advanced user environment. All potential arrays are visualized as 3D objects that the user can cut awaymore » to inspect ion trajectories and potential energy surfaces. User programs allow the user to customize the program for specific simulations. A geometry file option supports the definition of highly complex array geometry. Algorithm modifications have improved this version''s computational speed and accuracy.« less

  12. Surface delta interaction in the g7/2 - d5/2 model space

    NASA Astrophysics Data System (ADS)

    Yu, Xiaofei; Zamick, Larry

    2016-05-01

    Using an attractive surface delta interaction we obtain wave functions for 2 neutrons (or neutron holes) in the g7/2 -d5/2 model space. If we take the single particle energies to be degenerate we find that the g factors for I = 2 , 4 and 6 are all the same G (J) =gl, the orbital g factor of the nucleon. For a free neutron gl = 0, so in this case all 2 particles or 2 holes' g factors are equal to zero. Only the orbital part of the g-factors contributes - the spin part cancels out. We then consider the effects of introducing a single energy splitting between the 2 orbits. We make a linear approximation for all other n values.

  13. 3D prostate boundary segmentation from ultrasound images using 2D active shape models.

    PubMed

    Hodge, Adam C; Ladak, Hanif M

    2006-01-01

    Boundary outlining, or segmentation, of the prostate is an important task in diagnosis and treatment planning for prostate cancer. This paper describes an algorithm for semi-automatic, three-dimensional (3D) segmentation of the prostate boundary from ultrasound images based on two-dimensional (2D) active shape models (ASM) and rotation-based slicing. Evaluation of the algorithm used distance- and volume-based error metrics to compare algorithm generated boundary outlines to gold standard (manually generated) boundary outlines. The mean absolute distance between the algorithm and gold standard boundaries was 1.09+/-0.49 mm, the average percent absolute volume difference was 3.28+/-3.16%, and a 5x speed increase as compared manual planimetry was achieved. PMID:17946106

  14. Optical fiber poling by induction: analysis by 2D numerical modeling.

    PubMed

    De Lucia, F; Huang, D; Corbari, C; Healy, N; Sazio, P J A

    2016-04-15

    Since their first demonstration some 25 years ago, thermally poled silica fibers have been used to realize device functions such as electro-optic modulation, switching, polarization-entangled photons, and optical frequency conversion with a number of advantages over bulk free-space components. We have recently developed an innovative induction poling technique that could allow for the development of complex microstructured fiber geometries for highly efficient χ(2)-based device applications. To systematically implement these more advanced poled fiber designs, we report here the development of comprehensive numerical models of the induction poling mechanism itself via two-dimensional (2D) simulations of ion migration and space-charge region formation using finite element analysis. PMID:27082323

  15. Structure-approximating inverse protein folding problem in the 2D HP model.

    PubMed

    Gupta, Arvind; Manuch, Ján; Stacho, Ladislav

    2005-12-01

    The inverse protein folding problem is that of designing an amino acid sequence which has a particular native protein fold. This problem arises in drug design where a particular structure is necessary to ensure proper protein-protein interactions. In this paper, we show that in the 2D HP model of Dill it is possible to solve this problem for a broad class of structures. These structures can be used to closely approximate any given structure. One of the most important properties of a good protein (in drug design) is its stability--the aptitude not to fold simultaneously into other structures. We show that for a number of basic structures, our sequences have a unique fold. PMID:16379538

  16. Universality Class of the Nishimori Point in the 2D +/-J Random-Bond Ising Model

    NASA Astrophysics Data System (ADS)

    Honecker, A.; Picco, M.; Pujol, P.

    2001-07-01

    We study the universality class of the Nishimori point in the 2D +/-J random-bond Ising model by means of the numerical transfer-matrix method. Using the domain-wall free energy, we locate the position of the fixed point along the Nishimori line at the critical concentration value pc = 0.1094+/-0.0002 and estimate ν = 1.33+/-0.03. Then, we obtain the exponents for the moments of the spin-spin correlation functions as well as the value for the central charge c = 0.464+/-0.004. The main qualitative result is the fact that percolation is now excluded as a candidate for describing the universality class of this fixed point.

  17. Universality class of the Nishimori point in the 2D +/- J random-bond Ising model.

    PubMed

    Honecker, A; Picco, M; Pujol, P

    2001-07-23

    We study the universality class of the Nishimori point in the 2D +/- J random-bond Ising model by means of the numerical transfer-matrix method. Using the domain-wall free energy, we locate the position of the fixed point along the Nishimori line at the critical concentration value p(c) = 0.1094 +/- 0.0002 and estimate nu = 1.33 +/- 0.03. Then, we obtain the exponents for the moments of the spin-spin correlation functions as well as the value for the central charge c = 0.464 +/- 0.004. The main qualitative result is the fact that percolation is now excluded as a candidate for describing the universality class of this fixed point. PMID:11461639

  18. An application of the distributed hydrologic model CASC2D to a tropical montane watershed

    NASA Astrophysics Data System (ADS)

    Marsik, Matt; Waylen, Peter

    2006-11-01

    SummaryIncreased stormflow in the Quebrada Estero watershed (2.5 km 2), in the northwestern Central Valley tectonic depression of Costa Rica, reportedly has caused flooding of the city of San Ramón in recent decades. Although scientifically untested, urban expansion was deemed the cause and remedial measures were recommended by the Programa de Investigación en Desarrollo Humano Sostenible (ProDUS). CASC2D, a physically-based, spatially explicit hydrologic model, was constructed and calibrated to a June 10th 2002 storm that delivered 110.5 mm of precipitation in 4.5 h visibly exceeded the bankfull stage (0.9 m) of the Quebrada flooding portions of San Ramón. The calibrated hydrograph showed a peak discharge 16.68% (2.5 m 3 s -1) higher, an above flood stage duration 20% shorter, and time to peak discharge 11 min later than the same observed discharge hydrograph characteristics. Simulations of changing land cover conditions from 1979 to 1999 showed an increase also in the peak discharge, above flood stage duration, and time to peak discharge. Analysis using a modified location quotient identified increased urbanization in lower portions of the watershed over the time period studied. These results suggest that increased urbanization in the Quebrada Estero watershed have increased flooding peaks, and durations above threshold, confirming the ProDUS report. These results and the CASC2D model offer an easy-to-use, pragmatic planning tool for policymakers in San Ramón to assess future development scenarios and their potential flooding impacts to San Ramón.

  19. Assessment of the Impacts of Compensation Flow Changes Upon Instream Habitat Using 2D Modelling

    NASA Astrophysics Data System (ADS)

    Mould, D. C.; Lane, S. N.; Christmas, M.

    2004-05-01

    Many millstone-grit rivers in northern England are impounded. In such cases the water company in the area has to release compensation flows from the reservoirs, traditionally to meet industrial needs: these flows are rarely set with ecology in mind; and have commonly involved constant flow. Dam overtopping may create spates, but spawning in many fish species is prompted by a spate flow in the early autumn when dams are rarely full enough to overtop. Such flows are important for fine sediment flushing and controlling the wetted useable area for spawning. Classical physical habitat modelling for instream habitat has been largely reliant upon 1D approaches, such as the Instream Flow Incremental Methodology (IFIM). Here we use a 2D finite element model (FESWMS), to simulate changes in instream habitat with variations in the compensation flow regimes. The spatial resolution of 2D models can be adapted to the scale of fish habitats so providing better representation of the reach-scale flow processes (such as slack water in the margins, wetting and drying) than the 1D case. The model is applied to the Rivers Rivelin and Loxley in Sheffield, Northern England. At the confluence of the two rivers, the compensation flow level is set at 30.6 Thousand Cubic Metres per Day (TCMD). Due to historical reasons, the compensation is not divided equally, as the Loxley receives 28 TCMD whilst the Rivelin receives only 2.6 TCMD. The model is used to simulate a transfer of 6 TCMD from the Loxley to the Rivelin. After validation, model predictions are combined with available habitat requirement data (e.g. velocity and depth needs) to develop an index of change in habitat suitability in terms of first order variables (e.g. velocity, depth and wetted useable area). This suggests that the change in compensation may significantly improve instream ecology in relation to macroinvertebrates, brown trout (Salmo trutta) and bullhead (Cottus gobio) in the Rivelin without causing detrimental impacts

  20. 2D Distinct Element Method (DEM) models of the initiation, propagation and saturation of rock joints

    NASA Astrophysics Data System (ADS)

    Arslan, A.; Schöpfer, M. P.; Walsh, J. J.; Childs, C.

    2009-12-01

    In layered sequences, rock joints usually best develop within the more brittle layers and commonly display a regular spacing that scales with layer thickness. A variety of conceptual and mechanical models have been developed for these observations. A limitation of previous approaches, however, is that fracture initiation and associated interface slip are not explicitly simulated; instead, fractures were predefined and interfaces were welded. To surmount this problem, we have modelled the formation and growth of joints in layered sequences by using the two-dimensional Distinct Element Method (DEM) as implemented in the Particle Flow Code (PFC-2D). In PFC-2D, rock is represented by an assemblage of circular particles that are bonded at particle-particle contacts. Failure occurs if either the tensile or shear strength of a bond is exceeded. The models comprise a central brittle layer with high Young’s modulus, which is embedded in a low Young’s modulus matrix. The interfaces between the layers are defined by ‘smooth joint’ contacts, a modelling feature that eliminates interparticle bumpiness and associated interlocking friction. Consequently, this feature allows the user to assign macroscopic properties such as friction and cohesion along layer interfaces in a controlled manner. Layer parallel extension is applied by assigning a velocity to particles at the lateral boundaries of the model while maintaining a constant vertical confining pressure. Models were extended until joint saturation in the central layer was reached. We thereby explored the impact of confining pressure and interface properties (friction, cohesion) on joint spacing. A number of important conclusions can be drawn from our models: (i) The distributions of average horizontal normal stress within the layer and of shear stress at the interface are consistent with analytical solutions (stress-transfer theory). (ii) At low interfacial shear strength, new joints form preferentially midway between

  1. Field Evaluation of a Novel 2D Preferential Flow Snowpack Hydrology Model

    NASA Astrophysics Data System (ADS)

    Leroux, N.; Pomeroy, J. W.; Kinar, N. J.

    2015-12-01

    Accurate estimation of snowmelt flux is of primary importance for runoff hydrograph prediction, which is used for water management and flood forecasting. Lateral flows and preferential flow pathways in porous media flow have proven critical for improving soil and groundwater flow models, but though many physically-based layered snowmelt models have been developed, only 1D matrix flow is accounted for in these models. Therefore, there is a need for snowmelt models that include these processes so as to examine the potential to improve snowmelt hydrological modelling. A 2D model is proposed that enables an improved understanding of energy and water flows within deep heterogeneous snowpacks, including those on slopes. A dual pathway theory is presented that simulates the formation of preferential flow paths, vertical and lateral water flows through the snow matrix and flow fingers, internal energy fluxes, melt, wet snow metamorphism, and internal refreezing. The dual pathway model utilizes an explicit finite volume method to solve for the energy and water flux equations over a non-orthogonal grid. It was run and evaluated using in-situ data collected from snowpit - accessed gravimetric, thermometric, photographic, and dielectric observations and novel non-invasive acoustic observations of layering, temperature, flowpath geometry, density and wetness at the Fortress Mountain Snow Laboratory, Alberta, Canada. The melt of a natural snowpack was artificially generated after detailed observation of snowpack initial conditions such as snow layer properties, temperature, and liquid water content. Snowpack ablation and liquid water content distribution over time were then measured and used for model parameterization and validation. Energy available at the snow surface and soil slope angle were set as mondel inputs. Model verification was based on snowpack property evolution. The heterogeneous flow model can be an important tool to help understand snowmelt flow processes, how

  2. Thermochemical Nonequilibrium 2D Modeling of Nitrogen Inductively Coupled Plasma Flow

    NASA Astrophysics Data System (ADS)

    Yu, Minghao; Yusuke, Takahashi; Hisashi, Kihara; Ken-ichi, Abe; Kazuhiko, Yamada; Takashi, Abe; Satoshi, Miyatani

    2015-09-01

    Two-dimensional (2D) numerical simulations of thermochemical nonequilibrium inductively coupled plasma (ICP) flows inside a 10-kW inductively coupled plasma wind tunnel (ICPWT) were carried out with nitrogen as the working gas. Compressible axisymmetric Navier-Stokes (N-S) equations coupled with magnetic vector potential equations were solved. A four-temperature model including an improved electron-vibration relaxation time was used to model the internal energy exchange between electron and heavy particles. The third-order accuracy electron transport properties (3rd AETP) were applied to the simulations. A hybrid chemical kinetic model was adopted to model the chemical nonequilibrium process. The flow characteristics such as thermal nonequilibrium, inductive discharge, effects of Lorentz force were made clear through the present study. It was clarified that the thermal nonequilibrium model played an important role in properly predicting the temperature field. The prediction accuracy can be improved by applying the 3rd AETP to the simulation for this ICPWT. supported by Grant-in-Aid for Scientific Research (No. 23560954), sponsored by the Japan Society for the Promotion of Science

  3. Incorporating a Turbulence Transport Model into 2-D Hybrid Hall Thruster Simulations

    NASA Astrophysics Data System (ADS)

    Cha, Eunsun; Cappelli, Mark A.; Fernandez, Eduardo

    2014-10-01

    2-D hybrid simulations of Hall plasma thrusters that do not resolve cross-field transport-generating fluctuations require a model to capture how electrons migrate across the magnetic field. We describe the results of integrating a turbulent electron transport model into simulations of plasma behavior in a plane spanned by the E and B field vectors. The simulations treat the electrons as a fluid and the heavy species (ions/neutrals) as discrete particles. The transport model assumes that the turbulent eddy cascade in the electron fluid to smaller scales is the primary means of electron energy dissipation. Using this model, we compare simulations to experimental measurements made on a laboratory Hall discharge over a range of discharge voltage. Both the current-voltage trends as well as the plasma properties such as plasma temperature, electron density, and ion velocities seem agree favorably with experiments, where a simple Bohm transport model tends to perform poorly in capturing much of the discharge behavior.

  4. Spin Circuit Model for Spin Orbit Torques in 2D Channels

    NASA Astrophysics Data System (ADS)

    Hong, Seokmin

    2015-03-01

    Recently, the unique coupling between charge and spin in topological insulators has been explored through various types of electrical measurements, which could have interesting applications. In this talk, we present a spin circuit model for spin orbit torques in topological insulator surface states and other 2D channels. We show with a simple example that results from the circuit model agree well with those obtained from nonequilibrium Green's function (NEGF) based quantum transport simulation. Some predictions of our model have already received experimental support and we hope this model can provide a unifying framework that can be used to critically evaluate experimental results, to explore new types of devices as well as to answer fundamental questions regarding these materials. The model for spin-orbit torques described here can be incorporated into a broader spin-circuit approach which, we believe, provides a natural platform for multi-physics, multi-component spintronic devices. This work was supported by FAME, a Semiconductor Research Corporation program sponsored by MARCO and DARPA.

  5. Spin Circuit Model for 2D Channels with Spin-Orbit Coupling

    PubMed Central

    Hong, Seokmin; Sayed, Shehrin; Datta, Supriyo

    2016-01-01

    In this paper we present a general theory for an arbitrary 2D channel with “spin momentum locking” due to spin-orbit coupling. It is based on a semiclassical model that classifies all the channel electronic states into four groups based on the sign of the z-component of the spin (up (U), down (D)) and the sign of the x-component of the velocity (+, −). This could be viewed as an extension of the standard spin diffusion model which uses two separate electrochemical potentials for U and D states. Our model uses four: U+, D+, U−, and D−. We use this formulation to develop an equivalent spin circuit that is also benchmarked against a full non-equilibrium Green’s function (NEGF) model. The circuit representation can be used to interpret experiments and estimate important quantities of interest like the charge to spin conversion ratio or the maximum spin current that can be extracted. The model should be applicable to topological insulator surface states with parallel channels as well as to other layered structures with interfacial spin-orbit coupling. PMID:26932563

  6. Spin Circuit Model for 2D Channels with Spin-Orbit Coupling.

    PubMed

    Hong, Seokmin; Sayed, Shehrin; Datta, Supriyo

    2016-01-01

    In this paper we present a general theory for an arbitrary 2D channel with "spin momentum locking" due to spin-orbit coupling. It is based on a semiclassical model that classifies all the channel electronic states into four groups based on the sign of the z-component of the spin (up (U), down (D)) and the sign of the x-component of the velocity (+, -). This could be viewed as an extension of the standard spin diffusion model which uses two separate electrochemical potentials for U and D states. Our model uses four: U+, D+, U-, and D-. We use this formulation to develop an equivalent spin circuit that is also benchmarked against a full non-equilibrium Green's function (NEGF) model. The circuit representation can be used to interpret experiments and estimate important quantities of interest like the charge to spin conversion ratio or the maximum spin current that can be extracted. The model should be applicable to topological insulator surface states with parallel channels as well as to other layered structures with interfacial spin-orbit coupling. PMID:26932563

  7. Spin Circuit Model for 2D Channels with Spin-Orbit Coupling

    NASA Astrophysics Data System (ADS)

    Hong, Seokmin; Sayed, Shehrin; Datta, Supriyo

    2016-03-01

    In this paper we present a general theory for an arbitrary 2D channel with “spin momentum locking” due to spin-orbit coupling. It is based on a semiclassical model that classifies all the channel electronic states into four groups based on the sign of the z-component of the spin (up (U), down (D)) and the sign of the x-component of the velocity (+, -). This could be viewed as an extension of the standard spin diffusion model which uses two separate electrochemical potentials for U and D states. Our model uses four: U+, D+, U-, and D-. We use this formulation to develop an equivalent spin circuit that is also benchmarked against a full non-equilibrium Green’s function (NEGF) model. The circuit representation can be used to interpret experiments and estimate important quantities of interest like the charge to spin conversion ratio or the maximum spin current that can be extracted. The model should be applicable to topological insulator surface states with parallel channels as well as to other layered structures with interfacial spin-orbit coupling.

  8. Integration of 2-D hydraulic model and high-resolution LiDAR-derived DEM for floodplain flow modeling

    NASA Astrophysics Data System (ADS)

    Shen, D.; Wang, J.; Cheng, X.; Rui, Y.; Ye, S.

    2015-02-01

    The rapid progress of Light Detection And Ranging (LiDAR) technology has made acquirement and application of high-resolution digital elevation model (DEM) data increasingly popular, especially with regards to the study of floodplain flow modeling. High-resolution DEM data include many redundant interpolation points, needs a high amount of calculation, and does not match the size of computational mesh. These disadvantages are a common problem for floodplain flow modeling studies. Two-dimensional (2-D) hydraulic modeling, a popular method of analyzing floodplain flow, offers high precision of elevation parameterization for computational mesh while ignoring much micro-topographic information of the DEM data itself. We offer a flood simulation method that integrates 2-D hydraulic model results and high-resolution DEM data, enabling the calculation of flood water levels in DEM grid cells through local inverse distance weighted interpolation. To get rid of the false inundation areas during interpolation, it employs the run-length encoding method to mark the inundated DEM grid cells and determine the real inundation areas through the run-length boundary tracing technique, which solves the complicated problem of the connectivity between DEM grid cells. We constructed a 2-D hydraulic model for the Gongshuangcha polder, a flood storage area of Dongting Lake, using our integrated method to simulate the floodplain flow. The results demonstrate that this method can solve DEM associated problems efficiently and simulate flooding processes with greater accuracy than DEM only simulations.

  9. A new model for two-dimensional numerical simulation of pseudo-2D gas-solids fluidized beds

    SciTech Connect

    Li, Tingwen; Zhang, Yongmin

    2013-10-11

    Pseudo-two dimensional (pseudo-2D) fluidized beds, for which the thickness of the system is much smaller than the other two dimensions, is widely used to perform fundamental studies on bubble behavior, solids mixing, or clustering phenomenon in different gas-solids fluidization systems. The abundant data from such experimental systems are very useful for numerical model development and validation. However, it has been reported that two-dimensional (2D) computational fluid dynamic (CFD) simulations of pseudo-2D gas-solids fluidized beds usually predict poor quantitative agreement with the experimental data, especially for the solids velocity field. In this paper, a new model is proposed to improve the 2D numerical simulations of pseudo-2D gas-solids fluidized beds by properly accounting for the frictional effect of the front and back walls. Two previously reported pseudo-2D experimental systems were simulated with this model. Compared to the traditional 2D simulations, significant improvements in the numerical predictions have been observed and the predicted results are in better agreement with the available experimental data.

  10. A case study of fluid flow in fractured rock mass based on 2-D DFN modeling

    NASA Astrophysics Data System (ADS)

    Han, Jisu; Noh, Young-Hwan; Um, Jeong-Gi; Choi, Yosoon

    2014-05-01

    A two dimensional steady-state fluid flow through fractured rock mass of an abandoned copper mine in Korea is addressed based on discrete fracture network modeling. An injection well and three observation wells were installed at the field site to monitor the variations of total heads induced by injection of fresh water. A series of packer tests were performed to estimate the rock mass permeability. First, the two dimensional stochastic fracture network model was built and validated for a granitic rock mass using the geometrical and statistical data obtained from surface exposures and borehole logs. This validated fracture network model was combined with the fracture data observed on boreholes to generate a stochastic-deterministic fracture network system. Estimated apertures for each of the fracture sets using permeability data obtained from borehole packer tests were discussed next. Finally, a systematic procedure for fluid flow modeling in fractured rock mass in two dimensional domain was presented to estimate the conductance, flow quantity and nodal head in 2-D conceptual linear pipe channel network. The results obtained in this study clearly show that fracture geometry parameters (orientation, density and size) play an important role in the hydraulic behavior of fractured rock masses.

  11. Distributed and coupled 2D electro-thermal model of power semiconductor devices

    NASA Astrophysics Data System (ADS)

    Belkacem, Ghania; Lefebvre, Stéphane; Joubert, Pierre-Yves; Bouarroudj-Berkani, Mounira; Labrousse, Denis; Rostaing, Gilles

    2014-05-01

    The development of power electronics in the field of transportations (automotive, aeronautics) requires the use of power semiconductor devices providing protection and diagnostic functions. In the case of series protections power semiconductor devices which provide protection may operate in shortcircuit and act as a current limiting device. This mode of operations is very constraining due to the large dissipation of power. In these particular conditions of operation, electro-thermal models of power semiconductor devices are of key importance in order to optimize their thermal design and increase their reliability. The development of such an electro-thermal model for power MOSFET transistors based on the coupling between two computation softwares (Matlab and Cast3M) is described in this paper. The 2D electro-thermal model is able to predict (i) the temperature distribution on chip surface well as in the volume under short-circuit operations, (ii) the effect of the temperature on the distribution of the current flowing within the die and (iii) the effects of the ageing of the metallization layer on the current density and the temperature. In this paper, the electrical and thermal models are described as well as the implemented coupling scheme.

  12. Studies of Arctic stratospheric ozone in a 2-D model including some effects of zonal asymmetries

    SciTech Connect

    Isaksen, I.S.A.; Rognerud, B.; Stordal, F. ); Coffey, M.T.; Mankin, W.G. )

    1990-03-01

    A two-dimensional (2-D) zonally averaged chemistry-transport model of the stratosphere has been extended to include some zonally asymmetric effects to study the chemically disturbed conditions in the Arctic winter during the occurrence of polar stratospheric clouds (PSCs). The model allows air parcels that have been in PSCs in the polar night to be exposed to sunlight during the passage south through a wave trough. Large enhancements of ClO are estimated as well as significant ozone reductions, most pronounced around the 20 km height level. The ozone depletions maximize in late March, about one month after the cease in PSC activity in the model, and amount to 5-8% in column ozone at 70{degree}N. In agreement with column measurements made from the DC-8, the model estimates an increase in the columns of HNO{sub 3} and ClONO{sub 2}, and a decrease in the HCl column within the polar vortex.

  13. LBQ2D, Extending the Line Broadened Quasilinear Model to TAE-EP Interaction

    NASA Astrophysics Data System (ADS)

    Ghantous, Katy; Gorelenkov, Nikolai; Berk, Herbert

    2012-10-01

    The line broadened quasilinear model was proposed and tested on the one dimensional electrostatic case of the bump on tailfootnotetextH.L Berk, B. Breizman and J. Fitzpatrick, Nucl. Fusion, 35:1661, 1995 to study the wave particle interaction. In conventional quasilinear theory, the sea of overlapping modes evolve with time as the particle distribution function self consistently undergo diffusion in phase space. The line broadened quasilinear model is an extension to the conventional theory in a way that allows treatment of isolated modes as well as overlapping modes by broadening the resonant line in phase space. This makes it possible to treat the evolution of modes self consistently from onset to saturation in either case. We describe here the model denoted by LBQ2D which is an extension of the proposed one dimensional line broadened quasilinear model to the case of TAEs interacting with energetic particles in two dimensional phase space, energy as well as canonical angular momentum. We study the saturation of isolated modes in various regimes and present the analytical derivation and numerical results. Finally, we present, using ITER parameters, the case where multiple modes overlap and describe the techniques used for the numerical treatment.

  14. Obtaining instantaneous water levels relative to a geoid with a 2D storm surge model

    NASA Astrophysics Data System (ADS)

    Slobbe, D. C.; Verlaan, M.; Klees, R.; Gerritsen, H.

    2013-01-01

    Current and new applications of 2D storm surge models such as the Dutch Continental Shelf Model (DCSM) require that the models provide proper estimates of the instantaneous water levels expressed relative to a particular geoid, rather than only the tide and surge components expressed relative to the ill-defined model's zero height surface. For DCSM, this is realized by adding the depth-averaged horizontal baroclinic pressure gradients to the model equations, which are derived from 4D salinity and temperature fields provided by the Proudman Oceanographic Laboratory hydrodynamic model (POL's hindcast). The vertical datum of the extended model is fixed to that of the European Gravimetric Geoid model 2008 (EGG08). This is done by an adjustment of the model parameters that depend on the choice of the reference surface (e.g., bathymetry) and by referring the water levels along the open boundaries to this reference surface. Using different numerical experiments we investigate the effects on the water levels of several approximations we have made during the implementation. The ability of the model to reproduce both the mean sea level (MSL) and instantaneous water levels is assessed by a comparison with the MSL derived from POL's hindcast as well as with instantaneous water levels acquired by various radar altimeter satellites. From this comparison we conclude that our modeled MSL is in good agreement with the MSL derived from POL's hindcast; the standard deviation of the differences is below 2 cm. However, larger differences in MSL are observed when comparing the model output with the MSL derived from radar altimeter data. They are attributed to either geoid errors or errors in the used salinity and temperature fields. The root mean squared (rms) differences between observed and modeled instantaneous water levels over the entire model domain varies from 9 cm for data acquired by the TOPEX satellite to 11 cm for data acquired by the GFO-1 satellite. These numbers improve to

  15. Dynamic Linkages Between the Transition Zone & Surface Plate Motions in 2D Models of Subduction

    NASA Astrophysics Data System (ADS)

    Arredondo, K.; Billen, M. I.

    2013-12-01

    While slab pull is considered the dominant force controlling plate motion and speed, its magnitude is controlled by slab behavior in the mantle, where tomographic studies show a wide range of possibilities from direct penetration to folding, or stagnation directly above the lower mantle (e.g. Fukao et al., 2009). Geodynamic studies have investigated various parameters, such as plate age and two phase transitions, to recreate observed behavior (e.g. Běhounková and Cízková, 2008). However, past geodynamic models have left out known slab characteristics that may have a large impact on slab behavior and our understanding of subduction processes. Mineral experiments and seismic observations have indicated the existence of additional phase transitions in the mantle transition zone that may produce buoyancy forces large enough to affect the descent of a subducting slab (e.g. Ricard et al., 2005). The current study systematically tests different common assumptions used in geodynamic models: kinematic versus free-slip boundary conditions, the effects of adiabatic heating, viscous dissipation and latent heat, compositional layering and a more complete suite of phase transitions. Final models have a complete energy equation, with eclogite, harzburgite and pyrolite lithosphere compositional layers, and seven composition-dependent phase transitions within the olivine, pyroxene and garnet polymorph minerals. Results show important feedback loops between different assumptions and new behavior from the most complete models. Kinematic models show slab weakening or breaking above the 660 km boundary and between compositional layers. The behavior in dynamic models with a free-moving trench and overriding plate is compared to the more commonly found kinematic models. The new behavior may have important implications for the depth distribution of deep earthquakes within the slab. Though the thermodynamic parameters of certain phase transitions may be uncertain, their presence and

  16. A New 2D-Advection-Diffusion Model Simulating Trace Gas Distributions in the Lowermost Stratosphere

    NASA Astrophysics Data System (ADS)

    Hegglin, M. I.; Brunner, D.; Peter, T.; Wirth, V.; Fischer, H.; Hoor, P.

    2004-12-01

    Tracer distributions in the lowermost stratosphere are affected by both, transport (advective and non-advective) and in situ sources and sinks. They influence ozone photochemistry, radiative forcing, and heating budgets. In-situ measurements of long-lived species during eight measurement campaigns revealed relatively simple behavior of the tracers in the lowermost stratosphere when represented in an equivalent-latitude versus potential temperature framework. We here present a new 2D-advection-diffusion model that simulates the main transport pathways influencing the tracer distributions in the lowermost stratosphere. The model includes slow diabatic descent of aged stratospheric air and vertical and/or horizontal diffusion across the tropopause and within the lowermost stratosphere. The diffusion coefficients used in the model represent the combined effects of different processes with the potential of mixing tropospheric air into the lowermost stratosphere such as breaking Rossby and gravity waves, deep convection penetrating the tropopause, turbulent diffusion, radiatively driven upwelling etc. They were specified by matching model simulations to observed distributions of long-lived trace gases such as CO and N2O obtained during the project SPURT. The seasonally conducted campaigns allow us to study the seasonal dependency of the diffusion coefficients. Despite its simplicity the model yields a surprisingly good description of the small scale features of the measurements and in particular of the observed tracer gradients at the tropopause. The correlation coefficients between modeled and measured trace gas distributions were up to 0.95. Moreover, mixing across isentropes appears to be more important than mixing across surfaces of constant equivalent latitude (or PV). With the aid of the model, the distribution of the fraction of tropospheric air in the lowermost stratosphere can be determined.

  17. Transectional heat transfer in thermoregulating bigeye tuna (Thunnus obesus) - a 2D heat flux model.

    PubMed

    Boye, Jess; Musyl, Michael; Brill, Richard; Malte, Hans

    2009-11-01

    We developed a 2D heat flux model to elucidate routes and rates of heat transfer within bigeye tuna Thunnus obesus Lowe 1839 in both steady-state and time-dependent settings. In modeling the former situation, we adjusted the efficiencies of heat conservation in the red and the white muscle so as to make the output of the model agree as closely as possible with observed cross-sectional isotherms. In modeling the latter situation, we applied the heat exchanger efficiencies from the steady-state model to predict the distribution of temperature and heat fluxes in bigeye tuna during their extensive daily vertical excursions. The simulations yielded a close match to the data recorded in free-swimming fish and strongly point to the importance of the heat-producing and heat-conserving properties of the white muscle. The best correspondence between model output and observed data was obtained when the countercurrent heat exchangers in the blood flow pathways to the red and white muscle retained 99% and 96% (respectively) of the heat produced in these tissues. Our model confirms that the ability of bigeye tuna to maintain elevated muscle temperatures during their extensive daily vertical movements depends on their ability to rapidly modulate heating and cooling rates. This study shows that the differential cooling and heating rates could be fully accounted for by a mechanism where blood flow to the swimming muscles is either exclusively through the heat exchangers or completely shunted around them, depending on the ambient temperature relative to the body temperature. Our results therefore strongly suggest that such a mechanism is involved in the extensive physiological thermoregulatory abilities of endothermic bigeye tuna. PMID:19880733

  18. Complex 2D matrix model and geometrical map on the complex-Nc plane

    NASA Astrophysics Data System (ADS)

    Nawa, Kanabu; Ozaki, Sho; Nagahiro, Hideko; Jido, Daisuke; Hosaka, Atsushi

    2013-08-01

    We study the parameter dependence of the internal structure of resonance states by formulating a complex two-dimensional (2D) matrix model, where the two dimensions represent two levels of resonances. We calculate a critical value of the parameter at which a "nature transition" with character exchange occurs between two resonance states, from the viewpoint of geometry on complex-parameter space. Such a critical value is useful for identifying the internal structure of resonance states with variation of the parameter in the system. We apply the model to analyze the internal structure of hadrons with variation of the color number N_c from infty to a realistic value 3. By regarding 1/N_c as the variable parameter in our model, we calculate a critical color number of the nature transition between hadronic states in terms of a quark-antiquark pair and a mesonic molecule as exotics from the geometry on the complex-N_c plane. For large-N_c effective theory, we employ the chiral Lagrangian induced by holographic QCD with a D4/D8/overline {D8} multi-D brane system in type IIA superstring theory.

  19. Comparing a 2D fluid model of the DC planar magnetron cathode to experiments

    SciTech Connect

    Garcia, M.

    1996-05-01

    Planar magnetron cathodes have arching magnetic field lines which concentrate plasma density near the electrode surface. This enhances the ion bombardment of the surface and the yield of sputtered atoms. Magnetron cathodes are used in the Plasma Electrode Pockels Cell (PEPC) devices of the Laser Program because they provide for significantly higher conduction than do glow discharges. An essential feature of magnetron cathodes is that the vector product of the perpendicular electric field, E[sub y], with the parallel component of the magnetic field, B[sub x], forms a closed track with a circulating current along the cathode surface. An analytical, 2D, two component, quasi-neutral, continuum model yields formulas for the plasma density, the total and component current densities, the electric field, and the positive electrical potential, between the cathode surface and a distant, uniform plasma. For a specific gas, the free parameters are electron temperature, gas number density, and total current. The model is applied to the interpretation of experimental data from the PEPC device, as well as a small vacuum facility for testing magnetron cathodes. Finally, the model has been applied to generate cross sectional views of a PEPC magnetron cathode track.

  20. Doubled CO2 Effects on NO(y) in a Coupled 2D Model

    NASA Technical Reports Server (NTRS)

    Rosenfield, J. E.; Douglass, A. R.

    1998-01-01

    Changes in temperature and ozone have been the main focus of studies of the stratospheric impact of doubled CO2. Increased CO2 is expected to cool the stratosphere, which will result in increases in stratospheric ozone through temperature dependent loss rates. Less attention has been paid to changes in minor constituents which affect the O3 balance and which may provide additional feedbacks. Stratospheric NO(y) fields calculated using the GSFC 2D interactive chemistry-radiation-dynamics model show significant sensitivity to the model CO2. Modeled upper stratospheric NO(y) decreases by about 15% in response to CO2 doubling, mainly due to the temperature decrease calculated to result from increased cooling. The abundance of atomic nitrogen, N, increases because the rate of the strongly temperature dependent reaction N + O2 yields NO + O decreases at lower temperatures. Increased N leads to an increase in the loss of NO(y) which is controlled by the reaction N + NO yields N2 + O. The NO(y) reduction is shown to be sensitive to the NO photolysis rate. The decrease in the O3 loss rate due to the NO(y) changes is significant when compared to the decrease in the O3 loss rate due to the temperature changes.

  1. Simulation of abrasive flow machining process for 2D and 3D mixture models

    NASA Astrophysics Data System (ADS)

    Dash, Rupalika; Maity, Kalipada

    2015-12-01

    Improvement of surface finish and material removal has been quite a challenge in a finishing operation such as abrasive flow machining (AFM). Factors that affect the surface finish and material removal are media viscosity, extrusion pressure, piston velocity, and particle size in abrasive flow machining process. Performing experiments for all the parameters and accurately obtaining an optimized parameter in a short time are difficult to accomplish because the operation requires a precise finish. Computational fluid dynamics (CFD) simulation was employed to accurately determine optimum parameters. In the current work, a 2D model was designed, and the flow analysis, force calculation, and material removal prediction were performed and compared with the available experimental data. Another 3D model for a swaging die finishing using AFM was simulated at different viscosities of the media to study the effects on the controlling parameters. A CFD simulation was performed by using commercially available ANSYS FLUENT. Two phases were considered for the flow analysis, and multiphase mixture model was taken into account. The fluid was considered to be a

  2. Prominence fine-structure dynamics as inferred from 2D non-LTE models

    NASA Astrophysics Data System (ADS)

    Gunar, Stanislav; Schmieder, Brigitte; Mein, Pierre; Heinzel, Petr

    2012-07-01

    2D multi-thread prominence fine structure models are able to produce synthetic Lyman spectra in very good agreement with spectral observations by SOHO/SUMER including the spectral line asymmetries. The synthetic differential emission measure curves derived from these models are also in a good agreement with observations. Now we show that these models are also able to produce synthetic H-alpha line profiles in very good agreement with observations which allows us to analyze not only the physical parameters of the prominence fine-structure plasma but also some aspects of its dynamical behaviour. We compare the synthetic H-alpha spectra with the observed spectra of the April 26, 2007 prominence using three statistical parameters: the line integrated intensity, the line full-width at the half-maximum (FWHM), and the Doppler velocity derived from shifts of the line profiles. This statistical analysis allows us to conclude that the overall statistical distribution of the LOS velocities in the April 26, 2007 prominence at the time of the observations was below +/-15 km/s and in the prominence core was close to +/-10 km/s. In combination with the analysis of the Lyman spectra we determine several physical parameters of the observed prominence fine structures which show that the April 26, 2007 prominence was relatively less massive. We are also able to put some constrains on the prominence core temperature that might be relatively low, reaching values below 6000 K.

  3. 2d Affine XY-Spin Model/4d Gauge Theory Duality and Deconfinement

    SciTech Connect

    Anber, Mohamed M.; Poppitz, Erich; Unsal, Mithat; /SLAC /Stanford U., Phys. Dept. /San Francisco State U.

    2012-08-16

    We introduce a duality between two-dimensional XY-spin models with symmetry-breaking perturbations and certain four-dimensional SU(2) and SU(2) = Z{sub 2} gauge theories, compactified on a small spatial circle R{sup 1,2} x S{sup 1}, and considered at temperatures near the deconfinement transition. In a Euclidean set up, the theory is defined on R{sup 2} x T{sup 2}. Similarly, thermal gauge theories of higher rank are dual to new families of 'affine' XY-spin models with perturbations. For rank two, these are related to models used to describe the melting of a 2d crystal with a triangular lattice. The connection is made through a multi-component electric-magnetic Coulomb gas representation for both systems. Perturbations in the spin system map to topological defects in the gauge theory, such as monopole-instantons or magnetic bions, and the vortices in the spin system map to the electrically charged W-bosons in field theory (or vice versa, depending on the duality frame). The duality permits one to use the two-dimensional technology of spin systems to study the thermal deconfinement and discrete chiral transitions in four-dimensional SU(N{sub c}) gauge theories with n{sub f} {ge} 1 adjoint Weyl fermions.

  4. Field-induced magnetization jumps and quantum criticality in the 2D J-Q model

    NASA Astrophysics Data System (ADS)

    Iaizzi, Adam; Sandvik, Anders

    The J-Q model is a `designer hamiltonian' formed by adding a four spin `Q' term to the standard antiferromagnetic S = 1 / 2 Heisenberg model. The Q term drives a quantum phase transition to a valence-bond solid (VBS) state: a non-magnetic state with a pattern of local singlets which breaks lattice symmetries. The elementary excitations of the VBS are triplons, i.e. gapped S=1 quasiparticles. There is considerable interest in the quantum phase transition between the Néel and VBS states as an example of deconfined quantum criticality. Near the phase boundary, triplons deconfine into pairs of bosonic spin-1/2 excitations known as spinons. Using exact diagonalization and the stochastic series expansion quantum monte carlo method, we study the 2D J-Q model in the presence of an external magnetic field. We use the field to force a nonzero density of magnetic excitations at T=0 and look for signatures of Bose-Einstein condensation of spinons. At higher magnetic fields, there is a jump in the induced magnetization caused by the onset of an effective attractive interaction between magnons on a ferromagnetic background. We characterize the first order quantum phase transition and determine the minimum value of the coupling ratio q ≡ Q / J required to produce this jump. Funded by NSF DMR-1410126.

  5. Automatic 3D high-fidelity traffic interchange modeling using 2D road GIS data

    NASA Astrophysics Data System (ADS)

    Wang, Jie; Shen, Yuzhong

    2011-03-01

    3D road models are widely used in many computer applications such as racing games and driving simulations. However, almost all high-fidelity 3D road models were generated manually by professional artists at the expense of intensive labor. There are very few existing methods for automatically generating 3D high-fidelity road networks, especially for those existing in the real world. Real road network contains various elements such as road segments, road intersections and traffic interchanges. Among them, traffic interchanges present the most challenges to model due to their complexity and the lack of height information (vertical position) of traffic interchanges in existing road GIS data. This paper proposes a novel approach that can automatically produce 3D high-fidelity road network models, including traffic interchange models, from real 2D road GIS data that mainly contain road centerline information. The proposed method consists of several steps. The raw road GIS data are first preprocessed to extract road network topology, merge redundant links, and classify road types. Then overlapped points in the interchanges are detected and their elevations are determined based on a set of level estimation rules. Parametric representations of the road centerlines are then generated through link segmentation and fitting, and they have the advantages of arbitrary levels of detail with reduced memory usage. Finally a set of civil engineering rules for road design (e.g., cross slope, superelevation) are selected and used to generate realistic road surfaces. In addition to traffic interchange modeling, the proposed method also applies to other more general road elements. Preliminary results show that the proposed method is highly effective and useful in many applications.

  6. Generalized Mechanistic Model for the Chemical Vapor Deposition of 2D Transition Metal Dichalcogenide Monolayers.

    PubMed

    Govind Rajan, Ananth; Warner, Jamie H; Blankschtein, Daniel; Strano, Michael S

    2016-04-26

    Transition metal dichalcogenides (TMDs) like molybdenum disulfide (MoS2) and tungsten disulfide (WS2) are layered materials capable of growth to one monolayer thickness via chemical vapor deposition (CVD). Such CVD methods, while powerful, are notoriously difficult to extend across different reactor types and conditions, with subtle variations often confounding reproducibility, particularly for 2D TMD growth. In this work, we formulate the first generalized TMD synthetic theory by constructing a thermodynamic and kinetic growth mechanism linked to CVD reactor parameters that is predictive of specific geometric shape, size, and aspect ratio from triangular to hexagonal growth, depending on specific CVD reactor conditions. We validate our model using experimental data from Wang et al. (Chem. Mater. 2014, 26, 6371-6379) that demonstrate the systemic evolution of MoS2 morphology down the length of a flow CVD reactor where variations in gas phase concentrations can be accurately estimated using a transport model (CSulfur = 9-965 μmol/m(3); CMoO3 = 15-16 mmol/m(3)) under otherwise isothermal conditions (700 °C). A stochastic model which utilizes a site-dependent activation energy barrier based on the intrinsic TMD bond energies and a series of Evans-Polanyi relations leads to remarkable, quantitative agreement with both shape and size evolution along the reactor. The model is shown to extend to the growth of WS2 at 800 °C and MoS2 under varied process conditions. Finally, a simplified theory is developed to translate the model into a "kinetic phase diagram" of the growth process. The predictive capability of this model and its extension to other TMD systems promise to significantly increase the controlled synthesis of such materials. PMID:26937889

  7. 2D dynamical magma propagation modeling: application to the 2001 Mount Etna eruption

    NASA Astrophysics Data System (ADS)

    Pinel, Virginie; Carrara, Alexandre; Maccaferri, Francesco; Rivalta, Eleonora; Corbi, Fabio

    2016-04-01

    Numerical and analog studies of dike propagation in a stress field induced by volcanic edifice construction have shown that surface loading tends both to attract the magma and to reduce its velocity. Available numerical models can either calculate the trajectory or the velocity of the ascending dikes, but not both of them simultaneously. We developed a hybrid model of dyke propagation in two dimensions solving both for the magma trajectory and velocity as a function of the source overpressure, the magma physical properties (density and viscosity) as well as the crustal density and stress field. We first calculate a dyke trajectory in 2D and secondly run a 1D dynamical model of dyke propagation along this trajectory taken into account the influence of the stress field seen by the magma along this path. This model is used to characterize the influence of surface load on magma migration towards the surface and compared to previous results obtained by analog modeling.We find that the amplitude of dyke deflection and magma velocity variation depend on the ratio between the dyke driving pressure (source overpressure as well buoyancy) and the stress field perturbation. Our model is then applied to the July 2001 eruption of Etna, where the final dyke deflection had been previously interpreted as due to the topographic load by Bonaccorso et al. [2010]. We show that the velocity decrease observed during the last stage of the propagation can also be attributed to the local stress field. We use the dyke propagation duration to estimate the magma overpressure at the dyke bottom to be less than 4 MPa.

  8. Effects of CYP2D6 Status on Harmaline Metabolism, Pharmacokinetics and Pharmacodynamics, and a Pharmacogenetics-Based Pharmacokinetic Model

    PubMed Central

    Wu, Chao; Jiang, Xi-Ling; Shen, Hong-Wu; Yu, Ai-Ming

    2009-01-01

    Harmaline is a β-carboline alkaloid showing neuroprotective and neurotoxic properties. Our recent studies have revealed an important role for cytochrome P450 2D6 (CYP2D6) in harmaline O-demethylation. This study, therefore, aimed to delineate the effects of CYP2D6 phenotype/genotype on harmaline metabolism, pharmacokinetics (PK) and pharmacodynamics (PD), and to develop a pharmacogenetics mechanism-based compartmental PK model. In vitro kinetic studies on metabolite formation in human CYP2D6 extensive metabolizer (EM) and poor metabolizer (PM) hepatocytes indicated that harmaline O-demethylase activity (Vmax/Km) was about 9-fold higher in EM hepatocytes. Substrate depletion showed mono-exponential decay trait, and estimated in vitro harmaline clearance (CLint, μL/min/106 cells) was significantly lower in PM hepatocytes (28.5) than EM hepatocytes (71.1). In vivo studies in CYP2D6-humanized and wild-type mouse models showed that wild-type mice were subjected to higher and longer exposure to harmaline (5 and 15 mg/kg; i.v. and i.p.), and more severe hypothermic responses. The PK/PD data were nicely described by our pharmacogenetics-based PK model involving the clearance of drug by CYP2D6 (CLCYP2D6) and other mechanisms (CLother), and an indirect response PD model, respectively. Wild-type mice were also more sensitive to harmaline in marble-burying tests, as manifested by significantly lower ED50 and steeper Hill slope. These findings suggest that distinct CYP2D6 status may cause considerable variations in harmaline metabolism, PK and PD. In addition, the pharmacogenetics-based PK model may be extended to define PK difference caused by other polymorphic drug-metabolizing enzyme in different populations. PMID:19445902

  9. On the assimilation of SWOT type data into 2D shallow-water models

    NASA Astrophysics Data System (ADS)

    Frédéric, Couderc; Denis, Dartus; Pierre-André, Garambois; Ronan, Madec; Jérôme, Monnier; Jean-Paul, Villa

    2013-04-01

    In river hydraulics, assimilation of water level measurements at gauging stations is well controlled, while assimilation of images is still delicate. In the present talk, we address the richness of satellite mapped information to constrain a 2D shallow-water model, but also related difficulties. 2D shallow models may be necessary for small scale modelling in particular for low-water and flood plain flows. Since in both cases, the dynamics of the wet-dry front is essential, one has to elaborate robust and accurate solvers. In this contribution we introduce robust second order, stable finite volume scheme [CoMaMoViDaLa]. Comparisons of real like tests cases with more classical solvers highlight the importance of an accurate flood plain modelling. A preliminary inverse study is presented in a flood plain flow case, [LaMo] [HoLaMoPu]. As a first step, a 0th order data processing model improves observation operator and produces more reliable water level derived from rough measurements [PuRa]. Then, both model and flow behaviours can be better understood thanks to variational sensitivities based on a gradient computation and adjoint equations. It can reveal several difficulties that a model designer has to tackle. Next, a 4D-Var data assimilation algorithm used with spatialized data leads to improved model calibration and potentially leads to identify river discharges. All the algorithms are implemented into DassFlow software (Fortran, MPI, adjoint) [Da]. All these results and experiments (accurate wet-dry front dynamics, sensitivities analysis, identification of discharges and calibration of model) are currently performed in view to use data from the future SWOT mission. [CoMaMoViDaLa] F. Couderc, R. Madec, J. Monnier, J.-P. Vila, D. Dartus, K. Larnier. "Sensitivity analysis and variational data assimilation for geophysical shallow water flows". Submitted. [Da] DassFlow - Data Assimilation for Free Surface Flows. Computational software http

  10. MAST-2D diffusive model for flood prediction on domains with triangular Delaunay unstructured meshes

    NASA Astrophysics Data System (ADS)

    Aricò, C.; Sinagra, M.; Begnudelli, L.; Tucciarelli, T.

    2011-11-01

    A new methodology for the solution of the 2D diffusive shallow water equations over Delaunay unstructured triangular meshes is presented. Before developing the new algorithm, the following question is addressed: it is worth developing and using a simplified shallow water model, when well established algorithms for the solution of the complete one do exist? The governing Partial Differential Equations are discretized using a procedure similar to the linear conforming Finite Element Galerkin scheme, with a different flux formulation and a special flux treatment that requires Delaunay triangulation but entire solution monotonicity. A simple mesh adjustment is suggested, that attains the Delaunay condition for all the triangle sides without changing the original nodes location and also maintains the internal boundaries. The original governing system is solved applying a fractional time step procedure, that solves consecutively a convective prediction system and a diffusive correction system. The non linear components of the problem are concentrated in the prediction step, while the correction step leads to the solution of a linear system of the order of the number of computational cells. A semi-analytical procedure is applied for the solution of the prediction step. The discretized formulation of the governing equations allows to handle also wetting and drying processes without any additional specific treatment. Local energy dissipations, mainly the effect of vertical walls and hydraulic jumps, can be easily included in the model. Several numerical experiments have been carried out in order to test (1) the stability of the proposed model with regard to the size of the Courant number and to the mesh irregularity, (2) its computational performance, (3) the convergence order by means of mesh refinement. The model results are also compared with the results obtained by a fully dynamic model. Finally, the application to a real field case with a Venturi channel is presented.

  11. A computational model of the short-cut rule for 2D shape decomposition.

    PubMed

    Luo, Lei; Shen, Chunhua; Liu, Xinwang; Zhang, Chunyuan

    2015-01-01

    We propose a new 2D shape decomposition method based on the short-cut rule. The short-cut rule originates from cognition research, and states that the human visual system prefers to partition an object into parts using the shortest possible cuts. We propose and implement a computational model for the short-cut rule and apply it to the problem of shape decomposition. The model we proposed generates a set of cut hypotheses passing through the points on the silhouette, which represent the negative minima of curvature. We then show that most part-cut hypotheses can be eliminated by analysis of local properties of each. Finally, the remaining hypotheses are evaluated in ascending length order, which guarantees that of any pair of conflicting cuts only the shortest will be accepted. We demonstrate that, compared with state-of-the-art shape decomposition methods, the proposed approach achieves decomposition results, which better correspond to human intuition as revealed in psychological experiments. PMID:25438318

  12. Learning structured models for segmentation of 2-D and 3-D imagery.

    PubMed

    Lucchi, Aurelien; Marquez-Neila, Pablo; Becker, Carlos; Li, Yunpeng; Smith, Kevin; Knott, Graham; Fua, Pascal

    2015-05-01

    Efficient and accurate segmentation of cellular structures in microscopic data is an essential task in medical imaging. Many state-of-the-art approaches to image segmentation use structured models whose parameters must be carefully chosen for optimal performance. A popular choice is to learn them using a large-margin framework and more specifically structured support vector machines (SSVM). Although SSVMs are appealing, they suffer from certain limitations. First, they are restricted in practice to linear kernels because the more powerful nonlinear kernels cause the learning to become prohibitively expensive. Second, they require iteratively finding the most violated constraints, which is often intractable for the loopy graphical models used in image segmentation. This requires approximation that can lead to reduced quality of learning. In this paper, we propose three novel techniques to overcome these limitations. We first introduce a method to "kernelize" the features so that a linear SSVM framework can leverage the power of nonlinear kernels without incurring much additional computational cost. Moreover, we employ a working set of constraints to increase the reliability of approximate subgradient methods and introduce a new way to select a suitable step size at each iteration. We demonstrate the strength of our approach on both 2-D and 3-D electron microscopic (EM) image data and show consistent performance improvement over state-of-the-art approaches. PMID:25438309

  13. Basic Brackets of a 2D Model for the Hodge Theory Without its Canonical Conjugate Momenta

    NASA Astrophysics Data System (ADS)

    Kumar, R.; Gupta, S.; Malik, R. P.

    2016-06-01

    We deduce the canonical brackets for a two (1+1)-dimensional (2D) free Abelian 1-form gauge theory by exploiting the beauty and strength of the continuous symmetries of a Becchi-Rouet-Stora-Tyutin (BRST) invariant Lagrangian density that respects, in totality, six continuous symmetries. These symmetries entail upon this model to become a field theoretic example of Hodge theory. Taken together, these symmetries enforce the existence of exactly the same canonical brackets amongst the creation and annihilation operators that are found to exist within the standard canonical quantization scheme. These creation and annihilation operators appear in the normal mode expansion of the basic fields of this theory. In other words, we provide an alternative to the canonical method of quantization for our present model of Hodge theory where the continuous internal symmetries play a decisive role. We conjecture that our method of quantization is valid for a class of field theories that are tractable physical examples for the Hodge theory. This statement is true in any arbitrary dimension of spacetime.

  14. Optimal implicit 2-D finite differences to model wave propagation in poroelastic media

    NASA Astrophysics Data System (ADS)

    Itzá, Reymundo; Iturrarán-Viveros, Ursula; Parra, Jorge O.

    2016-05-01

    Numerical modeling of seismic waves in heterogeneous porous reservoir rocks is an important tool for the interpretation of seismic surveys in reservoir engineering. We apply globally optimal implicit staggered-grid finite-differences to model 2-D wave propagation in heterogeneous poroelastic media at a low-frequency range (<10kHz). We validate the numerical solution by comparing it to an analytical-transient solution obtaining clear seismic wavefields including fast P, slow P and S waves (for a porous media saturated with fluid). The numerical dispersion and stability conditions are derived using von Neumann analysis, showing that over a wide range of porous materials the Courant condition governs the stability and this optimal implicit scheme improves the stability of explicit schemes. High order explicit finite-differences (FD) can be replaced by some lower order optimal implicit FD so computational cost will not be as expensive while maintaining the accuracy. Here we compute weights for the optimal implicit FD scheme to attain an accuracy of γ = 10-8. The implicit spatial differentiation involves solving tridiagonal linear systems of equations through Thomas' algorithm.

  15. Optimal implicit 2-D finite differences to model wave propagation in poroelastic media

    NASA Astrophysics Data System (ADS)

    Itzá, Reymundo; Iturrarán-Viveros, Ursula; Parra, Jorge O.

    2016-08-01

    Numerical modeling of seismic waves in heterogeneous porous reservoir rocks is an important tool for the interpretation of seismic surveys in reservoir engineering. We apply globally optimal implicit staggered-grid finite differences (FD) to model 2-D wave propagation in heterogeneous poroelastic media at a low-frequency range (<10 kHz). We validate the numerical solution by comparing it to an analytical-transient solution obtaining clear seismic wavefields including fast P and slow P and S waves (for a porous media saturated with fluid). The numerical dispersion and stability conditions are derived using von Neumann analysis, showing that over a wide range of porous materials the Courant condition governs the stability and this optimal implicit scheme improves the stability of explicit schemes. High-order explicit FD can be replaced by some lower order optimal implicit FD so computational cost will not be as expensive while maintaining the accuracy. Here, we compute weights for the optimal implicit FD scheme to attain an accuracy of γ = 10-8. The implicit spatial differentiation involves solving tridiagonal linear systems of equations through Thomas' algorithm.

  16. T-duality without isometry via extended gauge symmetries of 2D sigma models

    NASA Astrophysics Data System (ADS)

    Chatzistavrakidis, Athanasios; Deser, Andreas; Jonke, Larisa

    2016-01-01

    Target space duality is one of the most profound properties of string theory. However it customarily requires that the background fields satisfy certain invariance conditions in order to perform it consistently; for instance the vector fields along the directions that T-duality is performed have to generate isometries. In the present paper we examine in detail the possibility to perform T-duality along non-isometric directions. In particular, based on a recent work of Kotov and Strobl, we study gauged 2D sigma models where gauge invariance for an extended set of gauge transformations imposes weaker constraints than in the standard case, notably the corresponding vector fields are not Killing. This formulation enables us to follow a procedure analogous to the derivation of the Buscher rules and obtain two dual models, by integrating out once the Lagrange multipliers and once the gauge fields. We show that this construction indeed works in non-trivial cases by examining an explicit class of examples based on step 2 nilmanifolds.

  17. Hydraulic Modeling of Alluvial Fans along the Truckee Canal using the 2-Dimensional Model SRH2D

    NASA Astrophysics Data System (ADS)

    Wright, J.; Kallio, R.; Sankovich, V.

    2013-12-01

    Alluvial fans are gently sloping, fan-shaped landforms created by sediment deposition at the ends of mountain valleys. Their gentle slopes and scenic vistas are attractive to developers. Unfortunately, alluvial fans are highly flood-prone, and the flow paths of flood events are highly variable, thereby placing human developments at risk. Many studies have been performed on alluvial fans in the arid west because of the uncertainty of their flow paths and flood extents. Most of these studies have been focused on flood elevations and mitigation. This study is not focused on the flood elevations. Rather, it is focused on the attenuation effects of alluvial fans on floods entering and potentially failing a Reclamation canal. The Truckee Canal diverts water from the Truckee River to Lahontan Reservoir. The drainage areas along the canal are alluvial fans with complex distributary channel networks . Ideally, in nature, the sediment grain-size distribution along the alluvial fan flow paths would provide enough infiltration and subsurface storage to attenuate floods entering the canal and reduce risk to low levels. Human development, however, can prevent the natural losses from occurring due to concentrated flows within the alluvial fan. While the concentrated flows might mitigate flood risk inside the fan, they do not lower the flood risk of the canal. A 2-dimensional hydraulic model, SRH-2D, was coupled to a 1-dimensional rainfall-runoff model to estimate the flood attenuation effects of the alluvial fan network surrounding an 11 mile stretch of the Truckee Canal near Fernley, Nevada. Floods having annual exceedance probabilities ranging from 1/10 to 1/100 were computed and analyzed. SRH-2D uses a zonal approach for modeling river systems, allowing areas to be divided into separate zones based on physical parameters such as surface roughness and infiltration. One of the major features of SRH-2D is the adoption of an unstructured hybrid mixed element mesh, which is based

  18. Combining multitracing and 2D-modelling to identify the dynamic of heavy metals during flooding.

    NASA Astrophysics Data System (ADS)

    Hissler, C.; Hostache, R.; Matgen, P.; Tosheva, Z.; David, E.; Bates, P.; Stille, P.

    2012-04-01

    to assess the risk of floodplain contamination in heavy metal due to river sediment deposition and to heavy metal partitioning between particulate and dissolved phases. We focus on a multidisciplinary approach combining environmental geochemistry (multitracing) and hydraulic modelling (using TELEMAC-2D). One important single flood event was selected to illustrate this innovative approach. During the entire flood, the river water was sampled every hour in order to collect the particulate and the dissolved fractions. All the tracers were analyzed in both fractions. An important set of hydrological and sedimentological data are used to reach a more efficient calibration of the TELEMAC modelling system. In addition to standard techniques of hydrochemistry, new approaches of in situ suspended sediment transport monitoring will help getting new insights on the hydraulic system behaviour.

  19. Water cycling beneath subduction zones in 2D and 3D numerical models (Invited)

    NASA Astrophysics Data System (ADS)

    Rupke, L.; Iyer, K. H.; Hasenclever, J.; Morgan, J.

    2013-12-01

    Tracing the cycling of fluids and volatiles through subduction zones continues to be a challenging task with budgets still having large error bars attached to them. In this contribution we show how numerical models can help to integrate various geological, geophysical, and geochemical datasets and how they can be used to put better bounds on the likely amounts of water being subducted, released into the arc and back-arc melting regions, and recycled to the deeper mantle. To achieve this task we use a suite of numerical models. Bending related faulting and hydration of the incoming lithosphere is resolved using a reactive flow model that solves for crustal scale fluid flow and mantle serpentinization using reaction kinetics. Seismic tomography studies from offshore Chile and Central America are used to evaluate and constrain the effective reaction rate. These rates are then used to assess the contribution of serpentinization to the water budget at subduction zones. The pattern of hydration is controlled by the reaction kinetics and serpentinization is most intense around the 270°C isotherm. The depth of this isotherm correlates well with the dominant spacing of double seismic zones observed globally. Comparison of the results with heat flow data suggests that observed seafloor temperature gradients in the bend-fault region are too low to be caused by ';one-pass' downward water flow into the serpentinizing lithosphere, but rather suggest that bend-faults are areas of active hydrothermal circulation. This implies that serpentine-sourced vents and chemosynthetic vent communities should be found in this deep-sea environment as well. Dehydration reactions are resolved with a 2D kinematic subduction zone model that computes the temperature field and the likely locations and volumes of slab fluid release due to metamorphic dehydration reactions. Here we find that up to 1/3 of the subducted water may be transported into the deeper mantle for the coldest subduction zones

  20. Modelling 2001 lahars at Popocatépetl volcano using FLO2D numerical code

    NASA Astrophysics Data System (ADS)

    Caballero, L.; Capra, L.

    2013-12-01

    Popocatépetl volcano is located on the central part of the Transmexican Volcanic Belt. It is one of the most active volcanoes in Mexico and endanger more than 25 million people that lives in its surroundings. In the last months, the renewal of its volcanic activity put into alert scientific community. One of the possible scenarios is the 2001 explosive activity, which was characterized by a 8 km eruptive column and the subsequent formation of pumice flows up to 4 km from the crater. Lahars were generated few hours after, remobilizing the new deposits towards NE flank of the volcano, along Huiloac Gorge, almost reaching Santiago Xalitzintla town (Capra et al., 2004). The occurrence of a similar scenario makes very important to reproduce this event to delimitate accurately lahar hazard zones. In this work, 2001 lahar deposit is modeled using FLO2D numerical code. Geophone data is used to reconstruct initial hydrograph and sediment concentration. Sensitivity study of most important parameters used by this code like Manning, and α and β coefficients was conducted in order to achieve a good simulation. Results obtained were compared with field data and demonstrated a good agreement in thickness and flow distribution. A comparison with previously published data with laharZ program (Muñoz-Salinas, 2009) is also made. Additionally, lahars with fluctuating sediment concentrations but with similar volume are simulated to observe the influence of the rheological behavior on lahar distribution.

  1. Modeling of two-storey precast school building using Ruaumoko 2D program

    SciTech Connect

    Hamid, N. H.; Tarmizi, L. H.; Ghani, K. D.

    2015-05-15

    The long-distant earthquake loading from Sumatra and Java Island had caused some slight damages to precast and reinforced concrete buildings in West Malaysia such as cracks on wall panels, columns and beams. Subsequently, the safety of existing precast concrete building is needed to be analyzed because these buildings were designed using BS 8110 which did not include the seismic loading in the design. Thus, this paper emphasizes on the seismic performance and dynamic behavior of precast school building constructed in Malaysia under three selected past earthquakes excitations ; El Centro 1940 North-South, El Centro East-West components and San Fernando 1971 using RUAUMOKO 2D program. This program is fully utilized by using prototype precast school model and dynamic non-linear time history analysis. From the results, it can be concluded that two-storey precast school building has experienced severe damage and partial collapse especially at beam-column joint under San Fernando and El Centro North-South Earthquake as its exceeds the allowable inter-storey drift and displacement as specified in Eurocode 8. The San Fernando earthquake has produced a massive destruction to the precast building under viscous damping, ξ = 5% and this building has generated maximum building displacement of 435mm, maximum building drift of 0.68% and maximum bending moment at 8458kNm.

  2. Uncertainty in 2D hydrodynamic models from errors in roughness parameterization based on aerial images

    NASA Astrophysics Data System (ADS)

    Straatsma, Menno; Huthoff, Fredrik

    2011-01-01

    In The Netherlands, 2D-hydrodynamic simulations are used to evaluate the effect of potential safety measures against river floods. In the investigated scenarios, the floodplains are completely inundated, thus requiring realistic representations of hydraulic roughness of floodplain vegetation. The current study aims at providing better insight into the uncertainty of flood water levels due to uncertain floodplain roughness parameterization. The study focuses on three key elements in the uncertainty of floodplain roughness: (1) classification error of the landcover map, (2), within class variation of vegetation structural characteristics, and (3) mapping scale. To assess the effect of the first error source, new realizations of ecotope maps were made based on the current floodplain ecotope map and an error matrix of the classification. For the second error source, field measurements of vegetation structure were used to obtain uncertainty ranges for each vegetation structural type. The scale error was investigated by reassigning roughness codes on a smaller spatial scale. It is shown that classification accuracy of 69% leads to an uncertainty range of predicted water levels in the order of decimeters. The other error sources are less relevant. The quantification of the uncertainty in water levels can help to make better decisions on suitable flood protection measures. Moreover, the relation between uncertain floodplain roughness and the error bands in water levels may serve as a guideline for the desired accuracy of floodplain characteristics in hydrodynamic models.

  3. Modeling of two-storey precast school building using Ruaumoko 2D program

    NASA Astrophysics Data System (ADS)

    Hamid, N. H.; Tarmizi, L. H.; Ghani, K. D.

    2015-05-01

    The long-distant earthquake loading from Sumatra and Java Island had caused some slight damages to precast and reinforced concrete buildings in West Malaysia such as cracks on wall panels, columns and beams. Subsequently, the safety of existing precast concrete building is needed to be analyzed because these buildings were designed using BS 8110 which did not include the seismic loading in the design. Thus, this paper emphasizes on the seismic performance and dynamic behavior of precast school building constructed in Malaysia under three selected past earthquakes excitations ; El Centro 1940 North-South, El Centro East-West components and San Fernando 1971 using RUAUMOKO 2D program. This program is fully utilized by using prototype precast school model and dynamic non-linear time history analysis. From the results, it can be concluded that two-storey precast school building has experienced severe damage and partial collapse especially at beam-column joint under San Fernando and El Centro North-South Earthquake as its exceeds the allowable inter-storey drift and displacement as specified in Eurocode 8. The San Fernando earthquake has produced a massive destruction to the precast building under viscous damping, ξ = 5% and this building has generated maximum building displacement of 435mm, maximum building drift of 0.68% and maximum bending moment at 8458kNm.

  4. 2D and 3D multipactor modeling in dielectric-loaded accelerator structures

    NASA Astrophysics Data System (ADS)

    Sinitsyn, Oleksandr; Nusinovich, Gregory; Antonsen, Thomas

    2010-11-01

    Multipactor (MP) is known as the avalanche growth of the number of secondary electrons emitted from a solid surface exposed to an RF electric field under vacuum conditions. MP is a severe problem in modern rf systems and, therefore, theoretical and experimental studies of MP are of great interest to the researchers working in various areas of physics and engineering. In this work we present results of MP studies in dielectric-loaded accelerator (DLA) structures. First, we show simulation results obtained with the use of the 2D self-consistent MP model (O. V. Sinitsyn, et. al., Phys. Plasmas, vol. 16, 073102 (2009)) and compare those to experimental ones obtained during recent extensive studies of DLA structures performed by Argonne National Laboratory, Naval Research Laboratory, SLAC National Accelerator Laboratory and Euclid TechLabs (C. Jing, et al., IEEE Trans. Plasma Sci., vol. 38, pp. 1354-1360 (2010)). Then we present some new results of 3D analysis of MP which include studies of particle trajectories and studies of MP development at the early stage.

  5. Spectral functions in the 1D and 2D Bose Hubbard model

    NASA Astrophysics Data System (ADS)

    Ivancic, Robert; Duchon, Eric; Trivedi, Nandini

    2014-03-01

    We use state of the art numerical techniques including quantum Monte Carlo and maximum entropy methods to obtain the low energy excitation spectra in the superfluid and Mott-insulator phases of the Bose Hubbard model. These results are checked in 1D against Bethe Ansatz and tDMRG results and extended to 2D where such approaches are impossible. In the superfluid, we find linearly dispersing Bogoliubov sound modes as well as additional gapped modes broadened by interaction effects. In the Mott insulator, we find evidence for a finite gap and well defined quasiparticle excitations. We examine properties such as the excitation lifetime, density of states, and speed of sound as the system is tuned across the quantum phase transition that separates the superfluid and Mott states. These results provide an important theoretical framework for upcoming ultracold atom experiments in one and two dimensions. We acknowledge support from the NSF DMR-0907275 (R.I., E.D. and N.T.).

  6. GPU computing with OpenCL to model 2D elastic wave propagation: exploring memory usage

    NASA Astrophysics Data System (ADS)

    Iturrarán-Viveros, Ursula; Molero-Armenta, Miguel

    2015-01-01

    Graphics processing units (GPUs) have become increasingly powerful in recent years. Programs exploring the advantages of this architecture could achieve large performance gains and this is the aim of new initiatives in high performance computing. The objective of this work is to develop an efficient tool to model 2D elastic wave propagation on parallel computing devices. To this end, we implement the elastodynamic finite integration technique, using the industry open standard open computing language (OpenCL) for cross-platform, parallel programming of modern processors, and an open-source toolkit called [Py]OpenCL. The code written with [Py]OpenCL can run on a wide variety of platforms; it can be used on AMD or NVIDIA GPUs as well as classical multicore CPUs, adapting to the underlying architecture. Our main contribution is its implementation with local and global memory and the performance analysis using five different computing devices (including Kepler, one of the fastest and most efficient high performance computing technologies) with various operating systems.

  7. 1D and 2D urban dam-break flood modelling in Istanbul, Turkey

    NASA Astrophysics Data System (ADS)

    Ozdemir, Hasan; Neal, Jeffrey; Bates, Paul; Döker, Fatih

    2014-05-01

    Urban flood events are increasing in frequency and severity as a consequence of several factors such as reduced infiltration capacities due to continued watershed development, increased construction in flood prone areas due to population growth, the possible amplification of rainfall intensity due to climate change, sea level rise which threatens coastal development, and poorly engineered flood control infrastructure (Gallegos et al., 2009). These factors will contribute to increased urban flood risk in the future, and as a result improved modelling of urban flooding according to different causative factor has been identified as a research priority (Gallegos et al., 2009; Ozdemir et al. 2013). The flooding disaster caused by dam failures is always a threat against lives and properties especially in urban environments. Therefore, the prediction of dynamics of dam-break flows plays a vital role in the forecast and evaluation of flooding disasters, and is of long-standing interest for researchers. Flooding occurred on the Ayamama River (Istanbul-Turkey) due to high intensity rainfall and dam-breaching of Ata Pond in 9th September 2009. The settlements, industrial areas and transportation system on the floodplain of the Ayamama River were inundated. Therefore, 32 people were dead and millions of Euros economic loses were occurred. The aim of this study is 1 and 2-Dimensional flood modelling of the Ata Pond breaching using HEC-RAS and LISFLOOD-Roe models and comparison of the model results using the real flood extent. The HEC-RAS model solves the full 1-D Saint Venant equations for unsteady open channel flow whereas LISFLOOD-Roe is the 2-D shallow water model which calculates the flow according to the complete Saint Venant formulation (Villanueva and Wright, 2006; Neal et al., 2011). The model consists a shock capturing Godunov-type scheme based on the Roe Riemann solver (Roe, 1981). 3 m high resolution Digital Surface Model (DSM), natural characteristics of the pond

  8. D Recording for 2d Delivering - the Employment of 3d Models for Studies and Analyses -

    NASA Astrophysics Data System (ADS)

    Rizzi, A.; Baratti, G.; Jiménez, B.; Girardi, S.; Remondino, F.

    2011-09-01

    In the last years, thanks to the advances of surveying sensors and techniques, many heritage sites could be accurately replicated in digital form with very detailed and impressive results. The actual limits are mainly related to hardware capabilities, computation time and low performance of personal computer. Often, the produced models are not visible on a normal computer and the only solution to easily visualized them is offline using rendered videos. This kind of 3D representations is useful for digital conservation, divulgation purposes or virtual tourism where people can visit places otherwise closed for preservation or security reasons. But many more potentialities and possible applications are available using a 3D model. The problem is the ability to handle 3D data as without adequate knowledge this information is reduced to standard 2D data. This article presents some surveying and 3D modeling experiences within the APSAT project ("Ambiente e Paesaggi dei Siti d'Altura Trentini", i.e. Environment and Landscapes of Upland Sites in Trentino). APSAT is a multidisciplinary project funded by the Autonomous Province of Trento (Italy) with the aim documenting, surveying, studying, analysing and preserving mountainous and hill-top heritage sites located in the region. The project focuses on theoretical, methodological and technological aspects of the archaeological investigation of mountain landscape, considered as the product of sequences of settlements, parcelling-outs, communication networks, resources, and symbolic places. The mountain environment preserves better than others the traces of hunting and gathering, breeding, agricultural, metallurgical, symbolic activities characterised by different lengths and environmental impacts, from Prehistory to the Modern Period. Therefore the correct surveying and documentation of this heritage sites and material is very important. Within the project, the 3DOM unit of FBK is delivering all the surveying and 3D material to

  9. Approaches to Modeling Coupled Flow and Reaction in a 2-D Cementation Experiment

    SciTech Connect

    Steefel, Carl; Cochepin, B.; Trotignon, L.; Bildstein, O.; Steefel, C.; Lagneau, V.; van der Lee, J.

    2008-04-01

    Porosity evolution at reactive interfaces is a key process that governs the evolution and performances of many engineered systems that have important applications in earth and environmental sciences. This is the case, for example, at the interface between cement structures and clays in deep geological nuclear waste disposals. Although in a different transport regime, similar questions arise for permeable reactive barriers used for biogeochemical remediation in surface environments. The COMEDIE project aims at investigating the coupling between transport, hydrodynamics and chemistry when significant variations of porosity occur. The present work focuses on a numerical benchmark used as a design exercise for the future COMEDIE-2D experiment. The use of reactive transport simulation tools like Hytec and Crunch provides predictions of the physico-chemical evolutions that are expected during the future experiments in laboratory. Focus is given in this paper on the evolution during the simulated experiment of precipitate, permeability and porosity fields. A first case is considered in which the porosity is constant. Results obtained with Crunch and Hytec are in relatively good agreement. Differences are attributable to the models of reactive surface area taken into account for dissolution/precipitation processes. Crunch and Hytec simulations taking into account porosity variations are then presented and compared. Results given by the two codes are in qualitative agreement, with differences attributable in part to the models of reactive surface area for dissolution/precipitation processes. As a consequence, the localization of secondary precipitates predicted by Crunch leads to lower local porosities than for predictions obtained by Hytec and thus to a stronger coupling between flow and chemistry. This benchmark highlights the importance of the surface area model employed to describe systems in which strong porosity variations occur as a result of dissolution

  10. Simulation of Subgrid Orographic Convection and Precipitation with 2-D Cloud-Resolving Models Embedded in a GCM Grid

    NASA Astrophysics Data System (ADS)

    Jung, J.; Arakawa, A.

    2015-12-01

    Through explicitly resolved cloud-scale processes by embedded 2-D cloud-resolving models (CRMs), the Multiscale Modeling Framework (MMF) known as the superparameterization has been reasonably successful to simulate various atmospheric events over a wide range of time scales. One thing to be justified is, however, if the influence of complex 3-D topography can be adequately represented by the embedded 2-D CRMs. In this study, simulations are performed in the presence of a variety of topography with embedded 3-D and 2-D CRMs in a single-column inactive GCM. Through the comparison between these simulations, it is demonstrated that the 2-D representation of topography is able to simulate the statistics of precipitation due to 3-D topography reasonably well as long as the topographic characteristics, such as the mean and standard deviation, are closely recognized. It is also shown that the use of two perpendicular sets of 2-D representations tends to reduce the error due to a 2-D representation.