Science.gov (United States)

terms of 11) numerical models, 2D SPH simulations, and analytic models. The 1D model was useful for determining the evolution over very long timescales. ...

Energy Technology Data Exchange (ETDEWEB)

There is an ever increasing demand to perform flow simulations that incorporate the complete details of geometry as well as sophisticated flow physics. This has led to the development of numerical algorithms that can simulate the actual flow phenomena with greater fidelity. However, the success of these algorithms hinges on the grid that models the geometry. Grid generation methods for 2-D models have long existed and the general lack of complexity of the simpler 2-D models has not quite challenged the efforts in this area. However, demands for generating better 3-D geometric models for flow simulations involving complex geometries have completely changed the perspective of grid generation strategies. As a consequence, grid generation efforts have earned equal significance as that of numerical solver ...

Energy Technology Data Exchange (ETDEWEB)

We demonstrate a two-dimensional device simulator for MOSFET structures that incorporates models for defect distributions and show predicted effects on device switching performance for various spatial distributions of defects in amorphous and polycrystalline silicon.

CERN Document Server

An effective computer program for three dimensional relativistic hydrodynamical model has been developed. It implements a new approach to the early hot phase of relativistic heavy-ion collisions. The computer program simulates time-space evolution of nuclear matter in terms of ideal-fluid dynamics. Equations of motions of hydrodynamics are solved making use of finite difference methods. Commonly-used algorithms of numerical relativistic hydrodynamics RHLLE and MUSTA-FORCE have been applied in simulations. To speed-up calculations, parallel processing has been made available for solving hydrodynamical equations. The test results of simulations for 3D, 2D and Bjorken expansion are reported in this paper. As a next step we plan to implement the hadronization algorithm by implementing the continuous particle emission for freeze-out and comparing it with Cooper-Frye ...

Science.gov (United States)

We show that accretion disks, both in the subcritical and supercritical accretion rate regime, may exhibit significant amplitude luminosity oscillations. The luminosity time behavior has been obtained by performing a set of time-dependent 2D SPH simulations of accretion disks with different values of ? and accretion rate. An explanation of this luminosity behavior is proposed in terms of limit-cycle instability: the disk oscillates between a radiation pressure dominated configuration (with a high luminosity value) and a gas pressure dominated one (with a low luminosity value). The origin of this instability is the difference between the heat produced by viscosity and the energy emitted as radiation from the disk surface (the well-known thermal instability mechanism). We support this hypothesis showing that the limit-cycle behavior produces a sequence of collapsing and refilling states of the innermost disk region.

Science.gov (United States)

This dissertation reports on the application of numerical optimization techniques as applied to fuel cell simulation and design. Due to the "multi-physics" inherent in a fuel cell, which results in a highly coupled and non-linear behavior, an experimental program to analyze and improve the performance of fuel cells is extremely difficult. This program applies new optimization techniques with computational methods from the field of aerospace engineering to the fuel cell design problem. After an overview of fuel cell history, importance, and classification, a mathematical model of solid oxide fuel cells (SOFC) is presented. The governing equations are discretized and solved with computational fluid dynamics (CFD) techniques including unstructured meshes, non-linear solution methods, numerical derivatives with complex variables, and sensitivity analysis with adjoint methods. Following the validation of the fuel cell model in ...

Energy Technology Data Exchange (ETDEWEB)

Porosity evolution at reactive interfaces is a key process that governs the evolution and performances of many engineered systems that have important applications in earth and environmental sciences. This is the case, for example, at the interface between cement structures and clays in deep geological nuclear waste disposals. Although in a different transport regime, similar questions arise for permeable reactive barriers used for biogeochemical remediation in surface environments. The COMEDIE project aims at investigating the coupling between transport, hydrodynamics and chemistry when significant variations of porosity occur. The present work focuses on a numerical benchmark used as a design exercise for the future COMEDIE-2D experiment. The use of reactive transport simulation tools like Hytec and Crunch provides predictions of the physico-chemical evolutions that are expected during the future ...

Science.gov (United States)

Turbulent flow in a 2D channel with repeated rectangular rib roughness was numerically simulated using a low Reynolds number form of the k-epsilon turbulence model. Friction factors and average Stanton numbers were calculated for various pitch to rib height ratios and bulk Reynolds numbers. Comparisons with experiment were generally adequate, with the predictions of friction superior to those for heat transfer. The effect of variable properties for channel flow was investigated, and the results showed a greater effect for friction than for heat transfer. Comparison with experiment yielded no clear conclusions. The turbulence model was also validated for a related problem, that of flow downstream of an abrupt pipe expansion.

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.

British Library Electronic Table of Contents (United Kingdom)

A 2-D model has been proposed to investigate the approximate estimation of the natural convection heat loss from modified cavity receiver of without insulation (WOI) and with insulation (WI) at the bottom of the aperture plane in our previous article. In this paper, a 3-D numerical model is presented to investigate the accurate estimation of natural convection heat loss from modified cavity receiver (WOI) of fuzzy focal solar dish concentrator. A comparison of 2-D and 3-D natural convection heat loss from a modified cavity receiver is carried out. A parametric study is carried out to develop separate Nusselt number correlations for 2-D and 3-D geometries of modified cavity receiver for estimation of convective heat loss from the receiver. The results show that the 2-D and 3-D are comparabl...

International Nuclear Information System (INIS)

Electrostatic plasma simulation code for 2D rectangular geometry is presented. Main distinguishing feature of the code is its orientation on the beam-plasma interaction. The code and its graphical interface were developed using MATLAB programming language. Simulation results of inhomogeneous plasma interaction with modulated electron beams of different width are compared. In case of wide beam the front of Langmuir waves generated in point of local plasma resonance is planar and in case of thin beam (or ribbon beam) the front has approximately half-circular form.

CERN Document Server

We point out that electromagnetic one-way edge modes analogous to quantum Hall edge states, originally predicted by Raghu and Haldane in 2D gyroelectric photonic crystals possessing Dirac point-derived bandgaps, can appear in more general settings. In particular, we show that the TM modes in a gyromagnetic photonic crystal can be formally mapped to electronic wavefunctions in a periodic electromagnetic field, so that the only requirement for the existence of one-way edge modes is that the Chern number for all bands below a gap is non-zero. In a square-lattice gyromagnetic Yttrium-Iron-Garnet photonic crystal operating at microwave frequencies, which lacks Dirac points, time-reversal breaking is strong enough that the effect should be easily observable. For realistic material parameters, the edge modes occupy a 10% band gap. Numerical simulations of a one-way waveguide incorporating this crystal show ...

Energy Technology Data Exchange (ETDEWEB)

2D simulation modeling of a river basin is based on notion of a River Basin as an opened non equilibrium nature system. All components of the system (elevation, water, soil properties and others) are linked to single multi - functional model. The objective of flood simulation is one the model functions. The objectives of the paper are as follows: a mean scale basin floods computing and the model validation; a large scale basin simulation and demands for its validation; results and demands for flood hazard assessment. (orig.)

Energy Technology Data Exchange (ETDEWEB)

A 2-D model has been proposed to investigate the approximate estimation of the natural convection heat loss from modified cavity receiver of without insulation (WOI) and with insulation (WI) at the bottom of the aperture plane in our previous article. In this paper, a 3-D numerical model is presented to investigate the accurate estimation of natural convection heat loss from modified cavity receiver (WOI) of fuzzy focal solar dish concentrator. A comparison of 2-D and 3-D natural convection heat loss from a modified cavity receiver is carried out. A parametric study is carried out to develop separate Nusselt number correlations for 2-D and 3-D geometries of modified cavity receiver for estimation of convective heat loss from the receiver. The results show that the 2-D and 3-D are comparable only at ...

International Nuclear Information System (INIS)

2D MEDICI simulator is used to investigate hardening solutions to single-event burnout (SEB). SEB parametric dependencies such as carrier lifetime reduction, base enlargement, and emitter doping decrease have been verified and a p"+ plug modification approach for SEB hardening of power MOSFETs is validated with simulations on actual device structures.

International Nuclear Information System (INIS)

The diffusion of ion beam injected self-interstitials (I) and their interaction with impurities in crystalline Si has been investigated and modeled. In particular, the I-substitutional carbon (C) interactions have been studied, using a molecular-beam-epitaxy grown Si_1_-_yC_y layer interposed between the shallow I-source and a deeper B-spike (marker for I-concentration). Substitutional C atoms are shown to trap I's, to be removed from their substitutional sites, and to form stable precipitates into the C-rich region. The I-trapping mechanism was quantitatively studied by a simulation code. The reactions causing trapping and deactivation are described. In addition, the boron markers approach was extended to the two dimensional (2D) diffusion. High resolution scanning capacitance microscopy was used for quantitative measurements of the 2D boron transient enhanced diffusion induced on ...

International Nuclear Information System (INIS)

Power devices such as MOSFETSs and IGBTs, include parasitic structures that can give rise to destructive failures such as breakdown and latch-up. To determine a suitable strategy for device radiation hardening, simulation software like MEDICI-2D can be used to model the effects of technological modifications and device parameters that are difficult to measure experimentally. (authors).

International Nuclear Information System (INIS)

With the objective to improve the reactor physics calculation on a 2D and 3D nuclear reactor via the Diffusion Equation, an adaptive automatic finite element remeshing method, based on the elementary area (2D) or volume (3D) constraints, has been developed. The adaptive remeshing technique, guided by a posteriori error estimator, makes use of two external mesh generator programs: Triangle and TetGen. The use of these free external finite element mesh generators and an adaptive remeshing technique based on the current field continuity show that they are powerful tools to improve the neutron flux distribution calculation and by consequence the power solution of the reactor core even though they have a minor influence on the critical coefficient of the calculated reactor core examples. Two numerical examples are presented: the 2D IAEA reactor ...

Energy Technology Data Exchange (ETDEWEB)

The WBURN (2-D, 2-group, coarse mesh) code is developed to analyze the equilibrium core characteristics of CANDU-PHWR. The equilibrium characteristics of Wolsung reactor computed by using WBURN are compared with the values given in the Wolsung FSR. The changes of equilibrium core characteristics caused by the variation of design parameters for operating conditions are also investigated. The numerical results indicate that the average discharge irradiation in the Wolsung reactor can be increased up to about 5%.

Energy Technology Data Exchange (ETDEWEB)

We present in this paper the computer code BACCHUS, to analyze the thermal-hydraulics in a rod bundle in single or two-phase flow regime. The model is 2-D and uses the porous body approach. The two-phase model is an extension of the classical homogeneous model, and includes a differential non-equilibrium equation. Results are shown for the extension of the boiling region in a 19-pin bundle.

Science.gov (United States)

Some authors have concluded that spiral structures and shocks do not develop if an adiabatic index gamma > 1.16 is adopted in accretion disc modelling, whilst others have claimed that they obtained well defined spirals and shocks adopting a gamma = 1.2 and a M_2/M_1 = 1 stellar mass ratio. In our opinion, it should be possible to develop spiral structures for low compressibility gas accretion discs if the primary component is a black hole. We considered a primary black hole of 8 solar mass and a small secondary component of 0.5 solar mass to favour spiral structures formations and possible spiral shocks via gas compression due to a strong gravitational attraction. We performed two 3D SPH simulations and two 2D SPH simulations and characterized a low compressibility model and a high compressibility model for each couple of simulations. 2D ...

CERN Document Server

We assess the accuracy and relevance of the numerical algorithms based on the principles of Geometrical Optics (GO) and Physical Optics (PO) in the analysis of reduced-size homogeneous dielectric lenses prone to behave as open resonators. As a benchmark solution, we use the Muller boundary integral equations discretized with trigonometric Galerkin scheme that has guaranteed and fast convergence as well as controllable accuracy. The lens cross-section is chosen typical for practical applications, namely an extended hemiellipse whose eccentricity satisfies the GO focusing condition. The analysis concerns homogeneous lenses made of rexolite, fused quartz, and silicon with the size varying between 3 and 20 wavelengths in free space. We consider the 2-D case with both E- and H-polarized plane waves under normal and oblique incidence, and compare characteristics of the near fields.

Science.gov (United States)

We develop a 3-D model to simulate the synthetic aperture radar (SAR) image formation process of an undulated vegetation canopy such as corn grown in fields with large periodic drainage reliefs. We explain how the simulated SAR image of undulated vegetation medium is obtained by the convolution of a 2-D slice of the 3-D simulated SAR system point spread function [(PSF), emulating the SAR beam modeled by a cosine modulated Gaussian], with the 2-D projection of the observed undulated vegetation canopy (modeled with scatterers randomly distributed in 3-D undulated space) followed by the extraction of each look envelope, the summation of looks, and sampling in azimuthal and range directions. Our model is useful to study the parameters involved in the formation and the analysis of SAR images of undulated vegetation medium. Validation of ...

Science.gov (United States)

The vortex-excited dynamics of a uniform pivoted cylinder in uniform and sheared flow was investigated experimentally. The experiments were numerically simulated using a diffusive Van der Pol oscillator model developed by Balasubramanian & Skop recently. Salient features of the experimental investigations and the numerical simulations are presented here. Comparisons between the experimentally recorded and numerically predicted structural response to vortex-excited vibrations, power spectral density measurements of near-wake velocity fluctuations and lock-in ranges are made. A comparison of the numerical predictions and the experimental data reveals good agreement.

Science.gov (United States)

... that range from relatively low frequency strong ground motion modeling [Vidale et al., 1985] to high frequency reverse time migration of seismic ...

UK PubMed Central (United Kingdom)

Many simulation studies in biomedicine are based on a similar sequence of processing steps, starting from images and running through geometric model generation, assignment of tissue properties, numerical...Full Text Available

Energy Technology Data Exchange (ETDEWEB)

