Chapter 14: Seismic imaging using 2- and 3-D migration
Energy Technology Data Exchange (ETDEWEB)
Development of a number of original 3-D migration methods and algorithms is described. The computational efficiency of these algorithms is demonstrated by implementing them on vector and parallel supercomputers. Two-dimensional post-stack and pre-stack depth migration algorithms in the frequency - space domain using implicit finite difference method were also developed and implemented on a parallel computer, and adapted for applications that range from deep crustal imaging of seismic wavefields that involve wide ranging travel times and frequencies. All implicit finite difference migration algorithms were highly prallelized . The one pass 3-D post-stack depth migration algorithm was extensively used for imaging of seismic monitoring data from Cold Lake, Alberta. Experience shows that it provides the right trade-off between accuracy and computational efficiency. A new formulation of the 3-D reverse-time ...
1999-09-01
T-HEMP3D, 3-D Time-Dependent Elastic Plastic Flow
International Nuclear Information System (INIS)
1 - Description of program or function: T-HEMP3D (Transportable HEMP3D) is a derivative of the STEALTH three-dimensional thermodynamics code developed by Science Applications, Inc. Changes were made primarily in three areas: material specification, coordinate generation, and the addition of sliding surface boundary conditions. STEALTH itself is based entirely on the HEMP3D code written at LLNL. HEMP3D solves problems in solid mechanics involving dynamic plasticity and time-dependent material behavior and problems in gas dynamics. 2 - Method of solution: The equations of motion, conservation equations, and constitutive relations for stresses, strains, hydrostatic pressure, the von Mises yield condition, and artificial viscosity are solved by finite difference methods following the format of the HEMP code (NESC Abstract 775). The physical object is divided into zones defined by 8 grid points. The grid (i,j,k) moves with the material, and the mass ...
Diffusion mechanism of implanted Be in GaAs
Energy Technology Data Exchange (ETDEWEB)
The transient enhanced diffusion of low and high dose implanted beryllium in undoped gallium arsenide during post-implant rapid thermal annealing in the temperature range of 700-900 C for 60-240 s has been studied and successfully simulated by the kick-out diffusion model, involving singly positively charged Be interstitials and doubly positively charged Ga self-interstitials. Using the ''plus one'' approach for Ga interstitial generation after implantation with the local Ga interstitial sink concept as well as the appropriate initial and boundary conditions for involved mobile species, and taking into account Fermi-level and built-in electric field effects, the obtained partial differential equations have been solved numerically by means of an explicit finite difference method. The thermal equilibrium concentrations and the diffusivities of Be and Ga interstitials, all as a function of ...
2008-01-15
An imaging comparison of three depth migration algorithms on Foothills datasets
Energy Technology Data Exchange (ETDEWEB)
Three popular imaging methods were used to carry out seismic migration from the Alberta foothills: Kirchhoff integral, reverse-time and explicit finite-difference migration methods. The imaging comparisons of reverse-time, Kirchhoff and f-x post-stack depth migration leads to the following conclusions. All of the three migration methods can provide acceptable and well-imaged intermediate results on the simplest model. Post-stack migration was not very sensitive to the accuracy of the velocity field once a reliable structural stack was available. With respect to the imaging accuracy and calculation cost, the Kirchhoff migration method seemed to be superior to reverse-time and f-x migration algorithms because its results on different kinds of velocity cases were basically kept coherent and consistent with the velocity model given by Mobil. Because of the lack of an accurate velocity field, reverse-time migration was less stable than the other two ...
1999-11-01
An imaging comparison of three depth migration algorithms on Foothills datasets
Energy Technology Data Exchange (ETDEWEB)
Three popular imaging methods were used to carry out seismic migration from the Alberta foothills: Kirchhoff integral, reverse-time and explicit finite-difference migration methods. The imaging comparisons of reverse-time, Kirchhoff and f-x post-stack depth migration leads to the following conclusions. All of the three migration methods can provide acceptable and well-imaged intermediate results on the simplest model. Post-stack migration was not very sensitive to the accuracy of the velocity field once a reliable structural stack was available. With respect to the imaging accuracy and calculation cost, the Kirchhoff migration method seemed to be superior to reverse-time and f-x migration algorithms because its results on different kinds of velocity cases were basically kept coherent and consistent with the velocity model given by Mobil. Because of the lack of an accurate velocity field, reverse-time migration was less stable than the other two ...
1999-01-01
Energy Technology Data Exchange (ETDEWEB)
Finite-difference acoustic-wave modeling and reverse-time depth migration based on the full wave equation are general approaches that can take into account arbitrary variations in velocity and density and can handle turning waves as well. However, conventional finite-difference methods for solving the acoustic- or elastic-wave equation suffer from numerical dispersion when too few samples per wavelength are used. The flux-corrected transport (FCT) algorithm, adapted from hydrodynamics, reduces the numerical dispersion in finite-difference wavefield continuation. The flux-correction procedure endeavors to incorporate diffusion into the wavefield continuation process only where needed to suppress the numerical dispersion. Incorporating the flux-correction procedure in conventional finite-difference modeling or reverse-time migration can provide finite-difference solutions with no ...
1995-11-01
Lagrange-Function Approach to Real-Space Order-N Electronic-Structure Calculations
Energy Technology Data Exchange (ETDEWEB)
The Lagrange functions are a family of analytical, complete, and orthonormal basis sets that are suitable for efficient, accurate, real-space, order-N electronic-structure calculations. Convergence is controlled by a single monotonic parameter, the dimension of the basis set, and computational complexity is lower than that of conventional approaches. In this paper we review their construction and applications in linearscaling electronic-structure calculations.
2006-04-01
Energy Technology Data Exchange (ETDEWEB)
MAGMA is a FORTRAN computer code designed to viscous flow in in situ vitrification melt pools. It models three-dimensional, incompressible, viscous flow and heat transfer. The momentum equation is coupled to the temperature field through the buoyancy force terms arising from the Boussinesq approximation. All fluid properties, except density, are assumed variable. Density is assumed constant except in the buoyancy force terms in the momentum equation. A simple melting model based on the enthalpy method allows the study of the melt front progression and latent heat effects. An indirect addressing scheme used in the numerical solution of the momentum equation voids unnecessary calculations in cells devoid of liquid. Two-dimensional calculations can be performed using either rectangular or cylindrical coordinates, while three-dimensional calculations use rectangular coordinates. All derivatives are approximated by finite differences. The ...
1991-11-01
Analytical method of theoretical simulation of collective hydrodynamic instabilities of intensive flows of discrete radiators, interacting with each other only through the coherent fields of their spontaneous radiation in corresponding media was suggested...
1989-01-01
The physics of tachyons. Part 4. Tachyon electrodynamics
Energy Technology Data Exchange (ETDEWEB)
A new formulation of the theory of tachyons using the same two postulates as in special relativity is applied to the electrodynamics of material media. A discussion of Lagrange`s equations and Hamilton`s equations for `classical` charged tachyons shows that such a formalism is a viable approach. An essay is included on why tachyons can be considered to be localised particles for the purpose of calculations. Tachyonic transformations of the electromagnetic fields D, P, H and M are shown to be the same as for bradyonic transformations. Examples discussed include the electric dipole moment of a tachyonic current loop, constitutive equations, polarisation in tachyonic dielectric materials and the velocity of light in tachyonic dielectric media. This is followed by discussions of the collision energy loss for charged tachyons interacting with a material medium and a mathematical proof that tachyons cannot emit Cherenkov radiation when passing through a bradyonic ...
1998-12-31
In this paper, the authors first describe a fourth order accurate finite difference discretization for both the Laplace equation and the heat equation with Dirichlet boundary conditions on irregular domains. In the case of the heat equation, they use an i...
2004-01-01
Ghost properties of generalized theories of gravitation
Energy Technology Data Exchange (ETDEWEB)
We investigate theories of gravitation, in which spacetime is non-Riemannian and the metric g/sub munu/ is nonsymmetric, for ghosts and tachyons, using a spin-projection operator formalism. Ghosts are removed not by gauge invariance but by a Lagrange multiplier W/sub ..mu../, which occurs due to the breaking of projective invariance in the theory. Unified theories based on a Lagrangian containing a term lambdag/sup munu/g/sub / are proved to contain ghosts or tachyons.
1982-10-15
JPF is an explicit state software model checker for Java bytecode. Today, JPF is a swiss army knife
2005-01-01
Energy Technology Data Exchange (ETDEWEB)
A wavelet transformation is performed over each of the spatial coordinates of the scalar wave equation. This transformed equation is solved directly with a finite-difference scheme for both homogeneous and smooth inhomogeneous media. Wavefield extrapolation is performed completely in the spatial wavelet domain without transforming back into the space domain at each time step. The wavelet coefficients are extrapolated, rather than the wavefield itself. The numerical solution of the scalar wave equation in the spatial wavelet domain is closely related to the finite-difference method because of the compact support of the wavelet bases. Poststack reverse-time migration is implemented as an application. The resolution spaces of the wavelet transform provide a natural framework for multigrid analysis. Migrated images are constructed from various resolution spaces.
1998-03-01
Comparison of fluid-dynamic modeling of flow with velocity-encoded MR imaging
International Nuclear Information System (INIS)
This paper develops a fluid dynamic model using finite difference methods characterizing flow in phantoms simulating in vivo conditions and to compare those results with velocity encoded MR images. The phantom consisted of 1-inch (2.5-cm) tubing with semicircular insert and fluid with viscosity, T1, and T2 comparable to blood. Numeric solutions to Navier-Stokes equations for this system were obtained using finite difference methods, with velocity input function of zero at walls and parabolic at both ends. In resulting color raster (CR) images, color temperature represented velocity value. In velocity-encoded MR images acquired under the same flow conditions, phase is proportional to average velocity during application of flow-encoding gradients. Because these gradients are applied along one direction per acquisition, magnitude and direction of velocity are obtained.
1990-11-25
A nodal integral method for neutron diffusion in hexagonal geometry
International Nuclear Information System (INIS)
A nodal integral method is derived for the monoenergetic, steady-state, fixed source neutron diffusion equation in hexagonal geometry based on a coordinate transformation that maps a parallelogram into a rectangle. The new hexagonal nodal diffusion method is implemented in the computer code HND where the discrete-variable equations are solved via an iterative scheme. Because the new method's equations are derived for a rhombus, they can be solved on a sequence of embedded meshes to study the method's error order. Indeed a preliminary numerical error analysis reveals a second-order error in the mesh size, and comparison with finite difference results obtained with the finite difference based BOLD-VENTURE code indicate the superior accuracy of our new nodal method.
1992-03-08
Comparing S duality to string/(D{minus}5)-brane duality on actions
Energy Technology Data Exchange (ETDEWEB)
We realize the string/(D{minus}5)-brane duality on the action at the classical level between the T{sup 10{minus}D}-compactified heterotic string effective action and the (D{minus}5)-brane effective action in D dimensions by managing a Lagrange multiplier field. In the same way the S duality is also reconstructed as a dual transformation on the T{sup 6}-compactified heterotic string effective action. {copyright} {ital 1997} {ital The American Physical Society}
1997-07-01
International Nuclear Information System (INIS)
The model of transient enhanced diffusion of ion-implanted As is formulated and the finite-difference method for numerical solution of the system of equations obtained is developed. The nonuniform distribution of point defects near the interface and more accurate description of arsenic clustering are simultaneously taken into account. Simulation of As diffusion during rapid annealing gives a reasonable agreement with the experimental data. (authors)
2005-09-01
3-D numerical simulation of the transient electromagnetic pulse coupling to double cylinder object
International Nuclear Information System (INIS)
In this paper, finite-difference time-domain method is used to research the transient electromagnetic pulse (EMP) coupling to the cylinder object which is shielded in another shielding cavity with a hole numerically. The different coupling effectiveness of three transient EMP is simulated. The distribution of the electric field in the shielding cavity and the current on the cavity are researched. The shielding effectiveness is evaluated. (authors)
2006-10-21
Utilization of fractional Brownian motion in constrained least-squares restoration of medical images
International Nuclear Information System (INIS)
Constrained least-squares techniques have been used to produce a well known class of image restoration algorithms. These techniques typically involve minimizing a linear operator on a vector representation of an image, subject to a constraint. For cases where an equality constraint is appropriate the method of Lagrange multipliers can be used to produce a restored image. In this work a fractal textural model, fractional Brownian motion, is used to represent images of interest. Using a variance fractal dimension estimator a non-linear operator, that represents the squared difference between the fractal dimension of the restored image and on a priori value is minimized, subject to the constraint that the norm of the residual between the restored image and available measurement equal the norm of the additive noise.
The physics of tachyons. Part 4. Tachyon electrodynamics
International Nuclear Information System (INIS)
A new formulation of the theory of tachyons using the same two postulates as in special relativity is applied to the electrodynamics of material media. A discussion of Lagrange's equations and Hamilton's equations for 'classical' charged tachyons shows that such a formalism is a viable approach. An essay is included on why tachyons can be considered to be localised particles for the purpose of calculations. Tachyonic transformations of the electromagnetic fields D, P, H and M are shown to be the same as for bradyonic transformations. Examples discussed include the electric dipole moment of a tachyonic current loop, constitutive equations, polarisation in tachyonic dielectric materials and the velocity of light in tachyonic dielectric media. This is followed by discussions of the collision energy loss for charged tachyons interacting with a material medium and a mathematical proof that tachyons cannot emit Cherenkov radiation when passing through a bradyonic ...
International Nuclear Information System (INIS)
A calculation program (URA 6.F4) was elaborated on FORTRAN IV language, that through finite differences solves the unidimensional scalar Helmholtz equation, assuming only one energy group, in spherical cylindrical or plane geometry. The purpose is the determination of the flow distribution in a reactor of spherical cylindrical or plane geometry and the critical dimensions. Feeding as entrance datas to the program the geometry, diffusion coefficients and macroscopic transversals cross sections of absorption and fission for each region. The differential diffusion equation is converted with its boundary conditions, to one system of homogeneous algebraic linear equations using the box integration technique. The investigation on criticality is converted then in a succession of eigenvalue problems for the critical eigenvalue. In general, only is necessary to solve the first eigenvalue and its corresponding eigenvector, employing the power method. The ...
1993-11-18
Energy Technology Data Exchange (ETDEWEB)
The standards mapping methods are limited by the lack of information between surface measures and wells profiles. The cross well seismic survey has proven being an effective high-resolution method for reservoirs characterization due its high frequency band. In this work, a 2-D finite differences algorithm for numerical cross well seismic simulation was developed, based in the approach of 2nd order for the time derivatives and 4th order for the space derivatives. In relation to the imaging techniques, it was introduced a reverse time migration method based in the solutions of the wave equation in time for cross well seismic data by the method of the finite differences exclusively. The results of reverse time migration has showed that layers with 90 cm thickness for P-waves could be visualized in migrated sections with excellent resolution. The high frequency cross well seismic data combined with reverse ...
2003-07-01
Energy Technology Data Exchange (ETDEWEB)
The reversion time migration, as one of reflection seismic prospecting migration, reverse propagates a wave equation by reverse promoting the differential time, and conducts the migration. The method has such merits as no limit to the inclined angle of a reflection, correspondence to a complex topography etc. And, there is a finite difference time domain (FDTD) that uses the Maxwell equation as it is and makes a difference approximation in the methods differentiating and resolving electro-magnetic fields. It is also possible to apply the reverse time migration onto it. In this speech, using an explosive reflection face, a slope 2 layer structure and models of its reflection points are synthesized by a ground penetrating radar using the FDTD, and reverse time migration is conducted. As a result, it is clarified that the reverse time migration is valid for the ground penetrating radar recording. And, the cases including a diffusion item in the ...
1999-02-01
B-spline methods for radial Dirac equations
International Nuclear Information System (INIS)
Although B-spline techniques have been used to solve two-point boundary value problems with Dirac Hamiltonians for more than 20 years, the treatment of boundary conditions is still a matter of controversy. Spurious, non-physical, solutions are endemic when boundary conditions are not handled correctly. These pathological problems are absent when traditional finite difference methods are used as in computer packages such as GRASP. Accurate approximation using both finite differences and B-splines depends on controlling local approximation errors, and this common property suggests no a priori reason to suppose that B-spline algorithms should be more prone to generate spurious solutions. The relativistic Bloch operators of [24], when added to the Dirac differential operator, permit the construction of a self-adjoint differential operator for the two-point boundary value problem on a finite interval. ...
2009-03-14
Reverse time migration in tilted transversely isotropic media
Energy Technology Data Exchange (ETDEWEB)
This paper presents a reverse time migration (RTM) method for the migration of shot records in tilted transversely isotropic (TTI) media. It is based on the tilted TI acoustic wave equation that was derived from the dispersion relation. The RTM is a full depth migration allowing for velocity to vary laterally as well as vertically and has no dip limitations. The wave equation is solved by a tenth-order finite difference scheme. Using 2D numerical models, we demonstrate that ignoring the tilt angle will introduce both lateral and vertical shifts in imaging. The shifts can be larger than 0.5 wavelength in the vertical direction and 1.5 wavelength in the lateral direction.
2004-07-01
Reflection of electromagnetic waves by a nonuniform plasma layer covering a metal surface
International Nuclear Information System (INIS)
Reflection coefficients of electromagnetic waves in a nonuniform plasma layer with electrons, positive ions and negative ions, covering a metal surface are investigated by using the finite-difference-time-domain method. It is shown that the reflection coefficients are influenced greatly by the density gradient on the layer edge, layer thickness and electron proportion, i.e., the effect of the negative ions. It is also found that low reflection or high attenuation can be reached by properly choosing high electron proportion, thick plasma layer, and smooth density gradient in the low frequency regime, but sharp density gradient in the high frequency regime. (authors)
2008-07-01
Energy Technology Data Exchange (ETDEWEB)
A pre-stack migration algorithm for elastic waves in two-dimensional variable-velocity media is developed, implemented, and tested. The algorithm operates in the time-space domain and is based on reverse-time finite-difference extrapolation of elastic waves. The algorithm is explained and demonstrated in the context of imaging of elastic vertical seismic profile data, but is applicable to any source-recorder geometry. Synthetic test examples include a point diffractor, laterally homogeneous layers, and the flank of a salt dome.
1986-03-01
A finite three-dimensional layer-by-layer photonic crystal with planar defects in a layer is shown to drastically modify the spontaneous emission rate of an embedded dipole. Finite-difference time-domain calculations with one quarter symmetric boundary condition and perfectly matched layer demonstrate the strong enhancement effects induced by the cavity resonance of defect modes and band-edge resonant modes. Simulation shows that the emission spectra are quite different when the position or polarization of the dipole is changed. Moreover, the extraction efficiency is calculated to observe the percentage of light leakage through a substrate.
2010-01-01
Hydrothermal coupling in a rough fracture
Heat exchange during laminar flow is studied at the fracture scale on the basis of the Stokes equation. We used a synthetic aperture model (a self-affine model) that has been shown to be a realistic geometrical description of the fracture morphology. We developed a numerical modelling using a finite difference scheme of the hydrodynamic flow and its coupling with an advection/conduction description of the fluid heat. As a first step, temperature within the surrounding rock is supposed to be constant. Influence of the fracture roughness on the heat flux through the wall, is estimated and a thermalization length is shown to emerge. Implications for the Soultz-sous-For\\^{e}ts geothermal project are discussed.
2006-01-01
Dynamics of a reverse osmosis unit with application to pulsating regimes for process optimization
British Library Electronic Table of Contents (United Kingdom)
We propose unsteady-state reverse osmosis cell modelling in two dimensions. The convection-diffusion equation describing the concentration of the relevant chemical species is solved by a finite difference technique, while the velocity field is described by empirical expressions for spiral-wound membrane cells. A non-constant permeability is introduced to take into account the effects of membrane compaction at high operating pressures. The role of concentration polarization is discussed for different values of the parameters describing the global process. Finally, the model is applied to predict the effects of a pulsating flow where a cyclic pressure feed is adopted to enhance the permeate flux. In this context, an experimental validation of the model is proposed.
2011-01-01
Analysis of laminar flow heat transfer in uniform temperature circular tubes with tape inserts
Constant property, laminar flow heat transfer in a semicircular tube with uniform wall temperature has been analyzed to define the lower bound of heat transfer augmentation in circular tubes with twisted-tape inserts. Two thermal boundary conditions, which correspond to the two extremes of the fin effect of twisted tapes encountered in practical applications, are considered. Numerical solutions, employing finite-difference formulations for the governing momentum and energy equations were carried out for the thermal entrance region and for fully developed flow.
1986-05-01
Supersymmetric homogeneous quantum cosmologies coupled to a scalar field
Energy Technology Data Exchange (ETDEWEB)
Recent work on [ital N]=2 supersymmetric Bianchi type IX cosmologies coupled to a scalar field is extended to a general treatment of homogeneous quantum cosmologies with explicitly solvable momentum constraints, i.e., Bianchi types I, II, VII, VIII in addition to the Bianchi type IX, and special cases, namely, the Freidmann universes, the Kantowski-Sachs space, and Taub-NUT space. In addition to the earlier explicit solution of the Wheeler-DeWitt equation for Bianchi type IX, describing a virtual wormhole fluctuation, an additional explicit solution is given and identified with the no-boundary state.''
1994-01-15
Supersymmetric homogeneous quantum cosmologies coupled to a scalar field
International Nuclear Information System (INIS)
Recent work on N=2 supersymmetric Bianchi type IX cosmologies coupled to a scalar field is extended to a general treatment of homogeneous quantum cosmologies with explicitly solvable momentum constraints, i.e., Bianchi types I, II, VII, VIII in addition to the Bianchi type IX, and special cases, namely, the Freidmann universes, the Kantowski-Sachs space, and Taub-NUT space. In addition to the earlier explicit solution of the Wheeler-DeWitt equation for Bianchi type IX, describing a virtual wormhole fluctuation, an additional explicit solution is given and identified with the ''no-boundary state.''.
International Nuclear Information System (INIS)
We develop the notions of fusion for representations of the WA_2 algebra along the lines of Feigin and Fuchs. We present some explicit calculations for a WA_2 minimal model. (orig.).
1995-01-01
The explicit Laplace transform for the Wishart process
We derive the explicit formula for the joint Laplace transform of the Wishart process and its time integral which extends the original approach of Bru. We compare our methodology with the alternative results given by the variation of constants method, the linearization of the Matrix Riccati ODE's and the Runge-Kutta algorithm. The new formula turns out to be fast, accurate and very useful for applications when dealing with stochastic volatility and stochastic correlation modelling.
2011-01-01
Energy Technology Data Exchange (ETDEWEB)
In this note we prove to all orders in the small scale expansion that all off-shell parameters which appear in the chiral effective Lagrangian with explicit DELTA(1232) isobar degrees of freedom can be absorbed into redefinitions of certain low-energy constants and are therefore redundant.
2010-01-18
A general theorem on the GBDT version of the B\\"acklund-Darboux transformation for systems rationally depending on the spectral parameter is treated and its applications to nonlinear equations are given. Explicit solutions of direct and inverse problems for Dirac-type systems, including systems with singularities, and for the system auxiliary to the $N$-wave equation are reviewed. New results on explicit construction of the wave functions for radial Dirac equation are obtained.
2009-01-01
Prejudice against overweight people is rife. However, there is a paucity of research on the underlying reasons for it. In two studies the relationship between body image, the tendency to make physical appearance-related comparisons (PACS), and both explicit and implicit anti-fat attitudes was examined. In Study 1 (n = 227) people with a high tendency to make physical appearance-related comparisons (high PACS scorers) reported lower self-appearance evaluation, but higher appearance orientation and explicit anti-fat attitudes. The PACS fully mediated the relationship between appearance orientation and explicit anti-fat attitudes. Study 2 (n = 134) found that the PACS also mediated the relationship between appearance orientation and implicit anti-fat attitudes. Thus, individual differences in factors such as body image and the tendency to make appearance-related comparisons, appear to play a central role in both ...
2007-07-24
Multi-scale modeling of fiber and fabric reinforced cement based composites
With an increased use of fiber reinforced concrete in structural applications, proper characterization techniques and development of design guides are needed. This dissertation presents a multi-scale modeling approach for fiber and fabric reinforced cement-based composites. A micromechanics-based model of the yarn pullout mechanism due to the failure of the interfacial zone is presented. The effect of mechanical anchorage of transverse yarns is simulated using nonlinear spring elements. The yarn pullout mechanism was used in a meso-scale modeling approach to simulate the yarn bridging force in the crack evolution process. The tensile stress-strain response of a tension specimen that experiences distributed cracking can be simulated using a generalized finite difference approach. The stiffness degradation, tension stiffening, crack spacing evolution, and crack width characteristics of cement composites can be derived using matrix, interface and ...
2008-01-01
Well Test Analysis in Naturally Fissured, Geothermal Reservoirs with Fracture Skin
Mineral deposition or alteration is commonly found at fracture-block interfaces is fissured, geothermal reservoirs. In response to pressure reduction in the fissures such mineralization, if less permeable than the matrix rock, will retard the flow of fluid from the blocks to the fissures and is termed fracture skin in this paper. The problem of fluid flow to a production well in a double-porosity reservoir with fracture skin was analyzed theoretically. One of the findings of the analysis was that fully transient block-to-fissure flow can be approximated by pseudo-steady state flow if fracture skin permeability is sufficiently low. Type curves generated by numerical inversion of Laplace transform solutions are used to cooroborate the results of a finite-difference model of steam transport to a well in a naturally fissured, geothermal reservoir with fracture skin.
1983-12-15
The new computer program for three dimensional relativistic hydrodynamical model
An effective computer program for three dimensional relativistic hydrodynamical model has been developed. It implements a new approach to the early hot phase of relativistic heavy-ion collisions. The computer program simulates time-space evolution of nuclear matter in terms of ideal-fluid dynamics. Equations of motions of hydrodynamics are solved making use of finite difference methods. Commonly-used algorithms of numerical relativistic hydrodynamics RHLLE and MUSTA-FORCE have been applied in simulations. To speed-up calculations, parallel processing has been made available for solving hydrodynamical equations. The test results of simulations for 3D, 2D and Bjorken expansion are reported in this paper. As a next step we plan to implement the hadronization algorithm by implementing the continuous particle emission for freeze-out and comparing it with Cooper-Frye formula.
2006-01-01
Resonant frequencies and Q factors of dielectric parallelepipeds by measurement and by FDTD
Energy Technology Data Exchange (ETDEWEB)
This paper describes the measurement and computation of the resonant frequencies and the associated Q factors of dielectric parallelepipeds made of high-permittivity, low-loss ceramic materials. Each resonance peak is measured separately with a fine frequency step. A curve-fitting method is used to accurately estimate the resonant frequency and 3 dB bandwidth from the somewhat noisy measured data. The finite-difference time-domain method is used to compute the initial portion of the backscattered field due to a Gaussian pulse plane wave. The time response is then extended to zero value by Prony`s method. The measured and computed data is compared for a parallelepiped resonator of permittivity 37.84.
1994-12-31
Energy Technology Data Exchange (ETDEWEB)
The paper presents a procedure for the calculation of an absorption-refrigeration plant that uses an LiBr/ZnBr{sub 2}/CH{sub 3}OH solution and is driven throughout the year by heat received from a solar pond. The calculations for the solar pond are based on a finite-difference solution of the transient heat conduction differential equation during the year, and the operation of the refrigeration unit is simulated by the use of functions describing the thermodynamic behaviour of the working medium. On the basis of the procedure devised, correlations and nomographs are developed that describe the seasonally optimum operation of the proposed composite system. The nomographs may be employed for a first estimation in the design of such systems. (Author).
1992-07-01
British Library Electronic Table of Contents (United Kingdom)
The effects of thermal radiation and viscous dissipation on magneto-hydrodynamic (MHD) unsteady free-convection flow over a semi-infinite vertical porous plate are analysed. The fluid considered is non-gray (absorption coefficient dependent on wave length). The Network Simulation Method is used to solve the boundary-layer equations based on the finite-difference formulation; only discretization of the spatial co-ordinates is necessary, while time remains as a real continuous variable. This method provides a solution for both transient and steady-state problems at the same time, and programming does not require manipulation of the sophisticated mathematical software that is inherent in other numerical methods. The velocity, temperature, local skin-friction and local Nusselt number are studi...
2007-01-01
Energy Technology Data Exchange (ETDEWEB)
The forward modelling and the prestack reverse time migration of seismic P-SV wave field was carried out in 2-D models of isotropic and anisotropic media which allow separation of P-SV and SH motion. The P-SV wave field can be described by a system of hyperbolic, first order differential equations in terms of particle velocity and stress. The system of five equations and five unknowns, namely horizontal (U) and vertical (V) velocity components, and three components of stress (T{sub xx}, T-z{sub z} and T{sub xz}) was solved numerically using second order space and forth order time finite differences operators. In order to attenuate numerical dispersion, a staggered grid was used. (author). 48 refs., 5 figs
1993-12-31
Mathematical modeling of the behavior of geothermal systems under exploitation
Energy Technology Data Exchange (ETDEWEB)
Analytical and numerical methods have been used in this investigation to model the behavior of geothermal systems under exploitation. The work is divided into three parts: (1) development of a numerical code, (2) theoretical studies of geothermal systems, and (3) field applications. A new single-phase three-dimensional simulator, capable of solving heat and mass flow problems in a saturated, heterogeneous porous or fractured medium has been developed. The simulator uses the integrated finite difference method for formulating the governing equations and an efficient sparse solver for the solution of the linearized equations. In the theoretical studies, various reservoir engineering problems have been examined. These include (a) well-test analysis, (b) exploitation strategies, (c) injection into fractured rocks, and (d) fault-charged geothermal reservoirs.
1982-01-01
Hydromagnetic free convection currents effects on boundary layer thickness
Energy Technology Data Exchange (ETDEWEB)
In this study we discuss an unsteady free convection MHD flow past semi-infinite vertical porous plate. We have considered the flow in the presence of a strong magnetic field and therefore the electromagnetic force is very large. This brings in the phenomenon of Hall and Ion-slip currents. The effects of these two parameters together with that of viscous dissipation and radiation absorption among others on velocity, temperature and concentration profiles are presented. The profiles are presented graphically. As the partial differential equations governing this problem are highly non-linear they are solved numerically by a finite difference method. It is found that in presence of heating of the plate by free convection current the velocity boundary layer thickness decreases.
2010-06-15
Hydromagnetic free convection currents effects on boundary layer thickness
British Library Electronic Table of Contents (United Kingdom)
In this study we discuss an unsteady free convection MHD flow past semi-infinite vertical porous plate. We have considered the flow in the presence of a strong magnetic field and therefore the electromagnetic force is very large. This brings in the phenomenon of Hall and Ion-slip currents. The effects of these two parameters together with that of viscous dissipation and radiation absorption among others on velocity, temperature and concentration profiles are presented. The profiles are presented graphically. As the partial differential equations governing this problem are highly non-linear they are solved numerically by a finite difference method. It is found that in presence of heating of the plate by free convection current the velocity boundary layer thickness decreases.
2010-01-01
Energy Technology Data Exchange (ETDEWEB)
The authors study, in forced convection, the evaporation of an Ostwaldian film flowing over an isothermal inclined plane surface to determine the influence of the behaviour index of the liquid on the dynamic and thermal characteristics of liquid-air system. The liquid flow is considered partially two-dimensional whereas for the air it is two-dimensional. The coupled equations with the interfacial conditions are solved using a fully implicit finite differences method. From the study, it appears that the behaviour index influences considerably the transfers which are more important for pseudo-plastic liquids than for dilatant ones. (authors)
2003-12-01
Estimation of plasma density by surface plasmons for surface-wave plasmas
International Nuclear Information System (INIS)
An estimation method of plasma density based on surface plasmons theory for surface-wave plasmas is proposed. The number of standing-wave is obtained directly from the discharge image, and the propagation constant is calculated with the trim size of the apparatus in this methods, then plasma density can be determined with the value of 9.1 x 1017 m-3. Plasma density is measured using a Langmuir probe, the value is 8.1 x 1017 m-3 which is very close to the predicted value of surface plasmons theory. Numerical simulation is used to check the number of standing-wave by the finite-difference time-domain (FDTD) method also. All results are compatible both of theoretical analysis and experimental measurement. (authors)
2008-12-01
Combined radiation and convection in absorbing, emitting, nongray gas-particulate tube flow
Energy Technology Data Exchange (ETDEWEB)
The interaction of thermal radiation with conduction and convection in thermally developing absorbing, emitting, nongray gas-particulate turbulent suspension flow through a circular tube is investigated. The contribution of thermal radiation is obtained through evaluation of the total hemispherical emittance of the particulate cloud and through evaluation of single band absorptances for molecular gases, modified to account for the interaction with the particles. The governing differential equation is derived as a (nonlinear) energy equation, coupled with integral equations to find the thermal radiation contributions. The energy equation is solved numerically by an implicit finite difference method with an iterative procedure. Qualitative results for Nusselt numbers are shown for a variety and range of parameters, such as optical thickness of particulates and single molecular gas bands, relative gas band position and band width, and temperature ...
1987-05-01
Combined radiation and convection in absorbing, emitting, non-Gray gas-particulate tube flow
Energy Technology Data Exchange (ETDEWEB)
The interaction of thermal radiation with conduction and convection in thermally developing absorbing, emitting, non-gray gas particulate turbulent suspension flow through a circular tube is investigated. The contribution of thermal radiation is obtained through evaluation of the total hemispherical emittance of the particulate cloud and through evaluation of single band absorptances for molecular gases, modified to account for the interaction with the particles. The governing differential equation is derived as a (nonlinear) energy equation, coupled with integral equations to find the thermal radiation contributions. The energy equation is solved numerically by an implicit finite difference with its iterative procedure. Qualitative results for Nusselt numbers are shown for a variety and range of parameters, such as optical thickness of particulates and single molecular gas bands, relative gas band position and band width, and temperature ...
1985-01-01
British Library Electronic Table of Contents (United Kingdom)
The heating of water layer using microwave oven with a rectangular waveguide has been studied both numerically and experimentally. The mathematical model is validated with the experimental data. The transient Maxwell's equations are solved by using the Finite Difference Time Domain (FDTD) method to describe the electromagnetic field inside the waveguide and sample. The temperature profile and velocity field within sample are determined by the solutions of the momentum, energy and Maxwell's equations. In this study, the effects of physical parameters, e.g. microwave power level, placement of sample inside the waveguide, volume of sample, are studied. The distribution of electric field, temperature profile and velocity field are presented in details. The results show good agreement between s...
2012-01-01
3-D prestack migration in anisotropic media
Energy Technology Data Exchange (ETDEWEB)
A three-dimensional (3-D) prestack reverse-time migration algorithm for common-source P-wave data from anisotropic media is developed and illustrated by application to synthetic data. Both extrapolation of the data and computation of the excitation-time imaging condition are implemented using a second-order finite-difference solution of the 3-D anisotropic scalar-wave equation. Poorly focused, distorted images are obtained if data from anisotropic media are migrated using isotropic extrapolation; well focused, clear images are obtained using anisotropic extrapolation. A priori estimation of the 3-D anisotropic velocity distribution is required. Zones of anomalous, directionally dependent reflectivity associated with anisotropic fracture zones are detectable in both the 3-D common-source data and the corresponding migrated images.
1993-01-01
3-D depth migration via McClellan transformations
Energy Technology Data Exchange (ETDEWEB)
Three-dimensional seismic wavefields may be extrapolated in depth, one frequency at a time, by two-dimensional convolution with a circularly symmetric, frequency- and velocity-dependent filter. This depth extrapolation, performed for each frequency independently, lies at the heart of 3-D finite-difference depth migration. The computational efficiency of 3-D depth migration depends directly on the efficiency of this depth extrapolation. McClellan transformations provide an efficient method for both designing and implementing two-dimensional digital filters that have a particular form of symmetry, such as the circularly symmetric depth extrapolation filters used in 3-D depth migration. Given the coefficients of one-dimensional, frequency- and velocity-dependent filters used to accomplish 2-D depth migration, McClellan transformations lead to a simple and efficient algorithm for 3-D depth migration. 21 refs., 12 figs.
1990-01-01
On cosmological type solutions in multidimensional model with Gauss-Bonnet term
A (n +1)-dimensional Einstein-Gauss-Bonnet (EGB) model is considered. For diagonal cosmological type metrics the equations of motion are reduced to a set of Lagrange equations. The effective Lagrangian contains two minisuperspace metrics on R^n. The first one is the well-known 2-metric of pseudo-Euclidean signature and the second one is the Finslerian 4-metric that is proportional to n-dimensional Berwald-Moor 4-metric. When a synchronous-like time gauge is considered the equations of motion are reduced to an autonomous system of first order differential equations. For the case of the pure Gauss-Bonnet model two exact solutions with power-law and exponential dependence of scale factors (w.r.t. synchronous-like variable) are obtained. In cosmological case the first solution was obtained earlier by A. Toporensky and P. Tretyakov for n = 5 (and verified for n = 6,7). A generalization of the effective Lagrangian to the Lowelock case is conjectured. This hypothesis ...
2009-01-01
Non-Riemannian geometrical asymmetrical damping stresses on the Lagrange instability of shear flows
It is shown that the physical interpretation of Elie Cartan three-dimensional space torsion as couple asymmetric stress, has the effect of damping, previously Riemannian unstable Couette planar shear flow, leading to stability of the flow in the Lagrangean sense. Actually, since the flow speed is inversely proportional to torsion, it has the effect of causing a damping in the planar flow atenuating the instability effect. In this sense we may say that Cartan torsion induces shear viscous asymmetric stresses in the fluid, which are able to damp the instability of the flow. The stability of the flow is computed from the sectional curvature in non-Riemannian three-dimensional manifold. Marginal stability is asssumed by making the sectional non-Riemannian curvature zero, which allows us to determine the speeds of flows able to induce this stability. The ideas discussed here show that torsion plays the geometrical role of magnetic field in hydromagnetic instability of Couette flows recently ...
2007-01-01
Geophysical remote sensing of water reservoirs suitable for desalinization.
Energy Technology Data Exchange (ETDEWEB)
In many parts of the United States, as well as other regions of the world, competing demands for fresh water or water suitable for desalination are outstripping sustainable supplies. In these areas, new water supplies are necessary to sustain economic development and agricultural uses, as well as support expanding populations, particularly in the Southwestern United States. Increasing the supply of water will more than likely come through desalinization of water reservoirs that are not suitable for present use. Surface-deployed seismic and electromagnetic (EM) methods have the potential for addressing these critical issues within large volumes of an aquifer at a lower cost than drilling and sampling. However, for detailed analysis of the water quality, some sampling utilizing boreholes would be required with geophysical methods being employed to extrapolate these sampled results to non-sampled regions of the aquifer. The research in this report addresses using seismic and EM methods in ...
2009-12-01
Energy Technology Data Exchange (ETDEWEB)
In this work the development of a third order scheme of finite differences centered in mesh is presented and it is applied in the numerical solution of those diffusion equations in multi groups in stationary state and X Y geometry. Originally this scheme was developed by Hennart and del Valle for the monoenergetic diffusion equation with a well-known source and they show that the one scheme is of third order when comparing the numerical solution with the analytical solution of a model problem using several mesh refinements and boundary conditions. The scheme by them developed it also introduces the application of numeric quadratures to evaluate the rigidity matrices and of mass that its appear when making use of the finite elements method of Galerkin. One of the used quadratures is the open quadrature of 4 points, no-standard, of Newton-Cotes to evaluate in approximate form the elements of the rigidity matrices. The other quadrature is that of ...
2003-07-01
A 3-D hydrodynamic dispersion model for modeling tracer transport in Geothermal Reservoirs
International Nuclear Information System (INIS)
A 3-D hydrodynamic dispersion model for tracer transport is developed and implemented into the TOUGH2 EOS3 (T2R3D) module. The model formulation incorporates a full dispersion tensor, based on a 3-D velocity field with a 3-D, irregular grid in a heterogeneous geological system. Two different weighting schemes are proposed for spatial average of 3-D velocity fields and concentration gradients to evaluate the mass flux by dispersion and diffusion of a tracer or a radionuclide. This new module of the TOUGH2 code is designed to simulate processes of tracer/radionuclide transport using an irregular, 3-D integral finite difference grid in non-isothermal, three-dimensional, multiphase, porous/fractured subsurface systems. The numerical method for this transport module is based on the integral finite difference scheme, as in the TOUGH2 code. The major assumptions of the tracer transport module are: (a) a tracer ...
Explicit Randomness is not Necessary when Modeling Probabilistic Encryption
Although good encryption functions are probabilistic, most symbolic models do not capture this aspect explicitly. A typical solution, recently used to prove the soundness of such models with respect to computational ones, is to explicitly represent the dependency of ciphertexts on random coins as labels. In order to make these label-based models useful, it seems natural to try to extend the underlying decision procedures and the implementation of existing tools. In this paper we put forth a more practical alternative based on the following soundness theorem. We prove that for a large class of security properties (that includes rather standard formulations for secrecy and authenticity properties), security of protocols in the simpler model implies security in the label-based model. Combined with the soundness result of (\\textbf{?}) our theorem enables the translation of security results in unlabeled symbolic models to computational security.
2006-01-01
Distinguishing three levels in explicit self-awareness
British Library Electronic Table of Contents (United Kingdom)
This paper focuses on the development of explicit self-awareness in children. Mirror self-recognition has been the most popular paradigm used to assess this ability in children. Nevertheless, according to Rochat (2003), there are, at least, three different levels of explicit self-awareness. We therefore designed three different self-recognition tasks, each corresponding to one of these levels (a mirror self-recognition task, a picture self-recognition task and a masked self-recognition task). We observed a decrease in performance across the three tasks. This supports a developmental scale in self-awareness. Besides, the masked self-recognition performance makes it possible to assess the final and the most sophisticated level of self-awareness, i.e. the external self. To our best knowledge,...
2011-01-01
Unified kinematics of bradyons and tachyons in six-dimensional space-time
International Nuclear Information System (INIS)
An explicit form of the transformation matrix in sin-dimensional space-time is given both for subluminal and superluminal cases, and it is shown how a single formalism can be used for both bradyons and tachyons. (author).
Toward a Unified Timestamp with explicit precision
UK PubMed Central (United Kingdom)
Demographic and health surveillance (DS) systems monitor and document individual and group-level processes in well-defined populations over long periods of time. The resulting data are complex...Full Text Available
2005-03-30
Tachyons and the second law of black hole physics
International Nuclear Information System (INIS)
It is shown that the usual proof of the second law of black hole physics breaks down if there are tachyons present in the vicinity of a black hole. Explicit cases are discussed where a tachyon of positive energy falling into the Kerr singularity actually decreases the area of the Kerr black hole. (author).
Some perturbative results for two-dimensional gravity
Energy Technology Data Exchange (ETDEWEB)
Perturbative approach to two-dimensional gravity and supergravity is considered. One-loop renormalization of the central charge of SL(2,R) Kac-Moody algebra is calculated perturbatively by functional integration and by explicit calculations of the Feynman diagrams. Also the wavefunction renormalization and the anomalous dimensions in the presence of gravity are calculated.
1990-04-20
UK PubMed Central (United Kingdom)
Four-dimensional (4D) radiotherapy is the explicit inclusion of the temporal changes in anatomy during the imaging, planning, and delivery of radiotherapy. One key component of 4D radiotherapy planning...Full Text Available
2008-04-01
UK PubMed Central (United Kingdom)
An explicit proposal for experiments leading to abrupt transitions in spatially extended bacterial populations in a Petri dish is presented on the basis of an exact formula obtained through an analytic...Full Text Available
2008-12-02
Hamiltonian formulation of the N=1 D=10 supergravity coupled to super Yang-Mills theory
Energy Technology Data Exchange (ETDEWEB)
The authors derive in full details the Hamiltonian formulation of the N = 1D = 10 supergravity coupled to super Yang-Mills theory. In particular, they present the explicit form of the first class constraints and compute the constraints gauge algebra.
1988-01-01
Hamiltonian BRST formalism for gauge-invariant systems with closed algebra
International Nuclear Information System (INIS)
The systems for which the algebra of gauge transformations in the lagrangian formalism is closed, are considered. The hamiltonian BRST charge and the BRST-invariant hamiltonian are found explicitly. Their expansions in powers of the ghost variables contain, in general, an infinite number of terms. (orig.).
1991-11-01
Fragment-Based Learning of Visual Object Categories in Non-Human Primates
UK PubMed Central (United Kingdom)
When we perceive a visual object, we implicitly or explicitly associate it with an object category we know. Recent research has shown that the visual system can use local, informative image fragments...Full Text Available
Explicit risk in acute coronary syndrome management
UK PubMed Central (United Kingdom)
At least implicitly, most clinical decisions represent an integration of disease and treatment-based risk assessments. Often, as is the case with acute coronary syndrome (ACS), these decisions need...Full Text Available
2009-06-01
Exactly solvable field-theoretical model with tachyons
Energy Technology Data Exchange (ETDEWEB)
Explicit soliton solutions describing the inelastic interaction between sub- and superluminal particles are found within the framework of a new integrable model of relativistic classical field theory. The corresponding energies are nonnegative irrespective of the choice of reference frame.
1988-09-01
UK PubMed Central (United Kingdom)
Diabetic pregnancies are characterized by chronic metabolic insults, including iron deficiency, that place the developing brain at risk and for memory impairment later in life. A behavioral...Full Text Available
2009-11-01
Accurate, stable, explicit, parabolized navier-stokes solver for high speed flows
Energy Technology Data Exchange (ETDEWEB)
A stable, accurate, and efficient implementation of MacCormack's explicit algorithm for the Parabolized Navier-Stokes equations is demonstrated. The familiar problem of decoding the conservative axial flux vector is solved, resulting in accurate, smooth dependent variable profiles through the viscous-layer sonic line. Source terms due to transformation of the parabolized governing equations into the computational plane and the equations into the computational plane and the resulting metric differencing have been identified and eliminated through inclusion of appropriate geometric conservation law terms. Test cases computed include two- and three-dimensional supersonic and hypersonic flow at laminar and turbulent Reynolds numbers. The computed results demonstrate very good agreement with experiment and with solutions of the full Navier-Stokes equations. Computational times required for the MacCormack explicit PNS code are approximately ...
1986-01-01
A comparison of the level of care predictions of six long-term care patient assessment systems.
UK PubMed Central (United Kingdom)
Six patient assessment systems that have explicit decision rules for replicating team judgments on level of care patient placement were selected for analysis. The six were selected because of their...Full Text Available
1980-11-01
Symmetry analysis and exact explicit solutions for Kadomtsev-Petviashvili-Burgers equation
British Library Electronic Table of Contents (United Kingdom)
We apply the group theory to Kadomtsev-Petviashvili-Burgers (KPBII) equation which is a natural model for the propagation of the two-dimensional damped waves. In correspondence with the generators of the symmetry group allowed by the equation, new types of symmetry reductions are performed. Some new exact solutions are obtained, which can be in the form of solitary waves and periodic waves. Specially, our solutions indicate that the equation may have time-dependent nonlinear shears. Such exact explicit solutions and symmetry reductions are important in both applications and the theory of nonlinear science.
2011-01-01
Matrix realization of string algebra axioms and conditions of invariance
International Nuclear Information System (INIS)
The matrix representations of Witten's and B-algebras of the field string theory in finite dimensional space of the ghost states are suggested for the case of Virasoro algebra truncated to its SU(1,1) subalgebra. In this case all algebraic operations of Witten's and B-algebras are realized in explicit form as some matrix operations in the graded complex vector space. The structure of string action coincides with the universal non-linear cubic matrix form of action for the gauge field theories. These representations lead to matrix conditions of theory invariance which can be used for finding of the explicit form of corresponding operators of the string algebras. (author).
Many-Body Rate Limit on Photoassociation of a Bose-Einstein Condensate
We briefly report on zero-temperature photoassociation of a Bose-Einstein condensate, focusing on the many-body rate limit for atom-molecule conversion. An upgraded model that explicitly includes spontaneous radiative decay leads to an unanticipated shift in the position of the photoassociation resonance, which affects whether the rate (constant) maximizes or saturates, as well as the limiting value itself. A simple analytical model agrees with numerical experiments, but only for high density. Finally, an explicit comparison with the two-body unitary limit, set by the size of the condensate, finds that the many-body rate limit is generally more strict.
2010-01-01
The thermal dissolver, the main reactor of the SRC unit, has suffered a recurring problem. Specifically, it has been observed that whenever the reactor vessel is cooled to below 400/sup 0/F, its bottom head gasket leaks. An analysis of the thermal stress induced in the gasket, owing to transients across the bottom head flange, was sought. The analysis was facilitated by judiciously dividing a symmetric section of the reactor into 79 differential elements. Heat balances have been developed around each element. A numerical technique, the backward finite-difference approach, was employed to obtain the thermal behavior across the bottom head flange as a function of reactor heat-up time. The analysis performed affords an explanation for the failure of the gasket. Based on results of this work, recommendations have been suggested to provide the gasket and bolt stress requirements that are necessary to avoid leaks due to temperature changes.
1983-08-01
Validation of WIMS-AECL reactivity device calculations for CANDU reactor
Energy Technology Data Exchange (ETDEWEB)
An important component of the overall program to validate WIMS-AECL for use with RFSP in the analysis of CANDU-6 reactors for design and safety analysis calculations is the validation of calculations of incremental cross sections used to represent reactivity devices. A method has been developed for the calculation of the three-dimensional neutron flux distribution in and around CANDU reactor fuel channels and reactivity control devices. The methods is based on one- and two dimensional transport calculations with the WIMS-AECL lattice cell code, SPH homogenization, and three-dimensional flux calculations with finite-difference diffusion theory using the MULTICELL code. Simulations of Wolsung 1 Phase-B commissioning measurements and Point Lepreau restart tests have been performed, as a part of the program to validate WIMS-AECL lattice cell calculations for application to CANDU reactor simulations in RFSP. The incremental cross section properties of the Wolsung 1 and ...
1997-06-01
International Nuclear Information System (INIS)
Experimental and theoretical work have been carried out to find out the effects of heat transfer augmentation on two-phase flow instabilities in a single channel system. The effect of mass flow rate, heat input and inlet subcooling on the system behavior is studied using Freon-11 as the test fluid and six different heater surfaces are tested at various experimental conditions. Experimental evaluation is done using the steady-state pressure drop versus mass flow rate curves along with the curves of additional inlet pressure drop required to stabilize the system during the oscillations, and tables generated using the experimental data. Homogeneous equilibrium flow model and finite differences are used in theoretical analysis. An empirical relationship is obtained to calculate the pressure drop across the exit restriction. Characteristic equation for the system is found by linearizing the dynamic equations of the overall system and analyzed to ...
1985-11-01
Solid oxide fuel cell simulation and design optimization with numerical adjoint techniques
This dissertation reports on the application of numerical optimization techniques as applied to fuel cell simulation and design. Due to the "multi-physics" inherent in a fuel cell, which results in a highly coupled and non-linear behavior, an experimental program to analyze and improve the performance of fuel cells is extremely difficult. This program applies new optimization techniques with computational methods from the field of aerospace engineering to the fuel cell design problem. After an overview of fuel cell history, importance, and classification, a mathematical model of solid oxide fuel cells (SOFC) is presented. The governing equations are discretized and solved with computational fluid dynamics (CFD) techniques including unstructured meshes, non-linear solution methods, numerical derivatives with complex variables, and sensitivity analysis with adjoint methods. Following the validation of the fuel cell model in 2-D and 3-D, the results of the sensitivity analysis are ...
2008-01-01
Reflection tomography using finite element method ray tracing
Energy Technology Data Exchange (ETDEWEB)
The seismic reflection exploration technique which is one of the geophysical methods for oil exploration became effectively to image the subsurface structure with rapid development of computer. As a tool to perform seismic inversion, seismic forward modeling program using ray tracing should be developed. In this study, we have developed the algorithm that is to calculate the travel time of the complex geological structure using ray tracing by subdividing the geologic model into triangular element (finite element) having the constant velocity. We can analytically calculate Jacobian with some information by this current ray tracing. With this Jacobian, we will develop new algorithm which is to obtain geological properties and to image the subsurface. Since the FEM (Finite Element Method) ray tracing we have developed goes well the inverse velocities structure, we can apply the inversion problem to complex geological model. This report consists of 3 articles. 1) Development of 2-D ...
1996-12-01
Prestack wave equation reverse-time migration of VSP data
Energy Technology Data Exchange (ETDEWEB)
Based on the wave equation theory, the authors use finite difference method to derive wave-field extrapolation formula, by which both wavefield modeling and migration can be made. The migration program VSP MIG that was developed using P-wave equation may be used to migrate both surface seismic data and VSP data. Besides, it may be applied to the migrations of P-wave data and P-SV wave data, in which case the original velocity model should be modified properly. The migrated image times for P-wave and P-SV wave respectively are absolutely equal. The migrations of upgoing P-wave and upgoing P-SV wave which have been separated can offer migrated P-wave section and migrated S-wave section. Thus, one can get more information than from simultaneous migration of P-wave and S-wave data.
1989-01-01
Energy Technology Data Exchange (ETDEWEB)
We intend to solve equations governing turbulent plane-vertical isotherm and non isotherm jets by taking into account inflow conditions at the exit of the nozzle. The analysis is focused on the influence of these conditions on this type of flow. Two cases are considered (uniform and parabolic velocity and temperature profiles). A finite difference scheme is developed to solve the governing equations. This numeric model allows us to show that the region of fully developed regime begins much nearer the nozzle for the turbulent case than for the laminar flow case. Indeed, the turbulence increases the mixing between the incoming gas from the nozzle and the ambient fluid, and consequently the size of the potential core zone decreases. The results are compared to other works introducing mathematical variables based on the energy conservation for the case of the mixed convection and the momentum conservation for the forced convection, which allows the ...
1999-11-01
New methods for electronic structure calculations
Energy Technology Data Exchange (ETDEWEB)
This thesis consists of five chapters, each of which is a self-contained unit. The first chapter overviews methods for electronic-structure calculations. Chapter 2 introduces a new method to generate a rapidly converging configuration expansion. The approach iteratively combines (1) a least-squares fitting of a configuration expansion to a many-body wave function with (2) a transfer-matrix method for projecting out the ground state. Results are shown to be equivalent to multiconfiguration Hartree-Fock. Results from test calculations are given for a simple finite difference model of the helium atom. In Chapter 3 the use of the finite-element method in electronic structure calculations is discussed. Chapters 4 and 5 discuss developments in Monte Carlo methods based on Hubbard-Stratonovich transformations. Chapter 4 introduces a canonical ensemble formulation of the method, which is more appropriate than the usual grand canonical formulation for ...
1988-01-01
Integrated EM & Thermal Simulations with Upgraded VORPAL Software
Energy Technology Data Exchange (ETDEWEB)
Nuclear physics accelerators are powered by microwaves which must travel in waveguides between room-temperature sources and the cryogenic accelerator structures. The ohmic heat load from the microwaves is affected by the temperature-dependent surface resistance and in turn affects the cryogenic thermal conduction problem. Integrated EM & thermal analysis of this difficult non-linear problem is now possible with the VORPAL finite-difference time-domain simulation tool. We highlight thermal benchmarking work with a complex HOM feed-through geometry, done in collaboration with researchers at the Thomas Jefferson National Accelerator Laboratory, and discuss upcoming design studies with this emerging tool. This work is part of an effort to generalize the VORPAL framework to include generalized PDE capabilities, for wider multi-physics capabilities in the accelerator, vacuum electronics, plasma processing and fusion R&D fields, and we will also discuss user ...
2011-03-01
Energy Technology Data Exchange (ETDEWEB)
The effect of transverse magnetic field parameter (Hartmann number, Ha) Reynolds number (Re) and Prandtl number (Pr) on the mixed convection flow past a semi-infinite vertical porous plate in a non-Darcian porous medium with variable viscosity and porosity, viscous dissipation and fluid-solid thermal conductivity ratio in the presence of plate transpiration (lateral mass flux) is investigated theoretically and numerically using Keller`s implicit finite difference scheme. It is shown that the Harmann number acts as a retarding force and increases the momentum boundary layer thickness, analogous to the flow against a positive pressure gradient, simultaneously decreasing local skin friction (shear stress). The heat transfer rate is however enhanced by the magnetic field (for positive values of the Eckert number) since the fluid is heated and temperature gradients become reduced between the fluid and the plate, with important potential applications ...
1997-01-01
Correlations for the yearly or seasonally optimum salt-gradient solar pond in Greece
Energy Technology Data Exchange (ETDEWEB)
Simple correlations and corresponding nomographs are presented, which express the maximum useful heat received from salt-gradient solar ponds throughout the year or during a specified season of the year, and the corresponding optimum depth of the nonconvective zone in terms of the thickness of the upper convective zone and the temperature under which the maximum useful heat is received. The correlations are valid for the Athens (Greece) area or for regions with a similar climate, because solar radiation and ambient temperature values for Athens have been employed, obtained by a statistical process of hourly measurements over a period of about 20 years. For other climates, it is easy to develop similar correlations using the same methodology, Development of the proposed correlations is based on a method, which simulates the transient operation of the salt-gradient pond using finite-differences, and calculates the useful heat received hourly along the typical year. ...
1993-05-01
Energy Technology Data Exchange (ETDEWEB)
Much of the current technology on 3D-prestack depth migration is based on the Kirchhoff method and the traveltime information needed for this purpose is obtained by the ray-theory approximation. The main motivation driving these options is to reduce computations but it can also have an adverse effect on the quality of results. Experience with real data indicates that the corresponding poststack data migrated by the reverse-time method often yields superior results. Considerable computational efficiency can be obtained by employing nonuniform coarse grid spacing during the reverse-time migration. We followed this strategy for developing a parallel scheme for doing 3D- prestack depth migration. However, we worked only on a small subset of shot gathers whose determined by the Quasi-Monte Carlo method. The results were indistinguishable from those obtained by using a much larger population of uniformly spaced shot gathers. We present examples of migration of 3-D data. (author)
1997-07-01
International Nuclear Information System (INIS)
The COSA II (computer codes for salt) benchmark problem has been pursued with the ADINA (Automatic Dynamic Incremental Nonlinear Analysis) program code. With the use of this, the code should be validated by means of experimental data and the ability to reproduce real-life calculation results of the KfK (Kernforschungszentrum Karlsruhe/Nuclear Research Center in Karlsruhe) should be proven. A successful validation of the code then forms the foundation stone for the ability to use different calculation problems in the final (ultimate) storage. This also accompanies the consequent reaction of replacing the STEALTH (Solids and Thermal Hydraulics Code for EPRI Adapted from LAGRANGE TOODY and HEMP) program which has a number of program-specific weaknesses compared to the ADINA computer code. In order to reproduce the approximate values from the KfK, the same values have been used. Differences were evident in the discretion and the selection of the initial values for ...
Space charge field in a FEL with axially symmetric electron beam
Energy Technology Data Exchange (ETDEWEB)
Nonlinear two-dimensional theory of the space charge of an axially symmetric electron beam propagating in combined right-hand polarized wiggler and uniform axial guide fields in a presence of high-frequency electromagnetic wave is presented. The well-known TE{sub 01} mode in a cylindrical waveguide for the model of radiation fields and paraxial approximation for the wiggler field are used. Space charge field components are written in the Lagrange coordinates by the twice averaged Green`s functions of two equally charged infinitely thin discs. For that {open_quotes}compensating charges{close_quotes} method is applied in which an electron ring model is substituted by one with two different radii and signs discs. On this approach the initial Green`s functions peculiarities are eliminated and all calculations are considerably simplified. Coefficients of a twice averaged Green`s function expansion into a Fourier series are obtained by use of corresponding expansion ...
1995-12-31
Energy Technology Data Exchange (ETDEWEB)
After having recalled the formal convergence of the semi-classical multi-species Boltzmann equations toward the multi-species Euler system (i.e. mixture of gases having the same velocity), we generalize to this system the closure relations proposed by B. Despres and by F. Lagoutiere for the multi-components Euler system (i.e. mixture of non miscible fluids having the same velocity). Then, we extend the energy relaxation schemes proposed by F. Coquel and by B. Perthame for the numerical resolution of the mono-species Euler system to the multi-species isothermal Euler system and to the multi-components isobar-isothermal Euler system. This allows to obtain a class of entropic schemes under a CFL criteria. In the multi-components case, this class of entropic schemes is perhaps a way for the treatment of interface problems and, then, for the treatment of the numerical mixture area by using a Lagrange + projection scheme. Nevertheless, we have to find a good projection ...
2001-07-01
Implementing Explicit and Finding Implicit Sharing in Embedded DSLs
Aliasing, or sharing, is prominent in many domains, denoting that two differently-named objects are in fact identical: a change in one object (memory cell, circuit terminal, disk block) is instantly reflected in the other. Languages for modelling such domains should let the programmer explicitly define the sharing among objects or expressions. A DSL compiler may find other identical expressions and share them, implicitly. Such common subexpression elimination is crucial to the efficient implementation of DSLs. Sharing is tricky in embedded DSL, since host aliasing may correspond to copying of the underlying objects rather than their sharing. This tutorial summarizes discussions of implementing sharing in Haskell DSLs for automotive embedded systems and hardware description languages. The technique has since been used in a Haskell SAT solver and the DSL for music synthesis. We demonstrate the embedding in pure Haskell of a simple DSL with a language form for ...
2011-01-01
Explicit representations of problem-solving strategies to support knowledge acquisition
Energy Technology Data Exchange (ETDEWEB)
Role-limiting approaches support knowledge acquisition (KA) by centering knowledge base construction on common types of tasks or domain-independent problem-solving strategies. Within a particular problem-solving strategy, domain-dependent knowledge plays specific roles. A KA tool then helps a user to fill these roles. Although role-limiting approaches are useful for guiding KA, they are limited because they only support users in filling knowledge roles that have been built in by the designers of the KA system. EXPECT takes a different approach to KA by representing problem-solving knowledge explicitly, and deriving from the current knowledge base the knowledge gaps that must be resolved by the user during KA. This paper contrasts role-limiting approaches and EXPECT`s approach, using the propose-and-revise strategy as an example. EXPECT not only supports users in filling knowledge roles, but also provides support in making other modifications to the knowledge base, ...
1996-12-31
British Library Electronic Table of Contents (United Kingdom)
This study quantifies the external water footprint of the Netherlands by partner country and import product and assesses the impact of this footprint by contrasting the geographically-explicit water footprint with water scarcity in the different parts of the world. The total water footprint of the Netherlands is estimated to be about 2300 m^3/year/cap, of which 67% relates to the consumption of agricultural goods, 31% to the consumption of industrial goods, and 2% to domestic water use. The Dutch water footprint related to the consumption of agricultural goods, is composed as follows: 46% related to livestock products; 17% oil crops and oil from oil crops; 12% coffee, tea, cocoa and tobacco; 8% cereals and beer; 6% cotton products; 5% fruits; and 6% other agricultural products. About 11% o...
2009-01-01
On the Critical Coupling for Kuramoto Oscillators
The celebrated Kuramoto model captures various synchronization phenomena in biological and man-made dynamical systems of coupled oscillators. It is well-known that there exists a critical coupling strength among the oscillators at which a phase transition from incoherency to synchronization occurs. This paper features three contributions. First, we characterize and distinguish the different notions of synchronization used throughout the literature and formally introduce the concept of phase cohesiveness as an analysis tool and performance index for synchronization. Second, we review the vast literature providing necessary, sufficient, implicit, and explicit estimates of the critical coupling strength in the finite and infinite-dimensional case. Finally, we present the first explicit necessary and sufficient condition on the critical coupling strength to achieve synchronization in the finite-dimensional Kuramoto model for an arbitrary ...
2010-01-01
International Nuclear Information System (INIS)
Casimir operators and the Cartan subalgebra are used to construct the scalar superfields in 10-dimensions. In massless case it is shown that the scalar superfield contains two irreducible pieces, one bosonic and one fermionic. The bosonic one contains the supergravity multiplet. Supersymmetric version of the Cartan subalgebra is used to obtain the explicit expressions of the irreducible superfields. In massive case the scalar superfield contains two bosonic and one fermionic irreducible components. It is shown explicitly that the one of the bosonic pieces reduces to the above mentioned massless bosonic piece containing the supergravity multiplet in the massless limit. Supersymmetric generators corresponding to the root vectors of the Lie algebra are found and used with the Cartan subalgebra to construct the irreducible scalar superfields. Finally this method is also applied to the 4-dimensional case and as a result the Transverse Vector ...
1764-01-01
Feynman rules in N=2 projective superspace. II. Massive hypermultiplets
International Nuclear Information System (INIS)
For pt.I see ibid., vol.516, p.426-48, 1998. Manifest N=2 supersymmetric hypermultiplet mass terms can be introduced in the projective N=2 superspace formalism. In the case of complex hypermultiplets, where the gauge covariantized spinor derivatives have an explicit representation in terms of gauge prepotentials, it is possible to interpret such masses as vacuum expectation values of an Abelian vector multiplet. The duality transformation that relates the N=2 off-shell projective description of the hypermultiplet to the on-shell description involving two N=1 chiral superfields allows us to obtain the massive propagators of the N=1 complex linear fields in the projective hypermultiplet. The N=1 massive propagators of the component superfields in the projective hypermultiplet suggest a possible ansatz for the N=2 massive propagator, which agrees with an explicit calculation in N=2 superspace. (orig.).
Advanced distillation curve measurements for corrosive fluids: Application to two crude oils
British Library Electronic Table of Contents (United Kingdom)
We have recently introduced several important improvements in the measurement of distillation curves for complex fluids. The modifications include a composition-explicit data channel for each distillate fraction (for both qualitative and quantitative analysis) and corrosivity assessment of each distillate fraction. The composition-explicit information is achieved with a new sampling approach that allows precise qualitative as well as quantitative analyses of each fraction, on the fly. We have applied the new method to a variety of fluids, including simple n-alkanes, rocket propellant, gasoline, jet fuels, and a hydrocarbon fluid made corrosive with dissolved hydrogen sulfide. In the current contribution, we present the application of the advanced distillation curve method to two samples of...
2008-01-01
The generating function of amplitudes with N twisted and M untwisted states
We show that the generating function of all amplitudes with N twisted and M untwisted states, i.e. the Reggeon vertex for magnetized branes on R^2 can be computed once the correlator of N non excited twisted states and the corresponding Green function are known and we give an explicit expression as a functional of the these objects
2011-01-01
Tachyons: may they have a role in elementary particle physics
International Nuclear Information System (INIS)
The possible role of space-like objects in elementary particle physics (and in quantum mechanics) is reviewed and discussed, mainly by exploiting the explicit consequences of the peculiar relativistic mechanics of Tachyons. Particular attention is paid: (i) to tachyons as the possible carriers of interactions; (ii) to the possibility of ''vacuum decays'' at the classical level; (iii) to a Lorentz-invariant bootstrap model; (iv) to the apparent shape of the tachyonic elementary particles and its possible connection with the de Broglie wave-particle dualism. (author).
Energy Technology Data Exchange (ETDEWEB)
We discuss condensations of closed string tachyons localized in compact spaces. Time evolution of an on-shell condensation is naturally related to the worldsheet RG flow. Some explicit tachyonic compactifications of Type II string theory is considered, and some of them are shown to decay into supersymmetric theories known as the little string theories.
2005-05-01
Surface brightness, galaxy evolution and cosmology
Energy Technology Data Exchange (ETDEWEB)
The Thomsen and Frandsen (1983) cosmological test employing the observed correlation between elliptical galaxy surface brightness and scale size in order to avoid problems due to dynamical evolution is presently extended to make use of arbitrary galaxy samples. The explicit dependence on stellar evolution is also demonstrated. On the assumption that this evolution is calculable, an equation for the deceleration parameter entirely in terms of observables is derived. The test is applied to two available samples. 25 references.
1985-01-01
International Nuclear Information System (INIS)
An explicite PN solution of the multi-dimensional homogeneous neutron transport equation is given by expanding the angular flux into a series of geometry-independent spherical harmonics operators. An algorithm is developed for representing the spherical harmonic operators in orthogonal curvilinear coordinates. The general formulae are applied to two-dimensional spherical geometry; detailed P3 formulae are given. (orig.).
Sources of X-rays in school; Zrodla promieniowania rentgenowskiego w szkole
Energy Technology Data Exchange (ETDEWEB)
Exposure of teachers of physics to X-radiation emitted by instruments which are used for demonstrating rarefied gas discharges during physics lessons at secondary schools is discussed. The measurements performed provide an explicit evidence that an effective, annual exposure dose under the most unfavorable conditions does not exceed admissible levels according to Polish regulations pertaining to persons non-occupationally exposed to ionizing radiation. (author). 3 refs,1 fig.
1996-12-31
Scalar fields in the dimensional reduction scheme for symmetric spaces
Energy Technology Data Exchange (ETDEWEB)
The authors study the general features of the dimensional reduction scheme for multi-dimensional spaces of the type M/sup 4/ x S/R, S/R being a symmetric coset space. The properties of the scalar potentials of the reduced theories are investigated and an effective method of explicit calculation of these potentials is elaborated. They consider also a wide class of embeddings of Lie subalgebras into simple Lie algebras resulting in reduced theories of physical interest.
1989-01-01
On energetics of a charged black hole relative to tachyons
International Nuclear Information System (INIS)
The author considers the extraction of the electromagnetic energy from a charged black hole through charged tachyons. The 'generalized ergosphere' relative to tachyons is defined. It turns out that it is possible to extract energy with 100% efficiency (through reversible process) if the charged source (whether black hole or otherwise) is of finite size. Two explicit examples of energy extraction are discussed. (Auth.).
Nuclear electromagnetic currents from chiral EFT
International Nuclear Information System (INIS)
Using the method of unitary transformation in combination with chiral effective field theory we derive the pion exchange contributions to the two-nucleon electromagnetic current. A formal definition of the current operator in this scheme and the power counting is presented. We discuss the implications of additional unitary transformations that have to be present to ensure the renormalizability of the one-pion exchange current. Further, we give explicit and compact results for the current in coordinate-space.
2010-03-15
Infinitesimal symmetries and conservation laws of the DNLSE hierarchy and the Noether's theorem
Energy Technology Data Exchange (ETDEWEB)
The hierarchy of integrable nonlinear equations associated with the quadratic bundle is considered. The expressions for the solution of linearization of these equations and their conservation law in the terms of solutions of corresponding Lax pairs are found. It is shown for the first member of the hierarchy that the conservation law is connected with the solution of linearized equation due to the Noether's theorem. The local hierarchy and three nonlocal ones of the infinitesimal symmetries and conservation laws explicitly expressed through the variables of the nonlinear equations are derived. (author)
2007-08-15
Geometrically relating momentum cut-off and dimensional regularization
The $\\beta$ function for a scalar field theory describes the dependence of the coupling constant on the renormalization mass scale. This dependence is affected by the choice of regularization scheme. I explicitly relate the $\\beta$-functions of momentum cut-off regularization and dimensional regularization on scalar field theories by a gauge transformation using the Hopf algebras of the Feynman diagrams of the theories.
2011-01-01
Explicit equation for particle settling velocities in solid-liquid systems
Energy Technology Data Exchange (ETDEWEB)
Zanker has recently presented nomographs for determining particle settling velocities in solid-liquid systems. These nomographs were based on the general correlations developed by Barnea and Mizrahi and Barnea and Mednick. This work presents an equation directly computing particle settling velocities, eliminating the uncertainty associated with nomographs.
1981-11-01
Effective fracture toughness of microcracked materials
Energy Technology Data Exchange (ETDEWEB)
Explicit analytical formulas are derived for the stress intensity factors at the tips of a main crack and of a microcrack for the two-dimensional case of a collinear microcrack. This configuration is used to derive an estimate of the toughness degradation due to microcracks linking up with an advancing main crack. The implications of this estimate for theoretical predictions of the toughening due to stress-induced microcracking are discussed.
1986-03-01
Depth migration in transversely isotropic media with explicit operators
Energy Technology Data Exchange (ETDEWEB)
The author presents and analyzes three approaches to calculating explicit two-dimensional (2D) depth-extrapolation filters for all propagation modes (P, SV, and SH) in transversely isotropic media with vertical and tilted axis of symmetry. These extrapolation filters are used to do 2D poststack depth migration, and also, just as for isotropic media, these 2D filters are used in the McClellan transformation to do poststack 3D depth migration. Furthermore, the same explicit filters can also be used to do depth-extrapolation of prestack data. The explicit filters are derived by generalizations of three different approaches: the modified Taylor series, least-squares, and minimax methods initially developed for isotropic media. The examples here show that the least-squares and minimax methods produce filters with accurate extrapolation (measured in the ability to position steep reflectors) for a wider range of propagation angles ...
1994-12-01
Cut diagrams for high energy scatterings
International Nuclear Information System (INIS)
A new approach is introduced to study QCD amplitudes at high energy and comparatively small momentum transfer. Novel cut diagrams, representing the resummation of Feynman diagrams, are used to simplify the calculation and to avoid delicate cancellations encountered in the usual approach. An explicit calculation to the sixth order is carried out to demonstrate the advantage of cut diagrams over Feynman diagrams. copyright 1996 The American Physical Society.
Canonical Gravity with Fermions
Canonical gravity in real Ashtekar-Barbero variables is generalized to allow for fermionic matter. The resulting torsion changes several expressions in Holst's original vacuum analysis, which are explicitly displayed here. This in turn requires adaptations to the known canonical (loop) quantization of gravity coupled to fermions, which is discussed on the basis of the classical analysis.
2007-01-01
Calculation of atomic spontaneous emission rate in 1D finite photonic crystal with defects
We derive the expression for spontaneous emission rate in finite one-dimensional photonic crystal with arbitrary defects using the effective resonator model to describe electromagnetic field distributions in the structure. We obtain explicit formulas for contributions of different types of modes, i.e. radiation, substrate and guided modes. Formal calculations are illustrated with a few numerical examples, which demonstrate that the application of effective resonator model simplifies interpretation of results.
2009-01-01
International Nuclear Information System (INIS)
The time-dependent probability distribution of neutrons, precursors and detectrons in a space-independent, low-power, multiplying assembly with a source is developed in the multigroup approximation as forward and backward Kolmogorov equations. The relationship between these as adjoint equations is made explicit in a tensor notation and the equations developed in the generating-function formalism. (author).
1982-01-01
A High-Wavenumber Viscosity for High-Resolution Numerical Methods
Energy Technology Data Exchange (ETDEWEB)
Numerical simulations of compressible flows are commonly based on the Euler equations when effects of viscosity are thought to be negligible. These equations admit singular solutions, even in cases where the initial and boundary conditions are smooth. So-called ''Euler solvers'' rely on numerical dissipation, explicitly or implicitly present in the scheme, to regularize the problem, such that physical solutions are selected.
2003-02-19
X-ray ablation measurements and modeling for ICF applications
Energy Technology Data Exchange (ETDEWEB)
X-ray ablation of material from the first wall and other components of an ICF (Inertial Confinement Fusion) chamber is a major threat to the laser final optics. Material condensing on these optics after a shot may cause damage with subsequent laser shots. To ensure the successful operation of the ICF facility, removal rates must be predicted accurately. The goal for this dissertation is to develop an experimentally validated x-ray response model, with particular application to the National Ignition Facility (NIF). Accurate knowledge of the x-ray and debris emissions from ICF targets is a critical first step in the process of predicting the performance of the target chamber system. A number of 1-D numerical simulations of NIF targets have been run to characterize target output in terms of energy, angular distribution, spectrum, and pulse shape. Scaling of output characteristics with variations of both target yield and hohlraum wall thickness are also described. Experiments have been ...
1996-09-01
Unsteady-state flow of methane and water in coal beds
Energy Technology Data Exchange (ETDEWEB)
A hypothesis on the transition of a fully water-saturated coal bed methane reservoir from single-phase flow through unsaturated flow to the two-phase flow regime was developed and substantiated using field data. The concept of parameter measurement windows was introduced to obtain the required reservoir parameters from flow testing. A new empirical relationship for capillary pressure vs. saturation, referred to as the Fermi distribution, was developed and shown to fit the capillary pressure-saturation data for coal, sandstone, and other consolidated and unconsolidated materials. New approximate analytic solutions for flow to a sink from an infinite coal bed reservoir with desorbing methane have been developed and verified using numerical models for the following cases: (1) single-phase gas flow with desorption; (2) single-phase water flow in fully water-saturated and unsaturated flow regimes. A numerical model based on finite-difference method for simultaneous flow ...
1987-01-01
Thermal performance of phase change wallboard for residential cooling application
Energy Technology Data Exchange (ETDEWEB)
Cooling of residential California buildings contributes significantly to electrical consumption and peak power demand mainly due to very poor load factors in milder climates. Thermal mass can be utilized to reduce the peak-power demand, downsize the cooling systems, and/or switch to low-energy cooling sources. Large thermal storage devices have been used in the past to overcome the shortcomings of alternative cooling sources, or to avoid high demand charges. The manufacturing of phase change material (PCM) implemented in gypsum board, plaster or other wall-covering material, would permit the thermal storage to become part of the building structure. PCMs have two important advantages as storage media: they can offer an order-of-magnitude increase in thermal storage capacity, and their discharge is almost isothermal. This allows the storage of high amounts of energy without significantly changing the temperature of the room envelope. As heat storage takes place inside the building, where ...
1997-04-01
International Nuclear Information System (INIS)
A theoretical model was developed to simulate this phenomenon, under some specific assumptions. The model simulates: the generation and decay of radon within the soil; its transport throughout the soil due to diffusion and convection induced by the pressure disturbance applied at a crack in the basement; its entrance into the house through the crack; and the resultant indoor radon concentration. The most important assumptions adopted in the model were: a steady-state condition; a house with a basement; a geometrically well-defined crack at the wall-floor joint in the basement; and a constant negative pressure applied at the crack in relation to the outside atmospheric pressure. Two three-dimensional finite-difference computer programs were written to solve the mathematical equations of the model. The first program, called PRESSU, was used to calculate: the pressure distribution within the soil as a result of the applied disturbance pressure at the crack; and the ...
1997-10-16
Energy Technology Data Exchange (ETDEWEB)
Many thermoluminescent materials has been developed and used for photon personal dosimetry but no one has all desired characteristics alone. These characteristics include robustness, high sensitivity, energy photon independence, large range of photon energy detection, good reproducibility, small fading and simple glow curve with peaks above 150 deg C. Calcium Sulfate Dysprosium doped (CaSO{sub 4}:Dy) phosphor Thermoluminescent Dosimeter (TLD) has been used by many laboratories, mainly in Brazil and India. Another interesting phosphor is Calcium Fluoride (CaF{sub 2}). These phosphor advantages begin to be more required and its disadvantages have became more apparent, in a global market more and more competitive. These phosphors are used in environmental and area monitoring, once they present more sensibility than other phosphors, like LiF:Mg. Theirs mainly disadvantage is a strong energetic dependence response, which must be corrected for theirs application in the field, where photon ...
2006-07-01
Energy Technology Data Exchange (ETDEWEB)
This report documents a one-dimensional numerical model, SAMFT1D, developed to simulate single-phase and multiphase fluid flow and solute transport in variably saturated porous media. The formulation of the governing equations and the numerical procedures used in the code for single-phase and multiphase flow and transport are presented. The code is constructed to handle single-phase as well as two or three-phase flow conditions using two integrated sets of computational modules. The fully implicit scheme is used in the code for both single-phase and multiphase flow simulations. Either the Crank-Nicholson scheme or the fully implicit scheme may be used in the transport simulation. The single-phase modules employ the Galerkin and upstream weighted residual finite element techniques to model flow and transport of water (aqueous phase) containing dissolved single-species contaminants concurrently or sequentially, and include the treatment of various boundary conditions and physical ...
1991-09-01
Energy Technology Data Exchange (ETDEWEB)
A theoretical background of crosshole reflection imaging with pre stack, multi-mode reverse time migration has been investigated. The crosshole migration has been carried out by solving two-way scalar wave equation with finite difference approximation. The wave field extrapolation was performed backward in time, with the imaging condition applied at each extrapolation time step to create partial image of the reflector. The complete migrated image was obtained by adding all the partial images from each time step. With multi-component crosshole data available, it is possible to separate the recorded data into primary P- and S-wave field. Separating wave field before migration is helpful in enhancing migrated image quality since only one wave mode can be dealt at a time with the scalar wave equation. It is possible to migrate mixed-mode primary reflections(P-S and S-P) as well as pure mode ones. These multi-mode migrations provide not only image ...
1993-08-31
Numerical simulation of the direct-injection diesel engine under motored and firing conditions
Energy Technology Data Exchange (ETDEWEB)
An Eulerian/Langrangian model was developed to study the two-phase flow in a two-dimensional, axisymmetric, direct-injection diesel engine, under motored and firing conditions. The liquid phase was tracked in a Lagrangian sense using discrete droplet packets (DDP). The gaseous phase was solved using an implicit, iterative, finite-difference, scheme of the control volume variety. Gaseous-phase turbulence was modeled using a compressible two-equation k/epsilon model. Turbulence effects on the liquid phase were modeled by superimposing a randomly oriented turbulent eddy velocity proportional to the root of the turbulent kinetic energy, on each packet's mean relative velocity. Droplet-droplet interactions were neglected. The general knowledge obtained during validation studies, was used to simulate the direct-injection two-stroke diesel engine under motored conditions. The simulation showed that by beginning of the injection, (15/sup 0/ before top dead center: ...
1986-01-01
Numerical ray tracing of electrons in different 3D fringing fields of spherical deflectors
This paper presents the results of numerical investigations on electron trajectories in different 3D fringing fields of spherical defectors. We start with a chosen set of geometrical values for the electrodes such as radius of curvature, electrode gap, angle between entrance and exit face. We do a first numerical calculation of the spatial potential distribution using a 3D electron optics program based on the finite difference method (FDM). Ray tracing through this potential distribution is then performed using another program, which also gives the potential and the electric field strength along the curved beam path. In the next step we estimate the position and the geometry of the field clamp using Herzog's formulae assuming the entrance of a parallel plate condenser to be a good approximation for the entrance of the spherical deflector in first order. For the consecutive numerical recalculations of the real potential distribution of the ...
1999-11-01
Multiphase treatment of ODTX in HMX spheres
Energy Technology Data Exchange (ETDEWEB)
What will be discussed in this report represents a framework upon which multiphase and other real physical effects can be built. Chemical models of increasing complexity are envisioned and this methodology can provide a tool for evaluating new ideas against known experimental data. The recent work to be reported here addresses the multiphase issue of temperature deviation between phases undergoing chemical and heat transport processes. Modeling of the LLNL ODTX experiment will be performed with FLUENT, a commercially available computational fluid dynamics (CFD) code. FLUENT solves flows in 2D or 3D in Cartesian, cylindrical, or general curvilinear coordinates, with steady-state of fully time-dependent analysis. Multiphase flows in which two or more continuous phases are present can be solved with arbitrary volumetric sources of heat, mass, momentum, and chemical species applied through user-defined FORTRAN subroutines. FLUENT models these of phenomena by solving the conservation ...
1997-12-22
Liquid-metal flow in a sharp elbow in a uniform transverse magnetic field
In the self-cooling blankets of the Tokamak fusion reactor, a liquid metal, namely liquid lithium, is pumped through a system of ducts to transfer heat and capture neutrons. One of the blanket designs proposed in Argonne National Laboratory's Blanket Comparison and Selection study uses a combination of poloidal and toroidal ducts in order to maximize heat transfer while minimizing net pressure drop. In the design, the poloidal and toroidal ducts meet at sharp, abrupt corners. They were modelled as two identical, straight, semi-infinite, thin-walled, rectangular ducts with 45{degree} miters and joined at a 90{degree} angle in the plane of a strong, uniform magnetic field. While in the toroidal containment vessel (i.e. the blanket), the liquid lithium is subjected to a large electromagnetic body force due to the presence of a strong magnetic field. This body force so dominates the flow as to make the inertial and viscous forces negligible everywhere, except in thin boundary or ...
1989-01-01
Imaging of offset VSP data with an elastic iterative migration scheme
Energy Technology Data Exchange (ETDEWEB)
VSP data are usually acquired in order to obtain high-resolution images of complex structures in reservoirs and near boreholes. The authors present an elastic iterative migration scheme which has few limitations regarding the complexity of the geology, and where the macromodel for both P- and S-wave velocities is automatically improved and updated at each iteration. They avoid wavefield separation (up/down and P/S) and the simplifying assumptions of small dips underlying most such methods. The migration scheme is based on elastic inversion theory. The wavefield extrapolation is based on a high-order, coarse-grid, finite-difference solution to the elastic two-way wave equation. At each iteration, the macromodel is updated using a gradient method, in which the gradient is computed by correlation of forward-modelled fields with back-propagated residual fields. The first iteration of the migration scheme is equivalent to elastic reverse-time migration with an imaging ...
1997-03-01
Emergency ampacities of direct buried three phase underground cable systems
Energy Technology Data Exchange (ETDEWEB)
A thermal model is formulated for the determination of the emergency transient ampacities of a three phase, horizontally spaced underground cable system. The model is based on a finite difference heat transfer analysis and the principle of superposition which permits the summation of the temperature rise of individual heat sources to eventually provide for the temperature rise of a multiple cable system. The thermal model is transformed into a complex computer program that is capable of calculating the real time cable temperature for any change in conductor current. The program is used to calculate the response of conductor temperatures caused by both ramp and step changes in currents for four conductor sizes between 4/0 and 1500 kcmil. The amount of time to reach conductor temperatures of 100/sup 0/C and 120/sup 0/C from a steady conductor temperature of 75/sup 0/C is calculated for the different conductor sizes and for different emergency ...
1983-07-01
ERRICCA radon model intercomparison exercise
Energy Technology Data Exchange (ETDEWEB)
Numerical models based on finite-difference or finite-element methods are used by various research groups in studies of radon-222 transport through soil and building materials. Applications range from design of radon remediation systems to more fundamental studies of radon transport. To ascertain that results obtained with these models are of good quality, it is necessary that such models are tested. This document reports on a benchmark test organized by the EU project ERRICCA: European Research into Radon in Construction Concerted Action. The test comprises the following cases: 1) Steady-state diffusive radon profiles in dry and wet soils, 2) steady-state entry of soil gas and radon into a house, 3) time-dependent radon exhalation from a building-material sample. These cases cover features such as: soil heterogeneity, anisotropy, 3D-effects, time dependency, combined advective and diffusive transport of radon, flux calculations, and partitioning of radon between ...
1999-04-01
Analysis of enclosed sodium pool fire scenario in sodium fire experimental facility
International Nuclear Information System (INIS)
Liquid sodium is used as coolant in Fast Breeder Reactors (FBR). There is a likelyhood of sodium spillage in ambient air in the Steam Generator Building (SGB) of the FBR plant. Due to high chemical reactivity with oxygen, especially at temperatures greater than 573 K, it catches fire very easily. In order to carryout safety related experimental studies for different modes of sodium fires and to develop suitable mathematical models for the assessment of their consequences, an experimental facility (SFEF, Sodium Fire Experimental Facility) is being setup a IGCAR, Kalpakkam. The SFEF is having a 540 m"3 volume experimental hall. Stainless steel linear will be provided on the inside surfaces of experimental hall walls, ceiling and floor. Analysis has been carried out for enclosed sodium pool fire scenarios in SFEF by using sodium pool fire code SOFIRE II, which estimates the thermal transients like pressure rise, gas temperature rise, cell wall temperature rise and sodium burning rate. ...
2007-04-22
Alloying element diffusion in alloy 800 H (UNS N08810) and alloy 617 (UNS N06617)
Energy Technology Data Exchange (ETDEWEB)
Samples of alloy 800 H (LJNS N08810) and alloy 617 (UNS N06617) were oxidized in a simulated process gas atmosphere. The principal elements of the gaseous environment were 50% H{sub 2}O, 35% H{sub 2} and 5% of CO{sub 2}, CO, and CH{sub 4}, respectively. The experiments were carried out at temperatures from 850 C to 950 C for exposures ranging from 600 to 5,000 hours. After the exposure, the alloy surfaces showed compact oxide layers, which resulted in a depleted zone of oxidizing elements in the substrate layer. Using Electron Probe Microanalysis (EPMA), concentration profiles and the extension of the depleted zone were measured. Calculated theoretical profiles were then fitted to the EPMA-measured profiles, using Fick`s second law of diffusion and resolved by a finite difference method and diffusion coefficients as fitting parameters. The diffusion coefficients for chromium may be represented as D = Do exp.(-Q/RT) with the values Do = 11.4 ...
1996-10-01
Energy Technology Data Exchange (ETDEWEB)
We present an efficient parallel algorithm and its implementation for computing the diagonal of $H^-1$ where $H$ is a 2D Kohn-Sham Hamiltonian discretized on a rectangular domain using a standard second order finite difference scheme. This type of calculation can be used to obtain an accurate approximation to the diagonal of a Fermi-Dirac function of $H$ through a recently developed pole-expansion technique \\cite{LinLuYingE2009}. The diagonal elements are needed in electronic structure calculations for quantum mechanical systems \\citeHohenbergKohn1964, KohnSham 1965,DreizlerGross1990. We show how elimination tree is used to organize the parallel computation and how synchronization overhead is reduced by passing data level by level along this tree using the technique of local buffers and relative indices. We analyze the performance of our implementation by examining its load balance and communication overhead. We show that our implementation ...
2009-09-25
We analytically work out the orbital effects caused by a Rindlertype extra-acceleration ARin which naturally arises in some recent models of modified gravity at large distances. In particular, we focus on the perturbations induced by it on the two-body range {\\rho} and range-rate {\\rho}\\cdot which are commonly used in satellite and planetary investigations as primary observable quantities. The constraints obtained for ARin by comparing our calculations with the currently available range and range-rate residuals for some of the major bodies of the solar system, obtained without explicitly modeling ARin, are 1 - 2 \\times 10-13 m s-2 (Mercury and Venus), 1 \\times 10-14 m s-2 (Saturn), 1 \\times 10-15 m s-2 (Mars), while for a terrestrial Rindler acceleration we have 5 \\times 10-16 m s-2 (Moon). Another approach which could be followed consists of taking into account ARin in re-processing all the available data sets with accordingly modified dynamical models, and ...
2010-01-01
Entanglement Spectrum of a Disordered Topological Chern Insulator
How much information is stored in the ground-state of a system without \\emph{any symmetry} and how can we extract it? This question is investigated by analyzing the behavior of a topological Chern Insulator (CI) in the presence of disorder, with a focus on its entanglement spectrum (EtS) constructed from the ground state. For systems with symmetries, the EtS was shown to contain explicit information revealed by sorting the EtS against the conserved quantum numbers. In the absence of any symmetry, we demonstrate that statistical methods such as the level statistics of the EtS can be equally insightful, allowing us to distinguish when an insulator is in a topological or trivial phase and to map the boundary between the two phases, where EtS becomes entirely delocalized. The phase diagram of a CI is explicitly computed as function of Fermi level ($E_F$) and disorder strength using the level statistics of the EtS and energy spectrum (EnS), ...
2010-01-01
Electroweak Supersymmetry with an Approximate U(1)_PQ
A new predictive framework for supersymmetry at the TeV scale is presented. The \\mu parameter of the MSSM is replaced by (\\lambda S), where S is a singlet field, and the axion becomes a heavy pseudoscalar, G, by adding a mass, m_G, by hand. The explicit breaking of Peccei--Quinn (PQ) symmetry is assumed to be sufficiently weak at the TeV scale that the only observable consequence is the mass m_G. Three models for the explicit PQ breaking are given; but the utility of this framework is that the predictions for all physics at the electroweak scale are independent of the particular model for PQ breaking. The MSSM parameters \\mu, B are replaced by \\lambda, A_\\lambda. Taking the scalar mass m_S^2 to be sufficiently small and \\lambda = A_\\lambda \\sin (2 \\beta) /2. Our framework leads to a theory similar to the MSSM, except that \\mu is predicted and there are light, weakly-coupled states that lie dominantly in the superfield S. The ...
2004-01-01
Energy Technology Data Exchange (ETDEWEB)
Calculations of refrigerant thermal properties are desired to be very fast and stable in cases of simulation of refrigeration system, etc. The traditional method based on equation of state cannot meet such requirement because of unavoidable iterations in calculation. In this paper, a new calculation method for refrigerant thermal properties is presented. Low order implicit polynomial equations are got by using curve-fitting method at first, and then explicit formulae for calculating refrigerant thermal properties quickly are obtained by getting the analytical solution of these implicit equations. Explicit fast calculation formulae for thermal properties of R22 and R407C, covering the saturated temperature of -60{approx}80 {sup o}C and superheat of 0-65 {sup o}C, are presented as examples. The calculation speeds of the formulae of R22 are about 140 times faster than those of REFPROP 6.01 while the formulae of R407C are about 1000 times faster. ...
2005-09-01
[Fritz hauschild (1908-1974) and drug research in the 'German Democratic Republic' (GDR)].
The chemist and pharmacologist Fritz Hauschild developed the sympathomimetic agent Pervitin (metamphetamin) in the 1930s. Not only because of the abuse of the stimulant during the Second World War ("pilot's chocolate") it is one of the most controversial substances in drug history. Nearly forgotten are Hauschild's contributions to build up the drug system in the GDR. Although he was a convinced communist, the skilful pharmacologist gave very early warning of the imminent lack of innovation in the GDR pharmaceutical industry. A letter which he addressed to the Minister of Health, Max Sefrin (born 1913), did not lack explicitness. PMID:15997840
2005-06-01
Vibrations of circular cylindrical shells
An exact solution method for the free vibration problem of thin circular cylindrical shells is presented. The differential equations of motion are solved directly with the use of simple Fourier series as the modal displacement functions. Stokes' transformation is exploited to obtain correct series expressions for the derivatives of the Fourier series. From this method an explicit expression of the exact frequency equation can be obtained for any kind of boundary conditions. The accuracy of the present method is checked against available data. The proposed method is then used to find the modal characteristics of the thermal liner of the Fast Test Reactor (FTR). The numerical results obtained are compared with finite element method solutions. (auth)
1975-07-01
Unconstrained energy functionals for electronic structure calculations
Energy Technology Data Exchange (ETDEWEB)
The performance of conjugate gradient schemes for minimizing unconstrained energy functionals in the context of electronic structure calculations is studied. The unconstrained functionals allow a straightforward application of conjugate gradients by removing the explicit orthonormality constraints on the quantum-mechanical wave functions. However, the removal of the constraints can lead to slow convergence, in particular when preconditioning is used. The convergence properties of two previously suggested energy functionals are analyzed, and a new functional is proposed, which unifies some of the advantages of the other functionals. A numerical example confirms the analysis.
1998-04-01
Towards the entropy of gravity time-dependent models via the Cardy-Verlinde formula
International Nuclear Information System (INIS)
For models with several time-dependent components, generalized entropies can be defined. This is shown for the Bianchi type IX model. We first derive the Cardy-Verlinde formula under the assumption that the first law of thermodynamics is valid. This leads to an explicit expression of the total entropy associated with this type of universe. Assuming the validity of the Cardy entropy formula, we obtain expressions for the corresponding Bekenstein, Bekenstein-Hawking and Hubble entropies. We discuss the validity of the Cardy-Verlinde formula and possible extensions of the outlined procedure to other time-dependent models.
2003-07-15
Toward a rule-based biome model
Energy Technology Data Exchange (ETDEWEB)
Current projections of the response of the biosphere to global climatic change indicate as much as 50% to 90% spatial displacement of extratropical biomes. The mechanism of spatial shift could be dominated by either (1) competitive displacement of northern biomes by southern biomes, or (2) drought-induced dieback of areas susceptible to change. The current suite of global biosphere models cannot distinguish between these two processes, thus determining the need for a mechanistically based biome model. The model is in an early stage of development and will require several enhancements, including explicit simulation of potential evapotranspiration, extension to boreal and tropical biomes, a shift from steady-state to transient dynamics, and validation on other continents.
1991-01-01
Theoretical considerations for SRAM total-dose hardening
International Nuclear Information System (INIS)
The theoretical hardness against total dose of the six-transistor SRAM cell is investigated in detail. An explicit analytical expression of the maximum tolerable threshold voltage shift is derived for two cross-coupled inverters. A numerical method is used to explore the hardness of the read and write operations. Both N- and P-channel access transistors designs are considered and their respective advantages are compared. The study points out that the radiation hardness mainly relies on the technology. Results obtained with the very robust Gate-All-Around process are finally presented.
Tachyons: may they have a role in elementary particle physics
Energy Technology Data Exchange (ETDEWEB)
The possible role of space-like objects in elementary particle physics (and in quantum mechanics) is reviewed and discussed, mainly by exploiting the explicit consequences of the peculiar relativistic mechanics of Tachyons. Particular attention is paid: (i) to tachyons as the possible carriers of interactions; (ii) to the possibility of ''vacuum decays'' at the classical level; (iii) to a Lorentz-invariant bootstrap model; (iv) to the apparent shape of the tachyonic elementary particles and its possible connection with the de Broglie wave-particle dualism. (author).
1985-01-01
Symmetrizers and antisymmetrizers for the BMW algebra
Let $n\\in\\mathds{N}$ and $B_n(r,q)$ be the generic Birman-Murakami-Wenzl algebra with respect to indeterminants $r$ and $q$. It is known that $B_n(r,q)$ has two distinct linear representations generated by two central elements of $B_n(r,q)$ called the symmetrizer and antisymmetrizer of $B_n(r,q)$. These generate for $n\\geq 3$ the only one dimensional one sided ideals of $B_n(r,q)$ and generalize the corresponding notion for Hecke algebras of type $A$. In this paper the coefficients of these elements with respect to the graphical basis of $B_n(r,q)$ are determined explicitly.
2011-01-01
Simulation of concrete perforation based on a continuum damage model
Energy Technology Data Exchange (ETDEWEB)
Numerical simulation of dynamic fracture of concrete slabs, impacted by steel projectiles, was carried out in this study. The concrete response was described by a continuum damage model. This continuum damage model was originally developed to study rock fragmentation and was modified in the present study with an emphasis on the post-limit structural response. The model was implemented into a transient dynamic explicit finite element code LS-DYNA2D and the code was then used for the numerical simulations. The specific impact configuration of this study follows the experiment series conducted by Hanchak et al. Comparisons between calculated results and measured data were made. Good agreements were found.
1994-10-01
Separation of variables for the quantum SL(2,R) spin chain
We construct representation of the Separated Variables (SoV) for the quantum SL(2,R) Heisenberg closed spin chain and obtain the integral representation for the eigenfunctions of the model. We calculate explicitly the Sklyanin measure defining the scalar product in the SoV representation and demonstrate that the language of Feynman diagrams is extremely useful in establishing various properties of the model. The kernel of the unitary transformation to the SoV representation is described by the same "pyramid diagram" as appeared before in the SoV representation for the SL(2,C) spin magnet. We argue that this kernel is given by the product of the Baxter Q-operators projected onto a special reference state.
2003-01-01
Quasilinear theory of the 2D Euler equation
Motivated by the numerical investigations of Laval, Dubrulle & Nazarenko (1999), we develop a quasilinear theory of the 2D Euler equation and derive an integro-differential equation for the evolution of the coarse-grained vorticity. This equation respects all the invariance properties of the Euler equation and conserves angular momentum in a circular domain and linear impulse in a channel (as well as in an infinite domain). The explicit energy is not rigorously conserved as it is partly transfered into fine-grained fluctuations but the total energy is conserved. We prove a H-theorem for the Fermi-Dirac entropy and make the connection with statistical theories of 2D turbulence.
2000-01-01
Quantum theory of the interaction of Josephson junctions with non-classical microwaves
Energy Technology Data Exchange (ETDEWEB)
We present a study of the interaction between Josephson junctions in circular superconducting rings and non-classical microwaves, treating both quantum mechanically. A Hamiltonian that describes both inductive and capacitive coupling between the two systems is derived within the external field approximation. Other Hamiltonians which go beyond the external field approximation, and describe explicitly the interaction of the quantum circuit that produces the non-classical microwaves with the Josephson junction circuit, are also presented. A comparison between current experiments which use classical electromagnetic fields and the proposed experiments that use non-classical microwaves, is made. (orig.) With 6 figs., 32 refs.
1997-01-01
On Financial Markets Based on Telegraph Processes
The paper develops a new class of financial market models. These models are based on generalized telegraph processes: Markov random flows with alternating velocities and jumps occurring when the velocities are switching. While such markets may admit an arbitrage opportunity, the model under consideration is arbitrage-free and complete if directions of jumps in stock prices are in a certain correspondence with their velocity and interest rate behaviour. An analog of the Black-Scholes fundamental differential equation is derived, but, in contrast with the Black-Scholes model, this equation is hyperbolic. Explicit formulas for prices of European options are obtained using perfect and quantile hedging.
2007-01-01
Numerical integration of the Feynman Path Integral for radiative transport
International Nuclear Information System (INIS)
The radiative transport problem is cast in integral form using a transport kernel. The transport kernel has an explicit representation in terms of a Feynman Path Integral over all paths between selected points in a volume. This representation is setup in detail. Numerical evaluation of this Path Integral is formulated with a Frenet-Serret based procedure for generating valid random paths, and with a numerical evaluation of the weight for each valid path. Very early sanity checks of a numerical implementation are reported. Approaches to optimization are identified. (authors)
2009-05-03
Modeling Blast and High-Velocity Impact of Composite Sandwich Panels
Analytical models for predicting the deformation and failure of composite sandwich panels subjected to blast and projectile impact loading are presented in this paper. The analytical predictions of the transient deformations and damage initiation in the composite sandwich panels were compared with finite element solutions using ABAQUS Explicit. For the blast model, the predicted transient deformation of the sandwich panel was within 7%of FEA results, while the predicted damage initiation using Hashin's composite failure criteria was about 15%higher than FEA results in most cases. For the high velocity impact model, the predicted transient deformations were within 20%of FEA results.
2009-01-01
Minimal Polynomials of Some Matrices Via Quaternions
This work provides explicit characterizations and formulae for the minimal polynomials of a wide variety of structured $4\\times 4$ matrices. These include symmetric, Hamiltonian and orthogonal matrices. Applications such as the complete determination of the Jordan structure of skew-Hamiltonian matrices and the computation of the Cayley transform are given. Some new classes of matrices are uncovered, whose behaviour insofar as minimal polynomials are concerned, is remarkably similar to those of skew-Hamiltonian and Hamiltonian matrices. The main technique is the invocation of the associative algebra isomorphism between the tensor product of the quaternions with themselves and the algebra of real $4\\times 4$ matrices.
2010-01-01
Linac based free-electron laser
International Nuclear Information System (INIS)
A basic treatment of the principle of the linac-driven free-electron laser (FEL) is given. The first part of the paper describes the FEL in low-gain approximation, and in the second part the high-gain FEL theory is given. The majority of the treatment describes FELs in one dimensional approximation, neglecting effects by diffraction of radiation and by electron beam emittance. Only in the final section a few remarks on these issues are given. The ambition of the paper is by no means any progress in FEL theory but a clear presentation of basic FEL theory concepts with explicit derivation of the formulae from first principles. (orig.)
Inhibition of Two-Photon Absorption in a Four-Level Atomic System with Closed-Loop Configuration
International Nuclear Information System (INIS)
We theoretically investigate the features of two-photon absorption in a coherently driven four-level atomic system with closed-loop configuration. It is found that two-photon absorption can be completely suppressed just by properly adjusting the relative phase of four coherent low-intensity driving fields and the atomic system becomes transparent against two-photon absorption. From a physical point of view, we explicitly explain these results in terms of quantum interference induced by two different two-photon excitation channels.
2007-05-15
Indirect stabilization of weakly coupled systems with hybrid boundary conditions
We investigate stability properties of indirectly damped systems of evolution equations in Hilbert spaces, under new compatibility assumptions. We prove polynomial decay for the energy of solutions and optimize our results by interpolation techniques, obtaining a full range of power-like decay rates. In particular, we give explicit estimates with respect to the initial data. We discuss several applications to hyperbolic systems with {\\em hybrid} boundary conditions, including the coupling of two wave equations subject to Dirichlet and Robin type boundary conditions, respectively.
2011-01-01
Higher spin symmetries and w_#infinity# algebra in the conformal affine Toda model
International Nuclear Information System (INIS)
As recently shown the conformal affine Toda models can be obtained via hamiltonian reduction from a two-loop Kac-Moody algebra. In this paper we propose a systematic procedure to analyze the higher spin symmetries of the conformal affine Toda models. The method is based on an explicit construction of infinite towers of extended conformal symmetry generators. Two fundamental building blocks of this construction are special spin-one and -two primary fields characterizing the conformal structure of these models. The connection to the algebra of area preserving diffeomorphisms on a two-manifold (w_#infinity# algebra) is established. (orig.).
1992-05-01
FORTUM Participation in BARCOM Round Robin pre-test simulation: mid-term analysis
Energy Technology Data Exchange (ETDEWEB)
In the study a preliminary mid-term analysis of the BARCOM test model is presented. The BARCOM test model is a 1:4 scale of an existing pressurized heavy water reactor (PHWR) pre-stressed concrete inner containment of 540 MW Tarapur Atomic Power Station 3 and 4 units in India. The goal of this midterm analysis is to illustrate the modeling approach and achieve a prediction of the failure mode. The analysis was carried out using ABAQUS/CAE and ABAQUS/EXPLICIT version 6.7-EF1 software
2009-07-01
Dynamics of a H(n) atom in time-dependent electric and magnetic fields
International Nuclear Information System (INIS)
The symmetry properties of the Coulomb potential allow for a dynamical spin-1/2 description of any fixed n level of a hydrogenlike atom in a time-dependent sufficiently weak electric and/or magnetic field. An explicit expression for the time dependence of the l,m amplitudes pertaining to a general n level is derived. The derivation follows on purely algebraic operations. Based on the derivation, we give analytical n-independent solutions to established and proposed schemes for driving the atom into a high angular-momentum state.
2002-05-01
Concentration of measures via size-biased couplings
British Library Electronic Table of Contents (United Kingdom)
Let Y be a nonnegative random variable with mean?? and finite positive variance ? 2, and let Y s , defined on the same space as Y, have the Y size-biased distribution, characterized by Formula Not Shown Under a variety of conditions on Y and the coupling of Y and Y s , including combinations of boundedness and monotonicity, one sided concentration of measure inequalities such as Formula Not Shown hold for some explicit A and B. The theorem is applied to the number of bulbs switched on at the terminal time in the so called lightbulb process of Rao et?al. (Sankhy? 69:137?161, 2007).
2011-01-01
Computer Algebra Solving of First Order ODEs Using Symmetry Methods
A set of Maple V R.3/4 computer algebra routines for the analytical solving of 1st. order ODEs, using Lie group symmetry methods, is presented. The set of commands includes a 1st. order ODE-solver and routines for, among other things: the explicit determination of the coefficients of the infinitesimal symmetry generator; the construction of the most general invariant 1st. order ODE under given symmetries; the determination of the canonical coordinates of the underlying invariant group; and the testing of the returned results.
1996-01-01
Bianchi type IX cosmological models with homogeneous spinor fields
Energy Technology Data Exchange (ETDEWEB)
The diagonal and symmetric Bianchi type IX models are coupled to a homogeneous spinor field. An action for the combined fields is constructed, where the orthonormal basis used is given explicitly in terms of the metric. This allows one to vary the action with respect to the metric and the spinor fields only. Next, a Hamiltonian formulation is given, and a qualitative solution for the problem is presented. We also show that the k = +1 FRW (Friedmann--Robertson--Walker) model is not compatible with a homogeneous spinor field, while the more complicated models are.
1981-03-01
Bianchi type IX cosmological models with homogeneous spinor fields
International Nuclear Information System (INIS)
The diagonal and symmetric Bianchi type IX models are coupled to a homogeneous spinor field. An action for the combined fields is constructed, where the orthonormal basis used is given explicitly in terms of the metric. This allows one to vary the action with respect to the metric and the spinor fields only. Next, a Hamiltonian formulation is given, and a qualitative solution for the problem is presented. We also show that the k = +1 FRW (Friedmann--Robertson--Walker) model is not compatible with a homogeneous spinor field, while the more complicated models are.
Let $\\A$ be an irreducible Coxeter arrangement and $\\bfk$ be a multiplicity of $\\A$. We study the derivation module $D(\\A, \\bfk)$. Any two-dimensional irreducible Coxeter arrangement with even number of lines is decomposed into two orbits under the action of the Coxeter group. In this paper, we will {explicitly} construct a basis for $D(\\A, \\bfk)$ assuming $\\bfk$ is constant on each orbit. Consequently we will determine the exponents of $(\\A, \\bfk)$ under this assumption. For this purpose we develop a theory of universal derivations and introduce a map to deal with our exceptional cases.
2010-01-01
Are classical tachyons slower-than-light quantum particles
Energy Technology Data Exchange (ETDEWEB)
After having studied the shape that a tachyon T (e.g., intrinsically spherical) would take up, we show in an explicit example that the characteristics of classical tachyons are similar to those of the ordinary (slower-than-light) quantum particles. In particular, a realistic tachyon is associated with a ''phase speed'' V(V/sup 2/>c/sup 2/), but with a ''group speed'' upsilon=c/sup 2//V (upsilon/sup 2/ 1983-07-02
Anisotropic mass density by two-dimensional acoustic metamaterials
Energy Technology Data Exchange (ETDEWEB)
We show that specially designed two-dimensional arrangements of full elastic cylinders embedded in a nonviscous fluid or gas define (in the homogenization limit) a new class of acoustic metamaterials characterized by a dynamical effective mass density that is anisotropic. Here, analytic expressions for the dynamical mass density and the effective sound velocity tensors are derived in the long wavelength limit. Both show an explicit dependence on the lattice filling fraction, the elastic properties of cylinders relative to the background, their positions in the unit cell, and their multiple scattering interactions. Several examples of these metamaterials are reported and discussed.
2008-02-15
A NEM diffusion code for fuel management and time average core calculation
International Nuclear Information System (INIS)
A computer code based on Nodal expansion method has been developed for solving two groups three dimensional diffusion equation. This code can be used for fuel management and time average core calculation. Explicit Xenon and fuel temperature estimation are also incorporated in this code. TAPP-4 phase-B physics experimental results were analyzed using this code and a code based on FD method. This paper gives the comparison of the observed data and the results obtained with this code and FD code. (author)
2005-11-01
Advanced distillation curve measurements for corrosive fluids: Application to two crude oils
Energy Technology Data Exchange (ETDEWEB)
We have recently introduced several important improvements in the measurement of distillation curves for complex fluids. The modifications include a composition-explicit data channel for each distillate fraction (for both qualitative and quantitative analysis) and corrosivity assessment of each distillate fraction. The composition-explicit information is achieved with a new sampling approach that allows precise qualitative as well as quantitative analyses of each fraction, on the fly. We have applied the new method to a variety of fluids, including simple n-alkanes, rocket propellant, gasoline, jet fuels, and a hydrocarbon fluid made corrosive with dissolved hydrogen sulfide. In the current contribution, we present the application of the advanced distillation curve method to two samples of crude oil. A primary motivation behind the work is to precisely measure the distillation curves of these oils using our advanced distillation apparatus; ...
2008-10-15
{delta}-excitations and the three-nucleon force
Energy Technology Data Exchange (ETDEWEB)
We study the three-nucleon force in chiral effective field theory with explicit {delta}-resonance degrees of freedom. We show that up to next-to-next-to-leading order, the only contribution to the isospin symmetric three-nucleon force involving the spin-3/2 degrees of freedom is given by the two-pion-exchange diagram with an intermediate delta, frequently called the Fujita-Miyazawa force. We also analyze the leading isospin-breaking corrections due to the delta. For that, we give the first quantitative analysis of the delta quartet mass splittings in chiral effective field theory including the leading electromagnetic corrections. The charge-symmetry breaking three-nucleon force due to an intermediate delta excitation is small, of the order of a few keV.
2008-06-15
The interaction of accretion disks with the magnetospheres of young stars can produce X-winds and funnel flows. With the assumption of axial symmetry and steady state flow, the problem can be formulated in terms of quantities that are conserved along streamlines, such as the Bernoulli integral (BI), plus a partial differential equation (PDE), called the Grad-Shafranov equation (GSE), that governs the distribution of streamlines in the meridional plane. The GSE plus BI yields a PDE of mixed type, elliptic before critical surfaces where the flow speed equals certain characteristic wave speeds are crossed and hyperbolic afterward. The computational difficulties are exacerbated by the locations of the critical surfaces not being known in advance. To overcome these obstacles, we consider a variational principle by which the GSE can be attacked by extremizing an action integral, with all other conserved quantities of the problem explicitly included as part of the overall ...
2007-01-01
British Library Electronic Table of Contents (United Kingdom)
In multi-organisational contexts, scenario building has been used to engage stakeholders in a critical discussion on issues of mutual importance, and to gain their support with regards to possible future responses. A review of existing literature suggests that much has been written regarding the process of scenario development and the benefits of the process, but the detailed analysis of scenario building outcomes, which encompass a large number of issues and their complex interconnections, has not been made explicit for studying and enhancing understanding of a complex societal problem. This paper presents a systematic method for analysing such complex outcomes in order to facilitate reflective thinking on important issues within the wider context for policy development. The method was em...
2011-01-01
Unitary constraints on Deeply Virtual Compton Scattering
Energy Technology Data Exchange (ETDEWEB)
At moderately low momentum transfer ($-t$ up to 1 GeV$^2$) the coupling to the vector meson production channels gives the dominant contribution to real Compton and deeply virtual Compton scattering (DVCS). Starting from a Regge Pole approach that successfully describes vector meson production, the singular part of the corresponding box diagrams (where the intermediate vector meson-baryon pair propagates on-shell) is evaluated without any further assumptions (unitarity). Such a treatment explains not only the unexpectedly large DVCS unpolarized cross section that has been recently measured at Jefferson Laboratory (JLab), but also all the beam spin and charge asymmetries that has been measured at JLab and Hermes, without explicit need of Generalized Parton Distributions (GPD). The issue of the relationship between the two approaches is addressed.
2007-11-01
Understanding collaborative filtering parameters for personalized recommendations in e-commerce
British Library Electronic Table of Contents (United Kingdom)
Collaborative Filtering (CF) is a popular method for personalizing product recommendations for e-Commerce and customer relationship management (CRM). CF utilizes the explicit or implicit product evaluation ratings of customers to develop personalized recommendations. However, there has been no in-depth investigation of the parameters of CF in relation to the number of ratings on the part of an individual customer and the total number of ratings for an item. We empirically investigated the relationships between these two parameters and CF performance, using two publicly available data sets, EachMovie and MovieLens. We conducted three experiments. The first two investigated the relationship between a particular customer?s number of ratings and CF recommendation performance. The third experim...
2007-01-01
We formulate a complete theory of Edge Radiation based on a novel method relying on Fourier Optics techniques. Similar types of radiation like Transition Undulator Radiation are addressed in the framework of the same formalism. Special attention is payed in discussing the validity of approximations upon which the theory is built. Our study makes consistent use of both similarity techniques and comparisons with numerical results from simulation. We discuss both near and far zone. Physical understanding of many asymptotes is discussed. Based on the solution of the field equation with a tensor Green's function technique, we also discuss an analytical model to describe the presence of a vacuum chamber. In particular, explicit calculations for a circular vacuum chamber are reported. Finally, we consider the use of Edge Radiation as a tool for electron beam diagnostics. We discuss Coherent Edge Radiation, Extraction of Edge Radiation by a mirror, and other issues ...
2008-01-01
The planetary N-body problem: symplectic?foliation, reductions and invariant tori
British Library Electronic Table of Contents (United Kingdom)
The 6n-dimensional phase space of the planetary (1+n)-body problem (after the classical reduction of the total linear momentum) is shown to be foliated by symplectic leaves of dimension (6n?2) invariant for the planetary Hamiltonian Formula Not Shown . Such foliation is described by means of a new global set of Darboux coordinates related to a symplectic (partial) reduction of rotations. On each symplectic leaf Formula Not Shown has the same form and it is shown to preserve classical symmetries. Further sets of Darboux coordinates may be introduced on the symplectic leaves so as to achieve a complete (total) reduction of rotations. Next, by explicit computations, it is shown that, in the reduced settings, certain degeneracies are removed. In particular, full torsion is checked both in the ...
2011-01-01
The number of cycles of specified normalized length in permutations
We compute the limiting distribution, as n approaches infinity, of the number of cycles of length between gamma n and delta n in a permutation of [n] chosen uniformly at random, for constants gamma, delta such that 1/(k+1) <= gamma < delta <= 1/k for some integer k. This distribution is supported on {0, 1, ... k} and has 0th, 1st, ..., kth moments equal to those of a Poisson distribution with parameter log (delta/gamma). For more general choices of gamma, delta we show that such a limiting distribution exists, which can be given explicitly in terms of certain integrals over intersections of hypercubes with half-spaces; these integrals are analytically intractable but a recurrence specifying them can be given. The results herein provide a basis of comparison for similar statistics on restricted classes of permutations.
2009-01-01
The kernel polynomial method for non-orthogonal electronic structure calculations
Energy Technology Data Exchange (ETDEWEB)
The Kernel Polynomial Method (KPM) has been successfully applied to tight-binding electronic structure calculations as an O(N) method. Here the authors extend this method to nonorthogonal basis sets with a sparse overlap matrix S and a sparse Hamiltonian H. Since the KPM method utilizes matrix vector multiplications it is necessary to apply S{sup {minus}1}H onto a vector. The multiplication of S{sup {minus}1} is performed using a preconditioned conjugate gradient method and does not involve the explicit inversion of S. Hence the method scales the same way as the original KPM method, i.e. O(N), although there is an overhead due to the additional conjugate gradient part. The authors show an application of this method to defects in a titanate/platinum interface and to a large scale electronic structure calculation of amorphous diamond.
1996-11-01
British Library Electronic Table of Contents (United Kingdom)
The analysis of complex fluids such as crude oils, fuels, vegetable oils and mixed waste streams poses significant challenges arising primarily from the multiplicity of components, the different properties of the components (polarity, polarizability, etc.) and matrix properties. We have recently introduced an analytical strategy that simplifies many of these analyses, and provides the added potential of linking compositional information with physical property information. This aspect can be used to facilitate equation of state development for the complex fluids. In addition to chemical characterization, the approach provides the ability to calculate thermodynamic properties for such complex heterogeneous streams. The technique is based on the advanced distillation curve (ADC) metrology, wh...
2010-01-01
British Library Electronic Table of Contents (United Kingdom)
Effective integration of nonproliferation management into the design process is key to the broad deployment of advanced nuclear energy systems, and is an explicit goal of the Laser Inertial Fusion Energy (LIFE) project at Lawrence Livermore National Laboratory. The nuclear explosives utility of a nuclear material to a state (proliferator) or sub-state (terrorist) is a critical factor to be assessed and is one aspect of material attractiveness. In this work, we approached nuclear explosives utility through the calculation of a ''figure of merit'' (FOM) that has recently been developed to capture the relative viability and difficulty of constructing nuclear explosives starting from various nuclear material forms and compositions. We discuss the integration of the figure of merit into an asse...
2010-01-01
The Minimal Scale Invariant Extension of the Standard Model
We perform a systematic analysis of an extension of the Standard Model that includes a complex singlet scalar field and is scale invariant at the tree level. We call such a model the Minimal Scale Invariant extension of the Standard Model (MSISM). The tree-level scale invariance of the model is explicitly broken by quantum corrections, which can trigger electroweak symmetry breaking and potentially provide a mechanism for solving the gauge hierarchy problem. Even though the scale invariant Standard Model is not a realistic scenario, the addition of a complex singlet scalar field may result in a perturbative and phenomenologically viable theory. We present a complete classification of the flat directions which may occur in the classical scalar potential of the MSISM. After calculating the one-loop effective potential of the MSISM, we investigate a number of representative scenarios and determine their scalar boson mass spectra, as well as their perturbatively ...
2010-01-01
We prove an analogue of the MacMahon Master Theorem for the right quantum superalgebras. In particular, we obtain a new and simple proof of this theorem for the right quantum algebras. In the super case the theorem is then used to construct higher order Sugawara operators for the affine Lie superalgebra \\hat gl(m|n) in an explicit form. The operators are elements of a completed universal enveloping algebra of \\hat gl(m|n) at the critical level. They occur as the coefficients in the expansion of a noncommutative Berezinian and as the traces of powers of generator matrices. The same construction yields higher Hamiltonians for the Gaudin model associated with the Lie superalgebra gl(m|n).
2009-01-01
The Exchange Value Embedded In A Transport System
This paper shows that a well designed transport system has an embedded exchange value by serving as a market for potential exchange between consumers. Under suitable conditions, one can improve the welfare of consumers in the system simply by allowing some exchange of goods between consumers during transportation without incurring additional transportation costs. We propose an explicit valuation formula to measure this exchange value for a given compatible transport system. This value is always nonnegative and bounded from above. Criteria based on transport structures, preferences and prices are provided to determine the existence of a positive exchange value. Finally, we study a new optimal transport problem with an objective taking into account of both transportation cost and exchange value.
2010-01-01
Teaching elementary thermodynamics and energy conversion: Opinions
British Library Electronic Table of Contents (United Kingdom)
This presentation deals with innovation in teaching and understanding of thermodynamic principles. Key features of the approach being advocated are: (a) postulation of the existence of entropy, (b) explicitly associating energy transfers with other transports, (c) stating the 2nd Law in terms of Gibbs' available-energy, (d) systematic use of software such as EES. The paper outlines and elaborates upon an introductory course. Major headings in the course are: basic concepts: properties, additive properties and balances, primitive properties, energy, 1st Law. entropy, elementary academic applications of balances, available-energy, second law, exergy, thermostatic property relations, EES. Applications to processes, fluid flow, Heat transfer, thermochemical. Applications to devices, single-pro...
2010-01-01
Study of azeotropic mixtures with the advanced distillation curve approach
British Library Electronic Table of Contents (United Kingdom)
Classical methods for the study of complex fluid phase behavior include static and dynamic equilibrium cells that usually require vapor and liquid recirculation. These are sophisticated, costly apparatus that require highly trained operators, usually months of labor-intensive work per mixture, and the data analysis is also rather complex. Simpler approaches to the fundamental study of azeotropes are highly desirable, even if they provide only selected cuts through the phase diagram. Recently, we introduced an advanced distillation curve measurement method featuring: (1) a composition explicit data channel for each distillate fraction (for both qualitative and quantitative analysis), (2) temperature measurements that are true thermodynamic state points that can be modeled with an equation o...
2009-01-01
Energy Technology Data Exchange (ETDEWEB)
The linearized collective Schroedinger equation for nuclear quadrupole surface vibrations incorporates a new spin degree of freeedom with a spin value of 3/2. We use this equation to describe the low energy spectrum of certain even-odd Ir nuclei which have a spin 3/2 in their ground state. For that purpose we explicitly introduce collective spin-dependent potentials which simulate the interaction of the valence nucleon with the core. The linearized Schroedinger equation is transformed into an effective Schroedinger equation with collective spin-dependent potentials. Already collective spin-orbit couplings of SO(3) and SO(5) type are sufficient to reproduce the lowest excited states of even-odd Ir nuclei. (orig.).
1990-06-01
International Nuclear Information System (INIS)
The linearized collective Schroedinger equation for nuclear quadrupole surface vibrations incorporates a new spin degree of freeedom with a spin value of 3/2. We use this equation to describe the low energy spectrum of certain even-odd Ir nuclei which have a spin 3/2 in their ground state. For that purpose we explicitly introduce collective spin-dependent potentials which simulate the interaction of the valence nucleon with the core. The linearized Schroedinger equation is transformed into an effective Schroedinger equation with collective spin-dependent potentials. Already collective spin-orbit couplings of SO(3) and SO(5) type are sufficient to reproduce the lowest excited states of even-odd Ir nuclei. (orig.).
Social-ecological science in the humane metropolis
British Library Electronic Table of Contents (United Kingdom)
The Humane metropolis is a rubric to summarize and promote environmental and social quality in contemporary urban mosaics. Because cities, suburbs, and exurbs, as spatially extensive and connected socio-ecological systems, exhibit many negative features, the humane metropolis identifies a strategy to combat the ills and instill more positive and sustainable features and processes in urban systems. Because the humane metropolis as a program has arisen primarily from social motivations, there is the opportunity to articulate more explicitly the role that science can play in addressing the humane metropolis program and evaluating its success. A humane metropolis can be summarized as one that 1) protects and restores ecological services in cities and suburbs, 2) promotes physical and mental he...
2011-01-01
Scattering of discrete states in two-dimensional open string field theory
Energy Technology Data Exchange (ETDEWEB)
This is the second in a series of papers devoted to open string field theory in two dimensions. In this paper we aim to clarify the origin and the role of discrete physical states in the theory. To this end, we study interactions of discrete states and generic tachyons. In particular, we discuss at length four point amplitudes. We show that the behavior of the correlation functions is governed by the number of generic tachyons involved and values of the kinematic invariants [ital s], [ital t], and [ital u]. The divergence of certain classes of correlators is shown to be the consequence of the fact that certain kinematic invariants are nonpositive integers in that case. Explicit examples are included. We check our results by the standard conformal technique.
1993-12-15
British Library Electronic Table of Contents (United Kingdom)
Abstract: The recent National Academy of Sciences report on forensic sciences states that the study of fire patterns and debris in arson fires is in need of additional work and eventual standardization. We discuss a recently introduced method that can provide predicted evaporation patterns for ignitable liquids as a function of temperature. The method is a complex fluid analysis protocol, the advanced distillation curve approach, featuring a composition explicit data channel for each distillate fraction (for qualitative, quantitative, and trace analysis), low uncertainty temperature measurements that are thermodynamic state points that can be modeled with an equation of state, consistency with a century of historical data, and an assessment of the energy content of each distillate fraction...
2011-01-01
Pinned fluxons in a Josephson junction with a finite-length inhomogeneity
We consider a Josephson junction system installed with a finite length inhomogeneity, either of microresistor or of microresonator type. The system can be modelled by a sine-Gordon equation with a piecewise-constant function to represent the varying Josephson tunneling critical current. The existence of pinned fluxons depends on the length of the inhomogeneity, the variation in the Josephson tunneling critical current and the applied bias current. We establish that a system may either not be able to sustain a pinned fluxon, or - for instance by varying the length of the inhomogeneity - may exhibit various different types of pinned fluxons. Our stability analysis shows that changes of stability can only occur at critical points of the length of the inhomogeneity as a function of the (Hamiltonian) energy density inside the inhomogeneity - a relation we determine explicitly. In combination with continuation arguments and Sturm-Liouville theory, we determine the ...
2011-01-01
A method for systematically including topological degrees of freedom in perturbation theory is developed. This is not bound by the restrictions of semi-classical techniques. The Yang-Mills theory in three Euclidean dimensions is considered here. A well-defined separation of the topological and the ``spin wave'' degrees of freedom is obtained, motivated by a singular gauge. This has ``photons'' distorting the spherically symmetric magnetic fields of Dirac monopoles, and massless charged vector bosons ``W'' scattering off the latter. It is explicitly shown that the Dirac string does not contribute. The mode of the charged vector bosons with total angular momentum J=0 provides precisely the core to give a finite energy to the monopole. The radial equation for W is remarkably simplified and only two polarization states survive exactly for the anomalous magnetic moment required by the Yang-Mills interaction.
2003-01-01
Partition Functions of Normal Factor Graphs
One of the most common types of functions in mathematics, physics, and engineering is a sum of products, sometimes called a partition function. After "normalization," a sum of products has a natural graphical representation, called a normal factor graph (NFG), in which vertices represent factors, edges represent internal variables, and half-edges represent the external variables of the partition function. In physics, so-called trace diagrams share similar features. We believe that the conceptual framework of representing sums of products as partition functions of NFGs is an important and intuitive paradigm that, surprisingly, does not seem to have been introduced explicitly in the previous factor graph literature. Of particular interest are NFG modifications that leave the partition function invariant. A simple subclass of such NFG modifications offers a unifying view of the Fourier transform, tree-based reparameterization, loop calculus, and the Legendre ...
2011-01-01
Noncentral extension of the $AdS_5 x S^5$ superalgebra supermultiplet of brane charges
We propose an extension of the su(2,2|4) superalgebra to incorporate the F1/D1 string charges in type IIB string theory on the AdS_5 X S^5 background, or the electro-magnetic charges in the dual super Yang-Mills theory. With the charges introduced, the superalgebra inevitably undergoes a noncentral extension, as noted recently in [1]. After developing a group theoretical method of obtaining the noncentral extension, we show that the charges form a certain nonunitary representation of the original unextended superalgebra, subject to some constraints. We solve the constraints completely and show that, apart from the su(2,2|4) generators, there exist 899 complex brane charges in the extended algebra. Explicitly we present all the super-commutators among them.
2004-01-01
Non-Gaussian gravitational clustering field statistics
In this work we investigate the multivariate statistical description of the matter distribution in the nonlinear regime. We introduce the multivariate Edgeworth expansion of the lognormal distribution to model the cosmological matter field. Such a technique could be useful to generate and reconstruct three-dimensional nonlinear cosmological density fields with the information of higher order correlation functions. We explicitly calculate the expansion up to third order in perturbation theory making use of the multivariate Hermite polynomials up to sixth order. The probability distribution function for the matter field includes at this level the two-point, the three-point and the four-point correlation functions. We use the hierarchical model to formulate the higher order correlation functions based on combinations of the two-point correlation function. This permits us to find compact expressions for the skewness and kurtosis terms of the expanded lognormal field ...
2010-01-01
Neutrinos and long-range weak forces in cosmology
International Nuclear Information System (INIS)
The Lorentz and coordinate covariant calculus of spinors in Riemannian spacetime, which is the mathematical model for the description of the quantum mechanics of elementary particles with spin interacting with the classical gravitation field, is explored. The Dirac equation describing the interaction of neutrinos with the gravitational fields of the Robertson-Walker cosmological world models is separated, and the spectrum of eigenfunctions and eigenvalues for particular choices of the set of quantum numbers is given explicitly for the k = 0 and k = +1 models, although only the radial equations determining the final quantum number are given for the k = -1 model. The mathematical theory of the motion of a perfect fluid whose elements interact via long-range neutrino-exchange forces, as well as gravitationally, is developed. The formalism for calculating, by calculating the Bogoliubov transformation of the Fock space operators that instantaneously diagonalize their ...
Modeling Credit Risk with Partial Information
This paper provides an alternative approach to Duffie and Lando [Econometrica 69 (2001) 633-664] for obtaining a reduced form credit risk model from a structural model. Duffie and Lando obtain a reduced form model by constructing an economy where the market sees the manager's information set plus noise. The noise makes default a surprise to the market. In contrast, we obtain a reduced form model by constructing an economy where the market sees a reduction of the manager's information set. The reduced information makes default a surprise to the market. We provide an explicit formula for the default intensity based on an Azema martingale, and we use excursion theory of Brownian motions to price risky debt.
2008-01-01
Modal Meinongianism and fiction: the best of three worlds
British Library Electronic Table of Contents (United Kingdom)
We outline a neo-Meinongian framework labeled as Modal Meinongian Metaphysics (MMM) to account for the ontology and semantics of fictional discourse. Several competing accounts of fictional objects are originated by the fact that our talking of them mirrors incoherent intuitions: mainstream theories of fiction privilege some such intuitions, but are forced to account for others via complicated paraphrases of the relevant sentences. An ideal theory should resort to as few paraphrases as possible. In Sect.?1, we make this explicit via two methodological principles, called the Minimal Revision and the Acceptability Constraint. In Sect.?2, we introduce the standard distinction between internal and external fictional discourse. In Sects.?3?5, we discuss the approaches of (traditional) Meinongia...
2011-01-01
Energy Technology Data Exchange (ETDEWEB)
We consider the dynamics towards the initial singularity of Bianchi type IX vacuum and orthogonal perfect fluid models with a linear equation of state. Surprisingly few facts are known about the 'Mixmaster' dynamics of these models, while at the same time most of the commonly held beliefs are rather vague. In this paper, we use Mixmaster facts as a base to build an infrastructure that makes it possible to sharpen the main Mixmaster beliefs. We formulate explicit conjectures concerning (i) the past asymptotic states of type IX solutions and (ii) the relevance of the Mixmaster/Kasner map for generic past asymptotic dynamics. The evidence for the conjectures is based on a study of the stochastic properties of this map in conjunction with dynamical systems techniques. We use a dynamical systems formulation, since this approach has so far been the only successful path to obtain theorems, but we also make comparisons with the 'metric' ...
2009-04-07
Measuring perceptual centers using the phase correction response
British Library Electronic Table of Contents (United Kingdom)
The perceptual center (P-center) is fundamental to the timing of heterogeneous event sequences, including music and speech. Unfortunately, there is currently no comprehensive and reliable model of P-centers in acoustic events, so P-centers must instead be measured empirically. This study reviews existing measurement methods and evaluates two methods in detail?the rhythm adjustment method and a new method based on the phase correction response (PCR) in a synchronous tapping task. The two methods yielded consistent P-center estimates and showed no evidence of P-center context dependence. The PCR method appears promising because it is accurate and efficient and does not require explicit perceptual judgments. As a secondary result, the magnitude of the PCR is shown to vary systematically with ...
2011-01-01
Market-based pricing and demand-side participation in the electricity pool of England and Wales
Energy Technology Data Exchange (ETDEWEB)
This paper focuses upon the workings of the electricity pool of England and Wales. Created as a generators` pool, with the twin objectives of preserving order-of-merit dispatch and setting marginal-cost pricing signals, it is now evident that the system may not provide an efficient market-clearing mechanism from the consumers` perspective. Thus, over the five years since its inception, there have been repeated calls to incorporate more explicitly demand-side pricing signals into the process. A number of official reviews and various consultants have looked at the issue and there have been trial periods of active demand-side bidding. This paper reviews that experience to date.
1996-03-01
Magnetotelluric inversion via reverse time migration algorithm of seismic data
We propose a new algorithm for two-dimensional magnetotelluric (MT) inversion. Our algorithm is an MT inversion based on the steepest descent method, borrowed from the backpropagation technique of seismic inversion or reverse time migration, introduced in the middle 1980s by Lailly and Tarantola. The steepest descent direction can be calculated efficiently by using the symmetry of numerical Green's function derived from a mixed finite element method proposed by Nedelec for Maxwell's equation, without calculating the Jacobian matrix explicitly. We construct three different objective functions by taking the logarithm of the complex apparent resistivity as introduced in the recent waveform inversion algorithm by Shin and Min. These objective functions can be naturally separated into amplitude inversion, phase inversion and simultaneous inversion. We demonstrate our algorithm by showing three inversion results for synthetic data.
2007-07-01
Land use change in a biofuels hotspot: The case of Iowa, USA
British Library Electronic Table of Contents (United Kingdom)
This study looks at the land use impact of the biofuels expansion on both the intensive and extensive margin, and its environmental consequences. We link economic, geographical and environmental models by using spatially explicit common units of analysis and use remote sensing crop cover maps and digitized soils data as inputs. Land use changes are predicted via economic analysis of crop rotation choice and tillage under alternative crop prices, and the Environmental Policy Integrated Climate (EPIC) model is used to predict corresponding environmental impacts. The study focuses on Iowa, which is the leading biofuels hotspot in the U.S. due to intensive corn production and the high concentration of ethanol plants that comprise 28% of total U.S. production. We consider the impact of the biof...
2011-01-01
Instanton fermionic zero mode at finite temperature and chemical potential
In QCD the spontaneous breaking of chiral symmetry and the U(1) axial anomaly can be understood considering instantons as the gauge configurations mediating quark-quark interaction. The existence of an exact zero mode solution of the Dirac equation in the field of a single instanton is the fundamental ingredient of this analysis. Explicit expressions for psi_0 are available for T different from 0 and mu=0, and mu different from 0 and T=0. In this paper we derive the solution for the most general case T different from 0 and mu different from 0. This new result opens the possibility of investigating the QCD dynamics associated with instantons in the full phase diagram. As a first step in this direction we will study the dependence of the instanton density from the thermodynamic coordinates.
2011-01-01
Generation of coherent states of photon-added type via pathway of eigenfunctions
International Nuclear Information System (INIS)
We obtain and investigate the regular eigenfunctions of simple differential operators xr dr+1/dxr+1, r = 1, 2, ..., with the eigenvalues equal to 1. With the help of these eigenfunctions, we construct a non-unitary analogue of a boson displacement operator which will be acting on the vacuum. In this way, we generate collective quantum states of the Fock space which are normalized and equipped with the resolution of unity with the positive weight functions that we obtain explicitly. These states are thus coherent states in the sense of Klauder. They span the truncated Fock space without first r lowest-lying basis states: |0), |1), ..., |r - 1). These states are squeezed, sub-Poissonian in nature and reminiscent of photon-added states in Agarwal and Tara (1991 Phys. Rev. A 43 492).
2010-09-17
Environmental Radiation Effects: A Need to Question Old Paradigms
Energy Technology Data Exchange (ETDEWEB)
A historical perspective is given of the current paradigm that does not explicitly protect the environment from radiation, but instead, relies on the concept that if dose limits are set to protect humans then the environment is automatically protected as well. We summarize recent international questioning of this paradigm and briefly present three different frameworks for protecting biota that are being considered by the U.S. DOE, the Canadian government and the International Commission on Radiological Protection. We emphasize that an enhanced collaboration is required between what has traditionally been separated disciplines of radiation biology and radiation ecology if we are going to properly address the current environmental radiation problems. We then summarize results generated from an EMSP grant that allowed us to develop a Low Dose Irradiation Facility that specifically addresses effects of low-level, chronic irradiation on multiple levels of biological ...
2003-03-27
Entropy and Hausdorff Dimension in Random Growing Trees
We investigate the limiting behavior of random tree growth in preferential attachment models. The tree stems from a root, and we add vertices to the system one-by-one at random, according to a rule which depends on the degree distribution of the already existing tree. The so-called weight function, in terms of which the rule of attachment is formulated, is such that each vertex in the tree can have at most K children. We define the concept of a certain random measure mu on the leaves of the limiting tree, which captures a global property of the tree growth in a natural way. We prove that the Hausdorff and the packing dimension of this limiting measure is equal and constant with probability one. Moreover, the local dimension of mu equals the Hausdorff dimension at mu-almost every point. We give an explicit formula for the dimension, given the rule of attachment.
2010-01-01
Energy and reserve dispatch in a multi-zone electricity market
Energy Technology Data Exchange (ETDEWEB)
Competitive bidding for ancillary services (A/S), in addition to bidding for energy, is increasingly recognized as an important part of electricity markets. In this paper, an LP-based approach is presented for solving the multi-product (i.e. energy and A/S products), multi-zone/area physical market dispatch problem. The approach can explicitly represent various constraints, such as reserve requirement constraints, network security limits, and sharing of constrained resource capacity (e.g. maximum unit generation) between energy and reserve. The joint dispatch method, which dispatches energy and reserve market concurrently, is extended to allow inter-zonal reserve trading in multi-zone competitive pool. Numerical examples are included to demonstrate the salient characteristics of the joint dispatch solution.
1999-08-01
Effect of aging on learning new names and descriptions for objects
British Library Electronic Table of Contents (United Kingdom)
The present study consisted of two experiments to compare new word learning in healthy young (N = 11) and older (N = 17) adults within an explicit learning paradigm. Experiment 1 investigated the new name learning for familiar objects, while Experiment 2 investigated learning names and descriptions for unfamiliar objects. Participants attended five learning sessions over 5 consecutive days, during which they viewed objects with novel names with/without descriptions. The older adults were as accurate as the young adults when recalling and recognizing new names during the learning sessions. With respect to response times, the older adults were as rapid as the young adults at recognizing the new names for the familiar objects, but were slower during the follow-up sessions. The older adults we...
2011-01-01
Ecological risk assessment of water environment for Luanhe River Basin based on relative risk model
British Library Electronic Table of Contents (United Kingdom)
The relative risk model (RRM) was applied in regional ecological risk assessments successfully. In this study, the RRM was developed through increasing the data of risk source and introducing the source?stressor?habitat exposure filter (SSH), the endpoint?habitat exposure filter (EH) and the stressor?endpoint effect filter (SE) to reflect the meaning of exposure and effect more explicit. Water environment which include water quality, water quantity and aquatic ecosystems was selected as the ecological risk assessment endpoints. The Luanhe River Basin located in the North China was selected as model case. The results showed that there were three low risk regions, one medium risk region and two high risk regions in the Luanhe River Basin. The results also indicated habitat destruction was th...
2010-01-01
Dissipative dynamics of circuit-QED in the mesoscopic regime
We investigate the behavior of a circuit QED device when the resonator is initially populated with a mesoscopic coherent field. The strong coupling between the cavity and the qubit produces an entangled state involving mesoscopic quasi-pointer states with respect to cavity dissipation. The overlap of the associated field components results in collapse and revivals for the Rabi oscillation. Although qubit relaxation and dephasing do not preserve these states, a simple analytical description of the dissipative dynamics of the circuit QED device including cavity relaxation as well as qubit dissipation is obtained from the Monte-Carlo approach. Explicit predictions for the spontaneous and induced Rabi oscillation signals are derived and sucessfully compared with exact calculations. We show that these interesting effects could be observed with a 10 photon field in forthcoming circuit QED experiments.
2007-01-01
Dirichlet mean identities and laws of a class of subordinators
An interesting line of research is the investigation of the laws of random variables known as Dirichlet means. However, there is not much information on interrelationships between different Dirichlet means. Here, we introduce two distributional operations, one of which consists of multiplying a mean functional by an independent beta random variable, the other being an operation involving an exponential change of measure. These operations identify relationships between different means and their densities. This allows one to use the often considerable analytic work on obtaining results for one Dirichlet mean to obtain results for an entire family of otherwise seemingly unrelated Dirichlet means. Additionally, it allows one to obtain explicit densities for the related class of random variables that have generalized gamma convolution distributions and the finite-dimensional distribution of their associated L\\'{e}vy processes. The importance of this latter statement is ...
2010-01-01
International Nuclear Information System (INIS)
Dimensioning of actinides waste packages for long duration storage has to take into account helium production from natural decay and release rates from the material. For the latter, we propose here an improved method for the determination of the helium diffusion coefficient in britholite, to be used for minor actinides storage. This work is based on results we previously published using the classical three steps method: "3He implantation on a Van de Graaff facility, "3He profile determination analysing the protons resulting from the "3He(d,p)"4He reaction in a nuclear microprobe, evolution of the helium profile during annealings. Taking explicitly into account the incident deuterons energy stragglings allows us to show that the implanted helium profiles are bimodal, each component leading to a different helium diffusion coefficient.
2005-02-01
Creating Virtual Collections in Digital Libraries: Benefits and Implementation Issues.
Explicitly defining virtual collections in a digital library can benefit both users and contributors. This paper outlines the costs and benefits for defining virtual collections and describes an implementation of collection-level metadata in two digital libraries. Our work suggests that component repositories of NSDL would profit by describing their collection, and sub-collections, using a schema such as the one we implemented. This would enable NSDL to provide a rich, cross-linked hierarchy of virtual collections. By sharing collection-level, rather than item-level, metadata, collections would provide detailed resource information, presumably keeping most of their viewership, while NSDL would provide a broad description of resources.
2003-03-18
Computational experience with an explicit difference scheme for a one phase Stefan problem
The results of a computational scheme of Rose (a method for calculating solutions of parabolic equations with free boundary, Math. Comp., 14: 249-256(1960)) are compared with the known analytic solution of a single-phase Stefan problem with constant boundary temperature. Rose's scheme was implemented in FORTRAN on a DEC PDP-10. Computations were performed for Stefan numbers St = 0.01, 0.1, 1, 10, and 100. Three relative mesh sizes (..gamma.. = DT/(DX)/sup 2/ = 1/2, 1/6, and 1/25) were used. For each St, the relative mesh size had very little effect on accuracy. However, taking ..gamma.. = 1/25 resulted in dramatically increased computing time. The relative errors of temperature were large near the phase change front. 11 figures, 5 tables.
1980-07-01
Color Octet Contribution in Exclusive P-Wave Charmonium Decay
Recent advances in our understanding of the higher-wave quarkonia have generated much interests in quarkonium physics. However most are devoted to inclusive decays and productions. Experimental data of several two-body exclusive decay channels of P-wave charmonia such as \\pi \\pi and p \\bar p are available and some have recently been re-measured by the BES collaboration. It is not clear from the outset that color octet is needed for these exclusive channels. Indeed only color singlet has been used in the past and reasonable agreement with data was found. Contrary to these old results, we provide theoretical arguments for the inclusion of color octet and perform explicit calculations to back this up.
2001-01-01
Chern classes of graph hypersurfaces and deletion-contraction
We study the behavior of the Chern classes of graph hypersurfaces under the operation of deletion-contraction of an edge of the corresponding graph. We obtain an explicit formula when the edge satisfies two technical conditions, and prove that both these conditions hold when the edge is multiple in the graph. This leads to recursions for the Chern classes of graph hypersurfaces for graphs obtained by adding parallel edges to a given (regular) edge. Analogous results for the case of Grothendieck classes of graph hypersurfaces were obtained in previous work. Both Grothendieck classes and Chern classes were used to define `algebro-geometric' Feynman rules. The results in this paper provide further evidence that the polynomial Feynman rule defined in terms of the Chern-Schwartz-MacPherson class of a graph hypersurface reflects closely the combinatorics of the corresponding graph. The key to the proof of the main result is a more general formula for the ...
2011-01-01
Radiation protection of the personnel who will perform interventions in the LHC Beam Cleaning Insertions is mandatory and includes the design of equipment and the establishment of work procedures. Residual dose rates due to activated equipment are expected to reach significant values such that any maintenance has to be planned and optimized in advance. Three-dimensional maps of dose equivalent rates at different cooling times after operation of the LHC have been calculated with FLUKA. The simulations are based on an explicit calculation of induced radioactivity and of the transport of the radiation from the radioactive decay. The paper summarizes the results for the Beam Cleaning Insertions and discusses the estimation of individual and collective doses received by personnel during critical interventions, such as the exchange of a collimator or the installation of Phase 2. The given examples outline the potential and the need to optimize, in an iterative way, the ...
2005-01-01
Bound state QED effects from the Schroedinger equation
International Nuclear Information System (INIS)
We present a new relativistic bound-state formalism for two interacting Fermi-Dirac particles. The kernel of the integral equation for the bound-state system is generated by summing Feynman scattering amplitudes and multiplying by a bound-state amplitude. The method is illustrated through calculations of the hyperfine and fine splittings of positronium up to order #alpha#"5. Our calculations of the one-loop contributions are carried out in the explicitly covariant Feynman gauge. We also present new results for the hyperfine and fine splittings in positronium to order #alpha#"5 for arbitrary principal quantum number n, which are easily obtained owing to the virtue of conceptual and calculational simplicity of our formalism. In addition, we present the one-loop renormalization scheme in our formalism. (author).
Asymmetric Totally-corrective Boosting for Real-time Object Detection
Real-time object detection is one of the core problems in computer vision. The cascade boosting framework proposed by Viola and Jones has become the standard for this problem. In this framework, the learning goal for each node is asymmetric, which is required to achieve a high detection rate and a moderate false positive rate. We develop new boosting algorithms to address this asymmetric learning problem. We show that our methods explicitly optimize asymmetric loss objectives in a totally corrective fashion. The methods are totally corrective in the sense that the coefficients of all selected weak classifiers are updated at each iteration. In contract, conventional boosting like AdaBoost is stage-wise in that only the current weak classifier's coefficient is updated. At the heart of the totally corrective boosting is the column generation technique. Experiments on face detection show that our methods outperform the state-of-the-art asymmetric boosting methods.
2010-01-01
Assurance and assessment techniques for nuclear weapon related software
Energy Technology Data Exchange (ETDEWEB)
Sandia National Laboratories has the qualification evaluation responsibility for the design of certain components intended for use in nuclear weapons. Specific techniques in assurance and assessment have been developed to provide the quality evidence that the software has been properly qualified for use. Qualification Evaluation is a process for assessing the suitability of either a process used to develop or manufacture the product, or the product itself The qualification process uses a team approach to evaluating a product or process, chaired by a Quality Assurance professional, with other members representing the design organization, the systems organization, and the production agency. Suitable for use implies that adequate and appropriate definition and documentation has been produced and formally released, adequate verification and validation activities have taken place to ensure proper operation, and the software product meets all requirements, explicitly or ...
1993-12-31
An improved noise-robust voice activity detector based on hidden semi-Markov models
British Library Electronic Table of Contents (United Kingdom)
To improve the performance of voice activity detector (VAD) in noisy environments, this paper concentrates on three critical aspects related to noise robustness including speech features, feature distributions and temporal dependence. Based on the statistic on TIMIT and NOIZEUS, Mel-frequency cepstrum coefficients (MFCCs) are selected as speech features, Gaussian Mixture distributions (GMD) are applied to associate the observations in MFCC domain with both speech and non-speech states, and Weibull and Gamma distributions are used to explicitly model noise and speech durations, respectively. To integrate these aspects into VAD, the hidden semi-Markov model (HSMM) as a generalized hidden Markov model (HMM) is introduced first. Then the VAD decision is made according to the likelihood ratio t...
2011-01-01
An analytical framework for capacity development in EIA ? The case of Yemen
British Library Electronic Table of Contents (United Kingdom)
Most countries worldwide nowadays apply Environmental Assessment (EA) as an ex ante tool to evaluate environmental impacts of policies, plans, programmes, and projects. However, the application and performance of EA differ significantly. Scientific analysis of how EA performs mainly focuses on two levels: the micro (or project) level and the macro (or system) level. Macro level analysis usually focuses on institutions for EA and the organisation of stakeholder interaction in EA. This article proposes a more comprehensive framework for analysing EA systems that combines other approaches with a capacity approach and an explicit consideration of the context in which EA systems are developed and performed. In order to illustrate the value of our framework, we apply it to the Republic of Yemen,...
2010-01-01
An Intrinsic Approach to Forces in Magnetoelectric Media
This paper offers a conceptually straightforward method for the calculation of stresses in polarisable media based on the notion of a drive form and its property of being closed in spacetimes with symmetry. After an outline of the notation required to exploit the powerful exterior calculus of differential forms, a discussion of the relation between Killing isometries and conservation laws for smooth and distributional drive forms is given. Instantaneous forces on isolated spacetime domains and regions with interfaces are defined, based on manifestly covariant equations of motion. The remaining sections apply these notions to media that sustain electromagnetic stresses, with emphasis on homogeneous magnetoelectric material. An explicit calculation of the average pressure exerted by a monochromatic wave normally incident on a homogeneous, magnetoelectric slab in vacuo is presented and the concluding section summarizes how this pressure depends on the parameters in ...
2008-01-01
Airbags to Martian Landers: Analyses at Sandia National Laboratories
Energy Technology Data Exchange (ETDEWEB)
A new direction for the national laboratories is to assist US business with research and development, primarily through cooperative research and development agreements (CRADAs). Technology transfer to the private sector has been very successful as over 200 CRADAs are in place at Sandia. Because of these cooperative efforts, technology has evolved into some new areas not commonly associated with the former mission of the national laboratories. An example of this is the analysis of fabric structures. Explicit analyses and expertise in constructing parachutes led to the development of a next generation automobile airbag; which led to the construction, testing, and analysis of the Jet Propulsion Laboratory Mars Environmental Survey Lander; and finally led to the development of CAD based custom garment designs using 3D scanned images of the human body. The structural analysis of these fabric structures is described as well as a more traditional example Sandia with the ...
1994-03-01
British Library Electronic Table of Contents (United Kingdom)
We have recently introduced several important improvements in the measurement of distillation curves for complex fluids. This new method is a significant improvement over current approaches, with features including a composition-explicit data channel for each distillate fraction (for both qualitative and quantitative analysis), trace chemical analysis of each distillate fraction, and a corrosivity assessment of each distillate fraction. In the present contribution, we present the application of the advanced distillation curve method to a bio-derived crude oil prepared from swine manure. One motivation behind the work was to precisely measure the distillation curve of the oil derived from swine manure with our advanced distillation apparatus and compare the bio-derived oil to traditional pe...
2008-01-01
Advanced Distillation Curve Measurement with a Model Predictive Temperature Controller
British Library Electronic Table of Contents (United Kingdom)
In previous work, several significant improvements in the measurement of distillation curves for complex fluids were introduced. The modifications to the classical measurement provide for (1) temperature and volume measurement(s) of low uncertainty, and most important, (2) a composition-explicit data channel in addition to the usual temperature?volume relationship. This latter modification is achieved with a new sampling approach that allows precise qualitative as well as quantitative analyses of each fraction, on the fly. In the new approach, the distillation temperature is measured in two locations. The temperature is measured in the usual location, at the bottom of the take-off in the distillation head, but it is also measured directly in the fluid. We have further modified our developm...
2006-01-01
Adiabatic matter effect with three generation neutrinos and the solar neutrino problem
We find an exact analytic solution for the time evolution of a three Dirac neutrino system adiabatically oscillating in matter, constructing explicitly the relevant 3\\times 3 mixing matrix in matter. Using this result we investigate the solar neutrino data in a scenario where the neutrino masses are such that m_1\\alt m_2\\ll m_3, taking into account several phenomenological constraints on neutrino mixing angles and masses. A solution of the solar neutrino problem for large values of the parameter \\delta m^2=m_2^2-m_1^2 which are not usually associated with a resonance is found. This is an essentially three-generation effect.
1994-01-01
A unified framework for biological evolution and stochastic quantization
We investigate the profound relation between the equations of biological evolution and quantum mechanics by writing a biologically inspired equation for the stochastic dynamics of an ensemble of particles. Interesting behavior is observed which is related to a new type of stochastic quantization. We find that the probability distribution of the ensemble of particles can be decomposed into eigenfunctions associated to a discrete spectrum of eigenvalues. In absence of interactions between the particles, the out-of-equilibrium dynamics asymptotically relaxes towards the fundamental state. This phenomenon can be related with the Fisher theorem in biology. On the contrary, in presence of scattering processes the evolution reaches a steady state in which the distribution of the ensemble of particles is characterized by a Bose-Einstein statistics. In order to show a concrete example of this stochastic quantization we have solved explicitly the case in which the potential ...
2010-01-01
International Nuclear Information System (INIS)
This study is concerned to the perforation phenomena of the oblique dual-plate by projectile. Experiment and simulation related to that was carried out. the variables considered in this phenomena include the electrolytic zinc coated steel sheet and carbon steel rod. In the former, the confirmation and projectile velocity possible phenomena of real phenomena is done, the latter, the effect of parameter such as time-step and grid space length is analyzed by using the three-dimensional Lagrangian explicit time-integration finite element code, HEMP. This code use the eight node hexahedral elements and in this study, Von-Mises Criteria is used as the strength model, Mie-Gruneisen is as the Equation of State. The simulation was performed by contrast with the experiment. Through the calibration of the parameter of Lagrangian code, reasonable result was approached.
2004-11-03
Energy Technology Data Exchange (ETDEWEB)
We discuss the renormalisation properties of the complete set of {delta}B=2 four-quark operators with the heavy quark treated in the static approximation. We elucidate the role of heavy quark symmetry and other symmetry transformations in constraining their mixing under renormalisation. By employing the Schroedinger functional, a set of non-perturbative renormalisation conditions can be defined in terms of suitable correlation functions. As a first step in a fully non-perturbative determination of the scale-dependent renormalisation factors, we evaluate these conditions in lattice perturbation theory at one loop. Thereby we verify the expected mixing patterns and determine the anomalous dimensions of the operators at NLO in the Schroedinger functional scheme. Finally, by employing twisted-mass QCD it is shown how finite subtractions arising from explicit chiral symmetry breaking can be avoided completely. (Orig.)
2006-04-15
A robust H{sub -} power system stabilizer design using reduced-order models
Energy Technology Data Exchange (ETDEWEB)
This paper deals with a robust H{sub -} power system stabilizer (HPSS) design using reduced-order models to improve the damping oscillation in power systems. The stabilizer is dynamic, low order and robust. In order to obtain a reduced-order controller, the method of balanced truncation is used. Sufficient conditions in the form of two algebraic Riccati equations (AREs) and an upper bound explicitly characterize an H{sub -} controller of lower dimensions. Furthermore, the bilinear transformation has been used to the design to prevent the pole-zero cancellation of the poorly damped poles and to improve the control system performance. The proposed technique is illustrated with applications to the design of stabilizer for a multi-machine power system. Simulation results under various operation conditions are given which show that the proposed HPSS damps the low-frequency oscillation in an efficient manner. (author)
2006-01-15
A Bayesian Bioeconometric Model of Invasive Species Control: The Case of the Hemlock Woolly Adelgid
British Library Electronic Table of Contents (United Kingdom)
In this article we evaluate a US Forest Service plan to mitigate damages from an invasive insect on public, forested land. We develop a dynamic model of infestation and control to explicitly account for biological interactions, baseline conditions, and uncertainty, thus creating a more complete picture of policy impacts than a static cost benefit analysis could provide. We combine the results of the dynamic model with an empirical study of nonmarket forest benefits to create a bioeconomic model of ecosystem management. Estimating the empirical model in a Bayesian framework allows us to treat the economic coefficients of the dynamic model as random variables. We specify distributions for the biological parameters and examine the effects of both biological and economic uncertainty on the pre...
2011-01-01
Systematics of high temperature perturbation theory: The two-loop electron self-energy in QED
International Nuclear Information System (INIS)
In order to investigate the systematics of the loop expansion in high temperature gauge theories beyond the leading order hard thermal loop (HTL) approximation, we calculate the two-loop electron proper self-energy #SIGMA# in high temperature QED. The two-loop bubble diagram of #SIGMA# contains a linear infrared divergence. Even if regulated with a nonzero photon mass M of order of the Debye mass, this infrared sensitivity implies that the two-loop self-energy contributes terms to the fermion dispersion relation that are comparable to or even larger than the next-to-leading order (NLO) contributions of the one-loop #SIGMA#. Additional evidence for the necessity of a systematic restructuring of the loop expansion comes from the explicit gauge-parameter dependence of the fermion damping rate at both one and two loops. The leading terms in the high temperature expansion of the two-loop self-energy for all topologies arise from an explicit ...
2010-01-15
Convolution/superposition using the Monte Carlo method
International Nuclear Information System (INIS)
The convolution/superposition calculations for radiotherapy dose distributions are traditionally performed by convolving polyenergetic energy deposition kernels with TERMA (total energy released per unit mass) precomputed in each voxel of the irradiated phantom. We propose an alternative method in which the TERMA calculation is replaced by random sampling of photon energy, direction and interaction point. Then, a direction is randomly sampled from the angular distribution of the monoenergetic kernel corresponding to the photon energy. The kernel ray is propagated across the phantom, and energy is deposited in each voxel traversed. An important advantage of the explicit sampling of energy is that spectral changes with depth are automatically accounted for. No spectral or kernel hardening corrections are needed. Furthermore, the continuous sampling of photon direction allows us to model sharp changes in fluence, such as those due to collimator tongue-and-groove. The ...
2003-07-21
Energy Technology Data Exchange (ETDEWEB)
We have recently introduced several important improvements in the measurement of distillation curves for complex fluids. This new method is a significant improvement over current approaches, with features including a composition-explicit data channel for each distillate fraction (for both qualitative and quantitative analysis), trace chemical analysis of each distillate fraction, and a corrosivity assessment of each distillate fraction. In the present contribution, we present the application of the advanced distillation curve method to a bio-derived crude oil prepared from swine manure. One motivation behind the work was to precisely measure the distillation curve of the oil derived from swine manure with our advanced distillation apparatus and compare the bio-derived oil to traditional petroleum crude oils. Then, the information content of the bio-derived oil distillation was extended further with the composition-explicit data channel: gas ...
2008-11-15
A revised dosimetric model of the adult head and brain
Energy Technology Data Exchange (ETDEWEB)
During the last decade, new radiopharmaceutical have been introduced for brain imaging. The marked differences of these tracers in tissue specificity within the brain and their increasing use for diagnostic studies support the need for a more anthropomorphic model of the human brain and head. Brain and head models developed in the past have been only simplistic representations of this anatomic region. For example, the brain within the phantom of MIRD Pamphlet No. 5 Revised is modeled simply as a single ellipsoid of tissue With no differentiation of its internal structures. To address this need, the MIRD Committee established a Task Group in 1992 to construct a more detailed brain model to include the cerebral cortex, the white matter, the cerebellum, the thalamus, the caudate nucleus, the lentiform nucleus, the cerebral spinal fluid, the lateral ventricles, and the third ventricle. This brain model has been included within a slightly modified version of the head model developed by ...
1996-06-01
A revised dosimetric model of the adult head and brain
International Nuclear Information System (INIS)
During the last decade, new radiopharmaceutical have been introduced for brain imaging. The marked differences of these tracers in tissue specificity within the brain and their increasing use for diagnostic studies support the need for a more anthropomorphic model of the human brain and head. Brain and head models developed in the past have been only simplistic representations of this anatomic region. For example, the brain within the phantom of MIRD Pamphlet No. 5 Revised is modeled simply as a single ellipsoid of tissue With no differentiation of its internal structures. To address this need, the MIRD Committee established a Task Group in 1992 to construct a more detailed brain model to include the cerebral cortex, the white matter, the cerebellum, the thalamus, the caudate nucleus, the lentiform nucleus, the cerebral spinal fluid, the lateral ventricles, and the third ventricle. This brain model has been included within a slightly modified version of the head model developed by ...
1996-07-21
{ital Ab Initio} Pseudopotential calculations of dopant diffusion in Si
Energy Technology Data Exchange (ETDEWEB)
The ab initio pseudopotential method is used to study transient-enhanced-diffusion (TED) related processes. The electronic degrees of freedom are included explicitly, together with the fully self-consistent treatment of the electron charge density. A large supercell and a fine k-point mesh are used to ensure numerical convergence. Such method has been demonstrated to give quantitative description of defect energetic. We will show that boron diffusion is significantly enhanced in the presence of the Si interstitial due to the substantial lowering of the migrational barrier through a kick-out mechanism. The resulting mobile boron can also be trapped by another substitutional boron, forming an immobile and elect rically inactive two-boron pair. Similarly, carbon diffusion is also enhanced significantly due to the pairing with Si interstitial. However, carbon binds to Si interstitial much more strongly than boron does, taking away most Si interstitial from boron at ...
1997-04-28
Vanishing integrals for Hall-Littlewood polynomials
It is well known that if one integrates a Schur function indexed by a partition $\\lambda$ over the symplectic (resp. orthogonal) group, the integral vanishes unless all parts of $\\lambda$ have even multiplicity (resp. all parts of $\\lambda$ are even). In a recent paper of Rains and Vazirani, Macdonald polynomial generalizations of these identities and several others were developed and proved using Hecke algebra techniques. However at $q=0$ (the Hall-Littlewood level), these approaches do not work, although one can obtain the results by taking the appropriate limit. In this paper, we develop a direct approach for dealing with this special case. This technique allows us to prove some identities that were not amenable to the Hecke algebra approach, as well as to explicitly control the nonzero values. Moreover, we are able to generalize some of the identities by introducing extra parameters. This leads us to a finite-dimensional analog of a recent result of Warnaar, ...
2010-01-01
Type of brown coal from North-Bohemian basin effect to quality of pyrolysis products
Energy Technology Data Exchange (ETDEWEB)
By increasing the price of crude oil and natural gas and the decline in their stock a possibility of more intense use of coal reserves as a feedstock not only for energy purposes but also as a feedstock source for other industries, mainly chemical and steel industry, is opening up. At present brown coal exploited in the Czech Republic is explicitly burned in the power sources of various outputs. Coal pyrolysis under the defined conditions may be an outstanding way of non-energetic brown coal processing. The target of our work was to determine the behaviour of various types of brown coal from the North Bohemian Basin under the pyrolysis up to the temperature 750{sup o}C. At the same time, qualitative parameters of the pyrolysis products in dependence on the features of the input brown coal were observed. Various types of brown coal were used for the pyrolysis tests differing mainly in the content of ash, sulphur and tar. Within the work mass balances of the incurred ...
2005-07-01
We solve the loop equations of the $\\beta$-ensemble model analogously to the solution found for the Hermitian matrices $\\beta=1$. For \\beta=1$, the solution was expressed using the algebraic spectral curve of equation $y^2=U(x)$. For arbitrary $\\beta$, the spectral curve converts into a Schr\\"odinger equation $((\\hbar\\partial)^2-U(x))\\psi(x)=0$ with $\\hbar\\propto (\\sqrt\\beta-1/\\sqrt\\beta)/N$. This paper is similar to the sister paper~I, in particular, all the main ingredients specific for the algebraic solution of the problem remain the same, but here we present the second approach to finding a solution of loop equations using sectorwise definition of resolvents. Being technically more involved, it allows defining consistently the B-cycle structure of the obtained quantum algebraic curve (a D-module of the form $y^2-U(x)$, where $[y,x]=\\hbar$) and to construct explicitly the correlation functions and the corresponding symplectic invariants $F_h$, or ...
2010-01-01
We use a superspin Hamiltonian defined on an infinite-dimensional Fock space with positive definite scalar product to study localization and delocalization of noninteracting spinless quasiparticles in quasi-one-dimensional quantum wires perturbed by weak quenched disorder. Past works using this approach have considered a single chain. Here, we extend the formalism to treat a quasi-one-dimensional system: a quantum wire with an arbitrary number of channels coupled by random hopping amplitudes. The computations are carried out explicitly for the case of a chiral quasi-one-dimensional wire with broken time-reversal symmetry (chiral-unitary symmetry class). By treating the space direction along the chains as imaginary time, the effects of the disorder are encoded in the time evolution induced by a single site superspin (non-Hermitian) Hamiltonian. We obtain the density of states near the band center of an infinitely long quantum wire. Our results agree with those based ...
2009-01-01
The Role of the UK Planning System in Protecting and Enhancing Soils
Environmental Research Database
Summary'THE ROLE OF THE UK PLANNING SYSTEM IN PROTECTING AND ENHANCING SOILS PROJECT REF: UKLQ01 APRIL 2004 INTRODUCTION Soil has long been overlooked in environmental planning. Whilst it has sometimes benefited from indirect protection, until recently its importance and vulnerability have not been reflected in more explicit strategies. This study explored the current and potential role of the planning system in protecting and enhancing soils. Soil is the physical material that covers much o [continued...]DescriptionSNIFFER has commissioned this study on behalf of SEPA, SNH, EHS and the Environment Agency for England and Wales. These UK agencies have an interest in protecting soil quality through a combination of direct regulation of emissions to land, protection of designated sites of special value or sensitivity and engaging other regulatory regimes such as land use planning. All UK environment agencies devote considerable resources to engaging the planning ...
2003-01-31
The Mass Spectra, Hierarchy and Cosmology of B-L MSSM Heterotic Compactifications
The matter spectrum of the MSSM, including three right-handed neutrino supermultiplets and one pair of Higgs-Higgs conjugate superfields, can be obtained by compactifying the E_{8} x E_{8} heterotic string and M-theory on Calabi-Yau manifolds with specific SU(4) vector bundles. These theories have the standard model gauge group augmented by an additional gauged U(1)_{B-L}. Their minimal content requires that the B-L gauge symmetry be spontaneously broken by a vacuum expectation value of at least one right-handed sneutrino. In previous papers, we presented the results of a quasi-analytic renormalization group analysis showing that B-L gauge symmetry is indeed radiatively broken with an appropriate B-L/electroweak hierarchy. In this paper, we extend these results by 1) enlarging the initial parameter space and 2) explicitly calculating all renormalization group equations numerically, without approximation. The regions of the initial parameter space leading to ...
2010-01-01
The Hydrodynamic Environment for the s Process in the He-Shell Flash of AGB Stars
The He-shell flash convection in AGB stars is the site for the high-temperature component of the s-process in low- and intermediate mass giants, driven by the Ne22 neutron source. [...] The upper convection boundary plays a critical role during the H-ingestion episode that may lead to neutron-bursts in the most metal-poor AGB stars. We address these problems through global 3-dimensional hydrodynamic simulations including the entire spherical He-shell flash convection zone (as oposed to the 3D box-in-a-star simulations). An important aspect of our current effort is to establish the feasibility of our appoach. We explain why we favour the explicit treatment over the anelastic approximation for this problem. The simulations presented in this paper use a Cartesian grid of 512^3 cells and have been run on four 8-core workstations for four days to simulate ~5000s, which corresponds to almost ten convective turn-over times. The convection layer extends radially at the ...
2009-01-01
Environmental Research Database
ObjectivesThe overall aim of the proposed research is explicit in the project title, i.e. the creation of a map of current vertical land movements in the UK based on an optimal combination of absolute gravity (AG) and continuous GPS (CGPS). This is consistent with specific objective (ii) of WP1.9 of the Oceans 2025 programme and is related to priority topic area 4 (application of satellite geodesy to sea level science) of the NERC Strategic Ocean Funding Initiative (SOFI). From long term geological an [continued...]DescriptionThe proposed research aims to create a map of current vertical land movements in the UK based on an optimal combination of estimates from two geodetic surveying and monitoring techniques; the measurement of absolute gravity (AG) and the use of high precision, continuous GPS (CGPS) observations. From long term geological and geophysical studies, vertical land movements in the UK are thought to be of the order of 1 to 2 millimetres per year, ...
2009-01-31
International Nuclear Information System (INIS)
We study the lipid and phase transferability of our recently developed systematically coarse-grained solvent-free membrane model. The force field was explicitly parameterized to describe a fluid 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayer at 310 K with correct structure and area per lipid, while gaining at least three orders of magnitude in computational efficiency (see Wang and Deserno 2010 J. Phys. Chem. B 114 11207-20). Here, we show that exchanging CG tails, without any subsequent re-parameterization, creates reliable models of 1,2-dioleoylphosphatidylcholine (DOPC) and 1,2-dipalmitoylphosphatidylcholine (DPPC) lipids in terms of structure and area per lipid. Furthermore, all CG lipids undergo a liquid-gel transition upon cooling, with characteristics like those observed in experiments and all-atom simulations during phase transformation. These studies suggest a promising transferability of our force field parameters to different lipid species and ...
2010-09-01
Structure activity relationships to assess new chemicals under TSCA
International Nuclear Information System (INIS)
Under Section 5 of the Toxic Substances Control Act (TSCA), manufacturers must notify the US Environmental Protection Agency (EPA) 90 days before manufacturing, processing, or importing a new chemical substance. This is referred to as a premanufacture notice (PMN). The PMN must contain certain information including chemical identity, production volume, proposed uses, estimates of exposure and release, and any health or environmental test data that are available to the submitter. Because there is no explicit statutory authority that requires testing of new chemicals prior to their entry into the market, most PMNs are submitted with little or no data. As a result, EPA has developed special techniques for hazard assessment of PMN chemicals. These include (1) evaluation of available data on the chemical itself, (2) evaluation of data on analogues of the PMN, or evaluation of data on metabolites or analogues of metabolites of the PMN, (3) use of quantitative structure ...
1990-06-27
Strategic principles workshops: Discussion drafts and workshop notes
Energy Technology Data Exchange (ETDEWEB)
The Office of Civilian Radioactive Waste-Management in the Department of Energy (DOE) is responsible for disposing of this nation`s spent fuel and high-level radioactive waste in a manner that protects the health and safety of the public and the quality of the environment. Although embodied in the Federal repository program that began with studies in the late 1950s, this mission was explicitly established by the Nuclear Waste Policy Act of 1982 and reaffirmed by the Nuclear Waste Policy Amendments Act of 1987. To fulfill our mission, we are developing a waste management system consisting of a geologic repository for permanent disposal deep beneath the surface of the earth, a facility for monitored retrievable storage, and a system for transporting the waste. This discussion draft was developed to help involve parties affected by or interested in the waste-management program in the formulation of the basic principles on which the program will be based. It reviews ...
1991-09-01
Stochastic-convective transport with nonlinear reaction: Mathematical framework
Energy Technology Data Exchange (ETDEWEB)
A stochastic-convective reactive (SCR) transport method is developed for one-dimensional steady transport in physically heterogeneous media with nonlinear degradation. The method is free of perturbation amplitude limitations and circumvents the difficulty of scale dependence of phenomenological parameters by avoiding volume-averaged specifications of diffusive/dispersive fluxes. The transport system is conceptualized as an ensemble of independent convective-reactive streamlines, each characterized by a randomized convective velocity (or travel time). Dispersive effects are treated as a component of the randomness in the streamline velocity ensemble, so no explicit expression for hydrodynamic dispersive flux is written in the streamline transport equation. The expected value of the transport over the stream tube ensemble is obtained as an average of solutions to the reactive convection equation according to the stream tube (travel time) probability distribution ...
1995-11-01
International Nuclear Information System (INIS)
In this paper, molecular quantum computation is numerically studied with the quantum search algorithm (Grover's algorithm) by means of optimal control simulation. Qubits are implemented in the vibronic states of I_2, while gate operations are realized by optimally designed laser pulses. The methodological aspects of the simulation are discussed in detail. We show that the algorithm for solving a gate pulse-design problem has the same mathematical form as a state-to-state control problem in the density matrix formalism, which provides monotonically convergent algorithms as an alternative to the Krotov method. The sequential irradiation of separately designed gate pulses leads to the population distribution predicted by Grover's algorithm. The computational accuracy is reduced by the imperfect quality of the pulse design and by the electronic decoherence processes that are modeled by the non-Markovian master equation. However, as long as we focus on the population distribution of the ...
2010-04-01
Semirelativistic technique for k#centre dot#p calculations: Optical properties of Pd and Pt
International Nuclear Information System (INIS)
A semirelativistic two-component extended linear augmented plane-wave k#centre dot#p method is described. In order to ensure a high accuracy of the k#centre dot#p method, it is necessary to include into the radial-basis set, which is used for the augmentation of the plane waves, functions that are neither solutions of the Schrodinger equation nor their energy derivatives. The usual scalar relativistic procedure, which is nonlinear in energy, is not applicable to such basis sets. As an alternative, we suggest an approximation to the Foldy-Wouthuysen Hamiltonian that produces an explicitly Hermitean matrix in the augmented plane wave representation. The technique is applied to the calculation of the full dielectric matrix and optical properties of palladium and platinum metals over the photon energy region up to 100 eV. Special attention is paid to the far ultraviolet absorption by the excitations of semicore Pd 4p and Pt 5p and 4f states. A strong effect of local ...
2001-06-15
We consider the (2, 0) supersymmetric theory of tensor multiplets and self-dual strings in six space-time dimensions. Space-time diffeomorphisms that leave the string world-sheet invariant appear as gauge transformations on the normal bundle of the world-sheet. The naive invariance of the model under such transformations is however explicitly broken by anomalies: The electromagnetic coupling of the string to the two-form gauge field of the tensor multiplet suffers from a classical anomaly, and there is also a one-loop quantum anomaly from the chiral fermions on the string world-sheet. Both of these contributions are proportional to the Euler class of the normal bundle of the string world-sheet, and consistency of the model requires that they cancel. This imposes strong constraints on possible models, which are found to obey an ADE-classification. We then consider the decoupled world-sheet theory that describes low-energy fluctuations (compared to the scale set by ...
2004-01-01
Relativistic mean-field approach to nuclear surface properties and spin-orbit effects
Energy Technology Data Exchange (ETDEWEB)
We treat symmetric semi-infinite nuclear matter in the relativistic mean-field approximation for the scalar-vector field theoretical model. Using special-type Dirac spinors the nucleonic Dirac equation is decoupled into two sets of differential equations for the spin-orientation dependent orbital nucleon Dirac spinors. We also rewrite the Dirac equation in terms of second-order differential equations with the spin-orbit interaction appearing explicitly. These equations can be solved if the spin-orbit part is left out. The spin-orbit effects thus can be isolated, and are shown to reduce the surface energy coefficient a/sub s/ as well as the surface thickness t in such a way that their ratio remains practically unchanged. For realistic lagrangians - in linear as well as non-linear forms - consistent with the empirical spin-orbit single-partial level splittings, a/sub s/ is reduced by nearly 10%. We also discuss nuclear surface properties and their relations to bulk ...
1989-02-20
We apply the permutation symmetry S3 to both charged-lepton and neutrino mass matrices, and suggest a useful symmetry-breaking scheme, in which the flavor symmetry is explicitly broken down via S3 -> Z3 -> nothing in the charged-lepton sector and via S3 -> Z2 -> nothing in the neutrino sector. Such a two-stage breaking scenario is reasonable in the sense that both Z3 and Z2 are the subgroups of S3, while Z3 and Z2 only have a trivial subgroup. In this scenario, we can naturally obtain a relatively large value of the smallest neutrino mixing angle, e.g., theta13 ~ 9 degrees, which is compatible with the recent result from T2K experiment and will be precisely measured in the ongoing Double Chooz and Daya Bay reactor neutrino experiments. Moreover, the maximal atmospheric mixing angle theta23 ~ 45 degrees can also be obtained while the best-fit value of solar mixing angle theta12 ~ 34 degrees is assumed, which cannot be achieved in previous S3 symmetry ...
2011-01-01
This study reports on an examination of the effectiveness of a reading intervention for adults with disabilities in a vocational rehabilitation setting. Participants were 57 adults with disabilities and low reading skills enrolled at the Reading Clinic at the Michigan Career and Technical Institute. As part of a 3-year research and demonstration grant evaluation project, participants received an individualized reading intervention that targeted phonological processing, orthographic pattern recognition, and comprehension for adults with disabilities and low reading levels. The particular program incorporated fluency exercises at every level of the structured, explicit, and systematic intervention curriculum. Results showed that participants made moderate to large gains in passage reading accuracy and comprehension during the course of the intervention. Gains in passage reading rate were not statistically significant. Results are discussed in relation to theories of ...
2011-07-01
We study an exactly solvable model where an uniformly accelerated detector is linearly coupled to a massless scalar field initially in the Minkowski vacuum. Using the exact correlation functions we show that as soon as the coupling is switched on one can see information flowing from the detector to the field and propagating with the radiation into null infinity. By expressing the reduced density matrix of the detector in terms of the two-point functions, we calculate the purity function in the detector and study the evolution of quantum entanglement between the detector and the field. Only in the ultraweak coupling regime could some degree of recoherence in the detector appear at late times, but never in full restoration, as an earlier work seems to suggest. We explicitly show that under the most general conditions the detector never recovers its quantum coherence and the entanglement between the detector and the field remains large at late times. To the extent ...
2007-01-01
PRA and Risk Informed Analysis
The Boiler and Pressure Vessel Code (BPVC) of the American Society of Mechanical Engineers (ASME) has introduced a risk based approach into Section XI that covers Rules for Inservice Inspection of Nuclear Power Plant Components. The risk based approach requires application of the probabilistic risk assessments (PRA). Because no industry consensus standard existed for PRAs, ASME has developed a standard to evaluate the quality level of an available PRA needed to support a given risk based application. The paper describes the PRA standard, Section XI application of PRAs, and plans for broader applications of PRAs to other ASME nuclear codes and standards. The paper addresses several specific topics of interest to Section XI. Important consideration are special methods (surrogate components) used to overcome the lack of PRA treatments of passive components in PRAs. The approach allows calculations of conditional core damage probabilities both for component failures that cause ...
2006-01-01
PRA and Risk Informed Analysis
Energy Technology Data Exchange (ETDEWEB)
The Boiler and Pressure Vessel Code (BPVC) of the American Society of Mechanical Engineers (ASME) has introduced a risk based approach into Section XI that covers Rules for Inservice Inspection of Nuclear Power Plant Components. The risk based approach requires application of the probabilistic risk assessments (PRA). Because no industry consensus standard existed for PRAs, ASME has developed a standard to evaluate the quality level of an available PRA needed to support a given risk based application. The paper describes the PRA standard, Section XI application of PRAs, and plans for broader applications of PRAs to other ASME nuclear codes and standards. The paper addresses several specific topics of interest to Section XI. Important consideration are special methods (surrogate components) used to overcome the lack of PRA treatments of passive components in PRAs. The approach allows calculations of conditional core damage probabilities both for component failures that cause initiating ...
2006-01-01
PRA and Risk Informed Analysis
International Nuclear Information System (INIS)
The Boiler and Pressure Vessel Code (BPVC) of the American Society of Mechanical Engineers (ASME) has introduced a risk based approach into Section XI that covers Rules for Inservice Inspection of Nuclear Power Plant Components. The risk based approach requires application of the probabilistic risk assessments (PRA). Because no industry consensus standard existed for PRAs, ASME has developed a standard to evaluate the quality level of an available PRA needed to support a given risk based application. The paper describes the PRA standard, Section XI application of PRAs, and plans for broader applications of PRAs to other ASME nuclear codes and standards. The paper addresses several specific topics of interest to Section XI. Important consideration are special methods (surrogate components) used to overcome the lack of PRA treatments of passive components in PRAs. The approach allows calculations of conditional core damage probabilities both for component failures that cause initiating ...
Energy Technology Data Exchange (ETDEWEB)
Removing noise from data is often the first step in data analysis. Denoising techniques should not only reduce the noise, but do so without blurring or changing the location of the edges. Many approaches have been proposed to accomplish this; in this paper, they focus on one such approach, namely the use of non-linear diffusion operators. This approach has been studied extensively from a theoretical viewpoint ever since the 1987 work of Perona and Malik showed that non-linear filters outperformed the more traditional linear Canny edge detector. They complement this theoretical work by investigating the performance of several isotropic diffusion operators on test images from scientific domains. They explore the effects of various parameters such as the choice of diffusivity function, explicit and implicit methods for the discretization of the PDE, and approaches for the spatial discretization of the non-linear operator etc. They also compare these schemes with ...
2001-12-20
Open vs Closed Access Femtocells in the Uplink
Femtocells are assuming an increasingly important role in the coverage and capacity of cellular networks. In contrast to existing cellular systems, femtocells are end-user deployed and controlled, randomly located, and rely on third party backhaul (e.g. DSL or cable modem). Femtocells can be configured to be either open access or closed access. Open access allows an arbitrary nearby cellular user to use the femtocell, whereas closed access restricts the use of the femtocell to users explicitly approved by the owner. Seemingly, the network operator would prefer an open access deployment since this provides an inexpensive way to expand their network capabilities, whereas the femtocell owner would prefer closed access, in order to keep the femtocell's capacity and backhaul to himself. We show mathematically and through simulations that the reality is more complicated for both parties, and that the best approach depends heavily on whether the multiple access scheme is ...
2010-01-01
On Sums of Generating Sets in (Z_2)^n
Let A and B be two affinely generating sets of (Z_2)^n. As usual, we denote their Minkowski sum by A+B. How small can A+B be, given the cardinalities of A and B? We give a fairly tight answer to this question. Our bound is attained when both A and B are unions of cosets of a certain subgroup of (Z_2)^n. These cosets are arranged as Hamming balls, the smaller of which has radius 1. By similar methods, we reprove the Freiman-Ruzsa theorem in (Z_2)^n, with an optimal upper bound. Denote by F(K) the maximal spanning constant || / |A|, over all subsets A of (Z_2)^n with doubling constant |A+A| / |A| < K. We explicitly calculate F(K), and in particular show that 4^K / 4K < F(K) (1+o(1)) < 4^K / 2K. This improves the estimate F(K) = poly(K) 4^K, found recently by Green and Tao and by Konyagin.
2011-01-01
Mining associations in text in the presence of background knowledge
Energy Technology Data Exchange (ETDEWEB)
This paper describes the FACT system for knowledge discovery from text. It discovers associations - patterns of co-occurrence - amongst keywords labeling the items in a collection of textual documents. In addition, FACT is able to use background knowledge about the keywords labeling the documents in its discovery process. FACT takes a query-centered view of knowledge discovery, in which a discovery request is viewed as a query over the implicit set of possible results supported by a collection of documents, and where background knowledge is used to specify constraints on the desired results of this query process. Execution of a knowledge-discovery query is structured so that these back-ground-knowledge constraints can be exploited in the search for possible results. Finally, rather than requiring a user to specify an explicit query expression in the knowledge-discovery query language, FACT presents the user with a simple-to-use graphical interface to the query ...
1996-12-31
Gastric dilatation-volvulus (GDV) is a life-threatening condition in dogs and other species in which the stomach dilates and rotates on itself. The etiology of the disease is multi-factorial, but explicit precipitating causes are unknown. This study sought to determine if there was a significant association between changes in hourly-measured temperature and/or atmospheric pressure and the occurrence of GDV in the population of high-risk working dogs in Texas. The odds of a day being a GDV day, given certain temperature and atmospheric pressure conditions for that day or the day before, was estimated using logistic regression models. There were 57 days in which GDV(s) occurred, representing 2.60% of the days in the 6-year study period. The months of November, December, and January collectively accounted for almost half (47%) of all cases. Disease risk was negatively associated with daily maximum temperature. An increased risk of GDV was weakly associated with the ...
2008-01-01
Metal Coupled Folding of Cys2His2 Zinc-Finger
Zinc-fingers, which widely exist in eukaryotic cell and play crucial roles in life processes, depend on the binding of zinc ion for their proper folding. To computationally study the zinc coupled folding of the zinc-fingers, charge transfer and metal induced protonation/deprotonation effects have to be considered. Here, by attempting to implicitly account for such effects in classical molecular dynamics and performing intensive simulations with explicit solvent for the peptides with and without zinc binding, we investigate the folding of the Cys2His2 type zinc-finger motif and the coupling between the peptide folding and zinc binding. We find that zinc ion not only stabilizes the native structure, but also participates in the whole folding process. It binds to the peptide at early stage of folding, and directs or modulates the folding and stabilizations of the component beta-hairpin and alpha-helix. Such a crucial role of zinc binding is mediated by the packing of ...
2008-01-01
International Nuclear Information System (INIS)
A one-dimensional numerical model for the expansion of impact-produced vapor clouds is used to investigate magnetic field generation mechanisms in events such as meteor collisions with the moon. The resulting cloud properties, such as ionization fraction, electrical conductivity, radial expansion velocity, mass density, and energy density are estimated. The model is initiated with the peak shock states and pressure thresholds for incipient and complete vaporization of anorthosite lunar surface materials by iron and GA composition meteorites. The expansion of the spherical gas cloud into a vacuum was traced with a one-dimensional explicit lagrangian hydrodynamic code. The hypervelocity impact plasmas produced are found to be significant in the amplitudes and orientations of the magnetic fields generated. An ambient magnetic field could have been provided by the core dynamo, which would have interacted with the expanding plasmas and formed induced paleomagnetic ...
1984-03-12
Linearized Schroedinger equation for nuclear quadrupole surface vibrations
Energy Technology Data Exchange (ETDEWEB)
The Schroedinger equation for nuclear quadruple surface vibrations is linearized with the consequence that a new spin degree of freedom appears in the wave function of the linearized equation. This spin is called collective spin and has a value of 3/2. The linearized Schroedinger equation for quadrupole vibrations is used for the description of certain collective aspects of even-odd {sup 187,189,191}Ir nuclei which have a spin 3/2 in their ground state. As a potential we use the {gamma}-soft collective potential of the neighboring even-even nuclei, which is inserted into the linearized Schroedinger equation via a scalar coupling. This leads to a collective spin-dependent fine structure splitting of the energy levels governed by a collective SO(5) spin-orbit coupling and a correction to the kinetic energy. Further, we consider explicitly spin-dependent potentials which effectively describe the interaction of the valence nucleon with the core of the even-odd Ir ...
1990-07-01
Linearized Schroedinger equation for nuclear quadrupole surface vibrations
International Nuclear Information System (INIS)
The Schroedinger equation for nuclear quadruple surface vibrations is linearized with the consequence that a new spin degree of freedom appears in the wave function of the linearized equation. This spin is called collective spin and has a value of 3/2. The linearized Schroedinger equation for quadrupole vibrations is used for the description of certain collective aspects of even-odd "1"8"7","1"8"9","1"9"1Ir nuclei which have a spin 3/2 in their ground state. As a potential we use the #gamma#-soft collective potential of the neighboring even-even nuclei, which is inserted into the linearized Schroedinger equation via a scalar coupling. This leads to a collective spin-dependent fine structure splitting of the energy levels governed by a collective SO(5) spin-orbit coupling and a correction to the kinetic energy. Further, we consider explicitly spin-dependent potentials which effectively describe the interaction of the valence nucleon with the core of the even-odd Ir ...
Inverse polynomial optimization
We consider the inverse optimization problem associated with the polynomial program f^*=\\min \\{f(x): x\\in K\\}$ and a given current feasible solution $y\\in K$. We provide a systematic numerical scheme to compute an inverse optimal solution. That is, we compute a polynomial $\\tilde{f}$ (which may be of same degree as $f$ if desired) with the following properties: (a) $y$ is a global minimizer of $\\tilde{f}$ on $K$ with a Putinar's certificate with an a priori degree bound $d$ fixed, and (b), $\\tilde{f}$ minimizes $\\Vert f-\\tilde{f}\\Vert$ (which can be the $\\ell_1$, $\\ell_2$ or $\\ell_\\infty$-norm of the coefficients) over all polynomials with such properties. Computing $\\tilde{f}_d$ reduces to solving a semidefinite program whose optimal value also provides a bound on how far is $f(\\y)$ from the unknown optimal value $f^*$. The size of the semidefinite program can be adapted to the computational capabilities available. Moreover, if one uses the $\\ell_1$-norm, then ...
2011-01-01
Intra-cellular transport of single-headed molecular motors KIF1A
Motivated by recent experiments on KIF1A, a representative member of single-headed kinesin motor proteins family, we develop a theoretical model of intra-cellular transport by mutually interacting molecular motors. The model explicitly accounts not only for the hydrolysis of ATP, but also for the ratchet mechanism which is believed to drive each individual KIF1A motor. We study the model by a combination of analytical and numerical techniques. A remarkable feature of this model is that all the parameters in it can be completely evaluated from experimental data. Our results in the dilute limit are in excellent quantitative agreement with the empirical data from single molecule experiments. In the high density regime the predictions of the model also agree qualitatively with the corresponding experimental observations. We derive a phase diagram that shows the influence of hydrolysis and Langmuir kinetics on the collective spatio-temporal organization of the motors. ...
2005-01-01
Interacting tachyons in classical and quantum physics
International Nuclear Information System (INIS)
It is demonstrated that tachyons do not violate the principles of relativity, and that, with the aid of a reinterpretation principle to eliminate negative energies, tachyons can be characterized as particles of real, spacelike 4-momentum. The classical, charged tachyon is treated within conventional electromagnetic theory, and in an explicitly Lorentz-invariant way. It is shown that a charged tachyon would not emit electromagnetic radiation in a vacuum regardless of its state of motion. A theory based on the real-energy solutions of the Klein-Gordon equation with imaginary mass is shown to provide the best opportunity for describing spinless tachyons in quantum field theory. The theory should be Lorentz-invariant, incorporate the reinterpretation principle to remove negative energies, and be as close as possible to conventional quantum theory. The proposal of Arons and Sudarshan is adopted as best fulfilling these requirements. A perturbation-type expansion for the ...
Hydrogen synthesis via combustion of fuel-rich natural gas/air mixtures at elevated pressure
Energy Technology Data Exchange (ETDEWEB)
Combustion of extremely fuel-rich ({phi}=4) methane/air mixtures at elevated pressures is investigated as a potential means to generate molecular hydrogen by non-catalytic partial oxidation. This system is investigated both computationally and experimentally. The computations use a perfectly-stirred reactor model and an explicit methane cool-flame mechanism to investigate the effects of reactor parameters on reaction time and product composition. Under adiabatic conditions, such mixtures are predicted to autoignite at low temperatures {approx}700 K for pressures exceeding 8.5 atm. Above 15 atm, conversion to products is complete in roughly 1 s. The dependence of reaction time and hydrogen yield is investigated as a function of inlet temperature, system pressure, and flame equivalence ratio. Actual product yields are measured in a tube reactor facility, and many of the predictions of the model, including long relaxation times to chemical equilibrium and temperatures ...
2005-07-01
Gluon Plasma Frequency -- the Next-to-Leading Order Term
The longitudinal-electric oscillations of the hot gluon system are studied beyond the well known leading order term at high temperature $T$ and small coupling $g$. The coefficient $\\eta$ in $\\omega^2 = m^2 \\, (1+ \\eta \\, g \\wu N \\, )$ is calculated, where \\hbox{$\\omega \\equiv \\omega (\\vc q =0)$} is the long-wavelength limit of the frequency spectrum, $N$ the number of colours and $m^2=g^2 N T^2/9$. In the course of this, for the real part of the gluon self-energy, the Braaten-Pisarski resummation programme is found to work well in all details. The coefficient $\\eta$ is explicitly seen to be gauge independent within the class of covariant gauges. Infrared singularities cancel as well as collinear singularities in the two-loop diagrams with both inner momenta hard. However, as it turns out, none of these two-loop contributions reaches the relative order $O(g)$ under study. The minus sign in our numerical result $\\; \\eta = -.18 \\; $ is in accord with ...
1994-01-01
Gauge invariant perturbation theory in spatially homogeneous cosmology
International Nuclear Information System (INIS)
The perturbations of the L.R.S. class A spatially homogeneous spacetimes are treated using Hamiltonian methods in conjunction with techniques from the theory of Lie group harmonic analysis. These latter techniques lead to a simple way of handling any set of tensor equations on these background spacetimes which has the same symmetry group as the space-time metric. Once this approach is developed, the Hamiltonian formulation is used to recover in a clean way the Bonanos equations for the perturbations of the perfect fluid models of the class of spacetimes under consideration. The conserved quantities associated with the four-dimensional symmetry group are evaluated and their role in the linearized Hamiltonian dynamics is discussed. The time-dependent linear canonical transformation of the linearized vacuum gravitational phase space adapted to Moncrief's gauge invariant decomposition is described in general for these models and evaluated explicitly for a class of ...
Fuzzy decision making in multiobjective long-term scheduling of hydrothermal system
Energy Technology Data Exchange (ETDEWEB)
In the paper, a fuzzy decision-making methodology is presented to decide the generation schedule of long-term hydrothermal problems with explicit recognition of statistical uncertainties in system production cost data, NO{sub x} emission data, system load demand and hydro reservoir water inflows. In deciding the optimal operation, three objectives operating cost, NO{sub x} emission and unsatisfied load demand over the whole of the planning period are simultaneously minimised. Specific technique is put forth to convert the stochastic models into their deterministic equivalents. The weighted minimax method is used to simulate the tradeoff relation between the conflicting objectives in the non-inferior domain. The fuzzy set theory is exploited to choose the best operating point over the tradeoff curve. An efficient decomposition technique is applied to reduce the complexity of the problem. In each subproblem, thermal generations are obtained by using simplified method ...
2001-01-01
The Bingham plastic constitutive model has been widely used to predict the post-yield behavior of electro- and magneto- rheological fluids (ER and MR fluids). However, if these fluids experience shear thinning or shear thickening, the Bingham plastic model may not be an accurate predictor of behavior, since the post-yield plastic viscosity is assumed to be constant. In a recent study, it was theoretically and experimentally demonstrated that the Herschel-Bulkley fluid model can be successfully employed when evaluating non- Newtonian post-yield behavior of ER and MR fluids. In this paper, the Herschel-Bulkley model is employed to include a detailed analysis of ER and MR fluid dynamics through pipes and parallel plates. Simplified explicit expressions for the exact formulation are also developed. It is shown that the proposed simplified model of the Herschel-Bulkley steady flow equations for pipes and parallel plates can be used as an accurate design tool while ...
2000-04-01
FastFlow: Efficient Parallel Streaming Applications on Multi-core
Shared memory multiprocessors come back to popularity thanks to rapid spreading of commodity multi-core architectures. As ever, shared memory programs are fairly easy to write and quite hard to optimise; providing multi-core programmers with optimising tools and programming frameworks is a nowadays challenge. Few efforts have been done to support effective streaming applications on these architectures. In this paper we introduce FastFlow, a low-level programming framework based on lock-free queues explicitly designed to support high-level languages for streaming applications. We compare FastFlow with state-of-the-art programming frameworks such as Cilk, OpenMP, and Intel TBB. We experimentally demonstrate that FastFlow is always more efficient than all of them in a set of micro-benchmarks and on a real world application; the speedup edge of FastFlow over other solutions might be bold for fine grain tasks, as an example +35% on OpenMP, +226% on Cilk, +96% on TBB for ...
2009-01-01
Extension of the EQ3/6 computer codes to geochemical modeling of brines
Energy Technology Data Exchange (ETDEWEB)
Recent modifications to the EQ3/6 geochemical modeling software package provide for the use of Pitzer's equations to calculate the activity coefficients of aqueous species and the activity of water. These changes extend the range of solute concentrations over which the codes can be used to dependably calculate equilibria in geochemical systems, and permit the inclusion of ion pairs, complexes, and undissociated acids and bases as explicit component species in the Pitzer model. Comparisons of calculations made by the EQ3NR and EQ6 compuer codes with experimental data confirm that the modifications not only allow the codes to accurately evaluate activity coefficients in concentrated solutions, but also permit prediction of solubility limits of evaporite minerals in brines at 25/sup 0/C and elevated temperatures. Calculations for a few salts can be made at temperatures up to approx. 300/sup 0/C, but the temperature range for most electrolytes is constrained ...
1984-10-23
We show how to obtain a fast component-by-component construction algorithm for higher order polynomial lattice rules. Such rules are useful for multivariate quadrature of high-dimensional smooth functions over the unit cube as they achieve the near optimal order of convergence. The main problem addressed in this paper is to find an efficient way of computing the worst-case error. A general algorithm is presented and explicit expressions for base~2 are given. To obtain an efficient component-by-component construction algorithm we exploit the structure of the underlying cyclic group. We compare our new higher order multivariate quadrature rules to existing quadrature rules based on higher order digital nets by computing their worst-case error. These numerical results show that the higher order polynomial lattice rules improve upon the known constructions of quasi-Monte Carlo rules based on higher order digital nets.
2011-01-01
Dose, exposure time, and resolution in Serial X-ray Crystallography
The resolution of X-ray diffraction microscopy is limited by the maximum dose that can be delivered prior to sample damage. In the proposed Serial Crystallography method, the damage problem is addressed by distributing the total dose over many identical hydrated macromolecules running continuously in a single-file train across a continuous X-ray beam, and resolution is then limited only by the available molecular and X-ray fluxes and molecular alignment. Orientation of the diffracting molecules is achieved by laser alignment. We evaluate the incident X-ray fluence (energy/area) required to obtain a given resolution from (1) an analytical model, giving the count rate at the maximum scattering angle for a model protein, (2) explicit simulation of diffraction patterns for a GroEL-GroES protein complex, and (3) the frequency cut off of the transfer function following iterative solution of the phase problem, and reconstruction of an electron density map in the ...
2007-01-01
Density functional calculations of 15N chemical shifts in solvated dipeptides
International Nuclear Information System (INIS)
We performed density functional calculations to examine the effects of solvation, hydrogen bonding, backbone conformation, and the side chain on 15N chemical shielding in proteins. We used N-methylacetamide (NMA) and N-formyl-alanyl-X (with X being one of the 19 naturally occurring amino acids excluding proline) as model systems. In addition, calculations were performed for selected fragments from protein GB3. The conducting polarizable continuum model was employed to include the effect of solvent in the density functional calculations. Our calculations for NMA show that the augmentation of the polarizable continuum model with the explicit water molecules in the first solvation shell has a significant influence on isotropic 15N chemical shift but not as much on the chemical shift anisotropy. The difference in the isotropic chemical shift between the standard ?-sheet and ?-helical conformations ranges from 0.8 to 6.2 ppm depending on the residue type, with the mean ...
2008-06-01
Collective behavior of stock price movements in an emerging market
To investigate the universality of the structure of interactions in different markets, we analyze the cross-correlation matrix ${\\mathbf C}$ of stock price fluctuations in the National Stock Exchange (NSE) of India. We find that this {\\em emerging} market exhibits strong correlations in the movement of stock prices compared to {\\em developed} markets, such as the New York Stock Exchange (NYSE). This is shown to be due to the dominant influence of a common market mode on the stock prices. By comparison, interactions between related stocks, e.g., those belonging to the same business sector, are much weaker. This lack of distinct sector identity in emerging markets is explicitly shown by reconstructing the network of mutually interacting stocks. Spectral analysis of ${\\mathbf C}$ for NSE reveals that, the few largest eigenvalues deviate from the bulk of the spectrum predicted by random matrix theory, but they are far fewer in number compared to, e.g., NYSE. We ...
2007-01-01
Coherent transport of matter waves in disordered optical potentials
Energy Technology Data Exchange (ETDEWEB)
The development of modern techniques for the cooling and the manipulation of atoms in recent years, and the possibility to create Bose-Einstein condensates and degenerate Fermi gases and to load them into regular optical lattices or disordered optical potentials, has evoked new interest for the disorder-induced localization of ultra-cold atoms. This work studies the transport properties of matter waves in disordered optical potentials, which are also known as speckle potentials. The effect of correlated disorder on localization is first studied numerically in the framework of the Anderson model. The relevant transport parameters in the configuration average over many different realizations of the speckle potential are then determined analytically, using self-consistent diagrammatic perturbation techniques. This allows to make predictions for a possible experimental observation of coherent transport phenomena for cold atoms in speckle potentials. Of particular importance are the spatial ...
2007-07-01
CSI: A Paradigm for Behavior-oriented Delivery Services in Mobile Human Networks
We propose behavior-oriented services as a new paradigm of communication in mobile human networks. Our study is motivated by the tight user-network coupling in future mobile societies. In such a paradigm, messages are sent to inferred behavioral profiles, instead of explicit IDs. Our paper provides a systematic framework in providing such services. First, user behavioral profiles are constructed based on traces collected from two large wireless networks, and their spatio-temporal stability is analyzed. The implicit relationship discovered between mobile users could be utilized to provide a service for message delivery and discovery in various network environments. As an example application, we provide a detailed design of such a service in challenged opportunistic network architecture, named CSI. We provide a fully distributed solution using behavioral profile space gradients and small world structures. Our analysis shows that user behavioral profiles are ...
2008-01-01
CFD Simulations of Pb-Bi Two-Phase Flow
International Nuclear Information System (INIS)
In a Pb-Bi cooled direct contact steam generation fast reactor water is injected directly above the core, the produced steam is separated at the top and is send to the turbine. Neither the direct contact phenomenon nor the two-phase flow simulations in CFD have been thoroughly described yet. A first attempt in simulating such two-phase flow in 2D using the CFD code Fluent is presented in this paper. The volume of fluid explicit model was used. Other important simulation parameters were: pressure velocity relation PISO, discretization scheme body force weighted for pressure, second order upwind for momentum and CISCAM for void fraction. Boundary conditions were mass flow inlet (Pb-Bi 0 kg/s and steam 0.07 kg/s) and pressure outlet. The effect of mesh size (0.5 mm and 0.2 mm cells) was investigated as well as the effect of the turbulent model. It was found that using a fine mesh is very important in order to achieve larger bubbles and the turbulent model (k-? ...
2008-09-21
CANDU 9 - the CANDU product to meet customer and regulator requirements now and in the future
Energy Technology Data Exchange (ETDEWEB)
CANDU reactors developed under Canadian licensing regulations that placed the primary responsibility for safety on the licensee. The Atomic Energy Control Board (AECB), Canada's nuclear regulatory agency, state in their regulations what is expected in terms of safety performance so that designers are free to propose the best means of meeting this performance. This goal-oriented approach, besides encouraging innovation, allowed CANDU to be licensed in other jurisdictions. The latest design - the large, single unit, CANDU 9 - explicitly incorporates licensability in Canada through a formal AECB review of the design; lessons learned from licensing CANDU 6 in Asian countries, particularly with Wolsong 2, 3 and 4 in Korea, and more recently with Qinshan in China; utility requirements for modem evolutionary plants; and emerging international standards for safety, sponsored or issued by the IAEA. By combining the assurance of acceptability in Canada with ...
1998-07-01
CANDU 9 - the CANDU product to meet customer and regulator requirements now and in the future
International Nuclear Information System (INIS)
CANDU reactors developed under Canadian licensing regulations that placed the primary responsibility for safety on the licensee. The Atomic Energy Control Board (AECB), Canada's nuclear regulatory agency, state in their regulations what is expected in terms of safety performance so that designers are free to propose the best means of meeting this performance. This goal-oriented approach, besides encouraging innovation, allowed CANDU to be licensed in other jurisdictions. The latest design - the large, single unit, CANDU 9 - explicitly incorporates licensability in Canada through a formal AECB review of the design; lessons learned from licensing CANDU 6 in Asian countries, particularly with Wolsong 2, 3 and 4 in Korea, and more recently with Qinshan in China; utility requirements for modem evolutionary plants; and emerging international standards for safety, sponsored or issued by the IAEA. By combining the assurance of acceptability in Canada with compliance with ...
1998-05-03
Energy Technology Data Exchange (ETDEWEB)
Debris coolability in the lower plenum of the reactor pressure vessel is an important factor for evaluation of debris in-vessel retention. The debris coolability analysis module is developed for the accurate prediction of the safety margin of the present reactor vessels in a severe accident. The module calculates debris spreading and cooling through melting and solidification in combination with a temperature distribution and failure evaluation of the vessel wall. Debris spreading is solved by the explicit method on a quasi-three-dimensional scheme and debris coolability is solved on the basis of natural convection analysis. The calculation for spreading is compared with a water spreading experiment on the floor and the calculation for coolability is compared with a n-octadecane melting experiment in a rectangular vessel. The comparisons show capabilities for predictions of spearhead transportation in the debris spreading process and of melting front transportation ...
1999-07-01
International Nuclear Information System (INIS)
Debris coolability in the lower plenum of the reactor pressure vessel is an important factor for evaluation of debris in-vessel retention. The debris coolability analysis module is developed for the accurate prediction of the safety margin of the present reactor vessels in a severe accident. The module calculates debris spreading and cooling through melting and solidification in combination with a temperature distribution and failure evaluation of the vessel wall. Debris spreading is solved by the explicit method on a quasi-three-dimensional scheme and debris coolability is solved on the basis of natural convection analysis. The calculation for spreading is compared with a water spreading experiment on the floor and the calculation for coolability is compared with a n-octadecane melting experiment in a rectangular vessel. The comparisons show capabilities for predictions of spearhead transportation in the debris spreading process and of melting front transportation ...
1999-04-19
Antideuteron fluxes from dark matter annihilation in diffusion models
Antideuterons are among the most promising galactic cosmic ray-related targets for dark matter indirect detection. Currently only upper limits exist on the flux, but the development of new experiments, such as GAPS and AMS-02, provides exciting perspectives for a positive measurement in the near future. In this Paper, we present a novel and updated calculation of both the secondary and primary antideuteron fluxes. We employ a two-zone diffusion model which successfully reproduces cosmic-ray nuclear data and the observed antiproton flux. We review the nuclear and astrophysical uncertainties and provide an up to date secondary (i.e. background) antideuteron flux. The primary (i.e. signal) contribution is calculated for generic WIMPs annihilating in the galactic halo: we explicitly consider and quantify the various sources of uncertainty in the theoretical evaluations. Propagation uncertainties, as is the case of antiprotons, are sizeable. Nevertheless, antideuterons ...
2008-01-01
International Nuclear Information System (INIS)
The connections between several studies of the off-shell amplitude, based on apparently different criteria, are first clarified and expressed through more coincise operator notion. In all cases the resolution of the underlying dynamical relations is reduced to a two-step procedure. Underconditions which are discussed, the latter implies only Volterra-like integral equations followed by simple quadratures and algebraic operations. Then it is shown that the off-shell generalizations of the Jost formalism which are defined by such approaches can as well be introduced without reference to any explicit dynamical framework. Examples are given specifying in such model independent ways the main properties of the associated Jost-type functions. The relative interest of different representations built with these functions is also examined, and on this occasion new three-separable-term approximations are a useful guideline for introducing similar formalisms in other ...
Energy Technology Data Exchange (ETDEWEB)
This paper presents a new technique--Integrated Bayesian Uncertainty Estimator (IBUNE) to account for the major uncertainties of hydrologic rainfall-runoff predictions explicitly. The uncertainties from the input (forcing) data--mainly the precipitation observations and from the model parameters are reduced through a Monte Carlo Markov Chain (MCMC) scheme named Shuffled Complex Evolution Metropolis (SCEM) algorithm which has been extended to include a precipitation error model. Afterwards, the Bayesian Model Averaging (BMA) scheme is employed to further improve the prediction skill and uncertainty estimation using multiple model output. A series of case studies using three rainfall-runoff models to predict the streamflow in the Leaf River basin, Mississippi are used to examine the necessity and usefulness of this technique. The results suggests that ignoring either input forcings error or model structural uncertainty will lead to unrealistic model simulations and ...
2006-05-05
ADIABATIC MASS LOSS AND THE OUTCOME OF THE COMMON ENVELOPE PHASE OF BINARY EVOLUTION
International Nuclear Information System (INIS)
We have developed a new method for calculating common envelope (CE) events based on explicit consideration of the donor star's structural response to adiabatic mass loss. In contrast to existing CE prescriptions, which specify a priori the donor's remnant mass, we determine this quantity self-consistently and find that it depends on binary and CE parameters. This aspect of our model is particularly important to realistic modeling for upper main-sequence star donors without strongly degenerate cores (and hence without a clear core/envelope boundary). We illustrate the central features of our method by considering CE events involving 10 M_s_u_n donors on or before their red giant branch. For such donors, the remnant core mass can be as much as 30% larger than the star's He-core mass. Applied across a population of such binaries, our methodology results in a significantly broader remnant mass and final orbital separation distribution and a 20% increase in CE survival ...
2010-08-10
A model of coherent fluctuations of nuclear density
International Nuclear Information System (INIS)
A model for coherent fluctuations of the nuclear density including all nucleons in the fluctuation movement is formulated. The wave function of the flucton satisfies a Schroedinger-type equation in which the role of potential energy is played by the energy of the nuclear matter for the A nucleons. The appearance of high-momentum components in the one-particle momentum distribution is related to the volume fluctuations of the density. These fluctuations are supposed to be of coherent nature, i.e. all nucleons are involved in the motion. An explicit expression for the momentum distribution of the nucleons in "1"2C and "1"8"1Ta nuclei is obtained and compared with the phenomenological momentum distribution as well as with the results given by the model of a non-relativistic harmonic oscillator and by the Dirac equation in the self-consistent model. The differential cross-section of 1 GeV proton scattering by the "2"8Si, "3"2S, "4"0Ca, "4"8Ca, "5"8Ni, "2"0"8Pb are ...
A detailed physical model for ion implant induced damage in silicon
Energy Technology Data Exchange (ETDEWEB)
A unified physically based ion implantation damage model has been developed which successfully predicts both the impurity profiles and the damage profiles for a wide range of implant conditions for arsenic, phosphorus, BF{sub 2}, and boron implants into single-crystal silicon. In addition, the amorphous layer thicknesses predicted by this new damage model are also in excellent agreement with experimental measurements. This damage model is based on the physics of point defects in silicon, and explicitly simulates the defect production, diffusion, and their interactions which include interstitial-vacancy recombination, clustering of same type of defects, defect-impurity complex formation, emission of mobile defects from clusters, and surface effects for the first time. New computationally efficient algorithms have been developed to overcome the barrier of the excessive computational requirements. In addition, the new model has been incorporated in the UT-MARLOWE ion ...
1998-06-01
A detailed physical model for ion implant induced damage in silicon
International Nuclear Information System (INIS)
A unified physically based ion implantation damage model has been developed which successfully predicts both the impurity profiles and the damage profiles for a wide range of implant conditions for arsenic, phosphorus, BF_2, and boron implants into single-crystal silicon. In addition, the amorphous layer thicknesses predicted by this new damage model are also in excellent agreement with experimental measurements. This damage model is based on the physics of point defects in silicon, and explicitly simulates the defect production, diffusion, and their interactions which include interstitial-vacancy recombination, clustering of same type of defects, defect-impurity complex formation, emission of mobile defects from clusters, and surface effects for the first time. New computationally efficient algorithms have been developed to overcome the barrier of the excessive computational requirements. In addition, the new model has been incorporated in the UT-MARLOWE ion ...
1998-06-01
A New Measurement of Kaonic Hydrogen X rays
The $\\bar{K}N$ system at threshold is a sensitive testing ground for low energy QCD, especially for the explicit chiral symmetry breaking. Therefore, we have measured the $K$-series x rays of kaonic hydrogen atoms at the DA$\\Phi$NE electron-positron collider of Laboratori Nazionali di Frascati, and have determined the most precise values of the strong-interaction energy-level shift and width of the $1s$ atomic state. As x-ray detectors, we used large-area silicon drift detectors having excellent energy and timing resolution, which were developed especially for the SIDDHARTA experiment. The shift and width were determined to be $\\epsilon_{1s} = -283 \\pm 36 \\pm 6 {(syst)}$ eV and $\\Gamma_{1s} = 541 \\pm 89 {(stat)} \\pm 22 {(syst)}$ eV, respectively. The new values will provide vital constraints on the theoretical description of the low-energy $\\bar{K}N$ interaction.
2011-01-01
On the Role of Convection and Turbulence for Tropospheric Ozone and its Precursors
International Nuclear Information System (INIS)
The aim of the work in this thesis is to investigate the convective and diffusive transport in the TM chemistry transport model, and to investigate some aspects of the consequences for NOx. The large inaccuracy and uncertainty in the description of processes like convection and turbulent diffusion, the strong dependence of the radiative forcing of ozone on its vertical distribution, and the strong dependence of the ozone production on the distribution of NOx, are the main motivation. The availability of the ERA-40 data, where convective data and vertical diffusion coefficients are archived, allows a study of the effect of different convective mass flux sets, and different vertical diffusion coefficients on the model-simulated distribution of tracers. In this thesis the following questions are addressed : (1) How large is the sensitivity of the (model simulated) distribution of ozone and nitrogen oxides on (the) convection (parameterisation)?; (2) What requirements should be fulfilled ...
The controversy about '1/m{sub Q} duality violation': a quark model point of view
Energy Technology Data Exchange (ETDEWEB)
Non relativistic quark models have been invoked to support the statement of a 1/m{sub Q} duality violation in semileptonic B decays. However, as we recall, we have very explicitly shown that no 1/m{sub Q} duality violation is present in totally integrated rates, in such quark models. Then: 1) it is shown that some contributions under discussion are misleadingly described as '1/m{sub Q} duality violation'; as they stand, they are in fact parametrically much smaller: they are O((1/m{sub Q}){sup 5/2}) when properly referred to the total rate under discussion; numerically this is below 10{sup -2}; 2) moreover, relying on our rigorous expansion of the harmonic oscillator model around the Shifman-Voloshin limit, it is shown that even such small terms are not present at all in the total rate, and must therefore merely cancel; 3) finally, for physical masses, by an exact treatment of this particular, but not too unrealistic model, one finds a departure ...
2006-06-15
The contribution of the DOE`s R&D budget in natural gas to energy price security
Energy Technology Data Exchange (ETDEWEB)
The energy price volatility model suggests that some of the proposed natural gas programs can contribute to energy price stability. The sector most vulnerable to fuel price variations is, of course, the transportation sector. The most effective strategy to achieve energy pace stability is to reduce petroleum consumption in this sector. The natural gas vehicle program is therefore recommended as potentially important and worthy of further consideration. At this point, distinguishing the merits of various subprograms is not feasible. This result farther supports the conclusion that the DOE`s energy R&D portfolio is not efficiently balanced and an increase in oil and gas research should be a high priority. The DOE has responded favorably and has significantly increased its proposed research with the explicit objective of displacing oil in the transportation sector. The enhanced research and development program for energy security, in the NES, proposes major ...
1992-12-01
The analysis of complex fluids such as crude oils, fuels, vegetable oils and mixed waste streams poses significant challenges arising primarily from the multiplicity of components, the different properties of the components (polarity, polarizability, etc.) and matrix properties. We have recently introduced an analytical strategy that simplifies many of these analyses, and provides the added potential of linking compositional information with physical property information. This aspect can be used to facilitate equation of state development for the complex fluids. In addition to chemical characterization, the approach provides the ability to calculate thermodynamic properties for such complex heterogeneous streams. The technique is based on the advanced distillation curve (ADC) metrology, which separates a complex fluid by distillation into fractions that are sampled, and for which thermodynamically consistent temperatures are measured at atmospheric pressure. The collected sample ...
2009-11-17
Second-law analyses applied to internal combustion engines operation
Energy Technology Data Exchange (ETDEWEB)
This paper surveys the publications available in the literature concerning the application of the second-law of thermodynamics to internal combustion engines. The availability (exergy) balance equations of the engine cylinder and subsystems are reviewed in detail providing also relations concerning the definition of state properties, chemical availability, flow and fuel availability, and dead state. Special attention is given to identification and quantification of second-law efficiencies and the irreversibilities of various processes and subsystems. The latter being particularly important since they are not identified in traditional first-law analysis. In identifying these processes and subsystems, the main differences between second- and first-law analyses are also highlighted. A detailed reference is made to the findings of various researchers in the field over the last 40 years concerning all types of internal combustion engines, i.e. spark ignition, compression ignition (direct or ...
2006-07-01
Production and decay of scalar top squarkonium bound states
International Nuclear Information System (INIS)
In this paper we discuss possible signatures for the production of scalar t_1t_1"* (top squarkonium) bound states #sigma#_t_1 at hadron colliders, where t_1 is the lighter scalar top eigenstate. We first study the decay of #sigma#_t_1; explicit expressions are given for all potentially important decay modes. If t_1 has unsuppressed two-body decays, they will always overwhelm the annihilation decays of #sigma#_t_1. Among the latter, we find that usually either the gg or hh final state dominates, depending on the size of the off-diagonal entry of the top squark mass matrix; h is the lighter neutral scalar Higgs boson of the minimal supersymmetric model. If m_#sigma#_t happens to be close to the mass of one of the neutral scalar Higgs bosons, Q bar Q final states dominate (Q=b or t). W"+W"- and ZZ final states are subdominant. We argue that #sigma#_t_1#->##gamma##gamma# decays offer the best signal for top squarkonium production at hadron colliders. The Fermilab ...
Performance objectives for the Hanford immobilized low-activity waste (ILAW) performance assessment
Energy Technology Data Exchange (ETDEWEB)
Before low-level waste may be disposed of, a performance assessment must be written and then approved by the DOE (DOE 1988a DOE 1999a). The performance assessment is to determine whether ''reasonable assurance'' exists that the performance objectives of the disposal facility will be met. The DOE requirements for waste disposal (DOE 1988a DOE 1999a) require the protection of public health and safety; and the protection of the environment. Although quantitative limits are sometimes stated (for example, the all-pathways exposure limit is 25 mrem/year), usually the requirements are stated in a general nature. Quantitative limits were established by: investigating all potentially applicable regulations as well as interpretations of the review panels which DOE has established to review performance assessments, interacting with program management to establish the additional requirements of the program, and interacting with the public (i.e., the Hanford ...
1999-09-09
Energy Technology Data Exchange (ETDEWEB)
This lectures aim at giving graduate students an introduction to a working knowledge of path integral methods in a wide variety of fields in physics. Consequently, the the lecture notes are organized in three main parts dealing with non-relativistic quantum mechanics, many-body physics and field theory. In the first part the basic concepts of path integrals are developed in the usual heuristic, non-mathematical way followed by the standard examples of quadratic Lagrangians for which the path integrals can be solved exactly. Applications include semi-classical expansions, scattering problems and the representation of Green functions as path integrals. In the last chapter of this part it is shown how (euclidean) path integrals can be treated numerically by Monte-Carlo methods with a program for the anharmonic oscillator as an explicit example. The second part deals with the application of path integrals in statistical mechanics and many-body problems. Various ...
1997-12-01
Optimizing an analytical dose calculation algorithm for fast 2D calculations
International Nuclear Information System (INIS)
Previously, an analytical dose calculation algorithm for MLC-based radiotherapy was developed and commissioned, which includes a detailed model of various MLC effects as a unique feature [1]. The algorithm was originally developed as an independent verification of the treatment planning system's dose calculation and it explicitly modeled spatial and depth dependent MLC effects such as interleaf transmission, the tongue-and-groove effect, rounded leaf ends, MLC scatter, beam hardening, and gradual MLC transmission fall-off with increasing off-axis distance. Originally the algorithm was implemented in Mathematica trademark (Wolfram). To speed up the calculation time and to be able to calculate high resolution 2D dose distributions within a reasonable time frame (<2 s) the algorithm needs to be optimized and to be embedded in a user friendly environment. To achieve this goal, the dose calculation model is implemented in Visual Basic 6.0, which decreases the ...
On foundation of the generalized Nambu mechanics
International Nuclear Information System (INIS)
We outline basic principles of a canonical formalism for the Nambu mechanics - a generalization of Hamiltonian mechanics proposed by Yoichiro Nambu in 1973. It is based on the notion of a Nambu bracket, which generalizes the Poisson bracket - a 'binary'' operation on classical observables on the phase space - to the 'multiple' operation of higher order n#>=#3. Nambu dynamics is described by the phase flow given by Nambu-Hamilton equations of motion - a system of ODE's which involves n-1 'Hamiltonians'. We introduce the fundamental identity for the Nambu bracket - a generalization of the Jacobi identity - as a consistency condition for the dynamics. We show that Nambu bracket structure defines a hierarchy of infinite families of 'subordinated' structures of lower order, including Poisson bracket structure, which satisfy certain matching conditions. The notion of Nambu bracket enables us to define Nambu-Poisson manifolds - phase spaces for the Nambu mechanics, which turn out to be ...
International Nuclear Information System (INIS)
How should lightweight electric-powered vehicles be promoted, after the large-scale fleet test being conducted at Mendrisio (southern Switzerland) is completed in 2001, and are there reasons to put question marks behind the current approach? The demand for electric vehicles, and particularly the one in the automobile category, has remained at a persistently low level. As it proved, any appreciable improvement of this situation is almost impossible, even with substantial financial incentives. However, the unsatisfactory sales figures have little to do with the nature of the fleet test itself or with the specific conditions at Mendrisio. The problem is rather of structural nature. For (battery-operated) electric cars the main problem at present is the lack of an expanding market which could become self-supporting with only a few additional incentives. Various strategies have been evaluated. Two alternatives were considered in particular: a strategy to promote ...
Introducing DeBRa: a detailed breast model for radiological studies
Energy Technology Data Exchange (ETDEWEB)
Currently, x-ray mammography is the method of choice in breast cancer screening programmes. As the mammography technology moves from 2D imaging modalities to 3D, conventional computational phantoms do not have sufficient detail to support the studies of these advanced imaging systems. Studies of these 3D imaging systems call for a realistic and sophisticated computational model of the breast. DeBRa (Detailed Breast model for Radiological studies) is the most advanced, detailed, 3D computational model of the breast developed recently for breast imaging studies. A DeBRa phantom can be constructed to model a compressed breast, as in film/screen, digital mammography and digital breast tomosynthesis studies, or a non-compressed breast as in positron emission mammography and breast CT studies. Both the cranial-caudal and mediolateral oblique views can be modelled. The anatomical details inside the phantom include the lactiferous duct system, the Cooper ligaments and the pectoral muscle. The ...
2009-07-21
Influence of anchor behaviour on the earthquake response of liquid storage tanks
International Nuclear Information System (INIS)
The dynamic response of thin liquid storage tanks to earthquakes is a very complicated phenomenon, because it can be highly non linear. Among others, one can meet material and geometric non linearities of the tank shell leading eventually to static or dynamic buckling non linear behavior of anchor bolts, contact non-linearities due to the uplift of the tank base and to the unilateral character of the fluid pressure on the shell and high amplitude fluid oscillations. Moreover, linear or non linear soil structure interaction affects considerably the response of the fluid structure system under consideration. In this paper we focus attention on problems related only to the base uplift and anchors plastification. We study a tank similar to the Hualien project tank, but we neglect the soil structure interaction. The studied tank is representative of medium height to radius ratio ratio tanks with relatively thick bottom plate. The contact is simulated via a simple discrete penalty method in ...
Infinite Dimensional Groups and Riemann Surface Field Theories
International Nuclear Information System (INIS)
We show how to obtain positive energy representations of the group G of smooth maps from a union of circles to U(N) from geometric data associated with a Riemann surface having these circles as boundary. Using covering spaces we can reduce to the case where N=1. Then our main result shows that Mackey induction may be applied and yields representations of the connected component of the identity of G which have the form of a Fock representation of an infinite dimensional Heisenberg group tensored with a finite dimensional representation of a subgroup isomorphic to the first cohomology group of the surface obtained by capping the boundary circles with discs. We give geometric sufficient conditions for the correlation functions to be positive definite and derive explicit formulae for them and for the vacuum (or cyclic) vector. (This gives a geometric construction of correlation functions which had been obtained earlier using tau functions.) By choosing particular ...
1996-01-01
Implementation of FFT convolution and multigrid superposition models in the FOCUS RTP system
International Nuclear Information System (INIS)
In radiotherapy treatment planning, convolution/superposition algorithms currently represent the best practical approach for accurate photon dose calculation in heterogeneous tissues. In this work, the implementation, accuracy and performance of the FFT convolution (FFTC) and multigrid superposition (MGS) algorithms are presented. The FFTC and MGS models use the same 'TERMA' calculation and are commissioned using the same parameters. Both models use the same spectra, incorporate the same off-axis softening and base incident lateral fluence on the same measurements. In addition, corrections are explicitly applied to the polyenergetic and parallel kernel approximations, and electron contamination is modelled. Spectra generated by Monte Carlo (MC) modelling of treatment heads are used. Calculations using the MC spectra were in excellent agreement with measurements for many linear accelerator types. To speed up the calculations, a number of calculation techniques were ...
2000-04-01
Energy Technology Data Exchange (ETDEWEB)
A new three-dimensional (3D) acoustic modelling method was developed using a first-order hyperbolic wave system which was solved with explicit finite dfferences. The numerical solution of the 3D wave system provides a useful method for simulating evolution of a pressure field corresponding to compressional type waves. Existing two-dimensional (2D) elastic modelling algorithms were modified and fine-tuned for computationally efficient and realistic wave propagation simulations in complex structures. An original formulation of the 3D reverse time migration method was developed which is very accurate, does not suffer from unwanted evenescent energy, can image dips beyond 90{degree}, and does not generate multiple energy. Two case studies were performed that involved steam stimulation projects in the Cold Lake deposit. Simulations were performed during different phases of the steam stimulation process to examine the relation between reservoir properties and conditions ...
1994-12-31
Fundamental Elements of Geologic C02 Sequestration in Saline Aquifers
Energy Technology Data Exchange (ETDEWEB)
Geologic sequestration represents a promising strategy for isolating CO{sub 2} waste streams from the atmosphere. Successful implementation of this approach hinges on our ability to predict the relative effectiveness of subsurface CO{sub 2} migration and sequestration as a function of key target-formation and cap-rock properties, which will enable us to identify optimal sites and evaluate their long-term isolation performance. Quantifying this functional relationship requires a modeling capability that explicitly couples multiphase flow and kinetically controlled geochemical processes. We have developed a unique computational package that meets these criteria, and used it to model CO{sub 2} injection at Statoil's North-Sea Sleipner facility, the world's first saline-aquifer storage site. The package integrates a state-of-the-art reactive transport simulator (NUFT) with supporting geochemical software and databases (SUPCRT92). In our Sleipner ...
2001-11-19
International Nuclear Information System (INIS)
The objective was to study the brain areas which are activated when normal subjects make moral judgments. Ten normal adults underwent BOLD functional magnetic resonance imaging (fMRI) during the auditory presentation of sentences that they were instructed to silently judge as either 'right' or 'wrong'. Half of the sentences had an explicit moral content ('We break the law when necessary'), the other half comprised factual statements devoid of moral connotation ('Stones are made of water'). After scanning, each subject rated the moral content, emotional valence, and judgment difficulty of each sentence on Likert-like scales. To exclude the effect of emotion on the activation results, individual responses were hemo dynamically modeled for event-related f MRI analysis. The general linear model was used to evaluate the brain areas activated by moral judgment. Regions activated during moral judgment included the frontopolar cortex (FPC), medial frontal gyrus, right ...
2001-09-01
Energy Technology Data Exchange (ETDEWEB)
We have constructed a Monte Carlo generator (the corresponding FORTRAN code can be obtained from the authors upon request) for lowest-order predictions for the processes {gamma}{gamma}{yields}4f and {gamma}{gamma}{yields}4f{gamma} in the standard model and extensions thereof by an effective {gamma}{gamma}H coupling as well as anomalous triple and quartic gauge-boson couplings. Polarization is fully supported, and a realistic photon beam spectrum can be taken into account. For the processes {gamma}{gamma}{yields}4f all helicity amplitudes are explicitly given in a compact form. The presented numerical results contain, in particular, a survey of cross sections for representative final states and their comparison to results obtained with the program package Whizard/Madgraph. The impact of a realistic beam spectrum on cross sections and distributions is illustrated. Moreover, the size of various contributions to cross sections, such as from weak charged- or ...
2004-08-01
International Nuclear Information System (INIS)
We have constructed a Monte Carlo generator (the corresponding FORTRAN code can be obtained from the authors upon request) for lowest-order predictions for the processes #gamma##gamma##->#4f and #gamma##gamma##->#4f#gamma# in the standard model and extensions thereof by an effective #gamma##gamma#H coupling as well as anomalous triple and quartic gauge-boson couplings. Polarization is fully supported, and a realistic photon beam spectrum can be taken into account. For the processes #gamma##gamma##->#4f all helicity amplitudes are explicitly given in a compact form. The presented numerical results contain, in particular, a survey of cross sections for representative final states and their comparison to results obtained with the program package Whizard/Madgraph. The impact of a realistic beam spectrum on cross sections and distributions is illustrated. Moreover, the size of various contributions to cross sections, such as from weak charged- or ...
2004-08-01
Field simulation of axisymmetric plasma screw pinches by alternating-direction-implicit methods
Energy Technology Data Exchange (ETDEWEB)
An axisymmetric plasma screw pinch is an axisymmetric column of ionized gaseous plasma radially confined by forces from axial and azimuthal currents driven in the plasma and its surroundings. This dissertation is a contribution to detailed, high resolution computer simulation of dynamic plasma screw pinches in 2-d {ital rz}-coordinates. The simulation algorithm combines electron fluid and particle-in-cell (PIC) ion models to represent the plasma in a hybrid fashion. The plasma is assumed to be quasineutral; along with the Darwin approximation to the Maxwell equations, this implies application of Ampere`s law without displacement current. Electron inertia is assumed negligible so that advective terms in the electron momentum equation are ignored. Electrons and ions have separate scalar temperatures, and a scalar plasma electrical resistivity is assumed. Altemating-direction-implicit (ADI) methods are used to advance the electron fluid drift velocity and the magnetic fields in the ...
1996-06-01
Electron-ion recombination of neutral iron
Energy Technology Data Exchange (ETDEWEB)
The total and state-specific electron-ion recombination rate coefficients are obtained for FeI. The calculations are carried out using a new ab initio method that incorporates both the radiative and the dielectronic recombination processes in an unified and self-consistent manner. The computations employ the close coupling approximation and the R-matrix method from atomic collision theory. A 52 state close coupling eigenfunction expansion dominated by the states of the ground 3d{sup 6}4s and excited 3d{sup 7}, 3d{sup 6}4p, 3d{sup 5}4s{sup 2}, and 3d{sup 5}4s4p configurations of FeII are used in the present calculations. The important electron correlation and radiation damping effects are included via explicit coupling of autoionization and radiative channels. This is the first detailed atomic calculation for the recombination rates for FeI. The present rates are considerably higher than the radiative recombination rates being used currently in the low-temperature ...
1997-04-01
Diffusion of adatoms on face-centered cubic transition metal surfaces
Energy Technology Data Exchange (ETDEWEB)
Mechanisms and associated energetics for adatom diffusion on the (100) and (110) surfaces of Ni, Cu, Rh, Pd, and Ag are investigated. Self-diffusion was studied on (100) and (I 10) surfaces of Ni, Cu, Pd and Ag using corrected effective medium method (CEM) and approximation to CEM used for molecular dynamics and Monte Carlo studies (MD/MC-CEM). Self-diffusion on Pd(100), Ag(100), Ni(110), Cu(110), Pd(110), and Ag(110) is accomplished by classical diffusion: the adatom hops from its equilibrium adsorption site over an intervening bridge site to an adjacent equilibrium site. Self-diffusion on Ni(100) and Cu(100) proceeds by atomic-exchange diffusion: the adatom on the surface displaces an atom in the first surface layer. Aside from explicit inclusion of the kinetic-exchange-correlation energy, it is critical to include enough movable atoms in the calculation to insure correct energetics. Distortions induced by these diffusion mechanisms, especially atomic exchange, ...
1994-05-10
DEFF Research Database (Denmark)
Executive summary 1. Few studies have to date explained consumer attitudes and purchase decisions with regard to genetically engineered food products. However, the increased marketing of genetically engineered food products and the considerable concern that consumers seem to express with regard to the technology call for the development of a theoretical basis for research into these issues. 2. The aim of the paper is to present three models which we have developed to explain consumer attitudes, buying behaviour and attitude change with regard to genetically engineered food products. All three models build on established consumer behaviour theory and on existing and comparable research in the field. 3. Consumer attitudes toward genetic engineering in food products are explained in an attitude model that builds on Fishbein's multiattribute attitude model. The model deviates from Fishbein's model in a number of ways: there is an explicit distinction between perceived ...
1998-01-01
Composing simulations using persistent software components
Energy Technology Data Exchange (ETDEWEB)
The traditional process for developing large-scale simulations is cumbersome, time consuming, costly, and in some cases, inadequate. The topics of software components and component-based software engineering are being explored by software professionals in academic and industrial settings. A component is a well-delineated, relatively independent, and replaceable part of a software system that performs a specific function. Many researchers have addressed the potential to derive a component-based approach to simulations in general, and a few have focused on military simulations in particular. In a component-based approach, functional or logical blocks of the simulation entities are represented as coherent collections of components satisfying explicitly defined interface requirements. A simulation is a top-level aggregate comprised of a collection of components that interact with each other in the context of a simulated environment. A component may represent a ...
1999-03-01
Comparison of helmet-mounted display designs in support of infantry wayfinding
The Canadian Soldier Information Requirements Technology Demonstration (SIREQ TD) soldier modernization research and development program has conducted experiments to help determine the types and amount of information needed to support wayfinding across a range of terrain environments, the most effective display modality for providing the information (visual, auditory or tactile) that will minimize conflict with other infantry tasks, and to optimize interface design. In this study, seven different visual helmet-mounted display (HMD) designs were developed based on soldier feedback from previous studies. The displays and an in-service compass condition were contrasted to investigate how the visual HMD interfaces influenced navigation performance. Displays varied with respect to their information content, frame of reference, point of view, and display features. Twelve male infantry soldiers used all eight experimental conditions to locate bearings to waypoints. From a constant location, ...
2003-09-01
Classical tachyons and possible applications
Energy Technology Data Exchange (ETDEWEB)
This article is divided in five parts, the first one having nothing to do with tachyons. In fact, to prepare the ground, in part I (sect.2) we shall merely show that special relativity - even without tachyons - can be given a form such to describe both particles and antiparticles. Part II is the largest one: initially, after some historical remarks and having revisited the postulates of special relativity, we presnt a review of the elegant ''model theory'' of tachyons in two dimensions; passing then to four dimensions, we review the main results of the classical theory of tachyons that do not depend on the existence of Superluminal reference frames (or that are at least independent of the explicit form of the Superluminal Lorentz ''transformations''). In particular, we discuss how tachyons would look like, i.e. their apparent ''shape''. Last but not least, ...
1986-01-01
Choice of energy data in environmental assessment of the built environment
Energy Technology Data Exchange (ETDEWEB)
Life cycle oriented methods are increasingly used for environmental assessments (EAs) of the built environment. However, many assumptions are made in such assessments, potentially influencing the results and making the assessment more ambiguous. To increase the reliability of EAs, consequences of the assumptions made have to be better understood. Since energy use in the operation and maintenance phase is an important factor decisive for the overall environmental performance of a building, the purpose of this study is to investigate how the selection of heat and electricity mix affects the assessed environmental performance of buildings. It also aims to suggest how to choose heat and electricity data in EAs of the built environment in general. Applying four different modes of electricity production and two different modes of heat production in a case study of three different buildings with different technical solutions for heat and electricity supply, the study show that choices of heat ...
2003-03-01
Analysis of Thermally Induced Permeability Enhancement in Geothermal Injection Wells
Energy Technology Data Exchange (ETDEWEB)
Reinjection of spent geothermal brine is a common means of disposing of geothermal effluents and maintaining reservoir pressures. Contrary to the predictions of two-fluid models (two-viscosity) of nonisothermal injection, an increase of injectivity, with continued injection, is often observed. Injectivity enhancement and thermally-affected pressure transients are particularly apparent in short-term injection tests at the Los Azufres Geothermal Field, Mexico. During an injection test, it is not uncommon to observe that after an initial pressure increase, the pressure decreases with time. As this typically occurs far below the pressure at which hydraulic fracturing is expected, some other mechanism for increasing the near-bore permeability must explain the observed behavior. This paper focuses on calculating the magnitude of the near-bore permeability changes observed in several nonisothermal injection tests conducted at the Los Azufres Geothermal Field. In order to evaluate the pressure ...
1987-01-20
The performance of intermolecular potential models on the adsorption of benzene on graphitized thermal carbon black at various temperatures is investigated. Two models contain only dispersive sites, whereas the other two models account explicitly for the dispersive and electrostatic sites. Using numerous data in the literature on benzene adsorption on graphitized thermal carbon black at various temperatures, we have found that the effect of surface mediation on interaction between adsorbed benzene molecules must be accounted for to describe correctly the adsorption isotherm as well as the isosteric heat. Among the two models with partial charges tested, the WSKS model of Wick et al. that has only six dispersive sites and three discrete partial charges is better than the very expensive all-atom model of Jorgensen and Severance. Adsorbed benzene molecules on graphitized thermal carbon black have a complex orientation with respect to distance from the surface and also ...
2006-01-31
A stochastic regime switching model for the failure process of a repairable system
Energy Technology Data Exchange (ETDEWEB)
This paper presents a stochastic model and estimation procedure for analyzing the failure process of a repairable system. We consider repairable systems whose successive interfailure times reveal a significant dependence while showing an insignificant trend. Neither the renewal process nor the non-homogeneous Poisson process are adequate for modeling such failure processes. Especially when the interfailure times show a cyclic pattern, we may consider a switching of the regimes (states) governing the lifetime distribution of the system. We propose a Markov switching model describing the failure process for such a case. The model postulates that a finite number of states governs the distinct lifetime distributions, and the state makes transitions according to a discrete-time Markov chain. Each of the distinct lifetime distributions represents a failure type that may change after successive repairs. Our model generalizes the mixture model by allowing the mixture probabilities to change ...
1998-02-01
Retrospective Monte Carlo dose calculations with limited beam weight information
International Nuclear Information System (INIS)
An important unresolved issue in outcomes analysis for lung complications is the effect of poor or completely lacking heterogeneity corrections in previously archived treatment plans. To estimate this effect, we developed a novel method based on Monte Carlo (MC) dose calculations which can be applied retrospectively to RTOG/AAPM-style archived treatment plans (ATP). We applied this method to 218 archived nonsmall cell lung cancer lung treatment plans that were originally calculated either without heterogeneity corrections or with primitive corrections. To retrospectively specify beam weights and wedges, beams were broken into Monte Carlo-generated beamlets, simulated using the VMC++ code, and mathematical optimization was used to match the archived water-based dose distributions. The derived beam weights (and any wedge effects) were then applied to Monte Carlo beamlets regenerated based on the patient computed tomography densities. Validation of the process was performed against five ...
2007-01-01
International Nuclear Information System (INIS)
Requirements on PSA for risk informed applications are expressed in different international documents. The ASME PRA standard published in spring 2002 is one such document, PSA requirements are also expressed in the European Utility Requirements (EUR) for new reactors. The Swedish PSA requirements are provided in the Swedish regulators (SKI) statutes SKIFS 2004:1. SKI also has a review handbook for PSA activities (SKI report 2003:48). The review handbook is a support during review of the utilities PSA activities and the PSAs themselves. The review handbook expresses SKIs expectations by providing so called important aspects for both the PSA work and the PSAs, A comparison of SKIFS requirements and the important aspects in the Review handbook, on one side, and the requirements on PSA in EUR and ASME on the other side, is presented. The comparison shows a large difference in the level of detail in the different documents, where ASME is most detailed and specific. This is expected since ...
1998-11-15
Apparent slip at the surface of a small rotating sphere in a dilute quiescent suspension
Energy Technology Data Exchange (ETDEWEB)
We consider the case of a test sphere (ball) of radius a{sub 1} rotating at constant angular velocity {omega} in an otherwise quiescent unbounded suspension of uniformly sized spheres of radii a{sub 2} dispersed in a Newtonian fluid of viscosity {mu}. To the first order in the volume fraction c of suspended spheres it is shown that when the ball is small compared with the suspended spheres the suspension does not behave as regards the hydrodynamic torque L exerted on the ball like a homogeneous Newtonian fluid characterized by the usual Einstein viscosity coefficient {mu}{sub s}={mu}(1+5/2c). Explicitly, the torque on the rotating sphere does not obey Kirchoff{close_quote}s law, L=8{pi}{mu}{sub s}a{sub 1}{sup 3}{omega} for no slip. Rather, a modified form of Kirchoff{close_quote}s law is obtained in which the Einstein coefficient of 5/2 is multiplied by a coefficient which is less than unity in magnitude and is functionally dependent only upon the ...
1998-03-01
Aespoe Task Force on modelling of groundwater flow and transport of solutes. Review of Task 6C
Energy Technology Data Exchange (ETDEWEB)
This report forms part of an independent review of the specifications, execution and results of Task 6 of the Aespoe Task Force on Modelling of Groundwater Flow and Transport of Solutes, which is seeking to provide a bridge between site characterization (SC) and performance assessment (PA) approaches to solute transport in fractured rock. The present report is concerned solely with Task 6C, which relates to the construction and parametrisation of a block-scale hydrostructural model of the TRUE Block Scale region of the Aespoe Hard Rock laboratory. The task objectives, specifications and outcome are summarised and reviewed. Also, consideration is given to how the hydrostructural model might affect the outcomes of Task 6D and 6E. The main conclusions of this review are summarised below: The Task 6C hydrostructural model is a more comprehensive approach to quantitatively describing a volume of fractured rock than has been achieved hitherto. The idea of including solute retention ...
2005-03-01
International Nuclear Information System (INIS)
Current methods to calculate dose distributions with organ motion can be broadly classified as 'dose convolution' and 'fluence convolution' methods. In the former, a static dose distribution is convolved with the probability distribution function (PDF) that characterizes the motion. However, artifacts are produced near the surface and around inhomogeneities because the method assumes shift invariance. Fluence convolution avoids these artifacts by convolving the PDF with the incident fluence instead of the patient dose. In this paper we present an alternative method that improves the accuracy, generality as well as the speed of dose calculation with organ motion. The algorithm starts by sampling an isocenter point from a parametrically defined space curve corresponding to the patient-specific motion trajectory. Then a photon is sampled in the linac head and propagated through the three-dimensional (3-D) collimator structure corresponding to a particular MLC segment chosen randomly from ...
2005-04-01
The Path to Fusion Energy for Concepts Currently at the Concept Exploration Level
Energy Technology Data Exchange (ETDEWEB)
Concept Exploration (CE) experiments within the Innovative Confinement Concept Program have a unique role which impacts their contributions to the development of fusion energy. As stated in the FESAC ''Report on Alternate Concepts:'' These [CE] programs are aimed at innovation and basic understanding of relevant scientific phenomena. The emphasis on innovation motivates their application to the search for a better fusion reactor configuration. In addition, because of their unique character the CE experiments offer excellent opportunities to couple fusion-plasma physics to other sciences. A recent example of coupling is the fusion self-organized plasmas to reconnection physics and extra-terrestrial plasmas. Perhaps of even greater importance is the education of the future scientists needed for developing fusion energy. The CE experiments, both at universities and national labs, are of a size students can ''get their hands ...
2003-01-09
Scale-up of two-phase flow in heterogeneous chalk. Matrix properties
Energy Technology Data Exchange (ETDEWEB)
This investigation presents scale-up of a detailed heterogeneous geostatistical model to a full field reservoir simulation model, considering both single and two-phase flow properties. The model represents a typical low permeability Danish North Sea chalk reservoir and includes capillary pressure and saturation end-point variations. Two new up-scaling methods has been investigated, all based on fine scale simulation on a cross section of the geomodel. The first methods assumes piston style behaviour and a coupled viscosity is introduced into the basic Darcy`s equations. The second method is a modification of the JBN method traditionally applied in analysing results from core flooding experiments, which emerged as the most successful and therefore also the recommended method. 1. In addition to the up scaling work we review the Equivalent Radius Method for capillary pressure normalisation with explicit derivation of type functions for Maastrichtian and Danian chalk ...
1998-02-01
Development of fire simulation models for radiative heat transfer and probabilistic risk assessment
Energy Technology Data Exchange (ETDEWEB)
An essential part of fire risk assessment is the analysis of fire hazards and fire propagation. In this work, models and tools for two different aspects of numerical fire simulation have been developed. The primary objectives have been firstly to investigate the possibility of exploiting state-of-the-art fire models within probabilistic fire risk assessments and secondly to develop a computationally efficient solver of thermal radiation for the Fire Dynamics Simulator (FDS) code. In the first part of the work, an engineering tool for probabilistic fire risk assessment has been developed. The tool can be used to perform Monte Carlo simulations of fires and is called the Probabilistic Fire Simulator (PFS). In Monte Carlo simulation, the simulations are repeated multiple times, covering the whole range of variability of the input parameters and thus resulting in a distribution of results covering what can be expected in reality. In practical applications, advanced simulation techniques ...
2008-07-01
Quantum information processing in nanostructures[Quantum optics; Quantum computing
Energy Technology Data Exchange (ETDEWEB)
Since information has been regarded os a physical entity, the field of quantum information theory has blossomed. This brings novel applications, such as quantum computation. This field has attracted the attention of numerous researchers with backgrounds ranging from computer science, mathematics and engineering, to the physical sciences. Thus, we now have an interdisciplinary field where great efforts are being made in order to build devices that should allow for the processing of information at a quantum level, and also in the understanding of the complex structure of some physical processes at a more basic level. This thesis is devoted to the theoretical study of structures at the nanometer-scale, 'nanostructures', through physical processes that mainly involve the solid-state and quantum optics, in order to propose reliable schemes for the processing of quantum information. Initially, the main results of quantum information theory and quantum computation are ...
2002-07-01
Contaminant transport in fracture networks with heterogeneous rock matrices. The Picnic code
Energy Technology Data Exchange (ETDEWEB)
In the context of safety assessment of radioactive waste repositories, complex radionuclide transport models covering key safety-relevant processes play a major role. In recent Swiss safety assessments, such as Kristallin-I, an important drawback was the limitation in geosphere modelling capability to account for geosphere heterogeneities. In marked contrast to this limitation in modelling capabilities, great effort has been put into investigating the heterogeneity of the geosphere as it impacts on hydrology. Structural geological methods have been used to look at the geometry of the flow paths on a small scale and the diffusion and sorption properties of different rock materials have been investigated. This huge amount of information could however be only partially applied in geosphere transport modelling. To make use of these investigations the 'PICNIC project' was established as a joint cooperation of PSI/Nagra and QuantiSci to provide a new geosphere transport ...
2001-02-01
Website Policies and Important Links Comments
WorldWideScience.org is maintained by the U.S. Department of Energy's
Office of Scientific and Technical Information as the Operating Agent
for the WorldWideScience Alliance.
