Estimation of 85Kr dispersion from the spent nuclear fuel reprocessing plant in Rokkasho, Japan, using an atmospheric dispersion model.

PubMed

Abe, K; Iyogi, T; Kawabata, H; Chiang, J H; Suwa, H; Hisamatsu, S

2015-11-01

The spent nuclear fuel reprocessing plant of Japan Nuclear Fuel Limited (JNFL) located in Rokkasho, Japan, discharged small amounts of (85)Kr into the atmosphere during final tests of the plant with actual spent fuel from 31 March 2006 to October 2008. During this period, the gamma-ray dose rates due to discharged (85)Kr were higher than the background rates measured at the Institute for Environmental Sciences and at seven monitoring stations of the Aomori prefectural government and JNFL. The dispersion of (85)Kr was simulated by means of the fifth-generation Penn State/NCAR Mesoscale Model and the CG-MATHEW/ADPIC models (ver. 5.0) with a vertical terrain-following height coordinate. Although the simulated gamma-ray dose rates due to discharged (85)Kr agreed fairly well with measured rates, the agreement between the estimated monthly mean (85)Kr concentrations and the observed concentrations was poor. Improvement of the vertical flow of air may lead to better estimation of (85)Kr dispersion. © The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.

Dispersion in Fractures with Ramified Dissolution Patterns

NASA Astrophysics Data System (ADS)

Xu, Le; Marks, Benjy; Toussaint, Renaud; Flekkøy, Eirik G.; Måløy, Knut J.

2018-04-01

The injection of a reactive fluid into an open fracture may modify the fracture surface locally and create a ramified structure around the injection point. This structure will have a significant impact on the dispersion of the injected fluid due to increased permeability, which will introduce large velocity fluctuations into the fluid. Here, we have injected a fluorescent tracer fluid into a transparent artificial fracture with such a ramified structure. The transparency of the model makes it possible to follow the detailed dispersion of the tracer concentration. The experiments have been compared to two dimensional (2D) computer simulations which include both convective motion and molecular diffusion. A comparison was also performed between the dispersion from an initially ramified dissolution structure and the dispersion from an initially circular region. A significant difference was seen both at small and large length scales. At large length scales, the persistence of the anisotropy of the concentration distribution far from the ramified structure is discussed with reference to some theoretical considerations and comparison with simulations.

Effect of numerical dispersion as a source of structural noise in the calibration of a highly parameterized saltwater intrusion model

USGS Publications Warehouse

Langevin, Christian D.; Hughes, Joseph D.

2010-01-01

A model with a small amount of numerical dispersion was used to represent saltwater 7 intrusion in a homogeneous aquifer for a 10-year historical calibration period with one 8 groundwater withdrawal location followed by a 10-year prediction period with two groundwater 9 withdrawal locations. Time-varying groundwater concentrations at arbitrary locations in this low-10 dispersion model were then used as observations to calibrate a model with a greater amount of 11 numerical dispersion. The low-dispersion model was solved using a Total Variation Diminishing 12 numerical scheme; an implicit finite difference scheme with upstream weighting was used for 13 the calibration simulations. Calibration focused on estimating a three-dimensional hydraulic 14 conductivity field that was parameterized using a regular grid of pilot points in each layer and a 15 smoothness constraint. Other model parameters (dispersivity, porosity, recharge, etc.) were 16 fixed at the known values. The discrepancy between observed and simulated concentrations 17 (due solely to numerical dispersion) was reduced by adjusting hydraulic conductivity through the 18 calibration process. Within the transition zone, hydraulic conductivity tended to be lower than 19 the true value for the calibration runs tested. The calibration process introduced lower hydraulic 20 conductivity values to compensate for numerical dispersion and improve the match between 21 observed and simulated concentration breakthrough curves at monitoring locations. 22 Concentrations were underpredicted at both groundwater withdrawal locations during the 10-23 year prediction period.

Irradiation defect dispersions and effective dislocation mobility in strained ferritic grains: A statistical analysis based on 3D dislocation dynamics simulations

NASA Astrophysics Data System (ADS)

Li, Y.; Robertson, C.

2018-06-01

The influence of irradiation defect dispersions on plastic strain spreading is investigated by means of three-dimensional dislocation dynamics (DD) simulations, accounting for thermally activated slip and cross-slip mechanisms in Fe-2.5%Cr grains. The defect-induced evolutions of the effective screw dislocation mobility are evaluated by means of statistical comparisons, for various defect number density and defect size cases. Each comparison is systematically associated with a quantitative Defect-Induced Apparent Straining Temperature shift (or «ΔDIAT»), calculated without any adjustable parameters. In the investigated cases, the ΔDIAT level associated with a given defect dispersion closely replicates the measured ductile to brittle transition temperature shift (ΔDBTT) due to the same, actual defect dispersion. The results are further analyzed in terms of dislocation-based plasticity mechanisms and their possible relations with the dose-dependent changes of the ductile to brittle transition temperature.

Effect of gravitational sedimentation on simulated aerosol dispersion in the human acinus.

PubMed

Darquenne, Chantal; Prisk, G Kim

2003-04-01

We studied the effect of gravitational sedimentation on the dispersion of 0.5 and 1 micrometer-diameter particle boluses within a two-dimensional symmetric six-generation model of the human acinus. Boluses were introduced at the beginning of a 2-s inspiration immediately followed by a 4-s expiration, in normal gravity (1 G) and in the absence of gravity (0 G). The flow corresponded to a flow rate at the mouth of 500 ml/s. In 0 G, simulated dispersion (Hsim) was 16 ml for both particle sizes. In 1 G, Hsim was 71 and 242 ml for 0.5 and 1 micrometer-diameter particles, respectively, showing the effect of gravitational sedimentation. The difference between experimental data (J. Appl. Physiol. 86 (1999) 1402) and simulations was independent of particle size. This suggests that the residual dispersion was independent of the intrinsic properties of the particles and was more likely due to other mechanisms such as ventilation inhomogeneities, cardiogenic oscillations and alveolar wall motion. c2003 Elsevier Science Ltd. All rights reserved.

Effect of gravitational sedimentation on simulated aerosol dispersion in the human acinus

NASA Technical Reports Server (NTRS)

Darquenne, Chantal; Prisk, G. Kim

2003-01-01

We studied the effect of gravitational sedimentation on the dispersion of 0.5 and 1 micrometer-diameter particle boluses within a two-dimensional symmetric six-generation model of the human acinus. Boluses were introduced at the beginning of a 2-s inspiration immediately followed by a 4-s expiration, in normal gravity (1 G) and in the absence of gravity (0 G). The flow corresponded to a flow rate at the mouth of 500 ml/s. In 0 G, simulated dispersion (Hsim) was 16 ml for both particle sizes. In 1 G, Hsim was 71 and 242 ml for 0.5 and 1 micrometer-diameter particles, respectively, showing the effect of gravitational sedimentation. The difference between experimental data (J. Appl. Physiol. 86 (1999) 1402) and simulations was independent of particle size. This suggests that the residual dispersion was independent of the intrinsic properties of the particles and was more likely due to other mechanisms such as ventilation inhomogeneities, cardiogenic oscillations and alveolar wall motion. c2003 Elsevier Science Ltd. All rights reserved.

Lagrangian and Eulerian statistics obtained from direct numerical simulations of homogeneous turbulence

NASA Technical Reports Server (NTRS)

Squires, Kyle D.; Eaton, John K.

1991-01-01

Direct numerical simulation is used to study dispersion in decaying isotropic turbulence and homogeneous shear flow. Both Lagrangian and Eulerian data are presented allowing direct comparison, but at fairly low Reynolds number. The quantities presented include properties of the dispersion tensor, isoprobability contours of particle displacement, Lagrangian and Eulerian velocity autocorrelations and time scale ratios, and the eddy diffusivity tensor. The Lagrangian time microscale is found to be consistently larger than the Eulerian microscale, presumably due to the advection of the small scales by the large scales in the Eulerian reference frame.

Oak habitat recovery on California's largest islands: Scenarios for the role of corvid seed dispersal

USGS Publications Warehouse

Pesendorfer, Mario B.; Baker, Christopher M.; Stringer, Martin; McDonald-Madden, Eve; Bode, Michael; McEachern, A. Kathryn; Morrison, Scott A.; Sillett, T. Scott

2018-01-01

Seed dispersal by birds is central to the passive restoration of many tree communities. Reintroduction of extinct seed dispersers can therefore restore degraded forests and woodlands. To test this, we constructed a spatially explicit simulation model, parameterized with field data, to consider the effect of different seed dispersal scenarios on the extent of oak populations. We applied the model to two islands in California's Channel Islands National Park (USA), one of which has lost a key seed disperser.We used an ensemble modelling approach to simulate island scrub oak (Quercus pacifica) demography. The model was developed and trained to recreate known population changes over a 20-year period on 250-km2 Santa Cruz Island, and incorporated acorn dispersal by island scrub-jays (Aphelocoma insularis), deer mice (Peromyscus maniculatus) and gravity, as well as seed predation. We applied the trained model to 215-km2 Santa Rosa Island to examine how reintroducing island scrub-jays would affect the rate and pattern of oak population expansion. Oak habitat on Santa Rosa Island has been greatly reduced from its historical extent due to past grazing by introduced ungulates, the last of which were removed by 2011.Our simulation model predicts that a seed dispersal scenario including island scrub-jays would increase the extent of the island scrub oak population on Santa Rosa Island by 281% over 100 years, and by 544% over 200 years. Scenarios without jays would result in little expansion. Simulated long-distance seed dispersal by jays also facilitates establishment of discontinuous patches of oaks, and increases their elevational distribution.Synthesis and applications. Scenario planning provides powerful decision support for conservation managers. We used ensemble modelling of plant demographic and seed dispersal processes to investigate whether the reintroduction of seed dispersers could provide cost-effective means of achieving broader ecosystem restoration goals on California's second-largest island. The simulation model, extensively parameterized with field data, suggests that re-establishing the mutualism with seed-hoarding jays would accelerate the expansion of island scrub oak, which could benefit myriad species of conservation concern.

Particle-In-Cell simulations of electron beam microbunching instability in three dimensions

NASA Astrophysics Data System (ADS)

Huang, Chengkun; Zeng, Y.; Meyers, M. D.; Yi, S.; Albright, B. J.; Kwan, T. J. T.

2013-10-01

Microbunching instability due to Coherent Synchrotron Radiation (CSR) in a magnetic chicane is one of the major effects that can degrade the electron beam quality in an X-ray Free Electron Laser. Self-consistent simulation using the Particle-In-Cell (PIC) method for the CSR fields of the beam and their effects on beam dynamics have been elusive due to the excessive dispersion error on the grid. We have implemented a high-order finite-volume PIC scheme that models the propagation of the CSR fields accurately. This new scheme is characterized and optimized through a detailed dispersion analysis. The CSR fields from our improved PIC calculation are compared to the extended CSR numerical model based on the Lienard-Wiechert formula in 2D/3D. We also conduct beam dynamics simulation of the microbunching instability using our new PIC capability. Detailed self-consistent PIC simulations of the CSR fields and beam dynamics will be presented and discussed. Work supported by the U.S. Department of Energy through the LDRD program at Los Alamos National Laboratory.

On the numerical dispersion of electromagnetic particle-in-cell code: Finite grid instability

NASA Astrophysics Data System (ADS)

Meyers, M. D.; Huang, C.-K.; Zeng, Y.; Yi, S. A.; Albright, B. J.

2015-09-01

The Particle-In-Cell (PIC) method is widely used in relativistic particle beam and laser plasma modeling. However, the PIC method exhibits numerical instabilities that can render unphysical simulation results or even destroy the simulation. For electromagnetic relativistic beam and plasma modeling, the most relevant numerical instabilities are the finite grid instability and the numerical Cherenkov instability. We review the numerical dispersion relation of the Electromagnetic PIC model. We rigorously derive the faithful 3-D numerical dispersion relation of the PIC model, for a simple, direct current deposition scheme, which does not conserve electric charge exactly. We then specialize to the Yee FDTD scheme. In particular, we clarify the presence of alias modes in an eigenmode analysis of the PIC model, which combines both discrete and continuous variables. The manner in which the PIC model updates and samples the fields and distribution function, together with the temporal and spatial phase factors from solving Maxwell's equations on the Yee grid with the leapfrog scheme, is explicitly accounted for. Numerical solutions to the electrostatic-like modes in the 1-D dispersion relation for a cold drifting plasma are obtained for parameters of interest. In the succeeding analysis, we investigate how the finite grid instability arises from the interaction of the numerical modes admitted in the system and their aliases. The most significant interaction is due critically to the correct representation of the operators in the dispersion relation. We obtain a simple analytic expression for the peak growth rate due to this interaction, which is then verified by simulation. We demonstrate that our analysis is readily extendable to charge conserving models.

Air-mediated pollen flow from genetically modified to conventional crops.

PubMed

Kuparinen, Anna; Schurr, Frank; Tackenberg, Oliver; O'Hara, Robert B

2007-03-01

Tools for estimating pollen dispersal and the resulting gene flow are necessary to assess the risk of gene flow from genetically modified (GM) to conventional fields, and to quantify the effectiveness of measures that may prevent such gene flow. A mechanistic simulation model is presented and used to simulate pollen dispersal by wind in different agricultural scenarios over realistic pollination periods. The relative importance of landscape-related variables such as isolation distance, topography, spatial configuration of the fields, GM field size and barrier, and environmental variation are examined in order to find ways to minimize gene flow and to detect possible risk factors. The simulations demonstrated a large variation in pollen dispersal and in the predicted amount of contamination between different pollination periods. This was largely due to variation in vertical wind. As this variation in wind conditions is difficult to control through management measures, it should be carefully considered when estimating the risk of gene flow from GM crops. On average, the predicted level of gene flow decreased with increasing isolation distance and with increasing depth of the conventional field, and increased with increasing GM field size. Therefore, at a national scale and over the long term these landscape properties should be accounted for when setting regulations for controlling gene flow. However, at the level of an individual field the level of gene flow may be dominated by uncontrollable variation. Due to the sensitivity of pollen dispersal to the wind, we conclude that gene flow cannot be summarized only by the mean contamination; information about the frequency of extreme events should also be considered. The modeling approach described in this paper offers a way to predict and compare pollen dispersal and gene flow in varying environmental conditions, and to assess the effectiveness of different management measures.

On the numerical dispersion of electromagnetic particle-in-cell code: Finite grid instability

DOE Office of Scientific and Technical Information (OSTI.GOV)

Meyers, M.D., E-mail: mdmeyers@physics.ucla.edu; Department of Physics and Astronomy, University of California Los Angeles, Los Angeles, CA 90095; Huang, C.-K., E-mail: huangck@lanl.gov

The Particle-In-Cell (PIC) method is widely used in relativistic particle beam and laser plasma modeling. However, the PIC method exhibits numerical instabilities that can render unphysical simulation results or even destroy the simulation. For electromagnetic relativistic beam and plasma modeling, the most relevant numerical instabilities are the finite grid instability and the numerical Cherenkov instability. We review the numerical dispersion relation of the Electromagnetic PIC model. We rigorously derive the faithful 3-D numerical dispersion relation of the PIC model, for a simple, direct current deposition scheme, which does not conserve electric charge exactly. We then specialize to the Yee FDTDmore » scheme. In particular, we clarify the presence of alias modes in an eigenmode analysis of the PIC model, which combines both discrete and continuous variables. The manner in which the PIC model updates and samples the fields and distribution function, together with the temporal and spatial phase factors from solving Maxwell's equations on the Yee grid with the leapfrog scheme, is explicitly accounted for. Numerical solutions to the electrostatic-like modes in the 1-D dispersion relation for a cold drifting plasma are obtained for parameters of interest. In the succeeding analysis, we investigate how the finite grid instability arises from the interaction of the numerical modes admitted in the system and their aliases. The most significant interaction is due critically to the correct representation of the operators in the dispersion relation. We obtain a simple analytic expression for the peak growth rate due to this interaction, which is then verified by simulation. We demonstrate that our analysis is readily extendable to charge conserving models.« less

CFD model simulation of LPG dispersion in urban areas

NASA Astrophysics Data System (ADS)

Pontiggia, Marco; Landucci, Gabriele; Busini, Valentina; Derudi, Marco; Alba, Mario; Scaioni, Marco; Bonvicini, Sarah; Cozzani, Valerio; Rota, Renato

2011-08-01

There is an increasing concern related to the releases of industrial hazardous materials (either toxic or flammable) due to terrorist attacks or accidental events in congested industrial or urban areas. In particular, a reliable estimation of the hazardous cloud footprint as a function of time is required to assist emergency response decision and planning as a primary element of any Decision Support System. Among the various hazardous materials, the hazard due to the road and rail transportation of liquefied petroleum gas (LPG) is well known since large quantities of LPG are commercialized and the rail or road transportation patterns are often close to downtown areas. Since it is well known that the widely-used dispersion models do not account for the effects of any obstacle like buildings, tanks, railcars, or trees, in this paper a CFD model has been applied to simulate the reported consequences of a recent major accident involving an LPG railcar rupture in a congested urban area (Viareggio town, in Italy), showing both the large influence of the obstacles on LPG dispersion as well as the potentials of CFD models to foresee such an influence.

Dispersive waves induced by self-defocusing temporal solitons in a beta-barium-borate crystal.

PubMed

Zhou, Binbin; Bache, Morten

2015-09-15

We experimentally observe dispersive waves in the anomalous dispersion regime of a beta-barium-borate (BBO) crystal, induced by a self-defocusing few-cycle temporal soliton. Together the soliton and dispersive waves form an energetic octave-spanning supercontinuum. The soliton was excited in the normal dispersion regime of BBO through a negative cascaded quadratic nonlinearity. Using pump wavelengths from 1.24 to 1.4 μm, dispersive waves are found from 1.9 to 2.2 μm, agreeing well with calculated resonant phase-matching wavelengths due to degenerate four-wave mixing to the soliton. We also observe resonant radiation from nondegenerate four-wave mixing between the soliton and a probe wave, which was formed by leaking part of the pump spectrum into the anomalous dispersion regime. We confirm the experimental results through simulations.

A guided wave dispersion compensation method based on compressed sensing

NASA Astrophysics Data System (ADS)

Xu, Cai-bin; Yang, Zhi-bo; Chen, Xue-feng; Tian, Shao-hua; Xie, Yong

2018-03-01

The ultrasonic guided wave has emerged as a promising tool for structural health monitoring (SHM) and nondestructive testing (NDT) due to their capability to propagate over long distances with minimal loss and sensitivity to both surface and subsurface defects. The dispersion effect degrades the temporal and spatial resolution of guided waves. A novel ultrasonic guided wave processing method for both single mode and multi-mode guided waves dispersion compensation is proposed in this work based on compressed sensing, in which a dispersion signal dictionary is built by utilizing the dispersion curves of the guided wave modes in order to sparsely decompose the recorded dispersive guided waves. Dispersion-compensated guided waves are obtained by utilizing a non-dispersion signal dictionary and the results of sparse decomposition. Numerical simulations and experiments are implemented to verify the effectiveness of the developed method for both single mode and multi-mode guided waves.

Contaminant attenuation by shallow aquifer systems under steady flow

NASA Astrophysics Data System (ADS)

Soltani, S. S.; Cvetkovic, V.

2017-10-01

We present a framework for analyzing advection-dominated solute transport and transformation in aquifer systems of boreal catchments that are typically shallow and rest on crystalline bedrock. A methodology is presented for estimating tracer discharge based on particle trajectories from recharge to discharge locations and computing their first passage times assuming that the flow pattern is approximately steady-state. Transformation processes can be included by solving one-dimensional reactive transport with randomized water travel time as the independent variable; the distribution of the travel times incorporates morphological dispersion (due to catchment geometry/topography) as well as macro-dispersion (due to heterogeneity of underlying hydraulic properties). The implementation of the framework is illustrated for the well characterized coastal catchment of Forsmark (Sweden). We find that macro-dispersion has a notable effect on attenuation even though the morphological dispersion is significantly larger. Preferential flow on the catchment scale is found to be considerable with only 5% of the Eulerian velocities contributing to transport over the simulation period of 375 years. Natural attenuation is illustrated as a simple (linear decay) transformation process. Simulated natural attenuation can be estimated analytically reasonably well by using basic hydrological and structural information, the latter being the pathway length distribution and average aquifer depth to the bedrock.

Pore-scale and Continuum Simulations of Solute Transport Micromodel Benchmark Experiments

DOE Office of Scientific and Technical Information (OSTI.GOV)

Oostrom, Martinus; Mehmani, Yashar; Romero Gomez, Pedro DJ

