Science.gov (United States)

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...

Energy Technology Data Exchange (ETDEWEB)

The entire Virasoro, Ramond and Neveu-Schwarz algebras can each be constructed from a finite number of well-chosen generators satisfying a small number of conditions. Our most economical sets consist of just two starting generators in all cases, subject to no more than six conditions for the Virasoro case, five conditions for the Ramond case, and nine conditions for the Neveu-Schwarz case. Consequently, the Virasoro algebra simply amounts to 6 equations in two operator unknowns, and correspondingly 5 to 9 equations for the foregoing superalgebras. 2 refs.

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.

Science.gov (United States)

Test firings of 120-mm rounds were performed to characterize effects due to the combustible cartridge cases used. Comparisons were made between measured ballistic parameters, especially gas pressures, and those predicted by computerized simulation calcula...

CERN Document Server

We investigate the existence of time-periodic solutions of the Dirac equation in the Kerr-Newman background metric. To this end, the solutions are expanded in a Fourier series with respect to the time variable $t$ and the Chandrasekhar separation ansatz is applied so that the question of existence of a time-periodic solution is reduced to the solvability of a certain coupled system of ordinary differential equations. First, we prove the already known result that there are no time-periodic solutions in the non-extreme case. Then it is shown that in the extreme case for fixed black hole data there is a sequence of particle masses $(m_N)_{N\\in\\mathbb N}$ for which a time-periodic solution of the Dirac equation does exist. The period of the solution depends only on the data of the black hole described by the Kerr-Newman metric.

Energy Technology Data Exchange (ETDEWEB)

Many methods have been proposed for the determination of the corrosion current density from polarization curves; the two most commonly used ones are the three-point method and the curve-fitting method. The basis of these methods is a polarization equation that was derived for the complete absence of mass-transport effect or considers the cathodic partial reaction to be under complete mass-transport control. The mixed-control case has been much less investigated, even though this case may quite frequently reflect practical situations. A completely generalized polarization equation is proposed that can be used as a basis for curve-fitting data evaluation under any conditions. This equation was used to determine the error caused by the neglect of mass-transport effect in conventional data evaluation.

International Nuclear Information System (INIS)

We analyze the stability of the structure equations of the vacuum in the brane world models, by using both the linear (Lyapunov) stability analysis, and the Jacobi stability analysis, the Kosambi-Cartan-Chern theory. In the brane world models the four-dimensional effective Einstein equations acquire extra terms, called dark radiation and dark pressure, respectively, which arise from the embedding of the three-brane in the bulk. Generally, the spherically symmetric vacuum solutions of the brane gravitational field equations have properties quite distinct as compared to the standard black hole solutions of general relativity. We close the structure equations by assuming a simple linear equation of state for the dark pressure. In this case the vacuum is Jacobi stable only for a small range of values of the proportionality constant relating the dark pressure and the ...

International Nuclear Information System (INIS)

In terms of Berezins's theory of symbols of operators, the integral formulation is suggested for the free differential algebra which gives rise to consistent equations of motion of interacting massless fields of all spins 0#<=#s<#infinity# in the frameworks of gravity. In the first nontrivial order of the expansion in powers of curvatures, Frobenius consistency conditions for higher-spin equations of motion are shown to reduce to the simple geometrical fast that there are two ways for splitting any quadrangle in two triangles. To clarify our construction, we illustrate how it works in the simplest case of pure gravity. (orig.).

CERN Document Server

In symmetrical nuclear matter the solutions of pion dispersion equation are investigated in the complex plane of the pion frequency $\\omega$. There are three well-known branches of solutions on the physical sheet : sound, pion and isobar, at the matter density less than the critical one $\\rho <\\rho_c$. At the condition $\\omega^2_c\\leq0$ (in general case Re$ \\omega^2_c\\leq0$) takes place. This points out the instability of the ground state which is possibly related to the pion condensation.

CERN Document Server

We study perturbations of a Schwarzschild black hole in Chern-Simons modified gravity. We begin by showing that Birkhoff's theorem holds for a wide family of Chern-Simons coupling functions, a scalar field present in the theory that controls the strength of the Chern-Simons correction to the Einstein-Hilbert action. After decomposing the perturbations in spherical harmonics, we study the linearized modified field equations and find that axial and polar modes are coupled, in contrast to general relativity. The divergence of the modified equations leads to the Pontryagin constraint, which forces the vanishing of the Cunningham-Price-Moncrief master function associated with axial modes. We analyze the structure of these equations and find that the appearance of the Pontryagin constraint yields an overconstrained system that does not allow for generic black hole oscillations. We illustrate this situation by studying the ...

International Nuclear Information System (INIS)

Within the framework of the Bianchi type-IX homogeneous space, we set up a system of coupled equations for the cosmic scale factors, scalar field, and radiative energy density. At the tree level, the equations are written in a self-consistent, Hartree-Fock form. For phi"4 theory, the system of nine first-order differential equations is solved numerically for a varying ratio of the energy of anisotropy to the vacuum energy. As the vacuum energy increases, there appears to be less reheating, since the energy of anisotropy is more efficiently converted into isotropic expansion. If the energy of anisotropy is large enough, the inflationary phase is prevented. In this case, a series of cosmological phase transitions will take place each time the square of the effective mass changes its sign.

Energy Technology Data Exchange (ETDEWEB)

Turbulent premixed methane-air combustion in a cylindrical chamber is numerically simulated considering radiation effect. Reaction rates are considered as minimum rates between Arrhenius rates and eddy break up rates. A five step reduced mechanism is used. Turbulent modeling is done via standard k-{epsilon} model imposed by empirical inlet boundary conditions. Source terms of energy equation consist of reaction rates and radiation effects. The discrete ordinate method (DOM) is employed to solve the radiative transfer equation (RTE) and the weighted sum of gray gas model (WSGGM) is imposed to consider radiation effect of non-gray gases. The results indicate that in the case of turbulent combusting flows, the effect of radiation of gases can affect the temperature and species concentrations. The numerical results obtained considering radiation effect are closer to the experimental data than that of the ...

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 ...

International Nuclear Information System (INIS)

We present a general procedure to construct the first law of thermodynamics on the apparent horizon and illustrate its validity by examining it in some extended gravity theories. Applying this procedure, we can describe the thermodynamics on the apparent horizon in Randall-Sundrum braneworld imbedded in a nontrivial bulk. We discuss the mass-like function which was used to link Friedmann equation to the first law of thermodynamics and obtain its special case which gives the generalized Misner-Sharp mass in Lovelock gravity.

Energy Technology Data Exchange (ETDEWEB)

A mathematical model is proposed for three-phase filtering which is a particular case of a more general model of filtering multiple component mixtures and makes it possible to describe three-phase filtering flows under conditions where the pressure of saturation is variable in time and in space. In the proposed model, the system of unknown functions is selected so that it makes it possible to trace fairly simply the pressure of saturation and to compute the corresponding equation coefficients.

Energy Technology Data Exchange (ETDEWEB)

Singlet oxygen generators are multiphase flow chemical reactors used to generate energetic oxygen to be used as a fuel for chemical oxygen iodine lasers. In this paper, a theoretical model of the generator is presented along with its solutions over ranges of parameter space and oxygen maximizing optimizations. The singlet oxygen generator (SOG) is a low-pressure, multiphase flow chemical reactor that is used to produce molecular oxygen in an electronically excited state, i.e. singlet delta oxygen. The primary product of the reactor, the energetic oxygen, is used in a stage immediately succeeding the SOG to dissociate and energize iodine. The gas mixture including the iodine is accelerated to a supersonic speed and lased. Thus the SOG is the fuel generator for the chemical oxygen iodine laser (COIL). The COIL has important application for both military purposes--it was developed by the US Air Force in the 1970s--and, as the infrared beam is readily absorbed by metals, industrial cutting ...

International Nuclear Information System (INIS)

British Nuclear Fuels plc owns and operates the near-surface Drigg disposal facility for low level radioactive waste. The long-term performance of the site is modelled by a suite of computer codes called MONDRIAN. One of the modules of MONDRIAN deals with the transport of radionuclides through the environment, and this paper reports on the current status of this module (BIOS). We derive the basic set of working equations from first principles and show clearly how the approximate nature of the final equations is arrived at. This is done by an averaging process leading to compartments, in and out of which radionuclides, solids and water can flow. The equations allow radioactive decay chains and an arbitrary number of compartments. There is also the facility to deal with changes in the rate coefficients, thereby simulating different environmental states. It is also possible to include the creation of new compartments arising ...

Energy Technology Data Exchange (ETDEWEB)

This paper is an introduction course in modelling turbulent thermohydraulics, aimed at computational fluid dynamics users. No specific knowledge other than the Navier Stokes equations is required beforehand. Chapter I (which those who are not beginners can skip) provides basic ideas on turbulence physics and is taken up in a textbook prepared by the teaching team of the ENPC (Benque, Viollet). Chapter II describes turbulent viscosity type modelling and the 2k-{epsilon} two equations model. It provides details of the channel flow case and the boundary conditions. Chapter III describes the `standard` (R{sub ij}-{epsilon}) Reynolds tensions transport model and introduces more recent models called `feasible`. A second paper deals with heat transfer and the effects of gravity, and returns to the Reynolds stress transport model. (author). 37 refs.

International Nuclear Information System (INIS)

The stability of a couple stress fluid saturated horizontal porous layer heated from below and cooled from above when the fluid and solid phases are not in local thermal equilibrium is investigated. The Darcy model is used for the momentum equation and a two-field model is used for energy equation each representing the solid and fluid phases separately. The linear stability theory is employed to obtain the condition for the onset of convection. The effect of thermal non-equilibrium on the onset of convection is discussed. It is shown that the results of the thermal non-equilibrium Darcy model for the Newtonian fluid case can be recovered in the limit as couple stress parameter C?0. We also present asymptotic analysis for both small and large values of the inter phase heat transfer coefficient H. We found an excellent agreement between the exact solutions and asymptotic solutions when H is very small.

Energy Technology Data Exchange (ETDEWEB)

The performance of five well-known photosynthesis-light equations is compared by presenting a wide range of solutions in the form of dimensionless nomographs for the case where photosynthetically available radiation (PAR) reduces exponentially down the water column and is distributed sinusoidally through the photoperiod. These provide a simple means of calculating daily photosynthesis at any depth (e.g. at a benthic layer), or through a water column, avoiding the need to perform complex integrations. An examination is made of the accuracy of common approximate methods for calculating daily photosynthesis, assuming constant PAR. For optically deep water a modification is proposed to Talling's planimetric solution, to enable daily photosynthesis to be calculated more accurately, yet simply, over the whole range of possible PAR values. The errors induced by approximating the daily PAR distribution as a half-sinusoid are also reported. 13 ...

CERN Document Server

In the paper we investigate a mathematical model describing the growth of tumor in the presence of immune response of a host organism. The dynamics of tumour and immune cells is based on the generic Michaelis-Menten kinetics describing interaction and competition between the tumour and the immune system. The appropriate phenomenological equation modeling cell-mediated immune surveillance against cancer is of the predator-prey form and within a given choice of parameters exhibits bistability. Under the influence of spontaneous weak fluctuations, the model may be analyzed in terms of a stochastic differential equation bearing the form of an overdamped Langevin-like dynamics in the external quasi-potential represented by a double well. We analyze properties of the system within the range of parameters for which the potential wells are of the same depth and when the additional perturbation describing a periodic treatment is insufficient to overcome ...

CERN Document Server

We review and extend in several directions recent results on the asymptotic safety approach to quantum gravity. The central issue in this approach is the search of a Fixed Point having suitable properties, and the tool that is used is a type of Wilsonian renormalization group equation. We begin by discussing various cutoff schemes, i.e. ways of implementing the Wilsonian cutoff procedure. We compare the beta functions of the gravitational couplings obtained with different schemes, studying first the contribution of matter fields and then the so-called Einstein-Hilbert truncation, where only the cosmological constant and Newton's constant are retained. In this context we make connection with old results, in particular we reproduce the results of the epsilon expansion and the perturbative one loop divergences. We then apply the Renormalization Group to higher derivative gravity. In the case of a general action quadratic in curvature we recover, ...

CERN Document Server

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 ...

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 ...

CERN Document Server