The seismic reflection exploration technique which is one of the geophysical methods for oil exploration became effectively to image the subsurface structure with rapid development of computer. As a tool to perform seismic inversion, seismic forward modeling program using ray tracing should be developed. In this study, we have developed the algorithm that is to calculate the travel time of the complex geological structure using ray tracing by subdividing the geologic model into triangular element (finite element) having the constant velocity. We can analytically calculate Jacobian with some information by this current ray tracing. With this Jacobian, we will develop new algorithm which is to obtain geological properties and to image the subsurface. Since the FEM (Finite Element Method) ray tracing we have developed goes well the inverse velocities structure, we can apply the inversion problem to complex geological model. For the convenience of velocity model input data which is needed ...

Science.gov (United States)

We present here a study based on the migration of protoplanets in an accretion disc of a forming star, as the mainly proposed scenario for the formation of planetary systems. Attention is here focused on the mutual interactions between two protoplanets, both embedded in the accretion disc, as a function of the protoplanets masses, their relative positions, the dynamic properties of the accretion disc particles. The study is performed through a 2D SPH code and preliminary results show an oscillation of the distance between the two protoplanets, together with a slow migration of the two planets towards the central star when two Jupiter-like planets are considered. Less correlated behaviour is observed when at least one of the two protoplanets has an Earth-like mass. The role played by the disc particles initial angular momentum is discussed.

Science.gov (United States)

In our previous paper (Lanzafame et al. 2000, PASJ 52, 515) we showed, through 2D SPH simulations, that the stellar mass ratio, M 2 / M 1, of a close binary system (that determines the position and then the initial specific angular momentum at L 1) plays a fundamental role in the formation and development of spiral structures and shock fronts in the radial flow of accretion discs. In that work only a quasi-sonic value of the injection velocity at L 1 was considered. In the present work we also carried out 2D SPH simulations with the aim to investigate the development of such structures, while keeping constant the mass of the compact primary (M 1 = 1 M odot) and the separation between the two components, and assuming as an initial condition of two different supersonic injection velocities at L 1, characterizing two sets of simulations. For ...

Energy Technology Data Exchange (ETDEWEB)

Piping in nuclear power plants is vital to the proper operation and safety of these facilities. To assure safety in the unlikely event of a pipe break, it is necessary to evaluate the consequences from the resulting whipping pipe on neighboring components and structures. Numerical simulations allow for rapid evaluation of these consequences. Before simulations can be accepted, however, the methodology and computer codes must be validated against experimental results. This paper uses a probabilistic approach to validate pipe whip simulations against limited experimental results. Probabilistic analysis software was developed and coupled to existing deterministic finite element software. An example of a whipping pipe impacting against a reinforced concrete slab was simulated. The described probabilistic approach was used to validate the numerical ...

International Nuclear Information System (INIS)

Core-collapse supernovae are among Nature's most energetic events. They mark the end of massive star evolution and pollute the interstellar medium with the life-enabling ashes of thermonuclear burning. Despite their importance for the evolution of galaxies and life in the universe, the details of the core-collapse supernova explosion mechanism remain in the dark and pose a daunting computational challenge. We outline the multi-dimensional, multi-scale, and multi-physics nature of the core-collapse supernova problem and discuss computational strategies and requirements for its solution. Specifically, we highlight the axisymmetric (2D) radiation-MHD code VULCAN/2D and present results obtained from the first full-2D angle-dependent neutrino radiation-hydrodynamics simulations of the post-core-bounce supernova evolution. We then go on to discuss ...

International Nuclear Information System (INIS)

In this paper, finite-difference time-domain method is used to research the transient electromagnetic pulse (EMP) coupling to the cylinder object which is shielded in another shielding cavity with a hole numerically. The different coupling effectiveness of three transient EMP is simulated. The distribution of the electric field in the shielding cavity and the current on the cavity are researched. The shielding effectiveness is evaluated. (authors)

Energy Technology Data Exchange (ETDEWEB)

The forward modelling and the prestack reverse time migration of seismic P-SV wave field was carried out in 2-D models of isotropic and anisotropic media which allow separation of P-SV and SH motion. The P-SV wave field can be described by a system of hyperbolic, first order differential equations in terms of particle velocity and stress. The system of five equations and five unknowns, namely horizontal (U) and vertical (V) velocity components, and three components of stress (T{sub xx}, T-z{sub z} and T{sub xz}) was solved numerically using second order space and forth order time finite differences operators. In order to attenuate numerical dispersion, a staggered grid was used. (author). 48 refs., 5 figs

Energy Technology Data Exchange (ETDEWEB)

Three-dimensional magnetohydrodynamic simulations are executed in a full toroidal geometry to clarify the physical mechanisms of the Internal Reconnection Event (IRE), which is observed in the spherical tokamak experiments. The simulation results reproduce several main properties of IRE. Comparison between the numerical results and experimental observation indicates fairly good agreements regarding nonlinear behavior, such as appearance of localized helical distortion, appearance of characteristic conical shape in the pressure profile during thermal quench, and subsequent appearance of the m=2/n=1 type helical distortion of the torus. (author)

CERN Document Server

Current theories on planetary formation establish that giant planet formation should be contextual to their quick migration towards the central star due to the protoplanets-disc interactions on a timescale of the order of $10^5$ years, for objects of nearly 10 terrestrial masses. Such a timescale should be smaller by an order of magnitude than that of gas accretion onto the protoplanet during the hierarchical growing-up of protoplanets by collisions with other minor objects. These arguments have recently been analysed using N-body and/or fluid-dynamics codes or a mixing of them. In this work, inviscid 2D simulations are performed, using the SPH method, to study the migration of one protoplanet, to evaluate the effectiveness of the accretion disc in the protoplanet dragging towards the central star, as a function of the mass of the planet itself, of disc tangential kinematics and of the presence of a planet ...

Energy Technology Data Exchange (ETDEWEB)

The mainly nonstationary operation of a solar-heated receiver can be simulated with sufficient accuracy only if data about the dynamic behavior are available. For this reason, the dynamic behavior of a solar cavity receiver with parabolic dish collector is investigated. The development of a mathematical simulation considering heat transfer and storage processes is presented and the procedure for a numerical solution is illustrated. The performance of the calculation method is finally demonstrated by simulating the passage of a cloud.

CERN Document Server

Steeghs et al. (1997) have found the first convincing evidence for spiral structure in the accretion disc in IP Pegasi. We perform two kinds of 2D hydrodynamic simulations, a SFS finite volume scheme and a SPH scheme, in the case of mass ratio of 0.5. Both results agree well each other. We construct the Doppler maps and line flux-binary phase relation based on the density distributions. Both of our results agree well with those obtained by the observation.

Energy Technology Data Exchange (ETDEWEB)

This report summarizes preliminary 2-D reactive-transport simulations on the injection, storage and transport of supercritical CO{sub 2} in dolostone formations in the Ordos Basin in China. The purpose of the simulations was to evaluate the role that basin heterogeneity, permeability, CO{sub 2} flux, and geochemical reactions between the carbonate geology and the CO{sub 2} equilibrated brines have on the evolution of porosity and permeability in the storage reservoir. The 2-D simulation of CO{sub 2} injection at 10{sup 3} ton/year corresponds to CO{sub 2} injection at a rate of 3 x 10{sup 5} ton/year in a 3-D, low permeable rock. An average permeability of 10 md was used in the simulation and reflects the upper range of permeability reported for the Ordos Basin Majiagou Group. Transport and distribution of CO{sub 2} ...

International Nuclear Information System (INIS)

This paper presents a new methodology to create realistic 3D microstructures of polycrystals. The virtual microstructures are based on statistical data describing the morphological and crystallographic textures of a sample, obtained from an EBSD analysis. In addition, the methodology can reproduce the observed surface on top of the simulated microstructure. This feature allows finite element calculations on these virtual aggregates to be compared to experimental results of mechanical tests. Such a comparison leads to the identification of the mechanical parameters of constitutive laws, such as critical resolved shear stress and strain hardening, using an optimization algorithm. Two materials were simulated in this study: TiAl and grade 702 zirconium. The first one presents twins inside the microstructure and the second one has an anisotropic texture. Based on 2D simulations, the ...

Energy Technology Data Exchange (ETDEWEB)

The research documented herein was funded by a research contract between the Research Reactors Division (RRD) of Oak Ridge National Laboratory (ORNL) and the University of Tennessee, Knoxville (UTK) Mechanical, Aerospace and Biomedical Engineering Department (MABE). The research was governed by a statement of work (SOW) which clearly defines nine specific tasks. This report is outlined to follow and document the results of each of these nine specific tasks. The primary goal of this phase of the research is to demonstrate, through verification and validation methods, that COMSOL is a viable simulation tool for thermal-hydraulic modeling of the High Flux Isotope Reactor (HFIR) core. A secondary goal of this two-dimensional phase of the research is to establish methodology and data base libraries that are also needed in the full three-dimensional COMSOL simulation to follow. COMSOL version 3.5a was used for all of the models presented throughout ...

Energy Technology Data Exchange (ETDEWEB)

The interaction between a pair of counter-rotating vortices and a lean premixed CH{sub 4}/O{sub 2}/N{sub 2} flame ({Phi} = + 0.55) has been studied by direct numerical simulations using detailed and reduced chemical reaction schemes. Results from the complex chemistry simulation are discussed with respect to earlier experiments and differences in the simulations using detailed and reduces chemistry are investigated. Transient evolutions of the flame surface and the total heat release rate are compared and modifications in the evolution of the local flame structure are displayed. (authors) 22 refs.

Science.gov (United States)

This paper reports that, to obtain better simulation results for a Canada deuterium uranium (CANDU) reactor operation, a new simulation method is developed that uses actual detector readings as a correction factor. Detector readings from a CANDU reactor are used to correct the calculated flux distribution during core calculation iterations. A suitable function is found to describe the relationship between the detector flux and the fluxes of mesh points around the detector. The new simulation method is tested by performing numerical calculations for the Wolsung reactor (a CANDU-600). The results show that the new method predicts the core state more accurately with fewer iterations.

Energy Technology Data Exchange (ETDEWEB)

This paper presents a numerical simulation with the EUROPLEXUS fast dynamics software of a pipe whipping phenomenon occurring in the thermal hydraulic conditions of a loss of coolant accident in a PWR primary circuit. Different physical phenomena take place simultaneously during the rupture and the whipping of the pipe such as plasticity, contact, large displacements, two-phase flow regime and fluid structure interaction. Two kinds of numerical models - a simplified pipeline model and a mixed 1D/3D model - are considered and compared throughout modelling and computation. Numerical results are compared with experimental data belonging to the Aquitaine II test campaign.

CERN Document Server

We study diffusion of Cosmic Rays (CRs) in turbulent magnetic fields using test particle simulations. Electromagnetic fields are produced in direct numerical MHD simulations of turbulence and used as an input for particle tracing, particle feedback on turbulence being ignored. Statistical transport coefficients from the test particle runs are compared with earlier analytical predictions. We find qualitative correspondence between them in various aspects of CR diffusion. In the incompressible case, that we consider in this paper, the dominant scattering mechanism occurs to be the non-resonant mirror interactions with the slow-mode perturbations. Perpendicular transport roughly agrees with being produced by magnetic field wandering.

Energy Technology Data Exchange (ETDEWEB)

We study by means of Quantum Monte Carlo simulations based on the Worm Algorithm the low temperature (down to T = 0.05 K) properties of parahydrogen clusters comprising up to 40 molecules. Three different intermolecular interactions are employed: the Silvera-Goldman, the Buck and the Lennard-Jones potential. Despite important discrepancies observed in the numerical estimates of energy and superfluid fraction, the mechanism by which clusters melt at low T is independent of the particular choice of the potential, whose only effect is to alter the temperature scale.

Science.gov (United States)

Results are given of numerical simulation of three-dimensional pattern of flow of a two-phase steam-water mixture in the house of a PGV-1000 horizontal steam generator obtained using the BAGIRA best-estimate thermohydrodynamic computer codes. The space distributions of velocities and local void fractions in the steam generator housing for different modes of operation of power-generating unit are calculated and compared with available experimental data.

Energy Technology Data Exchange (ETDEWEB)

This dissertation is concerned with the development, experimental diagnostics and mathematical modelling and simulation of polymer electrolyte fuel cells (PEFC). The central themes throughout this thesis are the closely interlinked phenomena of mass and charge transfer. In the face of developing a PEFC system for vehicle propulsion these phenomena are scrutinized on a broad range of relevant scales. Starting from the material related level of the membrane and the gas diffusion layer (GDL) we turn to length scales, where structural features of the cell additionally come into play. These are the scale of flow channels and ribs, the single cell and the cell stack followed by the cell, stack, and system development for an automotive power train. In Chapter 3 selected fundamental material models and properties, respectively, are explored that are crucial for the mathematical modelling and simulation of PEFC, as needed in some succeeding parts of ...

British Library Electronic Table of Contents (United Kingdom)

The ground state of a three-dimensional (3D) rotating trapped superfluid Fermi gas in the BCS-BEC crossover is mapped to finite N v -body vortex states by a simple ansatz. The total vortex energy is measured from the ground-state energy of the system in the absence of the vortices. The vortex state is stable since the vortex potential and rotation energies are attractive while the vortex kinetic energy and interaction between vortices are repulsive. By combining the analytical and numerical works for the minimal vortex energy, the 2D configurations of N v vortices are studied by taking into account of the finite size effects both on xy-plane and on z-direction. The calculated vortex numbers as a function of the interaction strength are appropriate to the renew experimental results by Zwier...

DEFF Research Database (Denmark)

We propose a new physical implementation of spin qubits for quantum information processing, namely defect states in antidot lattices defined in the two-dimensional electron gas (2DEG) at a semiconductor heterostructure. Calculations of the band structure of a periodic antidot lattice are presented. A point defect is created by removing a single antidot, and calculations show that localized states form within the defect, with an energy structure which is robust against thermal dephasing. The exchange coupling between two electrons residing in two tunnel-coupled defect states is calculated numerically. We find results reminiscent of double quantum dot structures, indicating that the suggested structure is a feasible physical implementation of spin qubits.