Four sets of micromodel nonreactive solute transport experiments were conducted with flow velocity, grain diameter, pore-aspect ratio, and flow focusing heterogeneity as the variables. The data sets were offered to pore-scale modeling groups to test their simulators. Each set consisted of two learning experiments, for which all results was made available, and a challenge experiment, for which only the experimental description and base input parameters were provided. The experimental results showed a nonlinear dependence of the dispersion coefficient on the Peclet number, a negligible effect of the pore-aspect ratio on transverse mixing, and considerably enhanced mixing due to flow focusing.more » Five pore-scale models and one continuum-scale model were used to simulate the experiments. Of the pore-scale models, two used a pore-network (PN) method, two others are based on a lattice-Boltzmann (LB) approach, and one employed a computational fluid dynamics (CFD) technique. The learning experiments were used by the PN models to modify the standard perfect mixing approach in pore bodies into approaches to simulate the observed incomplete mixing. The LB and CFD models used these experiments to appropriately discretize the grid representations. The continuum model use published non-linear relations between transverse dispersion coefficients and Peclet numbers to compute the required dispersivity input values. Comparisons between experimental and numerical results for the four challenge experiments show that all pore-scale models were all able to satisfactorily simulate the experiments. The continuum model underestimated the required dispersivity values and, resulting in less dispersion. The PN models were able to complete the simulations in a few minutes, whereas the direct models needed up to several days on supercomputers to resolve the more complex problems.« less

Detailed source term estimation of the atmospheric release for the Fukushima Daiichi Nuclear Power Station accident by coupling simulations of an atmospheric dispersion model with an improved deposition scheme and oceanic dispersion model

NASA Astrophysics Data System (ADS)

Katata, G.; Chino, M.; Kobayashi, T.; Terada, H.; Ota, M.; Nagai, H.; Kajino, M.; Draxler, R.; Hort, M. C.; Malo, A.; Torii, T.; Sanada, Y.

2015-01-01

Temporal variations in the amount of radionuclides released into the atmosphere during the Fukushima Daiichi Nuclear Power Station (FNPS1) accident and their atmospheric and marine dispersion are essential to evaluate the environmental impacts and resultant radiological doses to the public. In this paper, we estimate the detailed atmospheric releases during the accident using a reverse estimation method which calculates the release rates of radionuclides by comparing measurements of air concentration of a radionuclide or its dose rate in the environment with the ones calculated by atmospheric and oceanic transport, dispersion and deposition models. The atmospheric and oceanic models used are WSPEEDI-II (Worldwide version of System for Prediction of Environmental Emergency Dose Information) and SEA-GEARN-FDM (Finite difference oceanic dispersion model), both developed by the authors. A sophisticated deposition scheme, which deals with dry and fog-water depositions, cloud condensation nuclei (CCN) activation, and subsequent wet scavenging due to mixed-phase cloud microphysics (in-cloud scavenging) for radioactive iodine gas (I2 and CH3I) and other particles (CsI, Cs, and Te), was incorporated into WSPEEDI-II to improve the surface deposition calculations. The results revealed that the major releases of radionuclides due to the FNPS1 accident occurred in the following periods during March 2011: the afternoon of 12 March due to the wet venting and hydrogen explosion at Unit 1, midnight of 14 March when the SRV (safety relief valve) was opened three times at Unit 2, the morning and night of 15 March, and the morning of 16 March. According to the simulation results, the highest radioactive contamination areas around FNPS1 were created from 15 to 16 March by complicated interactions among rainfall, plume movements, and the temporal variation of release rates. The simulation by WSPEEDI-II using the new source term reproduced the local and regional patterns of cumulative surface deposition of total 131I and 137Cs and air dose rate obtained by airborne surveys. The new source term was also tested using three atmospheric dispersion models (Modèle Lagrangien de Dispersion de Particules d'ordre zéro: MLDP0, Hybrid Single Particle Lagrangian Integrated Trajectory Model: HYSPLIT, and Met Office's Numerical Atmospheric-dispersion Modelling Environment: NAME) for regional and global calculations, and the calculated results showed good agreement with observed air concentration and surface deposition of 137Cs in eastern Japan.

Non-Parabolic Hydrodynamic Formulations for the Simulation of Inhomogeneous Semiconductor Devices

NASA Technical Reports Server (NTRS)

Smith, A. W.; Brennan, K. F.

1996-01-01