The Cahn-Hilliard-Navier-Stokes system is based on a well-known diffuse interface model and describes the evolution of an incompressible isothermal mixture of binary fluids. A nonlocal variant consists of the Navier-Stokes equations suitably coupled with a nonlocal Cahn-Hilliard equation. The authors, jointly with P. Colli, have already proven the existence of a global weak solution to a nonlocal Cahn-Hilliard-Navier-Stokes system subject to no-slip and no-flux boundary conditions. Uniqueness is still an open issue even in dimension two. However, in this case, the energy identity holds. This property is exploited here to define, following J.M. Ball's approach, a generalized semiflow which has a global attractor. Through a similar argument, we can also show the existence of a (connected) global attractor for the convective nonlocal Cahn-Hilliard equation with a given velocity field, even in dimension ...

Energy Technology Data Exchange (ETDEWEB)

This thesis enters the field of highly compressed materials equation of state studies. In particular, it focuses on the case of laser shock compressed iron. This work indeed aims at getting to the conditions of the earth's core, comprising a solid inner core and a liquid outer core. The understanding of phenomena governing the core's thermodynamics and the geodynamic process requires the knowledge of iron melting line locus around the solid-liquid interface at 3.3 Mbar. Several experiments were performed to that extent. First, an absolute measurement of iron Hugoniot was obtained. Following is a study of partially released states of iron into a window material: lithium fluoride (LiF). This configuration enables direct access to compressed iron optical properties such as reflectivity and self-emission. Interface velocity measurement is dominated by compressed LiF optical properties and is used as a pressure gauge. Using a dual ...

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 ...

Energy Technology Data Exchange (ETDEWEB)

The high energy limit of Quantum Chromodynamics is one of the most fascinating areas in the theory of strong interactions. Over a decade ago the HERA experiment at DESY in Hamburg provided strong evidence for the rise of the proton structure function at small values of the Bjorken variable x. This behavior can be explained as an increase of the gluon density of the proton with energy or correspondingly with smaller values of x. This increase can be attributed on the other hand to the large probability of gluon splitting in QCD. The natural framework for describing the gluon dynamics at small x is the Balitskii-Fadin-Kuraev-Lipatov formalism developed some 30 years ago. It predicts that the gluon density grows very fast with increasing energy, as a power with a large intercept. This increase has to be tamed in order to satisfy the unitarily bound. Over two decades ago, Gribov, Levin and Ryskin proposed the mechanism called the parton saturation, which slows down the fast rise of the ...

Energy Technology Data Exchange (ETDEWEB)

Several polynomial finite elements of nodal type are introduced that should lead to convergence of O(h{sup 3}) in the L{sup 2} norm. Two of these methods are new and are expected to achieve the same orders of convergence with fewer parameters than the third method. They are applied to the one-group diffusion equation under different formulations, namely, several versions (with or without reduced and transverse integrations) of the primal and the mixed-hybrid formulations. Convergence rates are checked for a model problem with an analytical solution. Two of these methods exhibit superconvergence phenomena [O(h{sup 4}) instead of O(h{sup 3})], a fact that can be explained heuristically. The most promising method, with only five parameters per cell, turns out to yield only O(h{sup 2}) in its most algebraically efficient versions, while it has the potential of O(h{sup 3}) convergence rates. Again, an explanation is given for this behavior and a fully O(h{sup 3}) ...

International Nuclear Information System (INIS)

The onset of double-diffusive convection in a couple-stress fluid-saturated horizontal porous layer is studied using linear and weak nonlinear stability analyses. The modified Darcy equation that includes the time derivative term and the inertia term is used to model the momentum equation. The expressions for stationary, oscillatory and finite-amplitude Rayleigh number are obtained as a function of the governing parameters. The effect of couple-stress parameter, solute Rayleigh number, Vadasz number and diffusivity ratio on stationary, oscillatory and finite-amplitude convection is shown graphically. It is found that the couple-stress parameter and the solute Rayleigh number have a stabilizing effect on stationary, oscillatory and finite-amplitude convection. The diffusivity ratio has a destabilizing effect in the case of stationary and finite-amplitude modes, with a dual effect in the case of ...

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.''

CERN Document Server

It has been proposed that a hidden conformal field theory (CFT) governs the dynamics of low frequency scattering in a general Kerr black hole background. We further investigate this correspondence by mapping higher order corrections to the massless wave equations in a Kerr background to an expansion within the CFT in terms of higher dimension operators. This implies the presence of infinite towers of CFT primary operators with positive conformal dimensions compatible with unitarity. The exact Kerr background softly breaks the conformal symmetry and the scaling dimensions of these operators run with frequency. The scale-invariant fixed point is dual to a degenerate case of flat spacetime.

Energy Technology Data Exchange (ETDEWEB)

The diffusion coefficient of particles in radioactive gases has been calculated with account of random wandering of aerosols (occurrence of local fields affecting the particles; recoils accompanying radiation emitted by particles, etc.). To determine the diffusion coefficient, the method of Fokker-Planck equation derivation was used. A formula is presented for calculating the radiation-stimulated diffusion coefficient. A linear growth of the diffusion coefficient with radioactivity is noted according to the formula, the diffusion coefficient is mainly determined by the field in the radiation damage region. The aerosol radioactivity may result in a more rapid deposition of aerosols in the pipelines and aerosol purification systems. The diffusion rate grows not only in the presence of intrinsic radioactivity but in case of external radiation exposure as well.

International Nuclear Information System (INIS)

Open descendants extend conformal field theory to unoriented surfaces with boundaries. The construction rests on two types of generalizations of the fusion algebra. The first is needed even in the relatively simple case of diagonal models. It leads to a new tensor that satisfies the fusion algebra, but whose entries are signed integers. The second is needed when dealing with non-diagonal models, where Cardy's ansatz does not apply. It leads to a new tensor with positive integer entries, that satisfies a set of polynomial equations and encodes the classification of the allowed boundary operators. (orig.).

Energy Technology Data Exchange (ETDEWEB)

The characteristics of the temperature dependence of the upper critical magnetic field, Hc2(T), of structurally inhomogeneous superconductors (e.g., ternary molybdenum chalcogenides, A-15 compounds, transition metal alloys, and amorphous films) are investigated analytically. The MWGH equation for Hc2(T) is generalized to the case of weakly inhomogeneous systems with a characteristic inhomogeneity scale much smaller than the effective coherence length. It is shown that an increase in the dispersion of the diffusion coefficient leads to an increase in the slope and width of the linear section of the Hc2(T) curve. 11 references.

CERN Document Server

We describe dimensionally constrained symbolic regression which has been developed for mass measurement in certain classes of events in high-energy physics (HEP). With symbolic regression, we can derive equations that are well known in HEP. However, in problems with large number of variables, we find that by constraining the terms allowed in the symbolic regression, convergence behavior is improved. Dimensionally constrained symbolic regression (DCSR) finds solutions with much better fitness than is normally possible with symbolic regression. In some cases, novel solutions are found.

British Library Electronic Table of Contents (United Kingdom)

In this paper we present a new, accurate form of the heat balance integral method, termed the combined integral method (CIM). The application of this method to Stefan problems is discussed. For simple test cases the results are compared with exact and asymptotic limits. In particular, it is shown that the CIM is more accurate than the second order, large Stefan number, perturbation solution for a wide range of Stefan numbers. In the initial examples it is shown that the CIM reduces the standard problem, consisting of a PDE defined over a domain specified by an ODE, to the solution of one or two algebraic equations. The latter examples, where the boundary temperature varies with time, reduce to a set of three first order ODEs.

Energy Technology Data Exchange (ETDEWEB)

The influence of a periodic electromagnetic field on the radiative corrections to the atomic energy levels is studied for the case of ''strong fields'' for which the interaction between the atom and field is of the order or greater than the radiative effects. The analysis is carried out on the basis of the Schwinger-Dirac equation for the propagation function of a bound electron in the field and on the basis of the density matrix in the Furry representation. It is shown that in the strong field approximation the radiative shifts and widths are manifest as radiative corrections to the quasi-energies. In super-high resolution experiments intensity effects in the radiative corrections to the atomic levels are obtained in the case of single-photon resonance. Some multiphoton processes are condidered by taking into account the effect of the field on the radiative structure of the levels.

CERN Document Server

We review a new theoretical approach to the kinetics of surfactant adsorption at fluid-fluid interfaces. It yields a more complete description of the kinetics both in the aqueous solution and at the interface, deriving all equations from a free-energy functional. It also provides a general method to calculate dynamic surface tensions. For non-ionic surfactants the results coincide with previous models. Common non-ionic surfactants are shown to undergo diffusion-limited adsorption, in agreement with experiments. Strong electrostatic interactions in salt-free ionic surfactant solutions are found to lead to kinetically limited adsorption. In this case the theory accounts for unusual experimental results which could not be understood using previous approaches. Added salt screens the electrostatic interactions and makes the ionic surfactant adsorption similar to the non-ionic case. The departure from the non-ionic behavior as ...

Energy Technology Data Exchange (ETDEWEB)

There are various reasons for preferring local (e.g., three bump) orbit correction methods to global corrections. One is the difficulty of solving the mN equations for the required mN correcting bumps, where N is the number of superperiods and m is the number of bumps per superperiod. The latter is not a valid reason for avoiding global corrections, since, we can take advantage of the superperiod symmetry to reduce the mN simultaneous equations to N separate problems, each involving only m simultaneous equations. Previously, I have shown how to solve the general problem when the machine contains unknown magnet errors of known probability distribution; we made measurements of known precision of the orbit displacements at a set of points, and we wish to apply correcting bumps to minimize the weighted rms orbit deviations. In this report, we will consider two simpler problems, using similar methods. We consider the ...

Energy Technology Data Exchange (ETDEWEB)

In this paper, the wavefunction of the universe with a tunneling boundary condition is considered in the context of the Brans-Dicke-type scalar-tensor theory with matter. The matter may be interpreted as a D-particle (or D0-brane) in string theory when the Brans-Dicke parameter {omega} is -1. We study two simple examples. The first example, the {gamma}=0 (matter) case, has a scale factor duality even if the low energy string action is coupled to matter. The universe undergoes quantum transition from super-inflationary (pre-big-bang) to deflationary (post-big-bang) phase. We calculate the transition rate by solving the Wheeler-DeWitt equation and find that it is non-vanishing. The two phases are disconnected classically. The second example is the {gamma}=1/3(radiation) case. With the help of earlier work this matter can be identified with a D0-brane in string theory. In this case, due to the absence of ...

Energy Technology Data Exchange (ETDEWEB)

The stability of a couple stress fluid saturated horizontal porous layer heated from below and cooled from above when the fluid and solid phases are not in local thermal equilibrium is investigated. The Darcy model is used for the momentum equation and a two-field model is used for energy equation each representing the solid and fluid phases separately. The linear stability theory is employed to obtain the condition for the onset of convection. The effect of thermal non-equilibrium on the onset of convection is discussed. It is shown that the results of the thermal non-equilibrium Darcy model for the Newtonian fluid case can be recovered in the limit as couple stress parameter C{yields}0. We also present asymptotic analysis for both small and large values of the inter phase heat transfer coefficient H. We found an excellent agreement between the exact solutions and asymptotic solutions when H is very small.

International Nuclear Information System (INIS)

To systematically analyze the dynamical implications of the matter content in cosmology, we generalize earlier dynamical systems approaches so that perfect fluids with a general barotropic equation of state can be treated. We focus on locally rotationally symmetric Bianchi type IX and Kantowski-Sachs orthogonal perfect fluid models, since such models exhibit a particularly rich dynamical structure and also illustrate typical features of more general cases. For these models, we recast Einstein's field equations into a regular system on a compact state space, which is the basis for our analysis. We prove that models expand from a singularity and recollapse to a singularity when the perfect fluid satisfies the strong energy condition. When the matter source admits Einstein's static model, we present a comprehensive dynamical description, which includes the qualitative asymptotic behavior, of models in the neighborhood of the ...

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 ...

International Nuclear Information System (INIS)

Analyses of the date obtained from the instrumented impact test and of its correlation with static and dynamic J integral values and various tensile properties are made in this study. It is shown that a dynamic bending yield load recorded in the instrumented Charpy V test can be easily converted to a dynamic tensile yield stress using Server's equation. It is also possible to estimate a true fracture strain, COD value, and SZW easily if the load point displacement is recorded. Temperature dependency of Charpy V energy and J integral values is rather well simulated by a curve fitting method. However, a half value transition temperture obtained in that method has no physical meaning. Dynamic J sub(ID) is larger than a static J sub(IC) in the ductile fracture range in this study. This rise is proportional to the yield stress ratio #sigma#sub(yd)/#sigma#sub(y)(--1.2-1.3). It is found, therefore, that the following relation exists in the both cases; ...