Energy Technology Data Exchange (ETDEWEB)

Details of the interaction of photons with tissue phantoms are elucidated using Monte Carlo simulations. In particular, photon sampling volumes and photon pathlengths are determined for a variety of scattering and absorption parameters. The Monte Carlo simulations are specifically designed to model light delivery and collection geometries relevant to clinical applications of optical biopsy techniques. The Monte Carlo simulations assume that light is delivered and collected by two, nearly-adjacent optical fibers and take into account the numerical aperture of the fibers as well as reflectance and refraction at interfaces between different media. To determine the validity of the Monte Carlo simulations for modeling the interactions between the photons and the tissue phantom in these geometries, the simulations were compared to measurements of aqueous suspensions ...

Energy Technology Data Exchange (ETDEWEB)

The influence of the velocity profile at the inlet boundary on the simulation of air velocity distribution inside an electrostatic precipitator is presented in this study. Measurements and simulations were performed in a duct and an electrostatic precipitator (ESP). A four-hole cobra probe was used for the measurement of velocity distribution. The flow simulation was performed by using the computational fluid dynamics (CFD) code FLUENT. Numerical calculations for the air flow were carried out by solving the Reynolds-averaged Navier-Stokes equations coupled with the realizable k-{epsilon} turbulence model equations. Simulations were performed with two different velocity profiles at the inlet boundary - one with a uniform (ideal) velocity profile and the other with a non-uniform (real) velocity profile to demonstrate the effect of velocity inlet boundary condition on the flow ...

CERN Document Server

We present a theoretical framework for formal study of systematic effects in Supernovae Type Ia (SN Ia) that utilizes 2-d simulations to implement a form of the deflagration-detonation transition (DDT) explosion scenario. The framework is developed from a randomized initial condition that leads to a sample of simulated SN Ia whose Ni56 masses have a similar average and range to those observed, and have many other modestly realistic features such as the velocity extent of intermediate mass elements. The intended purpose is to enable statistically well-defined studies of both physical and theoretical parameters of the SN Ia explosion simulation. We present here a thorough description of the outcome of the SN Ia explosions produced by our current simulations. A first application of this framework is utilized to study the dependence of the SN Ia on the Ne22 content, ...

Energy Technology Data Exchange (ETDEWEB)

Boiling jet impingement cooling is currently being explored to cool power electronics components. In hybrid vehicles, inverters are used for DC-AC conversion. These inverters involve a number of insulated-gate bipolar transistors (IGBTs), which are used as on/off switches. The heat dissipated in these transistors can result in heat fluxes of up to 200 W/cm{sup 2}, which makes the thermal management problem quite important. In this paper, turbulent jet impingement involving nucleate boiling is explored numerically. The framework for these computations is the CFD code FLUENT. For nucleate boiling, the Eulerian multiphase model is used. The numerical results for boiling water and R113 jets (submerged) are validated against existing experimental data in the literature. Some representative IGBT package simulations that use R134a as the cooling fluid are also presented. (author)

CERN Document Server

We present a new set of 3.5 Post-Newtonian equations in which Newtonian hydrodynamics is coupled to the nonconservative effects of gravitational radiation emission. Our formalism differs in two significant ways from a similar 3.5 Post-Newtonian approach proposed by Blanchet (1993, 1997). Firstly we concentrate only on the radiation-reaction effects produced by a time-varying mass-current quadrupole $S_{ij}$. Secondly, we adopt a gauge in which the radiation-reaction force densities depend on the fourth time derivative of $S_{ij}$, rather than on the fifth, as in Blanchet's approach. This difference makes our formalism particularly well-suited to numerical implementation and could prove useful in performing fully numerical simulations of the recently discovered $r$-mode instability for rotating neutron stars subject to axial perturbations.

International Nuclear Information System (INIS)

A three-dimensional (3D) thermal-hydraulic model and a numerical procedure for the simulation and analysis of a steady-state, as well as transient operation of nuclear power plant components are presented. A two-fluid approach is applied to modelling of two-phase flow. Thermal-hydraulics of a horizontal steam generator in the WWER 1000 nuclear power plant has been simulated at the full load, steady-state operation. A comparison of the numerical results with data measured at the NPP Novovoronjezh shows good agreement. 3D numerical results can be used in plant design or retrofitting, in nuclear power plant operation and safety analysis and as improvement of existing one-dimensional thermal-hydraulics models of the horizontal steam generator which are assessed by system codes used for the nuclear power plant safety analyses. (author)

International Nuclear Information System (INIS)

The model of transient enhanced diffusion of ion-implanted As is formulated and the finite-difference method for numerical solution of the system of equations obtained is developed. The nonuniform distribution of point defects near the interface and more accurate description of arsenic clustering are simultaneously taken into account. Simulation of As diffusion during rapid annealing gives a reasonable agreement with the experimental data. (authors)

British Library Electronic Table of Contents (United Kingdom)

SummaryA meshless numerical model is proposed to investigate shallow-water dam break flows in 1D open channels. The numerical model is to solve the shallow water equations (SWE) based on smoothed particle hydrodynamics (SPH). The concept of slice water particles (SWP) is adopted in the SPH-SWE formulation. The numerical sensitivity analysis is first performed to study the appropriate SWP number and variable smoothing length through dam break flows in an idealized 1D channel with dry/wet beds. Extensive validation by comparison with laboratory and field data is next conducted for four benchmark problems, including dam break flows through a rough flat channel, a rough bumpy channel with various downstream boundary conditions, a nonprismatic channel, and a realistic scale model of the Toce ri...

International Nuclear Information System (INIS)

Full text of publication follows: As the coolant experiences no phase change in the core, SCWRs, unlike LWRs, cannot use design criteria based on the critical heat flux concept. The commonly accepted practice in SCWRs is to specify cladding temperature limits that must be met during transient and accident events. Therefore for the design of the SCWR, it is very important to predict the heat transfer coefficient to the supercritical water coolant with great accuracy. Our recent study focuses on the critical issue of measuring heat transfer to supercritical water at prototypical SCWR conditions and to develop the tools to predict the SCWR thermal behavior. A heat transfer test loop using a surrogate fluids, CO_2, is under construction. The reason of using CO_2 instead of water is that (i) valuable insight of the physical phenomena can be obtained with this fluid, and (ii) some existing facilities already used surrogate fluids, which in general have lower critical pressure and ...

Energy Technology Data Exchange (ETDEWEB)

The understanding of transport processes that occur in naturally fractured geothermal systems is far from being complete. Often, evaluation and numerical simulations of fractured geothermal reservoirs, are carried out by assuming equivalent porous media and homogeneous petrophysical properties within big matrix blocks. The purpose of this paper, is to present a comparison between results obtained from numerical studies of a naturally fractured reservoir treated as a simple porous medium and the simulation of some real aspects of the fractured reservoir. A general conclusion outlines the great practical importance of considering even approximately, the true nature of such systems. Our results show that the homogeneous simplified evaluation of the energy resource in a fractured system, could result in unrealistic estimates of the reservoir capacity to generate electricity.

CERN Document Server

Numerical results on the translocation of long biopolymers through mid-sized and wide pores are presented. The simulations are based on a novel methodology which couples molecular motion to a mesoscopic fluid solvent. Thousands of events of long polymers (up to 8000 monomers) are monitored as they pass through nanopores. Comparison between the different pore sizes shows that wide pores can host a larger number of multiple biopolymer segments, as compared to smaller pores. The simulations provide clear evidence of folding quantization in the translocation process as the biopolymers undertake multi-folded configurations, characterized by a well-defined integer number of folds. Accordingly, the translocation time is no longer represented by a single-exponent power law dependence on the length, as it is the case for single-file translocation through narrow pores. The folding quantization increases with the biopolymer length, ...

Energy Technology Data Exchange (ETDEWEB)

To represent a heterogeneous unsaturated fractured rock by its homogeneous equivalent, Monte Carlo simulations are used to obtain upscaled (effective) flow properties. In this study, we present a numerical procedure for upscaling the van Genuchten parameters of unsaturated fractured rocks by conducting Monte Carlo simulations of the unsaturated flow in a domain under gravity-dominated regime. The simulation domain can be chosen as the scale of block size in the field-scale modeling. The effective conductivity is computed from the steady-state flux at the lower boundary and plotted as a function of the averaging pressure head or saturation over the domain. The scatter plot is then fitted using van Genuchten model and three parameters, i.e., the saturated conductivity K{sub s}, the air-entry parameter {alpha}, the pore-size distribution parameter n, corresponding to this model are considered as the ...

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

CERN Document Server

We present a dynamical spectral model for Large Eddy Simulation of the incompressible magnetohydrodynamic (MHD) equations based on the Eddy Damped Quasi Normal Markovian approximation. This model extends classical spectral Large Eddy Simulations for the Navier-Stokes equations to incorporate general (non Kolmogorovian) spectra as well as eddy noise. We derive the model for MHD and show that introducing a new eddy-damping time for the dynamics of spectral tensors in the absence of equipartition between the velocity and magnetic fields leads to better agreement with direct numerical simulations, an important point for dynamo computations.

Energy Technology Data Exchange (ETDEWEB)

A new scheme to accurately retain kinetic electron effects in particle-in-cell (PIC) simulations for the case of electrostatic drift waves is presented. The splitting scheme, which is based on exact separation between adiabatic and on adiabatic electron responses, is shown to yield more accurate linear growth rates than the standard df scheme. The linear and nonlinear elliptic problems that arise in the splitting scheme are solved using a multi-grid solver. The multi-grid particle-in-cell approach offers an attractive path, both from the physics and numerical points of view, to simulate kinetic electron dynamics in global toroidal plasmas.

Energy Technology Data Exchange (ETDEWEB)

The results of this study show that Monte Carlo simulations of diffusion in homogeneous turbulence can be formulated in terms of the Eulerian space-time velocity autocorrelation function. Numerical results obtained using this approach agree with results obtained by Taylor (1921) using the Lagrangian autocorrelation function. Use was made of the equivalence of the Lagrangian and Eulerian Monte Carlo approaches to derive analytic relations between the Lagrangian integral time scale and the Eulerian integral space and time scales. These analytic results were verified by comparison to Monte Carlo simulations and to other theoretical results. They are in general agreement with many existing theories and semi-empirical relations.

British Library Electronic Table of Contents (United Kingdom)

Piping in nuclear power plants is vital to the proper operation and safety of these facilities. To assure safety in the unlikely event of a pipe break, it is necessary to evaluate the consequences from the resulting whipping pipe on neighboring components and structures. Numerical simulations allow for rapid evaluation of these consequences. Before simulations can be accepted, however, the methodology and computer codes must be validated against experimental results. This paper uses a probabilistic approach to validate pipe whip simulations against limited experimental results. Probabilistic analysis software was developed and coupled to existing deterministic finite element software. An example of a whipping pipe impacting against a reinforced concrete slab was simulated. The described pr...

Energy Technology Data Exchange (ETDEWEB)

A spectral element method embodying Large Eddy Simulation based on Re- Normalization Group theory for simulating Sub Grid Scale viscosity was chosen for this work. This method is embodied in a computer code called NEKTON. NEKTON solves the unsteady, 2D or 3D,incompressible Navier Stokes equations by a spectral element method. The code was later extended to include the variable density and multiple reactive species effects at low Mach numbers, and to compute transport of large particles governed by inertia. Transport of small particles is computed by treating them as trace species. Code computations were performed for a number of test conditions typical of flow past a deep tube bank in a boiler. Results indicate qualitatively correct behavior. Predictions of deposition rates and deposit shape evolution also show correct qualitative behavior. These simulations are the first attempts ...

Science.gov (United States)

We show that accretion disks, both in the subcritical and supercritical accretion rate regime, may exhibit significant amplitude luminosity oscillations. The luminosity time behavior has been obtained by performing a set of time-dependent 2D SPH simulations of accretion disks with different values of ? and accretion rate. An explanation of this luminosity behavior is proposed in terms of limit-cycle thermal instability: the disk oscillates between a radiation pressure dominated configuration (with a high luminosity value) and a gas pressure dominated one (with a low luminosity value). We support this hypothesis showing that the limit-cycle behavior produces a sequence of collapsing and refilling states of the innermost disk region.

Energy Technology Data Exchange (ETDEWEB)

X-ray diffraction enhanced imaging (DEI) is one of X-ray phase-contrast imaging methods, which is applied to inspect internal structures of weakly absorbing low-Z samples. The key problem of the DEI is how to extract phase information which is expressed by refraction-angle images from a series of DEI images measured in different positions of the rocking curve of the analyzer. Three effective extraction methods are presented in this paper: the statistical geometric-optics-approximation method, the maximum refraction-angle method and the Gaussian curve fitting method. They are compared with the existing methods, such as the D. Chapman's geometric optics approximation method and the multiple-images statistical method. A 2D computer simulation experiment is performed to draw comparisons of these methods. The experimental results prove that the above three methods have more precision of refraction-angle values than ...

Energy Technology Data Exchange (ETDEWEB)

This paper describes TENT, a system designed for integrating the various tools and parts of a CFD simulation system (CFD - Computational Fluid Dynamics) and the simulation code itself in one interactive simulation environment. Initially the system was intended to be a TEstbed for Numerical Turbine simulation but it turned out that the approach was general enough to be suitable for more complex simulation tasks in the field of aircraft simulation. The paper will give a brief description of a CFD simulation system and its component interactions, which are the basis for describing the Component Architecture based on CORBA. This includes the control architecture using an event model and data flow considerations. Because most CFD codes use parallelization techniques based on Message Passing, a so called MIMD Object ...

