WorldWideScience

Sample records for 1d model atmospheres

  1. What causes the large extensions of red-supergiant atmospheres? Comparisons of interferometric observations with 1-D hydrostatic, 3-D convection, and 1-D pulsating model atmospheres

    Arroyo-Torres, B; Chiavassa, A; Scholz, M; Freytag, B; Marcaide, J M; Hauschildt, P H; Wood, P R; Abellan, F J

    2015-01-01

    We present the atmospheric structure and the fundamental parameters of three red supergiants, increasing the sample of RSGs observed by near-infrared spectro-interferometry. Additionally, we test possible mechanisms that may explain the large observed atmospheric extensions of RSGs. We carried out spectro-interferometric observations of 3 RSGs in the near-infrared K-band with the VLTI/AMBER instrument at medium spectral resolution. To comprehend the extended atmospheres, we compared our observational results to predictions by available hydrostatic PHOENIX, available 3-D convection, and new 1-D self-excited pulsation models of RSGs. Our near-infrared flux spectra are well reproduced by the PHOENIX model atmospheres. The continuum visibility values are consistent with a limb-darkened disk as predicted by the PHOENIX models, allowing us to determine the angular diameter and the fundamental parameters of our sources. Nonetheless, in the case of V602 Car and HD 95686, the PHOENIX model visibilities do not predict ...

  2. FORest Canopy Atmosphere Transfer (FORCAsT) 1.0: a 1-D model of biosphere-atmosphere chemical exchange

    Ashworth, K.; Chung, S. H.; Griffin, R. J.; Chen, J.; Forkel, R.; Bryan, A. M.; Steiner, A. L.

    2015-11-01

    Biosphere-atmosphere interactions play a critical role in governing atmospheric composition, mediating the concentrations of key species such as ozone and aerosol, thereby influencing air quality and climate. The exchange of reactive trace gases and their oxidation products (both gas and particle phase) is of particular importance in this process. The FORCAsT (FORest Canopy Atmosphere Transfer) 1-D model is developed to study the emission, deposition, chemistry and transport of volatile organic compounds (VOCs) and their oxidation products in the atmosphere within and above the forest canopy. We include an equilibrium partitioning scheme, making FORCAsT one of the few canopy models currently capable of simulating the formation of secondary organic aerosols (SOAs) from VOC oxidation in a forest environment. We evaluate the capability of FORCAsT to reproduce observed concentrations of key gas-phase species and report modeled SOA concentrations within and above a mixed forest at the University of Michigan Biological Station (UMBS) during the Community Atmosphere-Biosphere Interactions Experiment (CABINEX) field campaign in the summer of 2009. We examine the impact of two different gas-phase chemical mechanisms on modelled concentrations of short-lived primary emissions, such as isoprene and monoterpenes, and their oxidation products. While the two chemistry schemes perform similarly under high-NOx conditions, they diverge at the low levels of NOx at UMBS. We identify peroxy radical and alkyl nitrate chemistry as the key causes of the differences, highlighting the importance of this chemistry in understanding the fate of biogenic VOCs (bVOCs) for both the modelling and measurement communities.

  3. Assessing the habitability of planets with Earth-like atmospheres with 1D and 3D climate modeling

    Godolt, M.; Grenfell, J. L.; Kitzmann, D.; Kunze, M.; Langematz, U.; Patzer, A. B. C.; Rauer, H.; Stracke, B.

    2016-07-01

    Context. The habitable zone (HZ) describes the range of orbital distances around a star where the existence of liquid water on the surface of an Earth-like planet is in principle possible. The applicability of one-dimensional (1D) climate models for the estimation of the HZ boundaries has been questioned by recent three-dimensional (3D) climate studies. While 3D studies can calculate the water vapor, ice albedo, and cloud feedback self-consistently and therefore allow for a deeper understanding and the identification of relevant climate processes, 1D model studies rely on fewer model assumptions and can be more easily applied to the large parameter space possible for extrasolar planets. Aims: We evaluate the applicability of 1D climate models to estimate the potential habitability of Earth-like extrasolar planets by comparing our 1D model results to those of 3D climate studies in the literature. We vary the two important planetary properties, surface albedo and relative humidity, in the 1D model. These depend on climate feedbacks that are not treated self-consistently in most 1D models. Methods: We applied a cloud-free 1D radiative-convective climate model to calculate the climate of Earth-like planets around different types of main-sequence stars with varying surface albedo and relative humidity profile. We compared the results to those of 3D model calculations available in the literature and investigated to what extent the 1D model can approximate the surface temperatures calculated by the 3D models. Results: The 1D parameter study results in a large range of climates possible for an Earth-sized planet with an Earth-like atmosphere and water reservoir at a certain stellar insolation. At some stellar insolations the full spectrum of climate states could be realized, i.e., uninhabitable conditions due to surface temperatures that are too high or too low as well as habitable surface conditions, depending only on the relative humidity and surface albedo assumed. When

  4. 1-D Air-snowpack modeling of atmospheric nitrous acid at South Pole during ANTCI 2003

    W. Liao

    2008-12-01

    Full Text Available A 1-D air-snowpack model of HONO has been developed and constrained by observed chemistry and meteorology data. The 1-D model includes molecular diffusion and mechanical dispersion, windpumping in snow, gas phase to quasi-liquid layer phase HONO transfer and quasi-liquid layer nitrate and interstitial air HONO photolysis. Photolysis of nitrate is important as a dominant HONO source inside the snowpack, however, the observed HONO emission from the snowpack was triggered mainly by the equilibrium between quasi liquid layer nitrite and firn air HONO deep down the snow surface (i.e. 30 cm below snow surface. The high concentration of HONO in the firn air is subsequently transported above the snowpack by diffusion and windpumping. The model uncertainties come mainly from lack of measurements and the interpretation of the QLL properties based on the bulk snow measurements. One critical factor is the ionic strength of QLL nitrite, which is estimated here by the bulk snow pH, nitrite concentration, and QLL to bulk snow volume ratio.

  5. FORest canopy atmosphere transfer (FORCAsT 1.0: a 1-D model of biosphere–atmosphere chemical exchange

    K. Ashworth

    2015-07-01

    Full Text Available Biosphere-atmosphere interactions play a critical role in governing atmospheric composition, mediating the concentration of key species such as ozone and aerosol, thereby influencing air quality and climate. The exchange of reactive trace gases and their oxidation products (both gas and particle phase is of particular importance in this process. The FORCAsT (FORest Canopy AtmoSphere Transfer one-dimensional model is developed to study the emission, deposition, chemistry and transport of volatile organic compounds (VOCs and their oxidation products in the atmosphere within and above the forest canopy. We include an equilibrium partitioning scheme, making FORCAsT one of the few canopy models currently capable of simulating the formation of secondary organic aerosols (SOA from VOC oxidation in a forest environment. We evaluate the capability of FORCAsT to reproduce observed concentrations of key gas-phase species and report modeled SOA concentrations within and above a mixed forest at the University of Michigan Biological Station (UMBS during the Community Atmosphere-Biosphere Interactions Experiment (CABINEX field campaign in summer 2009. We examine the impact of two different gas-phase chemical mechanisms on modelled concentrations of short-lived primary emissions, such as isoprene and monoterpenes, and their oxidation products. While the two chemistry schemes perform similarly under high-NOx conditions, they diverge at the low levels of NOx at UMBS. We identify peroxy radical and alkyl nitrate chemistry as the key causes of the differences, highlighting the importance of this chemistry in understanding the fate of biogenic VOCs (bVOCs for both the modelling and measurement communities.

  6. Mg line formation in late-type stellar atmospheres: II. Calculations in a grid of 1D models

    Osorio, Yeisson

    2015-01-01

    Mg is the alpha element of choice for Galactic population and chemical evolution studies as it is easily detectable in all late-type stars. Such studies require precise elemental abundances, and thus departures from LTE need to be accounted for. Our goal is to provide reliable departure coefficients and equivalent widths in non-LTE, and for reference in LTE, for diagnostic lines of Mg studied in late-type stars. These can be used e.g., to correct LTE spectra and abundances. Using the model atom built and tested in the preceding paper in this series, we performed non-LTE radiative transfer calculations in a grid of 3945 stellar 1D atmospheric models. We used a sub-grid of 86 models to explore the propagation of errors in the recent atomic collision calculations to the radiative transfer results. We obtained departure coefficients for all the levels and equivalent widths (in LTE and non-LTE) for all the radiative transitions included in the "final" model atom of Osorio et al.. We present and describe our result...

  7. Assessing the habitability of planets with Earth-like atmospheres with 1D and 3D climate modeling

    Godolt, M; Kitzmann, D; Kunze, M; Langematz, U; Patzer, A B C; Rauer, H; Stracke, B

    2016-01-01

    The habitable zone (HZ) describes the range of orbital distances around a star where the existence of liquid water on the surface of an Earth-like planet is in principle possible. While 3D climate studies can calculate the water vapor, ice albedo, and cloud feedback self-consistently and therefore allow for a deeper understanding and the identification of relevant climate processes, 1D model studies rely on fewer model assumptions and can be more easily applied to the large parameter space possible for exoplanets. We evaluate the applicability of 1D climate models to estimate the potential habitability of Earth-like exoplanets by comparing our 1D model results to those of 3D climate studies in the literature. We applied a cloud-free 1D radiative-convective climate model to calculate the climate of Earth-like planets around different types of main-sequence stars with varying surface albedo and relative humidity profile. These parameters depend on climate feedbacks that are not treated self-consistently in most...

  8. Mg line formation in late-type stellar atmospheres. II. Calculations in a grid of 1D models

    Osorio, Y.; Barklem, P. S.

    2016-02-01

    Context. Mg is the α element of choice for Galactic population and chemical evolution studies because it is easily detectable in all late-type stars. Such studies require precise elemental abundances, and thus departures from local thermodynamic equilibrium (LTE) need to be accounted for. Aims: Our goal is to provide reliable departure coefficients and equivalent widths in non-LTE, and for reference in LTE, for diagnostic lines of Mg studied in late-type stars. These can be used, for example, to correct LTE spectra and abundances. Methods: Using the model atom built and tested in the preceding paper in this series, we performed non-LTE radiative transfer calculations in a grid of 3945 stellar 1D atmospheric models. We used a sub-grid of 86 models to explore the propagation of errors in the recent atomic collision calculations to the radiative transfer results. Results: We obtained departure coefficients for all the levels and equivalent widths (in LTE and non-LTE) for all the radiative transitions included in the "final" model atom presented in Paper I. Here we present and describe our results and show some examples of applications of the data. The errors that result from uncertainties in the collisional data are investigated and tabulated. The results for equivalent widths and departure coefficients are made freely available. Conclusions: Giants tend to have negative abundance corrections while dwarfs have positive, though small, corrections. Error analysis results show that uncertainties related to the atomic collision data are typically on the order of 0.01 dex or less, although for few stellar models in specific lines uncertainties can be as large as 0.03 dex. As these errors are less than or on the same order as typical corrections, we expect that we can use these results to extract Mg abundances from high-quality spectra more reliably than from classical LTE analysis. Full Table 1 is only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130

  9. NLTE line formation of Fe for late-type stars. I. Standard stars with 1D and <3D> model atmospheres

    Bergemann, Maria; Collet, Remo; Magic, Zazralt; Asplund, Martin

    2012-01-01

    We investigate departures from LTE in the line formation of Fe for a number of well-studied late-type stars in different evolutionary stages. A new model of Fe atom was constructed from the most up-to-date theoretical and experimental atomic data available so far. Non-local thermodynamic equilibrium (NLTE) line formation calculations for Fe were performed using 1D hydrostatic MARCS and MAFAGS-OS model atmospheres, as well as the spatial and temporal average stratifications from full 3D hydrodynamical simulations of stellar convection computed using the Stagger code. It is shown that the Fe I/Fe II ionization balance can be well established with the 1D and mean 3D models under NLTE including calibrated inelastic collisions with H I calculated from the Drawin's (1969) formulae. Strong low-excitation Fe I lines are very sensitive to the atmospheric structure; classical 1D models fail to provide consistent excitation balance, particularly so for cool metal-poor stars. A better agreement between Fe I lines spannin...

  10. Diagnostics from a 1-D atmospheric column

    Flatley, J.M.; Mace, G. [Pennsylvania State Univ., University Park, PA (United States)

    1996-04-01

    Various diagnostics were computed from an array of radiosondes during an intensive field operation arranged by the Atmospheric Radiation Measurement Program. The network data was centered around the site at Lamont, Oklahoma. The apparent heat source and apparent moisture sink were computed and compared to the kinematic vertical velocity for both real data and the mesoscale analysis and prediction system. Three different case studies of various weathe regimes were examined.

  11. Near-infrared spectro-interferometry of Mira variables and comparisons to 1D dynamic model atmospheres and 3D convection simulations

    Wittkowski, M; Freytag, B; Scholz, M; Hoefner, S; Karovicova, I; Whitelock, P A

    2016-01-01

    We obtained a total of 20 near-infrared K-band spectro-interferometric snapshot observations of the Mira variables o Cet, R Leo, R Aqr, X Hya, W Vel, and R Cnc with a spectral resolution of about 1500. We compared observed flux and visibility spectra with predictions by CODEX 1D dynamic model atmospheres and with azimuthally averaged intensities based on CO5BOLD 3D dynamic model atmospheres including convection. Our visibility data confirm the presence of spatially extended molecular atmospheres located above the continuum radii with large-scale inhomogeneities or clumps that contribute a few percent of the total flux. The detailed structure of the inhomogeneities or clumps show a variability on time scales of 3 months and above. Both modeling attempts provided satisfactory fits to our data. In particular, they are both consistent with the observed decrease in the visibility function at molecular bands of water vapor and CO, indicating a spatially extended molecular atmosphere. Observational variability phase...

  12. A 1D radiative-convective model of H2O-CO2 atmospheres around young telluric planets: an update

    Marcq, Emmanuel; Salvador, Arnaud; Massol, Hélène; Chassefière, Éric

    2016-04-01

    The study of the early phases of the evolution of terrestrial planets has recently known significant progress [1,2]. It appears that their cooling phase during the magma ocean stage is first dominated by a radiative cooling stage through its atmosphere. If the planet is able to reach radiative balance during this stage, then its further evolution is dominated by the escape flux, and no large scale condensation of water occurs (Hamano-type II planets). On the other hand, if the planet is far enough from the sun, then radiative equilibrium cannot be reached until the outgoing flux has fallen below the runaway greenhouse limit, implying the condensation of most atmospheric water vapor into a global water ocean, thus sheltering most water from atmospheric escape (Hamano-type I planet). In the solar system, Earth is clearly a type-I planet, whereas Venus was most likely a type-II planet from quite early on in its history [1,2]. In this presentation, we will deal with the atmospheric radiative model used by [2] and first described in [3]. After describing its recent improvements since [3] (pressure grid enabling an arbitrary total volatile amount, correction of the k-correlated radiative transfer in the thermal radiation, improvement of the numerical stability and integration scheme) and their consequences on the detectability of extrasolar type-I or type-II planets, we will deal with the possible improvements and extensions to such models, such as but not limited to: (1) adopting a 1D-spherical geometry suited for larger atmospheres around smaller planets, (2) improvement of the visible albedo parameterization based on recent 3D-modelling GCM [4]. [1] : K. Hamano et al., Nature (2013) [2] : T. Lebrun et al. JGR (2013) [3] : E. Marcq, JGR (2012) [4] : J. Leconte et al. (2015)

  13. YORP torques with 1D thermal model

    Breiter, Slawomir; Czekaj, Maria

    2010-01-01

    A numerical model of the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect for objects defined in terms of a triangular mesh is described. The algorithm requires that each surface triangle can be handled independently, which implies the use of a 1D thermal model. Insolation of each triangle is determined by an optimized ray-triangle intersection search. Surface temperature is modeled with a spectral approach; imposing a quasi-periodic solution we replace heat conduction equation by the Helmholtz equation. Nonlinear boundary conditions are handled by an iterative, FFT based solver. The results resolve the question of the YORP effect in rotation rate independence on conductivity within the nonlinear 1D thermal model regardless of the accuracy issues and homogeneity assumptions. A seasonal YORP effect in attitude is revealed for objects moving on elliptic orbits when a nonlinear thermal model is used.

  14. Forest-atmosphere BVOC exchange in diverse and structurally complex canopies: 1-D modeling of a mid-successional forest in northern Michigan

    Bryan, Alexander M.; Cheng, Susan J.; Ashworth, Kirsti; Guenther, Alex B.; Hardiman, Brady; Bohrer, Gil; Steiner, A. L.

    2015-11-01

    Foliar emissions of biogenic volatile organic compounds (BVOC)dimportant precursors of tropospheric ozone and secondary organic aerosolsdvary widely by vegetation type. Modeling studies to date typi-cally represent the canopy as a single dominant tree type or a blend of tree types, yet many forests are diverse with trees of varying height. To assess the sensitivity of biogenic emissions to tree height vari-ation, we compare two 1-D canopy model simulations in which BVOC emission potentials are homo-geneous or heterogeneous with canopy depth. The heterogeneous canopy emulates the mid-successional forest at the University of Michigan Biological Station (UMBS). In this case, high-isoprene-emitting fo-liage (e.g., aspen and oak) is constrained to the upper canopy, where higher sunlight availability increases the light-dependent isoprene emission, leading to 34% more isoprene and its oxidation products as compared to the homogeneous simulation. Isoprene declines from aspen mortality are 10% larger when heterogeneity is considered. Overall, our results highlight the importance of adequately representing complexities of forest canopy structure when simulating light-dependent BVOC emissions and chemistry.

  15. Forest-atmosphere BVOC exchange in diverse and structurally complex canopies: 1-D modeling of a mid-successional forest in northern Michigan

    Bryan, Alexander M.; Cheng, Susan J.; Ashworth, Kirsti; Guenther, Alex B.; Hardiman, Brady S.; Bohrer, Gil; Steiner, Allison L.

    2015-11-01

    Foliar emissions of biogenic volatile organic compounds (BVOC)-important precursors of tropospheric ozone and secondary organic aerosols-vary widely by vegetation type. Modeling studies to date typically represent the canopy as a single dominant tree type or a blend of tree types, yet many forests are diverse with trees of varying height. To assess the sensitivity of biogenic emissions to tree height variation, we compare two 1-D canopy model simulations in which BVOC emission potentials are homogeneous or heterogeneous with canopy depth. The heterogeneous canopy emulates the mid-successional forest at the University of Michigan Biological Station (UMBS). In this case, high-isoprene-emitting foliage (e.g., aspen and oak) is constrained to the upper canopy, where higher sunlight availability increases the light-dependent isoprene emission, leading to 34% more isoprene and its oxidation products as compared to the homogeneous simulation. Isoprene declines from aspen mortality are 10% larger when heterogeneity is considered. Overall, our results highlight the importance of adequately representing complexities of forest canopy structure when simulating light-dependent BVOC emissions and chemistry.

  16. GIS-BASED 1-D DIFFUSIVE WAVE OVERLAND FLOW MODEL

    KALYANAPU, ALFRED [Los Alamos National Laboratory; MCPHERSON, TIMOTHY N. [Los Alamos National Laboratory; BURIAN, STEVEN J. [NON LANL

    2007-01-17

    This paper presents a GIS-based 1-d distributed overland flow model and summarizes an application to simulate a flood event. The model estimates infiltration using the Green-Ampt approach and routes excess rainfall using the 1-d diffusive wave approximation. The model was designed to use readily available topographic, soils, and land use/land cover data and rainfall predictions from a meteorological model. An assessment of model performance was performed for a small catchment and a large watershed, both in urban environments. Simulated runoff hydrographs were compared to observations for a selected set of validation events. Results confirmed the model provides reasonable predictions in a short period of time.

  17. Slug modeling with 1D two-fluid model

    Simulations of condensation-induced water hammer with one-dimensional two-fluid model requires explicit modeling of slug formation, slug propagation, and in some cases slug decay. Stratified flow correlations that are more or less well known in 1D two-fluid models, are crucial for accurate description of the initial phase of the slug formation and slug propagation. Slug formation means transition to other flow regime that requires different set of correlations. To use such two-fluid model for condensation induced water hammer simulations, a single slug must be explicitly recognized and captured. In the present work two cases of condensation-induced water hammer simulations performed with WAHA code, are described and discussed: injection of cold liquid into horizontal pipe filled with steam and injection of hot steam into horizontal pipe partially filled with cold liquid. (author)

  18. Development of a new 1D urban canopy model: coherences between surface parameterizations

    BLOND, Nadège; Mauree, Dasaraden; Kohler, Manon; Clappier, Alain

    2015-01-01

    A 1-D Canopy Interface Model (CIM) was developed in order to better simulate the effect of urban obstacles on the atmosphere in the boundary layer. The model solves the Navier-Stokes equations on a high-resolved gridded vertical column. The effect of the surface is simulated testing a set of theories and urban parameterizations. The final proposition guarantees its coherence with past theories in any atmospheric stability and terrain configuration. Obstacle characteristics are computed using...

  19. A 1-D morphodynamic model of postglacial valley incision

    Tunnicliffe, Jon F.; Church, Michael

    2015-11-01

    Chilliwack River is typical of many Cordilleran valley river systems that have undergone dramatic Holocene degradation of valley fills that built up over the course of Pleistocene glaciation. Downstream controls on base level, mainly blockage of valleys by glaciers, led to aggradation of significant glaciofluvial and glaciolacustrine valley fills and fan deposits, subsequently incised by fluvial action. Models of such large-scale, long-term degradation present a number of important challenges since the evolution of model parameters, such as the rate of bedload transport and grain size characteristics, are governed by the nature of the deposit. Sediment sampling in the Chilliwack Valley reveals a complex sequence of very coarse to fine textural modes. We present a 1-D numerical morphodynamic model for the river-floodplain system tailored to conditions in the valley. The model is adapted to dynamically adjust channel width to optimize sediment transporting capacity and to integrate relict valley fill material as the channel incises through valley deposits. Sensitivity to model parameters is studied using four principal criteria: profile concavity, rate of downstream grain size fining, bed surface sand content, and the timescale to equilibrium. Model results indicate that rates of abrasion and coarsening of the grain size distributions exert the strongest controls on all of the interrelated model performance criteria. While there are a number of difficulties in satisfying all model criteria simultaneously, results indicate that 1-D models of valley bottom sedimentary systems can provide a suitable framework for integrating results from sediment budget studies and chronologies of sediment evacuation established from dating.

  20. On the parameterization of 1D vertical mixing in planetary atmospheres: insights from 2D and 3D simulations

    Zhang, Xi; Showman, Adam P.

    2015-11-01

    Most of the current atmospheric chemistry models for planets (e.g., Krasnopolsky & Parshev 1981; Yung & Demore 1982; Yung, Allen & Pinto 1984; Lavvas et al. 2008; Zhang et al. 2012) and exoplanets (e.g., Line, Liang & Yung 2010; Moses et al. 2011; Hu & Seager 2014) adopt a one-dimensional (1D) chemical-diffusion approach in the vertical coordinate. Although only a crude approximation, these 1D models have succeeded in explaining the global-averaged vertical profiles of many chemical species in observations. One of the important assumptions of these models is that all chemical species are transported via the same eddy diffusion profile--that is, the assumption is made that the eddy diffusivity is a fundamental property of the dynamics alone, and does not depend on the chemistry. Here we show that, as also noticed in the Earth community (e.g., Holton 1986), this “homogenous eddy diffusion” assumption generally breaks down. We first show analytically why the 1D eddy diffusivity must generally depend both on the horizontal eddy mixing and the chemical lifetime of the species. This implies that the long-lived species and short-lived chemical species will generally exhibit different eddy diffusion profiles, even in a given atmosphere with identical dynamics. Next, we present tracer-transport simulations in a 2D chemical-diffusion-advection model (Shia et al. 1989; Zhang, Shia & Yung 2013) and a 3D general circulation model (MITgcm, e.g., Liu & Showman 2013), for both rapid-rotating planets and tidally-locked exoplanets, to further explore the effect of chemical timescales on the eddy diffusivity. From the 2D and 3D simulation outputs, we derive effective 1D eddy diffusivity profiles for chemical tracers exhibiting a range of chemical timescales. We show that the derived eddy diffusivity can depend strongly on the horizontal eddy mixing and chemistry, although the dependences are more complex than the analytic model predicts. Overall, these results suggest that

  1. Modeling atrazine transport in soil columns with HYDRUS-1D

    John Leju CELESTINO LADU

    2011-09-01

    Full Text Available Both physical and chemical processes affect the fate and transport of herbicides. It is useful to simulate these processes with computer programs to predict solute movement. Simulations were run with HYDRUS-1D to identify the sorption and degradation parameters of atrazine through calibration from the breakthrough curves (BTCs. Data from undisturbed and disturbed soil column experiments were compared and analyzed using the dual-porosity model. The study results show that the values of dispersivity are slightly lower in disturbed columns, suggesting that the more heterogeneous the structure is, the higher the dispersivity. Sorption parameters also show slight variability, which is attributed to the differences in soil properties, experimental conditions and methods, or other ecological factors. For both of the columns, the degradation rates were similar. Potassium bromide was used as a conservative non-reactive tracer to characterize the water movement in columns. Atrazine BTCs exhibited significant tailing and asymmetry, indicating non-equilibrium sorption during solute transport. The dual-porosity model was verified to best fit the BTCs of the column experiments. Greater or lesser concentration of atrazine spreading to the bottom of the columns indicated risk of groundwater contamination. Overall, HYDRUS-1D successfully simulated the atrazine transport in soil columns.

  2. Computation of Solar Radiative Fluxes by 1D and 3D Methods Using Cloudy Atmospheres Inferred from A-train Satellite Data

    Barker, Howard W.; Kato, Serji; Wehr, T.

    2012-01-01

    The main point of this study was to use realistic representations of cloudy atmospheres to assess errors in solar flux estimates associated with 1D radiative transfer models. A scene construction algorithm, developed for the EarthCARE satellite mission, was applied to CloudSat, CALIPSO, and MODIS satellite data thus producing 3D cloudy atmospheres measuring 60 km wide by 13,000 km long at 1 km grid-spacing. Broadband solar fluxes and radiances for each (1 km)2 column where then produced by a Monte Carlo photon transfer model run in both full 3D and independent column approximation mode (i.e., a 1D model).

  3. Simplified 1D modelling of the HGA test

    Document available in extended abstract form only. The HGA test is located in the Mont Terri Rock Laboratory (Switzerland). It consists of a horizontal borehole of 1.00 m of diameter and 13.00 m of length excavated in the ultra-low permeable Opalinus clay. During the tunnel drilling, the Opalinus clay near the tunnel wall was damaged, giving rise to an EDZ (Excavation Damaged Zone) around the tunnel. A steel liner was placed along the 6.00 m close to the tunnel mouth in order to guarantee the stability. The last 4.00 m at the tunnel end were backfilled with gravel. Along the remaining 3.00 m, an inflatable rubber packer of 1.00 m in diameter, was installed and inflated, thereby compressing the EDZ that was created during the tunnel excavation. The test section was filled with de-aired water and care was taken in order to eliminate the air from this tunnel section. Subsequently, a series of water and gas injection tests were carried out with varying mega-packer pressure, whereby water or gas was injected into the test section and, due to the very low permeability of the intact Opalinus clay, forced to flow back along the EDZ. In order to model the water and gas flow through the EDZ, we have followed a two-track approach. On the one hand, a 2D axisymmetric numerical model using code-bright has been made. On the other hand, a 1D analytical-numerical model has been developed and implemented in an Excel spreadsheet, whereby the field equations defined on a 1D geometrical domain are numerically solved using the finite element method. The 1D model has been used in order to calibrate the 2D axisymmetric model. Both the Opalinus clay and the EDZ will be considered to be porous media, with an incompressible solid phase (clay), an incompressible liquid phase (water and air) and a gas phase (water and air). The properties of the liquid phase will be assumed to be independent of the concentration of dissolved air and the gas phase will be assumed to be a mixture of dry air and

  4. Modelling turbulent vertical mixing sensitivity using a 1-D version of NEMO

    G. Reffray

    2014-08-01

    Full Text Available Through two numerical experiments, a 1-D vertical model called NEMO1D was used to investigate physical and numerical turbulent-mixing behaviour. The results show that all the turbulent closures tested (k + l from Blanke and Delecluse, 1993 and two equation models: Generic Lengh Scale closures from Umlauf and Burchard, 2003 are able to correctly reproduce the classical test of Kato and Phillips (1969 under favourable numerical conditions while some solutions may diverge depending on the degradation of the spatial and time discretization. The performances of turbulence models were then compared with data measured over a one-year period (mid-2010 to mid-2011 at the PAPA station, located in the North Pacific Ocean. The modelled temperature and salinity were in good agreement with the observations, with a maximum temperature error between −2 and 2 °C during the stratified period (June to October. However the results also depend on the numerical conditions. The vertical RMSE varied, for different turbulent closures, from 0.1 to 0.3 °C during the stratified period and from 0.03 to 0.15 °C during the homogeneous period. This 1-D configuration at the PAPA station (called PAPA1D is now available in NEMO as a reference configuration including the input files and atmospheric forcing set described in this paper. Thus, all the results described can be recovered by downloading and launching PAPA1D. The configuration is described on the NEMO site (http://www.nemo-ocean.eu/Using-NEMO/Configurations/C1D_PAPA. This package is a good starting point for further investigation of vertical processes.

  5. Benchmarking of a 1D scrape-off layer code SOLF1D with SOLPS and its use in modelling long-legged divertors

    Havlickova, E.; Fundamenski, W.; Subba, F.; Coster, D; Wischmeier, M; Fishpool, G.

    2013-01-01

    A 1D code modelling SOL transport parallel to the magnetic field (SOLF1D) is benchmarked with 2D simulations of MAST-U SOL performed via the SOLPS code for two different collisionalities. Based on this comparison, SOLF1D is then used to model the effects of divertor leg stretching in 1D, in support of the planned Super-X divertor on MAST. The aim is to separate magnetic flux expansion from volumetric power losses due to recycling neutrals by stretching the divertor leg either vertically or ra...

  6. Benchmarking of a 1D Scrape-off layer code SOLF1D with SOLPS and its use in modelling long-legged divertors

    Havlickova, E; Subba, F; Coster, D; Wischmeier, M; Fishpool, G

    2013-01-01

    A 1D code modelling SOL transport parallel to the magnetic field (SOLF1D) is benchmarked with 2D simulations of MAST-U SOL performed via the SOLPS code for two different collisionalities. Based on this comparison, SOLF1D is then used to model the effects of divertor leg stretching in 1D, in support of the planned Super-X divertor on MAST. The aim is to separate magnetic flux expansion from volumetric power losses due to recycling neutrals by stretching the divertor leg either vertically or radially.

  7. Testing the accuracy of a 1-D volcanic plume model in estimating mass eruption rate

    Mastin, Larry G.

    2014-01-01

    During volcanic eruptions, empirical relationships are used to estimate mass eruption rate from plume height. Although simple, such relationships can be inaccurate and can underestimate rates in windy conditions. One-dimensional plume models can incorporate atmospheric conditions and give potentially more accurate estimates. Here I present a 1-D model for plumes in crosswind and simulate 25 historical eruptions where plume height Hobs was well observed and mass eruption rate Mobs could be calculated from mapped deposit mass and observed duration. The simulations considered wind, temperature, and phase changes of water. Atmospheric conditions were obtained from the National Center for Atmospheric Research Reanalysis 2.5° model. Simulations calculate the minimum, maximum, and average values (Mmin, Mmax, and Mavg) that fit the plume height. Eruption rates were also estimated from the empirical formula Mempir = 140Hobs4.14 (Mempir is in kilogram per second, Hobs is in kilometer). For these eruptions, the standard error of the residual in log space is about 0.53 for Mavg and 0.50 for Mempir. Thus, for this data set, the model is slightly less accurate at predicting Mobs than the empirical curve. The inability of this model to improve eruption rate estimates may lie in the limited accuracy of even well-observed plume heights, inaccurate model formulation, or the fact that most eruptions examined were not highly influenced by wind. For the low, wind-blown plume of 14–18 April 2010 at Eyjafjallajökull, where an accurate plume height time series is available, modeled rates do agree better with Mobs than Mempir.

  8. Validation of 1-D transport and sawtooth models for ITER

    In this paper the authors describe progress on validating a number of local transport models by comparing their predictions with relevant experimental data from a range of tokamaks in the ITER profile database. This database, the testing procedure and results are discussed. In addition a model for sawtooth oscillations is used to investigate their effect in an ITER plasma with alpha-particles

  9. Analyzing and modeling /1+1d markets

    Challet, Damien; Stinchcombe, Robin

    2001-11-01

    We report on a statistical analysis of the Island ECN (NASDAQ) order book. We determine the static and dynamic properties of this system, and then analyze them from a physicist's viewpoint using an equivalent particle system obtained by treating orders as massive particles and price as position. We identify the fundamental dynamical processes, test existing particles models of such markets against our findings, and introduce a new model of limit order markets.

  10. Control oriented 1D electrochemical model of lithium ion battery

    Lithium ion (Li-ion) batteries provide high energy and power density energy storage for diverse applications ranging from cell phones to hybrid electric vehicles (HEVs). For efficient and reliable systems integration, low order dynamic battery models are needed. This paper introduces a general method to generate numerically a fully observable/controllable state variable model from electrochemical kinetic, species and charge partial differential equations that govern the discharge/charge behavior of a Li-ion battery. Validated against a 313th order nonlinear CFD model of a 6 Ah HEV cell, a 12th order state variable model predicts terminal voltage to within 1% for pulse and constant current profiles at rates up to 50 C. The state equation is constructed in modal form with constant negative real eigenvalues distributed in frequency space from 0 to 10 Hz. Open circuit potential, electrode surface concentration/reaction distribution coupling and electrolyte concentration/ionic conductivity nonlinearities are explicitly approximated in the model output equation on a local, electrode-averaged and distributed basis, respectively. The balanced realization controllability/observability gramian indicates that the fast electrode surface concentration dynamics are more observable/controllable than the electrode bulk concentration dynamics (i.e. state of charge)

  11. Kinetic and Stochastic Models of 1D yeast ``prions"

    Kunes, Kay

    2005-03-01

    Mammalian prion proteins (PrP) are of public health interest because of mad cow and chronic wasting diseases. Yeasts have proteins, which can undergo similar reconformation and aggregation processes to PrP; yeast ``prions" are simpler to experimentally study and model. Recent in vitro studies of the SUP35 protein (1), showed long aggregates and pure exponential growth of the misfolded form. To explain this data, we have extended a previous model of aggregation kinetics along with our own stochastic approach (2). Both models assume reconformation only upon aggregation, and include aggregate fissioning and an initial nucleation barrier. We find for sufficiently small nucleation rates or seeding by small dimer concentrations that we can achieve the requisite exponential growth and long aggregates.

  12. Application of particle trajectory model in 1D planar ejection

    刘坤; 柏劲松; 李平

    2008-01-01

    A simple one-dimensional planar model for ejection was set up based on experiments.And numerical simulation was performed on this model with particle trajectory model method.An Eulerian finite volume method was conducted to resolve gas field.And Lagrangian method was imposed to track each particle.The interaction between gas and particles was responded as source terms in governing equations which were induced by forces.The effects of total spraying mass,particle size and other factors on the mixture of particles and gas were investigated.The spatial distributions of particle mass and velocity at different time were presented.The result shows that the numerical results are qualitatively consistent to those of experiments.

  13. Using a 1-D model to reproduce diurnal SST signals

    Karagali, Ioanna; Høyer, Jacob L.

    2014-01-01

    profiles, along with the selection of the coefficients for the 2-band parametrisation of light’s penetration in the water column, hold a key role in the agreement of the modelled output with observations. To improve the surface heat budget and the distribution of heat, the code was modified to include an...

  14. Kinetic Model for 1D aggregation of yeast ``prions''

    Kunes, Kay; Cox, Daniel; Singh, Rajiv

    2004-03-01

    Mammalian prion proteins (PrP) are of public health interest because of mad cow and chronic wasting diseases. Yeast have proteins which can undergo similar reconformation and aggregation processes to PrP; yeast forms are simpler to experimentally study and model. Recent in vitro studies of the SUP35 protein(1), showed long aggregates and pure exponential growth of the misfolded form. To explain this data, we have extended a previous model of aggregation kinetics(2). The model assumes reconformation only upon aggregation, and includes aggregate fissioning and an initial nucleation barrier. We find for sufficiently small nucleation rates or seeding by small dimer concentrations that we can achieve the requisite exponential growth and long aggregates. We will compare to a more realistic stochastic kinetics model and present prelimary attempts to describe recent experiments on SUP35 strains. *-Supported by U.S. Army Congressionally Mandated Research Fund. 1) P. Chien and J.S. Weissman, Nature 410, 223 (2001); http://online.kitp.ucsb.edu/online/bionet03/collins/. 2) J. Masel, V.A.> Jansen, M.A. Nowak, Biophys. Chem. 77, 139 (1999).

  15. Nonisothermal Pluto atmosphere models

    The present thermal profile calculation for a Pluto atmosphere model characterized by a high number fraction of CH4 molecules encompasses atmospheric heating by solar UV flux absorption and conductive transport cooling to the surface of Pluto. The stellar occultation curve predicted for an atmosphere of several-microbar surface pressures (which entail the existence of a substantial temperature gradient close to the surface) agrees with observations and implies that the normal and tangential optical depth of the atmosphere is almost negligible. The minimum period for atmospheric methane depletion is calculated to be 30 years. 29 refs

  16. Testing the Early Mars H2-CO2 Greenhouse Hypothesis with a 1-D Photochemical Model

    Batalha, Natasha; Ramirez, Ramses; Kasting, James

    2015-01-01

    A recent study by Ramirez et al. (2014) demonstrated that an atmosphere with 1.3-4 bar of CO2 and H2O, in addition to 5-20% H2, could have raised the mean annual and global surface temperature of early Mars above the freezing point of water. Such warm temperatures appear necessary to generate the rainfall (or snowfall) amounts required to carve the ancient martian valleys. Here, we use our best estimates for early martian outgassing rates, along with a 1-D photochemical model, to assess the conversion efficiency of CO, CH4, and H2S to CO2, SO2, and H2. Our outgassing estimates assume that Mars was actively recycling volatiles between its crust and interior, as Earth does today. H2 production from serpentinization and deposition of banded iron-formations is also considered. Under these assumptions, maintaining an H2 concentration of ~1-2% by volume is achievable, but reaching 5% H2 requires additional H2 sources or a slowing of the hydrogen escape rate below the diffusion limit. If the early martian atmosphere...

  17. Modeling blood flow circulation in intracranial arterial networks: a comparative 3D/1D simulation study.

    Grinberg, L; Cheever, E; Anor, T; Madsen, J R; Karniadakis, G E

    2011-01-01

    We compare results from numerical simulations of pulsatile blood flow in two patient-specific intracranial arterial networks using one-dimensional (1D) and three-dimensional (3D) models. Specifically, we focus on the pressure and flowrate distribution at different segments of the network computed by the two models. Results obtained with 1D and 3D models with rigid walls show good agreement in massflow distribution at tens of arterial junctions and also in pressure drop along the arteries. The 3D simulations with the rigid walls predict higher amplitude of the flowrate and pressure temporal oscillations than the 1D simulations with compliant walls at various segments even for small time-variations in the arterial cross-sectional areas. Sensitivity of the flow and pressure with respect to variation in the elasticity parameters is investigated with the 1D model. PMID:20661645

  18. 1D model for the dynamics and expansion of elongated Bose-Einstein condensates

    Massignan, Pietro; Modugno, Michele

    2002-01-01

    We present a 1D effective model for the evolution of a cigar-shaped Bose-Einstein condensate in time dependent potentials whose radial component is harmonic. We apply this model to investigate the dynamics and expansion of condensates in 1D optical lattices, by comparing our predictions with recent experimental data and theoretical results. We also discuss negative-mass effects which could be probed during the expansion of a condensate moving in an optical lattice.

  19. Simple model of the density of states in 1D photonic crystal

    Rudziński, Adam; Tyszka-Zawadzka, Anna; Szczepański, Paweł

    2010-01-01

    In this paper, we present a simple, yet versatile, analytical model of one-dimensional photonic crystal (1D PC). In our theoretical model, we take into account direction of propagation and therefore do not neglect anisotropic nature of photonic crystals. We derive analytical expressions for mode spectrum and density of states in 1D photonic crystal. With those formulas, we obtain mode spectrum characteristics, which depict formation of photonic band gap and reveal properties of photonic cryst...

  20. How realistic are solar model atmospheres?

    Pereira, Tiago M D; Collet, Remo; Thaler, Irina; Trampedach, Regner; Leenaarts, Jorrit

    2013-01-01

    Recently, new solar model atmospheres have been developed to replace classical 1D LTE hydrostatic models and used to for example derive the solar chemical composition. We aim to test various models against key observational constraints. In particular, a 3D model used to derive the solar abundances, a 3D MHD model (with an imposed 10 mT vertical magnetic field), 1D models from the PHOENIX project, the 1D MARCS model, and the 1D semi-empirical model of Holweger & M\\"uller. We confront the models with observational diagnostics of the temperature profile: continuum centre-to-limb variations (CLV), absolute continuum fluxes, and the wings of hydrogen lines. We also test the 3D models for the intensity distribution of the granulation and spectral line shapes. The predictions from the 3D model are in excellent agreement with the continuum CLV observations, performing even better than the Holweger & M\\"uller model (constructed largely to fulfil such observations). The predictions of the 1D theoretical models ...

  1. REAL-TIME FLOOD FORECASTING METHOD WITH 1-D UNSTEADY FLOW MODEL

    MU Jin-bin; ZHANG Xiao-feng

    2007-01-01

    A real-time forecasting method coupled with the 1-D unsteady flow model with the recursive least-square method was developed. The 1-D unsteady flow model was modified by using the time-variant parameter and revising it dynamically through introducing a variable weighted forgetting factor, such that the output of the model could be adjusted for the real time forecasting of floods. The application of the new real time forecasting model in the reach from Yichang to Luoshan of the Yangtze River was demonstrated. Computational result shows that the forecasting accuracy of the new model is much higher than that of the original 1-D unsteady flow model. The method developed is effective for flood forecasting, and can be used for practical operation in the flood forecasting.

  2. DEVELOPMENT OF COUPLED 1D-2D MATHEMATICAL MODELS FOR TIDAL RIVERS

    XU Zu-xin; YIN Hai-long

    2004-01-01

    Some coupled 1D-2D hydrodynamic and water quality models depicting tidal water bodies with complex topography were presented. For the coupled models, finite element method was used to solve the governing equations so as to study tidal rivers with complex topography. Since the 1D and 2D models were coupled, the principle of model coupling was proposed to account appropriately for the factors of water level, flow and pollutant flux and the related dynamical behavior was simulated. Specifically the models were used to probe quantitative pollution contribution of receiving water from neighboring Jiangsu and Zhejiang Provinces to the pollution in the Huangpu River passing through Shanghai City. Numerical examples indicated that the developed coupled 1D-2D models are applicable in tidal river network region of Shanghai.

  3. Benchmarks and models for 1-D radiation transport in stochastic participating media

    Miller, D S

    2000-08-21

    Benchmark calculations for radiation transport coupled to a material temperature equation in a 1-D slab and 1-D spherical geometry binary random media are presented. The mixing statistics are taken to be homogeneous Markov statistics in the 1-D slab but only approximately Markov statistics in the 1-D sphere. The material chunk sizes are described by Poisson distribution functions. The material opacities are first taken to be constant and then allowed to vary as a strong function of material temperature. Benchmark values and variances for time evolution of the ensemble average of material temperature energy density and radiation transmission are computed via a Monte Carlo type method. These benchmarks are used as a basis for comparison with three other approximate methods of solution. One of these approximate methods is simple atomic mix. The second approximate model is an adaptation of what is commonly called the Levermore-Pomraning model and which is referred to here as the standard model. It is shown that recasting the temperature coupling as a type of effective scattering can be useful in formulating the third approximate model, an adaptation of a model due to Su and Pomraning which attempts to account for the effects of scattering in a stochastic context. This last adaptation shows consistent improvement over both the atomic mix and standard models when used in the 1-D slab geometry but shows limited improvement in the 1-D spherical geometry. Benchmark values are also computed for radiation transmission from the 1-D sphere without material heating present. This is to evaluate the performance of the standard model on this geometry--something which has never been done before. All of the various tests demonstrate the importance of stochastic structure on the solution. Also demonstrated are the range of usefulness and limitations of a simple atomic mix formulation.

  4. Deconvolution of Complex 1D NMR Spectra Using Objective Model Selection.

    Travis S Hughes

    Full Text Available Fluorine (19F NMR has emerged as a useful tool for characterization of slow dynamics in 19F-labeled proteins. One-dimensional (1D 19F NMR spectra of proteins can be broad, irregular and complex, due to exchange of probe nuclei between distinct electrostatic environments; and therefore cannot be deconvoluted and analyzed in an objective way using currently available software. We have developed a Python-based deconvolution program, decon1d, which uses Bayesian information criteria (BIC to objectively determine which model (number of peaks would most likely produce the experimentally obtained data. The method also allows for fitting of intermediate exchange spectra, which is not supported by current software in the absence of a specific kinetic model. In current methods, determination of the deconvolution model best supported by the data is done manually through comparison of residual error values, which can be time consuming and requires model selection by the user. In contrast, the BIC method used by decond1d provides a quantitative method for model comparison that penalizes for model complexity helping to prevent over-fitting of the data and allows identification of the most parsimonious model. The decon1d program is freely available as a downloadable Python script at the project website (https://github.com/hughests/decon1d/.

  5. Solar Atmosphere Models

    Rutten, R. J.

    2002-12-01

    This contribution honoring Kees de Jager's 80th birthday is a review of "one-dimensional" solar atmosphere modeling that followed on the initial "Utrecht Reference Photosphere" of Heintze, Hubenet & de Jager (1964). My starting point is the Bilderberg conference, convened by de Jager in 1967 at the time when NLTE radiative transfer theory became mature. The resulting Bilderberg model was quickly superseded by the HSRA and later by the VAL-FAL sequence of increasingly sophisticated NLTE continuum-fitting models from Harvard. They became the "standard models" of solar atmosphere physics, but Holweger's relatively simple LTE line-fitting model still persists as a favorite of solar abundance determiners. After a brief model inventory I discuss subsequent work on the major modeling issues (coherency, NLTE, dynamics) listed as to-do items by de Jager in 1968. The present conclusion is that one-dimensional modeling recovers Schwarzschild's (1906) finding that the lower solar atmosphere is grosso modo in radiative equilibrium. This is a boon for applications regarding the solar atmosphere as one-dimensional stellar example - but the real sun, including all the intricate phenomena that now constitute the mainstay of solar physics, is vastly more interesting.

  6. Technical Note: Sensitivity of 1-D smoke plume rise models to the inclusion of environmental wind drag

    S. R. Freitas

    2010-01-01

    Full Text Available Vegetation fires emit hot gases and particles which are rapidly transported upward by the positive buoyancy generated by the combustion process. In general, the final vertical height that the smoke plumes reach is controlled by the thermodynamic stability of the atmospheric environment and the surface heat flux released by the fire. However, the presence of a strong horizontal wind can enhance the lateral entrainment and induce additional drag, particularly for small fires, impacting the smoke injection height. In this paper, we revisit the parameterization of the vertical transport of hot gases and particles emitted from vegetation fires, described in Freitas et al. (2007, to include the effects of environmental wind on transport and dilution of the smoke plume at its scale. This process is quantitatively represented by introducing an additional entrainment term to account for organized inflow of a mass of cooler and drier ambient air into the plume and its drag by momentum transfer. An extended set of equations including the horizontal motion of the plume and the additional increase of the plume radius is solved to simulate the time evolution of the plume rise and the smoke injection height. One-dimensional (1-D model results are presented for two deforestation fires in the Amazon basin with sizes of 10 and 50 ha under calm and windy atmospheric environments. The results are compared to corresponding simulations generated by the complex non-hydrostatic three-dimensional (3-D Active Tracer High resolution Atmospheric Model (ATHAM. We show that the 1-D model results compare well with the full 3-D simulations. The 1-D model may thus be used in field situations where extensive computing facilities are not available, especially under conditions for which several optional cases must be studied.

  7. Technical Note: Sensitivity of 1-D smoke plume rise models to the inclusion of environmental wind drag

    S. R. Freitas

    2009-07-01

    Full Text Available We revisit the parameterization of the vertical transport of hot gases and particles emitted from biomass burning, described in Freitas et al. (2007, to include the effects of environmental wind on transport and dilution of the smoke plume at the cloud scale. Typically, the final vertical height that the smoke plumes reach is controlled by the thermodynamic stability of the atmospheric environment and the surface heat flux released by the fire. However, the presence of a strong horizontal wind can enhance the lateral entrainment and induce additional drag, particularly for small fires, impacting the smoke injection height. This process is quantitatively represented by introducing an additional entrainment term to account for organized inflow of a mass of cooler and drier ambient air into the plume and its drag by momentum transfer. An extended set of equations including the horizontal motion of the plume and the additional increase of the plume radius is solved to explicitly simulate the time evolution of the plume rise with the additional mass and momentum. One-dimensional (1-D model results are presented for two deforestation fires in the Amazon basin with sizes of 10 and 50 ha under calm and windy atmospheric environments. The results are compared to corresponding simulations generated by the complex non-hydrostatic three dimensional (3-D Active Tracer High resolution Atmospheric Model (ATHAM. We show that the 1-D model results compare well with the full 3-D simulations. The 1-D model may thus be used in field situations where extensive computing facilities are not available, especially under conditions for which several optional cases must be studied.

  8. Partitioning of evaporation into transpiration, soil evaporation and interception: a comparison between isotope measurements and a HYDRUS-1D model

    Sutanto, S. J.; Wenninger, J.; Coenders-Gerrits, A. M. J.; Uhlenbrook, S.

    2012-08-01

    Knowledge of the water fluxes within the soil-vegetation-atmosphere system is crucial to improve water use efficiency in irrigated land. Many studies have tried to quantify these fluxes, but they encountered difficulties in quantifying the relative contribution of evaporation and transpiration. In this study, we compared three different methods to estimate evaporation fluxes during simulated summer conditions in a grass-covered lysimeter in the laboratory. Only two of these methods can be used to partition total evaporation into transpiration, soil evaporation and interception. A water balance calculation (whereby rainfall, soil moisture and percolation were measured) was used for comparison as a benchmark. A HYDRUS-1D model and isotope measurements were used for the partitioning of total evaporation. The isotope mass balance method partitions total evaporation of 3.4 mm d-1 into 0.4 mm d-1 for soil evaporation, 0.3 mm d-1 for interception and 2.6 mm d-1 for transpiration, while the HYDRUS-1D partitions total evaporation of 3.7 mm d-1 into 1 mm d-1 for soil evaporation, 0.3 mm d-1 for interception and 2.3 mm d-1 for transpiration. From the comparison, we concluded that the isotope mass balance is better for low temporal resolution analysis than the HYDRUS-1D. On the other hand, HYDRUS-1D is better for high temporal resolution analysis than the isotope mass balance.

  9. Periodic Properties of 1D FE Discrete Models in High Frequency Dynamics

    A. Żak

    2016-01-01

    Full Text Available Finite element discrete models of various engineering 1D structures may be considered as structures of certain periodic characteristics. The source of this periodicity comes from the discontinuity of stress/strain field between the elements. This behaviour remains unnoticeable, when low frequency dynamics of these structures is investigated. At high frequency regimes, however, its influence may be strong enough to dominate calculated structural responses distorting or even falsifying them completely. In this paper, certain computational aspects of structural periodicity of 1D FE discrete models are discussed by the authors. In this discussion, the authors focus their attention on an exemplary problem of 1D rod modelled according to the elementary theory.

  10. A Mathematical Model of T1D Acceleration and Delay by Viral Infection.

    Moore, James R; Adler, Fred

    2016-03-01

    Type 1 diabetes (T1D) is often triggered by a viral infection, but the T1D prevalence is rising among populations that have a lower exposure to viral infection. In an animal model of T1D, the NOD mouse, viral infection at different ages may either accelerate or delay disease depending on the age of infection and the type of virus. Viral infection may affect the progression of T1D via multiple mechanisms: triggering inflammation, bystander activation of self-reactive T-cells, inducing a competitive immune response, or inducing a regulatory immune response. In this paper, we create mathematical models of the interaction of viral infection with T1D progression, incorporating each of these four mechanisms. Our goal is to understand how each viral mechanism interacts with the age of infection. The model predicts that each viral mechanism has a unique pattern of interaction with disease progression. Viral inflammation always accelerates disease, but the effect decreases with age of infection. Bystander activation has little effect at younger ages and actually decreases incidence at later ages while accelerating disease in mice that do get the disease. A competitive immune response to infection can decrease incidence at young ages and increase it at older ages, with the effect decreasing over time. Finally, an induced Treg response decreases incidence at any age of infection, but the effect decreases with age. Some of these patterns resemble those seen experimentally. PMID:27030351

  11. Column Testing and 1D Reactive Transport Modeling to Evaluate Uranium Plume Persistence Processes

    Johnson, Raymond H. [Navarro Research and Engineering, Inc.; Morrison, Stan [Navarro Research and Engineering, Inc.; Morris, Sarah [Navarro Research and Engineering, Inc.; Tigar, Aaron [Navarro Research and Engineering, Inc.; Dam, William [U.S. Department of Energy, Office of Legacy Management; Dayvault, Jalena [U.S. Department of Energy, Office of Legacy Management

    2016-04-26

    Motivation for Study: Natural flushing of contaminants at various U.S. Department of Energy Office of Legacy Management sites is not proceeding as quickly as predicted (plume persistence) Objectives: Help determine natural flushing rates using column tests. Use 1D reactive transport modeling to better understand the major processes that are creating plume persistence Approach: Core samples from under a former mill tailings area Tailings have been removed. Column leaching using lab-prepared water similar to nearby Gunnison River water. 1D reactive transport modeling to evaluate processes

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

    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

  13. HYDRUS-1D Modeling of an Irrigated Agricultural Plot with Application to Aquifer Recharge Estimation

    A variety of methods are available for estimating aquifer recharge in semi-arid regions, each with advantages and disadvantages. We are investigating a procedure for estimating recharge in an irrigated basin. The method involves computing irrigation return flows based on HYDRUS-1D modeling of root z...

  14. Minimal representations of supersymmetry and 1D N-extended σ-models

    We discuss the minimal representations of the 1D N-Extended Supersymmetry algebra (the Z2-graded symmetry algebra of the Supersymmetric Quantum Mechanics) linearly realized on a finite number of fields depending on a real parameter t, the time. Their knowledge allows to construct one-dimensional sigma-models with extended off-shell supersymmetries without using superfields (author)

  15. Thermodynamics of 1D N-Component Bariev Model Under Open Boundary Conditions

    WANG Chun; KE San-Min; YUE Rui-Hong

    2006-01-01

    The thermodynamic Bethe ansatz equations and free energy for 1D N-component Bariev model under open boundary conditions are derived based on the string hypothesis for both, a repulsive and an attractive interaction.These equations are discussed in some limiting cases, such as the ground state, weak and strong couplings.

  16. A 1D microphysical cloud model for Earth, and Earth-like exoplanets. Liquid water and water ice clouds in the convective troposphere

    Zsom, A; Goldblatt, C

    2012-01-01

    One significant difference between the atmospheres of stars and exoplanets is the presence of condensed particles (clouds or hazes) in the atmosphere of the latter. The main goal of this paper is to develop a self-consistent microphysical cloud model for 1D atmospheric codes, which can reproduce some observed properties of Earth, such as the average albedo, surface temperature, and global energy budget. The cloud model is designed to be computationally efficient, simple to implement, and applicable for a wide range of atmospheric parameters for planets in the habitable zone. We use a 1D, cloud-free, radiative-convective, and photochemical equilibrium code originally developed by Kasting, Pavlov, Segura, and collaborators as basis for our cloudy atmosphere model. The cloud model is based on models used by the meteorology community for Earth's clouds. The free parameters of the model are the relative humidity and number density of condensation nuclei, and the precipitation efficiency. In a 1D model, the cloud c...

  17. Chemical uncertainties in modeling hot Jupiters atmospheres

    Hebrard, Eric; Domagal-Goldman, Shawn

    2015-11-01

    Most predictions and interpretations of observations in beyond our Solar System have occurred through the use of 1D photo-thermo-chemical models. Their predicted atmospheric compositions are highly dependent on model parameters. Chemical reactions are based on empirical parameters that must be known at temperatures ranging from 100 K to above 2500 K and at pressures from millibars to hundreds of bars. Obtained from experiments, calculations and educated-guessed estimations, these parameters are always evaluated with substantial uncertainties. However, although of practical use, few models of exoplanetary atmospheres have considered these underlying chemical uncertainties and their consequences. Recent progress has been made recently that allow us to (1) evaluate the accuracy and precision of 1D models of planetary atmospheres, with quantifiable uncertainties on their predictions for the atmospheric composition and associated spectral features, (2) identify the ‘key parameters’ that contribute the most to the models predictivity and should therefore require further experimental or theoretical analysis, (3) reduce and optimize complex chemical networks for their inclusion in multidimensional atmospheric models.First, a global sampling approach based on low discrepancy sequences has been applied in order to propose error bars on simulations of the atmospheres HD 209458b and HD 189733b, using a detailed kinetic model derived from applied combustion models that was methodically validated over a range of temperatures and pressures typical for these hot Jupiters. A two-parameters temperature-dependent uncertainty factor has been assigned to each considered rate constant. Second, a global sensitivity approach based on high dimensional model representations (HDMR) has been applied in order to identify those reactions which make the largest contributions to the overall uncertainty of the simulated results. The HDMR analysis has been restricted to the most important

  18. Silva. EDF two-phase 1D annular model of a CFB boiler furnace

    Montat, D.; Fauquet, Ph. [Electricite de France (EDF), 78 - Chatou (France). Researckh and Development Div.; Lafanechere, L.; Bursi, J.M. [Electricite de France (EDF) (France). Construction Div.

    1997-01-01

    SILVA computer code is used for the modelling of the thermal-hydraulics and of the combustion of a coal-fired CFBC solid loop. In a first step, only the furnace is considered. The model is based on a 1D annular two phases description of the hydrodynamics. The model is based on particle mass balances and pressure drop calculations. A basic combustion model is incorporated into this model. The coal combustion is divided in two phases, the combustion of volatile matter and the heterogeneous combustion. The model has been developed within LEGO software and can be included into the global model of the solid loop developed by EDF. (author) 26 refs.

  19. A 1D model for the description of mixing-controlled reacting diesel sprays

    Desantesa, J.M.; Pastor, J.V.; Garcia-Oliver, J.M.; Pastor, J.M. [CMT - Motores Termicos, Universidad Politecnica de Valencia, Camino de Vera s/n, 46022, Valencia (Spain)

    2009-01-15

    The paper reports an investigation on the transient evolution of diesel flames in terms of fuel-air mixing, spray penetration and combustion rate. A one-dimensional (1D) spray model, which was previously validated for inert diesel sprays, is extended to reacting conditions. The main assumptions of the model are the mixing-controlled hypothesis and the validity of self-similarity for conservative properties. Validation is achieved by comparing model predictions with both CFD gas jet simulations and experimental diesel spray measurements. The 1D model provides valuable insight into the evolution of the flow within the spray (momentum and mass fluxes, tip penetration, etc.) when shifting from inert to reacting conditions. Results show that the transient diesel flame evolution is mainly governed by two combustion-induced effects, namely the reduction in local density and the increase in flame radial width. (author)

  20. Prediction of car cabin environment by means of 1D and 3D cabin model

    Jícha M.; Pokorný J.; Fišer J.

    2012-01-01

    Thermal comfort and also reduction of energy requirements of air-conditioning system in vehicle cabins are currently very intensively investigated and up-to-date issues. The article deals with two approaches of modelling of car cabin environment; the first model was created in simulation language Modelica (typical 1D approach without cabin geometry) and the second one was created in specialized software Theseus-FE (3D approach with cabin geometry). Performance and capabilities of this t...

  1. Local source tomography for Vrancea (Romania) region: minimum 1D model

    The Vrancea zone is characterized by a narrow, vertical focal volume of intense and persistent seismicity in the depth range of 60-220 km. The processes associated with the Vrancea earthquakes outline a final stage of continental subduction with plate detachment, still far from being understood. Our main purpose is to perform local earthquake tomography to image crustal structures in Vrancea and adjacent regions as initial reference model in 3 D tomography. The study includes updating of the existing catalogues of seismic parameter data, calculation of a minimum 1D model and correlation with the existing controlled-source data. The model is determined by a trial-and-error process, with internal non-linear (iterative) inversion procedure (program VELEST). The average velocity model parameters are obtained by minimizing the standard deviations of the travel time residuals. The P-wave travel times for 500 well-locatable events are used, merging travel time data recorded by two separate independent seismic networks: the Romanian telemetered network (17 vertical-component S-13 seismometers) and the digital accelerometers network (30 three-component K2 instruments) installed in 1996 within the joint Romanian-German cooperation programme CRC 461. Since the two networks largely overlap geographically, we calculate a common minimum 1D model that guarantees a uniform location quality and uniform phase identification for the two networks. Absolute mislocation errors for the combined networks using the minimum 1D model are calculated using quarry blast data obtained during the VRANCEA99 refraction experiment deployed in 1999 in Romania as a contribution to the CRC 461 programme. The resulted minimum 1D model and station corrections are used to relocate all the considered events and to improve the consistency in the phase identification. The station corrections obtained are correlated with the lateral velocity variations in the surface geology of the region. This paper is a

  2. Numerical Methods and Comparisons for 1D and Quasi 2D Streamer Propagation Models

    Huang, Mengmin; Guan, Huizhe; Zeng, Rong

    2016-01-01

    In this work, we propose four different strategies to simulate the one-dimensional (1D) and quasi two-dimensional (2D) model for streamer propagation. Each strategy involves of one numerical method for solving Poisson's equation and another method for solving continuity equations in the models, and a total variation diminishing three-stage Runge-Kutta method in temporal discretization. The numerical methods for Poisson's equation include finite volume method, discontinuous Galerkin methods, mixed finite element method and least-squared finite element method. The numerical method for continuity equations is chosen from the family of discontinuous Galerkin methods. The accuracy tests and comparisons show that all of these four strategies are suitable and competitive in streamer simulations from the aspects of accuracy and efficiency. By applying any strategy in real simulations, we can study the dynamics of streamer propagations and influences due to the change of parameters in both of 1D and quasi 2D models. T...

  3. Thermodynamic nature of vitrification in a 1D model of a structural glass former

    We propose a new spin-glass model with no positional quenched disorder which is regarded as a coarse-grained model of a structural glass-former. The model is analyzed in the 1D case when the number N of states of a primary cell is large. For N → ∞, the model exhibits a sharp freezing transition of the thermodynamic origin. It is shown both analytically and numerically that the glass transition is accompanied by a significant growth of a static length scale ξ pointing to the structural (equilibrium) nature of dynamical slowdown effects in supercooled liquids

  4. Assessment of core thermo-hydrodynamic models of REFLA-1D with CCTF data

    In order to assess the core thermo-hydrodynamic models of REFLA-1D/MODE3, which is the latest version of REFLA-1D, several calculations of the core thermo-hydrodynamics have been performed for the CCTF Core-I series tests. The measured initial and boundary conditions were used for these calculations. The calculational results showed that the water accumulation model of Case 2 could predict the CCTF results fairly well as it could for the JAERI small scale facility. The calculated results for the base case and the EM tests were in good agreement with the CCTF data. The parameter effects, such as system pressure, initial clad temperature, Acc injection rate, LPCI injection rate and initial down-comer wall temperature, were predicted correctly, except for the high system pressure and the high LPCI injection rate tests. (author)

  5. Partitioning of evaporation into transpiration, soil evaporation and interception: a comparison between isotope measurements and a HYDRUS-1D model

    S. J. Sutanto

    2012-08-01

    Full Text Available Knowledge of the water fluxes within the soil-vegetation-atmosphere system is crucial to improve water use efficiency in irrigated land. Many studies have tried to quantify these fluxes, but they encountered difficulties in quantifying the relative contribution of evaporation and transpiration. In this study, we compared three different methods to estimate evaporation fluxes during simulated summer conditions in a grass-covered lysimeter in the laboratory. Only two of these methods can be used to partition total evaporation into transpiration, soil evaporation and interception. A water balance calculation (whereby rainfall, soil moisture and percolation were measured was used for comparison as a benchmark. A HYDRUS-1D model and isotope measurements were used for the partitioning of total evaporation. The isotope mass balance method partitions total evaporation of 3.4 mm d−1 into 0.4 mm d−1 for soil evaporation, 0.3 mm d−1 for interception and 2.6 mm d−1 for transpiration, while the HYDRUS-1D partitions total evaporation of 3.7 mm d−1 into 1 mm d−1 for soil evaporation, 0.3 mm d−1 for interception and 2.3 mm d−1 for transpiration. From the comparison, we concluded that the isotope mass balance is better for low temporal resolution analysis than the HYDRUS-1D. On the other hand, HYDRUS-1D is better for high temporal resolution analysis than the isotope mass balance.

  6. Box model and 1D longitudinal model of flow and transport in Bosten Lake, China

    Li, Ning; Kinzelbach, Wolfgang; Li, WenPeng; Dong, XinGuang

    2015-05-01

    Bosten Lake in the southeast of Yanqi Catchment, China, supports the downstream agricultural and natural environments. Over the last few decades the intensive agricultural activities in Yanqi Catchment resulted in decreased lake levels and deteriorated lake water quality. A two-box model is constructed to understand the evolution of lake level and salinity between 1958 and 2008. The two-box model of the lake indicates that the evaporation does have the same trend as the observed lake area and the annual average evaporation agrees with the value obtained from the Penman-Monteith approach. To achieve a correct salt balance, the ratio of outflow concentration and average lake concentration has to be around 0.7. This is due to the incomplete mixing of the lake caused by short-circuiting between tributary inflow and the main outflow via the pump stations abstracting water from the lake. This short-circuiting is investigated in more detail by a 1D numerical flow and transport model of the lake calibrated with observations of lake level and lake concentrations. The distributed model reproduces the correct time-varying outflow concentration. It is used for the assessment of two basic management options: increasing river discharge (by water saving irrigation, reduction of phreatic evaporation or reduction of agricultural area) and diverting saline drainage water to the desert. Increasing river discharge to the lake by 20% reduces the east basin salt concentration by 0.55 kg/m3, while capturing all the drainage water and discharging it to depressions instead of the lake reduces the east basin salt concentration by 0.63 kg/m3. A combination of increasing river inflow and decreasing drainage salt flux is sufficient to bring future lake TDS below the required 1 kg/m3, to keep a lake level that sustains the lake ecosystem, and to supply more water for downstream development and ecosystem rehabilitation.

  7. Modelling activities of experimental facilities related to advanced reactors. Considerations on 1D/3D issues

    The state of art of modelling activities related to integral experimental facilities of advanced passive reactors show to date important open items. The main advantage of using 1D plant codes is the capability of simulating the full interaction between components traditionally correctly modelled (condensers, heat exchangers, pipes and vessels) and other components for which codes are not 100% suitable (pools and containments). Polytechnical University of Catalonia (UPC) and Polytechnical University of Valencia (UPV) cooperated with other European research organizations in the 'Technology Enhancement for Passive Safety Systems' (TEPSS) project, within the European Fourth Framework Programme. It was a task of both Universities to supply analytical support of PANDA tests. The paper deals with the 1D/3D discussion in the framework of modelling activities related to integral passive facilities like PANDA. It starts choosing reference tests among those corresponding to our participation in TEPSS project. The discrepancies observed in a 1D simulation of the selected tests will be shown and analyzed. An evaluation of how the 3D version can lead to a better agreement with data will be included. Disadvantages of 3D codes will be shown too. Combining the use of different codes, and considering analyst criteria, will make possible to establish suitable recommendations from both engineering and scientific point of view. (author)

  8. Analysis of Flash Flood Routing by Means of 1D - Hydraulic Modelling

    Tesfay Abraha, Zerisenay

    2013-01-01

    This study was conducted at the mountainous catchment part of Batinah Region of the Sultanate of Oman called Al-Awabi watershed which is about 260km2 in area with about 40 Km long Wadi main channel. The study paper presents a proposed modeling approach and possible scenario analysis which uses 1D - hydraulic modeling for flood routing analysis; and the main tasks of this study work are (1) Model setup for Al-Awabi watershed area, (2) Sensitivity Analysis, and (3) Scenario Analysis on impacts ...

  9. Modeling of the diffraction pattern of 1D-disordered silicon carbide

    A method for calculating the diffraction pattern of a 1D-disordered crystal structure is considered by the example of silicon carbide. One-dimensional disordering is described using a cell setting the mutual position of all close-packed crystal layers. Two models of structure disordering during the polytypic transformation of the silicon carbide cubic modification into hexagonal are discussed. The results of the calculation of the diffraction spectrum in different stages of polytypic transformation are reported. It is shown that 1D disordering leads to the formation of a set of weak diffraction reflections. The experimentally observed changes in the diffraction pattern can be interpreted within the hypothesis on crystal structure disordering through displacement of adjacent close-packed layers.

  10. Verification and comparison of four numerical schemes for a 1D viscoelastic blood flow model.

    Wang, Xiaofei; Fullana, Jose-Maria; Lagrée, Pierre-Yves

    2015-01-01

    A reliable and fast numerical scheme is crucial for the 1D simulation of blood flow in compliant vessels. In this paper, a 1D blood flow model is incorporated with a Kelvin-Voigt viscoelastic arterial wall. This leads to a nonlinear hyperbolic-parabolic system, which is then solved with four numerical schemes, namely: MacCormack, Taylor-Galerkin, monotonic upwind scheme for conservation law and local discontinuous Galerkin. The numerical schemes are tested on a single vessel, a simple bifurcation and a network with 55 arteries. The numerical solutions are checked favorably against analytical, semi-analytical solutions or clinical observations. Among the numerical schemes, comparisons are made in four important aspects: accuracy, ability to capture shock-like phenomena, computational speed and implementation complexity. The suitable conditions for the application of each scheme are discussed. PMID:25145651

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

    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

  12. 1D-3D Hybrid Modelling - From Multi-Compartment Models to Full Resolution Models in Space and Time

    Stephan Grein

    2014-07-01

    Full Text Available Investigation of cellular and network dynamics in the brain by means of modeling & simulation has evolved into a highly interdisciplinary field, that uses sophisticated modeling & simulation approaches to understand distinct areas of brain function. Depending on the underlying complexity, these models vary in level of detail to cope with the attached computational cost. Hence for large network simulations, single neurons are typically reduced to time-dependent signal processors, dismissing spatial aspects of the cells. For single cell or small-world networks, general purpose simulators allow for space and time-dependent simulations of electrical signal processing, based on the cable equation theory. An emerging field in Computational Neuroscience encompasses a new level of detail by incorporating the 3D morphology of cells and organelles into 3D space and time-dependent simulations. Every approach has its advantages and limitations, such as computational cost, integrated and methods-spanning simulation approaches, depending on the network size could establish new ways to investigate the brain. We present a hybrid simulation approach, that makes use of reduced 1D-models using e.g. the NEURON which couples to fully resolved models for simulating cellular and sub-cellular dynamics, including the detailed 3D-morphology of neurons and organelles. To couple 1D- & 3D-simulations, we present a geometry and membrane potential mapping framework, with which graph-based morphologies, e.g. in swc-/hoc-format, are mapped to full surface and volume representations of the neuron; membrane potential data from 1D-simulations are used as boundary conditions for full 3D simulations. Thus, established models and data, based on general purpose 1D-simulators, can be directly coupled to the emerging field of fully resolved highly detailed 3D-modeling approaches. The new framework is applied to investigate electrically active neurons and their intracellular spatio

  13. Numerical modeling of 1-D transient poroelastic waves in the low-frequency range

    Chiavassa, Guillaume; Piraux, Joël

    2007-01-01

    Propagation of transient mechanical waves in porous media is numerically investigated in 1D. The framework is the linear Biot's model with frequency-independant coefficients. The coexistence of a propagating fast wave and a diffusive slow wave makes numerical modeling tricky. A method combining three numerical tools is proposed: a fourth-order ADER scheme with time-splitting to deal with the time-marching, a space-time mesh refinement to account for the small-scale evolution of the slow wave, and an interface method to incorporate the jump conditions at interfaces. Comparisons with analytical solutions confirm the validity of this approach.

  14. GE SBWR stability analysis using TRAC-BF1 1-D kinetics model

    GE's simplified boiling water reactor, with its unique feature of using natural circulation to remove the heat from the reactor core, is a complicated dynamic system. Previous work by authors using the TRAC-BF1 code and a point kinetics model predicted that an SBWR may experience large amplitude power oscillation under certain low pressure and high power operating conditions. To further confirm the existence of this power oscillation and explore the dynamic spatial reactor power distribution, the TRAC-BF1 1-D kinetics model was used. The results show that an instability exists and the power oscillation starting time and maximum peak power are different from the point kinetics results

  15. Prediction of car cabin environment by means of 1D and 3D cabin model

    Fišer, J.; Pokorný, J.; Jícha, M.

    2012-04-01

    Thermal comfort and also reduction of energy requirements of air-conditioning system in vehicle cabins are currently very intensively investigated and up-to-date issues. The article deals with two approaches of modelling of car cabin environment; the first model was created in simulation language Modelica (typical 1D approach without cabin geometry) and the second one was created in specialized software Theseus-FE (3D approach with cabin geometry). Performance and capabilities of this tools are demonstrated on the example of the car cabin and the results from simulations are compared with the results from the real car cabin climate chamber measurements.

  16. Prediction of car cabin environment by means of 1D and 3D cabin model

    Jícha M.

    2012-04-01

    Full Text Available Thermal comfort and also reduction of energy requirements of air-conditioning system in vehicle cabins are currently very intensively investigated and up-to-date issues. The article deals with two approaches of modelling of car cabin environment; the first model was created in simulation language Modelica (typical 1D approach without cabin geometry and the second one was created in specialized software Theseus-FE (3D approach with cabin geometry. Performance and capabilities of this tools are demonstrated on the example of the car cabin and the results from simulations are compared with the results from the real car cabin climate chamber measurements.

  17. Assessing the impact of different sources of topographic data on 1-D hydraulic modelling of floods

    Ali, A. Md; Solomatine, D. P.; Di Baldassarre, G.

    2015-01-01

    Topographic data, such as digital elevation models (DEMs), are essential input in flood inundation modelling. DEMs can be derived from several sources either through remote sensing techniques (spaceborne or airborne imagery) or from traditional methods (ground survey). The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), the Shuttle Radar Topography Mission (SRTM), the light detection and ranging (lidar), and topographic contour maps are some of the most commonly used sources of data for DEMs. These DEMs are characterized by different precision and accuracy. On the one hand, the spatial resolution of low-cost DEMs from satellite imagery, such as ASTER and SRTM, is rather coarse (around 30 to 90 m). On the other hand, the lidar technique is able to produce high-resolution DEMs (at around 1 m), but at a much higher cost. Lastly, contour mapping based on ground survey is time consuming, particularly for higher scales, and may not be possible for some remote areas. The use of these different sources of DEM obviously affects the results of flood inundation models. This paper shows and compares a number of 1-D hydraulic models developed using HEC-RAS as model code and the aforementioned sources of DEM as geometric input. To test model selection, the outcomes of the 1-D models were also compared, in terms of flood water levels, to the results of 2-D models (LISFLOOD-FP). The study was carried out on a reach of the Johor River, in Malaysia. The effect of the different sources of DEMs (and different resolutions) was investigated by considering the performance of the hydraulic models in simulating flood water levels as well as inundation maps. The outcomes of our study show that the use of different DEMs has serious implications to the results of hydraulic models. The outcomes also indicate that the loss of model accuracy due to re-sampling the highest resolution DEM (i.e. lidar 1 m) to lower resolution is much less than the loss of model accuracy due

  18. Optimal modeling of 1D azimuth correlations in the context of Bayesian inference

    De Kock, Michiel B; Trainor, Thomas A

    2015-01-01

    Analysis and interpretation of spectrum and correlation data from high-energy nuclear collisions is currently controversial because two opposing physics narratives derive contradictory implications from the same data-one narrative claiming collision dynamics is dominated by dijet production and projectile-nucleon fragmentation, the other claiming collision dynamics is dominated by a dense, flowing QCD medium. Opposing interpretations seem to be supported by alternative data models, and current model-comparison schemes are unable to distinguish between them. There is clearly need for a convincing new methodology to break the deadlock. In this study we introduce Bayesian Inference (BI) methods applied to angular correlation data as a basis to evaluate competing data models. For simplicity the data considered are projections of 2D angular correlations onto 1D azimuth from three centrality classes of 200 GeV Au-Au collisions. We consider several data models typical of current model choices, including Fourier seri...

  19. Microscopic approach to a class of 1D quantum critical models

    Starting from the finite volume form factors of local operators, we show how and under which hypothesis the c = 1 free boson conformal field theory in two-dimensions emerges as an effective theory governing the large-distance regime of multi-point correlation functions for a large class of one-dimensional massless quantum Hamiltonians. In our approach, in the large-distance critical regime, the local operators of the initial model are represented by well suited vertex operators associated to the free boson model. This provides an effective field theoretic description of the large distance behaviour of correlation functions in 1D quantum critical models. We develop this description starting from the first principles and directly at the microscopic level, namely in terms of the properties of the finite volume matrix elements of local operators. (paper)

  20. A 1-D evolutionary model for icy satellites, applied to Enceladus

    Prialnik, Uri Malamud Dina

    2015-01-01

    We develop a long-term 1-D evolution model for icy satellites that couples multiple processes: water migration and differentiation, geochemical reactions and silicate phase transitions, compaction by self-gravity, and ablation. The model further considers the following energy sources and sinks: tidal heating, radiogenic heating, geochemical energy released by serpentinization or absorbed by mineral dehydration, gravitational energy and insolation, and heat transport by conduction, convection, and advection. We apply the model to Enceladus, by guessing the initial conditions that would render a structure compatible with present-day observations, assuming the initial structure to have been homogeneous. Assuming the satellite has been losing water continually along its evolution, we postulate that it was formed as a more massive, more icy and more porous satellite, and gradually transformed into its present day state due to sustained long-term tidal heating. We consider several initial compositions and evolution...

  1. REAL-TIME FLOOD FORECASTING MODELING OF 1D UNSTEADY CHANNEL FLOW AND KALMAN FILTER

    2001-01-01

    The model of 1D unsteady channel flow combined with the Kalmanfilter for real-time channel flood forecasting was attempted in this study. The suitable upstream and downstream boundary conditions were suggested. The system equation was given by the linearization of the finitedifference equations of the mass conservation and momentum equations as well as the boundary conditions. In the Kalman filter updating model, because the number of measurement variable is less then that of state-space variables, the measurement error covariance matrix could be estimated in real time through the innovation sequence, and the system error covariance matrix needs to be estimated preliminarily. A real example of flood forecasting in the Huaihe River was given to explain how the method works. The results show that the model is reasonable and effective.

  2. A 1-D modelling of streaming potential dependence on water content during drainage experiment in sand

    Allègre, Vincent; Ackerer, Philippe; Jouniaux, Laurence; Sailhac, Pascal; 10.1111/j.1365-246X.2012.05371.x

    2012-01-01

    The understanding of electrokinetics for unsaturated conditions is crucial for numerous of geophysical data interpretation. Nevertheless, the behaviour of the streaming potential coefficient C as a function of the water saturation Sw is still discussed. We propose here to model both the Richards' equation for hydrodynamics and the Poisson's equation for electrical potential for unsaturated conditions using 1-D finite element method. The equations are first presented and the numerical scheme is then detailed for the Poisson's equation. Then, computed streaming potentials (SPs) are compared to recently published SP measurements carried out during drainage experiment in a sand column. We show that the apparent measurement of DV / DP for the dipoles can provide the SP coefficient in these conditions. Two tests have been performed using existing models for the SP coefficient and a third one using a new relation. The results show that existing models of unsaturated SP coefficients C(Sw) provide poor results in term...

  3. Survey of Multi-Material Closure Models in 1D Lagrangian Hydrodynamics

    Maeng, Jungyeoul Brad [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Hyde, David Andrew Bulloch [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

    2015-07-28

    Accurately treating the coupled sub-cell thermodynamics of computational cells containing multiple materials is an inevitable problem in hydrodynamics simulations, whether due to initial configurations or evolutions of the materials and computational mesh. When solving the hydrodynamics equations within a multi-material cell, we make the assumption of a single velocity field for the entire computational domain, which necessitates the addition of a closure model to attempt to resolve the behavior of the multi-material cells’ constituents. In conjunction with a 1D Lagrangian hydrodynamics code, we present a variety of both the popular as well as more recently proposed multi-material closure models and survey their performances across a spectrum of examples. We consider standard verification tests as well as practical examples using combinations of fluid, solid, and composite constituents within multi-material mixtures. Our survey provides insights into the advantages and disadvantages of various multi-material closure models in different problem configurations.

  4. EFDC1D - A ONE DIMENSIONAL HYDRODYNAMIC AND SEDIMENT TRANSPORT MODEL FOR RIVER AND STREAM NETWORKS: MODEL THEORY AND USERS GUIDE

    This technical report describes the new one-dimensional (1D) hydrodynamic and sediment transport model EFDC1D. This model that can be applied to stream networks. The model code and two sample data sets are included on the distribution CD. EFDC1D can simulate bi-directional unstea...

  5. Fluid friction and wall viscosity of the 1D blood flow model.

    Wang, Xiao-Fei; Nishi, Shohei; Matsukawa, Mami; Ghigo, Arthur; Lagrée, Pierre-Yves; Fullana, Jose-Maria

    2016-02-29

    We study the behavior of the pulse waves of water into a flexible tube for application to blood flow simulations. In pulse waves both fluid friction and wall viscosity are damping factors, and difficult to evaluate separately. In this paper, the coefficients of fluid friction and wall viscosity are estimated by fitting a nonlinear 1D flow model to experimental data. In the experimental setup, a distensible tube is connected to a piston pump at one end and closed at another end. The pressure and wall displacements are measured simultaneously. A good agreement between model predictions and experiments was achieved. For amplitude decrease, the effect of wall viscosity on the pulse wave has been shown as important as that of fluid viscosity. PMID:26862041

  6. Fluid friction and wall viscosity of the 1D blood flow model

    Wang, Xiao-Fei; Matsukawa, Mami; Ghigo, Arthur; Lagrée, Pierre-Yves; Fullana, Jose-Maria

    2015-01-01

    We study the behavior of the pulse waves of water into a flexible tube for application to blood flow simulations. In pulse waves both fluid friction and wall viscosity are damping factors, and difficult to evaluate separately. In this paper, the coefficients of fluid friction and wall viscosity are estimated by fitting a nonlinear 1D flow model to experimental data. In the experimental setup, a distensible tube is connected to a piston pump at one end and closed at another end. The pressure and wall displacements are measured simultaneously. A good agreement between model predictions and experiments was achieved. For amplitude decrease, the effect of wall viscosity on the pulse wave has been shown as important as that of fluid viscosity.

  7. Full Waveform 3D Synthetic Seismic Algorithm for 1D Layered Anelastic Models

    Schwaiger, H. F.; Aldridge, D. F.; Haney, M. M.

    2007-12-01

    Numerical calculation of synthetic seismograms for 1D layered earth models remains a significant aspect of amplitude-offset investigations, surface wave studies, microseismic event location approaches, and reflection interpretation or inversion processes. Compared to 3D finite-difference algorithms, memory demand and execution time are greatly reduced, enabling rapid generation of seismic data within workstation or laptop computational environments. We have developed a frequency-wavenumber forward modeling algorithm adapted to realistic 1D geologic media, for the purpose of calculating seismograms accurately and efficiently. The earth model consists of N layers bounded by two halfspaces. Each layer/halfspace is a homogeneous and isotropic anelastic (attenuative and dispersive) solid, characterized by a rectangular relaxation spectrum of absorption mechanisms. Compressional and shear phase speeds and quality factors are specified at a particular reference frequency. Solution methodology involves 3D Fourier transforming the three coupled, second- order, integro-differential equations for particle displacements to the frequency-horizontal wavenumber domain. An analytic solution of the resulting ordinary differential system is obtained. Imposition of welded interface conditions (continuity of displacement and stress) at all interfaces, as well as radiation conditions in the two halfspaces, yields a system of 6(N+1) linear algebraic equations for the coefficients in the ODE solution. An optimized inverse 2D Fourier transform to the space domain gives the seismic wavefield on a horizontal plane. Finally, three-component seismograms are obtained by accumulating frequency spectra at designated receiver positions on this plane, followed by a 1D inverse FFT from angular frequency ω to time. Stress-free conditions may be applied at the top or bottom interfaces, and seismic waves are initiated by force or moment density sources. Examples reveal that including attenuation

  8. Initial Stage of the Microwave Ionization Wave Within a 1D Model

    Semenov, V. E.; Rakova, E. I.; Glyavin, M. Yu.; Nusinovich, G. S.

    2016-06-01

    The dynamics of the microwave breakdown in a gas is simulated numerically within a simple 1D model which takes into account such processes as the impact ionization of gas molecules, the attachment of electrons to neutral molecules, and plasma diffusion. Calculations are carried out for different spatial distributions of seed electrons with account for reflection of the incident electromagnetic wave from the plasma. The results reveal considerable dependence of the ionization wave evolution on the relation between the field frequency and gas pressure, as well as on the existence of extended rarefied halo of seed electrons. At relatively low gas pressures (or high field frequencies), the breakdown process is accompanied by the stationary ionization wave moving towards the incident electromagnetic wave. In the case of a high gas pressure (or a relatively low field frequency), the peculiarities of the breakdown are associated with the formation of repetitive jumps of the ionization front.

  9. Study of phase space structures in driven 1D Vlasov poisson model

    Electrostatic waves in a collisionless, unmagnetized plasma are known to interact with particles that stream with velocities close to the wave phase speed to produce damping effects, particle trapping and interesting nonlinear coherent structures. For example, it is well known that if the initial amplitude of the wave is large enough, the damping effects can be overcome to form BGK structures. In the present work, we consider a 1D driven Vlasov-Poisson plasma model. It is demonstrated that by a careful choice of drive phase and for drive amplitudes smaller than or comparable to the linear limit, it is possible to generate surprisingly large amplitude coherent structures in phase space. This and other details will be presented. (author)

  10. Significance of flow clustering and sequencing on sediment transport: 1D sediment transport modelling

    Hassan, Kazi; Allen, Deonie; Haynes, Heather

    2016-04-01

    This paper considers 1D hydraulic model data on the effect of high flow clusters and sequencing on sediment transport. Using observed flow gauge data from the River Caldew, England, a novel stochastic modelling approach was developed in order to create alternative 50 year flow sequences. Whilst the observed probability density of gauge data was preserved in all sequences, the order in which those flows occurred was varied using the output from a Hidden Markov Model (HMM) with generalised Pareto distribution (GP). In total, one hundred 50 year synthetic flow series were generated and used as the inflow boundary conditions for individual flow series model runs using the 1D sediment transport model HEC-RAS. The model routed graded sediment through the case study river reach to define the long-term morphological changes. Comparison of individual simulations provided a detailed understanding of the sensitivity of channel capacity to flow sequence. Specifically, each 50 year synthetic flow sequence was analysed using a 3-month, 6-month or 12-month rolling window approach and classified for clusters in peak discharge. As a cluster is described as a temporal grouping of flow events above a specified threshold, the threshold condition used herein is considered as a morphologically active channel forming discharge event. Thus, clusters were identified for peak discharges in excess of 10%, 20%, 50%, 100% and 150% of the 1 year Return Period (RP) event. The window of above-peak flows also required cluster definition and was tested for timeframes 1, 2, 10 and 30 days. Subsequently, clusters could be described in terms of the number of events, maximum peak flow discharge, cumulative flow discharge and skewness (i.e. a description of the flow sequence). The model output for each cluster was analysed for the cumulative flow volume and cumulative sediment transport (mass). This was then compared to the total sediment transport of a single flow event of equivalent flow volume

  11. Monte Carlo Uncertainty Quantification Using Quasi-1D SRM Ballistic Model

    Davide Viganò

    2016-01-01

    Full Text Available Compactness, reliability, readiness, and construction simplicity of solid rocket motors make them very appealing for commercial launcher missions and embarked systems. Solid propulsion grants high thrust-to-weight ratio, high volumetric specific impulse, and a Technology Readiness Level of 9. However, solid rocket systems are missing any throttling capability at run-time, since pressure-time evolution is defined at the design phase. This lack of mission flexibility makes their missions sensitive to deviations of performance from nominal behavior. For this reason, the reliability of predictions and reproducibility of performances represent a primary goal in this field. This paper presents an analysis of SRM performance uncertainties throughout the implementation of a quasi-1D numerical model of motor internal ballistics based on Shapiro’s equations. The code is coupled with a Monte Carlo algorithm to evaluate statistics and propagation of some peculiar uncertainties from design data to rocker performance parameters. The model has been set for the reproduction of a small-scale rocket motor, discussing a set of parametric investigations on uncertainty propagation across the ballistic model.

  12. A 1-D evolutionary model for icy satellites, applied to Enceladus

    Malamud, Uri; Prialnik, Dina

    2016-04-01

    We develop a long-term 1-D evolution model for icy satellites that couples multiple processes: water migration and differentiation, geochemical reactions and silicate phase transitions, compaction by self-gravity, and ablation. The model further considers the following energy sources and sinks: tidal heating, radiogenic heating, geochemical energy released by serpentinization or absorbed by mineral dehydration, gravitational energy and insolation, and heat transport by conduction, convection, and advection. We apply the model to Enceladus, by guessing the initial conditions that would render a structure compatible with present-day observations, assuming the initial structure to have been homogeneous. Assuming the satellite has been losing water continually along its evolution, we postulate that it was formed as a more massive, more icy and more porous satellite, and gradually transformed into its present day state due to sustained long-term tidal heating. We consider several initial compositions and evolution scenarios and follow the evolution for the age of the Solar System, testing the present day model results against the available observational constraints. Our model shows the present configuration to be differentiated into a pure icy mantle, several tens of km thick, overlying a rocky core, composed of dehydrated rock at the center and hydrated rock in the outer part. For Enceladus, it predicts a higher rock/ice mass ratio than previously assumed and a thinner ice mantle, compatible with recent estimates based on gravity field measurements. Although, obviously, the model cannot be used to explain local phenomena, it sheds light on the internal structure invoked in explanations of localized features and activities.

  13. Modelling land surface - atmosphere interactions

    Rasmussen, Søren Højmark

    The study is investigates modelling of land surface – atmosphere interactions in context of fully coupled climatehydrological model. With a special focus of under what condition a fully coupled model system is needed. Regional climate model inter-comparison projects as ENSEMBLES have shown bias...

  14. 1D Tight-Binding Models Render Quantum First Passage Time "Speakable"

    Ranjith, V.; Kumar, N.

    2015-12-01

    The calculation of First Passage Time (moreover, even its probability density in time) has so far been generally viewed as an ill-posed problem in the domain of quantum mechanics. The reasons can be summarily seen in the fact that the quantum probabilities in general do not satisfy the Kolmogorov sum rule: the probabilities for entering and non-entering of Feynman paths into a given region of space-time do not in general add up to unity, much owing to the interference of alternative paths. In the present work, it is pointed out that a special case exists (within quantum framework), in which, by design, there exists one and only one available path (i.e., door-way) to mediate the (first) passage -no alternative path to interfere with. Further, it is identified that a popular family of quantum systems - namely the 1d tight binding Hamiltonian systems - falls under this special category. For these model quantum systems, the first passage time distributions are obtained analytically by suitably applying a method originally devised for classical (stochastic) mechanics (by Schroedinger in 1915). This result is interesting especially given the fact that the tight binding models are extensively used in describing everyday phenomena in condense matter physics.

  15. Topological order in 1D super-lattice Bose-Hubbard models

    Fleischhauer, Michael; Grusdt, Fabian; Hoening, Michael

    2013-05-01

    After the discovery of topological insulators as a new state of matter and their consequent classification for free fermions, the question arises what kind of topological order can be supported by incompressible systems of interacting bosons. We consider a 1D super-lattice Hamiltonian with a non-trivial band structure (the Su-Schrieffer-Heeger model) and show that its Mott-insulating (MI) states can be classified by a quantized many-body winding number. This quantization is protected by sub-lattice and time-reversal symmetries, and it allows the implementation of a quantized cyclic pumping process (Thouless pump) in a simple super-lattice Bose-Hubbard model (BHM). For extended BHMs we discuss a connection of such a pump with the fractional quantum Hall effect. Furthermore we show that the quantization of the winding number leads to localized, protected edge states at sharp interfaces between topologically distinct MI phases which can be experimentally realized using Bose-Fermi mixtures in optical superlattices. DMRG simulations show that these edge states manifest themself either in localized density maxima or localized density minima, which can easily be detected. Supported by research center OPTIMAS and graduate school MAINZ.

  16. An Analytic Radiative-Convective Model for Planetary Atmospheres

    Robinson, Tyler D; 10.1088/0004-637X/757/1/104

    2012-01-01

    We present an analytic 1-D radiative-convective model of the thermal structure of planetary atmospheres. Our model assumes that thermal radiative transfer is gray and can be represented by the two-stream approximation. Model atmospheres are assumed to be in hydrostatic equilibrium, with a power law scaling between the atmospheric pressure and the gray thermal optical depth. The convective portions of our models are taken to follow adiabats that account for condensation of volatiles through a scaling parameter to the dry adiabat. By combining these assumptions, we produce simple, analytic expressions that allow calculations of the atmospheric pressure-temperature profile, as well as expressions for the profiles of thermal radiative flux and convective flux. We explore the general behaviors of our model. These investigations encompass (1) worlds where atmospheric attenuation of sunlight is weak, which we show tend to have relatively high radiative-convective boundaries, (2) worlds with some attenuation of sunli...

  17. Modelling hydrology of a single bioretention system with HYDRUS-1D.

    Meng, Yingying; Wang, Huixiao; Chen, Jiangang; Zhang, Shuhan

    2014-01-01

    A study was carried out on the effectiveness of bioretention systems to abate stormwater using computer simulation. The hydrologic performance was simulated for two bioretention cells using HYDRUS-1D, and the simulation results were verified by field data of nearly four years. Using the validated model, the optimization of design parameters of rainfall return period, filter media depth and type, and surface area was discussed. And the annual hydrologic performance of bioretention systems was further analyzed under the optimized parameters. The study reveals that bioretention systems with underdrains and impervious boundaries do have some detention capability, while their total water retention capability is extremely limited. Better detention capability is noted for smaller rainfall events, deeper filter media, and design storms with a return period smaller than 2 years, and a cost-effective filter media depth is recommended in bioretention design. Better hydrologic effectiveness is achieved with a higher hydraulic conductivity and ratio of the bioretention surface area to the catchment area, and filter media whose conductivity is between the conductivity of loamy sand and sandy loam, and a surface area of 10% of the catchment area is recommended. In the long-term simulation, both infiltration volume and evapotranspiration are critical for the total rainfall treatment in bioretention systems. PMID:25133240

  18. Modeling Soil Salt and Nitrogen Transport under Different Fertigation Practices with Hydrus-1D

    Zeng Wen-zhi

    2013-05-01

    Full Text Available In this study the effects of different fertigation practices on salt and nitrogen dynamics were analyzed in the Hetao District, China by using the Hydrus-1D model. The results indicated that the soil electrical conductivity increased gradually with depth after irrigation and the electrical conductivity of 0~60 cm depth changed faster than that of 60~100 cm depth. However, the soil ammonium nitrogen concentration decreased with depth and high irrigation intensity could promote the increase of ammonium nitrogen while reducing the differences of their distributions in soil profile. In addition, when the initial urea application was in a small amount (10 g, the nitrate nitrogen concentration increased with soil depth while decreased with irrigation intensity after irrigation. Furthermore, both ammonium and nitrate nitrogen content of soil profile rose with increasing initial urea application amount, which played a more important role in the changes of soil ammonium (0~100 cm and nitrate (0~80 cm nitrogen content than irrigation intensity.

  19. Self-assembling morphologies in a 1D model of two-inclusion-containing lipid membranes

    Zhou, Ling; Cheng, Mingfei; Fang, Jinghuai; Peng, Ju

    2016-08-01

    The self-assembling morphologies in a 1D model of two-inclusion-containing lipid membranes are investigated by using self-consistent field theory. It is found that the shape and overall volume fraction of lipids, the hydrophobic strength and the distance of inclusions play important roles in the morphology of lipid membrane. The membrane consisting of cylindrical lipids with a symmetrical head and tail only forms the well-known normal morphology. However, for the membrane consisting of cone-like lipids with a relatively big head, the increase of the hydrophobic strength of inclusions can realize the membrane transition from the normal morphology to the pore morphologies. With increasing distance between two inclusions, two pores, three pores and four pores appear in turn. Conversely, the increase of the overall volume fraction of lipids can make the membrane undergo a reentrant transition from pore morphologies to normal morphologies. The results may be helpful in our understanding of the pore-forming mechanism.

  20. Results and limits in the 1-D analytical modeling for the asymmetric DG SOI MOSFET

    O. Cobianu

    2008-05-01

    Full Text Available This paper presents the results and the limits of 1-D analytical modeling of electrostatic potential in the low-doped p type silicon body of the asymmetric n-channel DG SOI MOSFET, where the contribution to the asymmetry comes only from p- and n-type doping of polysilicon used as the gate electrodes. Solving Poisson's equation with boundary conditions based on the continuity of normal electrical displacement at interfaces and the presence of a minimum electrostatic potential by using the Matlab code we have obtained a minimum potential with a slow variation in the central zone of silicon with the value pinned around 0.46 V, where the applied VGS voltage varies from 0.45 V to 0.95 V. The paper states clearly the validity domain of the analytical solution and the important effect of the localization of the minimum electrostatic potential value on the potential variation at interfaces as a function of the applied VGS voltage.

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

    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

  2. Column Testing and 1D Reactive Transport Modeling to Evaluate Uranium Plume Persistence Processes

    Johnson, R. H.; Morrison, S.; Morris, S.; Tigar, A.; Dam, W. L.; Dayvault, J.

    2015-12-01

    At many U.S. Department of Energy Office of Legacy Management sites, 100 year natural flushing was selected as a remedial option for groundwater uranium plumes. However, current data indicate that natural flushing is not occurring as quickly as expected and solid-phase and aqueous uranium concentrations are persistent. At the Grand Junction, Colorado office site, column testing was completed on core collected below an area where uranium mill tailings have been removed. The total uranium concentration in this core was 13.2 mg/kg and the column was flushed with laboratory-created water with no uranium and chemistry similar to the nearby Gunnison River. The core was flushed for a total of 91 pore volumes producing a maximum effluent uranium concentration of 6,110 μg/L at 2.1 pore volumes and a minimum uranium concentration of 36.2 μg/L at the final pore volume. These results indicate complex geochemical reactions at small pore volumes and a long tailing affect at greater pore volumes. Stop flow data indicate the occurrence of non-equilibrium processes that create uranium concentration rebound. These data confirm the potential for plume persistence, which is occurring at the field scale. 1D reactive transport modeling was completed using PHREEQC (geochemical model) and calibrated to the column test data manually and using PEST (inverse modeling calibration routine). Processes of sorption, dual porosity with diffusion, mineral dissolution, dispersion, and cation exchange were evaluated separately and in combination. The calibration results indicate that sorption and dual porosity are major processes in explaining the column test data. These processes are also supported by fission track photographs that show solid-phase uranium residing in less mobile pore spaces. These procedures provide valuable information on plume persistence and secondary source processes that may be used to better inform and evaluate remedial strategies, including natural flushing.

  3. Titan atmospheric models intercomparison

    Pernot, P.

    2008-09-01

    Several groups over the world have developed independently models of the photochemistry of Titan. The Cassini mission reveals daily that the chemical complexity is beyond our expectations e. g. observation of heavy positive and negative ions..., and the models are updated accordingly. At this stage, there is no consensus on the various input parameters, and it becomes increasingly difficult to compare outputs form different models. An ISSI team of experts of those models will be gathered shortly to proceed to an intercomparison, i.e. to assess how the models behave, given identical sets of inputs (collectively defined). Expected discrepancies will have to be elucidated and reduced. This intercomparison will also be an occasion to estimate explicitly the importance of various physicalchemical processes on model predictions versus observations. More robust and validated models are expected from this study for the interpretation of Titanrelated data.

  4. 1D modelling of alteration of compacted bentonite due to cementitious leachate

    Document available in extended abstract form only. The interaction between the bentonite buffer and alkaline fluids deriving from cementitious materials used in the construction and operation of a KBS-3 spent-fuel repository is considered as a risk to the long-term chemical stability of bentonite. The associated long-term safety concerns include a loss of bentonite swelling pressure, increased hydraulic conductivity, and possibly fracturing of bentonite due to cementation. Consequently, it is important to identify and quantify the mechanisms controlling the degradation of the buffer by a cement leachate. This would subsequently help take adequate measures to mitigate the risk posed by alkaline attack. In this work, 1D reactive transport model (TOUGHREACT v. 1.0) calculations have been carried out to shed light on the possible extent of MX-80 bentonite alteration due to three leachates envisaged to represent different stages of cement degradation and subsequent mixing with the ground water. The reaction time for the calculations was set at 10,000 years or until the porosity was clogged in the compacted bentonite, whichever occurred earlier. A common feature of the model outcome for each leachate is the total consumption of gypsum and calcite in the bentonite and the clogging of the pore space by secondary mineral phases at the interface between the buffer and the rock fracture. Depending on the pH (12.17, 11.60 or 9.70), this clogging was calculated to occur after 10, 18 and 5,900 years of interaction, respectively. For each case, the heavily altered zone in bentonite is typically located within 1 cm from the interface, although for the lowest pH, the alteration extends slightly further into the bentonite due to the longer time it takes to clog the porosity. Overall, the model results are found sensitive to the chemical composition of the leachate interacting with the buffer. Regardless of the case studied, there is a significant smectite fraction left at the end of

  5. Simulation of a semiarid stream flow using the 1D model (Rubarbe) : case of Ksob Wadi in Algeria

    Hasbaia, M.; Adoui, H.; Paquier, A.

    2015-01-01

    This study aims to show the ability of 1D model to simulate the wadi flow during the flood. The wadi is a semiarid stream characterized by the non uniform geometry and an intermittent flow. These complexities can be modelled by a robust 1D model such as Rubarbe. In this article we study the 1994 flood of Ksob wadi that caused a large inundation of M’sila town in the center of Algeria. The simulated reach located downstream the Ksob dam, crosses M’sila over a length of 6800 m w...

  6. A Global Atmospheric Model of Meteoric Iron

    Feng, Wuhu; Marsh, Daniel R.; Chipperfield, Martyn P.; Janches, Diego; Hoffner, Josef; Yi, Fan; Plane, John M. C.

    2013-01-01

    The first global model of meteoric iron in the atmosphere (WACCM-Fe) has been developed by combining three components: the Whole Atmosphere Community Climate Model (WACCM), a description of the neutral and ion-molecule chemistry of iron in the mesosphere and lower thermosphere (MLT), and a treatment of the injection of meteoric constituents into the atmosphere. The iron chemistry treats seven neutral and four ionized iron containing species with 30 neutral and ion-molecule reactions. The meteoric input function (MIF), which describes the injection of Fe as a function of height, latitude, and day, is precalculated from an astronomical model coupled to a chemical meteoric ablation model (CABMOD). This newly developed WACCM-Fe model has been evaluated against a number of available ground-based lidar observations and performs well in simulating the mesospheric atomic Fe layer. The model reproduces the strong positive correlation of temperature and Fe density around the Fe layer peak and the large anticorrelation around 100 km. The diurnal tide has a significant effect in the middle of the layer, and the model also captures well the observed seasonal variations. However, the model overestimates the peak Fe+ concentration compared with the limited rocket-borne mass spectrometer data available, although good agreement on the ion layer underside can be obtained by adjusting the rate coefficients for dissociative recombination of Fe-molecular ions with electrons. Sensitivity experiments with the same chemistry in a 1-D model are used to highlight significant remaining uncertainties in reaction rate coefficients, and to explore the dependence of the total Fe abundance on the MIF and rate of vertical transport.

  7. Diesel Engine performance improvement in a 1-D engine model using Particle Swarm Optimization

    Karra, Prashanth

    2015-12-01

    A particle swarm optimization (PSO) technique was implemented to improve the engine development and optimization process to simultaneously reduce emissions and improve the fuel efficiency. The optimization was performed on a 4-stroke 4-cylinder GT-Power based 1-D diesel engine model. To achieve the multi-objective optimization, a merit function was defined which included the parameters to be optimized: Nitrogen Oxides (NOx), Nonmethyl hydro carbons (NMHC), Carbon Monoxide (CO), Brake Specific Fuel Consumption (BSFC). EPA Tier 3 emissions standards for non-road diesel engines between 37 and 75 kW of output were chosen as targets for the optimization. The combustion parameters analyzed in this study include: Start of main Injection, Start of Pilot Injection, Pilot fuel quantity, Swirl, and Tumble. The PSO was found to be very effective in quickly arriving at a solution that met the target criteria as defined in the merit function. The optimization took around 40-50 runs to find the most favourable engine operating condition under the constraints specified in the optimization. In a favourable case with a high merit function values, the NOx+NMHC and CO values were reduced to as low as 2.9 and 0.014 g/kWh, respectively. The operating conditions at this point were: 10 ATDC Main SOI, -25 ATDC Pilot SOI, 0.25 mg of pilot fuel, 0.45 Swirl and 0.85 tumble. These results indicate that late main injections preceded by a close, small pilot injection are most favourable conditions at the operating condition tested.

  8. MARCS model atmospheres

    Plez, Bertrand

    2008-01-01

    In this review presented at the Symposium A stellar journey in Uppsala, June 2008, I give my account of the historical development of the MARCS code from the first version published in 1975 and its premises to the 2008 grid. It is shown that the primary driver for the development team is the science that can be done with the models, and that they constantly strive to include the best possible physical data. A few preliminary comparisons of M star model spectra to spectrophotometric observations are presented. Particular results related to opacity effects are discussed. The size of errors in the spectral energy distribution (SED) and model thermal stratification are estimated for different densities of the wavelength sampling. The number of points used in the MARCS 2008 grid (108000) is large enough to ensure errors of only a few K in all models of the grid, except the optically very thin layers of metal-poor stars. Errors in SEDs may reach about 10% locally in the UV. The published sampled SEDs are thus appro...

  9. Comparison of 1D and 2D CSR Models with Application to the FERMI@ELETTRA Bunch Compressors

    Bassi, G.; Ellison, J.A.; Heinemann, K.

    2011-03-28

    We compare our 2D mean field (Vlasov-Maxwell) treatment of coherent synchrotron radiation (CSR) effects with 1D approximations of the CSR force which are commonly implemented in CSR codes. In our model we track particles in 4D phase space and calculate 2D forces [1]. The major cost in our calculation is the computation of the 2D force. To speed up the computation and improve 1D models we also investigate approximations to our exact 2D force. As an application, we present numerical results for the Fermi{at}Elettra first bunch compressor with the configuration described in [1].

  10. Comparison of 1D and 2D CSR Models with Application to the FERMI(at)ELETTRA Bunch Compressors

    We compare our 2D mean field (Vlasov-Maxwell) treatment of coherent synchrotron radiation (CSR) effects with 1D approximations of the CSR force which are commonly implemented in CSR codes. In our model we track particles in 4D phase space and calculate 2D forces (1). The major cost in our calculation is the computation of the 2D force. To speed up the computation and improve 1D models we also investigate approximations to our exact 2D force. As an application, we present numerical results for the Fermi(at)Elettra first bunch compressor with the configuration described in (1).

  11. 1D Runoff-runon stochastic model in the light of queueing theory : heterogeneity and connectivity

    Harel, M.-A.; Mouche, E.; Ledoux, E.

    2012-04-01

    Runoff production on a hillslope during a rainfall event may be simplified as follows. Given a soil of constant infiltrability I, which is the maximum amount of water that the soil can infiltrate, and a constant rainfall intensity R, runoff is observed where R is greater than I. The infiltration rate equals the infiltrability when runoff is produced, R otherwise. When ponding time, topography, and overall spatial and temporal variations of physical parameters, such as R and I, are neglected, the runoff equation remains simple. In this study, we consider soils of spatially variable infiltrability. As runoff can re-infiltrate on down-slope areas of higher infiltrabilities (runon), the resulting process is highly non-linear. The stationary runoff equation is: Qn+1 = max(Qn + (R - In)*Δx , 0) where Qn is the runoff arriving on pixel n of size Δx [L2/T], R and In the rainfall intensity and infiltrability on that same pixel [L/T]. The non-linearity is due to the dependence of infiltration on R and Qn, that is runon. This re-infiltration process generates patterns of runoff along the slope, patterns that organise and connect to each other differently depending on the rainfall intensity and the nature of the soil heterogeneity. The runoff connectivity, assessed using the connectivity function of Allard (1993), affects greatly the dynamics of the runoff hillslope. Our aim is to assess, in a stochastic framework, the runoff organization on 1D slopes with random infiltrabilities (log-normal, exponential, bimodal and uniform distributions) by means of theoretical developments and numerical simulations. This means linking the nature of soil heterogeneity with the resulting runoff organisation. In term of connectivity, we investigate the relations between structural (infiltrability) and functional (runoff) connectivity. A theoretical framework based on the queueing theory is developed. We implement the idea of Jones et al. (2009), who remarked that the above formulation is

  12. A One-Dimensional (1-D) Three-Region Model for a Bubbling Fluidized-Bed Adsorber

    Lee, Andrew; Miller, David C.

    2012-01-01

    A general one-dimensional (1-D), three-region model for a bubbling fluidized-bed adsorber with internal heat exchangers has been developed. The model can predict the hydrodynamics of the bed and provides axial profiles for all temperatures, concentrations, and velocities. The model is computationally fast and flexible and allows for any system of adsorption and desorption reactions to be modeled, making the model applicable to any adsorption process. The model has been implemented in both gPROMS and Aspen Custom Modeler, and the behavior of the model has been verified.

  13. 1-D seismic velocity model and hypocenter relocation using double difference method around West Papua region

    West Papua region has fairly high of seismicity activities due to tectonic setting and many inland faults. In addition, the region has a unique and complex tectonic conditions and this situation lead to high potency of seismic hazard in the region. The precise earthquake hypocenter location is very important, which could provide high quality of earthquake parameter information and the subsurface structure in this region to the society. We conducted 1-D P-wave velocity using earthquake data catalog from BMKG for April, 2009 up to March, 2014 around West Papua region. The obtained 1-D seismic velocity then was used as input for improving hypocenter location using double-difference method. The relocated hypocenter location shows fairly clearly the pattern of intraslab earthquake beneath New Guinea Trench (NGT). The relocated hypocenters related to the inland fault are also observed more focus in location around the fault

  14. INFIL1D: a quasi-analytical model for simulating one-dimensional, constant flux infiltration

    The program INFIL1D is designed to calculate approximate wetting-front advance into an unsaturated, uniformly moist, homogeneous soil profile, under constant surface-flux conditions. The code is based on a quasi-analytical method, which utilizes an assumed invariant functional relationship between reduced (normalized) flux and water content. The code uses general hydraulic property data in tabular form to simulate constant surface-flux infiltration. 10 references, 4 figures

  15. Crustal seismic velocity in the Marche region (Central Italy): computation of a minimum 1-D model with seismic station corrections.

    Scarfì, L.; Istituto Nazionale di Geofisica e Vulcanologia, Sezione Catania, Catania, Italia; Imposa, S.; Dipartimento di Scienze Geologiche, University of Catania, Italy; Raffaele, R.; Dipartimento di Scienze Geologiche, Università di Catania, Italy; Scaltrito, A.; Istituto Nazionale di Geofisica e Vulcanologia, Sezione Catania, Catania, Italia

    2008-01-01

    A 1-D velocity model for the Marche region (central Italy) was computed by inverting P- and S-wave arrival times of local earthquakes. A total of 160 seismic events with a minimum of ten observations, a travel time residual ≤ 0.8 s and an azimuthal gap lower than 180° have been selected. This “minimum 1-D velocity model” is complemented by station corrections, which can be used to take into account possible near-surface velocity heterogeneities beneath each station. Using this new P-wave ...

  16. 1D spin-1/2 XY models as a testing ground for spin systems theory methods

    Elementary excitation energy spectrum that gives thermodynamic properties is calculated for few partial non-random and random versions of 1D spin-1/2 XY model. The exact result obtained is compared with the results derived within some well known approximate approaches that permits to understand the region of their validity. (author). 6 refs, 6 figs

  17. Mapping temporal extent of Chiang Mai floods using coupled 1-D and quasi 2-D floodplain inundation models

    Kowit Boonrawd; Chatchai Jothityangkoon

    2015-01-01

    A coupling of a 1-D flood routing model and quasi 2-D floodplain inundation model is applied for mapping spacetime flood extent. The routing model is formulated based on a non-linear storage-discharge relationship which is converted from an observed and synthetic rating curve. To draw the rating curve, required parameters for each reaches are estimated from hydraulic properties, floodplain geometry and vegetation and building cover of compound channels. The shape of the floodplain...

  18. Revisions to Photochemical Data for Use in Atmospheric Modeling

    Shutter, Joshua D.; Willacy, Karen; Allen, Mark

    2012-01-01

    Solar and stellar flux incident on an atmosphere can cause molecules to dissociate into highly reactive species and allows for photochemical processes to play a fundamental role in atmospheric chemistry. While models have tried to simulate such processes, they are extremely sensitive to photoabsorption cross-sections and quantum yields: two parameters that are important in determining the photodissociation rate, and hence the lifetime, of atmospheric compounds. Obtaining high-resolution and current data for these parameters is therefore highly desirable. Due to this, database and literature searches for high-quality cross-sections and quantum yields were performed and compiled for KINETICS, a Caltech/JPL Chemical Transport Model that can be used in modeling planetary atmospheres. Furthermore, photodissociation rates determined by running a Titan 1-D model were used to verify the completeness of these latest revisions.

  19. Modeling of atmospheric pollutant transfers

    Modeling is today a common tool for the evaluation of the environmental impact of atmospheric pollution events, for the design of air monitoring networks or for the calculation of pollutant concentrations in the ambient air. It is even necessary for the a priori evaluation of the consequences of a pollution plume. A large choice of atmospheric transfer codes exist but no ideal tool is available which allows to model all kinds of situations. The present day approach consists in combining different types of modeling according to the requested results and simulations. The CEA has a solid experience in this domain and has developed independent tools for the impact and safety studies relative to industrial facilities and to the management of crisis situations. (J.S.)

  20. An evaluation of 1D loss model collections for the off-design performance prediction of automotive turbocharger compressors

    Harley, P.; Spence, S.; Early, J.; Filsinger, D.; Dietrich, M.

    2013-12-01

    Single-zone modelling is used to assess different collections of impeller 1D loss models. Three collections of loss models have been identified in literature, and the background to each of these collections is discussed. Each collection is evaluated using three modern automotive turbocharger style centrifugal compressors; comparisons of performance for each of the collections are made. An empirical data set taken from standard hot gas stand tests for each turbocharger is used as a baseline for comparison. Compressor range is predicted in this study; impeller diffusion ratio is shown to be a useful method of predicting compressor surge in 1D, and choke is predicted using basic compressible flow theory. The compressor designer can use this as a guide to identify the most compatible collection of losses for turbocharger compressor design applications. The analysis indicates the most appropriate collection for the design of automotive turbocharger centrifugal compressors.

  1. An evaluation of 1D loss model collections for the off-design performance prediction of automotive turbocharger compressors

    Single-zone modelling is used to assess different collections of impeller 1D loss models. Three collections of loss models have been identified in literature, and the background to each of these collections is discussed. Each collection is evaluated using three modern automotive turbocharger style centrifugal compressors; comparisons of performance for each of the collections are made. An empirical data set taken from standard hot gas stand tests for each turbocharger is used as a baseline for comparison. Compressor range is predicted in this study; impeller diffusion ratio is shown to be a useful method of predicting compressor surge in 1D, and choke is predicted using basic compressible flow theory. The compressor designer can use this as a guide to identify the most compatible collection of losses for turbocharger compressor design applications. The analysis indicates the most appropriate collection for the design of automotive turbocharger centrifugal compressors

  2. Model for Simulation Atmospheric Turbulence

    Lundtang Petersen, Erik

    1976-01-01

    A method that produces realistic simulations of atmospheric turbulence is developed and analyzed. The procedure makes use of a generalized spectral analysis, often called a proper orthogonal decomposition or the Karhunen-Loève expansion. A set of criteria, emphasizing a realistic appearance, a...... eigenfunctions and estimates of the distributions of the corresponding expansion coefficients. The simulation method utilizes the eigenfunction expansion procedure to produce preliminary time histories of the three velocity components simultaneously. As a final step, a spectral shaping procedure is then applied....... The method is unique in modeling the three velocity components simultaneously, and it is found that important cross-statistical features are reasonably well-behaved. It is concluded that the model provides a practical, operational simulator of atmospheric turbulence....

  3. Improving practical atmospheric dispersion models

    The new generation of practical atmospheric dispersion model (for short range ≤ 30 km) are based on dispersion science and boundary layer meteorology which have widespread international acceptance. In addition, recent improvements in computer skills and the widespread availability of small powerful computers make it possible to have new regulatory models which are more complex than the previous generation which were based on charts and simple formulae. This paper describes the basis of these models and how they have developed. Such models are needed to satisfy the urgent public demand for sound, justifiable and consistent environmental decisions. For example, it is preferable that the same models are used to simulate dispersion in different industries; in many countries at present different models are used for emissions from nuclear and fossil fuel power stations. The models should not be so simple as to be suspect but neither should they be too complex for widespread use; for example, at public inquiries in Germany, where simple models are mandatory, it is becoming usual to cite the results from highly complex computational models because the simple models are not credible. This paper is written in a schematic style with an emphasis on tables and diagrams. (au) (22 refs.)

  4. Improved 1D model for calculating hydraulic properties in meandering rivers: Comparisons with measurements and 3D numerical simulations

    Haji Mohammadi, M.; Kang, S.; Sotiropoulos, F.

    2011-12-01

    It is well-known that meander bends impose local losses of energy to the flow in rivers. These local losses should be added together with friction loss to get the total loss of energy. In this work, we strive to develop a framework that considers the effect of bends in meandering rivers for one-dimensional (1-D) homogenous equations of flow. Our objective is to develop a simple, yet physically sound, and efficient model for carrying out engineering computations of flow through meander bends. We consider several approaches for calculating 1-D hydraulic properties of meandering rivers such as friction factor and Manning coefficient. The method of Kasper et al. (2005), which is based on channel top width, aspect ratio and radius of curvature, is adopted for further calculations. In this method, a correction is implemented in terms of local energy loss, due to helical motion and secondary currents of fluid particles driven by centrifugal force, in meanders. To validate the model, several test cases are simulated and the computed results are compared with the reported data in the literature in terms of water surface elevation, shear velocity, etc. For all cases the computed results are in reasonable agreement with the experimental data. 3-D RANS turbulent flow simulations are also carried out, using the method of Kang et al. (Adv. In Water Res., vol. 34, 2011), for different geometrical parameters of Kinoshita Rivers to determine the spatial distribution of shear stress on river bed and banks, which is the key factor in scour/deposition patterns. The 3-D solutions are then cross-sectionally averaged and compared with the respective solutions from the 1-D model. The comparisons show that the improved 1D model, which incorporates the effect of local bend loss, captures key flow parameters with reasonable accuracy. Our results also underscore the range of validity and limitations of 1D models for meander bend simulations. This work was supported by NSF Grants (as part of

  5. Local and regional minimum 1D models for earthquake location and data quality assessment in complex tectonic regions: application to Switzerland

    One-dimensional (1D) velocity models are still widely used for computing earthquake locations at seismological centers or in regions where three-dimensional (3D) velocity models are not available due to the lack of data of sufficiently high quality. The concept of the minimum 1D model with appropriate station corrections provides a framework to compute initial hypocenter locations and seismic velocities for local earthquake tomography. Since a minimum 1D model represents a solution to the coupled hypocenter-velocity problem it also represents a suitable velocity model for earthquake location and data quality assessment, such as evaluating the consistency in assigning pre-defined weighting classes and average picking error. Nevertheless, the use of a simple 1D velocity structure in combination with station delays raises the question of how appropriate the minimum 1D model concept is when applied to complex tectonic regions with significant three-dimensional (3D) variations in seismic velocities. In this study we compute one regional minimum 1D model and three local minimum 1D models for selected subregions of the Swiss Alpine region, which exhibits a strongly varying Moho topography. We compare the regional and local minimum 1D models in terms of earthquake locations and data quality assessment to measure their performance. Our results show that the local minimum 1D models provide more realistic hypocenter locations and better data fits than a single model for the Alpine region. We attribute this to the fact that in a local minimum 1D model local and regional effects of the velocity structure can be better separated. Consequently, in tectonically complex regions, minimum 1D models should be computed in sub-regions defined by similar structure, if they are used for earthquake location and data quality assessment. (authors)

  6. Evaluating 3-D and 1-D mathematical models for mass transport in heterogeneous biofilms

    Morgenroth, Eberhard Friedrich; Eberl, H.; van Loosdrecht, M. C. M.

    Results from a three dimensional model for heterogeneous biofilms including the numerical solution of hydrodynamics were compared to simplified one dimensional models. A one dimensional model with a variable diffusion coefficient over the thickness of the biofilm was well suited to approximate av...

  7. Evaluating 3-D and 1-D mathematical models for mass transport in heterogeneous biofilms

    Morgenroth, Eberhard Friedrich; Eberl, H.; van Loosdrecht, M. C. M.

    2000-01-01

    Results from a three dimensional model for heterogeneous biofilms including the numerical solution of hydrodynamics were compared to simplified one dimensional models. A one dimensional model with a variable diffusion coefficient over the thickness of the biofilm was well suited to approximate av...... in a growing biofilm and in a mushroom type biofilm assuming different modes of detachment....

  8. 1D model of precursors to frictional stick-slip motion allowing for robust comparison with experiments

    Amundsen, David Skålid; Thøgersen, Kjetil; Trømborg, Jørgen; Malthe-Sørenssen, Anders

    2011-01-01

    We study the dynamic behaviour of 1D spring-block models of friction when the external loading is applied from a side, and not on all blocks like in the classical Burridge-Knopoff-like models. Such a change in the loading yields specific difficulties, both from numerical and physical viewpoints. To address some of these difficulties and clarify the precise role of a series of model parameters, we start with the minimalistic model by Maegawa et al. (Tribol. Lett. 38, 313, 2010) which was proposed to reproduce their experiments about precursors to frictional sliding in the stick-slip regime. By successively adding (i) an internal viscosity, (ii) an interfacial stiffness and (iii) an initial tangential force distribution at the interface, we manage to (i) avoid the model's unphysical stress fluctuations, (ii) avoid its unphysical dependence on the spatial resolution and (iii) improve its agreement with the experimental results, respectively. Based on the behaviour of this improved 1D model, we develop an analyti...

  9. 2D MHD and 1D HD models of a solar flare -- a comprehensive comparison of the results

    Falewicz, R; Murawski, K; Srivastava, A K

    2015-01-01

    Without any doubt solar flaring loops possess a multi-thread internal structure that is poorly resolved and there are no means to observe heating episodes and thermodynamic evolution of the individual threads. These limitations cause fundamental problems in numerical modelling of flaring loops, such as selection of a structure and a number of threads, and an implementation of a proper model of the energy deposition process. A set of 1D hydrodynamic and 2D magnetohydrodynamic models of a flaring loop are developed to compare energy redistribution and plasma dynamics in the course of a prototypical solar flare. Basic parameters of the modeled loop are set according to the progenitor M1.8 flare recorded in the AR10126 on September 20, 2002 between 09:21 UT and 09:50 UT. The non-ideal 1D models include thermal conduction and radiative losses of the optically thin plasma as energy loss mechanisms, while the non-ideal 2D models take into account viscosity and thermal conduction as energy loss mechanisms only. The 2...

  10. 1D Nonisothermal Fiber Spinning Models for Thermotropic Polymeric Liquid Crystals

    Zhou, Hong; Forest, M. Gregory; Wang, Qi

    1997-11-01

    Previous slender one-dimensional models for axisymmetric filaments of liquid crystalline polymers (LCPs) are extended to include temperature-dependent material behavior and an energy equation. A two-phase model is posited, where below the glass transition temperature the material is modeled as a rigid cooling LCP fiber. We present families of numerical steady boundary-value solutions for thermal spinning flows; effects of temperature-dependent viscosity, LCP relaxation, excluded-volume potential, and viscous heating are modeled and exhibited. The predictions focus on thermal influence on spun fiber performance properties, such as birefringence and axial force, and process stability. A cooling ambient clearly contributes to faster stable spinning speeds.

  11. Stochastic models for atmospheric dispersion

    Ditlevsen, Ove Dalager

    2003-01-01

    Simple stochastic differential equation models have been applied by several researchers to describe the dispersion of tracer particles in the planetary atmospheric boundary layer and to form the basis for computer simulations of particle paths. To obtain the drift coefficient, empirical vertical...... variation by height is adopted. A particular problem for simulation studies with finite time steps is the construction of a reflection rule different from the rule of perfect reflection at the boundaries such that the rule complies with the imposed skewness of the velocity distribution for particle...... positions close to the boundaries. Different rules have been suggested in the literature with justifications based on simulation studies. Herein the relevant stochastic differential equation model is formulated in a particular way. The formulation is based on the marginal transformation of the position...

  12. Fuzzy knowledge-bBased curve evaluation for 1-D river model calibration

    Vidal, J.P.; Moisan, S.

    2006-01-01

    Model calibration requires an evaluation of the agreement between model outputs and reference data. This article presents an automatic fuzzy knowledge-based approach to identify the relevant discrepancies. This evaluation module is intended to be integrated within an existing knowledge-based calibration support system.

  13. Introduction to Displacements, Strains and Stresses in a 1D CVM-model

    Frandsen, Jens Ole

    This lecture note contains an introduction to displacements, strains and stresses in an one-dimensional sg-FVM model of a tensile test bar.......This lecture note contains an introduction to displacements, strains and stresses in an one-dimensional sg-FVM model of a tensile test bar....

  14. Assessing the impact of different sources of topographic data on 1-D hydraulic modelling of floods

    A. Md Ali; D. P. Solomatine; G. Di Baldassarre

    2014-01-01

    Topographic data, such as digital elevation models (DEMs), are essential input in flood inundation modelling. DEMs can be derived from several sources either through remote sensing techniques (space-borne or air-borne imagery) or from traditional methods (ground survey). The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), the Shuttle Radar Topography Mission (SRTM), the Light Detection and Ranging (LiDAR), and topographic contour maps are some of the ...

  15. Assessing the impact of different sources of topographic data on 1-D hydraulic modelling of floods

    A. Md Ali; D. P. Solomatine; G. Di Baldassarre

    2015-01-01

    Topographic data, such as digital elevation models (DEMs), are essential input in flood inundation modelling. DEMs can be derived from several sources either through remote sensing techniques (spaceborne or airborne imagery) or from traditional methods (ground survey). The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), the Shuttle Radar Topography Mission (SRTM), the light detection and ranging (lidar), and topographic contour maps are some of the mo...

  16. Bottom Roughness and Bathymetry Estimation of 1-D Shallow Water Equations Model Using Ensemble Kalman Filter

    Hooshyar, M.; Hagen, S. C.; Wang, D.

    2014-12-01

    Hydrodynamic models are widely applied to coastal areas in order to predict water levels and flood inundation and typically involve solving a form of the Shallow Water Equations (SWE). The SWE are routinely discretized by applying numerical methods, such as the finite element method. Like other numerical models, hydrodynamic models include uncertainty. Uncertainties are generated due to errors in the discrete approximation of coastal geometry, bathymetry, bottom friction and forcing functions such as tides and wind fields. Methods to counteract these uncertainties should always begin with improvements to physical characterization of: the geometric description through increased resolution, parameters that describe land cover variations in the natural and urban environment, parameters that enhance transfer of surface forcings to the water surface, open boundary forcings, and the wetting/drying brought upon by flood and ebb cycles. When the best possible physical representation is achieved, we are left with calibration and data assimilation to reduce model uncertainty. Data assimilation has been applied to coastal hydrodynamic models to better estimate system states and/or system parameters by incorporating observed data into the model. Kalman Filter is one of the most studied data assimilation methods that minimizes the mean square errors between model state estimations and the observed data in linear systems (Kalman , 1960). For nonlinear systems, as with hydrodynamic models, a variation of Kalman filter called Ensemble Kalman Filter (EnKF), is applied to update the system state according to error statistics in the context of Monte Carlo simulations (Evensen , 2003) & (Hitoshi et. al, 2014). In this research, Kalman Filter is incorporated to simultaneously estimate an influential parameter used in the shallow water equations, bottom roughness, and to adjust the physical feature of bathymetry. Starting from an initial estimate of bottom roughness and bathymetry, and

  17. Study of the tritium production in a 1-D blanket model with Monte Carlo methods

    Cubí Ricart, Álvaro

    2015-01-01

    In this work a method to collapse a 3D geometry into a mono dimensional model of a fusion reactor blanket is developed and tested. Using this model, neutron and photon uxes and its energy deposition will be obtained with a Monte Carlo code. This results will allow to calculate the TBR and the thermal power of the blanket and will be able to be integrated in the AINA code.

  18. CATHARE Multi-1D Modeling of Coolant Mixing in VVER-1000 for RIA Analysis

    Sabotinov, L.; Kolev, N. P.; Spasov, I.; Donov, J.

    2010-01-01

    The paper presents validation results for multichannel vessel thermal-hydraulic models in CATHARE used in coupled 3D neutronic/thermal hydraulic calculations. The mixing is modeled with cross flows governed by local pressure drops. The test cases are from the OECD VVER-1000 coolant transient benchmark (V1000CT) and include asymmetric vessel flow transients and main steam line break (MSLB) transients. Plant data from flow mixing experiments are available for comparison. Sufficient mesh refinem...

  19. Test of 1-D transport models, and their predictions for ITER

    A number of proposed tokamak thermal transport models are tested by comparing their predictions with measurements from several tokamaks. The necessary data have been provided for a total of 75 discharges from C-mod, DIII-D, JET, JT-60U, T10, and TFTR. A standard prediction methodology has been developed, and three codes have been benchmarked; these 'standard' codes have been relied on for testing most of the transport models. While a wide range of physical transport processes has been tested, no single model has emerged as clearly superior to all competitors for simulating H-mode discharges. In order to winnow the field, further tests of the effect of sheared flows and of the 'stiffness' of transport are planned. Several of the models have been used to predict ITER performance, with widely varying results. With some transport models ITER's predicted fusion power depends strongly on the 'pedestal' temperature, but ∼ 1GW (Q=10) is predicted for most models if the pedestal temperature is at least 4 keV. (author)

  20. Assessment of TRAC-BF1 1D reflood model with CCTF and SCTF data

    Post test calculations for six selected Cylindrical Core Test Facility (CCTF) and Slab Core Test Facility (SCTF) tests were performed to assess the core thermal hydraulic models of the TRAC-BF1 code during the reflood in a PWR LOCA. A special version of the TRAC code was developed at JAERI by implementing the constitutive package of the TRAC-BF1 code into the TRAC-PF1 code for this assessment. The TRAC-BF1 model predicted well the void fraction at either bottom or top part of the core and overpredicted the void fraction at the center part of the core in the CCTF and SCTF tests performed under so-called licensing conditions. The TRAC-BF1 model overpredicted the clad temperatures at the center part of the core. The TRAC-BF1 model predicted a jump of void fraction where the flow pattern transition between the bubbly/slug flow and the annular/dispersed flow regimes occurred. The jump caused the water mass flow rate to be unstable and resulted in the overprediction of the void fraction at the center part of the core. It was also found that the TRAC-BF1 film boiling model underestimated the heat transfer coefficient in the vicinity of the quench front and caused the quench front propagation to be delayed. These assessment results suggest the following areas should be improved in future to apply the TRAC-BF1 code to the reflood in a PWR LOCA: (1) Core hydraulic model where flow pattern transition occurs, (2) Core heat transfer model in the film boiling regime, especially for the dependence on the distance from the quench front. (author)

  1. DCC Dynamics in (2+1)D-O(3) model

    Holzwarth, G

    2001-01-01

    The dynamics of symmetry-breaking after a quench is numerically simulated on a lattice for the (2+1)-dimensional O(3) model. In addition to the standard sigma-model with temperature-dependent Phi^4-potential the energy functional includes a four-derivative current-current coupling to stabilize the size of the emerging extended topological textures. The total winding number can be conserved by constraint. As a model for the chiral phase transition during the cooling phase after a hadronic collision this allows to investigate the interference of 'baryon-antibaryon' production with the developing disoriented aligned domains. The growth of angular correlations, condensate, average orientation is studied in dependence of texture size, quench rate, symmetry breaking. The classical dissipative dynamics determines the rate of energy emitted from the relaxing source for each component of the 3-vector field which provides a possible signature for domains of Disoriented Chiral Condensate. We find that the 'pions' are em...

  2. Biot-JKD model: simulation of 1D transient poroelastic waves with fractional derivatives

    Blanc, Emilie; Lombard, Bruno

    2012-01-01

    A time-domain numerical modeling of Biot poroelastic waves is presented. The viscous dissipation occurring in the pores is described using the dynamic permeability model developed by Johnson-Koplik-Dashen (JKD). Some of the coefficients in the Biot-JKD model are proportional to the square root of the frequency: in the time-domain, these coefficients introduce order 1/2 shifted fractional derivatives involving a convolution product. Based on a diffusive representation, the convolution kernel is replaced by a finite number of memory variables that satisfy local-in-time ordinary differential equations. Thanks to the dispersion relation, the coefficients in the diffusive representation are obtained by performing an optimization procedure in the frequency range of interest. A splitting strategy is then applied numerically: the propagative part of Biot-JKD equations is discretized using a fourth-order ADER scheme on a Cartesian grid, whereas the diffusive part is solved exactly. Comparisons with analytical solution...

  3. Development of a 1D-2D coupled hydrodynamic model for the Øyeren Delta in southern Norway

    2011-01-01

    In this study a coupled 1D-2D hydrodynamic model, MIKE FLOOD was used to simulate flood inundation extent, water levels and water velocities in the delta region of Lake Øyeren in southern Norway. The objective was to evaluate the improvement gained using a more complex framework. In addition, the credibility of existing flood zone maps made for Lillestrøm by Norges Vassdrag- og Energidirektorat (NVE) in 2005 was assessed. They were based on the assumption that the water levels predicted for F...

  4. The band-centre anomaly in the 1D Anderson model with correlated disorder

    We study the band-centre anomaly in the one-dimensional Anderson model with weak correlated disorder. Our analysis is based on the Hamiltonian map approach; the correspondence between the discrete model and its continuous counterpart is discussed in detail. We obtain analytical expressions of the localization length and of the invariant measure of the phase variable, valid for energies in a neighbourhood of the band centre. By applying these general results to specific forms of correlated disorder, we show how correlations can enhance or suppress the anomaly at the band centre. (paper)

  5. Modeling 1D structures on semiconductor surfaces: synergy of theory and experiment

    Atomic scale nanowires attract enormous interest in a wide range of fields. On the one hand, due to their quasi-one-dimensional nature, they can act as an experimental testbed for exotic physics: Peierls instability, charge density waves, and Luttinger liquid behavior. On the other hand, due to their small size, they are of interest not only for future device applications in the micro-electronics industry, but also for applications regarding molecular electronics. This versatile nature makes them interesting systems to produce and study, but their size and growth conditions push both experimental production and theoretical modeling to their limits. In this review, modeling of atomic scale nanowires on semiconductor surfaces is discussed, focusing on the interplay between theory and experiment. The current state of modeling efforts on Pt- and Au-induced nanowires on Ge(001) is presented, indicating their similarities and differences. Recently discovered nanowire systems (Ir, Co, Sr) on the Ge(001) surface are also touched upon. The importance of scanning tunneling microscopy as a tool for direct comparison of theoretical and experimental data is shown, as is the power of density functional theory as an atomistic simulation approach. It becomes clear that complementary strengths of theoretical and experimental investigations are required for successful modeling of the atomistic nanowires, due to their complexity. (topical review)

  6. Pushing 1D CCSNe to explosions: model and SN 1987A

    Perego, A; Fröhlich, C; Ebinger, K; Eichler, M; Casanova, J; Liebendoerfer, M; Thielemann, F -K

    2015-01-01

    We report on a method, PUSH, for triggering core-collapse supernova explosions of massive stars in spherical symmetry. We explore basic explosion properties and calibrate PUSH such that the observables of SN1987A are reproduced. Our simulations are based on the general relativistic hydrodynamics code AGILE combined with the detailed neutrino transport scheme IDSA for electron neutrinos and ALS for the muon and tau neutrinos. To trigger explosions in the otherwise non-exploding simulations, we rely on the neutrino-driven mechanism. The PUSH method locally increases the energy deposition in the gain region through energy deposition by the heavy neutrino flavors. Our setup allows us to model the explosion for several seconds after core bounce. We explore the progenitor range 18-21M$_{\\odot}$. Our studies reveal a distinction between high compactness (HC) and low compactness (LC) progenitor models, where LC models tend to explore earlier, with a lower explosion energy, and with a lower remnant mass. HC models are...

  7. Modeling Physiological Systems in the Human Body as Networks of Quasi-1D Fluid Flows

    Staples, Anne

    2008-11-01

    Extensive research has been done on modeling human physiology. Most of this work has been aimed at developing detailed, three-dimensional models of specific components of physiological systems, such as a cell, a vein, a molecule, or a heart valve. While efforts such as these are invaluable to our understanding of human biology, if we were to construct a global model of human physiology with this level of detail, computing even a nanosecond in this computational being's life would certainly be prohibitively expensive. With this in mind, we derive the Pulsed Flow Equations, a set of coupled one-dimensional partial differential equations, specifically designed to capture two-dimensional viscous, transport, and other effects, and aimed at providing accurate and fast-to-compute global models for physiological systems represented as networks of quasi one-dimensional fluid flows. Our goal is to be able to perform faster-than-real time simulations of global processes in the human body on desktop computers.

  8. ANISN-E, 1-D Transport Program ANISN with Exponential Model. ANISN-JR, 1-D Transport Program ANISN with ZZ JSD Data and Flux Plot

    . B - Method of solution: ANISN solves the one-dimensional Boltzmann transport equation for neutrons or gamma-rays in slab, sphere, or cylinder geometry. The source may be fixed, fission or a subcritical combination of the two. Criticality search may be performed on any one of several parameters. Cross sections may be weighted using the space and energy dependent flux generated in solving the transport equation. ANISN-E : Besides diamond and weighted difference supplementary equations, exponential supplementary equations are available. The new model: (1) always gives positive solutions, without using any 'fix up' technique provided that the source is non-negative; (2) allows, especially in deep penetration problems, the use of larger spatial meshes, hence requires shorter computer times than the ones requested by the diamond model combined with various types of fix up techniques or by weighted difference schemes to get the same accuracy; (3) supplies solutions that are always reasonable overestimates of the exact solution. In ANISN-JR, some optional functions are added to increase the utility of the code: (1) print the total fluxes at boundary points of all mesh intervals. (The original ANISN prints the total fluxes at midpoint only.) (2) calculate, print and plot the lethargy width spectra. (3) print the angular fluxes at only required mesh boundaries or midpoints (maximum 10 points). The original ANISN prints at mid- point of all meshes, and therefore the number of print pages becomes vast according to the number of spatial and angular meshes. (4) use the asymmetric quadrature set. (5) calculate and plot the reaction rates for neutron and gamma-ray detectors, and collapse the response functions of detectors. (6) generate volume-flux weighted cross sections for arbitrary zone or region. In the original ANISN, the cross sections can be collapsed only for a homogeneous zone or region. (7) collapse into few group cross sections in ANISN, DOT, or TWOTRAN format. (In

  9. Prediction of the expansion velocity of ultracold 1D quantum gases for integrable models

    Mei, Zhongtao; Vidmar, Lev; Heidrich-Meisner, Fabian; Bolech, Carlos

    In the theory of Bethe-ansatz integrable quantum systems, rapidities play an important role as they are used to specify many-body states. The physical interpretation of rapidities going back to Sutherland is that they are the asymptotic momenta after letting a quantum gas expand into a larger volume rendering it dilute and noninteracting. We exploit this picture to calculate the expansion velocity of a one-dimensional Fermi-Hubbard model by using the distribution of rapidities defined by the initial state. Our results are consistent with the ones from time-dependent density-matrix renormalization. We show in addition that an approximate Bethe-ansatz solution works well also for the Bose-Hubbard model. Our results are of interests for future sudden-expansion experiments with ultracold quantum gases.

  10. Medicanes in an ocean–atmosphere coupled regional climate model

    N. Akhtar

    2014-03-01

    Full Text Available So-called medicanes (Mediterranean hurricanes are meso-scale, marine, and warm-core Mediterranean cyclones that exhibit some similarities to tropical cyclones. The strong cyclonic winds associated with medicanes threaten the highly populated coastal areas around the Mediterranean basin. To reduce the risk of casualties and overall negative impacts, it is important to improve the understanding of medicanes with the use of numerical models. In this study, we employ an atmospheric limited-area model (COSMO-CLM coupled with a one-dimensional ocean model (1-D NEMO-MED12 to simulate medicanes. The aim of this study is to assess the robustness of the coupled model in simulating these extreme events. For this purpose, 11 historical medicane events are simulated using the atmosphere-only model, COSMO-CLM, and coupled model, with different setups (horizontal atmospheric grid-spacings of 0.44°, 0.22°, and 0.08°; with/without spectral nudging, and an ocean grid-spacing of 1/12°. The results show that at high-resolution, the coupled model is able to not only simulate most of medicane events but also improve the track length, core temperature, and wind speed of simulated medicanes compared to the atmosphere-only simulations. The results suggest that the coupled model is more proficient for systemic and detailed studies of historical medicane events, and that this model can be an effective tool for future projections.

  11. A 1-D Size Specific Numerical Model for Gravel Transport That Includes Sediment Exchange with a Floodplain

    Lauer, Wesley; Viparelli, Enrica; Piegay, Herve

    2014-05-01

    Sedimentary deposits adjacent to rivers can represent important sources and sinks for bed material sediment, particularly on decadal and longer timescales. The Morphodynamics and Sediment Tracers in 1-D model (MAST-1D) is a size-specific sediment transport model that allows for active exchange between channel and floodplain sediment on river reaches of tens to hundreds of kilometers in length. The model is intended to provide a mechanism for performing a first-order assessment of the likely importance of off-channel sediment exchange in controlling decadal-scale geomorphic trends, thereby helping plan and/or prioritize field data collection and higher resolution modeling work. The model develops a sediment budget for short segments of an alluvial valley. Each segment encompasses several active river bends. In each segment, a sediment transport capacity computation is performed to determine the downstream flux of bed material sediment, following the approach of most other 1-D sediment transport models. However, the model differs from most other bed evolution models in that sediment can be exchanged with the floodplain in each segment, and mass conservation is applied to both the active layer and floodplain sediment storage reservoirs. The potential for net imbalances in overall exchange as well as the size specific nature of the computations allows the model to simulate reach-scale aggradation/degradation and/or changes in bed texture. The inclusion of fine sediment in the model allows it to track geochemical tracer material and also provides a mechanism to simulate, to first order, the effects of changes in the supply of silt and clay on overall channel hydraulic capacity. The model is applied to a ~40 km reach of the Ain River, a tributary of the Rhône River in eastern France that has experienced a significant sediment deficit as a result of the construction of several dams between 1920 and 1970. MAST-1D simulations result in both incision and the formation of a

  12. Constraining the temporal evolution of a deep hypersaline anoxic basin by 1D geochemical modelling

    Goldhammer, Tobias; Aiello, Ivano; Zabel, Matthias

    2014-05-01

    Deep hypersaline anoxic basins (DHABs) are seafloor features of the accretionary prism of the Mediterranean Ridge. They have formed by the dissolution of exhumed shallow Messinian evaporites and subsequent concentration of the ultra-saline solutions in depressions on the seafloor. As an example, the horseshoe-shaped Urania basin is a DHAB south of the Peloponnese peninsula contains one of the most saline (about six times higher than Mediterranean seawater) and sulfidic (up to 15mM) water bodies of the Earth. Furthermore, its deepest part is underlain by a mud volcano that is responsible for the injection of fluid mud beneath the brine lake, with exceptionally sharp chemoclines between water column, brine, and mud layer. We here present a model approach to reconstruct the temporal aspects of the formation, dynamics and persistence of the brine-mud-system in the deep pit of the Urania Basin. Based on data from a sampling campaign with RV Meteor (Cruise M84/1 in February 2011), we set up a one-dimensional geochemical model that integrates diffusion, reaction and advective transport and mixing. Using a set of model preconditions, we aimed to answer (1) which processes are required to maintain the current situation of steep chemical gradients of the brine-mud-system, (2) how fast the current situation could have developed under different scenarios, and (3) how long such extraordinary conditions could have persisted through Earth's history. We further discuss the consequences of the temporal framework for the evolution of prokaryotic life in this extreme habitat.

  13. Transverse susceptibility of the 1D isotropic XY-model at zero temperature

    An exact expression is obtained for the dynamic transverse susceptibility Xxx (w,i,j) of the one-dimensional isotropic XY-model both on an open and on a closed chain with arbitrary number of sites at zero temperature, when the transverse field is greater than the absolute value of the exchange constant. The dynamic transverse wave-vector-dependent susceptibility for the closed infinite chain is also determined, and it is shown that in all cases the isothermal susceptibility is identical to the static susceptibility. (Author)

  14. 1-D transient numerical model of a regenerator in a novel sub Kelvin Active Magnetic Regenerative Refrigerator

    Jahromi, Amir E.; Miller, Franklin K.

    2016-03-01

    A sub Kelvin Active Magnetic Regenerative Refrigerator (AMRR) is being developed at the University of Wisconsin - Madison. This AMRR consists of two circulators, two regenerators, one superleak, one cold heat exchanger, and two warm heat exchangers. The circulators are novel non-moving part pumps that reciprocate a superfluid mixture of 4He-3He in the system. Heat from the mixture is removed within the two regenerators of this tandem system. An accurate model of the regenerators in this AMRR is necessary in order to predict the performance of these components, which in turn helps predicting the overall performance of the AMRR system. This work presents modeling methodology along with results from a 1-D transient numerical model of the regenerators of an AMRR capable of removing 2.5 mW at 850 mK at cyclic steady state.

  15. 1D and 2D Numerical Modeling for Solving Dam-Break Flow Problems Using Finite Volume Method

    Szu-Hsien Peng

    2012-01-01

    Full Text Available The purpose of this study is to model the flow movement in an idealized dam-break configuration. One-dimensional and two-dimensional motion of a shallow flow over a rigid inclined bed is considered. The resulting shallow water equations are solved by finite volumes using the Roe and HLL schemes. At first, the one-dimensional model is considered in the development process. With conservative finite volume method, splitting is applied to manage the combination of hyperbolic term and source term of the shallow water equation and then to promote 1D to 2D. The simulations are validated by the comparison with flume experiments. Unsteady dam-break flow movement is found to be reasonably well captured by the model. The proposed concept could be further developed to the numerical calculation of non-Newtonian fluid or multilayers fluid flow.

  16. Plasma Processes : A self-consistent kinetic modeling of a 1-D, bounded, plasma in equilibrium

    Monojoy Goswami; H Ramachandran

    2000-11-01

    A self-consistent kinetic treatment is presented here, where the Boltzmann equation is solved for a particle conserving Krook collision operator. The resulting equations have been implemented numerically. The treatment solves for the entire quasineutral column, making no assumptions about mfp/, where mfp is the ion-neutral collision mean free path and the size of the device. Coulomb collisions are neglected in favour of collisions with neutrals, and the particle source is modeled as a uniform Maxwellian. Electrons are treated as an inertialess but collisional fluid. The ion distribution function for the trapped and the transiting orbits is obtained. Interesting findings include the anomalous heating of ions as they approach the presheath, the development of strongly non-Maxwellian features near the last mfp, and strong modifications of the sheath criterion.

  17. On non-minimal N = 4 supermultiplets in 1D and their associated {sigma}-models

    Gonzales, Marcelo; Khodaee, Sadi; Toppan, Francesco, E-mail: marcbino@cbpf.b, E-mail: khodaee@cbpf.b, E-mail: toppan@cbpf.b

    2010-10-15

    We construct the non-minimal linear representations of the N = 4 Extended Supersymmetry in one-dimension. They act on 8 bosonic and 8 fermionic fields. Inequivalent representations are specified by the mass-dimension of the fields and the connectivity of the associated graphs. The oxidation to minimal N = 5 linear representations is given. Two types of N = 4 {sigma}-models based on non-minimal representations are obtained: the resulting off-shell actions are either manifestly invariant or depend on a constrained pre potential. The connectivity properties of the graphs play a decisive role in discriminating inequivalent actions. These results find application in partial breaking of supersymmetric theories. (author)

  18. Exact 1-D Model for Coherent Synchrotron Radiation with Shielding and Bunch Compression

    Mayes, Christopher

    2008-01-01

    Coherent Synchrotron Radiation has been studied effectively using a 1-dimensional model for the charge distribution in the realm of small angle approximations and high energies. Here we use Jefimenko's form of Maxwell's equations, without such approximations, to calculate the exact wake-fields due to this effect in multiple bends and drifts. It has been shown before that the influence of a drift can propagate well into a subsequent bend. We show, for reasonable parameters, that the influence of a previous bend can also propagate well into a subsequent bend, and that this is especially important at the beginning of a bend. Shielding by conducting parallel plates is simulated using the image charge method. We extend the formalism to situations with compressing and decompressing distributions, and conclude that simpler approximations to bunch compression usually overestimates the effect. Additionally, an exact formula for the coherent power radiated by a Gaussian bunch is derived by considering the coherent sync...

  19. Exact solution of the 1D Hubbard model with NN and NNN interactions in the narrow-band limit

    Mancini, Ferdinando; Plekhanov, Evgeny; Sica, Gerardo

    2013-10-01

    We present the exact solution, obtained by means of the Transfer Matrix (TM) method, of the 1D Hubbard model with nearest-neighbor (NN) and next-nearest-neighbor (NNN) Coulomb interactions in the atomic limit ( t = 0). The competition among the interactions ( U, V 1, and V 2) generates a plethora of T = 0 phases in the whole range of fillings. U, V 1, and V 2 are the intensities of the local, NN and NNN interactions, respectively. We report the T = 0 phase diagram, in which the phases are classified according to the behavior of the principal correlation functions, and reconstruct a representative electronic configuration for each phase. In order to do that, we make an analytic limit T → 0 in the transfer matrix, which allows us to obtain analytic expressions for the ground state energies even for extended transfer matrices. Such an extension of the standard TM technique can be easily applied to a wide class of 1D models with the interaction range beyond NN distance, allowing for a complete determination of the T = 0 phase diagrams.

  20. An analytical 1-D model for vertical momentum and energy flux through a fully developed wind farm

    Markfort, Corey D.; Zhang, Wei; Porté-Agel, Fernando

    2014-05-01

    Wind farms capture momentum from the atmospheric boundary layer (ABL) both at the leading edge and from the atmosphere above. Momentum is advected into the wind farm and wake turbulence draws excess momentum in from between turbines until momentum is only available from above the wind farm. This distance can be described by the so-called drag development length scale, which arises from the canopy drag force term in the momentum equation. At this point the flow can be considered fully developed. The horizontally-averaged velocity profile for a fully developed wind farm flow exhibits a characteristic inflection point near the top of the wind farm, similar to that of sparse canopy-type flows (Markfort et al., JoT, 2012). The inflected vertical velocity profile is associated with the presence of a dominant characteristic turbulence scale, which may be responsible for a significant portion of the vertical momentum flux. We evaluate an analytical canopy-type flow model for wind farm-atmosphere interaction by testing it against wind-tunnel experimental data of flow through a model wind farm. The model is adapted to predict the mean flow, vertical momentum flux, and the mean kinetic energy flux as well as kinetic energy dissipation within the wind farm. This model is particularly useful for wind farm configuration optimization, considering wind turbine spacing and surface roughness and can also be useful to represent wind farms in regional scale atmospheric simulations.

  1. A 1D pulse wave propagation model of the hemodynamics of calf muscle pump function.

    Keijsers, J M T; Leguy, C A D; Huberts, W; Narracott, A J; Rittweger, J; van de Vosse, F N

    2015-07-01

    The calf muscle pump is a mechanism which increases venous return and thereby compensates for the fluid shift towards the lower body during standing. During a muscle contraction, the embedded deep veins collapse and venous return increases. In the subsequent relaxation phase, muscle perfusion increases due to increased perfusion pressure, as the proximal venous valves temporarily reduce the distal venous pressure (shielding). The superficial and deep veins are connected via perforators, which contain valves allowing flow in the superficial-to-deep direction. The aim of this study is to investigate and quantify the physiological mechanisms of the calf muscle pump, including the effect of venous valves, hydrostatic pressure, and the superficial venous system. Using a one-dimensional pulse wave propagation model, a muscle contraction is simulated by increasing the extravascular pressure in the deep venous segments. The hemodynamics are studied in three different configurations: a single artery-vein configuration with and without valves and a more detailed configuration including a superficial vein. Proximal venous valves increase effective venous return by 53% by preventing reflux. Furthermore, the proximal valves shielding function increases perfusion following contraction. Finally, the superficial system aids in maintaining the perfusion during the contraction phase and reduces the refilling time by 37%. PMID:25766693

  2. Phase transitions at strong coupling in the 2+1-d abelian Higgs model

    We study, using numerical Monte-Carlo simulations, an effective description of the 2+1 dimensional Abelian Higgs model which is valid at strong coupling, in the broken symmetry sector. In this limit, the massive gauge boson and the massive neutral Higgs decouple leaving only the massive vortices. The vortices have no long range interactions. We find a phase transition as the mass of the vortices is made lighter and lighter. At the transition, the contributions to the functional integral come from a so-called infinite vortex anti-vortex loop. Adding the Chern-Simons term simply counts the linking number between the vortices. We find that the Wilson loop exhibits perimeter law behaviour in both phases, although the polarization cloud increases by an order of magnitude at the transition. We also study the 't Hooft loop. We find the 't Hooft loop exhibits perimeter law behaviour in the presence of the Chern-Simons term but is trivial in its absence. Thus we have a theory with perimeter law for both the Wilson loop and the 't Hooft loop, but contains no massless particles

  3. Assessment of phenol infiltration resilience in soil media by HYDRUS-1D transport model for a waste discharge site.

    Adhikari, K; Pal, S; Chakraborty, B; Mukherjee, S N; Gangopadhyay, A

    2014-10-01

    The movement of contaminants through soil imparts a variety of geo-environmental problem inclusive of lithospheric pollution. Near-surface aquifers are often vulnerable to contamination from surface source if overlying soil possesses poor resilience or contaminant attenuation capacity. The prediction of contaminant transport through soil is urged to protect groundwater from sources of pollutants. Using field simulation through column experiments and mathematical modeling like HYDRUS-1D, assessment of soil resilience and movement of contaminants through the subsurface to reach aquifers can be predicted. An outfall site of effluents of a coke oven plant comprising of alarming concentration of phenol (4-12.2 mg/L) have been considered for studying groundwater condition and quality, in situ soil characterization, and effluent characterization. Hydrogeological feature suggests the presence of near-surface aquifers at the effluent discharge site. Analysis of groundwater of nearby locality reveals the phenol concentration (0.11-0.75 mg/L) exceeded the prescribed limit of WHO specification (0.002 mg/L). The in situ soil, used in column experiment, possess higher saturated hydraulic conductivity (KS  = 5.25 × 10(-4) cm/s). The soil containing 47 % silt, 11 % clay, and 1.54% organic carbon content was found to be a poor absorber of phenol (24 mg/kg). The linear phenol adsorption isotherm model showed the best fit (R(2) = 0.977, RMSE = 1.057) to the test results. Column experiments revealed that the phenol removal percent and the length of the mass transfer zone increased with increasing bed heights. The overall phenol adsorption efficiency was found to be 42-49%. Breakthrough curves (BTCs) predicted by HYDRUS-1D model appears to be close fitting with the BTCs derived from the column experiments. The phenol BTC predicted by the HYDRUS-1D model for 1.2 m depth subsurface soil, i.e., up to the depth of groundwater in the study area, showed that the exhaustion

  4. Development of 1D and 2D coupled model to simulate urban inundation: An application to Beijing Olympic Village

    LI WeiFeng; CHEN QiuWen; MAO JingQiao

    2009-01-01

    Urban inundation due to anomalous storms is a serious problem for many cities worldwide. Therefore, it is important to accurately simulate urban hydrological processes and efficiently predict the potential risks of urban floods for the improvement of drainage designs and implementation of emergency ac-tions. However, the complexity of urban landforma and the diversity of hydraulic infrastructure pose particular challenges for the simulation and risk assessment of urban drainage processes. This study developed a methodology to comprehensively simulate inundation processes by dynamically coupling 1D and 2D hydrodynamic models. By allowing the simultaneous solution of the processes of rainfall and runoff, urban drainage, and flooding, this method can be used to estimate the potential inundation risks of any designed drainage system. Furthermore, a Geographical Information System (GIS) based platform was fully integrated with the model engine to effectively illustrate the context of the problem. The developed model was then demonstrated on the Beijing 2008 Olympic Village under the conditions of the 5-year and 50-year design storms. The sewer discharge, channel discharge, and flood propaga-tion (inundation initiation, extent, depths, and duration) were numerically validated and analyzed. The results identified the potential inundation risks. From the study, it is found that the coupled GIS and 1D and 2D hydrodynamic models have the potential to simulate urban inundation processes, and hence efficiently predict flood risks and support cost-effective drainage design and management. It also im-plies promising prospects about the wide availability of high quality digital data, GIS techniques, and well-developed monitoring infrastructure to develop online urban inundation forecasts.

  5. Models of Pluto's upper atmosphere

    Best guesses as to the thermal structure and composition of Pluto's atmosphere have led to speculations of substantial loss rates (∼1028 s-1) of methane from the planet over cosmogonic time scales. Results from recent stellar occultation measurements, and using a Parker-type hydrodynamic calculation, show that the loss rates may actually be lower by as much as a factor ∼5, depending upon the efficiency of heating of the atmosphere via the absorption of solar EUV and upon the true atmospheric composition, if the thermal structure of the upper atmosphere is properly taken into account. The loss rate may even be less (by another factor ∼10) if there is minimal heating of the upper atmosphere

  6. Analytical 1-D dual-porosity equivalent solutions to 3-D discrete single-continuum models. Application to karstic spring hydrograph modelling

    Cornaton, F

    2011-01-01

    One dimensional analytical porosity-weighted solutions of the dual-porosity model are derived, providing insights on how to relate exchange and storage coefficients to the volumetric density of the high-permeability medium. It is shown that porosity-weighted storage and exchange coefficients are needed when handling highly heterogeneous systems - such as karstic aquifers - using equivalent dual-porosity models. The sensitivity of these coefficients is illustrated by means of numerical experiments with theoretical karst systems. The presented 1-D dual-porosity analytical model is used to reproduce the hydraulic responses of reference 3-D karst aquifers, modelled by a discrete single-continuum approach. Under various stress conditions, simulation results show the relations between the dual-porosity model coefficients and the structural features of the discrete single-continuum model. The calibration of the equivalent 1-D analytical dual-porosity model on reference hydraulic responses confirms the dependence of ...

  7. Combinatorial variation of structure in considerations of compound lumping in one- and two-dimensional property representations of condensable atmospheric organic compounds. 1. Lumping by 1-D volatility with nC fixed

    Pankow, James F.; Niakan, Negar; Asher, William E.

    2013-12-01

    Many current models that aim to predict urban and regional levels of organic particulate matter (OPM) use either the 2 product (2p) framework for secondary organic aerosol (SOA) formation, or a static 1-D volatility basis set (1-D-VBS). These approaches assume that: 1) the compounds involved in OPM condensation/evaporation can be lumped simply by volatility with no specificity regarding carbon number nC, MW, or polar functionality; 2) water uptake does not occur; and 3) the compounds are non-ionizing. This work considers the consequences for uniphasic PM caused by the first two assumptions due to effects of the condensed-phase mean molecular weight MWbar and activity coefficients (ζi), including when RH (relative humidity) > 0. Setting nC = 10 for all bins, multiple chemical structures were developed for each bin of a 1-D-VBS for un-aged SOA in the α-pinene/ozone system. For each bin, a group-contribution vapor pressure (pLo) prediction method was used to find multiple structures such that the groups-based log pLo for nC = 10 and variable numbers of aldehyde, ketone, hydroxyl, and carboxylic acid groups agrees, within ±0.5, with the bin volatility. The number of possible combinations with one structure taken from each bin was 17,640. The Raster-Roulette Organic Aerosol (RROA) model was used to calculate the equilibrium mass concentrations (μg m-3) of OPM (Mo) and co-condensed water (Mw) at 25 °C for each combination for ranges of RH and ΔHC (change in parent hydrocarbon concentration). UNIFAC was used to determine the needed values of ζi. Frequency distributions from RROA for Mo, Mw, and the O:C ratio were developed. For Mo levels typical of the ambient atmosphere, then for the 1-D-VBS and all bins constrained at nC = 10, significant RH-induced enhancement of OPM condensation was observed in the distributions. The spread of the distributions was found to increase rapidly as the level of OPM decreased. The within-bin spread of ±0.5 log units in the groups

  8. Behaviour of tracer diffusion in simple atmospheric boundary layer models

    P. S. Anderson

    2006-12-01

    Full Text Available 1-D profiles and time series from an idealised atmospheric boundary layer model are presented, which show agreement with measurements of polar photogenic NO and NO2. Diffusion models are increasingly being used as the framework for studying tropospheric air chemistry dynamics. Models based on standard boundary layer diffusivity profiles have an intrinsic behaviour that is not necessarily intuitive, due to the variation of turbulent diffusivity with height. The relatively simple model provides both a programming and a conceptual tool in the analysis of observed trace gas evolution. A time scale inherent in the model can be tuned by fitting model time series to observations. This scale is then applicable to the more physically simple but chemically complex zeroth order or box models of chemical interactions.

  9. Behaviour of tracer diffusion in simple atmospheric boundary layer models

    P. S. Anderson

    2007-10-01

    Full Text Available 1-D profiles and time series from an idealised atmospheric boundary layer model are presented, which show agreement with boundary layer measurements of polar NOx. Diffusion models are increasingly being used as the framework for studying tropospheric air chemistry dynamics. Models based on standard boundary layer diffusivity profiles have an intrinsic behaviour that is not necessarily intuitive, due to the variation of turbulent diffusivity with height. The simple model presented captures the essence of the evolution of a trace gas released at the surface, and thereby provides both a programming and a conceptual tool in the analysis of observed trace gas evolution. A time scale inherent in the model can be tuned by fitting model time series to observations. This scale is then applicable to the more physically simple but chemically complex zeroth order or box models of chemical interactions.

  10. Relocation of Earthquakes in the Northeastern Tianshan Mountains Area and Improvement of Local 1-D Crustal Velocity ModelI

    Sun Anhui; Chen Qifu; Chen Yong; Li Gang; Zhang Zhiqiang

    2012-01-01

    We apply three methods to relocate 599 earthquake events that occurred from August 2004 to August 2005 in the northeastern Tianshan Mountains area (85~30'- 88~30'E, 43~00'- 44~40'N) by using travel times recorded by regional seismic network and 10 portable seismic stations deployed around the Urumqi city. By comparing the reliability of different results, we determined a suitable location method, and an improved 1-D crustal velocity model of the study area. The uncertainty of earthquake location is significantly reduced with combined data of seismic network and portable stations. The relocated events are clearly associated with regional tectonics of the northeastern Tianshan Mountains area, and are also in agreement with the existence of active faults imaged by deep seismic reflection profile. The relocated seismicity discovers some potential traces of buried active faults, which need to be validated further.

  11. Salient Features of High-Energy Multiparticle Distributions Learned from Exact Solutions of 1-d Ising Model

    Chau, L L; Chau, Ling-Lie; Huang, Ding-Wei

    1993-01-01

    We have derived explicit expressions in the 1-d Ising model for multiplicity distributions $P_{\\del\\xi}(n)$ and factorial moments $F_q(\\del\\xi)$. We identify the salient features of $P_{\\del\\xi}(n)$ that lead to scaling, $F_q(\\del\\xi)=\\F_q[F_2]$, and universality. These results compare well with the presently available high-energy data of $\\bar{p}p$ and $e^+e^-$ reactions. We point out the important features that should be studied in future higher-energy experiments of multiparticle productions in $pp$, $\\bar{p}p$, $ep$, $e^+e^-$, and $NN$. We also make comments on comparisons with KNO and negative-binomial distributions.

  12. 2D Axisymmetric vs 1D: A PIC/DSMC Model of Breakdown in Triggered Vacuum Spark Gaps

    Moore, Stan; Moore, Chris; Boerner, Jeremiah

    2015-09-01

    Last year at GEC14, we presented results of one-dimensional PIC/DSMC simulations of breakdown in triggered vacuum spark gaps. In this talk, we extend the model to two-dimensional axisymmetric and compare the results to the previous 1D case. Specially, we vary the fraction of the cathode that emits electrons and neutrals (holding the total injection rates over the cathode surface constant) and show the effects of the higher dimensionality on the time to breakdown. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U. S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

  13. Mapping temporal extent of Chiang Mai floods using coupled 1-D and quasi 2-D floodplain inundation models

    Kowit Boonrawd

    2015-04-01

    Full Text Available A coupling of a 1-D flood routing model and quasi 2-D floodplain inundation model is applied for mapping spacetime flood extent. The routing model is formulated based on a non-linear storage-discharge relationship which is converted from an observed and synthetic rating curve. To draw the rating curve, required parameters for each reaches are estimated from hydraulic properties, floodplain geometry and vegetation and building cover of compound channels. The shape of the floodplain is defined by using fitting exercise based on the reverse approach between past and simulated inundation flood extent, to solve the current problem of inadequate topographic input data for floodplain. Mapping of daily flood can be generated relying on flat water levels. The quasi 2-D raster model is tested and applied to generate more realistic water surface and is used to estimate flood extent. The model is applied to the floodplains of Chiang Mai, north of Thailand and used to estimate a time series of hourly flood maps. Extending from daily to hourly flood extent, mapping development provides more details of flood inundation extent and depth.

  14. A Dual EnKF for Estimating Water Level, Bottom Roughness, and Bathymetry in a 1-D Hydrodynamic Model

    Hooshyar, Milad; Wang, Dingbao; Hagen, Scott C

    2016-01-01

    Data assimilation has been applied to coastal hydrodynamic models to better estimate system states or parameters by incorporating observed data into the model. Kalman Filter (KF) is one of the most studied data assimilation methods whose application is limited to linear systems. For nonlinear systems such as hydrodynamic models a variation of the KF called Ensemble Kalman Filter (EnKF) is applied to update the system state in the context of Monte Carlo simulation. In this research, a dual EnKF approach is used to simultaneously estimate state (water surface elevation) and parameters (bottom roughness and bathymetry) of the shallow water models. The sensitivity of the filter to 1) the quantity and precision of the observations, and 2) the initial estimation of parameters is investigated in a 1-D shallow water problem located in the Gulf of Mexico. Results show that starting from an initial estimate of bottom roughness and bathymetry within a logical range and utilizing observations available at a limited numbe...

  15. Estimation of exhaust gas aerodynamic force on the variable geometry turbocharger actuator: 1D flow model approach

    Highlights: • Estimation of aerodynamic force on variable turbine geometry vanes and actuator. • Method based on exhaust gas flow modeling. • Simulation tool for integration of aerodynamic force in automotive simulation software. - Abstract: This paper provides a reliable tool for simulating the effects of exhaust gas flow through the variable turbine geometry section of a variable geometry turbocharger (VGT), on flow control mechanism. The main objective is to estimate the resistive aerodynamic force exerted by the flow upon the variable geometry vanes and the controlling actuator, in order to improve the control of vane angles. To achieve this, a 1D model of the exhaust flow is developed using Navier–Stokes equations. As the flow characteristics depend upon the volute geometry, impeller blade force and the existing viscous friction, the related source terms (losses) are also included in the model. In order to guarantee stability, an implicit numerical solver has been developed for the resolution of the Navier–Stokes problem. The resulting simulation tool has been validated through comparison with experimentally obtained values of turbine inlet pressure and the aerodynamic force as measured at the actuator shaft. The simulator shows good compliance with experimental results

  16. A Generic 1D Forward Modeling and Inversion Algorithm for TEM Sounding with an Arbitrary Horizontal Loop

    Li, Zhanhui; Huang, Qinghua; Xie, Xingbing; Tang, Xingong; Chang, Liao

    2016-08-01

    We present a generic 1D forward modeling and inversion algorithm for transient electromagnetic (TEM) data with an arbitrary horizontal transmitting loop and receivers at any depth in a layered earth. Both the Hankel and sine transforms required in the forward algorithm are calculated using the filter method. The adjoint-equation method is used to derive the formulation of data sensitivity at any depth in non-permeable media. The inversion algorithm based on this forward modeling algorithm and sensitivity formulation is developed using the Gauss-Newton iteration method combined with the Tikhonov regularization. We propose a new data-weighting method to minimize the initial model dependence that enhances the convergence stability. On a laptop with a CPU of i7-5700HQ@3.5 GHz, the inversion iteration of a 200 layered input model with a single receiver takes only 0.34 s, while it increases to only 0.53 s for the data from four receivers at a same depth. For the case of four receivers at different depths, the inversion iteration runtime increases to 1.3 s. Modeling the data with an irregular loop and an equal-area square loop indicates that the effect of the loop geometry is significant at early times and vanishes gradually along the diffusion of TEM field. For a stratified earth, inversion of data from more than one receiver is useful in noise reducing to get a more credible layered earth. However, for a resistive layer shielded below a conductive layer, increasing the number of receivers on the ground does not have significant improvement in recovering the resistive layer. Even with a down-hole TEM sounding, the shielded resistive layer cannot be recovered if all receivers are above the shielded resistive layer. However, our modeling demonstrates remarkable improvement in detecting the resistive layer with receivers in or under this layer.

  17. A COUPLED 1-D AND 2-D CHANNEL NETWORK MATHEMATICAL MODEL USED FOR FLOW CALCULATIONS IN THE MIDDLE REACHES OF THE YANGTZE RIVER

    HAN Dong; FANG Hong-wei; BAI Jing; HE Guo-jian

    2011-01-01

    A coupled one-dimensional(1-D)and two-dimensional(2-D)channel network mathematical model is proposed for flow calculations at nodes in a channel network system in this article.For the 1-D model,the finite difference method is used to discretize the Saint-Venant equations in all channels of a looped network.The Alternating Direction Implicit(ADI)method is adopted for the 2-D model at the nodes.In the coupled model,the 1-D model provides a good approximation with small computational effort,while the 2-D model is applied for complex topography to achieve a high accuracy.An Artificial Neural Network(ANN)method is used for the data exchange and the connectivity between the 1-D and 2-D models.The coupled model is applied to the Jingjiang-Dongting Lake region,to simulate the tremendous looped channel network system,and the results are compared with field data.The good agreement shows that the coupled hydraulic model is more effective than the conventional 1-D model.

  18. Comparison between a coupled 1D-2D model and a fully 2D model for supercritical flow simulation in crossroads

    Ghostine, Rabih

    2014-12-01

    In open channel networks, flow is usually approximated by the one-dimensional (1D) Saint-Venant equations coupled with an empirical junction model. In this work, a comparison in terms of accuracy and computational cost between a coupled 1D-2D shallow water model and a fully two-dimensional (2D) model is presented. The paper explores the ability of a coupled model to simulate the flow processes during supercritical flows in crossroads. This combination leads to a significant reduction in the computational time, as a 1D approach is used in branches and a 2D approach is employed in selected areas only where detailed flow information is essential. Overall, the numerical results suggest that the coupled model is able to accurately simulate the main flow processes. In particular, hydraulic jumps, recirculation zones, and discharge distribution are reasonably well reproduced and clearly identified. Overall, the proposed model leads to a 30% reduction in run times. © 2014 International Association for Hydro-Environment Engineering and Research.

  19. Modelling the angular momentum J of 1s, 1p, 1d, 2s and 1f nucleons

    By using the liquid drop model of 14 alpha particles representing a nickel 56 nuclide it can be shown that the mean distance of each of the 1d and 2s nucleons is r3 = 2.85 fm from the nuclide centre. It was found that the velocity of all nucleons is the same and is independent of the energy level. This implies that the de Broglie wavelength (w) of all nucleons is w h / m v = 6.3 fm ∼ 2π fm . Therefore for r1 ∼ 1 fm there is one w per orbit; for r2 ∼ 2 fm there are 2 w per orbit and so on. This implies that in the first magic number closed shell of nucleons there are 2 orbits each containing 2 standing wave maxima representing 1 proton and 1 neutron. The second closed shell consists of 3 orbits each containing 2 proton and 2 neutron standing wave maxima. While the third closed shell consists of 4 orbits each containing 3 protons and 3 neutrons the fourth closed shell consists of only 2 orbits each containing 4 protons and 4 neutrons. The Bernal liquid drop alpha particle models of nuclear structure appear to accord quite well with the quantum mechanical prescriptions of nucleon angular momentum and de Broglie wavelength

  20. High resolution transmission spectroscopy as a diagnostic for Jovian exoplanet atmospheres: constraints from theoretical models

    Kempton, Eliza M.-R. [Department of Physics, Grinnell College, Grinnell, IA 50112 (United States); Perna, Rosalba [Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY 11794 (United States); Heng, Kevin, E-mail: kemptone@grinnell.edu [University of Bern, Center for Space and Habitability, Sidlerstrasse 5, CH-3012 Bern (Switzerland)

    2014-11-01

    We present high resolution transmission spectra of giant planet atmospheres from a coupled three-dimensional (3D) atmospheric dynamics and transmission spectrum model that includes Doppler shifts which arise from winds and planetary motion. We model Jovian planets covering more than two orders of magnitude in incident flux, corresponding to planets with 0.9-55 day orbital periods around solar-type stars. The results of our 3D dynamical models reveal certain aspects of high resolution transmission spectra that are not present in simple one-dimensional (1D) models. We find that the hottest planets experience strong substellar to anti-stellar (SSAS) winds, resulting in transmission spectra with net blueshifts of up to 3 km s{sup –1}, whereas less irradiated planets show almost no net Doppler shifts. We find only minor differences between transmission spectra for atmospheres with temperature inversions and those without. Compared to 1D models, peak line strengths are significantly reduced for the hottest atmospheres owing to Doppler broadening from a combination of rotation (which is faster for close-in planets under the assumption of tidal locking) and atmospheric winds. Finally, high resolution transmission spectra may be useful in studying the atmospheres of exoplanets with optically thick clouds since line cores for very strong transitions should remain optically thick to very high altitude. High resolution transmission spectra are an excellent observational test for the validity of 3D atmospheric dynamics models, because they provide a direct probe of wind structures and heat circulation. Ground-based exoplanet spectroscopy is currently on the verge of being able to verify some of our modeling predictions, most notably the dependence of SSAS winds on insolation. We caution that interpretation of high resolution transmission spectra based on 1D atmospheric models may be inadequate, as 3D atmospheric motions can produce a noticeable effect on the absorption

  1. A 1D-ecosystem model for pelagic waters in the southern Baltic Sea. Numerical simulations (future decades)

    Dzierzbicka-Glowacka, L.; Maciejewska, A.; Osiński, R.; Jakacki, J.; Jędrasik, J.

    2009-04-01

    This paper presents a one-dimensional Ecosystem Model. Mathematically, the pelagic variables in the model are described by a second-order partial differential equation of the diffusion type with biogeochemical sources and sinks. The temporal changes in the phytoplankton biomass are caused by primary production, respiration, mortality, grazing by zooplankton and sinking. The zooplankton biomass is affected by ingestion, excretion, respiration, fecal production, mortality, and carnivorous grazing. The changes in the pelagic detritus concentration are determined by input of: dead phytoplankton and zooplankton, natural mortality of predators, fecal pellets, and sinks: sedimentation, zooplankton grazing and decomposition. The nutrient concentration is caused by nutrient release, zooplankton excretion, predator excretion, detritus decomposition and benthic regeneration as sources and by nutrient uptake by phytoplankton as sinks. However, the benthic detritus is described by phytoplankton sedimentation, detritus sedimentation and remineralisation. The particulate organic carbon concentration is determined as the sum of phytoplankton, zooplankton and dead organic matter (detritus) concentrations. The 1D ecosystem model was used to simulate the seasonal dynamics of pelagic variables (phytoplankton, zooplankton, pelagic detritus and POC) in the southern Baltic Sea (Gdańsk Deep, Bornholm Deep and Gotland Deep). The calculations were made assuming: 1) increase in the water temperature in the upper layer - 0.008oC per year, 2) increase in the available light - 0.2% per year. Based on this trend, daily, monthly and seasonal and annual variability of phytoplankton, zooplankton, pelagic detritus and particulate organic carbon in different areas of the southern Baltic Sea (Gdańsk Deep, Borrnholm Deep and Gotland Deep) in the euphotic layer was calculated for the years: 2000, 2010, 2020, 2030, 2040 and 2050.

  2. Assessing the impact of uncertainty on flood risk estimates with reliability analysis using 1-D and 2-D hydraulic models

    L. Altarejos-García

    2012-07-01

    Full Text Available This paper addresses the use of reliability techniques such as Rosenblueth's Point-Estimate Method (PEM as a practical alternative to more precise Monte Carlo approaches to get estimates of the mean and variance of uncertain flood parameters water depth and velocity. These parameters define the flood severity, which is a concept used for decision-making in the context of flood risk assessment. The method proposed is particularly useful when the degree of complexity of the hydraulic models makes Monte Carlo inapplicable in terms of computing time, but when a measure of the variability of these parameters is still needed. The capacity of PEM, which is a special case of numerical quadrature based on orthogonal polynomials, to evaluate the first two moments of performance functions such as the water depth and velocity is demonstrated in the case of a single river reach using a 1-D HEC-RAS model. It is shown that in some cases, using a simple variable transformation, statistical distributions of both water depth and velocity approximate the lognormal. As this distribution is fully defined by its mean and variance, PEM can be used to define the full probability distribution function of these flood parameters and so allowing for probability estimations of flood severity. Then, an application of the method to the same river reach using a 2-D Shallow Water Equations (SWE model is performed. Flood maps of mean and standard deviation of water depth and velocity are obtained, and uncertainty in the extension of flooded areas with different severity levels is assessed. It is recognized, though, that whenever application of Monte Carlo method is practically feasible, it is a preferred approach.

  3. Assessing the impact of uncertainty on flood risk estimates with reliability analysis using 1-D and 2-D hydraulic models

    L. Altarejos-García

    2012-01-01

    Full Text Available This paper addresses the use of reliability techniques such as Rosenblueth's Point-Estimate Method (PEM as a practical alternative to more precise Monte Carlo approaches to get estimates of the mean and variance of uncertain flood parameters water depth and velocity. These parameters define the flood severity, which is a concept used for decision-making in the context of flood risk assessment. The method proposed is particularly useful when the degree of complexity of the hydraulic models makes Monte Carlo inapplicable in terms of computing time, but when a measure of the variability of these parameters is still needed. The capacity of PEM, which is a special case of numerical quadrature based on orthogonal polynomials, to evaluate the first two moments of performance functions such as the water depth and velocity is demonstrated in the case of a single river reach using a 1-D HEC-RAS model. It is shown that in some cases, using a simple variable transformation, statistical distributions of both water depth and velocity approximate the lognormal. As this distribution is fully defined by its mean and variance, PEM can be used to define the full probability distribution function of these flood parameters and so allowing for probability estimations of flood severity. Then, an application of the method to the same river reach using a 2-D Shallow Water Equations (SWE model is performed. Flood maps of mean and standard deviation of water depth and velocity are obtained, and uncertainty in the extension of flooded areas with different severity levels is assessed. It is recognized, though, that whenever application of Monte Carlo method is practically feasible, it is a preferred approach.

  4. Development of simplified 1D and 2D models for studying a PWR lower head failure under severe accident conditions

    In the study of severe accidents of nuclear pressurized water reactors, the scenarios that describe the relocation of significant quantities of liquid corium at the bottom of the lower head are investigated from the mechanical point of view. In these scenarios, the risk of a breach and the possibility of a large quantity of corium being released from the lower head exist. This may lead to direct heating of the containment or outer vessel steam explosion. These issues are important due to their early containment failure potential. Since the TMI-2 accident, many theoretical and experimental investigations, relating to lower head mechanical behaviour under severe thermo-mechanical loading in the event of a core meltdown accident have been performed. IRSN participated actively in the one-fifth scale USNRC/SNL LHF and OECD LHF (OLHF) programs. Within the framework of these programs, two simplified models were developed by IRSN: the first is a simplified 1D approach based on the theory of pressurized spherical shells and the second is a simplified 2D model based on the theory of shells of revolution under symmetric loading. The mathematical formulation of both models and the creep constitutive equations used are presented in detail in this paper. The corresponding models were used to interpret some of the OLHF program experiments and the calculation results were quite consistent with the experimental data. The two simplified models have been used to simulate the thermo-mechanical behaviour of a 900 MWe pressurized water reactor lower head under severe accident conditions leading to failure. The average transient heat flux produced by the corium relocated at the bottom of the lower head has been determined using the IRSN HARAR code. Two different methods, both taking into account the ablation of the internal surface, are used to determine the temperature profiles across the lower head wall and their effect on the time to failure is discussed. Using these simplified models

  5. Soil-vegetation-atmosphere transfer modeling

    Ikonen, J.P.; Sucksdorff, Y. [Finnish Environment Agency, Helsinki (Finland)

    1996-12-31

    In this study the soil/vegetation/atmosphere-model based on the formulation of Deardorff was refined to hour basis and applied to a field in Vihti. The effect of model parameters on model results (energy fluxes, temperatures) was also studied as well as the effect of atmospheric conditions. The estimation of atmospheric conditions on the soil-vegetation system as well as an estimation of the effect of vegetation parameters on the atmospheric climate was estimated. Areal surface fluxes, temperatures and moistures were also modelled for some river basins in southern Finland. Land-use and soil parameterisation was developed to include properties and yearly variation of all vegetation and soil types. One classification was selected to describe the hydrothermal properties of the soils. Evapotranspiration was verified against the water balance method

  6. CD1d-dependent NKT cells play a protective role in acute and chronic arthritis models by ameliorating antigen-specific Th1 responses

    Teige, Anna; Bockermann, Robert; Hasan, Maruf;

    2010-01-01

    A protective and anti-inflammatory role for CD1d-dependent NKT cells (NKTs) has been reported in experimental and human autoimmune diseases. However, their role in arthritis has been unclear, with conflicting reports of CD1d-dependent NKTs acting both as regulatory and disease-promoting cells...... in arthritis. These differing modes of action might be due to genetic differences of inbred mice and incomplete backcrossing of gene-modified mice. We therefore put special emphasis on controlling the genetic backgrounds of the mice used. Additionally, we used two different murine arthritis models, Ag......-induced arthritis (AIA) and collagen-induced arthritis (CIA), to evaluate acute and chronic arthritis in CD1d knockout mice and mice depleted of NK1.1(+) cells. CD1d-deficient mice developed more severe AIA compared with wild-type littermates, with a higher degree of inflammation and proteoglycan depletion. Chronic...

  7. Modelling the atmospheric chemistry of volcanic plumes

    Surl, Luke

    2016-01-01

    Abstract Volcanoes are the principal way by which volatiles are transferred from the solid Earth to the atmosphere-hydrosphere system. Once released into the atmosphere, volcanic emissions rapidly undergo a complex series of chemical reactions. This thesis seeks to further the understanding of such processes by both observation and numerical modelling. I have adapted WRF-Chem to model passive degassing from Mount Etna, the chemistry of its plume, and its influence on the ...

  8. Atmospheric pollution. From processes to modelling

    Air quality, greenhouse effect, ozone hole, chemical or nuclear accidents.. All these phenomena are tightly linked to the chemical composition of atmosphere and to the atmospheric dispersion of pollutants. This book aims at supplying the main elements of understanding of 'atmospheric pollutions': stakes, physical processes involved, role of scientific expertise in decision making. Content: 1 - classifications and scales: chemical composition of the atmosphere, vertical structure, time scales (transport, residence); 2 - matter/light interaction: notions of radiative transfer, application to the Earth's atmosphere; 3 - some elements about the atmospheric boundary layer: notion of scales in meteorology, atmospheric boundary layer (ABL), thermal stratification and stability, description of ABL turbulence, elements of atmospheric dynamics, some elements about the urban climate; 4 - notions of atmospheric chemistry: characteristics, ozone stratospheric chemistry, ozone tropospheric chemistry, brief introduction to indoor air quality; 5 - aerosols, clouds and rains: aerosols and particulates, aerosols and clouds, acid rains and leaching; 6 - towards numerical simulation: equation of reactive dispersion, numerical methods for chemistry-transport models, numerical resolution of the general equation of aerosols dynamics (GDE), modern simulation chains, perspectives. (J.S.)

  9. Development of a 1 D hydrodynamic habitat model for the Hippopotamus amphibious as basis for sustainable exploitation of hydroelectric power

    Manful, D. Y.; Kaule, G.; Wieprecht, S.; Rees, J.; Hu, W.

    2009-12-01

    Hydroelectric Power (HEP) is proving to be a good alternative to carbon based energy. In the past hydropower especially large scale hydro attracted significant criticism as a result of its impact on the environment. A new breed of hydroelectric dam is in the offing. The aim is to have as little a footprint as possible on the environment in both pre and post construction phases and thus minimize impact on biodiversity whilst producing clean renewable energy. The Bui dam is 400 MW scheme currently under development on the Black Volta River in the Bui national park in Ghana. The reservoir created by the Bui barrage is expected to impact (through inundation) the habitat of two species of hippos know to exist in the park, the Hippopotamus amphibius and the Choeropsis liberiensis. Computer-based models present a unique opportunity to assess quantitatively the impact of the new reservoir on the habitat of the target species in this case the H. amphibious. Until this undertaking, there were very few studies documenting the habitat of the H. amphibious let alone model it. The work and subsequent presentation will show the development of a habitat model for the Hippopotamus amphibius. The Habitat Information retrieval Program based on Streamflow Analysis, in short HIPStrA, is a one dimensional (1D) in-stream, spatially explicit hybrid construct that combines physico-chemical evidence and expert knowledge to forecast river habitat suitability (Hs) for the Hippopotamus amphibius. The version of the model presented is specifically developed to assess the impact of a reservoir created by a hydroelectric dam on potential dwelling areas in the Bui gorge for hippos. Accordingly, this version of HIPStrA simulates a special reservoir suitability index (Rsi), a metric that captures the”hippo friendliness” of any lake or reservoir. The impact of measured and simulated flood events as well as low flows, representing extreme events is also assessed. Recommendations are made for the

  10. Regional forecasting with global atmospheric models

    The scope of the report is to present the results of the fourth year's work on the atmospheric modeling part of the global climate studies task. The development testing of computer models and initial results are discussed. The appendices contain studies that provide supporting information and guidance to the modeling work and further details on computer model development. Complete documentation of the models, including user information, will be prepared under separate reports and manuals

  11. Users of middle atmosphere models remarks

    Gamble, Joe

    1987-01-01

    The procedure followed for shuttle operations is to calculate descent trajectories for each potential shuttle landing site using the Global Reference Atmosphere Model (GRAM) to interactively compute density along the flight path 100 times to bound the statistics. The purpose is to analyze the flight dynamics, along with calculations of heat loads during reentry. The analysis program makes use of the modified version of the Jacchia-70 atmosphere, which includes He bulges over the poles and seasonal latitude variations at lower altitudes. For the troposphere, the 4-D Model is used up to 20 km, Groves from 30 km up to 90 km. It is extrapolated over the globe and faired into the Jacchia atmosphere between 90 and 115 km. Since data on the Southern Hemisphere was lacking, what was done was that the data was flipped over and lagged 6 months. Sometimes when winds are calculated from pressure data in the model there appear to be discontinuities. Modelers indicated that the GRAM was not designed to produce winds, but good wind data is needed for the landing phase of shuttle operations. Use of atmospheric models during reentry is one application where it is obvious that a single integrated atmosphere model is required.

  12. Stellar model atmospheres with magnetic line blanketing

    Kochukhov, O; Shulyak, D

    2004-01-01

    Model atmospheres of A and B stars are computed taking into account magnetic line blanketing. These calculations are based on the new stellar model atmosphere code LLModels which implements direct treatment of the opacities due to the bound-bound transitions and ensures an accurate and detailed description of the line absorption. The anomalous Zeeman effect was calculated for the field strengths between 1 and 40 kG and a field vector perpendicular to the line of sight. The model structure, high-resolution energy distribution, photometric colors, metallic line spectra and the hydrogen Balmer line profiles are computed for magnetic stars with different metallicities and are discussed with respect to those of non-magnetic reference models. The magnetically enhanced line blanketing changes the atmospheric structure and leads to a redistribution of energy in the stellar spectrum. The most noticeable feature in the optical region is the appearance of the 5200 A depression. However, this effect is prominent only in ...

  13. New atmospheric model of Epsilon Eridani

    Vieytes, Mariela; Fontenla, Juan; Buccino, Andrea; Mauas, Pablo

    2016-05-01

    We present a new semi-empirical model of the atmosphere of the widely studied K-dwarf Epsilon Eridani (HD 22049). The model is build to reproduce the visible spectral observations from 3800 to 6800 Angstrom and the h and k Mg II lines profiles. The computations were carried out using the Solar-Stellar Radiation Physical Modeling (SSRPM) tools, which calculate non-LTE population for the most important species in the stellar atmosphere. We show a comparison between the synthetic and observed spectrum, obtaining a good agreement in all the studied spectral range.

  14. Establishing the Capability of a 1D SVAT Modelling Scheme in Predicting Key Biophysical Vegetation Characterisation Parameters

    Ireland, Gareth; Petropoulos, George P.; Carlson, Toby N.; Purdy, Sarah

    2015-04-01

    Sensitivity analysis (SA) consists of an integral and important validatory check of a computer simulation model before it is used to perform any kind of analysis. In the present work, we present the results from a SA performed on the SimSphere Soil Vegetation Atmosphere Transfer (SVAT) model utilising a cutting edge and robust Global Sensitivity Analysis (GSA) approach, based on the use of the Gaussian Emulation Machine for Sensitivity Analysis (GEM-SA) tool. The sensitivity of the following model outputs was evaluated: the ambient CO2 concentration and the rate of CO2 uptake by the plant, the ambient O3 concentration, the flux of O3 from the air to the plant/soil boundary, and the flux of O3 taken up by the plant alone. The most sensitive model inputs for the majority of model outputs were related to the structural properties of vegetation, namely, the Leaf Area Index, Fractional Vegetation Cover, Cuticle Resistance and Vegetation Height. External CO2 in the leaf and the O3 concentration in the air input parameters also exhibited significant influence on model outputs. This work presents a very important step towards an all-inclusive evaluation of SimSphere. Indeed, results from this study contribute decisively towards establishing its capability as a useful teaching and research tool in modelling Earth's land surface interactions. This is of considerable importance in the light of the rapidly expanding use of this model worldwide, which also includes research conducted by various Space Agencies examining its synergistic use with Earth Observation data towards the development of operational products at a global scale. This research was supported by the European Commission Marie Curie Re-Integration Grant "TRANSFORM-EO". SimSphere is currently maintained and freely distributed by the Department of Geography and Earth Sciences at Aberystwyth University (http://www.aber.ac.uk/simsphere). Keywords: CO2 flux, ambient CO2, O3 flux, SimSphere, Gaussian process emulators

  15. Constitutive model for flake graphite cast iron automotive brake discs: from macroscopic multiscale models to a 1D rheological description

    Augustins, L.; Billardon, R.; Hild, F.

    2016-07-01

    One of the critical points of the thermomechanical fatigue design process is the correct description of the cyclic behavior of the material. This work focuses on the material of automotive brake discs, namely flake graphite cast iron. The specificity of this material is its asymmetric behavior under tensile and compressive loadings, which is due to the shape of graphite that acts as small cracks. Multiscale models inspired from the literature are first presented. They lead to a good description of the material behavior under cyclic loadings. An elastoviscoplastic constitutive model is then proposed in a one-dimensional setting in order to accurately describe cyclic tests from room temperature up to {600^{circ}{C}}.

  16. Model-Based Least Squares Reconstruction of Coded Source Neutron Radiographs: Integrating the ORNL HFIR CG1D Source Model

    Santos-Villalobos, Hector J [ORNL; Gregor, Jens [University of Tennessee, Knoxville (UTK); Bingham, Philip R [ORNL

    2014-01-01

    At the present, neutron sources cannot be fabricated small and powerful enough in order to achieve high resolution radiography while maintaining an adequate flux. One solution is to employ computational imaging techniques such as a Magnified Coded Source Imaging (CSI) system. A coded-mask is placed between the neutron source and the object. The system resolution is increased by reducing the size of the mask holes and the flux is increased by increasing the size of the coded-mask and/or the number of holes. One limitation of such system is that the resolution of current state-of-the-art scintillator-based detectors caps around 50um. To overcome this challenge, the coded-mask and object are magnified by making the distance from the coded-mask to the object much smaller than the distance from object to detector. In previous work, we have shown via synthetic experiments that our least squares method outperforms other methods in image quality and reconstruction precision because of the modeling of the CSI system components. However, the validation experiments were limited to simplistic neutron sources. In this work, we aim to model the flux distribution of a real neutron source and incorporate such a model in our least squares computational system. We provide a full description of the methodology used to characterize the neutron source and validate the method with synthetic experiments.

  17. Pose Estimation using 1D Fourier Transform and Euclidean Distance Matching of CAD Model and Inspected Model Part

    Zulkoffli, Zuliani; Abu Bakar, Elmi

    2016-02-01

    This paper present pose estimation relation of CAD model object and Projection Real Object (PRI). Image sequence of PRI and CAD model rotate on z axis at 10 degree interval in simulation and real scene used in this experiment. All this image is go through preprocessing stage to rescale object size and image size and transform all the image into silhouette. Correlation of CAD and PRI image is going through in this stage. Magnitude spectrum shows a reliable value in range 0.99 to 1.00 and Phase spectrum correlation shows a fluctuate graph in range 0.56 - 0.97. Euclidean distance correlation graph for CAD and PRI shows 2 zone of similar value due to almost symmetrical object shape. Processing stage of retrieval inspected PRI image in CAD database was carried out using range phase spectrum and maximum magnitude spectrum value within ±10% tolerance. Additional processing stage of retrieval inspected PRI image using Euclidean distance within ±5% tolerance also carried out. Euclidean matching shows a reliable result compared to range phase spectrum and maximum magnitude spectrum value by sacrificing more than 5 times processing time.

  18. Coupled atmosphere-wildland fire modelling

    Jacques Henri Balbi

    2009-10-01

    Full Text Available Simulating the interaction between fire and atmosphere is critical to the estimation of the rate of spread of the fire. Wildfire’s convection (i.e., entire plume can modify the local meteorology throughout the atmospheric boundary layer and consequently affect the fire propagation speed and behaviour. In this study, we use for the first time the Méso-NH meso-scale numerical model coupled to the point functional ForeFire simplified physical front-tracking wildfire model to investigate the differences introduced by the atmospheric feedback in propagation speed and behaviour. Both numerical models have been developed as research tools for operational models and are currently used to forecast localized extreme events. These models have been selected because they can be run coupled and support decisions in wildfire management in France and Europe. The main originalities of this combination reside in the fact that Méso-NH is run in a Large Eddy Simulation (LES configuration and that the rate of spread model used in ForeFire provides a physical formulation to take into account the effect of wind and slope. Simulations of typical experimental configurations show that the numerical atmospheric model is able to reproduce plausible convective effects of the heat produced by the fire. Numerical results are comparable to estimated values for fire-induced winds and present behaviour similar to other existing numerical approaches.

  19. Atmospheric neutrino flux calculation using the NRLMSISE00 atmospheric model

    Honda, M; Kajita, T; Kasahara, K; Midorikawa, S

    2015-01-01

    In this paper, we extend the calculation of the atmospheric neutrino flux~\\cite{hkkm2004,hkkms2006,hkkm2011} to the sites in polar and tropical regions. In our earliest full 3D-calculation~\\cite{hkkm2004}, we used DPMJET-III~\\cite{dpm} for the hadronic interaction model above 5~GeV, and NUCRIN~\\cite{nucrin} below 5~GeV. We modified DPMJET-III as in Ref.~\\cite{hkkms2006} to reproduce the experimental muon spectra better, mainly using the data observed by BESS group~\\cite{BESSTeVpHemu}. In a recent work~\\cite{hkkm2011}, we introduced JAM interaction model for the low energy hadronic interactions. JAM is a nuclear interaction model developed with PHITS (Particle and Heavy-Ion Transport code System)~\\cite{phits}. In Ref.~\\cite{hkkm2011}, we could reproduce the observed muon flux at the low energies at balloon altitude with DPMJET-III above 32 GeV and JAM below that better than the combination of DPMJET-III above 5~GeV and NUCRIN below that. Besides the interaction model, we have also improved the calculation sche...

  20. Portable University Model of the Atmosphere (PUMA)

    Fraedrich, K.; Kirk, E.; Lunkeit, F. [Hamburg Univ. (Germany). Meteorologisches Inst.

    1998-10-01

    The Portable University Model of the Atmosphere (PUMA) is based on the Reading multi-level spectral model SGCM (Simple Global Circulation Model) described by Hoskins and Simmons (1975) and James and Gray (1986). Originally developed as a numerical prediction model, it was changed to perform as a circulation model. For example, James and Gray (1986) studied the influence of surface friction on the circulation of a baroclinic atmosphere, James and James (1992), and James et al. (1994) investigated ultra-low-frequency variability, and Mole and James (1990) analyzed the baroclinic adjustment in the context of a zonally varying flow. Frisius et al. (1998) simulated an idealized storm track by embedding a dipole structure in a zonally symmetric forcing field and Lunkeit et al. (1998) investigated the sensitivity of GCM (General Circulation Model) scenarios by an adaption technique applicapable to SGCMs. (orig.)

  1. Coupling approaches used in atmospheric entry models

    Gritsevich, M. I.

    2012-09-01

    While a planet orbits the Sun, it is subject to impact by smaller objects, ranging from tiny dust particles and space debris to much larger asteroids and comets. Such collisions have taken place frequently over geological time and played an important role in the evolution of planets and the development of life on the Earth. Though the search for near-Earth objects addresses one of the main points of the Asteroid and Comet Hazard, one should not underestimate the useful information to be gleaned from smaller atmospheric encounters, known as meteors or fireballs. Not only do these events help determine the linkages between meteorites and their parent bodies; due to their relative regularity they provide a good statistical basis for analysis. For successful cases with found meteorites, the detailed atmospheric path record is an excellent tool to test and improve existing entry models assuring the robustness of their implementation. There are many more important scientific questions meteoroids help us to answer, among them: Where do these objects come from, what are their origins, physical properties and chemical composition? What are the shapes and bulk densities of the space objects which fully ablate in an atmosphere and do not reach the planetary surface? Which values are directly measured and which are initially assumed as input to various models? How to couple both fragmentation and ablation effects in the model, taking real size distribution of fragments into account? How to specify and speed up the recovery of a recently fallen meteorites, not letting weathering to affect samples too much? How big is the pre-atmospheric projectile to terminal body ratio in terms of their mass/volume? Which exact parameters beside initial mass define this ratio? More generally, how entering object affects Earth's atmosphere and (if applicable) Earth's surface? How to predict these impact consequences based on atmospheric trajectory data? How to describe atmospheric entry

  2. Radiation environment models and the atmospheric cutoff

    Konradi, Andrei; Hardy, Alva C.; Atwell, William

    1987-01-01

    The limitations of radiation environment models are examined by applying the model to the South Atlantic anomaly (SAA). The local magnetic-field-intensity (in gauss) and McIlwain (1961) drift-shell-parameter contours in the SAA are analyzed. It is noted that it is necessary to decouple the atmospheric absorption effects from the trapped radiation models in order to obtain accurate radiation dose predictions. Two methods for obtaining more accurate results are proposed.

  3. Fingering convection and cloudless models for cool brown dwarf atmospheres

    Tremblin, P; Mourier, P; Baraffe, I; Chabrier, G; Drummond, B; Homeier, D; Venot, O

    2015-01-01

    This work aims to improve the current understanding of the atmospheres of brown dwarfs, especially cold ones with spectral type T and Y, whose modeling is a current challenge. Silicate and iron clouds are believed to disappear at the photosphere at the L/T transition, but cloudless models fail to reproduce correctly the spectra of T dwarfs, advocating for the addition of more physics, e.g. other types of clouds or internal energy transport mechanisms. We use a one-dimensional (1D) radiative/convective equilibrium code ATMO to investigate this issue. This code includes both equilibrium and out-of-equilibrium chemistry and solves consistently the PT structure. Included opacity sources are H2-H2, H2-He, H2O, CO, CO2, CH4, NH3, K, Na, and TiO, VO if they are present in the atmosphere. We show that the spectra of Y dwarfs can be accurately reproduced with a cloudless model if vertical mixing and NH3 quenching are taken into account. T dwarf spectra still have some reddening in e.g. J - H compared to cloudless mode...

  4. ATMOSPHERIC HEALTH EFFECTS FRAMEWORK (AHEF) MODEL

    The Atmospheric and Health Effects Framework (AHEF) is used to assess theglobal impacts of substitutes for ozone-depleting substances (ODS). The AHEF is a series of FORTRAN modeling modules that collectively form a simulation framework for (a) translating ODS production into emi...

  5. Atmospheric statistical dynamic models. Model performance: the Lawrence Livermore Laboratoy Zonal Atmospheric Model

    Results from the zonal model indicate quite reasonable agreement with observation in terms of the parameters and processes that influence the radiation and energy balance calculations. The model produces zonal statistics similar to those from general circulation models, and has also been shown to produce similar responses in sensitivity studies. Further studies of model performance are planned, including: comparison with July data; comparison of temperature and moisture transport and wind fields for winter and summer months; and a tabulation of atmospheric energetics. Based on these preliminary performance studies, however, it appears that the zonal model can be used in conjunction with more complex models to help unravel the problems of understanding the processes governing present climate and climate change. As can be seen in the subsequent paper on model sensitivity studies, in addition to reduced cost of computation, the zonal model facilitates analysis of feedback mechanisms and simplifies analysis of the interactions between processes

  6. Fluctuation Theorem in an Atmospheric Circulation Model

    Schalge, Bernd; Wouters, Jeroen; Fraedrich, Klaus; Lunkeit, Frank

    2012-01-01

    Evidence for the validity of the Fluctuation Theorem (FT) in an atmospheric Global Circulation Model is found. The model is hydrostatic with variable numbers of vertical levels and different horizontal resolutions. For finite time intervals the largest local Lyapunov exponent (LLLE) is found to be negative consistent with predictions of the FT. The effect is present for resolutions up to wave numbers l=42 (~ 250km) and 10 levels.

  7. Development of a global 1-D chemically radiatively coupled model and an introduction to the development of a chemically coupled General Circulation Model

    A global one-dimensional, chemically and radiatively coupled model has been developed. The basic concept of the coupled model, definition of globally averaged zenith angles, the formulation of the model chemistry, radiation, the coupled processes, and profiles and diurnal variations of temperature and chemical species at a normal steady state are presented. Furthermore, a suddenly doubled CO2 experiment and a Pinatubo aerosol increase experiment were performed with the model. The time scales of variations in ozone and temperature in the lower stratosphere of the coupled system in the doubled CO2 experiment was long, due to a feedback process among ultra violet radiation, O(1D), NOy, NOx, and O3. From the Pinatubo aerosol experiment, a delay of maximum ozone decrease from the maximum aerosol loading is shown and discussed. Developments of 3-D chemical models with coupled processes are briefly described, and the ozone distribution from the first version of the 3-D model are presented. Chemical model development in National Institute for Environmental Studies (NIES) are briefly described. (author)

  8. Dynamics of ozone and nitrogen oxides at Summit, Greenland. II. Simulating snowpack chemistry during a spring high ozone event with a 1-D process-scale model

    Murray, Keenan A.; Kramer, Louisa J.; Doskey, Paul V.; Ganzeveld, Laurens; Seok, Brian; Van Dam, Brie; Helmig, Detlev

    2015-09-01

    Observed depth profiles of nitric oxide (NO), nitrogen dioxide (NO2), and ozone (O3) in snowpack interstitial air at Summit, Greenland were best replicated by a 1-D process-scale model, which included (1) geometrical representation of snow grains as spheres, (2) aqueous-phase chemistry confined to a quasi-liquid layer (QLL) on the surface of snow grains, and (3) initialization of the species concentrations in the QLL through equilibrium partitioning with mixing ratios in snowpack interstitial air. A comprehensive suite of measurements in and above snowpack during a high O3 event facilitated analysis of the relationship between the chemistry of snowpack and the overlying atmosphere. The model successfully reproduced 2 maxima (i.e., a peak near the surface of the snowpack at solar noon and a larger peak occurring in the evening that extended down from 0.5 to 2 m) in the diurnal profile of NO2 within snowpack interstitial air. The maximum production rate of NO2 by photolysis of nitrate (NO3-) was approximately 108 molec cm-3 s-1, which explained daily observations of maxima in NO2 mixing ratios near solar noon. Mixing ratios of NO2 in snowpack interstitial air were greatest in the deepest layers of the snowpack at night and were attributed to thermal decomposition of peroxynitric acid, which produced up to 106 molec NO2 cm-3 s-1. Highest levels of NO in snowpack interstitial air were confined to upper layers of the snowpack and observed profiles were consistent with photolysis of NO2. Production of nitrogen oxides (NOx) from NO3- photolysis was estimated to be two orders of magnitude larger than NO production and supports the hypothesis that NO3- photolysis is the primary source of NOx within sunlit snowpack in the Arctic. Aqueous-phase oxidation of formic acid by O3 resulted in a maximum consumption rate of ∼106-107 molec cm-3 s-1 and was the primary removal mechanism for O3.

  9. Regional transport model of atmospheric sulfates

    As part of the Sulfate Regional Experiment (SURE) Design Project, a regional transport model of atmospheric sulfates has been developed. This quasi-Lagrangian three-dimensional grid numerical model uses a detailed SO2 emission inventory of major anthropogenic sources in the Eastern U.S. region, and observed meteorological data during an episode as inputs. The model accounts for advective transport and turbulent diffusion of the pollutants. The chemical transformation of SO2 and SO4/sup =/ and the deposition of the species at the earth's surface are assumed to be linear processes at specified constant rates. The numerical model can predict the daily average concentrations of SO2 and SO4/sup =/ at all receptor locations in the grid region during the episode. Because of the spatial resolution of the grid, this model is particularly suited to investigate the effect of tall stacks in reducing the ambient concentration levels of sulfur pollutants. This paper presents the formulations and assumptions of the regional sulfate transport model. The model inputs and results are discussed. Isopleths of predicted SO2 and SO4/sup =/ concentrations are compared with the observed ground level values. The bulk of the information in this paper is directed to air pollution meteorologists and environmental engineers interested in the atmospheric transport modeling studies of sulfur oxide pollutants

  10. Effects of 1-D versus 3-D velocity models on moment tensor inversion in the Dobrá Voda area in the Little Carpathians region, Slovakia

    Jechumtálová, Zuzana; Bulant, P.

    2014-01-01

    Roč. 18, č. 3 (2014), s. 511-531. ISSN 1383-4649 R&D Projects: GA ČR GAP210/10/1728 EU Projects: European Commission(XE) 230669 - AIM Institutional support: RVO:67985530 Keywords : ray tracing * 1-D and 3-D velocity models * earthquake mechanism * amplitude inversion * Dobrá Voda earthquakes Subject RIV: DC - Siesmology, Volcanology, Earth Structure Impact factor: 1.386, year: 2014

  11. Net community production of oxygen derived from in vitro and in situ 1-D modeling techniques in a cyclonic mesoscale eddy in the Sargasso Sea

    B. Mouriño-Carballido; Anderson, L A

    2009-01-01

    It has been proposed that the disagreement traditionally reported between in vitro incubation and in situ estimates of oxygen net community production (NCP) could be explained, at least partially, by undersampling episodic pulses of net autotrophy associated with mesoscale dynamics. In this study we compare in vitro incubation estimates of net community production with in situ estimates, derived from oxygen profiles and a 1-D model, within a cyclonic eddy investigated in the Sargasso Sea in s...

  12. Centrifuge modeling of soil atmosphere interaction

    CAICEDO, B; TRISTANCHO, J; THOREL, Luc

    2010-01-01

    Atmosphere process of infiltration or evaporation affect the behavior of geotechnical structures located near the soil surface. This paper focuses on the drying process of soils due to evaporation. The scaling laws are analyzed and afterwards the results on applying two cycles of heating and cooling on a soil mass are presented. Based on these results, conclusions about the feasibility of reproducing evaporation on centrifuge models are recommended.

  13. On the thermally stratified atmospheric flow modeling

    Sládek, Ivo; Kozel, K.; Jaňour, Zbyněk

    Praha : Ústav termomechaniky, AV ČR, v. v. i., 2010 - (Příhoda, J.; Kozel, K.), s. 135-138 ISBN 978-80-87012-25-3. [Topical Problems of Fluid Mechanics 2010. Praha (CZ), 10.02.2010-11.02.2010] Institutional research plan: CEZ:AV0Z20760514 Keywords : atmospheric boundary layer * turbulence model * finite volume method Subject RIV: DG - Athmosphere Sciences, Meteorology

  14. Observations and Modeling of Tropical Planetary Atmospheres

    Laraia, Anne

    2016-01-01

    This thesis is a comprised of three different projects within the topic of tropical atmospheric dynamics. First, I analyze observations of thermal radiation from Saturn's atmosphere and from them, determine the latitudinal distribution of ammonia vapor near the 1.5-bar pressure level. The most prominent feature of the observations is the high brightness temperature of Saturn's subtropical latitudes on either side of the equator. After comparing the observations to a microwave radiative transfer model, I find that these subtropical bands require very low ammonia relative humidity below the ammonia cloud layer in order to achieve the high brightness temperatures observed. We suggest that these bright subtropical bands represent dry zones created by a meridionally overturning circulation. Second, I use a dry atmospheric general circulation model to study equatorial superrotation in terrestrial atmospheres. A wide range of atmospheres are simulated by varying three parameters: the pole-equator radiative equilibrium temperature contrast, the convective lapse rate, and the planetary rotation rate. A scaling theory is developed that establishes conditions under which superrotation occurs in terrestrial atmospheres. The scaling arguments show that superrotation is favored when the off-equatorial baroclinicity and planetary rotation rates are low. Similarly, superrotation is favored when the convective heating strengthens, which may account for the superrotation seen in extreme global-warming simulations. Third, I use a moist slab-ocean general circulation model to study the impact of a zonally-symmetric continent on the distribution of monsoonal precipitation. I show that adding a hemispheric asymmetry in surface heat capacity is sufficient to cause symmetry breaking in both the spatial and temporal distribution of precipitation. This spatial symmetry breaking can be understood from a large-scale energetic perspective, while the temporal symmetry breaking requires

  15. Regional forecasting with global atmospheric models

    This report was prepared by the Applied Research Corporation (ARC), College Station, Texas, under subcontract to Pacific Northwest Laboratory (PNL) as part of a global climate studies task. The task supports site characterization work required for the selection of a potential high-level nuclear waste repository and is part of the Performance Assessment Scientific Support (PASS) Program at PNL. The work is under the overall direction of the Office of Civilian Radioactive Waste Management (OCRWM), US Department of Energy Headquarters, Washington, DC. The scope of the report is to present the results of the third year's work on the atmospheric modeling part of the global climate studies task. The development testing of computer models and initial results are discussed. The appendices contain several studies that provide supporting information and guidance to the modeling work and further details on computer model development. Complete documentation of the models, including user information, will be prepared under separate reports and manuals

  16. Atmospheric dispersion modeling of radioactive effluents

    In case of a nuclear accident, which could lead to release of radioactive contaminants, fastest countermeasures are needed, relating to sheltering, iodine distribution, evacuation and interdiction of food and water consumption. All these decisions should be based either on estimation of inhaled dose and the dose due to external exposure for public or on the estimation of radioactive concentration in food (which will depend on the radioactive concentration in air and ground deposition). In order to perform any of these calculations of consequences in case of nuclear accident, which leads to release of radioactive contaminants in the atmosphere, we must start with atmospheric dispersion calculations. In the last few years, considerable efforts have been devoted in order to improve computer codes for dispersion in the atmosphere of the radioactive contaminants released in a nuclear accident. The paper presents the model used in computer codes for assessment of nuclear accident consequences and a special attention was paid to the dispersion model used in the Institute for Nuclear Research Pitesti. The values for the used parameters and the results for air and ground concentration are also presented. (authors)

  17. Atmospheric corrosion: statistical validation of models

    In this paper we discuss two different methods for validation of regression models, applied to corrosion data. One of them is based on the correlation coefficient and the other one is the statistical test of lack of fit. Both methods are used here to analyse fitting of bi logarithmic model in order to predict corrosion for very low carbon steel substrates in rural and urban-industrial atmospheres in Uruguay. Results for parameters A and n of the bi logarithmic model are reported here. For this purpose, all repeated values were used instead of using average values as usual. Modelling is carried out using experimental data corresponding to steel substrates under the same initial meteorological conditions ( in fact, they are put in the rack at the same time). Results of correlation coefficient are compared with the lack of it tested at two different signification levels (α=0.01 and α=0.05). Unexpected differences between them are explained and finally, it is possible to conclude, at least in the studied atmospheres, that the bi logarithmic model does not fit properly the experimental data. (Author) 18 refs

  18. Organic chemistry in the atmosphere. [laboratory modeling of Titan atmosphere

    Sagan, C.

    1974-01-01

    The existence of an at least moderately complex organic chemistry on Titan is stipulated based on clear evidence of methane, and at least presumptive evidence of hydrogen in its atmosphere. The ratio of methane to hydrogen is the highest of any atmosphere in the solar system. Irradiation of hydrogen/methane mixtures produces aromatic and aliphatic hydrocarbons. A very reasonable hypothesis assumes that the red cloud cover of Titan is made of organic chemicals. Two-carbon hydrocarbons experimentally produced from irradiated mixtures of methane, ammonia, water, and hydrogen bear out the possible organic chemistry of the Titanian environment.

  19. CIDGA - Coupling of Interior Dynamic models with Global Atmosphere models

    Noack, Lena; Plesa, Ana-Catalina; Breuer, Doris

    2010-05-01

    Atmosphere temperatures and in particular the surface temperatures mostly depend on the solar heat flux and the atmospheric composition. The latter can be influenced by interior processes of the planet, i.e. volcanism that releases greenhouse gases such as H2O, CO2 and methane into the atmosphere and plate tectonics through which atmospheric CO2 is recycled via carbonates into the mantle. An increasing concentration of greenhouse gases in the atmosphere results in an increase of the surface temperature. Changes in the surface temperature on the other hand may influence the cooling behaviour of the planet and hence influence its volcanic activity [Phillips et al., 2001]. This feedback relation between mantle convection and atmosphere is not very well understood, since until now mostly either the interior dynamic of a planet or its atmosphere was investigated separately. 2D or 3D mantle convection models to the authors' knowledge haven't been coupled to the atmosphere so far. We have used the 3D spherical simulation code GAIA [Hüttig et al., 2008] including partial melt production and coupled it with the atmosphere module CIDGA using a gray greenhouse model for varying H2O concentrations. This way, not only the influence of mantle dynamics on the atmosphere can be investigated, but also the recoupling effect, that the surface temperature has on the mantle dynamics. So far, we consider one-plate planets without crustal and thus volatile recycling. Phillips et al. [2001] already investigated the coupling effect of the surface temperature on mantle dynamics by using simple parameterized convection models for Venus. In their model a positive feedback mechanism has been observed, i.e., an increase of the surface temperature leads to an increase of partial melt and hence an increase of atmosphere density and surface temperature. Applying our model to Venus, we show that an increase of surface temperature leads not only to an increase of partial melt in the mantle; it also

  20. Atmospheric transmittance model for photosynthetically active radiation

    Paulescu, Marius; Stefu, Nicoleta; Gravila, Paul; Paulescu, Eugenia; Boata, Remus; Pacurar, Angel; Mares, Oana [Physics Department, West University of Timisoara, V Parvan 4, 300223 Timisoara (Romania); Pop, Nicolina [Department of Physical Foundations of Engineering, Politehnica University of Timisoara, V Parvan 2, 300223 Timisoara (Romania); Calinoiu, Delia [Mechanical Engineering Faculty, Politehnica University of Timisoara, Mihai Viteazu 1, 300222 Timisoara (Romania)

    2013-11-13

    A parametric model of the atmospheric transmittance in the PAR band is presented. The model can be straightforwardly applied for calculating the beam, diffuse and global components of the PAR solar irradiance. The required inputs are: air pressure, ozone, water vapor and nitrogen dioxide column content, Ångström's turbidity coefficient and single scattering albedo. Comparison with other models and ground measured data shows a reasonable level of accuracy for this model, making it suitable for practical applications. From the computational point of view the calculus is condensed into simple algebra which is a noticeable advantage. For users interested in speed-intensive computation of the effective PAR solar irradiance, a PC program based on the parametric equations along with a user guide are available online at http://solar.physics.uvt.ro/srms.

  1. Supermodeling With A Global Atmospheric Model

    Wiegerinck, Wim; Burgers, Willem; Selten, Frank

    2013-04-01

    In weather and climate prediction studies it often turns out to be the case that the multi-model ensemble mean prediction has the best prediction skill scores. One possible explanation is that the major part of the model error is random and is averaged out in the ensemble mean. In the standard multi-model ensemble approach, the models are integrated in time independently and the predicted states are combined a posteriori. Recently an alternative ensemble prediction approach has been proposed in which the models exchange information during the simulation and synchronize on a common solution that is closer to the truth than any of the individual model solutions in the standard multi-model ensemble approach or a weighted average of these. This approach is called the super modeling approach (SUMO). The potential of the SUMO approach has been demonstrated in the context of simple, low-order, chaotic dynamical systems. The information exchange takes the form of linear nudging terms in the dynamical equations that nudge the solution of each model to the solution of all other models in the ensemble. With a suitable choice of the connection strengths the models synchronize on a common solution that is indeed closer to the true system than any of the individual model solutions without nudging. This approach is called connected SUMO. An alternative approach is to integrate a weighted averaged model, weighted SUMO. At each time step all models in the ensemble calculate the tendency, these tendencies are weighted averaged and the state is integrated one time step into the future with this weighted averaged tendency. It was shown that in case the connected SUMO synchronizes perfectly, the connected SUMO follows the weighted averaged trajectory and both approaches yield the same solution. In this study we pioneer both approaches in the context of a global, quasi-geostrophic, three-level atmosphere model that is capable of simulating quite realistically the extra

  2. A new time-dependent analytic model for radiation-induced photocurrent in finite 1D epitaxial diodes.

    Verley, Jason C.; Axness, Carl L.; Hembree, Charles Edward; Keiter, Eric Richard; Kerr, Bert (New Mexico Institute of Mining and Technology, Socorro, NM)

    2012-04-01

    Photocurrent generated by ionizing radiation represents a threat to microelectronics in radiation environments. Circuit simulation tools such as SPICE [1] can be used to analyze these threats, and typically rely on compact models for individual electrical components such as transistors and diodes. Compact models consist of a handful of differential and/or algebraic equations, and are derived by making simplifying assumptions to any of the many semiconductor transport equations. Historically, many photocurrent compact models have suffered from accuracy issues due to the use of qualitative approximation, rather than mathematically correct solutions to the ambipolar diffusion equation. A practical consequence of this inaccuracy is that a given model calibration is trustworthy over only a narrow range of operating conditions. This report describes work to produce improved compact models for photocurrent. Specifically, an analytic model is developed for epitaxial diode structures that have a highly doped subcollector. The analytic model is compared with both numerical TCAD calculations, as well as the compact model described in reference [2]. The new analytic model compares well against TCAD over a wide range of operating conditions, and is shown to be superior to the compact model from reference [2].

  3. Numerical modeling of atmospheric washout processes

    For the washout of particles from the atmosphere by clouds and rain one has to distinguish between processes which work in the first phase of cloud development, when condensation nuclei build up in saturated air (Nucleation Aerosol Scavenging, NAS) and those processes which work at the following cloud development. In the second case particles are taken off by cloud droplets or by falling rain drops via collision (Collision Aerosol Scavenging, CAS). The physics of both processes is described. For the CAS process a numerical model is presented. The report contains a documentation of the mathematical equations and the computer programs (FORTRAN). (KW)

  4. Analysis of software for modeling atmospheric dispersion

    During last few years, a number software packages for microcomputes have appeared with the aim to simulate diffusion of atmospheric pollutants. These codes, simplifying the models used for safety analyses of industrial plants are becoming more useful, and are even used for post-accidental conditions. The report presents for the first time in a critical manner, principal models available up to this date. The problem arises in adapting the models to the demanded post-accidental interventions. In parallel to this action an analysis of performance was performed. It means, identifying the need of forecasting the most appropriate actions to be performed having in mind short available time and lack of information. Because of these difficulties, it is possible to simplify the software, which will not include all the options but could deal with a specific situation. This would enable minimisation of data to be collected on the site

  5. Differences in Water Vapor Radiative Transfer among 1D Models Can Significantly Affect the Inner Edge of the Habitable Zone

    Yang, Jun; Leconte, Jérémy; Wolf, Eric T.; Goldblatt, Colin; Feldl, Nicole; Merlis, Timothy; Wang, Yuwei; Koll, Daniel D. B.; Ding, Feng; Forget, François; Abbot, Dorian S.

    2016-08-01

    An accurate estimate of the inner edge of the habitable zone is critical for determining which exoplanets are potentially habitable and for designing future telescopes to observe them. Here, we explore differences in estimating the inner edge among seven one-dimensional radiative transfer models: two line-by-line codes (SMART and LBLRTM) as well as five band codes (CAM3, CAM4_Wolf, LMDG, SBDART, and AM2) that are currently being used in global climate models. We compare radiative fluxes and spectra in clear-sky conditions around G and M stars, with fixed moist adiabatic profiles for surface temperatures from 250 to 360 K. We find that divergences among the models arise mainly from large uncertainties in water vapor absorption in the window region (10 μm) and in the region between 0.2 and 1.5 μm. Differences in outgoing longwave radiation increase with surface temperature and reach 10–20 W m‑2 differences in shortwave reach up to 60 W m‑2, especially at the surface and in the troposphere, and are larger for an M-dwarf spectrum than a solar spectrum. Differences between the two line-by-line models are significant, although smaller than among the band models. Our results imply that the uncertainty in estimating the insolation threshold of the inner edge (the runaway greenhouse limit) due only to clear-sky radiative transfer is ≈10% of modern Earth’s solar constant (i.e., ≈34 W m‑2 in global mean) among band models and ≈3% between the two line-by-line models. These comparisons show that future work is needed that focuses on improving water vapor absorption coefficients in both shortwave and longwave, as well as on increasing the resolution of stellar spectra in broadband models.

  6. Using 1 -D and 2-D modelling of ground motion for seismic zonation criteria: results for the city of Rome

    A. Caserta; L. Malagnini; A. Rovelli; Marra, F

    1995-01-01

    The geological information collected in the last years by the Istituto Nazionale di Geofisica for the city of Rome is used to construct 1- and 2-D models of the nearsurface structure. These models are the basis for the numerical generation of synthetic accelerograms which can simulate the horizontal ground motion (SH waves) produced in the different areas of the city by a large (M ? 7) potential earthquake 100 km away in Central Apennines. The proposed methodology yields earthquake engineerin...

  7. Status of the solar and infrared radiation submodels in the LLNL 1-D and 2-D chemical-transport models

    The authors have implemented a series of state of the art radiation transport submodels in previously developed one dimensional and two dimensional chemical transport models of the troposphere and stratosphere. These submodels provide the capability of calculating accurate solar and infrared heating rates. They are a firm basis for further radiation submodel development as well as for studying interactions between radiation and model dynamics under varying conditions of clear sky, clouds, and aerosols. 37 refs., 3 figs

  8. Numerical evaluation of inverse modelling methods for 1D and 3D water infiltration experiments in homogeneous soils

    Lassabatère, L.; Yilmaz, D.; Angulo-Jaramillo, R.; Soria Ugalde, J.; Braud, I.; Simunek, J.

    2010-01-01

    International audience Modelling and understanding water fluxes in the vadose zone are important in regards to water management and require appropriate characterization methods of soil hydraulic properties. The presented work studies three common methods for characterization of soil hydraulic properties based on the inverse modelling of Beerkan water infiltration experiments: the CI method for Cumulative Information method and two BEST methods for Beerkan Estimation Soil pedotransfer metho...

  9. Swift observations of the 2006 outburst of the recurrent nova RS Ophiuchi: II. 1D hydrodynamical models of wind driven shocks

    Vaytet, N M H; Bode, M F

    2007-01-01

    Following the early Swift X-ray observations of the latest outburst of the recurrent nova RS Ophiuchi in February 2006 (Paper I), we present new 1D hydrodynamical models of the system which take into account all three phases of the remnant evolution. The models suggest a novel way of modelling the system by treating the outburst as a sudden increase then decrease in wind mass-loss rate and velocity. The differences between this wind model and previous Primakoff-type simulations are described. A more complex structure, even in 1D, is revealed through the presence of both forward and reverse shocks, with a separating contact discontinuity. The effects of radiative cooling are investigated and key outburst parameters such as mass-loss rate, ejecta velocity and mass are varied. The shock velocities as a function of time are compared to the ones derived in Paper I. We show how the manner in which the matter is ejected controls the evolution of the shock and that for a well-cooled remnant, the shock deceleration ra...

  10. The Atmospheric Radionuclide Transport Model (ARTM) - Validation of a long-term atmospheric dispersion model

    Hettrich, Sebastian; Wildermuth, Hans; Strobl, Christopher; Wenig, Mark

    2016-04-01

    In the last couple of years, the Atmospheric Radionuclide Transport Model (ARTM) has been developed by the German Federal Office for Radiation Protection (BfS) and the Society for Plant and Reactor Security (GRS). ARTM is an atmospheric dispersion model for continuous long-term releases of radionuclides into the atmosphere, based on the Lagrangian particle model. This model, developed in the first place as a more realistic replacement for the out-dated Gaussian plume models, is currently being optimised for further scientific purposes to study atmospheric dispersion in short-range scenarios. It includes a diagnostic wind field model, allows for the application of building structures and multiple sources (including linear, 2-and 3-dimensional source geometries), and considers orography and surface roughness. As an output it calculates the activity concentration, dry and wet deposition and can model also the radioactive decay of Rn-222. As such, ARTM requires to undergo an intense validation process. While for short-term and short-range models, which were mainly developed for examining nuclear accidents or explosions, a few measurement data-sets are available for validation, data-sets for validating long-term models are very sparse and the existing ones mostly prove to be not applicable for validation. Here we present a strategy for the validation of long-term Lagrangian particle models based on the work with ARTM. In our validation study, the first part we present is a comprehensive analysis of the model sensitivities on different parameters like e.g. (simulation grid size resolution, starting random number, amount of simulation particles, etc.). This study provides a good estimation for the uncertainties of the simulation results and consequently can be used to generate model outputs comparable to the available measurements data at various distances from the emission source. This comparison between measurement data from selected scenarios and simulation results

  11. Improved Large-Scale Inundation Modelling by 1D-2D Coupling and Consideration of Hydrologic and Hydrodynamic Processes - a Case Study in the Amazon

    Hoch, J. M.; Bierkens, M. F.; Van Beek, R.; Winsemius, H.; Haag, A.

    2015-12-01

    Understanding the dynamics of fluvial floods is paramount to accurate flood hazard and risk modeling. Currently, economic losses due to flooding constitute about one third of all damage resulting from natural hazards. Given future projections of climate change, the anticipated increase in the World's population and the associated implications, sound knowledge of flood hazard and related risk is crucial. Fluvial floods are cross-border phenomena that need to be addressed accordingly. Yet, only few studies model floods at the large-scale which is preferable to tiling the output of small-scale models. Most models cannot realistically model flood wave propagation due to a lack of either detailed channel and floodplain geometry or the absence of hydrologic processes. This study aims to develop a large-scale modeling tool that accounts for both hydrologic and hydrodynamic processes, to find and understand possible sources of errors and improvements and to assess how the added hydrodynamics affect flood wave propagation. Flood wave propagation is simulated by DELFT3D-FM (FM), a hydrodynamic model using a flexible mesh to schematize the study area. It is coupled to PCR-GLOBWB (PCR), a macro-scale hydrological model, that has its own simpler 1D routing scheme (DynRout) which has already been used for global inundation modeling and flood risk assessments (GLOFRIS; Winsemius et al., 2013). A number of model set-ups are compared and benchmarked for the simulation period 1986-1996: (0) PCR with DynRout; (1) using a FM 2D flexible mesh forced with PCR output and (2) as in (1) but discriminating between 1D channels and 2D floodplains, and, for comparison, (3) and (4) the same set-ups as (1) and (2) but forced with observed GRDC discharge values. Outputs are subsequently validated against observed GRDC data at Óbidos and flood extent maps from the Dartmouth Flood Observatory. The present research constitutes a first step into a globally applicable approach to fully couple

  12. Effects of a space modulation on the behavior of a 1D alternating Heisenberg spin-1/2 model.

    Mahdavifar, Saeed; Abouie, Jahanfar

    2011-06-22

    The effects of a magnetic field (h) and a space modulation (δ) on the magnetic properties of a one-dimensional antiferromagnetic-ferromagnetic Heisenberg spin-1/2 model have been studied by means of numerical exact diagonalization of finite size systems, the nonlinear σ model, and a bosonization approach. The space modulation is considered on the antiferromagnetic couplings. At δ = 0, the model is mapped to a gapless Lüttinger liquid phase by increasing the magnetic field. However, the space modulation induces a new gap in the spectrum of the system and the system experiences different quantum phases which are separated by four critical fields. By opening the new gap, a magnetization plateau appears at ½M(sat). The effects of the space modulation are reflected in the emergence of a plateau in other physical functions such as the F-dimer and the bond-dimer order parameters, and the pair-wise entanglement. PMID:21613724

  13. Butyrylcholinesterase is Associated with β-Amyloid Plaques in the Transgenic APPSWE/PSEN1dE9 Mouse Model of Alzheimer Disease

    Darvesh, Sultan; Cash, Meghan K.; Reid, G. Andrew; Martin, Earl; Mitnitski, Arnold; Geula, Changiz

    2012-01-01

    Histochemical analysis of Alzheimer disease (AD) brain tissues indicates that butyrylcholinesterase (BuChE) is present in β-amyloid (Aβ) plaques. The role of BuChE in AD pathology is unknown but an animal model developing similar BuChE-associated Aβ plaques could provide insights. The APPSWE/PSEN1dE9 mouse (ADTg), which develops Aβ plaques, was examined to determine if BuChE associates with these plaques, as in AD. We found that in mature ADTg mice, BuChE activity associated with Aβ plaques. ...

  14. Pulse wave propagation in a model human arterial network: Assessment of 1-D visco-elastic simulations against in vitro measurements

    Alastruey, Jordi; Khir, Ashraf W.; Matthys, Koen S.; Segers, Patrick; Sherwin, Spencer J.; Verdonck, Pascal R.; Parker, Kim H.; Peiro, Joaquim

    2011-01-01

    The accuracy of the nonlinear one-dimensional (1-D) equations of pressure and flow wave propagation in Voigt-type visco-elastic arteries was tested against measurements in a well-defined experimental 1:1 replica of the 37 largest conduit arteries in the human systemic circulation. The parameters required by the numerical algorithm were directly measured in the in vitro setup and no data fitting was involved. The inclusion of wall visco-elasticity in the numerical model reduced the underdamped...

  15. Analytical solution to the 1D Lemaitre's isotropic damage model and plane stress projected implicit integration procedure

    Andriollo, Tito; Thorborg, Jesper; Hattel, Jesper Henri

    2016-01-01

    In the present paper, for the first time in literature an exact analytical solution to Lemaitre's isotropic damage model is developed for the special case of uniaxial tensile testing. This is achieved by taking advantage of a convenient formulation of the isotropic hardening function, which allows...... optimization, as all issues associated with classical numerical solution procedures of the constitutive equations are eliminated. In addition, an implicit implementation of the plane stress projected version of Lemaitre's model is discussed, showing that the resulting algebraic system can be reduced to a...

  16. Integrating models to simulate emergent behaviour: effects of organic matter on soil hydraulics in the ICZ-1D soil-vegetation model

    Valstar, Johan; Rowe, Ed; Konstantina, Moirogiorgou; Giannakis, Giorgos; Nikolaidis, Nikolaos

    2014-05-01

    explore the complex interactions involved in soil development and change. We were unable to identify appropriately-detailed existing models for plant productivity and for the dynamics of soil aggregation and porosity, and so developed the PROSUM and CAST models, respectively, to simulate these subsystems. Moreover, we applied the BRNS generator to obtain a chemical equilibrium model. These were combined with HYDRUS-1D (water and solute transport), a weathering model (derived from the SAFE model) and a simple bioturbation model. The model includes several feedbacks, such as the effect of soil organic matter on water retention and hydraulic conductivity. We encountered several important challenges when building the integrated model. First, a mechanism was developed that initiates the execution of a single time step for an individual sub-model and accounts for the relevant mass transfers between sub-models. This allows for different and sometimes variable time step duration in the submodels. Secondly, we removed duplicated processes and identified and included relevant solute production terms that had been neglected. The model is being tested against datasets obtained from several Soil Critical Zone Observatories in Europe. This contribution focuses on the design strategy for the model.

  17. PCBs in the Arctic atmosphere: determining important driving forces using a global atmospheric transport model

    Friedman, Carey L.; Selin, Noelle E

    2016-01-01

    We present a spatially and temporally resolved global atmospheric polychlorinated biphenyl (PCB) model, driven by meteorological data, that is skilled at simulating mean atmospheric PCB concentrations and seasonal cycles in the Northern Hemisphere midlatitudes and mean Arctic concentrations. However, the model does not capture the observed Arctic summer maximum in atmospheric PCBs. We use the model to estimate global budgets for seven PCB congeners, and we demonstrate that c...

  18. Using 1 -D and 2-D modelling of ground motion for seismic zonation criteria: results for the city of Rome

    A. Caserta

    1995-06-01

    Full Text Available The geological information collected in the last years by the Istituto Nazionale di Geofisica for the city of Rome is used to construct 1- and 2-D models of the nearsurface structure. These models are the basis for the numerical generation of synthetic accelerograms which can simulate the horizontal ground motion (SH waves produced in the different areas of the city by a large (M ? 7 potential earthquake 100 km away in Central Apennines. The proposed methodology yields earthquake engineering parameters (peak ground acceleration and velocity, Arias intensity, energy flux, response spectra whose spatial variations are consistent with the damage distribution caused by the strongest earthquakes felt in Rome during its long history. Based on the macroseismic inforination and the results of the numerical simulations, general criteria for seismic zonation of the city of Rome are proposed.

  19. Estimation of the Origin of a Speaker's and singer's Formant Cluster using an Optimization of 1D Vocal Tract Model

    Laukkanen, A. M.; Radolf, Vojtěch; Horáček, Jaromír; Leino, T.

    Madrid : EUIT Telecomunicación, 2009, s. 1-4. ISBN 84-95227-64-9. [Advanced Voice Function Assessment International Workshop /3./. Madrid (ES), 18.05.2009-20.05.2009] R&D Projects: GA ČR GA101/08/1155 Institutional research plan: CEZ:AV0Z20760514 Keywords : vocal tract acoustics * mathematical modelling * operatic singing * good speaking voice Subject RIV: BI - Acoustics

  20. Minimum 1D velocity models in Central and Southern Italy: a contribution to better constrain hypocentral determinations

    A. Frepoli

    1997-06-01

    Full Text Available We computed one-dimensional ( I D velocity models and station corrections for Centrai and Southern Italy, in- verting re-picked P-wave alTival times recorded by the Istituto Nazionale di Geofisica seismic network. The re-picked data yield resolved P-wave velocity results and proved to be more suited than bulletin data for de- tailed tomographic studies. Using the improved velocity models, we relocated the most significant earthquakes which occurt.ed in the Apennines in the past 7 years, achieving constrained hypocentral determinations for events within most of the Apenninic belt. The interpretation of the obtained lD velocity models allows us to infer interesting features on the deep structure of the Apennines. Smooth velocity gradients with depth and low P-wave velocities are ob,'ierved beneath the Apennines. We believe that our results are effective to constrain hypocentral locations in Italy and may represent a first step towards more detailed seismotectonic analyses.

  1. Comparing the Degree of Land-Atmosphere Interaction in Four Atmospheric General Circulation Models

    Koster, Randal D.; Dirmeyer, Paul A.; Hahmann, Andrea N.; Ijpelaar, Ruben; Tyahla, Lori; Cox, Peter; Suarez, Max J.; Houser, Paul R. (Technical Monitor)

    2001-01-01

    Land-atmosphere feedback, by which (for example) precipitation-induced moisture anomalies at the land surface affect the overlying atmosphere and thereby the subsequent generation of precipitation, has been examined and quantified with many atmospheric general circulation models (AGCMs). Generally missing from such studies, however, is an indication of the extent to which the simulated feedback strength is model dependent. Four modeling groups have recently performed a highly controlled numerical experiment that allows an objective inter-model comparison of land-atmosphere feedback strength. The experiment essentially consists of an ensemble of simulations in which each member simulation artificially maintains the same time series of surface prognostic variables. Differences in atmospheric behavior between the ensemble members then indicates the degree to which the state of the land surface controls atmospheric processes in that model. A comparison of the four sets of experimental results shows that feedback strength does indeed vary significantly between the AGCMs.

  2. Non-trivial exponents in the zero temperature dynamics of the 1D Ising and Potts models

    Derrida, B.; Bray, A. J.; Godrèche, C.

    1994-01-01

    We consider the Glauber dynamics of the $ q $-state Potts model in one dimension at zero temperature. Starting with a random initial configuration, we measure the density $ r_t $ of spins which have never flipped from the beginning of the simulation until time $ t. $ We find that for large $ t, $ the density $ r_t $ has a power law decay $ \\left(r_t \\sim t^{-\\theta} \\right) $ where the exponent $ \\theta $ varies with $ q. $ Our simulations lead to $ \\theta \\simeq .37 $ for $ q=2, $ $ \\theta \\...

  3. High Resolution Transmission Spectroscopy as a Diagnostic for Jovian Exoplanet Atmospheres: Constraints from Theoretical Models

    Kempton, Eliza M -R; Heng, Kevin

    2014-01-01

    We present high resolution transmission spectra of giant planet atmospheres from a coupled 3-D atmospheric dynamics and transmission spectrum model that includes Doppler shifts which arise from winds and planetary motion. We model jovian planets covering more than two orders of magnitude in incident flux, corresponding to planets with 0.9 to 55 day orbital periods around solar-type stars. The results of our 3-D dynamical models reveal certain aspects of high resolution transmission spectra that are not present in simple 1-D models. We find that the hottest planets experience strong substellar to anti-stellar (SSAS) winds, resulting in transmission spectra with net blue shifts of up to 3 km s$^{-1}$, whereas less irradiated planets show almost no net Doppler shifts. Compared to 1-D models, peak line strengths are significantly reduced for the hottest atmospheres owing to Doppler broadening from a combination of rotation (which is faster for close-in planets under the assumption of tidal locking) and atmospheri...

  4. Simulation of Marine Boundary Layer characteristics using a 1-D PBL model over the Bay of Bengal during BOBMEX-99

    N V Sam; U C Mohanty; A N V Satyanarayana

    2003-06-01

    The characteristic features of the marine boundary layer (MBL) over the Bay of Bengal during the southwest monsoon and the factors influencing it are investigated. The Bay of Bengal and Monsoon Experiment (BOBMEX) carried out during July-August 1999 is the first observational experiment under the Indian Climate Research Programme (ICRP). A very high-resolution data in the vertical was obtained during this experiment, which was used to study the MBL characteristics off the east coast of India in the north and south Bay of Bengal. Spells of active and suppressed convection over the Bay were observed, of which, three representative convective episodes were considered for the study. For this purpose a one-dimensional multi-level PBL model with a TKE- closure scheme was used. The soundings, viz., the vertical profiles of temperature, humidity, zonal and meridional component of wind, obtained onboard ORV Sagar Kanya and from coastal stations along the east coast are used for the study. The temporal evolution of turbulent kinetic energy, marine boundary layer height (MBLH), sensible and latent heat fluxes and drag coefficient of momentum are simulated for different epochs of monsoon and monsoon depressions during BOBMEX-99.The model also generates the vertical profiles of potential temperature, specific humidity, zonal and meridional wind. These simulated values compared reasonably well with the observations available from BOBMEX.

  5. Evaluation of Bulk Charging in Geostationary Transfer Orbit and Earth Escape Trajectories Using the Numit 1-D Charging Model

    Minow, Joseph I.; Coffey, Victoria N.; Parker, Linda N.; Blackwell, William C., Jr.; Jun, Insoo; Garrett, Henry B.

    2007-01-01

    The NUMIT 1-dimensional bulk charging model is used as a screening to ol for evaluating time-dependent bulk internal or deep dielectric) ch arging of dielectrics exposed to penetrating electron environments. T he code is modified to accept time dependent electron flux time serie s along satellite orbits for the electron environment inputs instead of using the static electron flux environment input originally used b y the code and widely adopted in bulk charging models. Application of the screening technique ts demonstrated for three cases of spacecraf t exposure within the Earth's radiation belts including a geostationa ry transfer orbit and an Earth-Moon transit trajectory for a range of orbit inclinations. Electric fields and charge densities are compute d for dielectric materials with varying electrical properties exposed to relativistic electron environments along the orbits. Our objectiv e is to demonstrate a preliminary application of the time-dependent e nvironments input to the NUMIT code for evaluating charging risks to exposed dielectrics used on spacecraft when exposed to the Earth's ra diation belts. The results demonstrate that the NUMIT electric field values in GTO orbits with multiple encounters with the Earth's radiat ion belts are consistent with previous studies of charging in GTO orb its and that potential threat conditions for electrostatic discharge exist on lunar transit trajectories depending on the electrical proper ties of the materials exposed to the radiation environment.

  6. Transient runoff-runon model for a 1-D slope with random infiltrability: flow statistics and connectivity

    Harel, Marie-Alice; Mouche, Emmanuel

    2015-04-01

    Despite the recent research focused on runoff pattern connectivity in hydrology, there is a surprising lack of theoretical knowledge regarding hillslope runoff generation and dynamics during a rainfall event. The transient problem is especially unaddressed. In this paper we propose a model based on queueing theory formalism for the infiltration-excess overland flow generation on soils with random infiltration properties. The influence of rainfall intensity and duration on runoff dynamics and connectivity is studied thanks to this model, numerical simulation and available steady-state results. We limit our study to a rainfall intensity that is a rectangular function of time. Exact solutions for the case of spatially random exponential distributions of soil infiltrability and rainfall intensity are developed. Simulations validate these analytical results and allow for the study the rising and recession limbs of the hydrograph for different rainfall characteristics. The case of a deterministic uniform rainfall rate and different infiltrability distributions is also discussed in light of runoff connectivity. We show that the connectivity framework contributes to a better understanding and prediction of runoff pattern formation and evolution with time. A fragmented overland flow is shown to have shorter charge and discharge periods after the onset and offset of rainfall compared to well connected runoff fields. These results demonstrate that the transient regime characteristics are linked with connectivity parameters, rainstorm properties and scale issues.

  7. A time series generalized functional model based method for vibration-based damage precise localization in structures consisting of 1D, 2D, and 3D elements

    Sakaris, C. S.; Sakellariou, J. S.; Fassois, S. D.

    2016-06-01

    This study focuses on the problem of vibration-based damage precise localization via data-based, time series type, methods for structures consisting of 1D, 2D, or 3D elements. A Generalized Functional Model Based method is postulated based on an expanded Vector-dependent Functionally Pooled ARX (VFP-ARX) model form, capable of accounting for an arbitrary structural topology. The FP model's operating parameter vector elements are properly constrained to reflect any given topology. Damage localization is based on operating parameter vector estimation within the specified topology, so that the location estimate and its uncertainty bounds are statistically optimal. The method's effectiveness is experimentally demonstrated through damage precise localization on a laboratory spatial truss structure using various damage scenarios and a single pair of random excitation - vibration response signals in a low and limited frequency bandwidth.

  8. Aerial Radioactivity Monitoring Using Atmospheric Dispersion Model

    Since North Korea announced the underground nuclear test on last October 9th, 2006, many countries including South Korea have worried about the atmospheric dispersion and pollution of radioactive material by nuclear test. To verify the existence of nuclear test by detecting radioactive materials such as xenon and krypton at the early stage, to locate the position of test site, and to chase the trajectory of radioactivity have been heavily issued. And radioactivity detection and radiation monitoring technology using an aircraft have been recently examined by an authority concerned in South Korea. Although various techniques of aerial radioactivity monitoring are developed and operated in the world such as United States of America, Japan, Germany, etc., the relevant technical development or research is wholly lacking in our country. In this study, we performed some case studies on North Korea's nuclear test and accidental releases from nuclear power plant (NPP) using HYSPLIT (HYbrid Single Particle Lagrangian Integrated Trajectory) model developed by National Oceanic and Atmospheric Administration (NOAA) of U.S. Department of Commerce. We also investigated a feasibility of HYSPLIT to the aerial radioactivity monitoring system in terms of deciding potential measuring location and time

  9. Modeling seismic wave propagation and amplification in 1D/2D/3D linear and nonlinear unbounded media

    Semblat, Jean-François

    2011-01-01

    To analyze seismic wave propagation in geological structures, it is possible to consider various numerical approaches: the finite difference method, the spectral element method, the boundary element method, the finite element method, the finite volume method, etc. All these methods have various advantages and drawbacks. The amplification of seismic waves in surface soil layers is mainly due to the velocity contrast between these layers and, possibly, to topographic effects around crests and hills. The influence of the geometry of alluvial basins on the amplification process is also know to be large. Nevertheless, strong heterogeneities and complex geometries are not easy to take into account with all numerical methods. 2D/3D models are needed in many situations and the efficiency/accuracy of the numerical methods in such cases is in question. Furthermore, the radiation conditions at infinity are not easy to handle with finite differences or finite/spectral elements whereas it is explicitely accounted in the B...

  10. 1-d model for propagation and absorption of h.f. waves near ion cyclotron resonances in large tokamak plasmas

    The behaviour of h.f. waves near ion cyclotron and ion ion hybrid resonances in a tokamak is investigated by means of an one-dimensional finite element code. Our model takes into account: - strength and orientation of the poloidal component of the static magnetic field, - finite larmor radius corrections to the dielectric tensor, - ion cyclotron damping at the fundamental and first harmonic resonance, - electron Landau damping to lowest order msub(e)/msub(i). We assume that an incoming fast wave approaches the singular layer from the low or high field side making an arbitrary angle relative to the local magnetic flux surface and to the resonance layer. These initial conditions might be provided by ray tracing from the antenna. Then we calculate the electromagnetic wavefield and the power fluxes of the transmitted or reflected fast and slow waves as well as the power absorbed by ions and electrons. (orig.)