International Nuclear Information System (INIS)

The present paper is concerned with development and application of a so-called Effective Convection Model (ECM), which aims to provide a detailed, mechanistic description of heat transfer processes in a BWR lower plenum. The ECM is a Computational Fluid Dynamics (CFD)-like tool which employs a simpler and more effective approach to compute heat transfer by solving only energy conservation equation instead of solving the full set of Navier-Stokes and energy equations by a CFD code. We implement the ECM in a CFD code (Fluent), with detailed description of the ECM development, implementation and validation. A dual approach is used to validate the ECM, namely validation against experimental data and against heat transfer results obtained by CFD predictions in the same geometries and conditions. Insights gained from CFD simulations are also used to improve ECM. The ECM capability as an effective tool to simulate heat transfer of an internally heated ...

International Nuclear Information System (INIS)

High-energy ionizing radiation is a prominent modality for the treatment of many cancers. The approaches to electron dose calculation can be categorized into semi-empirical models (e.g. Fermi-Eyges, convolution-superposition) and probabilistic methods (e.g. Monte Carlo). A third approach to dose calculation has only recently attracted attention in the medical physics community. This approach is based on the deterministic kinetic equations of radiative transfer. We derive a macroscopic partial differential equation model for electron transport in tissue. This model involves an angular closure in the phase space. It is exact for the free streaming and the isotropic regime. We solve it numerically by a newly developed HLLC scheme based on Berthon et al (2007 J. Sci. Comput. 31 347-89) that exactly preserves the key properties of the analytical solution on the discrete level. We discuss several test cases taken from the medical ...

CERN Document Server

We study the dynamics of states perturbatively expanded about a harmonic system of loop quantum cosmology, exhibiting a bounce. In particular, the evolution equations for the first and second order moments of the system are analyzed. These moments back-react on the trajectories of the expectation values of the state and hence alter the energy density at the bounce. This analysis is performed for isotropic loop quantum cosmology coupled to a scalar field with a small but non-zero constant potential, hence in a regime in which the kinetic energy of matter dominates. Analytic restrictions on the existence of dynamical coherent states and the meaning of semi-classicality within these systems are discussed. A numerical investigation of the trajectories of states that remain semi-classical across the bounce demonstrates that, at least for such states, the bounce persists and that its properties are similar to the standard case, in which the moments ...

Energy Technology Data Exchange (ETDEWEB)

In this article, the ability of artificial neural networks in prediction of separation in steady two dimensional boundary layer flows is studied. Data for network training is extracted from numerical solution of an ODE obtained from Von Karman integral equation with approximate one parameter Pohlhousen velocity profile. As an appropriate neural network, a two layer radial basis generalized regression artificial neural network is used. The results shows good agreements between the overall behavior of the flow fields predicted by the artificial neural network and the actual flow fields for some cases. The method easily can be extended to unsteady separation and turbulent as well as compressible boundary layer flows. (author)

CERN Document Server

In this paper, we discuss the large--time behavior of solution of a simple kinetic model of Boltzmann--Maxwell type, such that the temperature is time decreasing and/or time increasing. We show that, under the combined effects of the nonlinearity and of the time--monotonicity of the temperature, the kinetic model has non trivial quasi-stationary states with power law tails. In order to do this we consider a suitable asymptotic limit of the model yielding a Fokker-Planck equation for the distribution. The same idea is applied to investigate the large-time behavior of an elementary kinetic model of economy involving both exchanges between agents and increasing and/or decreasing of the mean wealth. In this last case, the large-time behavior of the solution shows a Pareto power law tail. Numerical results confirm the previous analysis.

CERN Document Server

Within a fully relativistic framework, we derive and solve numerically the perturbation equations of relativistic stars, including the stresses produced by a non-vanishing shear viscosity in the stress-energy tensor. With this approach, the real and imaginary parts of the frequency of the modes are consistently obtained. We find that, approaching the inviscid limit from the finite viscosity case, the continuous spectrum is regularized and we can calculate the quasi-normal modes for stellar models that do not admit solutions at first order in perturbation theory when the coupling between the polar and axial perturbations is neglected. The viscous damping time is found to agree within factor 2 with the usual estimate obtained by using the eigenfunctions of the inviscid limit and some approximation for the energy dissipation integrals. We find that the frequencies and viscous damping times for relativistic $r-$modes lie between the Newtonian and ...

British Library Electronic Table of Contents (United Kingdom)

We have discovered a previously unidentified thermal explosion mode using the Los Alamos deflagration cylinder test (DFCT). The DFCT is a "pipe bomb"-style test similar to the detonation cylinder test (DTCT), which has been used for many years to calibrate detonation product equations of state. The shot is heated in an oven to a uniform test temperature. The pre-heated high explosive (HE) is triggered by a hot wire initiator on one end. The tube is back-illuminated by a bright light source, and its combustion-driven deformation and subsequent break-up are observed by a high-speed framing camera. Like the DTCT, the DFCT tube wall motion provides the primary diagnostic. A variety of reactive responses are possible, including quasi-steady deflagration and deflagration-to-detonation transition...

British Library Electronic Table of Contents (United Kingdom)

There are two difficulties encountered in modeling valveless micropumps using lumped-element methods. The pressure loss coefficient for fluidic diodes used in valveless pumps to rectify flow depends on the flow direction. A problem arises in choosing the proper loss correlation because the flow direction is not known a priori. Another problem is the quadratic form of the equation for the flow through the fluidic diodes, which brings about multiple solutions. The above problems become even more serious in multi-chamber cases. They are overcome in this study by suitably formulating the flow resistance. In addition, the flow inertia is accounted for in the unsteady model. The steady and unsteady models are evaluated by comparing with CFD simulations, which also serve to illustrate the flow fi...

Energy Technology Data Exchange (ETDEWEB)

The paper presents an alternative approach for the studies of high voltage transmission lines based on artificial intelligence and more specifically artificial neural networks (ANNs). In contrast to the existing conventional-analytical techniques and simulations which are using in the calculations empirical and/or approximating equations, this approach is based only on actual field data and actual measurements. The proposed approach is applied on high voltage transmission lines in order to calculate the lightning outages, on grounding systems in order to assess the grounding resistance and on high voltage transmission lines' polluted insulators in order to estimate the critical flashover voltage. The obtained results are very close to the actual ones for all three case studies, something which clearly implies that the ANN approach is well working and has an acceptable accuracy, constituting an additional tool of electric engineers. ...

Energy Technology Data Exchange (ETDEWEB)

Whereas in part 1 and 2 the detailed observation equations of terrestrial type for an integrated geodetic model are presented this paper will outline how the theory is transferred in an operational program. The main lines of functions of the FORTRAN IV program are discussed. Since the input specifications as well as the data requirements are described in detail and illustrated by example of input and output part 3 can serve as an user-guide for OPERA. OPERA is an acronym for operational adjustment. Besides the integrated determination of 3D-geocentric coordinates and the gravity disturbing potential the program can handle all cases of traditional geodesy (in total 13 variants of solution are provided) equivalent also with adjustments in ellipsoidal coordinates B, L or H, as well as constrained adjustments, pure prediction of the gravity disturbing field, etc.

British Library Electronic Table of Contents (United Kingdom)

Forced harmonic vibration of a non-uniform elastic beam with attached dynamic vibration absorbers (DVA) is studied. Analytical approximation of the solution is obtained by the functional perturbation method (FPM). The problem has application to cutting tools operations where the resistance of the tool holder against regenerative chatter can be enhanced by optimizing the real part of the frequency response function (FRF). A test case of a beam with step-like heterogeneity and single DVA at the tip shows that the FPM solution is very accurate for up to ?40 percent deviation in both stiffness and mass density. Using the analytical results and Sims approach, optimal DVA tuning is found for each set of beam heterogeneity parameters by solving a set of nonlinear algebraic equations numer...

International Nuclear Information System (INIS)

A numerical analysis is performed to clarify the heat transfer characteristics of combined convection around the tubes which form a single row perpendicularly arranged to a forced flow in the direction of gravity. A boundary-fitted coordinate transformation technique is adopted to solve the governing equations numerically. It is found in the case of liquid sodium that gravity-induced flow gives no substantial effect to heat transfer augmentation in a low Reynolds number region while it can be sufficiently expected in the ordinary fluids in the same Reynolds number region. This means that such heat transfer augmentation as expected in ordinary fluids can not be realized when crossflow-type heat exchangers installed for decay-heat removal operates in a low Reynolds number region. (author).

International Nuclear Information System (INIS)

A one-dimensional tachyon Klein-Gordon equation is reduced to a nonrelativistic-tachyon equation of motion. The interpretation of this reduced equation leads to the following conclusions: 1) tachyons can be localized in time instead of in space as compared with bradyons, 2) space representation and momentum representation of bradyonic quantum equation of motion are replaced by time representation and energy representation in tachyon quantum equation of motion and 3) with the aid of these results, it has been found that the solutions of the tachyon Klein-Gordon equation of motion form a complete set. (author).

International Nuclear Information System (INIS)

The linearized-augmented-plane-wave (LAPW) method for thin films is generalized by removing the remaining shape approximation to the potential inside the atomic spheres. A new technique for solving Poisson's equation for a general charge density and potential is described and implemented in the film LAPW method. In the resulting full-potential LAPW method (FLAPW), all contributions to the potential are completely taken into account in the Hamiltonian matrix elements. The accuracy of the method: already well known for clean metal surfaces: is demonstrated for the case of a nearly free (noninteracting) O_2 molecule which is a severe test case of the method because of its large anisotropic charge distribution. Detailed comparisons show that the accuracy of the FLAPW results for O_2 exceeds that of existing state-of-the-art local-density linear-combination-of-atomic-orbitals (LCAO)-type calculations, and that taking the full ...

International Nuclear Information System (INIS)

We study the magnetic braking and viscous damping of differential rotation in incompressible, uniform density stars in general relativity. Differentially rotating stars can support significantly more mass in equilibrium than nonrotating or uniformly rotating stars, according to general relativity. The remnant of a binary neutron star merger or supernova core collapse may produce such a 'hypermassive' neutron star. Although a hypermassive neutron star may be stable on a dynamical time scale, magnetic braking and viscous damping of differential rotation will ultimately alter the equilibrium structure, possibly leading to delayed catastrophic collapse. Here we treat the slow-rotation, weak-magnetic field limit in which E_r_o_t<

Energy Technology Data Exchange (ETDEWEB)

The Kadomtsev-Petviashvili equation describes nonlinear dispersive waves which travel mainly in one direction, generalizing the Korteweg-de Vries equation for purely uni-directional waves. In this Letter we derive an improved KP-equation that has exact dispersion in the main propagation direction and that is accurate in second order of the wave height. Moreover, different from the KP-equation, this new equation is also valid for waves on deep water. These properties are inherited from the AB-equation (E. van Groesen, Andonowati, 2007 ) which is the unidirectional improvement of the KdV equation. The derivation of the equation uses the variational formulation of surface water waves, and inherits the basic Hamiltonian structure.

International Nuclear Information System (INIS)

In the procedure of the hydraulic control system analysis, a linearized approximate equation described by the first order term of Taylor's series has been widely used. Such a linearized equation is effective just near the operating point. In this study, the authors estimate computational errors in the process of applying the existing linearized equation stated above. For evaluating the computational accuracy in practical applications of the linearized equations, dynamic behaviors of hydraulic control systems are investigated through simulations with several kinds of representative hydraulic systems and the linearized equations suggested in this study.

Science.gov (United States)

A method for calculating flow through an airfoil cascade drawn on a surface of revolution is discussed. The three-dimensional flow was assumed to be represented by part-channels of varying width. The basic equations are the equation of continuity and the ...

Science.gov (United States)

... 4 II. INERTIAL NAVIGATION ERROR EQUATIONS ..... 5 ... 4- Page 20. Chapter II INERTIAL NAVIGATION ERROR EQUATIONS ...

Science.gov (United States)

A rigorous solution for the spectrum of a quasioptical cylindrical cavity resonator with a randomly rough side boundary has been obtained. To accomplish this task, we have developed a method for the separation of variables in a wave equation, which enables one, in principle, to rigorously examine any limiting case-from negligibly weak to arbitrarily strong disorder at the resonator boundary. It is shown that the effect of disorder-induced scattering can be properly described in terms of two geometric potentials, specifically, the "amplitude" and the "gradient" potentials, which appear in wave equations in the course of conformal smoothing of the resonator boundaries. The scattering resulting from the gradient potential appears to be dominant, and its impact on the whole spectrum is governed by the unique sharpness parameter ?, the mean tangent of the asperity slope. As opposed to the resonator with bulk disorder, the ...

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 ...