CERN Document Server

We examine the effect of a threshold bias on the power spectrum and the bispectrum in an ensemble of numerical simulations (Gaussian initial perturbations with power law spectra P(k) \\sim k^n, n=+1, 0, -1, -2) and compare our results with theoretical predictions. Our simulations are evolved sufficiently that on the scale where we apply the threshold the rms fluctuation has developed significantly into the nonlinear regime. Thus, predictions based on perturbation theory do not necessarily apply. Nevertheless, we find our results for the power spectrum, biased power simply amplified by a numerical factor, follow predicted trends, far beyond the regime where perturbation theory is expected to be valid. We find that the biased bispectrum continues to follow the so-called hierarchical form, with reduced three-point amplitude Q \\approx 1 in the strongly nonlinear regime, independent of initial spectrum. In ...

Energy Technology Data Exchange (ETDEWEB)

This investigation concerns the nonlinear dynamics of heat transfer from a surface using an upstream eddy promoter. A numerical model is developed for the coupled fluid flow and heat transfer process based on a modified MacCormack scheme. Numerical simulations are carried out to determine the response and heat transfer enhancement due to the promoter. The average heat transfer from a cavity floor is seen to be increased by a factor of approximately five over the unpromoted'' flow. Another interesting feature of the study is the nonlinear viscous flow dynamics from the cylinder-wall interaction which differ significantly from the familiar cylinder-free stream patterns.

Energy Technology Data Exchange (ETDEWEB)

Stabilized formation temperatures were estimated at different depths in 40 wells from the Los Humeros geothermal field, Mexico, using the Horner and the spherical radial flow (SRF) methods. The results showed that the Horner method underestimates formation temperatures, while the SRF method gives temperatures that are closer to the true formation temperatures. This was supported by numerical simulation of a combined circulation and shut-in period in several wells, and results for well H-26 are presented. Numerical reproduction of logged temperature is more feasible if an initial temperature profile based on the SRF method is employed instead of using an initial temperature profile based on the Horner method. (Author)

Energy Technology Data Exchange (ETDEWEB)

The impact of the azimuthal wave refraction in the middle atmosphere on the distribution of gravity wave amplitudes, propagation azimuths, and other wave parameters is investigated using a numerical ray-tracing model of gravity wave propagation through a representative zonal mean reference model of geostrophic winds and temperature in the middle atmosphere. Simulations are first performed with only a single type of gravity wave to help explain some important effects which occur during the refraction process. Then, a multiray simulation is performed which traces a crude spectrum of waves from different altitudes through the atmosphere for every month of a climatological year. The simulated wave climatologies are compared with observations. 108 refs.

CERN Document Server

We examine the observable properties of simulated barred galaxies including radial mass profiles, edge-on structure and kinematics, bar lengths and pattern speed evolution for detailed comparison to real systems. We have run several simulations in which bars are created through inherent instabilities in self-consistent simulations of a realistic disc+halo galaxy model with a disc-dominated, flat rotation curve. These simulations were run at high (N=20M particles) and low (N=500K) resolution to test numerical convergence. We determine the pattern speeds in simulations directly from the phase angle of the bar versus time and the Tremaine-Weinberg method. Fundamental dynamics do not change between the high and low resolution, suggesting that convergence has been reached in this case. We find the higher resolution is needed to simulate ...

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 case without radiation ...

International Nuclear Information System (INIS)

Short communication.

International Nuclear Information System (INIS)

Results of three-dimensional laminar and standard K-#epsilon# turbulent numerical simulations of natural convection cooling of ten cubic aluminum blocks mounted on an insulated plate, facing a shrouding wall, are presented. This geometry is chosen so that comparison with experimental results is possible. The considered problem is of great practical importance because it simulates the case of heated electronic chips, mounted on printed board assemblies, which are frequently encountered in electronic industry applications. The problem is mathematically modeled by the three-dimensional conservation differential equations of mass, momentum, energy and turbulent kinetic energy and dissipation (for the turbulent flow model). IN this paper, these equations are numerically solved by a finite volume method and the laminar and turbulent results are compared to the experimental results obtained with similar ...

CERN Document Server

Dark matter particles form halos that contribute the major part of the mass of galaxy clusters. The formation of these cosmological structures have been investigated both observationally and in numerical simulations, which have confirmed the existence of a universal mass profile. However, the dynamic behaviour of dark matter in halos is not as well understood. We have used observations of 16 equilibrated galaxy clusters to show that the random velocities of dark matter particles are larger on average along the radial direction than along the tangential, and that the magnitude of this velocity anisotropy is radially varying. Our measurement implies that the collective behaviour of dark matter particles is fundamentally different from that of normal particles and the radial variation of the anisotropy velocity agrees with the predictions of numerical simulation.

British Library Electronic Table of Contents (United Kingdom)

In this research, the interactive effect of grain and specimen sizes on the flow stress of sheet metal in microforming is investigated via the tensile test of pure copper and numerical modeling. Models based on different assumptions are proposed to analyze the size effect phenomenon. It is found that the flow stress decreases linearly with the decrease of the ratio of specimen to grain sizes. The grain boundary thickness decreases and its volume fraction increases with the decrease of grain size. The variation of grain boundary thickness is not proportional to the variation of grain size. Furthermore, the fraction of grain boundary increases with the strain and the ratio of specimen to grain sizes. Based on the FE simulation, it is found that the simulated flow stress, which is modeled bas...

British Library Electronic Table of Contents (United Kingdom)

The FVM-LES-acoustic analogy method (FVM-LES-AAM), which is a hybrid prediction technique for the acoustical property computation, is presented and performed in this paper. The FVM-LES-AAM was developed by combining the finite volume method (FVM), the large eddy simulation (LES), and the Ffowcs Williams-Hawkings analogy algorithm (FWH-AA). To predict the acoustical properties of induction cookers, the FVM is used for discretizing the calculation field and building numerical equations, and the LES and FWH-AA are performed for computing the sound sources and predicting the far-field sound, respectively. Using the FVM with the unstructured grids method to discretize the control equation of Navier-Stokes was introduced for illuminating the above numerical simulation procedure. To prove the FVM...

Energy Technology Data Exchange (ETDEWEB)

Functional imaging of the pharynx used to be the domain of cineradiography, CT and ultrafast CT. The development of modern MRI techniques led to new access to functional disorders of the pharynx. The aim of this study was to implement a new MRI technique to examine oropharyngeal obstructive mechanisms in patients with obstructive sleep apnea (OSA). Sixteen patients suffering from OSA and 6 healthy volunteers were examined on a 1.5 T whole-body imager (`Vision`, Siemens, Erlangen Medical Engineering, Germany) using a circular polarized head coil. Imaging was performed with 2D flash sequences in midsagittal and axial planes. Patients and volunteers were asked to breathe normally through the nose and to simulate snoring and the Mueller maneuver during magnetic resonance imaging (MRI). Prior to MRI, all patients underwent an ear, nose and throat (ENT) examination, functional fiberoptic nasopharyngoscopy and polysomnography. A ...

CERN Document Server

We show that accretion disks, both in the subcritical and supercritical accretion rate regime, may exhibit significant amplitude luminosity oscillations. The luminosity time behavior has been obtained by performing a set of time-dependent 2D SPH simulations of accretion disks with different values of alpha and accretion rate. In this study, to avoid any influence of the initial disk configuration, we produced the disks injecting matter from an outer edge far from the central object. The period of oscillations is 2 - 50 s respectively for the two cases, and the variation amplitude of the disc luminosity is 10^38 - 10^39 erg/s. An explanation of this luminosity behavior is proposed in terms of limit cycle instability: the disk oscillates between a radiation pressure dominated configuration (with a high luminosity value) and a gas pressure dominated one (with a low luminosity value). The origin of this instability is the ...

International Nuclear Information System (INIS)

Fluid dynamic analysis of a commercial, counter-flow Ranque-Hilsch Vortex Tube (Rh vt), Ex air 25 s cfm, has been performed in this work both experimentally and numerically; in particular Rh vt cooling power and temperature separation performances have been tested in both direct cooling employment (jet impingement) and indirect cooling employment (supplying cold plates). Experimental techniques, used in this work, revealed several difficulties to produce detailed information about velocity and temperature fields inside the tube and at both the exits. Hence numerical simulation of the flow inside the tube has been conducted using the commercial Cfd code Fluent 6.3.26. Compressible, turbulent, high swirling flow inside Rh vt has been simulated by using both Rans and Les approaches. In particular several turbulence closures have been used in the Rans simulations and results have been ...

International Nuclear Information System (INIS)

Previously, an analytical dose calculation algorithm for MLC-based radiotherapy was developed and commissioned, which includes a detailed model of various MLC effects as a unique feature [1]. The algorithm was originally developed as an independent verification of the treatment planning system's dose calculation and it explicitly modeled spatial and depth dependent MLC effects such as interleaf transmission, the tongue-and-groove effect, rounded leaf ends, MLC scatter, beam hardening, and gradual MLC transmission fall-off with increasing off-axis distance. Originally the algorithm was implemented in Mathematica trademark (Wolfram). To speed up the calculation time and to be able to calculate high resolution 2D dose distributions within a reasonable time frame (<2 s) the algorithm needs to be optimized and to be embedded in a user friendly environment. To achieve this goal, the dose calculation model is implemented in Visual Basic 6.0, which ...

Energy Technology Data Exchange (ETDEWEB)

The Idaho National Engineering Laboratory (INEL) has been conducting geothermal reservoir research and testing sponsored by the US Department of Energy (DOE) since 1983. The INEL research program is primarily aimed at the development of reservoir engineering techniques for fractured geothermal reservoirs. Numerical methods have been developed which allow the simulation of fluid flow and heat transfer in complex fractured reservoirs. Sensitivity studies have illustrated the importance of incorporating the influence of fractures in reservoir simulations. Related efforts include fracture characterization, geochemical reaction kinetics and field testing.

International Nuclear Information System (INIS)

A theoretical model is presented to describe the electromagnetic, heat transfer and fluid flow phenomena within a magnetron plasma torch and in the resultant plume, by using a commercial computational fluid dynamics (CFD) code FLUENT. Specific calculations are presented for a pure argon system (i.e., an argon plasma discharging into an argon environment), operated in a turbulent mode. An important finding of this work is that the external axial magnetic field (AMF) may have a significant effect on the behavior of arc plasma and thus affects the resulting plume. The AMF impels the plasma to retract axially and expand radially. As a result, the plasma intensity distribution on the cross section of torch seems to be more uniform. Numerical results also show that with AMF, the highest plasma temperature decreases and the anode arc root moves upstream significantly, while the current density distribution at the anode is more concentrated with a higher peak value. In ...

Energy Technology Data Exchange (ETDEWEB)

The two-dimensional Eulerian hydrodynamic code 2DE with the shock initiation of heterogeneous explosive burn model called Forest Fire, is used to model numerically the interaction of jets of steel, copper, tantalum, aluminum, and water with steel, water, and explosive targets. The calculated and experimental critical condition for propagating detonation may be described by the Held V/sup 2/d expression (jet velocity squared times the jet diameter). In PBX 9502, jets initiate an overdriven detonation smaller than the critical diameter, which either fails or enlarges to greater than the critical diameter while the overdriven detonation decays to the C-J state. In PBX 9404, the jet initiates a detonation that propagates only if it is maintained by the jet for an interval sufficient to establish a stable curved detonation front. The calculated penetration velocities into explosives, initiated by a low-velocity jet, are ...

CERN Document Server

Without using any moduli, sheaves, stacks, nor any analytic, nor category-type arguments, we exhibit an analogue to Geometric Langlands Theory in an entirely model-independent, non-perturbative,purely smooth topological context in Artin Presentation Theory. A basic initial feature is that AP Theory, as a whole, is already, ab initio, a universal canonical 2D sigma-model, targeting smooth, compact, simply-connected 4-manifolds with a connected boundary, and its topological Planckian quantum starting point, as well as its cone-like, infinitely-generated at each stage, graded group of homology-preserving, but topology-changing transitions/interactions, exhibit the most general qualitative S-duality. We first point out the numerous mathematically rigorous, model-free, (i.e., intrinsic), topological AP analogues with the heuristic Kapustin-Witten version of Geometric Langlands theory, as well as the crucial differences between ...

International Nuclear Information System (INIS)

The mechanism proposed by Kormendy (1984) for the formation of counterrotating cores in elliptical galaxies is investigated using self-consistent numerical simulations of mergers between a high- and a low-luminosity elliptical galaxies. The conditions for a counterrotation to appear are determined, observational properties of the remnants are described, and the evolution of the structural and kinematic parameters of the larger galaxy is analyzed. It is shown that a counterrotation results only when the merging orbits are retrograde, due to a large change in the secondary spin during the merger. 36 refs.

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.

Energy Technology Data Exchange (ETDEWEB)

This paper suggests how nonlinear periodic optimal control of a pumped storage plant can be realized. The control problem consists in maximization of the plant benefits over an operational period. An optimal control law is proposed, yielding a bang-off-bang mode of operation. An algorithm for numerical solution of the problem was developed, and its effectiveness was demonstrated by simulation experiments.

Energy Technology Data Exchange (ETDEWEB)

In this letter we report a clear and unambiguous observation of the out-of-plane quadrupole magnetic field suggested by numerical simulations in the reconnecting current sheet in the Magnetic Reconnection Experiment (MRX). Measurements show that the Hall effect is large in collisionless regime and becomes small as the collisionality increases, indicating that the Hall effect plays an important role in collisionless reconnection.

International Nuclear Information System (INIS)