Hydrodynamic models are becoming prevalent design tools for small scale devices and other devices in which high energy effects can dominate transport. Most current hydrodynamic models use a parabolic band approximation to obtain fairly simple conservation equations. Interest in accounting for band structure effects in hydrodynamic device simulation has begun to grow since parabolic models cannot fully describe the transport in state of the art devices due to the distribution populating non-parabolic states within the band. This paper presents two different non-parabolic formulations or the hydrodynamic model suitable for the simulation of inhomogeneous semiconductor devices. The first formulation uses the Kane dispersion relationship ((hk)(exp 2)/2m = W(1 + alphaW). The second formulation makes use of a power law ((hk)(exp 2)/2m = xW(exp y)) for the dispersion relation. Hydrodynamic models which use the first formulation rely on the binomial expansion to obtain moment equations with closed form coefficients. This limits the energy range over which the model is valid. The power law formulation readily produces closed form coefficients similar to those obtained using the parabolic band approximation. However, the fitting parameters (x,y) are only valid over a limited energy range. The physical significance of the band non-parabolicity is discussed as well as the advantages/disadvantages and approximations of the two non-parabolic models. A companion paper describes device simulations based on the three dispersion relationships; parabolic, Kane dispersion and power law dispersion.

Non-parabolic hydrodynamic formulations for the simulation of inhomogeneous semiconductor devices

NASA Technical Reports Server (NTRS)

Smith, Arlynn W.; Brennan, Kevin F.

1995-01-01

Hydrodynamic models are becoming prevalent design tools for small scale devices and other devices in which high energy effects can dominate transport. Most current hydrodynamic models use a parabolic band approximation to obtain fairly simple conservation equations. Interest in accounting for band structure effects in hydrodynamic device simulation has begun to grow since parabolic models can not fully describe the transport in state of the art devices due to the distribution populating non-parabolic states within the band. This paper presents two different non-parabolic formulations of the hydrodynamic model suitable for the simulation of inhomogeneous semiconductor devices. The first formulation uses the Kane dispersion relationship (hk)(exp 2)/2m = W(1 + alpha(W)). The second formulation makes use of a power law ((hk)(exp 2)/2m = xW(sup y)) for the dispersion relation. Hydrodynamic models which use the first formulation rely on the binomial expansion to obtain moment equations with closed form coefficients. This limits the energy range over which the model is valid. The power law formulation readily produces closed form coefficients similar to those obtained using the parabolic band approximation. However, the fitting parameters (x,y) are only valid over a limited energy range. The physical significance of the band non-parabolicity is discussed as well as the advantages/disadvantages and approximations of the two non-parabolic models. A companion paper describes device simulations based on the three dispersion relationships: parabolic, Kane dispersion, and power low dispersion.

Assessment of sub-grid scale dispersion closure with regularized deconvolution method in a particle-laden turbulent jet

NASA Astrophysics Data System (ADS)

Wang, Qing; Zhao, Xinyu; Ihme, Matthias

2017-11-01

Particle-laden turbulent flows are important in numerous industrial applications, such as spray combustion engines, solar energy collectors etc. It is of interests to study this type of flows numerically, especially using large-eddy simulations (LES). However, capturing the turbulence-particle interaction in LES remains challenging due to the insufficient representation of the effect of sub-grid scale (SGS) dispersion. In the present work, a closure technique for the SGS dispersion using regularized deconvolution method (RDM) is assessed. RDM was proposed as the closure for the SGS dispersion in a counterflow spray that is studied numerically using finite difference method on a structured mesh. A presumed form of LES filter is used in the simulations. In the present study, this technique has been extended to finite volume method with an unstructured mesh, where no presumption on the filter form is required. The method is applied to a series of particle-laden turbulent jets. Parametric analyses of the model performance are conducted for flows with different Stokes numbers and Reynolds numbers. The results from LES will be compared against experiments and direct numerical simulations (DNS).

Dispersion and dispersion slope compensation impact on high channel bit rate optical signal transmission degradation

NASA Astrophysics Data System (ADS)

Hamidine, Mahamadou; Yuan, Xiuhua

2011-11-01

In this article a numerical simulation is carried out on a single channel optical transmission system with channel bit rate greater than 40 Gb/s to investigate optical signal degradation due to the impact of dispersion and dispersion slope of both transmitting and dispersion compensating fibers. By independently varying the input signal power and the dispersion slope of both transmitting and dispersion compensating fibers of an optical link utilizing a channel bit rate of 86 Gb/s, a good quality factor (Q factor) is obtained with a dispersion slope compensation ratio change of +/-10% for a faithful transmission. With this ratio change a minimum Q factor of 16 dB is obtained in the presence of amplifier noise figure of 5 dB and fiber nonlinearities effects at input signal power of 5 dBm and 3 spans of 100 km standard single mode fiber with a dispersion (D) value of 17 ps/nm.km.

Statistical mechanics of two-dimensional shuffled foams: Geometry-topology correlation in small or large disorder limits

NASA Astrophysics Data System (ADS)

Durand, Marc; Kraynik, Andrew M.; van Swol, Frank; Käfer, Jos; Quilliet, Catherine; Cox, Simon; Ataei Talebi, Shirin; Graner, François

2014-06-01

Bubble monolayers are model systems for experiments and simulations of two-dimensional packing problems of deformable objects. We explore the relation between the distributions of the number of bubble sides (topology) and the bubble areas (geometry) in the low liquid fraction limit. We use a statistical model [M. Durand, Europhys. Lett. 90, 60002 (2010), 10.1209/0295-5075/90/60002] which takes into account Plateau laws. We predict the correlation between geometrical disorder (bubble size dispersity) and topological disorder (width of bubble side number distribution) over an extended range of bubble size dispersities. Extensive data sets arising from shuffled foam experiments, surface evolver simulations, and cellular Potts model simulations all collapse surprisingly well and coincide with the model predictions, even at extremely high size dispersity. At moderate size dispersity, we recover our earlier approximate predictions [M. Durand, J. Kafer, C. Quilliet, S. Cox, S. A. Talebi, and F. Graner, Phys. Rev. Lett. 107, 168304 (2011), 10.1103/PhysRevLett.107.168304]. At extremely low dispersity, when approaching the perfectly regular honeycomb pattern, we study how both geometrical and topological disorders vanish. We identify a crystallization mechanism and explore it quantitatively in the case of bidisperse foams. Due to the deformability of the bubbles, foams can crystallize over a larger range of size dispersities than hard disks. The model predicts that the crystallization transition occurs when the ratio of largest to smallest bubble radii is 1.4.

Computation of turbulence and dispersion of cork in the NETL riser

DOE Office of Scientific and Technical Information (OSTI.GOV)

Jiradilok, Veeraya; Gidaspow, Dimitri; Breault, R.W.

The knowledge of dispersion coefficients is essential for reliable design of gasifiers. However, a literature review had shown that dispersion coefficients in fluidized beds differ by more than five orders of magnitude. This study presents a comparison of the computed axial solids dispersion coefficients for cork particles to the NETL riser cork data. The turbulence properties, the Reynolds stresses, the granular temperature spectra and the radial and axial gas and solids dispersion coefficients are computed. The standard kinetic theory model described in Gidaspow’s 1994 book, Multiphase Flow and Fluidization, Academic Press and the IIT and Fluent codes were used tomore » compute the measured axial solids volume fraction profiles for flow of cork particles in the NETL riser. The Johnson–Jackson boundary conditions were used. Standard drag correlations were used. This study shows that the computed solids volume fractions for the low flux flow are within the experimental error of those measured, using a two-dimensional model. At higher solids fluxes the simulated solids volume fractions are close to the experimental measurements, but deviate significantly at the top of the riser. This disagreement is due to use of simplified geometry in the two-dimensional simulation. There is a good agreement between the experiment and the three-dimensional simulation for a high flux condition. This study concludes that the axial and radial gas and solids dispersion coefficients in risers operating in the turbulent flow regime can be computed using a multiphase computational fluid dynamics model.« less

Direct exfoliation and dispersion of two-dimensional materials in pure water via temperature control.

PubMed

Kim, Jinseon; Kwon, Sanghyuk; Cho, Dae-Hyun; Kang, Byunggil; Kwon, Hyukjoon; Kim, Youngchan; Park, Sung O; Jung, Gwan Yeong; Shin, Eunhye; Kim, Wan-Gu; Lee, Hyungdong; Ryu, Gyeong Hee; Choi, Minseok; Kim, Tae Hyeong; Oh, Junghoon; Park, Sungjin; Kwak, Sang Kyu; Yoon, Suk Wang; Byun, Doyoung; Lee, Zonghoon; Lee, Changgu

2015-09-15

The high-volume synthesis of two-dimensional (2D) materials in the form of platelets is desirable for various applications. While water is considered an ideal dispersion medium, due to its abundance and low cost, the hydrophobicity of platelet surfaces has prohibited its widespread use. Here we exfoliate 2D materials directly in pure water without using any chemicals or surfactants. In order to exfoliate and disperse the materials in water, we elevate the temperature of the sonication bath, and introduce energy via the dissipation of sonic waves. Storage stability greater than one month is achieved through the maintenance of high temperatures, and through atomic and molecular level simulations, we further discover that good solubility in water is maintained due to the presence of platelet surface charges as a result of edge functionalization or intrinsic polarity. Finally, we demonstrate inkjet printing on hard and flexible substrates as a potential application of water-dispersed 2D materials.

Inversion of Surface-wave Dispersion Curves due to Low-velocity-layer Models

NASA Astrophysics Data System (ADS)

Shen, C.; Xia, J.; Mi, B.

2016-12-01

A successful inversion relies on exact forward modeling methods. It is a key step to accurately calculate multi-mode dispersion curves of a given model in high-frequency surface-wave (Rayleigh wave and Love wave) methods. For normal models (shear (S)-wave velocity increasing with depth), their theoretical dispersion curves completely match the dispersion spectrum that is generated based on wave equation. For models containing a low-velocity-layer, however, phase velocities calculated by existing forward-modeling algorithms (e.g. Thomson-Haskell algorithm, Knopoff algorithm, fast vector-transfer algorithm and so on) fail to be consistent with the dispersion spectrum at a high frequency range. They will approach a value that close to the surface-wave velocity of the low-velocity-layer under the surface layer, rather than that of the surface layer when their corresponding wavelengths are short enough. This phenomenon conflicts with the characteristics of surface waves, which results in an erroneous inverted model. By comparing the theoretical dispersion curves with simulated dispersion energy, we proposed a direct and essential solution to accurately compute surface-wave phase velocities due to low-velocity-layer models. Based on the proposed forward modeling technique, we can achieve correct inversion for these types of models. Several synthetic data proved the effectiveness of our method.

Performance of Dispersed Fringe Sensor in the Presence of Segmented Mirror Aberrations: Modeling and Simulation

NASA Technical Reports Server (NTRS)

Shi, Fang; Basinger, Scott A.; Redding, David C.

2006-01-01

Dispersed Fringe Sensing (DFS) is an efficient and robust method for coarse phasing of a segmented primary mirror such as the James Webb Space Telescope (JWST). In this paper, modeling and simulations are used to study the effect of segmented mirror aberrations on the fringe image, DFS signals and DFS detection accuracy. The study has shown due to the pixilation spatial filter effect from DFS signal extraction the effect of wavefront error is reduced and DFS algorithm will be more robust against wavefront aberration by using multi-trace DFS approach. We also studied the JWST Dispersed Hartmann Sensor (DHS) performance in presence of wavefront aberrations caused by the gravity sag and we use the scaled gravity sag to explore the JWST DHS performance relationship with the level of the wavefront aberration. This also includes the effect from line-of-sight jitter.

Comparison of Non-Parabolic Hydrodynamic Simulations for Semiconductor Devices

NASA Technical Reports Server (NTRS)

Smith, A. W.; Brennan, K. F.

1996-01-01

Parabolic drift-diffusion simulators are common engineering level design tools for semiconductor devices. Hydrodynamic simulators, based on the parabolic band approximation, are becoming more prevalent as device dimensions shrink and energy transport effects begin to dominate device characteristic. However, band structure effects present in state-of-the-art devices necessitate relaxing the parabolic band approximation. This paper presents simulations of ballistic diodes, a benchmark device, of Si and GaAs using two different non-parabolic hydrodynamic formulations. The first formulation uses the Kane dispersion relationship in the derivation of the conservation equations. The second model uses a power law dispersion relation {(hk)(exp 2)/2m = xW(exp Y)}. Current-voltage relations show that for the ballistic diodes considered. the non-parabolic formulations predict less current than the parabolic case. Explanations of this will be provided by examination of velocity and energy profiles. At low bias, the simulations based on the Kane formulation predict greater current flow than the power law formulation. As the bias is increased this trend changes and the power law predicts greater current than the Kane formulation. It will be shown that the non-parabolicity and energy range of the hydrodynamic model based on the Kane dispersion relation are limited due to the binomial approximation which was utilized in the derivation.

Development of a tool for calculating early internal doses in the Fukushima Daiichi nuclear power plant accident based on atmospheric dispersion simulation

NASA Astrophysics Data System (ADS)

Kurihara, Osamu; Kim, Eunjoo; Kunishima, Naoaki; Tani, Kotaro; Ishikawa, Tetsuo; Furuyama, Kazuo; Hashimoto, Shozo; Akashi, Makoto

2017-09-01

A tool was developed to facilitate the calculation of the early internal doses to residents involved in the Fukushima Nuclear Disaster based on atmospheric transport and dispersion model (ATDM) simulations performed using Worldwide version of System for Prediction of Environmental Emergency Information 2nd version (WSPEEDI-II) together with personal behavior data containing the history of the whereabouts of individul's after the accident. The tool generates hourly-averaged air concentration data for the simulation grids nearest to an individual's whereabouts using WSPEEDI-II datasets for the subsequent calculation of internal doses due to inhalation. This paper presents an overview of the developed tool and provides tentative comparisons between direct measurement-based and ATDM-based results regarding the internal doses received by 421 persons from whom personal behavior data available.

Jeans that fit: weighing the mass of the Milky Way analogues in the ΛCDM universe

NASA Astrophysics Data System (ADS)

Kafle, Prajwal R.; Sharma, Sanjib; Robotham, Aaron S. G.; Elahi, Pascal J.; Driver, Simon P.

2018-04-01

The spherical Jeans equation is a widely used tool for dynamical study of gravitating systems in astronomy. Here, we test its efficacy in robustly weighing the mass of Milky Way analogues, given they need not be in equilibrium or even spherical. Utilizing Milky Way stellar haloes simulated in accordance with Λ cold dark matter (ΛCDM) cosmology by Bullock and Johnston and analysing them under the Jeans formalism, we recover the underlying mass distribution of the parent galaxy, within distance r/kpc ∈ [10, 100], with a bias of ˜ 12 per cent and a dispersion of ˜ 14 per cent. Additionally, the mass profiles of triaxial dark matter haloes taken from the SURFS simulation, within scaled radius 0.2 < r/rmax < 3, are measured with a bias of ˜ - 2.4 per cent and a dispersion of ˜ 10 per cent. The obtained dispersion is not because of Poisson noise due to small particle numbers as it is twice the later. We interpret the dispersion to be due to the inherent nature of the ΛCDM haloes, for example being aspherical and out-of-equilibrium. Hence, the dispersion obtained for stellar haloes sets a limit of about 12 per cent (after adjusting for random uncertainty) on the accuracy with which the mass profiles of the Milky Way-like galaxies can be reconstructed using the spherical Jeans equation. This limit is independent of the quantity and quality of the observational data. The reason for a non-zero bias is not clear, hence its interpretation is not obvious at this stage.

Computational investigation of longitudinal diffusion, eddy dispersion, and trans-particle mass transfer in bulk, random packings of core-shell particles with varied shell thickness and shell diffusion coefficient.

PubMed

Daneyko, Anton; Hlushkou, Dzmitry; Baranau, Vasili; Khirevich, Siarhei; Seidel-Morgenstern, Andreas; Tallarek, Ulrich

2015-08-14

In recent years, chromatographic columns packed with core-shell particles have been widely used for efficient and fast separations at comparatively low operating pressure. However, the influence of the porous shell properties on the mass transfer kinetics in core-shell packings is still not fully understood. We report on results obtained with a modeling approach to simulate three-dimensional advective-diffusive transport in bulk random packings of monosized core-shell particles, covering a range of reduced mobile phase flow velocities from 0.5 up to 1000. The impact of the effective diffusivity of analyte molecules in the porous shell and the shell thickness on the resulting plate height was investigated. An extension of Giddings' theory of coupled eddy dispersion to account for retention of analyte molecules due to stagnant regions in porous shells with zero mobile phase flow velocity is presented. The plate height equation involving a modified eddy dispersion term excellently describes simulated data obtained for particle-packings with varied shell thickness and shell diffusion coefficient. It is confirmed that the model of trans-particle mass transfer resistance of core-shell particles by Kaczmarski and Guiochon [42] is applicable up to a constant factor. We analyze individual contributions to the plate height from different mass transfer mechanisms in dependence of the shell parameters. The simulations demonstrate that a reduction of plate height in packings of core-shell relative to fully porous particles arises mainly due to reduced trans-particle mass transfer resistance and transchannel eddy dispersion. Copyright © 2015 Elsevier B.V. All rights reserved.

Highly-dispersive electromagnetic induced transparency in planar symmetric metamaterials.

PubMed

Lu, Xiqun; Shi, Jinhui; Liu, Ran; Guan, Chunying

2012-07-30

We propose, design and experimentally demonstrate highly-dispersive electromagnetically induced transparency (EIT) in planar symmetric metamaterials actively switched and controlled by angles of incidence. Full-wave simulation and measurement results show EIT phenomena, trapped-mode excitations and the associated local field enhancement of two symmetric metamaterials consisting of symmetrically split rings (SSR) and a fishscale (FS) metamaterial pattern, respectively, strongly depend on angles of incidence. The FS metamaterial shows much broader spectral splitting than the SSR metamaterial due to the surface current distribution variation.

Pore-scale and continuum simulations of solute transport micromodel benchmark experiments

DOE PAGES

Oostrom, M.; Mehmani, Y.; Romero-Gomez, P.; ...

2014-06-18

Four sets of nonreactive solute transport experiments were conducted with micromodels. Three experiments with one variable, i.e., flow velocity, grain diameter, pore-aspect ratio, and flow-focusing heterogeneity were in each set. The data sets were offered to pore-scale modeling groups to test their numerical simulators. Each set consisted of two learning experiments, for which our results were made available, and one challenge experiment, for which only the experimental description and base input parameters were provided. The experimental results showed a nonlinear dependence of the transverse dispersion coefficient on the Peclet number, a negligible effect of the pore-aspect ratio on transverse mixing,more » and considerably enhanced mixing due to flow focusing. Five pore-scale models and one continuum-scale model were used to simulate the experiments. Of the pore-scale models, two used a pore-network (PN) method, two others are based on a lattice Boltzmann (LB) approach, and one used a computational fluid dynamics (CFD) technique. Furthermore, we used the learning experiments, by the PN models, to modify the standard perfect mixing approach in pore bodies into approaches to simulate the observed incomplete mixing. The LB and CFD models used the learning experiments to appropriately discretize the spatial grid representations. For the continuum modeling, the required dispersivity input values were estimated based on published nonlinear relations between transverse dispersion coefficients and Peclet number. Comparisons between experimental and numerical results for the four challenge experiments show that all pore-scale models were all able to satisfactorily simulate the experiments. The continuum model underestimated the required dispersivity values, resulting in reduced dispersion. The PN models were able to complete the simulations in a few minutes, whereas the direct models, which account for the micromodel geometry and underlying flow and transport physics, needed up to several days on supercomputers to resolve the more complex problems.« less

Solvation structures of water in trihexyltetradecylphosphonium-orthoborate ionic liquids

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wang, Yong-Lei, E-mail: wangyonl@gmail.com; System and Component Design, Department of Machine Design, KTH Royal Institute of Technology, SE-100 44 Stockholm; Sarman, Sten

2016-08-14

Atomistic molecular dynamics simulations have been performed to investigate effective interactions of isolated water molecules dispersed in trihexyltetradecylphosphonium-orthoborate ionic liquids (ILs). The intrinsic free energy changes in solvating one water molecule from gas phase into bulk IL matrices were estimated as a function of temperature, and thereafter, the calculations of potential of mean force between two dispersed water molecules within different IL matrices were performed using umbrella sampling simulations. The systematic analyses of local ionic microstructures, orientational preferences, probability and spatial distributions of dispersed water molecules around neighboring ionic species indicate their preferential coordinations to central polar segments in orthoboratemore » anions. The effective interactions between two dispersed water molecules are partially or totally screened as their separation distance increases due to interference of ionic species in between. These computational results connect microscopic anionic structures with macroscopically and experimentally observed difficulty in completely removing water from synthesized IL samples and suggest that the introduction of hydrophobic groups to central polar segments and the formation of conjugated ionic structures in orthoborate anions can effectively reduce residual water content in the corresponding IL samples.« less

Detailed source term estimation of the atmospheric release for the Fukushima Daiichi Nuclear Power Station accident by coupling simulations of atmospheric dispersion model with improved deposition scheme and oceanic dispersion model

NASA Astrophysics Data System (ADS)

Katata, G.; Chino, M.; Kobayashi, T.; Terada, H.; Ota, M.; Nagai, H.; Kajino, M.; Draxler, R.; Hort, M. C.; Malo, A.; Torii, T.; Sanada, Y.

2014-06-01

Temporal variations in the amount of radionuclides released into the atmosphere during the Fukushima Dai-ichi Nuclear Power Station (FNPS1) accident and their atmospheric and marine dispersion are essential to evaluate the environmental impacts and resultant radiological doses to the public. In this paper, we estimate a detailed time trend of atmospheric releases during the accident by combining environmental monitoring data with atmospheric model simulations from WSPEEDI-II (Worldwide version of System for Prediction of Environmental Emergency Dose Information), and simulations from the oceanic dispersion model SEA-GEARN-FDM, both developed by the authors. A sophisticated deposition scheme, which deals with dry and fogwater depositions, cloud condensation nuclei (CCN) activation and subsequent wet scavenging due to mixed-phase cloud microphysics (in-cloud scavenging) for radioactive iodine gas (I2 and CH3I) and other particles (CsI, Cs, and Te), was incorporated into WSPEEDI-II to improve the surface deposition calculations. The fallout to the ocean surface calculated by WSPEEDI-II was used as input data for the SEA-GEARN-FDM calculations. Reverse and inverse source-term estimation methods based on coupling the simulations from both models was adopted using air dose rates and concentrations, and sea surface concentrations. The results revealed that the major releases of radionuclides due to FNPS1 accident occurred in the following periods during March 2011: the afternoon of 12 March due to the wet venting and hydrogen explosion at Unit 1, the morning of 13 March after the venting event at Unit 3, midnight of 14 March when the SRV (Safely Relief Valve) at Unit 2 was opened three times, the morning and night of 15 March, and the morning of 16 March. According to the simulation results, the highest radioactive contamination areas around FNPS1 were created from 15 to 16 March by complicated interactions among rainfall, plume movements, and the temporal variation of release rates associated with reactor pressure changes in Units 2 and 3. The modified WSPEEDI-II simulation using the new source term reproduced local and regional patterns of cumulative surface deposition of total 131I and 137Cs and air dose rate obtained by airborne surveys. The new source term was also tested using three atmospheric dispersion models (MLDP0, HYSPLIT, and NAME) for regional and global calculations and showed good agreement between calculated and observed air concentration and surface deposition of 137Cs in East Japan. Moreover, HYSPLIT model using the new source term also reproduced the plume arrivals at several countries abroad showing a good correlation with measured air concentration data. A large part of deposition pattern of total 131I and 137Cs in East Japan was explained by in-cloud particulate scavenging. However, for the regional scale contaminated areas, there were large uncertainties due to the overestimation of rainfall amounts and the underestimation of fogwater and drizzle depositions. The computations showed that approximately 27% of 137Cs discharged from FNPS1 deposited to the land in East Japan, mostly in forest areas.

Ion sound instability driven by the ion flows

DOE Office of Scientific and Technical Information (OSTI.GOV)

Koshkarov, O., E-mail: koshkarov.alexandr@usask.ca; Smolyakov, A. I.; National Research Centre

2015-05-15

Ion sound instabilities driven by the ion flow in a system of a finite length are considered by analytical and numerical methods. The ion sound waves are modified by the presence of stationary ion flow resulting in negative and positive energy modes. The instability develops due to coupling of negative and positive energy modes mediated by reflections from the boundary. It is shown that the wave dispersion due to deviation from quasineutrality is crucial for the stability. In finite length system, the dispersion is characterized by the length of the system measured in units of the Debye length. The instabilitymore » is studied analytically and the results are compared with direct, initial value numerical simulations.« less

CORRECTING FOR INTERSTELLAR SCATTERING DELAY IN HIGH-PRECISION PULSAR TIMING: SIMULATION RESULTS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Palliyaguru, Nipuni; McLaughlin, Maura; Stinebring, Daniel

2015-12-20

Light travel time changes due to gravitational waves (GWs) may be detected within the next decade through precision timing of millisecond pulsars. Removal of frequency-dependent interstellar medium (ISM) delays due to dispersion and scattering is a key issue in the detection process. Current timing algorithms routinely correct pulse times of arrival (TOAs) for time-variable delays due to cold plasma dispersion. However, none of the major pulsar timing groups correct for delays due to scattering from multi-path propagation in the ISM. Scattering introduces a frequency-dependent phase change in the signal that results in pulse broadening and arrival time delays. Any methodmore » to correct the TOA for interstellar propagation effects must be based on multi-frequency measurements that can effectively separate dispersion and scattering delay terms from frequency-independent perturbations such as those due to a GW. Cyclic spectroscopy, first described in an astronomical context by Demorest (2011), is a potentially powerful tool to assist in this multi-frequency decomposition. As a step toward a more comprehensive ISM propagation delay correction, we demonstrate through a simulation that we can accurately recover impulse response functions (IRFs), such as those that would be introduced by multi-path scattering, with a realistic signal-to-noise ratio (S/N). We demonstrate that timing precision is improved when scatter-corrected TOAs are used, under the assumptions of a high S/N and highly scattered signal. We also show that the effect of pulse-to-pulse “jitter” is not a serious problem for IRF reconstruction, at least for jitter levels comparable to those observed in several bright pulsars.« less

Automated electric valve for electrokinetic separation in a networked microfluidic chip.

PubMed

Cui, Huanchun; Huang, Zheng; Dutta, Prashanta; Ivory, Cornelius F

2007-02-15

This paper describes an automated electric valve system designed to reduce dispersion and sample loss into a side channel when an electrokinetically mobilized concentration zone passes a T-junction in a networked microfluidic chip. One way to reduce dispersion is to control current streamlines since charged species are driven along them in the absence of electroosmotic flow. Computer simulations demonstrate that dispersion and sample loss can be reduced by applying a constant additional electric field in the side channel to straighten current streamlines in linear electrokinetic flow (zone electrophoresis). This additional electric field was provided by a pair of platinum microelectrodes integrated into the chip in the vicinity of the T-junction. Both simulations and experiments of this electric valve with constant valve voltages were shown to provide unsatisfactory valve performance during nonlinear electrophoresis (isotachophoresis). On the basis of these results, however, an automated electric valve system was developed with improved valve performance. Experiments conducted with this system showed decreased dispersion and increased reproducibility as protein zones isotachophoretically passed the T-junction. Simulations of the automated electric valve offer further support that the desired shape of current streamlines was maintained at the T-junction during isotachophoresis. Valve performance was evaluated at different valve currents based on statistical variance due to dispersion. With the automated control system, two integrated microelectrodes provide an effective way to manipulate current streamlines, thus acting as an electric valve for charged species in electrokinetic separations.

MODEL -- BIOPLUME III VERSION 1.0 - SEPTEMBER 1997 (SUBSURFACE PROTECTION AND REMEDIATION DIVISION, NRMRL)

EPA Science Inventory

BIOPLUME III is a 2D, finite difference model for simulating the natural attenuation of organic contaminants in groundwater due to the processes of advection, dispersion, sorption, and biodegradation. Biotransformation processes are potentially important in the restoration of aq...

Life history trade-off moderates model predictions of diversity loss from climate change.

PubMed

Moor, Helen

2017-01-01

Climate change can trigger species range shifts, local extinctions and changes in diversity. Species interactions and dispersal capacity are important mediators of community responses to climate change. The interaction between multispecies competition and variation in dispersal capacity has recently been shown to exacerbate the effects of climate change on diversity and to increase predictions of extinction risk dramatically. Dispersal capacity, however, is part of a species' overall ecological strategy and are likely to trade off with other aspects of its life history that influence population growth and persistence. In plants, a well-known example is the trade-off between seed mass and seed number. The presence of such a trade-off might buffer the diversity loss predicted by models with random but neutral (i.e. not impacting fitness otherwise) differences in dispersal capacity. Using a trait-based metacommunity model along a warming climatic gradient the effect of three different dispersal scenarios on model predictions of diversity change were compared. Adding random variation in species dispersal capacity caused extinctions by the introduction of strong fitness differences due an inherent property of the dispersal kernel. Simulations including a fitness-equalising trade-off based on empirical relationships between seed mass (here affecting dispersal distance, establishment probability, and seedling biomass) and seed number (fecundity) maintained higher initial species diversity and predicted lower extinction risk and diversity loss during climate change than simulations with variable dispersal capacity. Large seeded species persisted during climate change, but developed lags behind their climate niche that may cause extinction debts. Small seeded species were more extinction-prone during climate change but tracked their niches through dispersal and colonisation, despite competitive resistance from residents. Life history trade-offs involved in coexistence mechanisms may increase community resilience to future climate change and are useful guides for model development.

Impact of subgrid fluid turbulence on inertial particles subject to gravity

NASA Astrophysics Data System (ADS)

Rosa, Bogdan; Pozorski, Jacek

2017-07-01

Two-phase turbulent flows with the dispersed phase in the form of small, spherical particles are increasingly often computed with the large-eddy simulation (LES) of the carrier fluid phase, coupled to the Lagrangian tracking of particles. To enable further model development for LES with inertial particles subject to gravity, we consider direct numerical simulations of homogeneous isotropic turbulence with a large-scale forcing. Simulation results, both without filtering and in the a priori LES setting, are reported and discussed. A full (i.e. a posteriori) LES is also performed with the spectral eddy viscosity. Effects of gravity on the dispersed phase include changes in the average settling velocity due to preferential sweeping, impact on the radial distribution function and radial relative velocity, as well as direction-dependent modification of the particle velocity variance. The filtering of the fluid velocity, performed in spectral space, is shown to have a non-trivial impact on these quantities.

Evaluation of a micro-scale wind model's performance over realistic building clusters using wind tunnel experiments

NASA Astrophysics Data System (ADS)

Zhang, Ning; Du, Yunsong; Miao, Shiguang; Fang, Xiaoyi

2016-08-01

The simulation performance over complex building clusters of a wind simulation model (Wind Information Field Fast Analysis model, WIFFA) in a micro-scale air pollutant dispersion model system (Urban Microscale Air Pollution dispersion Simulation model, UMAPS) is evaluated using various wind tunnel experimental data including the CEDVAL (Compilation of Experimental Data for Validation of Micro-Scale Dispersion Models) wind tunnel experiment data and the NJU-FZ experiment data (Nanjing University-Fang Zhuang neighborhood wind tunnel experiment data). The results show that the wind model can reproduce the vortexes triggered by urban buildings well, and the flow patterns in urban street canyons and building clusters can also be represented. Due to the complex shapes of buildings and their distributions, the simulation deviations/discrepancies from the measurements are usually caused by the simplification of the building shapes and the determination of the key zone sizes. The computational efficiencies of different cases are also discussed in this paper. The model has a high computational efficiency compared to traditional numerical models that solve the Navier-Stokes equations, and can produce very high-resolution (1-5 m) wind fields of a complex neighborhood scale urban building canopy (~ 1 km ×1 km) in less than 3 min when run on a personal computer.

The Atacama Cosmology Telescope: Dynamical Masses for 44 SZ-Selected Galaxy Clusters over 755 Square Degrees

NASA Technical Reports Server (NTRS)

Sifon, Cristobal; Battaglia, Nick; Hasselfield, Matthew; Menanteau, Felipe; Barrientos, L. Felipe; Bond, J. Richard; Crichton, Devin; Devlin, Mark J.; Dunner, Rolando; Hilton, Matt;

Carbon Nanotube Dispersion in Solvents and Polymer Solutions: Mechanisms, Assembly, and Preferences.

PubMed

Pramanik, Chandrani; Gissinger, Jacob R; Kumar, Satish; Heinz, Hendrik

2017-12-26

Debundling and dispersion of carbon nanotubes (CNTs) in polymer solutions play a major role in the preparation of carbon nanofibers due to early effects on interfacial ordering and mechanical properties. A roadblock toward ultrastrong fibers is the difficulty to achieve homogeneous dispersions of CNTs in polyacrylonitrile (PAN) and poly(methyl methacrylate) (PMMA) precursor solutions in solvents such as dimethyl sulfoxide (DMSO), N,N-dimethylacetamide (DMAc), and N,N-dimethylformamide (DMF). In this contribution, molecular dynamics simulations with accurate interatomic potentials for graphitic materials that include virtual π electrons are reported to analyze the interaction of pristine single wall CNTs with the solvents and polymer solutions at 25 °C. The results explain the barriers toward dispersion of SWCNTs and quantify CNT-solvent, polymer-solvent, as well as CNT-polymer interactions in atomic detail. Debundling of CNTs is overall endothermic and unfavorable with dispersion energies of +20 to +30 mJ/m 2 in the pure solvents, + 20 to +40 mJ/m 2 in PAN solutions, and +20 to +60 mJ/m 2 in PMMA solutions. Differences arise due to molecular geometry, polar, van der Waals, and CH-π interactions. Among the pure solvents, DMF restricts CNT dispersion less due to the planar geometry and stronger van der Waals interactions. PAN and PMMA interact favorably with the pure solvents with dissolution energies of -0.7 to -1.1 kcal per mole monomer and -1.5 to -2.2 kcal per mole monomer, respectively. Adsorption of PMMA onto CNTs is stronger than that of PAN in all solvents as the molecular geometry enables more van der Waals contacts between alkyl groups and the CNT surface. Polar side groups in both polymers prefer interactions with the polar solvents. Higher polymer concentrations in solution lead to polymer aggregation via alkyl groups and reduce adsorption onto CNTs. PAN and PMMA solutions in DMSO and dilute solutions in DMF support CNT dispersion more than other combinations whereby the polymers significantly adsorb onto CNTs in DMSO solution. The observations by molecular simulations are consistent with available experimental data and solubility parameters and aid in the design of carbon nanofibers. The methods can be applied to other multiphase graphitic materials.

One-dimensional Analytical Modelling of Floating Seed Dispersal in Tidal Channels

NASA Astrophysics Data System (ADS)

Shi, W.; Purnama, A.; Shao, D.; Cui, B.; Gao, W.

2017-12-01

Seed dispersal is a primary factor influencing plant community development, and thus plays a critical role in maintaining wetland ecosystem functioning. However, compared with fluvial seed dispersal of riparian plants, dispersal of saltmarsh plant seeds in tidal channels is much less studied due to its complex behavior, and relevant mathematical modelling is particularly lacking. In this study, we developed a one-dimensional advection-dispersion model to explore the patterns of tidal seed dispersal. Oscillatory tidal current and water depth were assumed to represent the tidal effects. An exponential decay coefficient λ was introduced to account for seed deposition and retention. Analytical solution in integral form was derived using Green's function and further evaluated using numerical integration. The developed model was applied to simulate Spartina densiflora seed dispersal in a tidal channel located at the Mad River Slough in North Humboldt Bay, California, USA, to demonstrate its practical applicability. Model predictions agree satisfactorily with field observation and simulation results from Delft3D numerical model. Sensitivity analyses were also conducted to evaluate the effects of varying calibrated parameters on model predictions. The range of the seed dispersion as well as the distribution of the seed concentration were further analyzed through statistical parameters such as centroid displacement and variance of the seed cloud together with seed concentration contours. Implications of the modelling results on tidal marsh restoration and protection, e.g., revegetation through seed addition, were also discussed through scenario analysis. The developed analytical model provides a useful tool for ecological management of tidal marshes.

Solubilization of single-walled carbon nanotubes using a peptide aptamer in water below the critical micelle concentration.

PubMed

Li, Zha; Kameda, Tomoshi; Isoshima, Takashi; Kobatake, Eiry; Tanaka, Takeshi; Ito, Yoshihiro; Kawamoto, Masuki

2015-03-24

The solubilizing ability of single-walled carbon nanotubes (SWCNTs) in water with several dispersants was investigated. Among the dispersants, including low-molecular-weight surfactants, peptides, DNA, and a water-soluble polymer, the peptide aptamer, A2 (IFRLSWGTYFS), exhibited the highest dispersion capability below the critical micelle concentration at a concentration of 0.02 w/v%. The dispersion of supernatant aqueous solution of SWCNTs containing aptamer A2 was essentially unchanged for several months after high-speed ultracentrifugation and gave rise to an efficient and stable dispersion of the SWCNTs in water. From the results of isothermal titration calorimetry and molecular dynamics simulations, the effective binding capability of A2 was due to π-π interaction between aromatic groups in the peptide aptamer and the side walls of SWCNTs. Interestingly, the peptide aptamer showed the possibility of diameter separation of semiconducting SWCNTs using a uniform density gradient ultracentrifuge. These phenomena are encouraging results toward an effective approach to the dispersion and separation of SWCNTs.

Modeling Intrajunction Dispersion at a Well-Mixed Tidal River Junction

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wolfram, Phillip J.; Fringer, Oliver B.; Monsen, Nancy E.

In this paper, the relative importance of small-scale, intrajunction flow features such as shear layers, separation zones, and secondary flows on dispersion in a well-mixed tidal river junction is explored. A fully nonlinear, nonhydrostatic, and unstructured three-dimensional (3D) model is used to resolve supertidal dispersion via scalar transport at a well-mixed tidal river junction. Mass transport simulated in the junction is compared against predictions using a simple node-channel model to quantify the effects of small-scale, 3D intrajunction flow features on mixing and dispersion. The effects of three-dimensionality are demonstrated by quantifying the difference between two-dimensional (2D) and 3D model results.more » An intermediate 3D model that does not resolve the secondary circulation or the recirculating flow at the junction is also compared to the 3D model to quantify the relative sensitivity of mixing on intrajunction flow features. Resolution of complex flow features simulated by the full 3D model is not always necessary because mixing is primarily governed by bulk flow splitting due to the confluence–diffluence cycle. Finally, results in 3D are comparable to the 2D case for many flow pathways simulated, suggesting that 2D modeling may be reasonable for nonstratified and predominantly hydrostatic flows through relatively straight junctions, but not necessarily for the full junction network.« less

Modeling Intrajunction Dispersion at a Well-Mixed Tidal River Junction

DOE PAGES

Wolfram, Phillip J.; Fringer, Oliver B.; Monsen, Nancy E.; ...

2016-08-01

In this paper, the relative importance of small-scale, intrajunction flow features such as shear layers, separation zones, and secondary flows on dispersion in a well-mixed tidal river junction is explored. A fully nonlinear, nonhydrostatic, and unstructured three-dimensional (3D) model is used to resolve supertidal dispersion via scalar transport at a well-mixed tidal river junction. Mass transport simulated in the junction is compared against predictions using a simple node-channel model to quantify the effects of small-scale, 3D intrajunction flow features on mixing and dispersion. The effects of three-dimensionality are demonstrated by quantifying the difference between two-dimensional (2D) and 3D model results.more » An intermediate 3D model that does not resolve the secondary circulation or the recirculating flow at the junction is also compared to the 3D model to quantify the relative sensitivity of mixing on intrajunction flow features. Resolution of complex flow features simulated by the full 3D model is not always necessary because mixing is primarily governed by bulk flow splitting due to the confluence–diffluence cycle. Finally, results in 3D are comparable to the 2D case for many flow pathways simulated, suggesting that 2D modeling may be reasonable for nonstratified and predominantly hydrostatic flows through relatively straight junctions, but not necessarily for the full junction network.« less

Numerical simulation of tidal dispersion around a coastal headland

USGS Publications Warehouse

Signell, R.P.; Geyer, W. Rockwell; Cheng, Ralph T.

1990-01-01

Tidal flows around headlands can exhibit strong spatial gradients in the Eulerian currents, resulting in complex Lagrangian trajectories and dispersion of the vertically integrated flow. This typically occurs when the horizontal length scale of the headland is comparable to or smaller than the tidal excursion. The effects of these headlands on dispersion are investigated using a depthaveraged hydrodynamic model combined with a particle tracking model. The dispersion of patches of fluid is found to vary by more than an order of magnitude, depending both on position and tidal phase at the time of release. This is due to the infrequent interaction of material with the strongly sheared flow at the tip of the headland, where flow separation occurs during times of maximum tidal flow. Spreading of these patches over many tidal cycles is not Gaussian, but rather shows a patchy, streaky structure.

Dispersive stresses in wind farms

NASA Astrophysics Data System (ADS)

Segalini, Antonio; Braunbehrens, Robert; Hyvarinen, Ann

2017-11-01

One of the most famous models of wind farms is provided by the assumption that the farm can be approximated as a horizontally-homogeneous forest canopy with vertically-varying force intensity. By means of this approximation, the flow-motion equations become drastically simpler, as many of the three-dimensional effects are gone. However, the application of the horizontal average operator to the RANS equations leads to the appearance of new transport terms (called dispersive stresses) originating from the horizontal (small-scale) variation of the mean velocity field. Since these terms are related to the individual turbine signature, they are expected to vanish outside the roughness sublayer, providing a definition for the latter. In the present work, an assessment of the dispersive stresses is performed by means of a wake-model approach and through the linearised code ORFEUS developed at KTH. Both approaches are very fast and enable the characterization of a large number of wind-farm layouts. The dispersive stress tensor and its effect on the turbulence closure models are investigated, providing guidelines for those simulations where it is impossible to resolve the farm at a turbine scale due to grid requirements (as, for instance, mesoscale simulations).

CFD Modeling of LNG Spill: Humidity Effect on Vapor Dispersion

NASA Astrophysics Data System (ADS)

Giannissi, S. G.; Venetsanos, A. G.; Markatos, N.

2015-09-01

The risks entailed by an accidental spill of Liquefied Natural Gas (LNG) should be indentified and evaluated, in order to design measures for prevention and mitigation in LNG terminals. For this purpose, simulations are considered a useful tool to study LNG spills and to understand the mechanisms that influence the vapor dispersion. In the present study, the ADREA-HF CFD code is employed to simulate the TEEX1 experiment. The experiment was carried out at the Brayton Fire Training Field, which is affiliated with the Texas A&M University system and involves LNG release and dispersion over water surface in open- obstructed environment. In the simulation the source was modeled as a two-phase jet enabling the prediction of both the vapor dispersion and the liquid pool spreading. The conservation equations for the mixture are solved along with the mass fraction for natural gas. Due to the low prevailing temperatures during the spill ambient humidity condenses and this might affect the vapor dispersion. This effect was examined in this work by solving an additional conservation equation for the water mass fraction. Two different models were tested: the hydrodynamic equilibrium model which assumes kinetic equilibrium between the phases and the non hydrodynamic equilibrium model, in order to assess the effect of slip velocity on the prediction. The slip velocity is defined as the difference between the liquid phase and the vapor phase and is calculated using the algebraic slip model. Constant droplet diameter of three different sizes and a lognormal distribution of the droplet diameter were applied and the results are discussed and compared with the measurements.

Simulating dispersion in porous media and the influence of segmentation on stagnancy in carbonates

NASA Astrophysics Data System (ADS)

Gray, F.; Cen, J.; Shah, S. M.; Crawshaw, J. P.; Boek, E. S.

2016-11-01

Understanding the transport of chemical components in porous media is fundamentally important to many reservoir processes such as contaminant transport and reactive flows involved in CO2 sequestration. Carbonate rocks in particular present difficulties for pore-scale simulations because they contain large amounts of sub-micron porosity. In this work, we introduce a new hybrid simulation model to calculate hydrodynamic dispersion in pore-scale images of real porous media and use this to elucidate the origins and behaviour of stagnant zones arising in transport simulations using micro-CT images of carbonates. For this purpose a stochastic particle model for simulating the transport of a solute is coupled to a Lattice-Boltzmann algorithm to calculate the flow field. The particle method incorporates second order spatial and temporal resolution to resolve finer features of the domain. We demonstrate how dispersion coefficients can be accurately obtained in capillaries, where corresponding analytical solutions are available, even when these are resolved to just a few lattice units. Then we compute molecular displacement distributions for pore-spaces of varying complexity: a pack of beads; a Bentheimer sandstone; and a Portland carbonate. Our calculated propagator distributions are compared directly with recent experimental PFG-NMR propagator distributions (Scheven et al., 2005; Mitchell et al., 2008), the latter excluding spin relaxation mechanisms. We observe that the calculated transport propagators can be quantitatively compared with the experimental distribution, provided that spin relaxations in the experiment are excluded, and good agreement is found for both the sandstone and the carbonate. However, due to the absence of explicit micro-porosity from the carbonate pore space image used for flow field simulations we note that there are fundamental differences in the physical origins of the stagnant zones for micro-porous rocks between simulation and experiment. We show that for a given micro-CT image of a carbonate, small variations in the parameters chosen for the segmentation process lead to different amounts of stagnancy which diffuse away at different rates. Finally, we use a filtering method to show that this is due to the presence of spurious isolated pores which arise from the segmentation process and suggest an approach to overcome this limitation.

Neutral Evolution and Dispersal Limitation Produce Biogeographic Patterns in Microcystis aeruginosa Populations of Lake Systems.

PubMed

Shirani, Sahar; Hellweger, Ferdi L

2017-08-01

Molecular observations reveal substantial biogeographic patterns of cyanobacteria within systems of connected lakes. An important question is the relative role of environmental selection and neutral processes in the biogeography of these systems. Here, we quantify the effect of genetic drift and dispersal limitation by simulating individual cyanobacteria cells using an agent-based model (ABM). In the model, cells grow (divide), die, and migrate between lakes. Each cell has a full genome that is subject to neutral mutation (i.e., the growth rate is independent of the genome). The model is verified by simulating simplified lake systems, for which theoretical solutions are available. Then, it is used to simulate the biogeography of the cyanobacterium Microcystis aeruginosa in a number of real systems, including the Great Lakes, Klamath River, Yahara River, and Chattahoochee River. Model output is analyzed using standard bioinformatics tools (BLAST, MAFFT). The emergent patterns of nucleotide divergence between lakes are dynamic, including gradual increases due to accumulation of mutations and abrupt changes due to population takeovers by migrant cells (coalescence events). The model predicted nucleotide divergence is heterogeneous within systems, and for weakly connected lakes, it can be substantial. For example, Lakes Superior and Michigan are predicted to have an average genomic nucleotide divergence of 8200 bp or 0.14%. The divergence between more strongly connected lakes is much lower. Our results provide a quantitative baseline for future biogeography studies. They show that dispersal limitation can be an important factor in microbe biogeography, which is contrary to the common belief, and could affect how a system responds to environmental change.

Comparison of ground motions from hybrid simulations to nga prediction equations

USGS Publications Warehouse

Star, L.M.; Stewart, J.P.; Graves, R.W.

2011-01-01

We compare simulated motions for a Mw 7.8 rupture scenario on the San Andreas Fault known as the ShakeOut event, two permutations with different hypocenter locations, and a Mw 7.15 Puente Hills blind thrust scenario, to median and dispersion predictions from empirical NGA ground motion prediction equations. We find the simulated motions attenuate faster with distance than is predicted by the NGA models for periods less than about 5.0 s After removing this distance attenuation bias, the average residuals of the simulated events (i.e., event terms) are generally within the scatter of empirical event terms, although the ShakeOut simulation appears to be a high static stress drop event. The intraevent dispersion in the simulations is lower than NGA values at short periods and abruptly increases at 1.0 s due to different simulation procedures at short and long periods. The simulated motions have a depth-dependent basin response similar to the NGA models, and also show complex effects in which stronger basin response occurs when the fault rupture transmits energy into a basin at low angle, which is not predicted by the NGA models. Rupture directivity effects are found to scale with the isochrone parameter ?? 2011, Earthquake Engineering Research Institute.

Delineating the role of ripples on the thermal expansion of 2D honeycomb materials: graphene, 2D h-BN and monolayer (ML)-MoS2.

PubMed

Anees, P; Valsakumar, M C; Panigrahi, B K

2017-04-19

We delineated the role of thermally excited ripples on the thermal expansion properties of 2D honeycomb materials (free-standing graphene, 2D h-BN, and ML-MoS 2 ), by explicitly carrying out three-dimensional (3D) and two-dimensional (2D) molecular dynamics simulations. In 3D simulations, the in-plane lattice parameter (a-lattice) of graphene and 2D h-BN shows thermal contraction over a wide range of temperatures and exhibits a strong system size dependence. The 2D simulations of the very same system show a reverse trend, where the a-lattice expands in the whole computed temperature range. In contrast to graphene and 2D h-BN, the a-lattice of ML-MoS 2 shows thermal expansion in both 2D and 3D simulations and their system size dependence is marginal. By analyzing the phonon dispersion at 300 K, we found that the discrepancy between 2D and 3D simulations of graphene and 2D h-BN is due to the absence of out-of-plane bending modes (ZA) in 2D simulations, which is responsible for the thermal contraction of the a-lattice at low temperature. Meanwhile, all the phonon modes are present in the 2D phonon dispersion of ML-MoS 2 , which indicates that the origin of the ZA mode is not purely due to the out-of-plane movement of atoms and also its effect on thermal expansion is not significant as found in graphene and 2D h-BN.

Examination of climatological wind patterns and simulated pollen dispersion in a complex island environment

NASA Astrophysics Data System (ADS)

Viner, Brian J.; Arritt, Raymond W.; Westgate, Mark E.

2017-08-01

Complex terrain creates small-scale circulations which affect pollen dispersion but may be missed by meteorological observing networks and coarse-grid meteorological models. On volcanic islands, these circulations result from differing rates of surface heating between land and sea as well as rugged terrain. We simulated the transport of bentgrass, ryegrass, and maize pollen from 30 sources within the agricultural regions of the Hawaiian island Kaua'i during climatological conditions spanning season conditions and the La Niña, El Niño, and neutral phases of the El Niño-Southern Oscillation. Both pollen size and source location had major effects on predicted dispersion over and near the island. Three patterns of pollen dispersion were identified in response to prevailing wind conditions: southwest winds transported pollen inland, funneling pollen grains through valleys; east winds transported pollen over the ocean, with dispersive tails for the smallest pollen grains following the mean wind and extending as far as the island of Ni'ihau 35 km away; and northeast winds moved pollen inland counter to the prevailing flow due to a sea breeze circulation that formed over the source region. These results are the first to predict the interactions between complex island terrain and local climatology on grass pollen dispersion. They demonstrate how numerical modeling can provide guidance for field trials by illustrating the common flow regimes present in complex terrain, allowing field trials to focus on areas where successful sampling is more likely to occur.

Dispersion-convolution model for simulating peaks in a flow injection system.

PubMed

Pai, Su-Cheng; Lai, Yee-Hwong; Chiao, Ling-Yun; Yu, Tiing

2007-01-12

A dispersion-convolution model is proposed for simulating peak shapes in a single-line flow injection system. It is based on the assumption that an injected sample plug is expanded due to a "bulk" dispersion mechanism along the length coordinate, and that after traveling over a distance or a period of time, the sample zone will develop into a Gaussian-like distribution. This spatial pattern is further transformed to a temporal coordinate by a convolution process, and finally a temporal peak image is generated. The feasibility of the proposed model has been examined by experiments with various coil lengths, sample sizes and pumping rates. An empirical dispersion coefficient (D*) can be estimated by using the observed peak position, height and area (tp*, h* and At*) from a recorder. An empirical temporal shift (Phi*) can be further approximated by Phi*=D*/u2, which becomes an important parameter in the restoration of experimental peaks. Also, the dispersion coefficient can be expressed as a second-order polynomial function of the pumping rate Q, for which D*(Q)=delta0+delta1Q+delta2Q2. The optimal dispersion occurs at a pumping rate of Qopt=sqrt[delta0/delta2]. This explains the interesting "Nike-swoosh" relationship between the peak height and pumping rate. The excellent coherence of theoretical and experimental peak shapes confirms that the temporal distortion effect is the dominating reason to explain the peak asymmetry in flow injection analysis.

THE EFFECT OF UNRESOLVED BINARIES ON GLOBULAR CLUSTER PROPER-MOTION DISPERSION PROFILES

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bianchini, P.; Norris, M. A.; Ven, G. van de

2016-03-20

High-precision kinematic studies of globular clusters (GCs) require an accurate knowledge of all possible sources of contamination. Among other sources, binary stars can introduce systematic biases in the kinematics. Using a set of Monte Carlo cluster simulations with different concentrations and binary fractions, we investigate the effect of unresolved binaries on proper-motion dispersion profiles, treating the simulations like Hubble Space Telescope proper-motion samples. Since GCs evolve toward a state of partial energy equipartition, more-massive stars lose energy and decrease their velocity dispersion. As a consequence, on average, binaries have a lower velocity dispersion, since they are more-massive kinematic tracers. Wemore » show that, in the case of clusters with high binary fractions (initial binary fractions of 50%) and high concentrations (i.e., closer to energy equipartition), unresolved binaries introduce a color-dependent bias in the velocity dispersion of main-sequence stars of the order of 0.1–0.3 km s{sup −1} (corresponding to 1%−6% of the velocity dispersion), with the reddest stars having a lower velocity dispersion, due to the higher fraction of contaminating binaries. This bias depends on the ability to distinguish binaries from single stars, on the details of the color–magnitude diagram and the photometric errors. We apply our analysis to the HSTPROMO data set of NGC 7078 (M15) and show that no effect ascribable to binaries is observed, consistent with the low binary fraction of the cluster. Our work indicates that binaries do not significantly bias proper-motion velocity-dispersion profiles, but should be taken into account in the error budget of kinematic analyses.« less

Study on longitudinal dispersion relation in one-dimensional relativistic plasma: Linear theory and Vlasov simulation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhang, H.; Wu, S. Z.; Zhou, C. T.

2013-09-15

The dispersion relation of one-dimensional longitudinal plasma waves in relativistic homogeneous plasmas is investigated with both linear theory and Vlasov simulation in this paper. From the Vlasov-Poisson equations, the linear dispersion relation is derived for the proper one-dimensional Jüttner distribution. Numerically obtained linear dispersion relation as well as an approximate formula for plasma wave frequency in the long wavelength limit is given. The dispersion of longitudinal wave is also simulated with a relativistic Vlasov code. The real and imaginary parts of dispersion relation are well studied by varying wave number and plasma temperature. Simulation results are in agreement with establishedmore » linear theory.« less

A numerical study on flow and pollutant transport in Singapore coastal waters.

PubMed

Xu, Ming; Chua, Vivien P

2016-10-15

Intensive economic and shipping activities in Singapore Strait have caused Singapore coastal waters to be under high risk of water pollution. A nested three-dimensional unstructured-grid SUNTANS model is applied to Singapore coastal waters to simulate flow and pollutant transport. The small domain (~50m resolution) Singapore coastal model is nested within a large domain (~200m resolution) regional model. The nested model is able to predict water surface elevations and velocities with high R(2) values of 0.96 and 0.91, respectively. Model results delineate the characteristics of circulation pattern in Singapore coastal waters during the Northeast and Southwest monsoons. The pollutants are modeled as passive tracers, and are released at six key sailing locations Points 1-6 in Singapore coastal waters and are named as Passive Tracers 1-6, respectively. Our results show that the rate of dispersion is twice as large for the Northeast monsoon compared to the Southwest monsoon due to differences in large-scale monsoons and small-scale local winds. The volume averaged concentration (VAC) diminishes faster and the local flushing time is shorter during the Northeast monsoon than the Southwest monsoon. Dispersion coefficients K and the VAC decreasing rate are maximum for Tracers 2 and 3 with shortest local flushing time due to the strong surrounding currents and abrupt bathymetry changes near Senang and St. John Islands. Dispersion coefficients K and the VAC decreasing rate are minimum for Tracer 1 due to weak currents induced by the semi-enclosed coastline near Tuas. It is found that both the lateral dispersion coefficient Ky and the compound dispersion coefficient K obey a "4/3-law", which defines a linear correlation between dispersion coefficients and 4/3-power of selected length scale. Copyright © 2016 Elsevier Ltd. All rights reserved.

Water quality modeling in the dead end sections of drinking water distribution networks.

PubMed

Abokifa, Ahmed A; Yang, Y Jeffrey; Lo, Cynthia S; Biswas, Pratim

2016-02-01

Dead-end sections of drinking water distribution networks are known to be problematic zones in terms of water quality degradation. Extended residence time due to water stagnation leads to rapid reduction of disinfectant residuals allowing the regrowth of microbial pathogens. Water quality models developed so far apply spatial aggregation and temporal averaging techniques for hydraulic parameters by assigning hourly averaged water demands to the main nodes of the network. Although this practice has generally resulted in minimal loss of accuracy for the predicted disinfectant concentrations in main water transmission lines, this is not the case for the peripheries of the distribution network. This study proposes a new approach for simulating disinfectant residuals in dead end pipes while accounting for both spatial and temporal variability in hydraulic and transport parameters. A stochastic demand generator was developed to represent residential water pulses based on a non-homogenous Poisson process. Dispersive solute transport was considered using highly dynamic dispersion rates. A genetic algorithm was used to calibrate the axial hydraulic profile of the dead-end pipe based on the different demand shares of the withdrawal nodes. A parametric sensitivity analysis was done to assess the model performance under variation of different simulation parameters. A group of Monte-Carlo ensembles was carried out to investigate the influence of spatial and temporal variations in flow demands on the simulation accuracy. A set of three correction factors were analytically derived to adjust residence time, dispersion rate and wall demand to overcome simulation error caused by spatial aggregation approximation. The current model results show better agreement with field-measured concentrations of conservative fluoride tracer and free chlorine disinfectant than the simulations of recent advection dispersion reaction models published in the literature. Accuracy of the simulated concentration profiles showed significant dependence on the spatial distribution of the flow demands compared to temporal variation. Copyright © 2015 Elsevier Ltd. All rights reserved.

A fuzzy-theory-based method for studying the effect of information transmission on nonlinear crowd dispersion dynamics

NASA Astrophysics Data System (ADS)

Fu, Libi; Song, Weiguo; Lo, Siuming

2017-01-01

Emergencies involved in mass events are related to a variety of factors and processes. An important factor is the transmission of information on danger that has an influence on nonlinear crowd dynamics during the process of crowd dispersion. Due to much uncertainty in this process, there is an urgent need to propose a method to investigate the influence. In this paper, a novel fuzzy-theory-based method is presented to study crowd dynamics under the influence of information transmission. Fuzzy functions and rules are designed for the ambiguous description of human states. Reasonable inference is employed to decide the output values of decision making such as pedestrian movement speed and directions. Through simulation under four-way pedestrian situations, good crowd dispersion phenomena are achieved. Simulation results under different conditions demonstrate that information transmission cannot always induce successful crowd dispersion in all situations. This depends on whether decision strategies in response to information on danger are unified and effective, especially in dense crowds. Results also suggest that an increase in drift strength at low density and the percentage of pedestrians, who choose one of the furthest unoccupied Von Neumann neighbors from the dangerous source as the drift direction at high density, is helpful in crowd dispersion. Compared with previous work, our comprehensive study improves an in-depth understanding of nonlinear crowd dynamics under the effect of information on danger.

Numerical studies of dispersion due to tidal flow through Moskstraumen, northern Norway

NASA Astrophysics Data System (ADS)

Lynge, Birgit Kjoss; Berntsen, Jarle; Gjevik, Bjørn

2010-08-01

The effect of horizontal grid resolution on the horizontal relative dispersion of particle pairs has been investigated on a short time scale, i.e. one tidal M 2 cycle. Of particular interest is the tidal effect on dispersion and transports in coastal waters where small-scale flow features are important. A three-dimensional ocean model has been applied to simulate the tidal flow through the Moskstraumen Maelstrom outside Lofoten in northern Norway, well known for its strong current and whirlpools (Gjevik et al., Nature 388(6645):837-838, 1997; Moe et al., Cont Shelf Res 22(3):485-504, 2002). Simulations with spatial resolution down to 50 m have been carried out. Lagrangian tracers were passively advected with the flow, and Lyapunov exponents and power law exponents have been calculated to analyse the separation statistics. It is found that the relative dispersion of particles on a short time scale (12-24 h) is very sensitive to the grid size and that the spatial variability is also very large, ranging from 0 to 100 km2 over a distance of 100 m. This means that models for prediction of transport and dispersion of oil spills, fish eggs, sea lice etc. using a single diffusion coefficient will be of limited value, unless the models actually resolves the small-scale eddies of the tidal current.

Dispersive effects on multicomponent transport through porous media

NASA Astrophysics Data System (ADS)

Dutta, Sourav; Daripa, Prabir

2017-11-01

We use a hybrid numerical method to solve a global pressure based porous media flow model of chemical enhanced oil recovery. This is an extension of our recent work. The numerical method is based on the use of a discontinuous finite element method and the modified method of characteristics. The impact of molecular diffusion and mechanical dispersion on the evolution of scalar concentration distributions are studied through numerical simulations of various flooding schemes. The relative importance of the advective, capillary diffusive and dispersive fluxes are compared over different flow regimes defined in the parameter space of Capillary number, Peclet number, longitudinal and transverse dispersion coefficients. Such studies are relevant for the design of effective injection policies and determining optimal combinations of chemical components for improving recovery. This work has been possible due to financial support from the U.S. National Science Foundation Grant DMS-1522782.

Effects of optical layer impairments on 2.5 Gb/s optical CDMA transmission.

PubMed

Feng, H; Mendez, A; Heritage, J; Lennon, W

2000-07-03

We conducted a computer simulation study to assess the effects of optical layer impairments on optical CDMA (O-CDMA) transmission of 8 asynchronous users at 2.5 Gb/s each user over a 214-km link. It was found that with group velocity dispersion compensation, two other residual effects, namely, the nonzero chromatic dispersion slope of the single mode fiber (which causes skew) and the non-uniform EDFA gain (which causes interference power level to exceed signal power level of some codes) degrade the signal to multi-access interference (MAI) ratio. In contrast, four wave mixing and modulation due to the Kerr and Raman contributions to the fiber nonlinear refractive index are less important. Current wavelength-division multiplexing (WDM) technologies, including dispersion management, EDFA gain flattening, and 3 rd order dispersion compensation, are sufficient to overcome the impairments to the O-CDMA transmission system that we considered.

Effect of tidal fluctuations on transient dispersion of simulated contaminant concentrations in coastal aquifers

USGS Publications Warehouse

La Licata, Ivana; Langevin, Christian D.; Dausman, Alyssa M.; Alberti, Luca

2011-01-01

Variable-density groundwater models require extensive computational resources, particularly for simulations representing short-term hydrologic variability such as tidal fluctuations. Saltwater-intrusion models usually neglect tidal fluctuations and this may introduce errors in simulated concentrations. The effects of tides on simulated concentrations in a coastal aquifer were assessed. Three analyses are reported: in the first, simulations with and without tides were compared for three different dispersivity values. Tides do not significantly affect the transfer of a hypothetical contaminant into the ocean; however, the concentration difference between tidal and non-tidal simulations could be as much as 15%. In the second analysis, the dispersivity value for the model without tides was increased in a zone near the ocean boundary. By slightly increasing dispersivity in this zone, the maximum concentration difference between the simulations with and without tides was reduced to as low as 7%. In the last analysis, an apparent dispersivity value was calculated for each model cell using the simulated velocity variations from the model with tides. Use of apparent dispersivity values in models with a constant ocean boundary seems to provide a reasonable approach for approximating tidal effects in simulations where explicit representation of tidal fluctuations is not feasible.

Effect of tidal fluctuations on transient dispersion of simulated contaminant concentrations in coastal aquifers

USGS Publications Warehouse

La Licata, Ivana; Langevin, Christian D.; Dausman, Alyssa M.; Alberti, Luca

2013-01-01

Variable-density groundwater models require extensive computational resources, particularly for simulations representing short-term hydrologic variability such as tidal fluctuations. Saltwater-intrusion models usually neglect tidal fluctuations and this may introduce errors in simulated concentrations. The effects of tides on simulated concentrations in a coastal aquifer were assessed. Three analyses are reported: in the first, simulations with and without tides were compared for three different dispersivity values. Tides do not significantly affect the transfer of a hypothetical contaminant into the ocean; however, the concentration difference between tidal and non-tidal simulations could be as much as 15%. In the second analysis, the dispersivity value for the model without tides was increased in a zone near the ocean boundary. By slightly increasing dispersivity in this zone, the maximum concentration difference between the simulations with and without tides was reduced to as low as 7%. In the last analysis, an apparent dispersivity value was calculated for each model cell using the simulated velocity variations from the model with tides. Use of apparent dispersivity values in models with a constant ocean boundary seems to provide a reasonable approach for approximating tidal effects in simulations where explicit representation of tidal fluctuations is not feasible.

Analysis of the coherent and turbulent stresses of a numerically simulated rough wall pipe

NASA Astrophysics Data System (ADS)

Chan, L.; MacDonald, M.; Chung, D.; Hutchins, N.; Ooi, A.

2017-04-01

A turbulent rough wall flow in a pipe is simulated using direct numerical simulation (DNS) where the roughness elements consist of explicitly gridded three-dimensional sinusoids. Two groups of simulations were conducted where the roughness semi-amplitude h+ and the roughness wavelength λ+ are systematically varied. The triple decomposition is applied to the velocity to separate the coherent and turbulent components. The coherent or dispersive component arises due to the roughness and depends on the topological features of the surface. The turbulent stress on the other hand, scales with the friction Reynolds number. For the case with the largest roughness wavelength, large secondary flows are observed which are similar to that of duct flows. The occurrence of these large secondary flows is due to the spanwise heterogeneity of the roughness which has a spacing approximately equal to the boundary layer thickness δ.

Life history trade-off moderates model predictions of diversity loss from climate change

PubMed Central

2017-01-01

Climate change can trigger species range shifts, local extinctions and changes in diversity. Species interactions and dispersal capacity are important mediators of community responses to climate change. The interaction between multispecies competition and variation in dispersal capacity has recently been shown to exacerbate the effects of climate change on diversity and to increase predictions of extinction risk dramatically. Dispersal capacity, however, is part of a species’ overall ecological strategy and are likely to trade off with other aspects of its life history that influence population growth and persistence. In plants, a well-known example is the trade-off between seed mass and seed number. The presence of such a trade-off might buffer the diversity loss predicted by models with random but neutral (i.e. not impacting fitness otherwise) differences in dispersal capacity. Using a trait-based metacommunity model along a warming climatic gradient the effect of three different dispersal scenarios on model predictions of diversity change were compared. Adding random variation in species dispersal capacity caused extinctions by the introduction of strong fitness differences due an inherent property of the dispersal kernel. Simulations including a fitness-equalising trade-off based on empirical relationships between seed mass (here affecting dispersal distance, establishment probability, and seedling biomass) and seed number (fecundity) maintained higher initial species diversity and predicted lower extinction risk and diversity loss during climate change than simulations with variable dispersal capacity. Large seeded species persisted during climate change, but developed lags behind their climate niche that may cause extinction debts. Small seeded species were more extinction-prone during climate change but tracked their niches through dispersal and colonisation, despite competitive resistance from residents. Life history trade-offs involved in coexistence mechanisms may increase community resilience to future climate change and are useful guides for model development. PMID:28520770

Classical density functional theory and Monte Carlo simulation study of electric double layer in the vicinity of a cylindrical electrode

NASA Astrophysics Data System (ADS)

Zhou, Shiqi; Lamperski, Stanisław; Sokołowska, Marta

2017-07-01

We have performed extensive Monte-Carlo simulations and classical density functional theory (DFT) calculations of the electrical double layer (EDL) near a cylindrical electrode in a primitive model (PM) modified by incorporating interionic dispersion interactions. It is concluded that (i) in general, an unsophisticated use of the mean field (MF) approximation for the interionic dispersion interactions does not distinctly worsen the classical DFT performance, even if the salt ions considered are highly asymmetrical in size (3:1) and charge (5:1), the bulk molar concentration considered is high up to a total bulk ion packing fraction of 0.314, and the surface charge density of up to 0.5 C m-2. (ii) More specifically, considering the possible noises in the simulation, the local volume charge density profiles are the most accurately predicted by the classical DFT in all situations, and the co- and counter-ion singlet distributions are also rather accurately predicted; whereas the mean electrostatic potential profile is relatively less accurately predicted due to an integral amplification of minor inaccuracy of the singlet distributions. (iii) It is found that the layered structure of the co-ion distribution is abnormally possible only if the surface charge density is high enough (for example 0.5 C m-2) moreover, the co-ion valence abnormally influences the peak height of the first counter-ion layer, which decreases with the former. (iv) Even if both the simulation and DFT indicate an insignificant contribution of the interionic dispersion interaction to the above three ‘local’ quantities, it is clearly shown by the classical DFT that the interionic dispersion interaction does significantly influence a ‘global’ quantity like the cylinder surface-aqueous electrolyte interfacial tension, and this may imply the role of the interionic dispersion interaction in explaining the specific Hofmeister effects. We elucidate all of the above observations based on the arguments from the liquid state theory and at the molecular scale.

Quantum Drude oscillator model of atoms and molecules: Many-body polarization and dispersion interactions for atomistic simulation

NASA Astrophysics Data System (ADS)

Jones, Andrew P.; Crain, Jason; Sokhan, Vlad P.; Whitfield, Troy W.; Martyna, Glenn J.

2013-04-01

Treating both many-body polarization and dispersion interactions is now recognized as a key element in achieving the level of atomistic modeling required to reveal novel physics in complex systems. The quantum Drude oscillator (QDO), a Gaussian-based, coarse grained electronic structure model, captures both many-body polarization and dispersion and has linear scale computational complexity with system size, hence it is a leading candidate next-generation simulation method. Here, we investigate the extent to which the QDO treatment reproduces the desired long-range atomic and molecular properties. We present closed form expressions for leading order polarizabilities and dispersion coefficients and derive invariant (parameter-free) scaling relationships among multipole polarizability and many-body dispersion coefficients that arise due to the Gaussian nature of the model. We show that these “combining rules” hold to within a few percent for noble gas atoms, alkali metals, and simple (first-row hydride) molecules such as water; this is consistent with the surprising success that models with underlying Gaussian statistics often exhibit in physics. We present a diagrammatic Jastrow-type perturbation theory tailored to the QDO model that serves to illustrate the rich types of responses that the QDO approach engenders. QDO models for neon, argon, krypton, and xenon, designed to reproduce gas phase properties, are constructed and their condensed phase properties explored via linear scale diffusion Monte Carlo (DMC) and path integral molecular dynamics (PIMD) simulations. Good agreement with experimental data for structure, cohesive energy, and bulk modulus is found, demonstrating a degree of transferability that cannot be achieved using current empirical models or fully ab initio descriptions.

Analysis and iterative equalization of transient and adiabatic chirp effects in DML-based OFDM transmission systems.

PubMed

Wei, Chia-Chien

2012-11-05

This work theoretically studies the transmission performance of a DML-based OFDM system by small-signal approximation, and the model considers both the transient and adiabatic chirps. The dispersion-induced distortion is modeled as subcarrier-to-subcarrier intermixing interference (SSII), and the theoretical SSII agrees with the distortion obtained from large-signal simulation statistically and deterministically. The analysis shows that the presence of the adiabatic chirp will ease power fading or even provide gain, but will increase the SSII to deteriorate OFDM signals after dispersive transmission. Furthermore, this work also proposes a novel iterative equalization to eliminate the SSII. From the simulation, the distortion could be effectively mitigated by the proposed equalization such that the maximum transmission distance of the DML-based OFDM signal is significantly improved. For instance, the transmission distance of a 30-Gbps DML-based OFDM signal can be extended from 10 km to more than 100 km. Besides, since the dispersion-induced distortion could be effectively mitigated by the equalization, negative power penalties are observed at some distances due to chirp-induced power gain.

The edge of galaxy formation - II. Evolution of Milky Way satellite analogues after infall

NASA Astrophysics Data System (ADS)

Frings, Jonas; Macciò, Andrea; Buck, Tobias; Penzo, Camilla; Dutton, Aaron; Blank, Marvin; Obreja, Aura

2017-12-01

In the first paper, we presented 27 hydrodynamical cosmological simulations of galaxies with total masses between 5 × 108 and 1010 M⊙. In this second paper, we use a subset of these cosmological simulations as initial conditions (ICs) for more than 40 hydrodynamical simulations of satellite and host galaxy interaction. Our cosmological ICs seem to suggest that galaxies on these mass scales have very little rotational support and are velocity dispersion (σ) dominated. Accretion and environmental effects increase the scatter in the galaxy scaling relations (e.g. size-velocity dispersion) in very good agreement with observations. Star formation is substantially quenched after accretion. Mass removal due to tidal forces has several effects: it creates a very flat stellar velocity dispersion profile, and it reduces the dark matter content at all scales (even in the centre), which in turn lowers the stellar velocity on scales around 0.5 kpc even when the galaxy does not lose stellar mass. Satellites which start with a cored dark matter profile are more prone to either be destroyed or to end up in a very dark matter poor galaxy. Finally, we found that tidal effects always increase the 'cuspyness' of the dark matter profile, even for haloes that infall with a core.

Dispersive Phase in the L-band InSAR Image Associated with Heavy Rain Episodes

NASA Astrophysics Data System (ADS)

Furuya, M.; Kinoshita, Y.

2017-12-01

Interferometric synthetic aperture radar (InSAR) is a powerful geodetic technique that allows us to detect ground displacements with unprecedented spatial resolution, and has been used to detect displacements due to earthquakes, volcanic eruptions, and glacier motion. In the meantime, due to the microwave propagation through ionosphere and troposphere, we often encounter non-negligible phase anomaly in InSAR data. Correcting for the ionsphere and troposphere is therefore a long-standing issue for high-precision geodetic measurements. However, if ground displacements are negligible, InSAR image can tell us the details of the atmosphere.Kinoshita and Furuya (2017, SOLA) detected phase anomaly in ALOS/PALSAR InSAR data associated with heavy rain over Niigata area, Japan, and performed numerical weathr model simulation to reproduce the anomaly; ALOS/PALSAR is a satellite-based L-band SAR sensor launched by JAXA in 2006 and terminated in 2011. The phase anomaly could be largely reproduced, using the output data from the weather model. However, we should note that numerical weather model outputs can only account for the non-dispersive effect in the phase anomaly. In case of severe weather event, we may expect dispersive effect that could be caused by the presence of free-electrons.In Global Navigation Satellite System (GNSS) positioning, dual frequency measurements allow us to separate the ionospheric dispersive component from tropospheric non-dispersive components. In contrast, SAR imaging is based on a single carrier frequency, and thus no operational ionospheric corrections have been performed in InSAR data analyses. Recently, Gomba et al (2016) detailed the processing strategy of split spectrum method (SSM) for InSAR, which splits the finite bandwidth of the range spectrum and virtually allows for dual-frequency measurements.We apply the L-band InSAR SSM to the heavy rain episodes, in which more than 50 mm/hour precipitations were reported. We report the presence of phase anomaly in both dispersive and non-dispersive components. While the original phase anomaly turns out to be mostly due to the non-dispersive effect, we could recognize local anomalies in the dispersive component as well. We will discuss its geophysical implications, and may show several case studies.

Theoretic Study on Dispersion Mechanism of Boron Nitride Nanotubes by Polynucleotides

PubMed Central

Liang, Lijun; Hu, Wei; Zhang, Zhisen; Shen, Jia-Wei

2016-01-01

Due to the unique electrical and mechanical properties of boron nitride nanotubes (BNNT), BNNT has been a promising material for many potential applications, especially in biomedical field. Understanding the dispersion of BNNT in aqueous solution by biomolecules is essential for its use in biomedical applications. In this study, BNNT wrapped by polynucleotides in aqueous solution was investigated by molecular dynamics (MD) simulations. Our results demonstrated that the BNNT wrapped by polynucleotides could greatly hinder the aggregation of BNNTs and improve the dispersion of BNNTs in aqueous solution. Dispersion of BNNTs with the assistance of polynucleotides is greatly affected by the wrapping manner of polynucleotides on BNNT, which mainly depends on two factors: the type of polynucleotides and the radius of BNNT. The interaction between polynucleotides and BNNT(9, 9) is larger than that between polynucleotides and BNNT(5, 5), which leads to the fact that dispersion of BNNT(9, 9) is better than that of BNNT(5, 5) with the assistance of polynucleotides in aqueous solution. Our study revealed the molecular-level dispersion mechanism of BNNT with the assistance of polynucleotides in aqueous solution. It shades a light on the understanding of dispersion of single wall nanotubes by biomolecules. PMID:28004832

Flow and axial dispersion in a sinusoidal-walled tube: Effects of inertial and unsteady flows

DOE Office of Scientific and Technical Information (OSTI.GOV)

Richmond, Marshall C.; Perkins, William A.; Scheibe, Timothy D.

2013-12-01

Dispersion in porous media flows has been the subject of much experimental, theoretical and numerical study. Here we consider a wavy-walled tube (a three-dimensional tube with sinusoidally-varying diameter) as a simplified conceptualization of flow in porous media, where constrictions represent pore throats and expansions pore bodies. A theoretical model for effective (macroscopic) longitudinal dispersion in this system has been developed by volume averaging the microscale velocity field. Direct numerical simulation using computational fluid dynamics (CFD) methods was used to compute velocity fields by solving the Navier-Stokes equations, and also to numerically solve the volume averaging closure problem, for a rangemore » of Reynolds numbers (Re) spanning the low-Re to inertial flow regimes, including one simulation at Re = 449 for which unsteady flow was observed. Dispersion values were computed using both the volume averaging solution and a random walk particle tracking method, and results of the two methods were shown to be consistent. Our results are compared to experimental measurements of dispersion in porous media and to previous theoretical results for the low-Re, Stokes flow regime. In the steady inertial regime we observe an power-law increase in effective longitudinal dispersion (DL) with Re, consistent with previous results. This rapid rate of increase is caused by trapping of solute in expansions due to flow separation (eddies). For the unsteady case (Re = 449), the rate of increase of DL with Re was smaller than that observed at lower Re. Velocity fluctuations in this regime lead to increased rates of solute mass transfer between the core flow and separated flow regions, thus diminishing the amount of tailing caused by solute trapping in eddies and thereby reducing longitudinal dispersion.« less

Assessment of Fiber Chromatic Dispersion Based on Elimination of Second-Order Harmonics in Optical OFDM Single Sideband Modulation Using Mach Zehnder Modulator

NASA Astrophysics Data System (ADS)

Patel, Dhananjay; Singh, Vinay Kumar; Dalal, U. D.

2016-07-01

This work addresses the analytical and numerical investigations of the transmission performance of an optical Single Sideband (SSB) modulation technique generated by a Mach Zehnder Modulator (MZM) with a 90° and 120° hybrid coupler. It takes into account the problem of chromatic dispersion in single mode fibers in Passive Optical Networks (PON), which severely degrades the performance of the system. Considering the transmission length of the fiber, the SSB modulation generated by maintaining a phase shift of π/2 between the two electrodes of the MZM provides better receiver sensitivity. However, the power of higher-order harmonics generated due to the nonlinearity of the MZM is directly proportional to the modulation index, making the SSB look like a quasi-double sideband (DSB) and causing power fading due to chromatic dispersion. To eliminate one of the second-order harmonics, the SSB signal based on an MZM with a 120° hybrid coupler is simulated. An analytical model of conventional SSB using 90° and 120° hybrid couplers is established. The latter suppresses unwanted (upper/lower) first-order and second-order (lower/upper) sidebands. For the analysis, a varying quadrature amplitude modulation (QAM) Orthogonal Frequency Division Multiplexing (OFDM) signal with a data rate of 5 Gb/s is upconverted using both of the SSB techniques and is transmitted over a distance of 75 km in Single Mode Fiber (SMF). The simulation results show that the SSB with 120° hybrid coupler proves to be more immune to chromatic dispersion as compared to the conventional SSB technique. This is in tandem with the theoretical analysis presented in the article.

Mate-finding as an overlooked critical determinant of dispersal variation in sexually-reproducing animals.

PubMed

Gilroy, James J; Lockwood, Julie L

2012-01-01

Dispersal is a critically important process in ecology, but robust predictive models of animal dispersal remain elusive. We identify a potentially ubiquitous component of variation in animal dispersal that has been largely overlooked until now: the influence of mate encounters on settlement probability. We use an individual-based model to simulate dispersal in sexually-reproducing organisms that follow a simple set of movement rules based on conspecific encounters, within an environment lacking spatial habitat heterogeneity. We show that dispersal distances vary dramatically with fluctuations in population density in such a model, even in the absence of variation in dispersive traits between individuals. In a simple random-walk model with promiscuous mating, dispersal distributions become increasingly 'fat-tailed' at low population densities due to the increasing scarcity of mates. Similar variation arises in models incorporating territoriality. In a model with polygynous mating, we show that patterns of sex-biased dispersal can even be reversed across a gradient of population density, despite underlying dispersal mechanisms remaining unchanged. We show that some widespread dispersal patterns found in nature (e.g. fat tailed distributions) can arise as a result of demographic variability in the absence of heterogeneity in dispersive traits across the population. This implies that models in which individual dispersal distances are considered to be fixed traits might be unrealistic, as dispersal distances vary widely under a single dispersal mechanism when settlement is influenced by mate encounters. Mechanistic models offer a promising means of advancing our understanding of dispersal in sexually-reproducing organisms.

Dispersal and fallout simulations for urban consequences management (u)

DOE Office of Scientific and Technical Information (OSTI.GOV)

Grinstein, Fernando F; Wachtor, Adam J; Nelson, Matt

2010-01-01

Hazardous chemical, biological, or radioactive releases from leaks, spills, fires, or blasts, may occur (intentionally or accidentally) in urban environments during warfare or as part of terrorist attacks on military bases or other facilities. The associated contaminant dispersion is complex and semi-chaotic. Urban predictive simulation capabilities can have direct impact in many threat-reduction areas of interest, including, urban sensor placement and threat analysis, contaminant transport (CT) effects on surrounding civilian population (dosages, evacuation, shelter-in-place), education and training of rescue teams and services. Detailed simulations for the various processes involved are in principle possible, but generally not fast. Predicting urban airflowmore » accompanied by CT presents extremely challenging requirements. Crucial technical issues include, simulating turbulent fluid and particulate transport, initial and boundary condition modeling incorporating a consistent stratified urban boundary layer with realistic wind fluctuations, and post-processing of the simulation results for practical consequences management. Relevant fluid dynamic processes to be simulated include, detailed energetic and contaminant sources, complex building vortex shedding and flows in recirculation zones, and modeling of particle distributions, including particulate fallout, as well as deposition, re-suspension and evaporation. Other issues include, modeling building damage effects due to eventual blasts, addressing appropriate regional and atmospheric data reduction.« less

Influence of regional climate change on meteorological characteristics and their subsequent effect on ozone dispersion in Taiwan

NASA Astrophysics Data System (ADS)

Cheng, Fang-Yi; Jian, Shan-Ping; Yang, Zhih-Min; Yen, Ming-Cheng; Tsuang, Ben-Jei

2015-02-01

The objective of this study is to understand the influence of regional climate change on local meteorological conditions and their subsequent effect on local ozone (O3) dispersion in Taiwan. The 33-year NCEP-DOE Reanalysis 2 (NNR2) data set (1979-2011) was analyzed to understand the variations in regional-scale atmospheric conditions in East Asia and the western North Pacific. To save computational processing time, two scenarios representative of past (1979-86) and current (2004-11) atmospheric conditions were selected but only targeting the autumn season (September, October and November) when the O3 concentrations were at high levels. Numerical simulations were performed using weather research and forecasting (WRF) model and Community Multiscale Air Quality (CMAQ) model for the past and current scenarios individually but only for the month of October because of limited computational resources. Analysis of NNR2 data exhibited increased air temperature, weakened Asian continental anticyclone, enhanced northeasterly monsoonal flow, and a deepened low-pressure system forming near Taiwan. With enhanced evaporation from oceans along with a deepened low-pressure system, precipitation amounts increased in Taiwan in the current scenario. As demonstrated in the WRF simulation, the land surface physical process responded to the enhanced precipitation resulting in damper soil conditions, and reduced ground temperatures that in turn restricted the development of boundary layer height. The weakened land-sea breeze flow was simulated in the current scenario. With reduced dispersion capability, air pollutants would tend to accumulate near the emission source leading to a degradation of air quality in this region. The conditions would be even worse in southwestern Taiwan due to the fact that stagnant wind fields would occur more frequently in the current scenario. On the other hand, in northern Taiwan, the simulated O3 concentrations are lower during the day in the current scenario due to the enhanced cloud conditions and reduced solar radiation.

Examination of climatological wind patterns and simulated pollen dispersion in a complex island environment

DOE PAGES

Viner, Brian J.; Arritt, Raymond W.; Westgate, Mark E.

2017-03-29

Complex terrain creates small-scale circulations which affect pollen dispersion but may be missed by meteorological observing networks and coarse-grid meteorological models. On volcanic islands, these circulations result from differing rates of surface heating between land and sea as well as rugged terrain. We simulated the transport of bentgrass, ryegrass, and maize pollen from 30 sources within the agricultural regions of the Hawaiian island Kaua’i during climatological conditions spanning season conditions and the La Niña, El Niño, and neutral phases of the El Niño-Southern Oscillation. Both pollen size and source location had major effects on predicted dispersion over and near themore » island. Three patterns of pollen dispersion were identified in response to prevailing wind conditions: southwest winds transported pollen inland, funneling pollen grains through valleys; east winds transported pollen over the ocean, with dispersive tails for the smallest pollen grains following the mean wind and extending as far as the island of Ni’ihau 35 km away; and northeast winds moved pollen inland counter to the prevailing flow due to a sea breeze circulation that formed over the source region. These results are the first to predict the interactions between complex island terrain and local climatology on grass pollen dispersion. As a result, they demonstrate how numerical modeling can provide guidance for field trials by illustrating the common flow regimes present in complex terrain, allowing field trials to focus on areas where successful sampling is more likely to occur.« less

Examination of climatological wind patterns and simulated pollen dispersion in a complex island environment

DOE Office of Scientific and Technical Information (OSTI.GOV)

Viner, Brian J.; Arritt, Raymond W.; Westgate, Mark E.

Complex terrain creates small-scale circulations which affect pollen dispersion but may be missed by meteorological observing networks and coarse-grid meteorological models. On volcanic islands, these circulations result from differing rates of surface heating between land and sea as well as rugged terrain. We simulated the transport of bentgrass, ryegrass, and maize pollen from 30 sources within the agricultural regions of the Hawaiian island Kaua’i during climatological conditions spanning season conditions and the La Niña, El Niño, and neutral phases of the El Niño-Southern Oscillation. Both pollen size and source location had major effects on predicted dispersion over and near themore » island. Three patterns of pollen dispersion were identified in response to prevailing wind conditions: southwest winds transported pollen inland, funneling pollen grains through valleys; east winds transported pollen over the ocean, with dispersive tails for the smallest pollen grains following the mean wind and extending as far as the island of Ni’ihau 35 km away; and northeast winds moved pollen inland counter to the prevailing flow due to a sea breeze circulation that formed over the source region. These results are the first to predict the interactions between complex island terrain and local climatology on grass pollen dispersion. As a result, they demonstrate how numerical modeling can provide guidance for field trials by illustrating the common flow regimes present in complex terrain, allowing field trials to focus on areas where successful sampling is more likely to occur.« less

Gravity or turbulence? IV. Collapsing cores in out-of-virial disguise

NASA Astrophysics Data System (ADS)

Ballesteros-Paredes, Javier; Vázquez-Semadeni, Enrique; Palau, Aina; Klessen, Ralf S.

2018-06-01

We study the dynamical state of massive cores by using a simple analytical model, an observational sample, and numerical simulations of collapsing massive cores. From the analytical model, we find that cores increase their column density and velocity dispersion as they collapse, resulting in a time evolution path in the Larson velocity dispersion-size diagram from large sizes and small velocity dispersions to small sizes and large velocity dispersions, while they tend to equipartition between gravity and kinetic energy. From the observational sample, we find that: (a) cores with substantially different column densities in the sample do not follow a Larson-like linewidth-size relation. Instead, cores with higher column densities tend to be located in the upper-left corner of the Larson velocity dispersion σv, 3D-size R diagram, a result explained in the hierarchical and chaotic collapse scenario. (b) Cores appear to have overvirial values. Finally, our numerical simulations reproduce the behavior predicted by the analytical model and depicted in the observational sample: collapsing cores evolve towards larger velocity dispersions and smaller sizes as they collapse and increase their column density. More importantly, however, they exhibit overvirial states. This apparent excess is due to the assumption that the gravitational energy is given by the energy of an isolated homogeneous sphere. However, such excess disappears when the gravitational energy is correctly calculated from the actual spatial mass distribution. We conclude that the observed energy budget of cores is consistent with their non-thermal motions being driven by their self-gravity and in the process of dynamical collapse.

Composition measurement in substitutionally disordered materials by atomic resolution energy dispersive X-ray spectroscopy in scanning transmission electron microscopy.

PubMed

Chen, Z; Taplin, D J; Weyland, M; Allen, L J; Findlay, S D

2017-05-01

The increasing use of energy dispersive X-ray spectroscopy in atomic resolution scanning transmission electron microscopy invites the question of whether its success in precision composition determination at lower magnifications can be replicated in the atomic resolution regime. In this paper, we explore, through simulation, the prospects for composition measurement via the model system of Al x Ga 1-x As, discussing the approximations used in the modelling, the variability in the signal due to changes in configuration at constant composition, and the ability to distinguish between different compositions. Results are presented in such a way that the number of X-ray counts, and thus the expected variation due to counting statistics, can be gauged for a range of operating conditions. Copyright © 2016 Elsevier B.V. All rights reserved.

Simulating oil droplet dispersal from the Deepwater Horizon spill with a Lagrangian approach

USGS Publications Warehouse

North, Elizabeth W.; Schlag, Zachary; Adams, E. Eric; Sherwood, Christopher R.; He, Ruoying; Hyun, Hoon; Socolofsky, Scott A.

2011-01-01

An analytical multiphase plume model, combined with time-varying flow and hydrographic fields generated by the 3-D South Atlantic Bight and Gulf of Mexico model (SABGOM) hydrodynamic model, were used as input to a Lagrangian transport model (LTRANS), to simulate transport of oil droplets dispersed at depth from the recent Deepwater Horizon MC 252 oil spill. The plume model predicts a stratification-dominated near field, in which small oil droplets detrain from the central plume containing faster rising large oil droplets and gas bubbles and become trapped by density stratification. Simulated intrusion (trap) heights of ∼ 310–370 m agree well with the midrange of conductivity-temperature-depth observations, though the simulated variation in trap height was lower than observed, presumably in part due to unresolved variability in source composition (percentage oil versus gas) and location (multiple leaks during first half of spill). Simulated droplet trajectories by the SABGOM-LTRANS modeling system showed that droplets with diameters between 10 and 50 μm formed a distinct subsurface plume, which was transported horizontally and remained in the subsurface for >1 month. In contrast, droplets with diameters ≥90 μm rose rapidly to the surface. Simulated trajectories of droplets ≤50 μm in diameter were found to be consistent with field observations of a southwest-tending subsurface plume in late June 2010 reported by Camilli et al. [2010]. Model results suggest that the subsurface plume looped around to the east, with potential subsurface oil transport to the northeast and southeast. Ongoing work is focusing on adding degradation processes to the model to constrain droplet dispersal.

Simulating dispersal of reintroduced species within heterogeneous landscapes

Treesearch

Robert H. Gardner; Eric J. Gustafson

2004-01-01

This paper describes the development and application of a spatially explicit, individual based model of animal dispersal (J-walk) to determine the relative effects of landscape heterogeneity, prey availability, predation risk, and the energy requirements and behavior of dispersing organisms on dispersal success. Significant unknowns exist for the simulation of complex...

The problem of complex eigensystems in the semianalytical solution for advancement of time in solute transport simulations: a new method using real arithmetic

USGS Publications Warehouse

Umari, Amjad M.J.; Gorelick, Steven M.

1986-01-01

In the numerical modeling of groundwater solute transport, explicit solutions may be obtained for the concentration field at any future time without computing concentrations at intermediate times. The spatial variables are discretized and time is left continuous in the governing differential equation. These semianalytical solutions have been presented in the literature and involve the eigensystem of a coefficient matrix. This eigensystem may be complex (i.e., have imaginary components) due to the asymmetry created by the advection term in the governing advection-dispersion equation. Previous investigators have either used complex arithmetic to represent a complex eigensystem or chosen large dispersivity values for which the imaginary components of the complex eigenvalues may be ignored without significant error. It is shown here that the error due to ignoring the imaginary components of complex eigenvalues is large for small dispersivity values. A new algorithm that represents the complex eigensystem by converting it to a real eigensystem is presented. The method requires only real arithmetic.

Wetting-Dewetting and Dispersion-Aggregation Transitions Are Distinct for Polymer Grafted Nanoparticles in Chemically Dissimilar Polymer Matrix.

PubMed

Martin, Tyler B; Mongcopa, Katrina Irene S; Ashkar, Rana; Butler, Paul; Krishnamoorti, Ramanan; Jayaraman, Arthi

2015-08-26

Simulations and experiments are conducted on mixtures containing polymer grafted nanoparticles in a chemically distinct polymer matrix, where the graft and matrix polymers exhibit attractive enthalpic interactions at low temperatures that become progressively repulsive as temperature is increased. Both coarse-grained molecular dynamics simulations, and X-ray scattering and neutron scattering experiments with deuterated polystyrene (dPS) grafted silica and poly(vinyl methyl ether) PVME matrix show that the sharp phase transition from (mixed) dispersed to (demixed) aggregated morphologies due to the increasingly repulsive effective interactions between the blend components is distinct from the continuous wetting-dewetting transition. Strikingly, this is unlike the extensively studied chemically identical graft-matrix composites, where the two transitions have been considered to be synonymous, and is also unlike the free (ungrafted) blends of the same graft and matrix homopolymers, where the wetting-dewetting is a sharp transition coinciding with the macrophase separation.

Simulating the effects of upstream turbulence on dispersion around a building

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhang, Y.Q.; Arya, S.P.S.; Huber, A.H.

The effects of high turbulence versus no turbulence in a sheared boundary-layer flow approaching a building are being investigated by a turbulent kinetic energy/dissipation model (TEMPEST). The effects on both the mean flow and the concentration field around a cubical building are presented. The numerical simulations demonstrate significant effects due to the differences in the incident flow. The addition of upstream turbulence results in a reduced size of the cavity directly behind the building. The velocity deficits in the wake strongly depend on the upstream turbulence intensities. The accuracy of numerical simulations is verified by comparing the predicted mean flowmore » and concentration fields with the wind tunnel measurements of Castro and Robins (1977) and Robins and Castro (1977, 1975). Comparing the results with experimental data, the authors show that the TEMPEST model can reasonably simulate the mean flow. The numerical simulations of the concentration fields due to a source on the roof-top of the building are presented. Both the value and the position of the maximum ground-level concentration are changed dramatically due to the effects of the upstream level of turblence.« less

Energy Budget of Forming Clumps in Numerical Simulations of Collapsing Clouds

NASA Astrophysics Data System (ADS)

Camacho, Vianey; Vázquez-Semadeni, Enrique; Ballesteros-Paredes, Javier; Gómez, Gilberto C.; Fall, S. Michael; Mata-Chávez, M. Dolores

2016-12-01

We analyze the physical properties and energy balance of density enhancements in two SPH simulations of the formation, evolution, and collapse of giant molecular clouds. In the simulations, no feedback is included, so all motions are due either to the initial decaying turbulence or to gravitational contraction. We define clumps as connected regions above a series of density thresholds. The resulting full set of clumps follows the generalized energy equipartition relation, {σ }v/{R}1/2\\propto {{{Σ }}}1/2, where {σ }v is the velocity dispersion, R is the “radius,” and Σ is the column density. We interpret this as a natural consequence of gravitational contraction at all scales rather than virial equilibrium. Nevertheless, clumps with low Σ tend to show a large scatter around equipartition. In more than half of the cases, this scatter is dominated by external turbulent compressions that assemble the clumps rather than by small-scale random motions that would disperse them. The other half does actually disperse. Moreover, clump sub-samples selected by means of different criteria exhibit different scalings. Sub-samples with narrow Σ ranges follow Larson-like relations, although characterized by their respective values of Σ. Finally, we find that (I) clumps lying in filaments tend to appear sub-virial, (II) high-density cores (n≥slant {10}5 cm3) that exhibit moderate kinetic energy excesses often contain sink (“stellar”) particles and the excess disappears when the stellar mass is taken into account in the energy balance, and (III) cores with kinetic energy excess but no stellar particles are truly in a state of dispersal.

Flow and axial dispersion in a sinusoidal-walled tube: Effects of inertial and unsteady flows

NASA Astrophysics Data System (ADS)

Richmond, Marshall C.; Perkins, William A.; Scheibe, Timothy D.; Lambert, Adam; Wood, Brian D.

2013-12-01

In this work, we consider a sinusoidal-walled tube (a three-dimensional tube with sinusoidally-varying diameter) as a simplified conceptualization of flow in porous media. Direct numerical simulation using computational fluid dynamics (CFD) methods was used to compute velocity fields by solving the Navier-Stokes equations, and also to numerically solve the volume averaging closure problem, for a range of Reynolds numbers (Re) spanning the low-Re to inertial flow regimes, including one simulation at Re=449 for which unsteady flow was observed. The longitudinal dispersion observed for the flow was computed using a random walk particle tracking method, and this was compared to the longitudinal dispersion predicted from a volume-averaged macroscopic mass balance using the method of volume averaging; the results of the two methods were consistent. Our results are compared to experimental measurements of dispersion in porous media and to previous theoretical results for both the low-Re, Stokes flow regime and for values of Re representing the steady inertial regime. In the steady inertial regime, a power-law increase in the effective longitudinal dispersion (DL) with Re was found, and this is consistent with previous results. This rapid rate of increase is caused by trapping of solute in expansions due to flow separation (eddies). One unsteady (but non-turbulent) flow case (Re=449) was also examined. For this case, the rate of increase of DL with Re was smaller than that observed at lower Re. Velocity fluctuations in this regime lead to increased rates of solute mass transfer between the core flow and separated flow regions, thus diminishing the amount of tailing caused by solute trapping in eddies and thereby reducing longitudinal dispersion. The observed tailing was further explored through analysis of concentration skewness (third moment) and its assymptotic convergence to conventional advection-dispersion behavior (skewness = 0). The method of volume averaging was applied to develop a skewness model, and demonstrated that the skewness decreases as a function of inverse square root of time. Our particle tracking simulation results were shown to conform to this theoretical result in most of the cases considered.

Evaluation of a Eulerian and Lagrangian air quality model using perfluorocarbon tracers released in Texas for the BRAVO haze study

NASA Astrophysics Data System (ADS)

Schichtel, Bret A.; Barna, Michael G.; Gebhart, Kristi A.; Malm, William C.

The Big Bend Regional Aerosol and Visibility Observational (BRAVO) study was designed to determine the sources of haze at Big Bend National Park, Texas, using a combination of source and receptor models. BRAVO included an intensive monitoring campaign from July to October 1999 that included the release of perfluorocarbon tracers from four locations at distances 230-750 km from Big Bend and measured at 24 sites. The tracer measurements near Big Bend were used to evaluate the dispersion mechanisms in the REMSAD Eulerian model and the CAPITA Monte Carlo (CMC) Lagrangian model used in BRAVO. Both models used 36 km MM5 wind fields as input. The CMC model also used a combination of routinely available 80 and 190 km wind fields from the National Weather Service's National Centers for Environmental Prediction (NCEP) as input. A model's performance is limited by inherent uncertainties due to errors in the tracer concentrations and a model's inability to simulate sub-resolution variability. A range in the inherent uncertainty was estimated by comparing tracer data at nearby monitoring sites. It was found that the REMSAD and CMC models, using the MM5 wind field, produced performance statistics generally within this inherent uncertainty. The CMC simulation using the NCEP wind fields could reproduce the timing of tracer impacts at Big Bend, but not the concentration values, due to a systematic underestimation. It appears that the underestimation was partly due to excessive vertical dilution from high mixing depths. The model simulations were more sensitive to the input wind fields than the models' different dispersion mechanisms. Comparisons of REMSAD to CMC tracer simulations using the MM5 wind fields had correlations between 0.75 and 0.82, depending on the tracer, but the tracer simulations using the two wind fields in the CMC model had correlations between 0.37 and 0.5.

Simulated distribution and ecotoxicity-based assessment of chemically-dispersed oil in Tokyo Bay.

PubMed

Koyama, Jiro; Imakado, Chie; Uno, Seiichi; Kuroda, Takako; Hara, Shouichi; Majima, Takahiro; Shirota, Hideyuki; Añasco, Nathaniel C

2014-08-30

To assess risks of chemically-dispersed oil to marine organisms, oil concentrations in the water were simulated using a hypothetical spill accident in Tokyo Bay. Simulated oil concentrations were then compared with the short-term no-observed effect concentration (NOEC), 0.01 mg/L, obtained through toxicity tests using marine diatoms, amphipod and fish. Area of oil concentrations higher than the NOEC were compared with respect to use and non-use of dispersant. Results of the simulation show relatively faster dispersion near the mouth of the bay compared to its inner sections which is basically related to its stronger water currents. Interestingly, in the inner bay, a large area of chemically-dispersed oil has concentrations higher than the NOEC. It seems emulsifying oil by dispersant increases oil concentrations, which could lead to higher toxicity to aquatic organisms. When stronger winds occur, however, the difference in toxic areas between use and non-use of dispersant is quite small. Copyright © 2014 Elsevier Ltd. All rights reserved.

Amorphous Solid Dispersion of Epigallocatechin Gallate for Enhanced Physical Stability and Controlled Release.

PubMed

Cao, Yizheng; Teng, Jing; Selbo, Jon

2017-11-09

Epigallocatechin gallate (EGCG) has been recognized as the most prominent green tea extract due to its healthy influences. The high instability and low bioavailability, however, strongly limit its utilization in food and drug industries. This work, for the first time, develops amorphous solid dispersion of EGCG to enhance its bioavailability and physical stability. Four commonly used polymeric excipients are found to be compatible with EGCG in water-dioxane mixtures via a stepwise mixing method aided by vigorous mechanical interference. The dispersions are successfully generated by lyophilization. The physical stability of the dispersions is significantly improved compared to pure amorphous EGCG in stress condition (elevated temperature and relative humidity) and simulated gastrointestinal tract environment. From the drug release tests, one of the dispersions, EGCG-Soluplus ® 50:50 ( w / w ) shows a dissolution profile that only 50% EGCG is released in the first 20 min, and the remains are slowly released in 24 h. This sustained release profile may open up new possibilities to increase EGCG bioavailability via extending its elimination time in plasma.

Amorphous Solid Dispersion of Epigallocatechin Gallate for Enhanced Physical Stability and Controlled Release

PubMed Central

Cao, Yizheng; Teng, Jing; Selbo, Jon

2017-01-01

Epigallocatechin gallate (EGCG) has been recognized as the most prominent green tea extract due to its healthy influences. The high instability and low bioavailability, however, strongly limit its utilization in food and drug industries. This work, for the first time, develops amorphous solid dispersion of EGCG to enhance its bioavailability and physical stability. Four commonly used polymeric excipients are found to be compatible with EGCG in water-dioxane mixtures via a stepwise mixing method aided by vigorous mechanical interference. The dispersions are successfully generated by lyophilization. The physical stability of the dispersions is significantly improved compared to pure amorphous EGCG in stress condition (elevated temperature and relative humidity) and simulated gastrointestinal tract environment. From the drug release tests, one of the dispersions, EGCG-Soluplus® 50:50 (w/w) shows a dissolution profile that only 50% EGCG is released in the first 20 min, and the remains are slowly released in 24 h. This sustained release profile may open up new possibilities to increase EGCG bioavailability via extending its elimination time in plasma. PMID:29120370

Simulation of the Processes of Formation of a Dust Cloud in a Vacuum and in the Absence of Gravitation

NASA Astrophysics Data System (ADS)

Avdeev, A. V.; Boreisho, A. S.; Ivakin, S. V.; Moiseev, A. A.; Savin, A. V.; Sokolov, E. I.; Smirnov, P. G.

2018-01-01

This article is devoted to the simulation of the processes of formation of dust clouds in the absence of gravitation, which is necessary for understanding the processes proceeding in dust clusters in outer space, upper planetary atmosphere, and on the surface of space objects, as well as for evaluating the possibilities of creating disperse structures with given properties. The chief aim of the simulation is to determine the general laws of the dynamics of the dust cloud at the initial stage of its formation. With the use of the original approach based on the particle-in-cell method that permits investigating the mechanics of large ensembles of particles on contemporary computational platforms, we consider the mechanics of a dusty medium in the process of its excitation in a closed container due to the vibration of the walls, and then in the process of particle scattering when the container opens in outer space. The main formation mechanisms of a dust formation have been elucidated, and the possibilities of mathematical simulation for predicting spatial and time characteristics of disperse structures have been shown.

Simulation of windblown dust transport from a mine tailings impoundment using a computational fluid dynamics model

NASA Astrophysics Data System (ADS)

Stovern, Michael; Felix, Omar; Csavina, Janae; Rine, Kyle P.; Russell, MacKenzie R.; Jones, Robert M.; King, Matt; Betterton, Eric A.; Sáez, A. Eduardo

2014-09-01

Mining operations are potential sources of airborne particulate metal and metalloid contaminants through both direct smelter emissions and wind erosion of mine tailings. The warmer, drier conditions predicted for the Southwestern US by climate models may make contaminated atmospheric dust and aerosols increasingly important, due to potential deleterious effects on human health and ecology. Dust emissions and dispersion of dust and aerosol from the Iron King Mine tailings in Dewey-Humboldt, Arizona, a Superfund site, are currently being investigated through in situ field measurements and computational fluid dynamics modeling. These tailings are heavily contaminated with lead and arsenic. Using a computational fluid dynamics model, we model dust transport from the mine tailings to the surrounding region. The model includes gaseous plume dispersion to simulate the transport of the fine aerosols, while individual particle transport is used to track the trajectories of larger particles and to monitor their deposition locations. In order to improve the accuracy of the dust transport simulations, both regional topographical features and local weather patterns have been incorporated into the model simulations. Results show that local topography and wind velocity profiles are the major factors that control deposition.

Simulation of windblown dust transport from a mine tailings impoundment using a computational fluid dynamics model.

PubMed

Stovern, Michael; Felix, Omar; Csavina, Janae; Rine, Kyle P; Russell, MacKenzie R; Jones, Robert M; King, Matt; Betterton, Eric A; Sáez, A Eduardo

2014-09-01

Mining operations are potential sources of airborne particulate metal and metalloid contaminants through both direct smelter emissions and wind erosion of mine tailings. The warmer, drier conditions predicted for the Southwestern US by climate models may make contaminated atmospheric dust and aerosols increasingly important, due to potential deleterious effects on human health and ecology. Dust emissions and dispersion of dust and aerosol from the Iron King Mine tailings in Dewey-Humboldt, Arizona, a Superfund site, are currently being investigated through in situ field measurements and computational fluid dynamics modeling. These tailings are heavily contaminated with lead and arsenic. Using a computational fluid dynamics model, we model dust transport from the mine tailings to the surrounding region. The model includes gaseous plume dispersion to simulate the transport of the fine aerosols, while individual particle transport is used to track the trajectories of larger particles and to monitor their deposition locations. In order to improve the accuracy of the dust transport simulations, both regional topographical features and local weather patterns have been incorporated into the model simulations. Results show that local topography and wind velocity profiles are the major factors that control deposition.

Simulation of windblown dust transport from a mine tailings impoundment using a computational fluid dynamics model

PubMed Central

Stovern, Michael; Felix, Omar; Csavina, Janae; Rine, Kyle P.; Russell, MacKenzie R.; Jones, Robert M.; King, Matt; Betterton, Eric A.; Sáez, A. Eduardo

2014-01-01

Mining operations are potential sources of airborne particulate metal and metalloid contaminants through both direct smelter emissions and wind erosion of mine tailings. The warmer, drier conditions predicted for the Southwestern US by climate models may make contaminated atmospheric dust and aerosols increasingly important, due to potential deleterious effects on human health and ecology. Dust emissions and dispersion of dust and aerosol from the Iron King Mine tailings in Dewey-Humboldt, Arizona, a Superfund site, are currently being investigated through in situ field measurements and computational fluid dynamics modeling. These tailings are heavily contaminated with lead and arsenic. Using a computational fluid dynamics model, we model dust transport from the mine tailings to the surrounding region. The model includes gaseous plume dispersion to simulate the transport of the fine aerosols, while individual particle transport is used to track the trajectories of larger particles and to monitor their deposition locations. In order to improve the accuracy of the dust transport simulations, both regional topographical features and local weather patterns have been incorporated into the model simulations. Results show that local topography and wind velocity profiles are the major factors that control deposition. PMID:25621085

Ensemble Simulation of the Atmospheric Radionuclides Discharged by the Fukushima Nuclear Accident

NASA Astrophysics Data System (ADS)

Sekiyama, Thomas; Kajino, Mizuo; Kunii, Masaru

2013-04-01

Enormous amounts of radionuclides were discharged into the atmosphere by a nuclear accident at the Fukushima Daiichi nuclear power plant (FDNPP) after the earthquake and tsunami on 11 March 2011. The radionuclides were dispersed from the power plant and deposited mainly over eastern Japan and the North Pacific Ocean. A lot of numerical simulations of the radionuclide dispersion and deposition had been attempted repeatedly since the nuclear accident. However, none of them were able to perfectly simulate the distribution of dose rates observed after the accident over eastern Japan. This was partly due to the error of the wind vectors and precipitations used in the numerical simulations; unfortunately, their deterministic simulations could not deal with the probability distribution of the simulation results and errors. Therefore, an ensemble simulation of the atmospheric radionuclides was performed using the ensemble Kalman filter (EnKF) data assimilation system coupled with the Japan Meteorological Agency (JMA) non-hydrostatic mesoscale model (NHM); this mesoscale model has been used operationally for daily weather forecasts by JMA. Meteorological observations were provided to the EnKF data assimilation system from the JMA operational-weather-forecast dataset. Through this ensemble data assimilation, twenty members of the meteorological analysis over eastern Japan from 11 to 31 March 2011 were successfully obtained. Using these meteorological ensemble analysis members, the radionuclide behavior in the atmosphere such as advection, convection, diffusion, dry deposition, and wet deposition was simulated. This ensemble simulation provided the multiple results of the radionuclide dispersion and distribution. Because a large ensemble deviation indicates the low accuracy of the numerical simulation, the probabilistic information is obtainable from the ensemble simulation results. For example, the uncertainty of precipitation triggered the uncertainty of wet deposition; the uncertainty of wet deposition triggered the uncertainty of atmospheric radionuclide amounts. Then the remained radionuclides were transported downwind; consequently the uncertainty signal of the radionuclide amounts was propagated downwind. The signal propagation was seen in the ensemble simulation by the tracking of the large deviation areas of radionuclide concentration and deposition. These statistics are able to provide information useful for the probabilistic prediction of radionuclides.

PIC simulations of a three component plasma described by Kappa distribution functions as observed in Saturn's magnetosphere

NASA Astrophysics Data System (ADS)

Barbosa, Marcos; Alves, Maria Virginia; Simões Junior, Fernando

2016-04-01

In plasmas out of thermodynamic equilibrium the particle velocity distribution can be described by the so called Kappa distribution. These velocity distribution functions are a generalization of the Maxwellian distribution. Since 1960, Kappa velocity distributions were observed in several regions of interplanetary space and astrophysical plasmas. Using KEMPO1 particle simulation code, modified to introduce Kappa distribution functions as initial conditions for particle velocities, the normal modes of propagation were analyzed in a plasma containing two species of electrons with different temperatures and densities and ions as a third specie.This type of plasma is usually found in magnetospheres such as in Saturn. Numerical solutions for the dispersion relation for such a plasma predict the presence of an electron-acoustic mode, besides the Langmuir and ion-acoustic modes. In the presence of an ambient magnetic field, the perpendicular propagation (Bernstein mode) also changes, as compared to a Maxwellian plasma, due to the Kappa distribution function. Here results for simulations with and without external magnetic field are presented. The parameters for the initial conditions in the simulations were obtained from the Cassini spacecraft data. Simulation results are compared with numerical solutions of the dispersion relation obtained in the literature and they are in good agreement.

Enhanced Bioaccessibility of Curcuminoids in Buttermilk Yogurt in Comparison to Curcuminoids in Aqueous Dispersions.

PubMed

Fu, Shishan; Augustin, Mary Ann; Sanguansri, Luz; Shen, Zhiping; Ng, Ken; Ajlouni, Said

2016-03-01

Curcuminoids have low bioavailability due to low aqueous solubility. We compared the bioaccessibility of curcuminoids delivered in buttermilk yogurt to that of curcuminoid powder in an aqueous dispersion. Buttermilk containing added curcuminoids (300 mg/100 g, 0.3% w/w) was used for yogurt manufacture. We measured percentage of curcuminoids remaining in yogurts after manufacture and after exposure to simulated gastrointestinal fluids, and the in vitro bioaccessibility of the curcuminoids. Curcuminoids were stable during yogurt manufacture. At the end of in vitro digestion, approximately 11% of the curcuminoids delivered in yogurt was degraded compared to <1% for curcuminoids in an aqueous dispersion. However, curcuminoids delivered in yogurt was 15-fold more bioaccessible than curcuminoids in aqueous dispersion. The small change in yogurt properties (decrease in total lactic acid bacteria counts of <1 log and increased viscosity) on addition of curcuminoids has to be balanced against the benefits of increased bioaccessibility of curcuminoids when delivered in yogurts. © 2016 Institute of Food Technologists®

Oceanic dispersion of Fukushima-derived Cs-137 in the coastal, offshore, and open oceans simulated by multiple oceanic general circulation models

NASA Astrophysics Data System (ADS)

Kawamura, H.; Furuno, A.; Kobayashi, T.; In, T.; Nakayama, T.; Ishikawa, Y.; Miyazawa, Y.; Usui, N.

2017-12-01

To understand the concentration and amount of Fukushima-derived Cs-137 in the ocean, this study simulates the oceanic dispersion of Cs-137 by an oceanic dispersion model SEA-GEARN-FDM developed at Japan Atomic Energy Agency (JAEA) and multiple oceanic general circulation models. The Cs-137 deposition amounts at the sea surface were used as the source term in oceanic dispersion simulations, which were estimated by atmospheric dispersion simulations with a Worldwide version of System for Prediction of Environmental Emergency Dose Information version II (WSPEEDI-II) developed at JAEA. The direct release from the Fukushima Daiichi Nuclear Power Plant into the ocean based on in situ Cs-137 measurements was used as the other source term in oceanic dispersion simulations. The simulated air Cs-137 concentrations qualitatively replicated those measured around the North Pacific. The accumulated Cs-137 ground deposition amount in the eastern Japanese Islands was consistent with that estimated by aircraft measurements. The oceanic dispersion simulations relatively well reproduced the measured Cs-137 concentrations in the coastal and offshore oceans during the first few months after the Fukushima disaster, and in the open ocean during the first year post-disaster. It was suggested that Cs-137 dispersed along the coast in the north-south direction during the first few months post-disaster, and were subsequently dispersed offshore by the Kuroshio Current and Kuroshio Extension. Mesoscale eddies accompanied by the Kuroshio Current and Kuroshio Extension played an important role in dilution of Cs-137. The Cs-137 amounts were quantified in the coastal, offshore, and open oceans during the first year post-disaster. It was demonstrated that Cs-137 actively dispersed from the coastal and offshore oceans to the open ocean, and from the surface layer to the deeper layer in the North Pacific.

Chromatography modelling to describe protein adsorption at bead level.

PubMed

Gerontas, Spyridon; Shapiro, Michael S; Bracewell, Daniel G

2013-04-05

Chromatographic modelling can be used to describe and further understand the behaviour of biological species during their chromatography separation on adsorption resins. Current modelling approaches assume uniform rate parameters throughout the column. Software and hardware advances now allow us to consider what can be learnt from modelling at bead level, enabling simulation of heterogeneity in bead and packed bed structure due to design or due to changes during operation. In this paper, a model has been developed to simulate at bead level protein loading in 1.5 μl microfluidic columns. This model takes into account the heterogeneity in bead sizes and the spatial variations of the characteristics of a packed bed, such as bed void fraction and dispersion, thus offering a detailed description of the flow field and mass transfer phenomena. Simulations were shown to be in good agreement with published experimental data. Copyright © 2013 Elsevier B.V. All rights reserved.

Material parameters from frequency dispersion simulation of floating gate memory with Ge nanocrystals in HfO2

NASA Astrophysics Data System (ADS)

Palade, C.; Lepadatu, A. M.; Slav, A.; Lazanu, S.; Teodorescu, V. S.; Stoica, T.; Ciurea, M. L.

2018-01-01

Trilayer memory capacitors with Ge nanocrystals (NCs) floating gate in HfO2 were obtained by magnetron sputtering deposition on p-type Si substrate followed by rapid thermal annealing at relatively low temperature of 600 °C. The frequency dispersion of capacitance and resistance was measured in accumulation regime of Al/HfO2 gate oxide/Ge NCs in HfO2 floating gate/HfO2 tunnel oxide/SiOx/p-Si/Al memory capacitors. For simulation of the frequency dispersion a complex circuit model was used considering an equivalent parallel RC circuit for each layer of the trilayer structure. A series resistance due to metallic contacts and Si substrate was necessary to be included in the model. A very good fit to the experimental data was obtained and the parameters of each layer in the memory capacitor, i.e. capacitances and resistances were determined and in turn the intrinsic material parameters, i.e. dielectric constants and resistivities of layers were evaluated. The results are very important for the study and optimization of the hysteresis behaviour of floating gate memories based on NCs embedded in oxide.

Explosive particle soil surface dispersion model for detonated military munitions.

PubMed

Hathaway, John E; Rishel, Jeremy P; Walsh, Marianne E; Walsh, Michael R; Taylor, Susan

2015-07-01

The accumulation of high explosive mass residue from the detonation of military munitions on training ranges is of environmental concern because of its potential to contaminate the soil, surface water, and groundwater. The US Department of Defense wants to quantify, understand, and remediate high explosive mass residue loadings that might be observed on active firing ranges. Previously, efforts using various sampling methods and techniques have resulted in limited success, due in part to the complicated dispersion pattern of the explosive particle residues upon detonation. In our efforts to simulate particle dispersal for high- and low-order explosions on hypothetical firing ranges, we use experimental particle data from detonations of munitions from a 155-mm howitzer, which are common military munitions. The mass loadings resulting from these simulations provide a previously unattained level of detail to quantify the explosive residue source-term for use in soil and water transport models. In addition, the resulting particle placements can be used to test, validate, and optimize particle sampling methods and statistical models as applied to firing ranges. Although the presented results are for a hypothetical 155-mm howitzer firing range, the method can be used for other munition types once the explosive particle characteristics are known.

An analytic description of electrodynamic dispersion in free-flow zone electrophoresis.

PubMed

Dutta, Debashis

2015-07-24

The present work analyzes the electrodynamic dispersion of sample streams in a free-flow zone electrophoresis (FFZE) chamber resulting due to partial or complete blockage of electroosmotic flow (EOF) across the channel width by the sidewalls of the conduit. This blockage of EOF has been assumed to generate a pressure-driven backflow in the transverse direction for maintaining flow balance in the system. A parallel-plate based FFZE device with the analyte stream located far away from the channel side regions has been considered to simplify the current analysis. Applying a method-of-moments formulation, an analytic expression was derived for the variance of the sample zone at steady state as a function of its position in the separation chamber under these conditions. It has been shown that the increase in stream broadening due to the electrodynamic dispersion phenomenon is additive to the contributions from molecular diffusion and sample injection, and simply modifies the coefficient for the hydrodynamic dispersion term for a fixed lateral migration distance of the sample stream. Moreover, this dispersion mechanism can dominate the overall spatial variance of analyte zones when a significant fraction of the EOF is blocked by the channel sidewalls. The analysis also shows that analyte streams do not undergo any hydrodynamic broadening due to unwanted pressure-driven cross-flows in an FFZE chamber in the absence of a transverse electric field. The noted results have been validated using Monte Carlo simulations which further demonstrate that while the sample concentration profile at the channel outlet approaches a Gaussian distribution only in FFZE chambers substantially longer than the product of the axial pressure-driven velocity and the characteristic diffusion time in the system, the spatial variance of the exiting analyte stream is well described by the Taylor-Aris dispersion limit even in analysis ducts much shorter than this length scale. Copyright © 2015 Elsevier B.V. All rights reserved.

How are quasars fueled? Simulating interstellar gas in tidally disturbed galaxies

NASA Technical Reports Server (NTRS)

Byrd, Gene G.

1986-01-01

Whether gravitational tides from companions trigger global instabilities in spiral galaxy disks and thus rapid flows of gas into the nucleus to fuel activity is investigated. An n-body computer program is used to simulate the disk of the spiral galaxy within a much more stable, high-velocity dispersion spherical halo. Under sufficient perturbation, the disk undergoes violent distortions due to the disturber and its self-gravitation. The tidal action of companions was simulated and the tidal strengths at which the instabilities appear to match those of the observed companions of Seyferts and quasars was shown. With the additional modifications planned, the gas flow will be more realistically simulated to compare with observations (e.g., colors, velocity fields) of active galaxies.

Road traffic impact on urban water quality: a step towards integrated traffic, air and stormwater modelling.

PubMed

Fallah Shorshani, Masoud; Bonhomme, Céline; Petrucci, Guido; André, Michel; Seigneur, Christian

2014-04-01

Methods for simulating air pollution due to road traffic and the associated effects on stormwater runoff quality in an urban environment are examined with particular emphasis on the integration of the various simulation models into a consistent modelling chain. To that end, the models for traffic, pollutant emissions, atmospheric dispersion and deposition, and stormwater contamination are reviewed. The present study focuses on the implementation of a modelling chain for an actual urban case study, which is the contamination of water runoff by cadmium (Cd), lead (Pb), and zinc (Zn) in the Grigny urban catchment near Paris, France. First, traffic emissions are calculated with traffic inputs using the COPERT4 methodology. Next, the atmospheric dispersion of pollutants is simulated with the Polyphemus line source model and pollutant deposition fluxes in different subcatchment areas are calculated. Finally, the SWMM water quantity and quality model is used to estimate the concentrations of pollutants in stormwater runoff. The simulation results are compared to mass flow rates and concentrations of Cd, Pb and Zn measured at the catchment outlet. The contribution of local traffic to stormwater contamination is estimated to be significant for Pb and, to a lesser extent, for Zn and Cd; however, Pb is most likely overestimated due to outdated emissions factors. The results demonstrate the importance of treating distributed traffic emissions from major roadways explicitly since the impact of these sources on concentrations in the catchment outlet is underestimated when those traffic emissions are spatially averaged over the catchment area.

Assessment of risk due to the use of carbon fiber composites in commercial and general aviation

NASA Technical Reports Server (NTRS)

Fiksel, J.; Rosenfield, D.; Kalelkar, A.

1980-01-01

The development of a national risk profile for the total annual aircraft losses due to carbon fiber composite (CFC) usage through 1993 is discussed. The profile was developed using separate simulation methods for commercial and general aviation aircraft. A Monte Carlo method which was used to assess the risk in commercial aircraft is described. The method projects the potential usage of CFC through 1993, investigates the incidence of commercial aircraft fires, models the potential release and dispersion of carbon fibers from a fire, and estimates potential economic losses due to CFC damaging electronic equipment. The simulation model for the general aviation aircraft is described. The model emphasizes variations in facility locations and release conditions, estimates distribution of CFC released in general aviation aircraft accidents, and tabulates the failure probabilities and aggregate economic losses in the accidents.

Oil viscosity limitation on dispersibility of crude oil under simulated at-sea conditions in a large wave tank.

PubMed

Trudel, Ken; Belore, Randy C; Mullin, Joseph V; Guarino, Alan

2010-09-01

This study determined the limiting oil viscosity for chemical dispersion of oil spills under simulated sea conditions in the large outdoor wave tank at the US National Oil Spill Response Test Facility in New Jersey. Dispersant effectiveness tests were completed using crude oils with viscosities ranging from 67 to 40,100 cP at test temperature. Tests produced an effectiveness-viscosity curve with three phases when oil was treated with Corexit 9500 at a dispersant-to-oil ratio of 1:20. The oil viscosity that limited chemical dispersion under simulated at-sea conditions was in the range of 18,690 cP to 33,400 cP. Visual observations and measurements of oil concentrations and droplet size distributions in the water under treated and control slicks correlated well with direct measurements of effectiveness. The dispersant effectiveness versus oil viscosity relationship under simulated at sea conditions at Ohmsett was most similar to those from similar tests made using the Institut Francais du Pétrole and Exxon Dispersant Effectiveness (EXDET) test methods. Copyright 2010 Elsevier Ltd. All rights reserved.

Kinetic theory-based numerical modeling and analysis of bi-disperse segregated mixture fluidized bed

DOE Office of Scientific and Technical Information (OSTI.GOV)

Konan, N. A.; Huckaby, E. D.

We discuss a series of continuum Euler-Euler simulations of an initially mixed bi-disperse fluidized bed which segregates under certain operating conditions. The simulations use the multi-phase kinetic theory-based description of the momentum and energy exchanges between the phases by Simonin’s Group [see e.g. Gourdel, Simonin and Brunier (1999). Proceedings of 6th International Conference on Circulating Fluidized Beds, Germany, pp. 205-210]. The discussion and analysis of the results focus on the fluid-particle momentum exchange (i.e. drag). Simulations using mono- and poly-disperse fluid-particle drag correlations are analyzed for the Geldart D-type size bi-disperse gas-solid experiments performed by Goldschmidt et al. [Powder Tech.,more » pp. 135-159 (2003)]. The poly-disperse gas-particle drag correlations account for the local particle size distribution by using an effective mixture diameter when calculating the Reynolds number and then correcting the resulting force coefficient. Simulation results show very good predictions of the segregation index for bidisperse beds with the mono-disperse drag correlations contrary to the poly-disperse drag correlations for which the segregation rate is systematically under-predicted. The statistical analysis of the results shows a clear separation in the distribution of the gas-particle mean relaxation times of the small and large particles with simulations using the mono-disperse drag. In contrast, the poly-disperse drag simulations have a significant overlap and also a smaller difference in the mean particle relaxation times. This results in the small and large particles in the bed to respond to the gas similarly without enough relative time lag. The results suggest that the difference in the particle response time induce flow dynamics favorable to a force imbalance which results in the segregation.« less

Kinetic theory-based numerical modeling and analysis of bi-disperse segregated mixture fluidized bed

DOE PAGES

Konan, N. A.; Huckaby, E. D.

2017-06-21

We discuss a series of continuum Euler-Euler simulations of an initially mixed bi-disperse fluidized bed which segregates under certain operating conditions. The simulations use the multi-phase kinetic theory-based description of the momentum and energy exchanges between the phases by Simonin’s Group [see e.g. Gourdel, Simonin and Brunier (1999). Proceedings of 6th International Conference on Circulating Fluidized Beds, Germany, pp. 205-210]. The discussion and analysis of the results focus on the fluid-particle momentum exchange (i.e. drag). Simulations using mono- and poly-disperse fluid-particle drag correlations are analyzed for the Geldart D-type size bi-disperse gas-solid experiments performed by Goldschmidt et al. [Powder Tech.,more » pp. 135-159 (2003)]. The poly-disperse gas-particle drag correlations account for the local particle size distribution by using an effective mixture diameter when calculating the Reynolds number and then correcting the resulting force coefficient. Simulation results show very good predictions of the segregation index for bidisperse beds with the mono-disperse drag correlations contrary to the poly-disperse drag correlations for which the segregation rate is systematically under-predicted. The statistical analysis of the results shows a clear separation in the distribution of the gas-particle mean relaxation times of the small and large particles with simulations using the mono-disperse drag. In contrast, the poly-disperse drag simulations have a significant overlap and also a smaller difference in the mean particle relaxation times. This results in the small and large particles in the bed to respond to the gas similarly without enough relative time lag. The results suggest that the difference in the particle response time induce flow dynamics favorable to a force imbalance which results in the segregation.« less

Attenuation, dispersion and nonlinearity effects in graphene-based waveguides

PubMed Central

Mota, João Cesar Moura; Sombra, Antonio Sergio Bezerra

2015-01-01

Summary We simulated and analyzed in detail the behavior of ultrashort optical pulses, which are typically used in telecommunications, propagating through graphene-based nanoribbon waveguides. In this work, we showed the changes that occur in the Gaussian and hyperbolic secant input pulses due to the attenuation, high-order dispersive effects and nonlinear effects. We concluded that it is possible to control the shape of the output pulses with the value of the input signal power and the chemical potential of the graphene nanoribbon. We believe that the obtained results will be highly relevant since they can be applied to other nanophotonic devices, for example, filters, modulators, antennas, switches and other devices. PMID:26171299

Development of an analytical Lagrangian model for passive scalar dispersion in low-wind speed meandering conditions

NASA Astrophysics Data System (ADS)

Stefanello, M. B.; Degrazia, G. A.; Mortarini, L.; Buligon, L.; Maldaner, S.; Carvalho, J. C.; Acevedo, O. C.; Martins, L. G. N.; Anfossi, D.; Buriol, C.; Roberti, D.

2018-02-01

Describing the effects of wind meandering motions on the dispersion of scalars is a challenging task, since this type of flow represents a physical state characterized by multiple scales. In this study, a Lagrangian stochastic diffusion model is derived to describe scalar transport during the horizontal wind meandering phenomenon that occurs within a planetary boundary layer. The model is derived from the linearization of the Langevin equation, and it employs a heuristic functional form that represents the autocorrelation function of meandering motion. The new solutions, which describe the longitudinal and lateral wind components, were used to simulate tracer experiments that were performed in low-wind speed conditions. The results of the comparison indicate that the new model can effectively reproduce the observed concentrations of the contaminants, and therefore, it can satisfactorily describe enhanced dispersion effects due to the presence of meandering.

SIRANERISK: Modelling dispersion of steady and unsteady pollutant releases in the urban canopy

NASA Astrophysics Data System (ADS)

Soulhac, L.; Lamaison, G.; Cierco, F.-X.; Ben Salem, N.; Salizzoni, P.; Mejean, P.; Armand, P.; Patryl, L.

2016-09-01

SIRANERISK is an operational model for the simulation of the dispersion of unsteady atmospheric releases of pollutant within and above an urban area. SIRANERISK is built on the same principles as the SIRANE model, and couples a street network model for the pollutant transfers within the urban canopy with a Gaussian puff model for the transfers above it. The performance of the model are here analysed by a detailed comparisons with wind-tunnel experiments. These experiments concern the dispersion of steady and unsteady pollutant releases within and above obstacle arrays with varying geometrical configurations, representing different topologies of idealised urban districts. The overall good agreement between numerical and experimental data demonstrates the reliability of SIRANERISK as an operational tool for the assessment of risk analysis and for the management of crises due to the accidental release of harmful airborne pollutants within a built environment.

Compact first and second order polarization mode dispersion emulator

NASA Astrophysics Data System (ADS)

Zhang, Yang; Li, Shiguang; Yang, Changxi

2005-08-01

We propose a 1st and 2nd order polarization mode dispersion emulator (PMDE) with one variable differential group delay (DGD) element using birefringence crystals and four polarization controllers (PCs). Monte Carlo simulations demonstrate that the output 1st and 2nd order polarization mode dispersion (PMD) generated by the PMDE consists with statistic theory. Compared with former PMDEs, this design is tunable, lower-cost, and more integrated for fabrication, which shows response time of 150 ?s, response frequency of 3.8 kHz, working wavelength of 1550 nm, total power consumption of less than 3 W, working range of 0---84 ps and 0---3600 ps2 for 1st and 2nd order PMD emulation, respectively. Also, it is programmable and can be controlled by either singlechip or computer. It can be applied to study the outage probability of optical communication systems due to PMD effect and the effectiveness of PMD compensation.

Simulations of the propagation of multiple-FM smoothing by spectral dispersion on OMEGA EP

DOE PAGES

Kelly, J. H.; Shvydky, A.; Marozas, J. A.; ...

2013-02-18

A one-dimensional (1-D) smoothing by spectral dispersion (SSD) system for smoothing focal-spot nonuniformities using multiple modulation frequencies has been commissioned on one long-pulse beamline of OMEGA EP, the first use of such a system in a high-energy laser. Frequency modulation (FM) to amplitude modulation (AM) conversion in the infrared (IR) output, frequency conversion, and final optics affected the accumulation of B-integral in that beamline. Modeling of this FM-to-AM conversion using the code Miró. was used as input to set the beamline performance limits for picket (short) pulses with multi-FM SSD applied. This article first describes that modeling. The 1-D SSDmore » analytical model of Chuang is first extended to the case of multiple modulators and then used to benchmark Miró simulations. Comparison is also made to an alternative analytic model developed by Hocquet et al. With the confidence engendered by this benchmarking, Miró results for multi-FM SSD applied on OMEGA EP are then presented. The relevant output section(s) of the OMEGA EP Laser System are described. The additional B-integral in OMEGA EP IR components upstream of the frequency converters due to AM is modeled. The importance of locating the image of the SSD dispersion grating at the frequency converters is demonstrated. In conclusion, since frequency conversion is not performed in OMEGA EP’s target chamber, the additional AM due to propagation to the target chamber’s vacuum window is modeled.« less

Potential influence of wildfire in modulating climate-induced forest redistribution in a central Rocky Mountain landscape

USGS Publications Warehouse

Campbell, John L.; Shinneman, Douglas

2017-01-01

IntroductionClimate change is expected to impose significant tension on the geographic distribution of tree species. Yet, tree species range shifts may be delayed by their long life spans, capacity to withstand long periods of physiological stress, and dispersal limitations. Wildfire could theoretically break this biological inertia by killing forest canopies and facilitating species redistribution under changing climate. We investigated the capacity of wildfire to modulate climate-induced tree redistribution across a montane landscape in the central Rocky Mountains under three climate scenarios (contemporary and two warmer future climates) and three wildfire scenarios (representing historical, suppressed, and future fire regimes).MethodsDistributions of four common tree species were projected over 90 years by pairing a climate niche model with a forest landscape simulation model that simulates species dispersal, establishment, and mortality under alternative disturbance regimes and climate scenarios.ResultsThree species (Douglas-fir, lodgepole pine, subalpine fir) declined in abundance over time, due to climate-driven contraction in area suitable for establishment, while one species (ponderosa pine) was unable to exploit climate-driven expansion of area suitable for establishment. Increased fire frequency accelerated declines in area occupied by Douglas-fir, lodgepole pine, and subalpine fir, and it maintained local abundance but not range expansion of ponderosa pine.ConclusionsWildfire may play a larger role in eliminating these conifer species along trailing edges of their distributions than facilitating establishment along leading edges, in part due to dispersal limitations and interspecific competition, and future populations may increasingly depend on persistence in locations unfavorable for their establishment.

Stable Stratification Effects on Flow and Pollutant Dispersion in Boundary Layers Entering a Generic Urban Environment

NASA Astrophysics Data System (ADS)

Tomas, J. M.; Pourquie, M. J. B. M.; Jonker, H. J. J.

2016-05-01

Large-eddy simulations (LES) are used to investigate the effect of stable stratification on rural-to-urban roughness transitions. Smooth-wall turbulent boundary layers are subjected to a generic urban roughness consisting of cubes in an in-line arrangement. Two line sources of pollutant are added to investigate the effect on pollutant dispersion. Firstly, the LES method is validated with data from wind-tunnel experiments on fully-developed flow over cubical roughness. Good agreement is found for the vertical profiles of the mean streamwise velocity component and mean Reynolds stress. Subsequently, roughness transition simulations are done for both neutral and stable conditions. Results are compared with fully-developed simulations with conventional double-periodic boundary conditions. In stable conditions, at the end of the domain the streamwise velocity component has not yet reached the fully-developed state even though the surface forces are nearly constant. Moreover, the internal boundary layer is shallower than in the neutral case. Furthermore, an investigation of the turbulence kinetic energy budget shows that the buoyancy destruction term is reduced in the internal boundary layer, above which it is equal to the undisturbed (smooth wall) value. In addition, in stable conditions pollutants emitted above the urban canopy enter the canopy farther downstream due to decreased vertical mixing. Pollutants emitted below the top of the urban canopy are 85 % higher in concentration in stable conditions mostly due to decreased advection. If this is taken into account concentrations remain 17 % greater in stable conditions due to less rapid internal boundary-layer growth. Finally, it is concluded that in the first seven streets the vertical advective pollutant flux is significant, in contrast to the fully-developed case.

Passive advection-dispersion in networks of pipes: Effect of connectivity and relationship to permeability

NASA Astrophysics Data System (ADS)

Bernabé, Y.; Wang, Y.; Qi, T.; Li, M.

2016-02-01

The main purpose of this work is to investigate the relationship between passive advection-dispersion and permeability in porous materials presumed to be statistically homogeneous at scales larger than the pore scale but smaller than the reservoir scale. We simulated fluid flow through pipe network realizations with different pipe radius distributions and different levels of connectivity. The flow simulations used periodic boundary conditions, allowing monitoring of the advective motion of solute particles in a large periodic array of identical network realizations. In order to simulate dispersion, we assumed that the solute particles obeyed Taylor dispersion in individual pipes. When a particle entered a pipe, a residence time consistent with local Taylor dispersion was randomly assigned to it. When exiting the pipe, the particle randomly proceeded into one of the pipes connected to the original one according to probabilities proportional to the outgoing volumetric flow in each pipe. For each simulation we tracked the motion of at least 6000 solute particles. The mean fluid velocity was 10-3 ms-1, and the distance traveled was on the order of 10 m. Macroscopic dispersion was quantified using the method of moments. Despite differences arising from using different types of lattices (simple cubic, body-centered cubic, and face-centered cubic), a number of general observations were made. Longitudinal dispersion was at least 1 order of magnitude greater than transverse dispersion, and both strongly increased with decreasing pore connectivity and/or pore size variability. In conditions of variable hydraulic radius and fixed pore connectivity and pore size variability, the simulated dispersivities increased as power laws of the hydraulic radius and, consequently, of permeability, in agreement with previously published experimental results. Based on these observations, we were able to resolve some of the complexity of the relationship between dispersivity and permeability.

Formation of multiple energy dispersion of H+, He+, and O+ ions in the inner magnetosphere in response to interplanetary shock

NASA Astrophysics Data System (ADS)

Tsuji, H.; Ebihara, Y.; Tanaka, T.

2017-04-01

An interplanetary (IP) shock has a large impact on magnetospheric ions. Satellite observations have shown that soon after arrival of the IP shock, overall intensity of the ions rapidly increases and multiple energy dispersion appears in an energy-time spectrogram of the ions. In order to understand the response of the magnetospheric ions to IP shock, we have performed test particle simulation under the electric and magnetic fields provided by the global magnetohydrodynamic simulation. We reconstructed the differential flux of H+, He+, and O+ ions at (7, 0, 0) Re in GSM coordinates by means of the semi-Lagrangian (phase space mapping) method. Simulation results show that the ions respond to the IP shock in two different ways. First, overall intensity of the flux gradually increases at all pitch angles. As the compressional wave propagates tailward, the magnetic field increases, which accelerates the ions due to the gyrobetatron. Second, multiple energy-time dispersion appears in the reconstructed spectrograms of the ion flux. The energy-time dispersion is caused by the ion moving toward mirror point together with tailward propagating compressional wave at off-equator. The ions are primarily accelerated by the drift betatron under the strong electric field looking dawnward. The dispersion is absent in the spectrogram of equatorially mirroring ions. The dispersion appears at higher energy for heavier ions. These features are consistent with the satellite observations. Because the acceleration depends on bounce phase, the bounce-averaged approximation is probably invalid for the ions during the interval of geomagnetic sudden commencement.