International Nuclear Information System (INIS)

It was shown by Gibbons and Tsarev (1996 Phys. Lett. A 211 19, 1999 Phys. Lett. A 258 263) that N-parameter reductions of the Benney equations correspond to particular N-parameter families of conformal maps. In recent papers (Baldwin and Gibbons 2003 J. Phys. A: Math. Gen. 36 8393-417, Baldwin and Gibbons 2004 J. Phys. A: Math. Gen. 37 5341-54), the present authors have constructed examples of such reductions where the mappings take the upper half p-plane to a polygonal slit domain in the #lambda#-plane. In those cases, the mapping function was expressed in terms of the derivatives of Kleinian #sigma# functions of hyperelliptic curves, restricted to the one-dimensional stratum #THETA#_1 of the #THETA#-divisor. This was done using an extension of the method given in Enolskii et al (2003 J. Nonlinear Sci. 13 157) extended to a genus 3 curve (Enolski V Z and Gibbons J Addition theorems on the strata of the theta divisor of genus three ...

Energy Technology Data Exchange (ETDEWEB)

The boundary-integral equation method is well suited for the calculation of the dynamic-stiffness matrix of foundations embedded in a layered visco-elastic halfspace (or a transmitting boundary of arbitrary shape), which represents an unbounded domain. It also allows pile groups to be analyzed, taking pile-soil-pile interaction into account. The discretization of this boundary-element method is restricted to the structure-soil interface. All trial functions satisfy exactly the field equations and the radiation condition at infinity. In the indirect boundary-element method distributed source loads of initially unknown intensities act on a source line located in the excavated part of the soil and are determined such that the prescribed boundary conditions on the structure-soil interface are satisfied in an average sense. In the two-dimensional case the variables are expanded in a Fourier integral in the wave number domain, ...

Energy Technology Data Exchange (ETDEWEB)

The temporal evolution of linear toroidal ion temperature gradient (ITG) modes is studied based on a kinetic integral equation including an initial condition. It is shown how to evaluate the analytic continuation of the integral kernel as a function of a complex-valued frequency, which is useful for analytical and numerical calculations of the asymptotic damping behavior of the ITG mode. In the presence of the toroidal {nabla}B-curvature drift, the temporal dependence of the density and potential perturbations consists of normal modes and a continuum mode, which correspond to contributions from poles and from an integral along a branch cut, respectively, of the Laplace-transformed potential function of the complex-valued frequency. The normal modes have exponential time dependence with frequencies and growth rates determined by the dispersion relation while the continuum mode, which has a ballooning structure, shows a power law decay {proportional_to} t{sup -2} in ...

Energy Technology Data Exchange (ETDEWEB)

The purpose of this study is to investigate the feasibility of visualization and void fraction measurement of air-water two-phase flow in a small diameter tube (inner diameter; 4.08mm) by using the real-time neutron radiography and image processing techniques. Video images of two-phase flow were taken by using the real-time neutron radiography system (thermal neutron radiography facility No.2) installed at the Japan Research Reactor 3 M of the Japan Atomic Energy Research Institute. The shape of bubbles and its moving behavior were clearly observed from the video images. The image corrections for dark current, shading, field intensity fluctuation and electrical system drift were examined in order to measure the void fraction from the video images. Though, generally speaking, the effect of the scattered neutron could not be ignored for quantification of the images taken by the neutron radiography, the scattered neutron could not affect the final results of void fraction in the ...

CERN Document Server

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.

Science.gov (United States)

... slenderness ratio. Equations which determine the ultimate compressive strength of a composite were established. Several ...

KoreaScience (Korea)

Full Text Available

British Library Electronic Table of Contents (United Kingdom)

A method for transforming the Euler and Navier-Stokes equations and a complete system of fluid dynamics equations in three dimensions to a closed system on any moving surface is proposed. As a result, for an arbitrary geometric configuration, the dimension of the equations is reduced by one, which makes them convenient for numerical simulation. The general principles of the method are described, and verifying examples are presented.

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).

International Nuclear Information System (INIS)

The report addresses the application of reverse osmosis to the treatment of radioactive liquid waste, focusing on the effects of non-radioactive cations and anions, existing in radioactive waste, on the removal of radioactive nuclides. The bivalent cations are generally about ten times as high in decontamination factor as monovalent ones. In the case of anions, the factor of SO_4"2"- is about ten times that of Cl"- while that of NO_3"- is somewhat smaller than that of Cl"-. During permeation through membrane, anions tend to accompany cations, and the decontamination factor of accompanying ions is found to increase with that of the paring ions. These results indicate that if one can select components that will remain in the liquid waste, one should consider their effects on the efficiency of its treatment. Furthermore, the interaction between two anions or between two cations can be approximated by a linear equation consisting of the logarithm ...

International Nuclear Information System (INIS)

The Flory model is applied to predict the isobaric expansion coefficients, ?P, isentropic, ?S, and isothermal, ?T, compressibilities and speeds of sound, u, of the highly complex mixtures: hydroxyether + alkane, + dibutylether, + 1-butanol or + 2-methoxyethanol, 1-alkanol + alkane, and 1-alkanol + dibutylether. Predictions were obtained using the energetic parameter, ?12, determined from values of excess enthalpies, HE, and from values of the excess isochoric internal energies, UVE, at equimolar composition. No meaningful differences exist between such predictions. Deviations between experimental and calculated values are lower than 2% for mixtures containing alkoxyethanols. Poorer results are obtained for 1-alkanol + dibutylether mixtures with deviations up to 5%. It is shown that predictions on ?P, ?S, ?T, and u essentially depend on structural effects and not on the orientational effects present in the studied mixtures. Results are improved using ?12 values fitted to molar excess ...

International Nuclear Information System (INIS)

In order to study the boundary conditions for the occurrence of ratcheting, a testing rig was erected. In this rig, small-scale mock-ups of the first wall were tested under mechanical conditions similar to those expected during disruptions. These mock-ups were made of stainless steel AISI 316. They were heated on the surface and water cooled from inside. Disruption forces were simulated by a high-speed hydraulic bending system by impact loads. The mock-ups were heated by a high-frequency generator and cooled by water through two cooling channels.In addition to the experimental work, simplified analytical analysis and numerical calculations were carried out using the finite element program ABAQUS. In these calculations, the temperature and stress fields in the mock-ups were investigated and the influence of the different constitutive equations was studied.The experimental results confirm the ratcheting behaviour predicted by the calculations. It has been ...

Energy Technology Data Exchange (ETDEWEB)

In order to design roof bolting systems safely and economically, it is essential to understand the flexural behaviour of the immediate roof. Based on the strata sequence, the strata in the immediate roof are divided into three types. The flexural behaviour of the three strata types are investigated in terms of the following effects: roof span, horizontal stress, thickness and Young's modulus of the lowest strata. The suspension reinforcement mechanism is analysed using beam-column theory. The equations for the maximum bending stress, deflection and transferred bolt load for the bolted strata are derived. In the analysis, the bolt load is assumed to be a point load and a horizontal stress is uniformly applied to each stratum. The friction reinforcement mechanism is also investigated. The major function of roof bolting in this case is to create frictional resistance by tensioning the roof bolts so that the individual layers are combined ...

Science.gov (United States)

In response to the consultation from the National Council on Radiation Protection and Measurements (NCRP) to the International Commission on Radiological Protection (ICRP) committee 2, retention and excretion of 239Pu deposited at wound site were calculated by coupling together the proposed NCRP wound model and the current ICRP systemic model of Pu. The physicochemical forms considered were the soluble form categorized into 'Strong Retention', and the colloidal, particulate and fragmentary forms. The results are summarized as follows. If in soluble form, immediate medical intervention is needed to prevent uptake of radionuclides to body tissues, and prompt wound monitoring is essential for an accurate estimation of the initially deposited radioactivity. If in particulate form, a multi-component exponential equation leads to an overestimation of the absorption rate to blood because of significant lymph node drainage. The committed doses in the organs for direct ...

International Nuclear Information System (INIS)

The electromagnetic field due to an air explosion is here studied: neutron and the gamma-rays, generated by neutronic reactions, deposit energy which ionizes atmosphere; recoil electrons are also created by gamma collisions. This data set allows to solve the Maxwell equations which manage the electromagnetic field. The TRIPOLI-2 code studies the coupled neutron-gamma transport in 3D- geometries by the Monte Carlo method. The code has been modified to calculate the photon energy deposited in matter and the recoil electron current created by Compton effect. The method is tested wiht a simple case; then neutron and gamma transport is studied in air kerma, deposited photon energy, electron current are calculated as functions of space and time and the contributions of the different neutronic reactions are separately evaluated. The calculations presented here are only part of studies about this subject. Developments will relate three dimensional ...

Energy Technology Data Exchange (ETDEWEB)

Nowadays, the floating production, storage, and offloading (FPSO) units are used to marginal oil fields and deep seas, because of mobility and reduced lead time from discovery to first production. One of the methods for offloading is to moor a tanker alongside the FPSO unit. In this case motion response of the system becomes highly complex, and will be affected by many factors such as hydrodynamic interaction and mechanical connections between the vessels. In this study, a quite general method is developed which has taken into account the hydrodynamic interaction as well as the effects of connectors and mooring lines, in motion response analysis of a multi-body floating system. For this purpose, 3-D source distribution method is applied for hydrodynamic analysis, and linear stiffness matrices are introduced to represent the effect of connectors and mooring lines on the motion equations. The method has been used for motion analysis of a ...

Energy Technology Data Exchange (ETDEWEB)

A flexible and computationally economical model of the internal combustion engine was developed for use on large digital computer systems. It is based on a system of ordinary differential equations for cylinder-averaged properties. The computer program is capable of multicycle calculations, with some parameters varying from cycle to cycle, and has restart capabilities. It can accommodate a broad spectrum of reactants, permits changes in physical properties, and offers a wide selection of alternative modeling functions without any reprogramming. It readily adapts to the amount of information available in a particular case because the model is in fact a hierarchy of five models. The models range from a simple model requiring only thermodynamic properties to a complex model demanding full combustion kinetics, transport properties, and poppet valve flow characteristics. Among its many features the model includes heat transfer, valve timing, ...

CERN Document Server

A risk-neutral method is always used to price and hedge contingent claims in complete market, but another method based on utility maximization or risk minimization is wildly used in more general case. One can find all kinds of special risk measure in literature. In this paper, instead of using market modified risk measure, we use a kind of risk measure induced by g_\\Gamma-solution or the minimal solution of a Constrained Backward Stochastic Differential Equation (CBSDE) directly when constraints on wealth and portfolio process comes to our consideration. Such g_\\Gamma-solution and the risk measure generated by it is well defined on appropriate space under suitable conditions. We adopt the inf-convolution of convex risk measures to solve some optimization problem. A dynamic version risk measures defined through g_\\Gamma-solution and some similar results about optimal problem can be got in our new framework and by our new approach.

Energy Technology Data Exchange (ETDEWEB)

An experimental study was conducted to investigate the heat transfer and friction loss characteristics in a horizontal rectangular channel having attachments of hollow rectangular profile fins over one of its heated surface. The Reynolds number based on the flow averaged inlet velocity and the hydraulic diameter, ranged from 3000 to 32,000. The hollow rectangular profile fins in 10cm height and axb=2cmx4cm dimensions with a thickness of 0.2cm were mounted on a heating surface vertically. Reynolds number, fin arrangement and fin pitch in the flow direction were the experimental parameters. Both in-line and staggered fin arrangements were studied for one-fixed spanwise (S{sub x}/a=3) and four different streamwise (S{sub y}/b=1.5, 1.875, 2.5 and 3.75) distances. Correlation equations for Nu, f and thermal performances were determined for fin configurations and the straight channel case without fins. (author)

CERN Document Server

Complete factorization scheme analysis of $F_2^{\\gamma}(x,Q^2)$ in the next--to--leading order QCD is performed. It is shown that the presence of the inhomogenous term in the evolution equations for parton distribution functions of the photon implies qualitative change in the way factorization mechanism works in photon--hadron and photon--photon collisions as compared to hadron--hadron ones. The incompletness of the existing NLO analyses of $F_2^{\\gamma}(x,Q^2)$ is pointed out and its origins traced back to the misinterpretation of the behaviour of $q^{\\gamma}(x,M)$ as a function of in the LO photonic coefficient function $C_{\\gamma}^{(0)}(x)$ is given and smooth transition between the properties of virtual and real photon discussed. A simple way of avoiding the problems with the term $\\ln(1-x)$ appearing in the standard expression for $C_{\\gamma}^{(0)}(x)$ in the case of the reaal photon is proposed.