A distinguishing design feature of CANDU nuclear reactors is the use of horizontal fuel channels housed in a horizontal vessel called a calandria, which is made of stainless steel 304L. Each channel consists of a Zr-2.5%Nb alloy pressure tube and an externally concentric Zr-2 calandria tube. The calandria tubes are joined to the end plates (tubesheets) of the calandria vessel by joints formed by roller expansion. The bores in the tubesheets are grooved. Roller expanded joints provide a cost effective means of joining dissimilar materials, require minimal space and no maintenance. The quality of these roller expanded joints is important from a sealing, strength and stress corrosion point of view. The roller expansion process consists of expanding the calandria tubes to deform them plastically against the bores and into the grooves of the tubesheets. Therefore, understanding the effect of the number, geometry and the pitch of the grooves on the quality of a roller expanded joint is very ...

Energy Technology Data Exchange (ETDEWEB)

There is a crescent interest in the scientific community in the study of natural circulation phenomenon. New generation of compact nuclear reactors uses the natural circulation of the fluid as a system of cooling and of residual heat removal in case of accident or shutdown. The objective of this paper is to present a study through the comparison of experimental data and numerical simulation for the natural circulation phenomenon in one and two-phase flow regime. An experimental circuit built with glass tubes is used for the experiments. Thus, it allows the thermal hydraulic phenomena visualization. There is an electric heater as the heat source, a heat exchanger as the heat sink and an expansion tank to accommodate fluid density excursions. The circuit instrumentation consists of thermocouples and pressure meters to better keep track of the flow and heat transfer phenomena. Instrumentation data acquisition is performed through a computer ...

International Nuclear Information System (INIS)

The influence that electromagnetic fields have on mechanical and electronic equipment was discussed. A three dimensional numerical simulation model of electric fields around high voltage power transmission lines was described. 4 refs., 1 tab., 5 figs.

Energy Technology Data Exchange (ETDEWEB)

Geologic sequestration is currently being practiced and scientifically evaluated as a critical component in a broad strategy, comprising new practices and technologies, for mitigating global climate change due to anthropogenic emissions of CO2. Demonstrating that geologic sequestration of CO2 is safe and effective, and gaining public acceptance of sequestration technologies are critically important in meeting these global climate change challenges. Monitored field-scale demonstrations of geologic sequestration of carbon dioxide will contribute greatly toward growing trust and confidence in the technology; however, pilot demonstrations ultimately will not be the norm for new geological sequestration deployments. Instead, scientists, engineers, regulators, and ultimately the public will rely on numerical simulations to predict the performance of geologic repositories for carbon dioxide sequestration. The U.S. Department of Energy (DOE), through ...

Energy Technology Data Exchange (ETDEWEB)

This contribution is an informal essay based on a talk delivered at the Institute for Mathematics and its Applications (IMA) in Minneapolis, under the summer program in molecular biology, July 18-22, 1994. I exclude many technical details, which can be found elsewhere, and instead focus on the basic ideas of molecular dynamics simulations, with the goal of conveying to students and non-specialists the key concepts of the theory and practice of large-scale simulations. Following a description of the basic idea in molecular dynamics, I discuss some of the practical details involved in simulations of large biological molecules, the numerical timestep problem, and approaches to this problem based on implicit-integration techniques. I end with a perspective of open challenges in the field and directions for future research. 79 refs., 6 figs., 1 tab.

Energy Technology Data Exchange (ETDEWEB)

Incomplete convergence in numerical simulation such as computational physics simulations and/or Monte Carlo simulations can enter into the calculation of the objective function in an optimization problem, producing noise, bias, and topo- graphical inaccuracy in the objective function. These affect accuracy and convergence rate in the optimization problem. This paper is concerned with global searching of a diverse parameter space, graduating to accelerated local convergence to a (hopefully) global optimum, in a framework that acknowledges convergence uncertainty and manages model resolu- tion to efficiently reduce uncertainty in the final optimum. In its own right, the global-to-local optimization engine employed here (devised for noise tolerance) performs better than other classical and contemporary optimization approaches tried individually and in combination on the "industrial" test problem to be ...

International Nuclear Information System (INIS)

The energy deposited by the passage of a single, energetic, heavy-ion through a semiconductor produces dense electron-hole (eh) pair concentrations near the ion trajectory. The size, shape, and charge density of an ion track represent critical parameters for many models of single event phenomena. The authors describe the design and uses of possible semiconductor test structures for measuring the initial radial distribution of charge and subsequent charge transport in a high energy, heavy-ion track. Numerical simulations show how the test structure can resolve different radial distributions of charge within an ion track. The test structure simulations also show the importance of accurately representing ion track structure in single event effects simulations.

Energy Technology Data Exchange (ETDEWEB)

The aim of this study concerns the use of numerical methods for the resolution of the Reynolds Averaged Navier Stokes equations adapted to the simulation of the cooling of the trailing edge of a stator in a high pressure turbine. These methods, based on the elsA solver developed at ONERA, use a four steps Runge Kutta time discretization scheme and a Jameson centered space discretization scheme. The scheme is applied through a finite volume approach on control volume centered on the cells of a multi-block structured mesh. Turbulence is simulated either through the algebraic Michel model, or through the one-transport-equation Spalart-Allmaras model, or through the two-transport-equations k 1, k {omega} and k {epsilon} models, and through ASM model. A simulation of the flow in a bidimensional stator, without cooling, is carried out. The cooling, which is realized with trailing edge slots, is then ...

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 turbulence are modeled by ...

International Nuclear Information System (INIS)

This thesis was carried out at Forsmark Nuclear Power Plant. The power plant in Forsmark consists of three boiling water reactors (BWR) which produce about 17% of Swedish electricity. In a BWR the nuclear reactions are used to boil water inside the reactor vessel. The water works both as a coolant and as a moderator and the resulting steam is used directly to run the turbines. A problem when running a BWR at low flow conditions is the density wave oscillations that might occur to the water flow inside the fuel assemblies. These oscillations arise due to the connection between power and flow rate in a heated channel with two-phase flow. In order to improve the stability performance of the channel an orifice plate is placed at the inlet of each fuel assembly. Today these orifice plates have sharp edges and a constant resistance coefficient. Experimental work has been done with progressive orifices, the edge of which is half-oval in shape. The advantage of progressive orifices is the ...

CERN Document Server

Numerical simulations predict that a large number of baryons reside in intergalactic space at temperatures between 10^5-10^7 K. Highly-ionized metals, such as O VII and O VIII, are good tracers of this ``warm-hot intergalactic medium'', or WHIM. For collisionally-ionized gas, the ionization fraction of each ion peaks at some particular temperature (``peak temperatures''), so different ions can therefore trace the IGM at different temperatures. We performed a hydrodynamic simulation to study the metal distributions in the IGM. We then draw random lines-of-sight across the simulated region and synthesize resonance absorption line spectra in a similar way to simulating the Ly-alpha forest. By studying the distribution functions of H- and He-like O, Si and Fe in a collisionally-ionized IGM and comparing with semi-analytic results based on the Press-Schechter formalism, we find: (1) ions ...

UK PubMed Central (United Kingdom)

BackgroundRecessive mutations in guanylate cyclase-1 (Gucy2d) are associated with severe, early onset Leber congenital amaurosis-1(LCA1). Gucy2d...Full Text Available

UK PubMed Central (United Kingdom)

BackgroundIn current comparative proteomics studies, the large number of images generated by 2D gels is currently compared using spot matching algorithms. Unfortunately, differences...Full Text Available

CERN Document Server

The Sunyaev-Zeldovich (SZ) effect gives a measure of the thermal energy and electron pressure in groups and clusters of galaxies. In the near future SZ surveys will map hundreds of systems, shedding light on the pressure distribution in the systems. The thermal energy is related to the total mass of a system of galaxies, but it is only a projection that is observed through the SZ effect. A model for the 3D distribution of pressure is needed to link the SZ signal to the total mass of the system. In this work we construct an empirical model for the 2D and 3D SZ profile, and compare it to a set of realistic high resolution SPH simulations of galaxy clusters and groups, and to a stacked SZ profile for massive clusters derived from WMAP data. Furthermore, we combine observed temperature profiles with dark matter potentials to yield an additional constraint, under the assumption of hydrostatic equilibrium. We find a very tight ...

International Nuclear Information System (INIS)

IMRT optimization requires a fast yet relatively accurate algorithm to calculate the iteration dose with small memory demand. In this paper, we present a dose calculation algorithm that approaches these goals. By decomposing the infinitesimal pencil beam (IPB) kernel into the central axis (CAX) component and lateral spread function (LSF) and taking the beam's eye view (BEV), we established a non-voxel and non-beamlet-based dose calculation formula. Both LSF and CAX are determined by a commissioning procedure using the collapsed-cone convolution/superposition (CCCS) method as the standard dose engine. The proposed dose calculation involves a 2D convolution of a fluence map with LSF followed by ray tracing based on the CAX lookup table with radiological distance and divergence correction, resulting in complexity of O(N"3) both spatially and temporally. This simple algorithm is orders of magnitude faster than the CCCS method. Without ...

International Nuclear Information System (INIS)

The present study sought a biochemical explanation for retarded brain development in the heterozygous offspring of the phenylketonuric (PKU) mother. Two rat models of simulated maternal PKU, one induced by p-chloropheylalanine and phenylalanine and the other by phenylacetate, were employed in this investigation. Maternal PKU had no influence on cerebral concentrations of DNA, protein, and cholesterol, which were normal in the 2 d old pup. However, there was a noticeable disruption of the normal ganglioside pattern and a significant reduction of sialoglycoproteins. Concomitant with a delayed drop in the gangliosides Q/sub 1b/ and D_3, was a slower rise in M_1 and D/sub 1a/. At least 66% of sialoglycoproteins located on SDS-PAGE gel chromatograms, by radioactivity incorporated in vivo from radiolabeled N-acetylmannosamine and by ("3H) sialic acid released by neuraminidase from periodate-("3H) borohydride labeled ...

International Nuclear Information System (INIS)

Nuclear Engineering Laboratory of the Technical Research Centre of Finland has widely used RELAP5/MOD2 and -MOD3 codes to simulate horizontal steam generators. Several models have been developed and successfully used in the VVER-safety analysis. Nevertheless, the models developed have included only rather few nodes in the steam generator secondary side. The secondary side has normally been divided into about 10 to 15 nodes. Since the secondary side at the steam generators of VVER-440 type reactors consists of a rather large water pool, these models were only roughly capable to predict secondary side flows. The paper describes an attempt to use RELAP5/MOD3 code to predict secondary side flows in a steam generator of a VVER-440 reactor. A 2D/3D model has been developed using RELAP5/MOD3 codes cross-flow junctions. The model includes 90 volumes on the steam generator secondary side. The model has been used to calculate steady ...

International Nuclear Information System (INIS)

A detailed treatment is introduced to measure the dynamic stability of the relativistic electrons in a self-amplified spontaneous emission free-electron laser (FEL) system, which includes the numerical approach of the Kolmogorov entropy (entropy-like quantity), the general equations of motion for a charged particle and the method of monitoring the simulation accuracy. Numerical experiments reveal a new phenomenon that there exists the possibility of the transition from chaotic to non-chaotic phase-space trajectories of the strongly relativistic electrons due to the effect of their self-fields. The adiabatic magnetic field of a one-dimensional wiggler may have a slight influence on the electron transportation in the absence of the FEL fields, but substantially affects the dynamic stability of the electrons in the process of the FEL interaction. Moreover, the laser fields diminish the dynamic stability of the electrons as the ...

CERN Document Server

The Fourier transform of cosmological density perturbations can be represented in terms of amplitudes and phases for each Fourier mode. We investigate the phase evolution of these modes using a mixture of analytical and numerical techniques. Using a toy model of one-dimensional perturbations evolving under the Zel'dovich approximation as an initial motivation, we develop a statistic that quantifies the information content of the distribution of phases. Using numerical simulations beginning with more realistic Gaussian random-phase initial conditions, we show that the information content of the phases grows from zero in the initial conditions, first slowly and then rapidly when structures become non-linear. This growth of phase information can be expressed in terms of an effective entropy: Gaussian initial conditions are a maximum entropy realisation of the initial power spectrum, gravitational evolution decreases the phase ...

Energy Technology Data Exchange (ETDEWEB)

Combined gas and steam turbine processes based on direct coal firing show a high thermal efficiency. At RWTH Aachen, University of Technology, an experimental test furnace has been built to investigate the pressurized pulverized coal combustion (PPCC). The PPCC-facility has been constructed as a slag tap furnace. Particles hitting the walls at temperatures above the melting point cause slagging depositions and create a film flowing down the reactor walls. As a part of the PPCC-program different mathematical models have been developed and implemented into the CFD-code FLUENT to predict the behavior of slag films at the furnace walls. Numerical strategies and the mathematical models used are described in detail. 12 refs., 9 figs.

Energy Technology Data Exchange (ETDEWEB)

The scope of the present article is two-fold. Firstly, to conduct an experiment to provide the temperature-time history of the cooling of a hot ball bearing in quiescent ambient air. Secondly, to predict the temporal variation of the bearing under the hypothesis of natural convection, radiation or natural convection coexists with radiation for a non-vanishing total hemispherical emissivity of the surface of the bearing. Numerical solutions of the three governing nonlinear lumped heat equations were carried out with a Runge-Kutta-Fehlberg (RKF45) algorithm accounting for automatic step size control. The experimental data was obtained with chrome steel ball bearings of diameter 0.953 cm (7/16 in) heated in an electric oven to a pre-set temperature. The heated bearing was exposed later to ambient air at atmospheric temperature and pressure. (orig.)

British Library Electronic Table of Contents (United Kingdom)

A test bench was designed and assembled to carry out impact tests on samples and components. The system allows simple and rapid adjustment of the test parameters, such as the shot size and air pressure, with good repeatability of the results. Tests on steel and light alloys were carried out under both as-produced condition and on thin hard-coated samples. Significant reductions in dimple dimensions were seen after coating. FE models simulating the experiments overestimated the dimple depths, although the parameter trend was satisfactorily captured. The residual stresses from coating and shot peening determined numerically are believed to have been proven effective against fatigue.