Science.gov (United States)

Metabolic carbon labelling experiments enable a large amount of extracellular fluxes and intracellular carbon isotope enrichments to be measured. Since the relation between the measured quantities and the unknown intracellular metabolic fluxes is given by bilinear balance equations, flux determination from this data set requires the numerical solution of a nonlinear inverse problem. To this end, a general algorithm for flux estimation from metabolic carbon labelling experiments based on the least squares approach is developed in this contribution and complemented by appropriate tools for statistical analysis. The linearization technique usually applied for the computation of nonlinear confidence regions is shown to be inappropriate in the case of large exchange fluxes. For this reason a sophisticated compactification transformation technique for nonlinear statistical analysis is developed. Statistical analysis is then performed by computing ...

CERN Document Server

The evolution of FRI jets has been long studied in the framework of the FRI-FRII dichotomy. In this paper, we test the present theoretical and observational models via a relativistic numerical simulation of the jets in the radio galaxy 3C 31. We use the parameters derived from the modelling presented by \\cite{lb02a,lb02b} as input parameters for the simulation of the evolution of the source, thus assuming that they have not varied over the lifetime of the source. We simulate about 10 % of the total lifetime of the jets in 3C 31. Realistic density and pressure gradients for the atmosphere are used. The simulation includes an equation of state for a two-component relativistic gas that allows a separate treatment of leptonic and baryonic matter. We compare our results with the modelling of the observational data of the source. Our results show that the bow shock evolves self-similarly at a quasi-constant speed, with slight deceleration by the end of the simulation, ...

Energy Technology Data Exchange (ETDEWEB)

We consider the supercooled Stefan problem with a general anisotropic curvature- and velocity-dependent boundary condition on the moving interface. We present numerical methods, based on an integral equation formulation and including a new algorithm for moving curves with curvature-dependent velocity. These methods compute a periodic interface with {ital O}({Delta}{ital t}) accuracy, where {Delta}{ital t} is the time step. Previous work has been limited to short time spans and achieved slightly less than {ital O}({Delta}{ital t}{sup 1/2}) accuracy. Accurate numerical results are seen to agree with the predictions of linear stability theory. This agreement has eluded previous authors, because their numerical methods suffered from grid effects and their linear stability theory was incorrect. We study the long-time evolution of an unstable interface. Our computations exhibit the beginnings of a sidebranching instability when the boundary condition includes anisotropy ...

CERN Document Server

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.

International Nuclear Information System (INIS)

The top quark antiquark production system in the dilepton decay channel is described by a set of equations which is nonlinear in the unknown neutrino momenta. Its most precise and least time consuming solution is of major importance for measurements of top quark properties like the top quark mass and tt spin correlations. The initial system of equations can be transformed into two polynomial equations with two unknowns by means of elementary algebraic operations. These two polynomials of multidegree two can be reduced to one univariate polynomial of degree four by means of resultants. The obtained quartic equation is solved analytically.

International Nuclear Information System (INIS)

The lakes formed in the Amazon floodplain play a key role in sustaining the highest populational density in the Amazon Basin because of its high primary productivity. The water volume in these lakes vary greatly during the year and possible water sources are the feeding river during high water state of the hydrographs and local inputs from its drainage basin. As these sources are chemically distinct, to understand the nutrient dynamics in these lakes, it is very important to know the relative water contribution of different water sources. In this paper we use "1"8O as a conservative to know the water sources for the lakes and our hypothesis is that: 1) In the high water stage lake water is derived predominantly from its feeding river. Therefore their #delta#"1"8O should be similar. 2) In the low water stage lake water is a mixture of watershed and last flood event waters. In this case, lake #delta#"1"8O should be indicative of the relative contribution of each ...

International Nuclear Information System (INIS)

The present situation of the researches that have been made for predicting the process and range of diffusion of warm drain is reviewed. This review is divided into eight sections. The first section deals with the present situation of warm drain from power plants. For the establishment of drainage standard, there are many difficult problems to be solved because water temperature differs in its nature from other regulation items. In the second section, the process of diffusion and cooling of warm drain is explained. The third section deals with the diffusion characteristics of warm drain in Japanese coastal sea due to water temperature. Two types of diffusion are known. One is dominant irregular current, and the other is periodical reciprocating stream. The fourth section deals with the methods of prediction of diffusion. Research methods and simulation models are described. The fifth section deals with the prediction with single layer model. Fundamental consideration on the simulation ...

International Nuclear Information System (INIS)

We investigate the possibility of kaon condensation in the dense interior of neutron stars through the s-wave interaction of kaons with nucleons. We include nucleon-nucleon interactions by using simple parametrizations of realistic forces, and include electrons and muons in #beta#-equilibrium. The equation of state above the condensate threshold is derived in the mean field approximation. The conditions under which kaon condensed cores undergo a transition to quark matter containing strange quarks are also established.The critical density for kaon condensation lies in the range (2.3-5.0)#rho#_0, where #rho#_0=0.16 fm"-"3 is the equilibrium density of nuclear matter. The critical density depends largely on the value of the strangeness content of the proton, the size of which is controversial. For too large a value of the strangeness content, matter with a kaon condensate is not sufficiently stiff to support the lower limit of 1.44 M_o_e_d_o_t for a neutron star. ...

International Nuclear Information System (INIS)

The next relativistic correction to #alpha# to for bound state mass of two charged scalar particles is calculated in the quantum scalar electrodynamics by the functional integral method. Contribution of the ''nonphysical'' time variable turned out to be important and leads to nonanalytic dependence of the bound state mass on #alpha#. In conclusion, one can say that the functional approach is the best mathematical representation to preserve the gauge invariance. The lowest approximation of this functional representation is the pure nonrelativistic Feynman path integral representation of the nonrelativistic Schroedinger equation. The functional integral representation shows that any regular series for next corrections to #alpha# does not exist and these corrections cannot be reduced to some terms of the nonrelativistic potential in the Schroedinger picture. In other words, the ''nonphysical'' time coordinate is important and leads to corrections which are not ...

International Nuclear Information System (INIS)

To develop the advanced subchannel analysis code, the dominant factors that influence the boiling transitional process must be taken into account in the mechanistic constitutive equations based on the flow geometries and the fluid properties. The dominant factors that influence the boiling transitional processes are (1) the gas-liquid re-distribution by cross flow, (2) the liquid film dryout, (3) the two-phase flow regime transition, (4) the droplet deposition, and (5) the spacer-droplet interaction. At first, we indicated the strategy for the development of the constitutive equations for the five dominant factors based on the experimental database by the latest measurement technique and the latest computational fluid dynamics method. Then, the problems of the present constitutive equations and the improvement plan of the constitutive equations were indicated. Finally, the layered structure for the ...

International Nuclear Information System (INIS)

The static thermophysical properties of the molten salt system like LiF-NaF-BeF_2 influence the design and construction of the fuel salt and coolant in the Molten Salt Reactor for the new generation. In this paper, the equation of state of the ternary system 0.15LiF-0.58NaF-0.27BeF_2, over the temperature range from 873.15K to 1073.15K at one atmosphere pressure, is described by using modified Peng-Robinson equation. The density of the ternary system is evaluated by this equation directly, and compared with the experimental data. Base on the equation of state, the other static thermophysical properties such as the enthalpy, entropy and heat capacity at constant pressure are evaluated by the fugacity coefficient and residual function methods respectively. The density calculated by Peng-Robinson equation is in highly agreement with the experimental data, and the enthalpy, entropy and ...

Science.gov (United States)

The static thermophysical properties of the molten salt system like LiF-NaF-BeF2 influence the design and construction of the fuel salt and coolant in the Molten Salt Reactor for the new generation. In this paper, the equation of state of the ternary system 0.15LiF-0.58NaF-0.27BeF2, over the temperature range from 873.15K to 1073.15K at one atmosphere pressure, is described by using modified Peng-Robinson equation. The density of the ternary system is evaluated by this equation directly, and compared with the experimental data. Base on the equation of state, the other static thermophysical properties such as the enthalpy, entropy and heat capacity at constant pressure are evaluated by the fugacity coefficient and residual function methods respectively. The density calculated by Peng-Robinson equation is in highly agreement with the experimental data, and the enthalpy, entropy and ...

UK PubMed Central (United Kingdom)

Allometric equations are often used to extrapolate traits in animals for which only body mass estimates are known, such as dinosaurs. One important decision can be whether these equations should be...Full Text Available

Science.gov (United States)

... However, the solution of the Stefan problem in this context becomes abruptly complicated, since within each zone, the differential equations ...

Science.gov (United States)

... analyze importance sampling schemes for stochastic networks ... an interesting mixed open/closed queueing network ... with a suitable terminal condition ...

International Nuclear Information System (INIS)

... aromatics barium dienes heat exchangers ion exchange ions kinetic equations

Science.gov (United States)

... Abstract : The purpose of this study was to evaluate the tooling cost estimating equation of the DAPCA III model and determine if more accurate ...

Science.gov (United States)

the progressive simplification of composite sandwich panels. The equations ... The thermal and structural behavior of composite sandwich panels with a ...

Science.gov (United States)

... Following the application of the error covariance update equation, we apply the following operation to the parameter covariance matrix P+: ...

CERN Document Server

In this Letter the bound states of (2+1) Dirac equation with the cylindrically symmetric $\\delta (r-r_{0})$-potential are discussed. It is surprisingly found that the relation between the radial functions at two sides of $r_{0}$ can be established by an SO(2) transformation. We obtain a transcendental equation for calculating the energy of the bound state from the matching condition in the configuration space. The condition for existence of bound states is determined by the Sturm-Liouville theorem.

Energy Technology Data Exchange (ETDEWEB)

Sandia National Laboratories maintains several libraries of equation of state tables, in a modified Sesame format, for use in hydrocode calculations and other applications. This report discusses one of those libraries, the seslan file, which contains 78 tables from the Los Alamos equation of state library. Minor changes have been made to these tables, making them more convenient for code users and reducing numerical difficulties that occasionally arise in hydrocode calculations.

Energy Technology Data Exchange (ETDEWEB)

We report on the calculation of multi-loop Feynman integrals for single-scale problems by means of difference equations in Mellin space. The solution to these difference equations in terms of harmonic sums can be constructed algorithmically over difference fields, the so-called {pi}{sigma}{sup *}-fields. We test the implementation of the Mathematica package Sigma on examples from recent higher order perturbative calculations in Quantum Chromodynamics. (orig.)

Energy Technology Data Exchange (ETDEWEB)

The purpose of this study was to evaluate the ability to enhance distribution of eelgrass (Zostera marina) in the Columbia River Estuary to serve as refuge and feeding habitat for juvenile salmon, Dungeness crab, and other fish and wildlife. We strongly suspected that limited eelgrass seed dispersal has resulted in the present distribution of eelgrass meadows, and that there are other suitable places for eelgrass to survive and form functional meadows. Funded as part of the Bonneville Power Administration's call for Innovative Projects, we initiated a multistage study in 2008 that combined modeling, remote sensing, and field experimentation to: (1) Spatially predict habitat quality for eelgrass; (2) Conduct experimental plantings; and (3) Evaluate restoration potential. Baseline in-situ measurements and remote satellite observations were acquired for locations in the Lower Columbia River Estuary (LCRE) to determine ambient habitat ...

Science.gov (United States)

The aim of this paper is to illustrate the validity and efficiency of iterative methods for solving large linear systems arising from the finite element discretization of the equation governing conduction-controlled solidification processes. Starting from the basic enthalpy equation, two alternative formulations are obtained and fixed-grid finite element discretizations are developed. These discretizations yield a set of nonlinear equations that are linearized using the Newton-Raphson scheme. The linearized equations are used as a basis for evaluating different iterative methods of the conjugated gradient type. Symmetric scaling and incomplete factorization preconditioning of the linear equations are used to improve the convergence properties of the iterative methods. Vectorization and parallelization are also employed to make full use of the CRAY-2 supercomputer. The results ...

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.).

British Library Electronic Table of Contents (United Kingdom)

The molten salt reactor (MSR), which is one of the generation IV reactors, can meet the demand of transmutation and breeding. The thermodynamic properties of the molten salt system like LiF-NaF-BeF2 influence the design and construction of the fuel salt and coolant in the MSR for the new generation. In this paper, the equation of state of the ternary system 15%LiF-58%NaF-27%BeF2, over the temperature range from 873.15 to 1 073.15 K at one atmosphere pressure, is described using a modified Peng-Robinson (PR) equation. The densities of the ternary system and its components are estimated by this equation directly, and compared with the experimental data. Based on the equation of state, the other thermodynamic properties such as the enthalpy, entropy and heat capacity at constant pressure are ...

CERN Document Server

Single-case and Small-n Experimental Designs

Science.gov (United States)

... projectile. Descriptors : *Combustible cartridge cases, *Patents, Laminates, Ammunition propellants, Nitrocellulose. Subject ...

International Nuclear Information System (INIS)

The interfacial area transport equation dynamically models the two-phase flow regime transitions and predicts continuous change of the interfacial area concentration along the flow field. Hence, when employed in the numerical thermal-hydraulic system analysis codes, it eliminates artificial bifurcations stemming from the use of the static flow regime transition criteria. Accounting for the substantial differences in the transport phenomena of various sizes of bubbles, the two-group interfacial area transport equations have been developed. The group 1 equation describes the transport of small-dispersed bubbles that are either distorted or spherical in shapes, and the group 2 equation describes the transport of large cap, slug or churn-turbulent bubbles. The source and sink terms in the right-hand-side of the transport equations have been established by mechanistically modeling the ...

Energy Technology Data Exchange (ETDEWEB)

Full text of publication follows: Different to the very simple class of homogeneous non-equilibrium models (HEM) an one dimensional thermal-hydraulic theoretical drift-flux based and thus non-homogeneous coolant channel model and, as a result, an in itself complete thermal-hydraulic coolant channel module CCM have been established allowing to simulate in a very general way the steady state and transient behaviour of the most important parameters of a single- or two-phase fluid flowing within any type of heated or non-heated coolant channel (with an eventually varying cross flow area). To avoid mathematical discontinuities at the transition from single- to two-phase flow the coolant channel will, in its general form, be split into different regions, i.e. be looked as a basic channel (BC) which can consist of a number of different flow regimes and can, accordingly, be subdivided into a number of sub-channels (SC-s). All of them belong, obviously, to only two types of SC-s, a SC with an ...

Energy Technology Data Exchange (ETDEWEB)

The productivity characteristics of perforated gas well are investigated for the first time under real conditions considering drilling damage, perforation compacted zone and anisotropy in this paper. A non-linear finite element model for non-Darcy flow in perforated gas reservoir is built to consider the effects of twelve factors on productivity ratio (PR) of perforated gas well. The nomographs for two cases are presented in this paper. Compared with previous works, the following results are found: (1) the effects of drilling damage and perforation diameter on PR of gas well are more remarkable than that in oil well; (2) any perforation phase may be the best or the worst as to perforating gas well, which depends on anisotropy, drawdown and that whether or not perforations are beyond drilling damage zone; (3) the increase in perforation length hardly improve PR of gas well when perforations have been beyond drilling damage zone. The new method is suggested to ...

Energy Technology Data Exchange (ETDEWEB)

Shielding calculations for neutron-gamma radiation are usually done by using the full Theory of Transport or the Monte Carlo Techniques. After some works based on the Albedo Method, the shielding calculations for neutron-gamma radiation have a reliable tool with great didactical value which shows its clarity and simplicity for the resolution of cases that involve neutrons and photon shielding in nonmultiplying media. The excellent results of these works have motivated the elaboration and the development of this study that will be presented in this dissertation. The balance of a neutronic current entering a shield of two layers considering the coupling neutron-gamma will be determined by the Albedo Method. The shield will be composed of a layer of iron and another one of manganese with 10 cm of thickness each. The arrays of the materials coefficients will be obtained from the ANISN code. ANISN is a one dimensional deterministic code that is based on transport ...

Energy Technology Data Exchange (ETDEWEB)

Aging degradation is a phenomenon we all experience throughout life, both on a personal basis and in business. Many industries have been successful in postponing the inevitable impact on their related systems and components through programs to maintain long-term reliability, maintainability and safety. However, this has not always been the case for nuclear power. While all power plants are experiencing the world trend of increasing operating costs with age, few (if any) have been able to fully define the parameters that solve the aging equation, particularly in relation to major components. Inspection and preventive maintenance have not been effective in predicting life-limiting degradation and failure. In CANDU nuclear plants, utilities are taking a comprehensive approach in dealing with the aging problem. Programs have been established to identify the current condition and degradation mechanisms of critical components, the failure of which ...

Energy Technology Data Exchange (ETDEWEB)

In Chapter 2 the evaluation methods of input ground motions for seismic design are summarized and especially the empirical and/or stochastic Green's function method is introduced in order to calculate the broadband strong motions based on the physically reliable fault model. Then the evaluation procedure is applied to the Ulsan fault system near which Wolsung nuclear power plant is located. In chapter3, the PSHA case studies for Wolsung NPP site are performed, in which multi attenuation equations of earthquake motions are adopted. Moreover, the results of seismic hazard evaluation for eight sites in Korea are shown. In addition, the contents of two literatures on seismic hazard evaluation are introduced. In chapter 4, isolation devices for emergency diesel generator have been studied. The models for friction-pendulum system, high damping rubber bearing and natural rubber bearing are developed and implemented in the TDAP-III. Then, the ...

Energy Technology Data Exchange (ETDEWEB)

Ferrofluids have promising potential for heat transfer applications, since advective transport in a ferrofluid can be readily controlled by using an external magnetic field. However, unlike conventional free or forced convection, ferrohydrodynamic convection is not yet well characterized. A full understanding of the relationship between an imposed magnetic field, the resulting ferrofluid flow, and the temperature distribution is a prerequisite for the proper design and implementation of applications involving thermomagnetic convection. The literature variously assumes constant magnetic fields, does not completely represent the variation in the imposed field, or its descriptions are inaccurate, since the fields do not comply with the Maxwell's equations of electromagnetism. We address this by simulating two-dimensional forced convection heat transfer in a channel with a ferrofluid that is under the influence of a two-dimensional magnetic field created by a ...

Energy Technology Data Exchange (ETDEWEB)

By the use of thermoluminescence (TL) technique we have investigated the effects of gamma irradiation on the TL response of pulverised natural zircon (ZrSiO{sub 4}) crystals from the Jos Plateau area of Nigeria. Samples of the mineral that have been previously annealed were artificially irradiated to moderate gamma doses ranging from 0.7 to 4.9 Gy using {sup 60}Co gamma-source of the CERD. The natural TL glow curve of the sample showed a broad peak around 250 deg. C, the peak was however, found to be completely erased following heat treatment for 2 h at 400 deg. C. The TL glow curves of the artificially irradiated samples showed two overlapping glow peaks at about 112 deg. C (peak I) and at 156 deg. C (peak II) within the dose range utilised. We note that the natural TL peak appears to be different from the two observed in the artificially irradiated sample. The TL response with dose was investigated for the two glow peaks found in the artificially irradiated sample and in both ...

Energy Technology Data Exchange (ETDEWEB)

A series of stainless austenitic and ferritic materials was exposed for 100 days to boiling nitric acid which contained no corrosion products; the corrosion rates and depths of the grain boundary attack were observed. - Provided the structure is precipitation-free, the following are suitable for long-term exposure; the austenitic steels X 2 CrNi 19 11, X 1 CrNi 25 21, X 1 CrNiMoN 25 22 2 and X 1 NiCrMoCu 31 27 4, the practically Mo-free and Cu-free development steel X 1 NiCr 31 27, and the highly Mo-alloyed variant X 1 NiCrMoCu 31 27 5. - In the case of alloy NiCr21Mo it is advisable to limit the concentration and/or temperature of the nitric acid. - The 'superferrite' X 1 CrNiMoNb 28 4 2, the Japanese steel X 1 CrNiNb 30 2 and the austenitic steels X 2 CrNiMoN 17 13 3 and X 1 CrNiMoN 25 22 2 in the version with high nickel content are unsuitable. The decisive factor affecting the surface-removal rates is the chromium content of the material. The ...

Energy Technology Data Exchange (ETDEWEB)

External radiation dose exceeded 5,000 rad in three cases, and intravaginal radiation dose was 5,000 rad in one case. Radiation damage including perforation was seen in the end of ileus in one case, in the sigmoid and rectum in two cases, and in the end of ileus, sigmoid and rectum in the last case. Satisfactory results were obtained by the removal of the ileocecum in the case of the damage in the end of the ileus. However, only colostomy was performed for the damage in the sigmoid and rectum.

Science.gov (United States)

In this dissertation, a two-phase, air-droplets, dilute, turbulent, and steady state flow in a horizontal rectangular channel, is modeled and numerically simulated using a modified KIVA-3V code. The deposition of different sizes of droplets on the walls of the channel is also studied. In this model, the interaction effects between the phases (two-way coupling) are considered by source terms in the momentum and energy equations for the continuous phase and by the instantaneous local velocity of the air in the droplet equation of motion, which includes the aerodynamic and gravitational forces. The turbulence is modeled by a k-? model. The interaction effects between the turbulence and the dispersed droplets are also taken into account. The effects of the turbulence on the droplets are modeled by a fluctuating component added to the local air velocity in the droplet equation of motion. The effects of the droplets on the ...

Energy Technology Data Exchange (ETDEWEB)

Formulation is introduced for discretizing a boundary integral equation into an indirect boundary element method for the solution of 3-dimensional topographic problems. Yokoi and Takenaka propose an analytical solution-capable reference solution (solution for the half space elastic body with flat free surface) to problems of topographic response to seismic motion in a 2-dimensional in-plane field. That is to say, they propose a boundary integral equation capable of effectively suppressing the non-physical waves that emerge in the result of computation in the wake of the truncation of the discretized ground surface making use of the wave field in a semi-infinite elastic body with flat free surface. They apply the proposed boundary integral equation discretized into the indirect boundary element method to solve some examples, and succeed in proving its validity. In this report, the equation is expanded to ...

Energy Technology Data Exchange (ETDEWEB)

Installation of new Shortwave Spectrometer for permanent operation at SGP - In May 2006 the new ShortWave Spectrometer (SWS) was installed in the Optical Trailer at the Southern Great Plains Central Facility SGP on 27 April 2006. The SWS began full operation 28 April 2006 and has run continuously to the present. Over 25 GB of spectra has been collected, calibrated and archived. 3-D radiative transfer simulations - Retrieved fields of cloud optical thickness and effective radius to from the MODIS Airborne Simulator were used to reproduce 3D cloud fields that were used a input to 3D radiative transfer simulations and then compared with simultaneous Solar Spectral Flux Radiometer (SSFR) spectral irradiance measurements. The influence of both horizontal and vertical cloud structure, using accurate versus approximated optical properties in the radiative transfer model on the modeled irradiance was examined, as was the influence of using the full phase function versus using approximations of ...

Energy Technology Data Exchange (ETDEWEB)

Solving a coupled set of differential equations formed by a Fokker-Planck equation and a set of master equations, we can calculate the number of 'prefission' particles (n,p,a), which are emitted from a hot, rotating compound nucleus. Furthermore we get their energy spectrum and their angular momentum distribution. Our results depend sensitively on the friction parameter {gamma}, the inertia M and the fission potential. Comparing the results with experimental data, we may learn something about these quantities. The investigations have been made for {sup 187}Ir and {sup 185}Os. (orig.).

CERN Document Server

Constrained quantum dynamics is used to propose a nonlinear dynamical equation for pure states of a generalized coarse-grained system. The relevant constraint is given either by the generalized purity or by the generalized invariant fluctuation, and the coarse-grained pure states correspond to the generalized coherent i.e. generalized nonentangled states. Open system model of the coarse-graining is discussed. It is shown that in this model and in the weak coupling limit the constrained dynamical equations coincide with an equation for pointer states, based on Hilbert-Schmidt distance, that was previously suggested in the context of the decoherence theory.

CERN Document Server

Big Bang nucleosynthesis requires a fine balance between equations of state for photons and relativistic fermions. Several corrections to equation of state parameters arise from classical and quantum physics, which are derived here from a canonical perspective. In particular, loop quantum gravity allows one to compute quantum gravity corrections for Maxwell and Dirac fields. Although the classical actions are very different, quantum corrections to the equation of state are remarkably similar. To lowest order, these corrections take the form of an overall expansion-dependent multiplicative factor in the total density. We use these results, along with the predictions of Big Bang nucleosynthesis, to place bounds on these corrections.

British Library Electronic Table of Contents (United Kingdom)