DEFF Research Database (Denmark)

We present a theoretical study of the resonant interaction between dynamical localized states (discrete breathers) and linear electromagnetic excitations (EE's) in Josephson junction ladders. By making use of direct numerical simulations we find that such an interaction manifests itself by resonant steps and various sharp switchings (voltage jumps) in the current-voltage characteristics. Moreover, the power of ac oscillations away from the breather center (the breather tail) displays singularities as the externally applied dc bias decreases. All these features may be mapped to the spectrum of EE's that has been derived analytically and numerically. Using an improved analysis of the breather tail, a spectroscopy of the EE's is developed, The nature of breather instability driven by localized EE's is established.

Energy Technology Data Exchange (ETDEWEB)

The bulk removal of CO/sub 2/ from industrial gases is a frequently applied process in the field of gas-treating. Usually chemical solvents are used in order to increase both the solvent-loading capacity and the absorption rate. These solvents often are aqueous solutions of alkanolamines. The loaded-liquid is regenerated at high temperatures with the aid steam. In the present study absorption rates are experimentally determined in various mixtures of amines over a wide range of concentrations. Simultaneously, a model was developed for the description of mass transfer accompanied by parallel reversible chemical reactions. The equations of this mass transfer model could not be solved analytically and therefore numerical techniques were used to obtain an exact description of this process. This model also includes the asymptotic situations of homogeneous catalysts and shuttle mechanism for which in the literature several approximate solutions have been presented. The ...

CERN Document Server

Cosmological hydrodynamical simulations as well as observations indicate that spiral galaxies are comprised of five different components: dark matter halo, stellar disc, stellar bulge, gaseous disc and gaseous halo. While the first four components have been extensively considered in numerical simulations of binary galaxy mergers, the effect of a hot gaseous halo has usually been neglected even though it can contain up to 80% of the total gas within the galaxy virial radius. We present a series of hydrodynamic simulations of major mergers of disc galaxies, that for the first time include a diffuse, rotating, hot gaseous halo. Through cooling and accretion, the hot halo can dissipate and refuel the cold gas disc before and after a merger. This cold gas can subsequently form stars, thus impacting the morphology and kinematics of the remnant. Simulations of isolated systems with total ...

CERN Document Server

We investigate weak lensing by large-scale structure using ray tracing through N-body simulations. Photon trajectories are followed through high resolution simulations of structure formation to make simulated maps of shear and convergence on the sky. Tests with varying numerical parameters are used to calibrate the accuracy of computed lensing statistics on angular scales from about 1 arcminute to a few degrees. Various aspects of the weak lensing approximation are also tested. For fields a few degrees on a side the shear power spectrum is almost entirely in the nonlinear regime and agrees well with nonlinear analytical predictions. Sampling fluctuations in power spectrum estimates are investigated by comparing several ray tracing realizations of a given model. For survey areas smaller than a degree on a side the main source of scatter is nonlinear coupling to modes larger than the survey. We develop a ...

British Library Electronic Table of Contents (United Kingdom)

A parameterized surface wave-induced vertical mixing (B"v) is incorporated into the Modular Ocean Model version 4 (MOM4). A comparison of the two numerical experiments with and without B"v shows that B"v can significantly improve the upper-ocean (20-100m) simulation in summer. The simulated upper-ocean temperature errors are reduced in summer due to the surface wave-induced vertical diffusive heat flux. The non-breaking-wave-induced vertical mixing can increase the probability of the simulated SST biases between -1^oC and 1^oC from 64% to 76% in the Southern Hemisphere (60^oS-10^oS) in January, and from 66% to 75% in the Northern Hemisphere (10^oN-60^oN) in July. The averaged mixed layer depth (MLD) simulated by the MOM4 without B"v is 14.4m shallower than the observations in 10^oS-60^oS i...

International Nuclear Information System (INIS)

In this paper, we calculated the radioactive concentration distribution of radioactive waste water, the temperature distribution of drained cooling water and the effect of implement from the Daya Bay Nuclear Power Plant on nearby waters range, discussed and analysed some problems of computational results and computation with Alternating Direction Implicit Method (ADI). The contents of the article included: the establishment of two-dimension tidal current equation, radioactive waste water pollutant dispersion equation and cooling water heat convection diffusion equation, the numerical difference calculation model of tidal current field, concentration field as well as temperature field, effect impingement with ADI method, numerical calculation results. The result of research showed that: when the Daya Bay Nuclear Power Plant is on normal operation and after the low level radioactive waste water and low temperature cooling water have been drained ...

CERN Document Server

A Langmuir film is a molecularly thin film on the surface of a fluid; we study the evolution of a Langmuir film with two co-existing fluid phases driven by an inter-phase line tension and damped by the viscous drag of the underlying subfluid. Experimentally, we study an 8CB Langmuir film via digitally-imaged Brewster Angle Microscopy (BAM) in a four-roll mill setup which applies a transient strain and images the response. When a compact domain is stretched by the imposed strain, it first assumes a bola shape with two tear-drop shaped reservoirs connected by a thin tether which then slowly relaxes to a circular domain which minimizes the interfacial energy of the system. We process the digital images of the experiment to extract the domain shapes. We then use one of these shapes as an initial condition for the numerical solution of a boundary-integral model of the underlying hydrodynamics and compare the subsequent images of the experiment to the ...

International Nuclear Information System (INIS)

The RELAP5 code has been developed for best-estimate transient simulation of light water reactor coolant systems during postulated accidents. The code models the coupled behavior of the reactor coolant system and the core for loss-of-coolant accidents and operational transients such as anticipated transient without scram, loss of offsite power, loss of feedwater, and loss of flow. A generic modeling approach is used that permits simulating a variety of thermal hydraulic systems. Control system and secondary system components are included to permit modeling of plant controls, turbines, condensers, and secondary feedwater systems. RELAP5/MOD3 code documentation is divided into seven volumes: Volume I presents modeling theory and associated numerical schemes; Volume II details instructions for code application and input data preparation; Volume III presents the results of developmental assessment cases that demonstrate and ...

Energy Technology Data Exchange (ETDEWEB)

A numerical model was developed to simulate the combined effects of heat transfer, magnetite oxidation, and NO{sub x} formation in a grate-kiln furnace for taconite pellet induration. Heat transfer from the flame in the kiln was described by the net radiation method. The shrinking core model was used to account for magnetite oxidation on the grate. A novel approach to oxidation of tumbling pellets in a kiln was derived. The Zeldovich mechanism was used to predict thermal NO generation. Temperature fluctuations in the gas streams were estimated with a clipped Gaussian probability density function. The thermal energy balances and mass balances resulted in coupled systems of first-order differential equations, which were solved numerically. The model is capable of predicting NO production and magnetite oxidation in agreement with observation of plant performance. Although the design of the grate-kiln system is for efficient ...

Energy Technology Data Exchange (ETDEWEB)

The results of an experimental investigation of the performance of a small bore engine with a swirl chamber when operating as a dual fuel engine with commercial methane as the gaseous fuel were presented in this paper. The experiment involved using a 3-dimensional computational fluid dynamics model to predict the performance of the engine. A detailed chemical kinetics for the gaseous fuel component, consisting primarily of methane and a reduced detailed chemical kinetics for the diesel fuel while considering the turbulent combustion processes an associated performance of a dual fuel engine with a swirl chamber were incorporated in the simulation. The study experimentally and numerically investigated the effects of changes in the quantities of the liquid fuel pilot and gaseous fuels on the combustion processes, engine performance, cyclic variations, and emissions. The paper discussed the experimental approach and results. It also discussed the ...

CERN Document Server

This work investigates the impingement of a liquid microdroplet onto a glass substrate at different temperatures. A finite-element model is applied to simulate the transient fluid dynamics and heat transfer during the process. Results for impingement under both isothermal and non-isothermal conditions are presented for four liquids: isopropanol, water, dielectric fluid (FC-72) and eutectic tin-lead solder (63Sn-37Pb). The objective of the work is to select liquids for a combined numerical and experimental study involving a high resolution, laser-based interfacial temperature measurement to measure interfacial heat transfer during microdroplet deposition. Applications include spray cooling, micro-manufacturing and coating processes, and electronics packaging. The initial droplet diameter and impact velocity are 80 {\\mu}m and 5 m/s, respectively. For isothermal impact, our simulations with water and isopropanol show very ...

Energy Technology Data Exchange (ETDEWEB)

Although incompressible fluid flows can be regarded as a particular case of a general problem, numerical methods and the mathematical formulation aimed to solve compressible and incompressible flows have their own peculiarities, in such a way, that it is generally not possible to attain both regimes with a single approach. In this work, we start from a typically compressible formulation, slightly modified to make use of pressure variables and, through augmenting the stabilising parameters, we end up with a simplified model which is able to deal with a wide range of flow regimes, from supersonic to low speed gas flows. The resulting methodology is flexible enough to allow for the simulation of liquid flows as well. Examples using conservative and pressure variables are shown and the results are compared to those published in the literature, in order to validate the method. (author)

British Library Electronic Table of Contents (United Kingdom)

This work aims to contribute to the understanding of the influence of the ionospheric layer height (ILH) on the thin layer ionospheric model (TLIM) used to retrieve ionospheric information from the GNSS observations. Particular attention is paid to the errors caused on the estimation of the vertical total electron content (vTEC) and the GNSS satellites and receivers inter-frequency biases (IFB), by the use of an inappropriate ILH. The work relies upon numerical simulations performed with an empirical model of the Earth?s ionosphere: the model is used to create realistic but controlled ionospheric scenarios and the errors are evaluated after recovering those scenarios with the TLIM. The error assessment is performed in the Central and the northern part of the South American continents, a re...

CERN Document Server

After two decades of direct dynamical simulation of large-scale structure in the universe, it is safe to say the subject is now mature. Still, there are parts of the problem that are less well developed than others. In general, the collisionless dynamics of the dark matter component is better understood than the collisional gas dynamics of the baryonic component. In situations where the gas dynamics is relatively simple, such as the Lyman-$\\alpha$ forest and the intracluster medium in X-ray clusters, our ability to reproduce observational data has evolved rapidly, and the interpretive and predictive power of such experiments should now be taken seriously. A comparison of twelve gas dynamic codes to the problem of forming a single X-ray cluster shows that numerical inaccuracies are modest (typically below ten percent), leaving missing physics as the main source for large systematic differences between theory and observation. Galaxy formation, ...

International Nuclear Information System (INIS)

In this paper, a time-domain equivalent circuit method is applied to solve dispersion of coupled-cavity travelling-wave tube (CCTWT). First, the time-domain circuit equations of CCTWT coupled-cavity chain are deduced from the equivalent circuit model. Then, the equations are solved numerically by fourth-order Runge-Kutta method and a program CTTDCP is developed using MATLAB. Last, a L-band CCTWT is calculated using CTTDCP and the cavity pass-band of this tube is computed to be 1.08-1.48 GHz, which is consistent with the experimental results and the simulation results of electromagnetic code and demonstrates the validity of the time-domain equivalent circuit method. In addition, a new design method which uses the equivalent circuit method and electromagnetic simulation together to optimize the cold cavity characteristics of CCTWT is proposed. (authors)

Energy Technology Data Exchange (ETDEWEB)

Recently, we showed (Yu, Phys. Rev. E 58 (1998) 4991) that the widely used simulation code TDA3D, even though a single-frequency code, can be used to determine the power spectrum in the SASE process with excellent approximation in the exponential growth regime. In this paper, we apply this method to the BNL Cornell Wiggler A SASE experiment as an example. When the gain is not very high, there are many modes in the radiation, which seems to make the analytical calculation very difficult. However, we show that the increment of the radiation due to SASE over the spontaneous radiation can be expanded in terms of guided modes with rapid convergence. Thus when the spontaneous radiation is subtracted from the SASE power during the calculation, there is a good agreement between the analytical theory and the numerical simulation.

CERN Document Server

We investigate the influence of the turbulence forcing on the mass distributions of gravitationally unstable cores by postprocessing data from simulations of non-selfgravitating isothermal supersonic turbulence with varying resolution. In one set of simulations solenoidal forcing is applied, while the second set uses purely compressive forcing to excite turbulent motions. From the resulting density field, we compute the mass distribution of gravitationally unstable cores by means of a clump-finding algorithm. Using the time-averaged probability density functions of the mass density, semi-analytic mass distributions are calculated from analytical theories. We apply stability criteria that are based on the Bonnor-Ebert mass resulting from the thermal pressure and from the sum of thermal and turbulent pressure. Although there are uncertainties in the application of the clump-finding algorithm, we find systematic differences in the mass ...

British Library Electronic Table of Contents (United Kingdom)

In this paper linear and nonlinear models of spacecraft attitude dynamics equations and gravity gradient moments are investigated. In addition, effects of gravity gradient moments on attitude dynamics of the satellite are studied. The purpose of this paper is to present a comparison between nonlinear and linear models of spacecraft attitude dynamics and gravity gradient moments in order to determine divergence of linear approximation from the nonlinear model. Simulation results indicate that designer of spacecraft attitude control subsystem should be meticulous in applying linear approximation of equations especially in low earth orbits. Consequently, finding an upper bound for small angle to keep the linear model valid and precise enough would be a vital part of using linear approximation...

British Library Electronic Table of Contents (United Kingdom)

A suite of numerical experiments is conducted to extend our previous studies to explore the effects of linear, ambient wind shear on simulated mammatus-like clouds. Increasing values of unidirectional shear result in banded cloud bases, and for the strongest shears roll-type quasi-two-dimensional convection results rather than three-dimensional lobe morphology. Two observed soundings were used, and the sounding with a drier, shallower sub-cloud layer appeared to be most strongly affected by the presence of the ambient shear. The drier sounding also had less snow in the mammatus lobes owing to sublimation. Copyright Copyright 2009 Royal Meteorological Society

International Nuclear Information System (INIS)

There has been increased theoretical and experimental research interest in autonomous mobile robots exhibiting cooperative behaviour. This paper provides consistent quantitative measures of organizational degree of a two-dimensional environment. We proved, by the way of numerical simulations, that the theoretically derived values of the feature are reliable measures of aggregation degree. The slope of the feature's dependence on memory radius leads to an optimization criterion for stochastic functional self-organization. We also described the intellectual heritages that have guided our research, as well as possible future developments. (author)

International Nuclear Information System (INIS)

Axial offset anomaly (AOA) is caused by the deposition of crud on the fuel cladding of a PWR. When significant levels of crud build up on the cladding, boron can accumulate in the pores of the crud as a concentrated solution or solid phase, and cause the flux depression. Numerous studies have been conducted on the primary water chemistry to reduce the amount of crud in the primary circuit to avoid radioactivity buildup and unexpected power transition in the plant. However, experiments on the crud are restricted in the laboratory because the crud is a highly radioactive material. The objective of this study is to develop a test method for simulating the deposition of crud in a nuclear power plant

Energy Technology Data Exchange (ETDEWEB)

This report details an investigation into the efficacy of two approaches to solving the radiation diffusion equation within a radiation hydrodynamic simulation. Because leading-edge scientific computing platforms have evolved from large single-node vector processors to parallel aggregates containing tens to thousands of individual CPU's, the ability of an algorithm to maintain high compute efficiency when distributed over a large array of nodes is critically important. The viability of an algorithm thus hinges upon the tripartite question of numerical accuracy, total time to solution, and parallel efficiency.

International Nuclear Information System (INIS)

Hung and Quyen's model (analysis of strategic interactions between players in the game theoretic framework) is first recapitulated. Solutions to the embargo game with the seller acting as the Stackelberg leader in both periods. Then the timing coordination issue is discussed and perfect equilibrium outcomes under different structures of leader-follower in the game are compared. Numerical simulations show that the structure corresponding to alternated leadership from one period to the other yields the perfect equilibrium outcome that is Pareto improving with respect to the structure where the buyer is first mover in both periods and Pareto dominant for some specific values of the parameters embedded in the game.

DEFF Research Database (Denmark)

Nonlinear distortion added by the loudspeaker in a hearing aid lowers the signal-to-noise ratio and may degrade the hearing aid user's ability to understand speech. The balancedarmature- type loudspeakers, predominantly used in hearing aids, are inherently nonlinear devices, as any displacement of the loudspeaker diaphragm inevitably changes the magnetic and electrical characteristics of the loudspeaker. A numerical time-domain model capable of describing these nonlinearities is presented. By simulation it is demonstrated how the output distortion could potentially be reduced significantly through careful design of the mechanical properties of the armature.

Energy Technology Data Exchange (ETDEWEB)

A numerical model has been developed to simulate the charging and discharge behaviour of an electrically rechargeable alkaline zinc-air battery. Further a galvanostatic experiment including three charge/discharge cycles has been performed. The cell voltages, the Zn electrode potentials versus a Zn reference, and the O{sub 2} electrode potentials versus a Zn reference calculated with the model are in fairly good agreement with the corresponding experimental data. The model is expected to be useful for zinc-air battery design and for analysis of experimental data. (author)

Science.gov (United States)

The electronic transitions and photodissociation of the bromine molecule were studied in the visible-near UV continuum using dynamic simulation. The molar extinction coefficients in this study were obtained in numerical calculations. The quantum yields of the spin-orbit Br*(2 P 1/2) product at different photon frequencies were determined. Time-dependent density functional theory was used to analyze the highest five occupied and lowest five unoccupied Br2 orbitals. The transition to the 1? u state was found to be most probable in the visible-near UV absorption range.

Energy Technology Data Exchange (ETDEWEB)

The primary objective is to perform a fundamental study of filtration of nanoparticles, and to obtain filtration knowledge necessary to design particle collection devices/systems for nanoparticle processing and for preventing nanoparticle emissions into the environment. The research covered a wide area relevant to nanoparticle filtration, under these main topics: (1) nanoparticle filtration and molecular dynamics simulation, (2) nanoparticle virtual impactor, (3) particle transport under low pressure, and (4) development of a high-throughput nanoparticle generator. A number of novel tools and numerical models have been developed under the DOE support.

International Nuclear Information System (INIS)

Environmental computer science is a new partial discipline of applied computer science, which makes use of methods and techniques of information processing in environmental protection. Thanks to the inter-disciplinary nature of environmental problems, computer science acts as a mediator between numerous disciplines and institutions in this sector. The handbook reflects the broad spectrum of state-of-the art environmental computer science. The following important subjects are dealt with: Environmental databases and information systems, environmental monitoring, modelling and simulation, visualization of environmental data and knowledge-based systems in the environmental sector. (orig.).

Energy Technology Data Exchange (ETDEWEB)

Analytical and numerical calculations are presented for a reflexing electron beam type of collective ion accelerator. These results are then compared to those obtained through experiment. By constraining one free parameter to experimental conditions, the self-similar solution of the ion energy distribution agrees closely with the experimental distribution. Hence the reflexing beam model appears to be a valid model for explaining the experimental data. Simulation shows in addition to the agreement with the experimental ion distribution that synchronization between accelerated ions and electric field is phase unstable. This instability seems to further restrict the maximum ion energy to several times the electron energy.

Energy Technology Data Exchange (ETDEWEB)

This conference is composed of 21 communications and 21 posters in the domain of building ventilation and indoor air quality; the main themes are: indoor air quality assessment and optimization; performance enhancement and optimization of ventilation systems and equipment; ventilation systems for renovated and rehabilitated buildings; French and European regulations, standardizations and certifications; experimental and numerical simulation studies concerning ventilation systems, air flow, temperature distribution, air quality, thermal comfort and ventilation-induced acoustic levels

International Nuclear Information System (INIS)

The asymmetric fingered structure for polysilicon thin-film transistors (AF-TFTs) is analysed in detail by combining experimental characteristics and two-dimensional numerical simulations. This structure allows an effective reduction of the kink effect and off-current, without introducing any additional series resistance. In addition, a substantial improvement in the device stability is also observed when compared to conventional TFT. The AF-TFT characteristics have been explained by considering a two transistor model.

Energy Technology Data Exchange (ETDEWEB)

Longitudinal space charge (LSC) force can be a main effect driving the microbunching instability in the linac for an x-ray free-electron laser (FEL). In this paper, the LSC-induced beam modulation is studied using an integral equation approach that takes into account the transverse (radial) variation of LSC field for both the coasting beam limit and bunched beam. Changes of beam energy and the transverse beam size can be also incorporated. We discuss the validity of this approach and compare it with other analyses as well as numerical simulations.

CERN Document Server

An algorithm for pose and motion estimation using corresponding features in images and a digital terrain map is proposed. Using a Digital Terrain (or Digital Elevation) Map (DTM/DEM) as a global reference enables recovering the absolute position and orientation of the camera. In order to do this, the DTM is used to formulate a constraint between corresponding features in two consecutive frames. The utilization of data is shown to improve the robustness and accuracy of the inertial navigation algorithm. Extended Kalman filter was used to combine results of inertial navigation algorithm and proposed vision-based navigation algorithm. The feasibility of this algorithms is established through numerical simulations.

Energy Technology Data Exchange (ETDEWEB)

This paper deals with the efficient simulation of the dynamical behaviour of molten carbonate fuel cells (MCFCs). MCFCs allow an efficient and environmentally friendly energy production via electrochemical reactions. Their dynamics can be described by large scale systems of up to currently 22 nonlinear partial differential algebraic equations (PDAE). The paper also serves as a basis for later parameter identification and optimal control purposes. Therefore, the numerical simulations are particularly based on hierarchically embedded systems of PDAE, first of all in one space dimension. The PDAE are of mixed parabolic-hyperbolic type and are completed by nonlinear initial and boundary conditions of mixed type. For a series of embedded models in one space dimension, the vertical method of lines (MOL) is used throughout this paper. For the semi-discretization in space appropriate difference schemes are applied depending on the ...

Energy Technology Data Exchange (ETDEWEB)

This paper is dedicated to the numerical study of natural ventilation in a room through a large external opening using the CFD code Fluent. Bidimensional numerical simulations are performed for wind speeds up to 5.55 m-s{sup -1} (20 km-h{sup -1}). We propose a two-step solution procedure with grid refinement. When combined with appropriate thermal boundary conditions, this technique appears very efficient at limiting local convergence problems. Through an analysis of the flow pattern in the cavity, we qualitatively explain those results and define a critical Archimedes number. Calculated air change rates are compared to values deduced from empirical correlations of the literature. Agreement is in general poor, which can be explained by the differences between the conditions of our simulations and those of the experiments that lead to those correlations. However, from our ...

Energy Technology Data Exchange (ETDEWEB)

This manual describes the use of the Xyce Parallel Electronic Simulator. Xyce has been designed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: (1) Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). Note that this includes support for most popular parallel and serial computers. (2) Improved performance for all numerical kernels (e.g., time integrator, nonlinear and linear solvers) through state-of-the-art algorithms and novel techniques. (3) Device models which are specifically tailored to meet Sandia's needs, including some radiation-aware devices (for Sandia users ...

Energy Technology Data Exchange (ETDEWEB)

Evaluation and interpretation of CT colonography is based on both 2D and 3D techniques. The 2D techniques are popular mainly because the time needed for evaluation is short. The 3D techniques allow better definition than the 2D techniques, especially of polyps close to folds or at the base of a fold. The evaluation strategies generally accepted so far (primarily 2D, with 3D for problem cases, or vice versa) demand knowledge of both 2D and 3D techniques. Newer 3D visualization techniques help make it possible to acquire more complete and faster recording particularly of areas that are not easily accessible to endoscopic examination. These user-friendly developments are thus well suited to improving the detection and the security of detection of polyps. It must be remembered that ...

Energy Technology Data Exchange (ETDEWEB)

Currently, x-ray mammography is the method of choice in breast cancer screening programmes. As the mammography technology moves from 2D imaging modalities to 3D, conventional computational phantoms do not have sufficient detail to support the studies of these advanced imaging systems. Studies of these 3D imaging systems call for a realistic and sophisticated computational model of the breast. DeBRa (Detailed Breast model for Radiological studies) is the most advanced, detailed, 3D computational model of the breast developed recently for breast imaging studies. A DeBRa phantom can be constructed to model a compressed breast, as in film/screen, digital mammography and digital breast tomosynthesis studies, or a non-compressed breast as in positron emission mammography and breast CT studies. Both the cranial-caudal and mediolateral oblique views can be modelled. The anatomical details inside the phantom include the lactiferous duct system, the ...

Science.gov (United States)

Cosmological N-body and hydrodynamic simulations start with a realization of a random density fluctuation field representing a cosmological model at an early epoch. The density field is often replaced by a set of particles whose positions and velocities are set to conform to the desired density field. Each particle represents a cloud of huge number of real particles. Positions and velocities of particles are subsequently integrated by various numerical codes. We have simulated a set of collisionless collapses of Gaussian density peaks by using the PM and P(3) M codes. We find that in cosmological simulations the physics at scales below the mean particle separation(MPS) is dominated by inaccuracies in describing the initial density field, and cannot be studied even by the high force-resolution codes. Since density fluctuations are ill-defined at scales smaller than MPS, it is desirable not to amplify ...

Energy Technology Data Exchange (ETDEWEB)

An experimental and analytical study was performed to improve understanding of the dynamic impact behavior of carbon steel pipes. The test program addressed two types of pipe impact scenarios using both 2- and 4-in. Sch-80 pipes and elbows. Projectile-on-pipe tests simulated the behavior of a stationary target pipe which is impacted at its center by a larger, more rigid whipping pipe. These target pipes, which contained non-flowing water at about 290{degree}C temperature and ca 8.5 megapascals pressure, exhibited a peak deformation of up to 45% reduction in their diameter. For each test condition, the local deformation at the impact zone is a function of the peak impact force and impact velocity. Pipe-on-wall tests simulated the impact of an elbow at the free end of a cantilevered whipping pipe with a rigid barrier. The peak crush of the impact zone of the elbows was a function of peak impact force. Numerical analyses of ...

Science.gov (United States)

High precision, fast computation speed, as well as a good capability of fault tolerant and reconstruction are required more and more for spacecraft attitude determination system. To realize the above requirement, an approach was presented to the synthesis of federated filters using sigma point technique. In this algorithm, the sigma point technique brought the algorithm a high precision, while the federated structure significantly enhanced the filters' capability of multi-rate information fusion, fault tolerance, and system modularity. Within consideration of computation consumption, a simple information-sharing formulation was derived to adapt to the special property of sigma point distribution, and a dynamical information sharing strategy for multi-rate fusion was developed. A numerical simulation example was employed to give the algorithm a test, where the simulated system contained a suit of gyroscopes; a three-axis ...

Science.gov (United States)

Although remarkable progress has been made in developing technologies for the clean and efficient utilization of coal, the biggest challenge in the utilization of coal is still the protection of the environment. Specifically, electric utilities face increasingly stringent restriction on the emissions of NO{sub x} and SO{sub x}, new mercury emission standards, and mounting pressure for the mitigation of CO{sub 2} emissions, an environmental challenge that is greater than any they have previously faced. The Utah Clean Coal Program addressed issues related to innovations for existing power plants including retrofit technologies for carbon capture and sequestration (CCS) or green field plants with CCS. The Program focused on the following areas: simulation, mercury control, oxycoal combustion, gasification, sequestration, chemical looping combustion, materials investigations and student research experiences. The goal of this program was to begin to integrate the ...

Energy Technology Data Exchange (ETDEWEB)

In technical combustors, flow and combustion processes inside the combustion chamber determine the design and also the safety of the plant. A modelling technique was developed at TU Dresden which comprises both numerical simulations and experimental data and thus achieves a high degree of transparency in its visualisation of the processes inside the combustion chambers under investigation. With the experience gained, numerical models of commercial simulation software (FLUENT, CFX) can be validated by measured values, and potential applications of the programs can be investigated. The first stage of experimental modelling consisted in isothermal (cold) modelling in consideration of flow-relevant similarity characteristics. In the reactive, non-isothermal models, the fuel/air ratio and temperature conditions inside the original combustion chamber were modelled as well. In order to obtain detailed data on ...

Energy Technology Data Exchange (ETDEWEB)

A three-dimensional, three-phase, multiple-component numerical simulator was developed to investigate transport and growth of microorganisms in porous media and the impacts of microbial activities on oil recovery. The microbial activities modeled in this study included: (1) growth, retention, chemotaxis, and end product inhibition of growth, (2) the formation of metabolic products, and (3) the consumption of nutrients. Major mechanisms for microbial enhanced oil recovery (MEOR) processes were modeled as follows: (1) improvement in sweep efficiency of a displacement process due to in situ plugging of highly-permeable production zones by cell mass or due to improved mobility control achieved by increasing the viscosity of the displacing fluid with a biopolymer, and (2) solubilization and mobilization of residual oil in porous media due to the reduction of the interfacial tension between oleic and aqueous phases by the production of a ...

Science.gov (United States)

Subcommittee of the Joint Committee on Atomic Energy of the United States, Outer Space Propulsion by Nuclear Energy, hearings, 85th Cong., 2d sess., 22, ...

UK PubMed Central (United Kingdom)

The processes of neuronal outgrowth and guidance have typically been studied in classic 2D cell culture systems that do not recapitulate topographical cues present in the in vivo extracellular matrix...Full Text Available

Energy Technology Data Exchange (ETDEWEB)

No abstract prepared.

KoreaScience (Korea)

Full Text Available

Energy Technology Data Exchange (ETDEWEB)

Blade fatigue life is an important element in determining the economic viability of the Vertical-Axis Wind Turbine (VAWT). A principal source of blade fatigue is thought to be the stochastic (i.e., random) aerodynamic loads created by atmospheric turbulence. This report describes the theoretical background of the VAWT Stochastic Aerodynamic Loads (VAWT-SAL) computer code, whose purpose is to numerically simulate these random loads, given the rotor geometry, operating conditions, and assumed turbulence properties. A Double-Multiple-Stream Tube (DMST) analysis is employed to model the rotor's aerodynamic response. The analysis includes the effects of Reynolds number variations, different airfoil sections and chord lengths along the blade span, and an empirical model for dynamic stall effects. The mean ambient wind is assumed to have a shear profile which is described by either a power law or a logarithmic variation with height above ...

Energy Technology Data Exchange (ETDEWEB)

Simulation and experimental results on the transport of microbes and nutrients in one-dimensional cores are presented, and the development of a three-dimensional, three-phase, multiple-component numerical model to describe the microbial transport and oil recovery in porous media is described. The change of rock's wettability and associated relative permeability values after microbial treatments were accounted for in the model for additional oil recovery. Porosity and permeability reductions due to cell clogging have been considered and the production of gas by microbial metabolism has been incorporated. Governing equations for microbial and nutrient transport are coupled with continuity and flow equations under conditions appropriate for a black oil reservoir. The computer simulator has been used to determine the effects of various transport parameters on microbial transport phenomena. The model can accurately ...

Energy Technology Data Exchange (ETDEWEB)

In the groundwater flow modelling part of this work the effective dilution volume in the well scenario was estimated by means of transient simulations of groundwater flow and transport, which are coupled due to the varying salinity. Both deep, drilled wells and shallow surface wells in the vicinity of the repository were considered. The simulations covered the time period from the present to 1000 years after the present. Conceptually the fractured bedrock consists of planar fracture zones (with a high fracture density and a greater ability to conduct water) and the intact rock (in which the fracture density and the hydraulic conductivity are low). For them the equivalent-continuum model was applied separately. Thus, the fractured bedrock was considered as piecewise homogeneous (except for the depth dependence) and isotropic continuum with representative average characteristics. A generic simulation model for groundwater ...

Energy Technology Data Exchange (ETDEWEB)

A numerical investigation of a humidification dehumidification desalination (HDD) process using solar energy is presented. The HDD system consists mainly of a concentrating solar water heating collector, flat plate solar air heating collector, humidifying tower and dehumidifying exchanger. Two separate circulating loops constitute the HDD system, the first for heating the feed water and the second for heating air. A mathematical model is developed, simulating the HDD system, to study the influence of the different system configurations, weather and operating conditions on the system productivity. The model validity is examined by comparing the theoretical and experimental results of the same authors. It is found that the results of the developed mathematical model are in good agreement with the experimental results and other published works. The results show also that the productivity of the unit is strongly influenced by the air flow rate, ...

Energy Technology Data Exchange (ETDEWEB)

Increase usage of magnesium components in the portable electronics industries has propelled considerable research interest on the development of magnesium forging. There have been some studies done in the area of sheet metal forming and forging of axis-symmetrical parts of magnesium alloys but there are very little work done on the plane strain backward extrusion of magnesium alloy. This paper deals with the analysis of plane strain extrusion of magnesium alloy, AZ31, using numerical modeling at different extrusion conditions of temperatures and friction coefficients. A commercially available finite element package, ANSYS/LSDYNA, was used in the finite element analysis. FE simulated results, material deformation, punch force and stress/strain distributions, were compared with actual experimental data. From the numerical analysis, a decrease in forging temperature from 300 C to 250 C will lead to a moderate increase in ...

Science.gov (United States)

The PLANCK SURVEYOR is a European Space Agency satellite mission to image the very faint anisotropies in the temperature of the Cosmic Microwave Background (CMB) radiation. Maynooth is actively participating in an international collaboration of scientists involved with the optical design of the High Frequency Instrument (HFI). This paper outlines research which has been undertaken in Maynooth concerned with numerical modelling of the optical characteristics of the multi-frequency array of detectors making up the HFI instrument. In the study the commercial software package ZEMAX was used to model the coupling of the focal plane HFI detectors to the PLANCK telescope. This package is particularly useful in the optical design of PLANCK because of the powerful optimisation features of the software. It is thus possible to readily determine the optimum positioning of the detectors in the focal plane of the telescope. Although the package is based on a ...

Energy Technology Data Exchange (ETDEWEB)

Moving out of the scientific community research laboratories, computational fluid dynamics (CFD) software packages are now allowing industrials to analyse and optimize industrial processes involving the use of gases, liquids and even some two-phase fluids. Their attractiveness and their impact stems out from the opportunity they offer to bring insight into an existing unit, or even at the design stage, by displaying the spatial distribution of process relevant variables such as temperature, concentration. The filling of the spacing in between a two-layer window is a simple example. This new opportunity of visualisation is at times an unique way, when the process environment is an opaque one, such as liquid metal flowing into a tundish or when measurements of flows may be a long and tedious work, such as flows within water treatment basins. This environment we are to investigate in order to optimize can also be a harsh one, due to its high temperature level for example. Such are ...

Energy Technology Data Exchange (ETDEWEB)

Excessive erosion of the labyrinth seal of a 100 MW geothermal turbine has been investigated. This study used computational fluid dynamics (CFD) and aims to identify one cause of erosion and a possible solution for substantially reducing it. The predictions were based upon a numerical calculation using a CFD model of the labyrinth seal with a water/steam flow containing hard solid particles and solved with a commercial CFD code: Fluent V5.0. The results confirmed the existence of flow conditions that play a major role in the rotor labyrinth seal erosion. Afterwards, the flow path was simulated with changes of rotor labyrinth seal geometry, which are indeed feasible of being implemented. The results confirmed that it is possible to reduce the erosion process by approximately 80% by incorporating a steam flow deflector in the fourth stage diaphragm, which changes the steam flow direction in the inlet zone to the rotor labyrinth seal channel, ...

Energy Technology Data Exchange (ETDEWEB)

This is a numerical verification of two groups of pipe whip tests sponsored or cosponsored by EPRI. Experimental data of the two pipe whip tests, one by Tennessee Valley Authority (TVA) and by FRAMATOME/CEA, were provided by EPRI. A nonlinear finite element code, ABAQUS-EPGEN, developed under partial sponsorship by EPRI was used for modeling the pipe whip tests. Beam elements together with an equivalent nonlinear spring element or a partial shell mesh were used to model pipes and elbow in the pipe whip tests. Material nonlinearity due to plasticity, strain rate effects, and temperature, as well as geometric nonlinearity due to large rotation and boundary conditions were included in the study. Effects of strain rate and modeling techniques were assessed. Results by current industry approach were also included as a reference solution. This report can be used as a guideline for numerical simulation of pipe whip phenomena. 74 ...

Energy Technology Data Exchange (ETDEWEB)

This paper describes the comprehensive modeling method of fluid and heat flows in the hydrological system. Based upon the concept that the hydrological system is composed of the air-water 2-phase fluids, and the rock phase, various natural processes are modeled including surface/subsurface 2-phase mass transfer, heat exchanges between fluid phases, between fluid and solid phases, and sensible/latent heat exchanges on the ground surface. A field-oriented numerical simulator is developed, in which a set of governing equations is solved for different variables on the surface and in subsurface regions, respectively. Results of an experimental study are presented, in which transient formulation of heat exchange between fluid/solid phases is examined through matching of the observation and calculated performances. (author)

International Nuclear Information System (INIS)

As a three-year joint university-industry effort, development of a generalized boiling transition analysis method has been started in 2002 aiming at enhanced capabilities of subchannel analysis for a wide variety of BWR-type fuel bundle geometry from ordinary BWR to tight lattice fuel bundles. For this purpose, five dominant factors affecting boiling transition phenomena have been identified on which our efforts of experimentation and numerical analyses are focused. In this report, as the first-year achievement, we will describe a master plan of the development and contents for experimental approaches to construct thermal-hydraulic databases. The databases will be utilized for the developments of constitutive equations to describe the basic characteristics of the elementary processes. The planned experiments are divided into two groups. One is air-water experiments at atmospheric pressure, and the other is steam-water experiments up to 1 MPa. The former group of ...

International Nuclear Information System (INIS)

DF (decontamination factor) of two types of air filter was investigated against "1"9"8Au aerosol particles with AMAD (activity median aerodynamic diameter) = 1.2 #mu#m and #sigma#g (geometric standard deviation) = 2.0. One of the filters, SO filter, is HEPA (high efficiency particulate air) filter media with 46 mm#phi#, and the other is HE-40T (cellulose-glass fiber filter for dust sampling) filter media with 40 mm#phi#. DF of the SO filter was found to be over 5 x 10"4 at air flow rate of 1 l/min. And DF of HE-40 T filter was determined to be 94#+-#4 at the same airflow rate. In the multistage HE-40 T filter system, DF decreased at each succeeding stage. It was confirmed by numerical simulation that the decrease of DF was not inconsistent with the decrease of collection performance in each successive filter. It was also confirmed that DF could be predicted by the numerical simulaiton, if particle penetration function of ...

Energy Technology Data Exchange (ETDEWEB)

The convection heat loss from cavity receiver in parabolic dish solar thermal power system can significantly reduce the efficiency and consequently the cost effectiveness of the system. It is important to assess this heat loss and subsequently improve the thermal performance of the receiver. This paper aims to present a comprehensive review and systematic summarization of the state of the art in the research and progress in this area. The efforts include the convection heat loss mechanism, experimental and numerical investigations on the cavity receivers with varied shapes that have been considered up to date, and the Nusselt number correlations developed for convection heat loss prediction as well as the wind effect. One of the most important features of this paper is that it has covered numerous cavity literatures encountered in various other engineering systems, such as those in electronic cooling devices and buildings. The studies related ...

Science.gov (United States)

The main purpose of this paper is to explore a numerical algorithm for determining the contact stress when a circular crowned roller is compressed between two plates. To start with, the deformation curve on a plate surface will be derived by using the contact mechanical model. Then, the contact stress distribution along the roller which occurs on the plate surface is divided into three parts: from the center of contact to the edge, the edge and apart from the contact line. The first part is calculated by the elastic contact theorem for the contact subjected to nominal stress between non-crowned parts of roller and plates, the second part is obtained by the classical Hertzian contact solution for the contact between crowned parts of roller and plates, and the third part is simulated as exponential decay. In order to overcome the defect of the half space theorem, in which a plate with infinite thickness is assumed initially, a weighting method is ...

International Nuclear Information System (INIS)

It is well known that one of the main advantages of the high speed forming (HSF) processes is the improvement in the forming limits of the used materials.Using the Electromagnetic Forming (EMF) technology two materials have been tested with different mechanical and physical properties: the AA5754 aluminium and the AZ31B magnesium alloys.The EMF process principle can be described as follows: A significant amount of electrical energy is stored in a bank of capacitors which are suddenly discharged releasing all the stored energy. This electric discharge runs through a coil which generates an intense transient magnetic field. At the same time transient Eddy currents are induced in the electrically conductive part placed some millimetres far from the coil. Another intense magnetic field is generated due to those Eddy currents but on the opposite direction as the one generated by the coil. A big magnetic repulsion force is created between the part and the coil. This magnetic repulsion ...

UK PubMed Central (United Kingdom)

The halogenated benzimidazoles BDCRB (2-bromo-5,6-dichloro-1-;2;-d-riborfuranosyl benzimidazole riboside) and TCRB (2,5,6-trichloro-1-;2;-d-riborfuranosyl benzimidazole...Full Text Available