Pippenger (2011) recently proposed a solution to the longstanding forward-bias puzzle. He argues that the puzzling estimates obtained using the standard equation for the efficient markets hypothesis are due to omitted variable bias. He identifies the missing variables as the future change in the forward exchange rate and the future interest differential. When these are added to the standard equation, he finds a one-to-one relationship between the future change in the spot rate and the forward premium. However, we argue that his equation can only test covered interest parity and offers no insight into the forward-bias puzzle.

International Nuclear Information System (INIS)

HAMILTON is a computer code performing all algebraic operations necessary for an analytic determination of the power series of the Hamiltonian equations of motion in the electromagnetic fields with at least one plane of symmetry. It is written entirely in FORTRAN in order to achieve fast machine performance, a requirement which is essential due to the complexity of the equations of motion in higher orders. HAMILTON is considerably faster than common more versatile formula manipulators and uses noticeably less storage. Besides the mere solution of the equations of motion, HAMILTON also produces FORTRAN code compatible with the program COSY 5.0 allowing the computation of matrix elements of individual optical elements and their concatenation. The produced FORTRAN code is highly optimized and on average requires only 30% of the execution time of a handwritten comparable code. (orig.).

Energy Technology Data Exchange (ETDEWEB)

In this paper, research on pass schedule of tandem cold mill (TCM) is carried out. According to load (reduction, rolling force, motor power) balance, non-linear equations set with variables of inter-stand thickness is constructed. The pass schedule optimization is carried out by solving the non-linear equations set. As the traditional method, the Newton-Raphson method is used for solving the non-linear equations set. In this paper a new simple method is brought up. On basis of the monotone relations between thickness and load, the inter-stands thickness is adjusted dynamically. The solution of non-linear equations set can be converged by iterative calculation. This method can avoid the derivative calculation used by traditional method. So, this method is simple and calculation speed is high. It is suitable for on-line control. (author)

UK PubMed Central (United Kingdom)

The geographical arrangement of people on the surface of the earth is approximated by a mathematical equation of 361 terms. This is a convenient form for comparison with other distributions and for...Full Text Available

International Nuclear Information System (INIS)

Electrodynamic properties of processes involving charged tachyons are considered in connection with the real superluminal Lorentz transformations. These transformations are interpreted without making use of a tachon corridor and so that the interaction of tachyons from subluminal sources with ordinary matter is Lorentz invariant. Transformations of the electromagnetic fields are deduced from the extended principle of relativity. They are necessary in order to obtain the field equations that are satisfied by fields obeying Maxwell equations relative to superluminal inertial frames with respect to subluminal inertial frames. The field equations allow one to get the dependence of the field energy density and its current on the field strengths. The resulting equations are applied in a discussion of plane electromegnetic waves, sent out by a superluminal source. The obtained electrodynamics allows the ...

Science.gov (United States)

a large number of polycrystalline materials, has the fot_ of the. Manson-Coffin- Basquin equation, wherein the exponents of the two power law terms are given ...

UK PubMed Central (United Kingdom)

The solution of the continuity equation in practical applications is examined, and the values needed for approximate solutions are indicated. Models are adequate for investigating what could happen,...Full Text Available

International Nuclear Information System (INIS)

The objective of this research is to replace the current flow regime dependent interfacial area correlations in the thermal-hydraulic system analysis codes with an interfacial area transport equation, which dynamically models two-phase flow regime transitions. The interfacial area transport equation applicable to bubbly flows has been developed based on mechanistic modeling of bubble interaction. Detailed experiments have been performed for bubbly conditions in 2.54, 5.08 and 10.16-cm ID pipes to benchmark the model. This new approach predicts the continuous changes of the interfacial area and will eliminate artificial bifurcations stemming from the use of the static flow regime transition criteria. This paper presents the preliminary results of incorporating the one-group transport equation into the USNRC consolidated code, currently under development. The framework for the two-group transport equation ...

CERN Document Server

A well-known diffuse interface model consists of the Navier-Stokes equations nonlinearly coupled with a convective Cahn-Hilliard type equation. This system describes the evolution of an incompressible isothermal mixture of binary-fluids and it has been investigated by many authors. Here we consider a variant of this model where the standard Cahn-Hilliard equation is replaced by its nonlocal version. More precisely, the gradient term in the free energy functional is replaced by a spatial con- volution operator acting on the order parameter phi. Therefore the coupling with the Navier-Stokes equations is difficult to handle even in two spatial dimensions because of the lack of regularity of phi. We establish the global existence of a weak solution.

CERN Document Server

The possible physical linkage between galactic cosmic rays intensity and the Earth's cloud cover is discussed using the analysis of the first indirect aerosol effect (Twomey effect) and its experimental representation as the dependence of average cloud droplet effective radius on aerosol index characterizing the aerosol concentration in the atmospheric air column of unit section. It is shown that the basic kinetic equation of the Earth's climate energy-balance model is described by the bifurcation equation (with respect to the temperature of the Earth's surface) in the form of fold catastrophe with two governing parameters defining the variations of insolation and Earth's magnetic field (or galactic cosmic rays intensity in the atmosphere), respectively. The principle of hierarchical climatic models construction, which consists in the structural invariance of balance equations of these models evolving on the different time ...

British Library Electronic Table of Contents (United Kingdom)

The DN-tracking method is used to prove the existence of a closed trajectory in a quadratic system of ordinary differential equations in three dimensions.

CERN Document Server

A new mathematical framework is formulated to derive the effective equations of motion for the constrained quantum system which possesses an internal clock. In the realm close to classical behavior, the quantum evolution is approximated by a finite system of coupled but ordinary differential equations adhered to the weakly imposed Hamiltonian constraint. For the simplified version of loop quantum cosmology in the Bianchi I model with a free massless scalar filed, the resulting effective equations of motion affirm the bouncing scenario predicted by the previous studies: The big bang singularity is resolved and replaced by the big bounces, which take place up to three times, once in each diagonal direction, whenever the directional density approaches the critical value in the regime of Planckian density. It is also revealed that back-reaction arises from the quantum corrections and modifies the precise value of the ...

UK PubMed Central (United Kingdom)

Diffusion in the extracellular space (ECS) of the brain is constrained by the volume fraction and the tortuosity and a modified diffusion equation represents the transport behavior of many molecules...Full Text Available

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 ...

International Nuclear Information System (INIS)

The presence of an extraneous light gas must be taken into account in consideration of centrifugal separation of uranium isotopes, when there is inevitable leakage of such a gas through the gas sealant and/or leaking in from the atmosphere. Consideration is first given to the influence of the presence of the light gas on the maximum separative power. Then the basic equation for isotope separation containing a light gas is derived from Hirschfelder's diffusion equations. This equation is solved and the separative performance is expressed in terms of the shape factor and reflux parameter. The formulas for expressing the flow configurations of the gases are obtained for a simple model in which inflow and outflow prevail throughout the centrifuge. The corresponding equation for a model in which the gases flow in two concentric thin streams is also derived. It is concluded that the influence provided by the ...

International Nuclear Information System (INIS)

The set of nonlinear equations describing the standard model kinematics of the top quark antiquark production system in the dilepton decay channel has at most a fourfold ambiguity due to two not fully reconstructed neutrinos. Its most precise solution is of major importance for measurements of top quark properties like the top quark mass and tt spin correlations. Simple algebraic operations allow one to transform the nonlinear equations into a system of two polynomial equations with two unknowns. These two polynomials of multidegree eight can in turn be analytically reduced to one polynomial with one unknown by means of resultants. The obtained univariate polynomial is of degree 16. The number of its real solutions is determined analytically by means of Sturm's theorem, which is as well used to isolate each real solution into a unique pairwise disjoint interval. The solutions are polished by seeking the sign change of the ...

International Nuclear Information System (INIS)

We consider particle evaporation and fission of an ensemble of hot, rotating compound nuclei as a stochastic process. We derive a set of coupled differential equations formed by a Fokker-Planck equation describing fission, and master equations for calculating particle evaporation. From these equations, we are able to determine multiplicities of prefission neutrons, protons and #alpha#-particles, their energy spectra and their angular momentum distributions. A comparison of our results with experimental data provides us with information regarding the reduced friction coefficient #beta#, the fission barrier height and the level density parameter. For different iridium isotopes, ("1"8"1","1"8"5","1"8"7Ir), "1"8"5Os and "1"5"8Er, we obtain as an upper limit #beta##<=#8.0x10"2"1 s"-"1. (orig.).

Science.gov (United States)

... 3. Critical constants, and triple point 4. Compressibility isothermE S. Density, molar volume (of a condensed phase) 6. Equations of state, general ...

DEFF Research Database (Denmark)

Udgivelsesdato: March

KoreaScience (Korea)

Full Text Available

Science.gov (United States)

... Descriptors : *AMMUNITION, *CONTAINERS, *CARTRIDGES, *COMBUSTIBLE CARTRIDGE CASES, *CARTRIDGE CASES, *ANTITANK ...

Science.gov (United States)

A production method was developed to fibrillate commercial acrylic fiber suitable for the manufacture of combustible cartridge cases.

Science.gov (United States)

... initiation) Pressed e Combustible Cartridge Cases Quantity-Distance Table ... Combustible Cartridge Cases, and Electrically Exploded Devices in ...

Science.gov (United States)

... ADA185428. Title : Combustible Cartridge Cases, an Account of the Current Technology and Proposals for Future Development. ...

Science.gov (United States)

... Abstract : The beater additive process for manufacturing combustible cartridge cases is superior to the solvent resin, post-impregnation process. ...

Science.gov (United States)

... reverse aeb 9 nescesmsy amI identify by block number) Provides procedures for evaluating metal, consumable, and combustible cartridge cases. ...

British Library Electronic Table of Contents (United Kingdom)

An operator of FE-closure is introduced on the set of functions of a multivalued logic based on the systems of functional equations. It is proved that, for every k ? 2, the FE-closure operator generates a finite classification on the set P k of functions of k-valued logic. The least class in this classification is shown to be the class H k of all homogeneous functions. Also a series of corollaries are obtained concerning the finite FE-generating sets in the FE-closed classes.

Energy Technology Data Exchange (ETDEWEB)

The spin-up problem of He II in a cylindrical container of infinite and finite height is considered on the basis of the Hall-Vinen-Bekarevich-Khalatnikov equations. The boundary conditions include the assumption of the generation of superfluid vorticity at the walls of the container. In the lowest order approximation with respect to viscosity, the problem is reduced to a Stefan problem for a single diffusion-type equation. Examples of solutions are given which show that the method used gives a good insight into the physical development of the spin-up flow.

International Nuclear Information System (INIS)

For a system of noninteracting electrons scattered by resonant levels of neutral impurities, we show that virial and quasiparticle corrections have nearly equal magnitudes. We propose a modification of the Boltzmann equation that includes quasiparticle and virial corrections and discuss their interplay on a dielectric function. copyright 1997 The American Physical Society.

Energy Technology Data Exchange (ETDEWEB)

To simulate the effects of the breakup channel on rearrangement amplitudes, the conventional coupled-reaction-channel (CRC) expansion is augmented by pseudoreaction channels. The construction of the projector for the extended CRC space is discussed, and transition-operator equations on this space are given. By solving the full and post-approximation forms of the CRC equations for a model three-particle problem, the crucial role played by the nonorthogonality terms is demonstrated.

CERN Document Server

We study the dynamics of the scalar tensor cosmological model in the presence of tachyon field. In an alternative approach, in two exponential and power law form of the scalar field functions in the model, field equations are solved by simultaneously best fitting the model parameters with the most recent observational data. This approach gives us an observationally verified interpretation of the dynamics of the universe. We then discuss the best fitted of equation of state parameter, the statefinder parameters and the reconstructed scalar field in the model.

International Nuclear Information System (INIS)

A method of auxiliary spectrum is modified so that matrix elements of the reaction matrix in final nuclei are determined by means of algebraic operations alone (inversion of matrices). No differential equations need to be solved; Pauli's exclusion principle is accurately taken into account. A single-particle potential may be of any kind, but a two-particle interaction must have no solid core.

International Nuclear Information System (INIS)

The loop quantum cosmology 'improved dynamics' of the Bianchi type IX model are studied. The action of the Hamiltonian constraint operator is obtained via techniques developed for the Bianchi type I and type II models, no new input is required. It is shown that the big bang and big crunch singularities are resolved by quantum gravity effects. We also present effective equations which provide quantum geometry corrections to the classical equations of motion.

Energy Technology Data Exchange (ETDEWEB)

The optimization of ethylene furnace operations using cfd-based simulations has been addressed. The optimization problems have been cast into various formulations: the Multidisciplinary feasible (MDF) approach, the All-At-Once (AAO) approach and the Individual discipline feasible (IDF) approach. These approaches mainly differ in their handling of the state equations as constraints, and hence some of the formulations place restrictions on the methods used to solve the state equations.

CERN Document Server

In this paper we first establish global pointwise time-space estimates for a class of oscillatory integrals. Then, we use them to establish $L^p-L^q$ estimates for a class of higher order wave-type equations of the form $\\partial_{tt}u+P(D_{x})u=0$, where the symbol $P(\\xi)$ is a real non-degenerate elliptic polynomial of ${\\bf R}^n$.

Energy Technology Data Exchange (ETDEWEB)

An equation is described to estimate the value of the liquid-phase mass-transfer coefficient of aqueous solutions of surfactants from its value for pure water. The form of the equation is obtained from Levich's theory, and it is adjusted to experiments of absorption of CO/sub 2/ in aqueous solutions of five pure surfactants. The reduction of the absorption rate is correlated with the interfacial pressure of the solution.

CERN Document Server

We consider the numerical discretization of the time-domain Maxwell's equations with an energy-conserving discontinuous Galerkin finite element formulation. This particular formulation allows for higher order approximations of the electric and magnetic field. Special emphasis is placed on an efficient implementation which is achieved by taking advantage of recurrence properties and the tensor-product structure of the chosen shape functions. These recurrences have been derived symbolically with computer algebra methods reminiscent of the holonomic systems approach.

Energy Technology Data Exchange (ETDEWEB)

Equations are compiled for thermal balance in which for simplification, no consideration is made for heat conductivity along the axis of the cable and dependence of losses, heat capacitance and heat conductivity on temperature. Equations are modeled on a transistor analog calculator 42 TA. The solution to the task on the computer produced values of maximum temperature on the cable and coordinates of the point of maximum overheating. Using the analog model, one can study other parameters of the thermal mode.

UK PubMed Central (United Kingdom)

AimsTwo-dimensional echocardiographic (2DE) continuity-equation derived aortic valve area (AVA) in aortic stenosis (AS) relies on non-simultaneous measurement of left ventricular...Full Text Available

Energy Technology Data Exchange (ETDEWEB)

A formula for calculating the Dirac S-matrix for central Lorentz scalar and vector potentials is derived by use of a new amplitude-phase method. The derivation also makes use of certain invariants of the reduced 2-spinor radial Dirac equations.

Energy Technology Data Exchange (ETDEWEB)

We consider the problem of rotation in a homogeneous Bianchi type IX cosmological model. Studying the Wheeler-De Witt equation corresponding to this minisuperspace, and adopting a particular choice of the factor ordering, we are able to find a particular solution which is strongly peaked about isotropy. This result confirms all the previous investigations in this field, and suggests the conclusion that the machian nature of our universe is the natural, i.e. most likely, outcome of the Planck epoch. (orig.).

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 ...

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 ...

International Nuclear Information System (INIS)

For good image quality using ultrasound inverse scattering, one alternately solves the well-posed forward scattering equation for an estimated total field and the ill-posed inverse scattering equation for the desired object property function. In estimating the total field, error or noise contaminates the coefficients of both matrix and data of the inverse scattering equation. Previous work on ill-posed inverse ultrasonic scattering commonly used Tikhonov regularization, which considers error only in the data. The solution so obtained is not precise enough to reconstruct the quantitative internal structure of a large or high-contrast object. This paper adopts the truncated total least squares method, simultaneously considering error and noise on both sides of the inverse scattering equation, and compares it with the classical Tikhonov regularization. We show that it can substantially improve ...

Energy Technology Data Exchange (ETDEWEB)

Background: The difficulty of directly measuring cellular dose is a significant obstacle to application of target tissue dosimetry for nanoparticle and microparticle toxicity assessment. As a consequence, the target tissue paradigm for dosimetry and hazard assessment of nanoparticles has largely been ignored in favor of using metrics of exposure (e.g. ?g particle/mL culture medium, particle surface area/mL, particle number/mL). We have developed a computational model of solution particokinetics (sedimentation, diffusion) and dosimetry for non-interacting spherical particles and their agglomerates in monolayer cell culture systems. Particle transport to cells is calculated by simultaneous solution of Stokes Law (sedimentation) and the Stokes-Einstein equation (diffusion). Results: The In vitro Sedimentation, Diffusion and Dosimetry model (ISDD) was tested against measured transport rates or cellular doses for multiple sizes of polystyrene spheres (20-1100 nm), 35 nm ...

International Nuclear Information System (INIS)

The Forsmark and the Laxemar investigation areas are examined by the Swedish Nuclear Fuel and Waste Management Co. for a possible construction of a deep repository for nuclear waste. In the case of a future leakage of waste, the radioactive isotopes could end up in the ecosystems above the repository. The fate of the radionuclides and their possible radiological impacts are then highly determined by ecosystem carbon cycling. An important part of the carbon cycling is the soil carbon effluxes, and in the investigation areas soil carbon effluxes have been examined with the closed chamber technique. This paper is divided into two parts. Firstly, there were problems with the equipment measuring the soil carbon dioxide efflux, and the first part is a description of the problem, how it was corrected and its possible causes. The second part is a manual in how to analyse data and calculate annual estimates of soil carbon efflux. The field measurement by the EGM-4 is just ...

International Nuclear Information System (INIS)

Computed tomograms (CT) from 204 cases of premature and full term infants were studied. 1) In 70 infants of hyaline membrane disease, primary apnea and asymptomatic transient hypoglycemia and hypocalcemia, without any abnormal CT findings such as intracranial hemorrhage, periventricular lucency was found in 65 cases (93%) and a wide extracerebral space of the temporal lobes was found in 60 cases (86%) before 44 weeks of gestation. A wide interhemispheric fissure was found in 11 cases of extremely premature infants before term. Periventricular lucency and a wide extracerebral space of temporal lobes may not be the result of the intracranial pathological changes, but they may represent a stage of brain development. 2) In 204 cases, intracranial hemorrhage was found in 39 cases; intra ventricular 14 cases, subependymal 2 ...

Energy Technology Data Exchange (ETDEWEB)

A copyrighted computational fluid dynamics computer code, ICOMFLO, has been developed for the simulation of multiphase reacting flows. The code solves conservation equations for gaseous species and droplets (or solid particles) of various sizes. General conservation laws, expressed by elliptic type partial differential equations, are used in conjunction with rate equations governing the mass, momentum, enthalpy, species, turbulent kinetic energy, and turbulent dissipation. Associated phenomenological submodels of the code include integral combustion, two parameter turbulence, particle evaporation, and interfacial submodels. A newly developed integral combustion submodel replacing an Arrhenius type differential reaction submodel has been implemented to improve numerical convergence and enhance numerical stability. A two parameter turbulence submodel is modified for both gas and solid phases. An evaporation submodel treats ...

Energy Technology Data Exchange (ETDEWEB)

VSP reverse-time migration is a well adaptable wave equation migration method. Its control equation not only describes all-direction propagation of seismic wave but also removes interbed multiples. Clearbout's image principle is generalized to determine image conditions, real VSP data are used to determine boundary condition, and two way reflection-free wave equation is solved by making reverse-time extrapolation. In each step of extrapolation, the migration value at relevant image point is obtained by using the image condition. The complete migration of a seismic section is achieved when reverse-time extrapolation reaches the minimum image time. In this paper it is proved theoretically and practically that this method is applicable to any velocity variation and makes the migrated section have both good resolution and high S/N ratio. Besides, this method results in high processing efficiency.

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 ...

Energy Technology Data Exchange (ETDEWEB)

The space derivation term of a wave equation is accurately calculated using Fourier transform method, and the wave equation only relating to time derivative in time-space domain is derived. Solving this equation with the aid of central difference method brings the numerical recursion formula for forward modeling or reverse time migration of seismic wave field in the medium in which seismic velocity varies. The key to this method lies in introducing two auxiliary wave fields relating to velocity and spacial frequency of Fourier transform respectively after multi-dimensional Fourier transform respectively after multi-dimensional Fourier transform of space vector is made. Theoretically, this method is suitable to the forward modeling and migration of seismic wave field in complex area where seismic velocity and structure shape vary arbitrarily. Theoretical results proved this method satisfactory.

Energy Technology Data Exchange (ETDEWEB)

Equilibrium and energy equations have been developed in describing the solid-fluid transition zone during the melting and solidification of the binary alloys. Due to the existence of the transition region exhibited both solid and fluid phases at the same material point from continuum point of view, mixture theory was utilized to analyze the region. Unlike the Stefan problem, the latent heat due to the phase change appears as a source term in the heat equation. The molten fluid is treated as a thermoviscous and incompressible fluid, whereas the solid is thermoviscoplastic described by the Bodner-Partom/Walker type of constitutive equations. Thermal mechanical behaviors of the solid and the fluid phases are determined separately because of insignificant mechanical interactions between them. Volume fractions of the phases are obtained according to the equilibrium phase diagram. The simulation process of the transition zone and ...

Energy Technology Data Exchange (ETDEWEB)

Lateral plume dispersion in deep river valleys during neutral and stable lapse rate conditions can be exceptionally high due to the intense horizontal turbulence generated by prominent variations in the width, orientation and surface roughness of the valley walls. Use of the standard Gaussian dispersion equation to represent this enhanced lateral plume spread in a narrow valley may be inappropriate since consideration should also be given to the restriction of the horizontal plume spread due to impingement against the valley walls. The basic concept employed in the modification of the Gaussian dispersion equation was to assume that multiple eddy reflections occur between the valley walls in a manner similar to the vertical eddy reflections between the ground plane and an inversion layer aloft. The expressions are developed, but no attempt has been made to validate them with actual measured field data. 4 references, 3 figures, 1 table.

Energy Technology Data Exchange (ETDEWEB)

The code is being developed starting from the steady-state finite element code FENDER for the solution of the diffusion equation by extending it to become time-dependent. The numerical solution of the time-dependent multigroup diffusion equations within TRANSFUSION is performed at the present stage of development by using a backward difference scheme for the time variable, leading to a rearrangement of FENDER by adding a new loop over time steps. The code retains the multigroup coupled neutron-gamma features of FENDER, and provides a consistent two-, and quasi three-dimensional numerical solution of both static and time-dependent multigroup diffusion equations. (orig./DG)

International Nuclear Information System (INIS)

We study the thermodynamics of a one-dimensional attractive Fermi gas (the Gaudin-Yang model) with spin imbalance. The exact solution has been known from the thermodynamic Bethe ansatz for decades, but it involves an infinite number of coupled nonlinear integral equations whose physics is difficult to extract. Here the solution is analytically reduced to a simple, powerful set of four algebraic equations. The simplified equations become universal and exact in the experimental regime of strong interaction and relatively low temperature. Using the new formulation, we discuss the qualitative features of finite-temperature crossover and make quantitative predictions on the density profiles in traps. We propose a practical two-stage scheme to achieve accurate thermometry for a trapped spin-imbalanced Fermi gas.

British Library Electronic Table of Contents (United Kingdom)

In this paper, the unsteady MHD free convection heat and mass transfer of viscous fluid flowing through a Darcian porous regime adjacent to a moving vertical semi-infinite plate under Soret and Dufour effect have been examined. Viscous dissipation effects are included in the energy equation. A uniform magnetic field is applied transversely to the direction of the flow. The differential equations governing the problem have been transformed by a similarity transformation into a system of non-dimensional differential equations which are solved numerically by element free Galerkin method. The influence of Grashof number (Gr), magnetic parameter (M), heat absorption parameter (Q), permeability parameter (K), Schmidt number (Sc), Soret number (Sr), and Dufour number (Du) on the velocity, tempera...

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.

CERN Document Server

This contribution deals with identification of fractional-order dynamical systems. System identification, which refers to estimation of process parameters, is a necessity in control theory. Real processes are usually of fractional order as opposed to the ideal integral order models. A simple and elegant scheme of estimating the parameters for such a fractional order process is proposed. This method employs fractional calculus theory to find equations relating the parameters that are to be estimated, and then estimates the process parameters after solving the simultaneous equations. The said simultaneous equations are generated and updated using particle swarm optimization (PSO) technique, the fitness function being the sum of squared deviations from the actual set of observations. The data used for the calculations are intentionally corrupted to simulate real-life conditions. Results show that the proposed scheme offers a ...