CERN Document Server

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

Science.gov (United States)

Apparent evolution of relativistic flows as traced by radio emission results from a combination of several factors related to propagation of relativistic blobs or shocks, velocity, density and pressure stratification of the underlying flow, plasma instability and (possibly also) phase and time travel effect. This combination can create an intricate and chaotic patterns of the observed morphological changes in radio emission, which complicates the analysis and interpretation of kinematic and physical properties of the jet plasma. Recent studies have indicated that slow and quasi-stationary patterns in jets are most likely formed by plasma instabilities while faster, superluminally moving patterns are related to highly relativistic plasma condensations produced by the nuclear flares. Some of the stationary patterns may also be related to recollimation shocks or locations where strong non-thermal continuum is produced in jets. Similarities and differences of the AGN ...

International Nuclear Information System (INIS)

We develop a continuum description for the axial segregation of granular materials in a long rotating drum based on the dynamics of the thin near-surface granular flow coupled to bulk flow. The equations of motion are reduced to the one-dimensional system for two local variables only, the concentration difference and the dynamic angle of repose, or the average slope of the free surface. The parameters of the system are established from comparison with experimental data. The resulting system describes both initial transient traveling wave dynamics and the formation of quasi-stationary bands of segregated materials. A long-term evolution proceeds through slow logarithmic coarsening of the band structure which is analogous to the spinoidal decomposition described by the Cahn-Hilliard equation. copyright 1999 The American Physical Society

Energy Technology Data Exchange (ETDEWEB)

In this thesis a study is made of the of voltage characteristics and reactive power control in the transmission network of the Electric Power Systems (EPS) with static compensators of VARs (SCVs). The necessity of voltage and reactive power control and the application of SCVs as a solution means for three representative operation stages of the EPS, considering balanced operation conditions and fundamental frequency is analyzed: Operation in the presence of incremental voltage changes, quasi-stationary operation or steady state and behavior at the occurrence of big disturbances. The application of SCVs to systems of integral transmission is in general a closed loop control system, reason why the analysis of the interaction of the system EPS-SCV is considered as an essential part of the voltage and reactive power control study. Firstly, the response analysis of the control loop is made for incremental voltage changes, representing the power system by means of its ...

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

Science.gov (United States)

This viewgraph presentation provides an overview of research designed to provide a numerical

British Library Electronic Table of Contents (United Kingdom)

The shock stand-off distance of a spherical model flying with transonic speeds is determined through numerical simulations. The model decelerates due to drag forces caused by the pressure and viscous shear stress at the model surface. Two-dimensional axisymmetric numerical codes with numerical grids fixed to the flying spherical model are used in the simulations. Numerically determined shock stand-off distances are compared with experimental data obtained in a previous study as well as with those obtained in our ballistic-range experiments. The numerical results and the experimental data are found to be in good agreement. In addition, the time-dependent shock stand-off distance of a decelerating model is investigated.

International Nuclear Information System (INIS)

An integrated approach to the modelling of hydrogen assisted failure in 316L steel is presented. The approach includes experimental, fractographic, numerical and theoretical analysis of the phenomenon. The physical adequacy of the mechanical models of hydrogen embrittlement (notch extension model and notch cracking model) is discussed by comparing the virtual damage depth (theoretical) predicted by the models with the embrittled zone (microphysical) measured in the fractographic analysis by scanning electron microscopy. In addition, a numerical modelling of hydrogen diffusion is performed, concluding that bulk diffusion is not important in hydrogen embrittlement of 316L steel, so that hydrogen transport accelerated by the microdamage itself should be taken into account. (orig.)

Energy Technology Data Exchange (ETDEWEB)

Collisional-Radiative numerical models are commonly used to design or interpret experiments in atomic physics of laser-created plasmas, including X-ray laser studies. We describe our new code containing several options: average ion, more or less detailed configurations. It consists of an atomic data base coupled to subroutines evaluating ionic populations and emission and absorption coefficients. Numerical results are given to illustrate the capabilities of the code and to compare different models and types of approximation.

UK PubMed Central (United Kingdom)

The ras oncogenes function by indirectly controlling expression of a subset of yet-undefined genes that are crucial for cell growth and differentiation. In a differential display strategy, numerous...Full Text Available

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

Science.gov (United States)

... these features and yields results which are two orders of magnitude faster than those available using conventional numerical ray tracing techniques ...

Energy Technology Data Exchange (ETDEWEB)

It is shown numerically that on a time scale of order 1 msec a neutron star with pion condensation will undergo separation into a core and an envelope with a sharp interface. The envelope will not, however, be ejected.

International Nuclear Information System (INIS)

The aim of the safety calculation is to quantify through numerical modelling the radiological impact (molar flow, human dose) of a potential repository for radioactive waste on the Meuse/Haute Marne site at Bure. A selection process is underway for safety scenarios and their phenomenological and numerical conceptual models upstream from the safety calculation. This involves defining and quantifying the geometric and dimensional representations of the repository for each scenario plus the physical, mathematical and numerical models that reflect its behaviour and the uncertainties associated with all the parameters required to quantify the impact. A summary will be given of these various aspects. The numerical simulations are then performed on the Alliances platform which integrates the various computer codes required for the physical ...

International Nuclear Information System (INIS)

An attempt is taken to explain anomalies in "9"2Mo(d, n)"9"3Tc, "9"2Mo(d, #alpha#)"9"0Nb, "9"4Mo(d, n)"9"5Tc, "9"8Mo(d, n)"9"9Tc, "9"8Mo(d, p)"9"9Mo, "9"8Mo(d, #alpha#)"9"6Nb, "1"0"0Mo(d, p)"1"0"1Mo and "1"0"0Mo(d, n)"1"0"1Tc reactions with input states having a one-particle nature. Thin films saturated with molybdenum isotopes at the approximately 1 mgxcm"-"2 surface density are used as targets. The targets are irradiated by the extracted cyclotron beam. The deuteron energy is 5-12 MeV. The reaction cross sections are determined by the activation analysis method. Quasi-stationary levels of the nucleus-deuteron system are calculated. Weak anomalies revealing in a smooth (d, #alpha#) reaction cross section on sup(92, 98)Mo nuclei, which do not necessarily correlate with anomalies in the (d, n) and (d, p) channels, are observed. The ground states of the (d, #alpha#) reaction products "9"0Nb and "9"6Nb have (8"+) and (6"+) spins, respec,.ively, that testifies to a ...

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

International Nuclear Information System (INIS)

This paper presents a methodology for modeling slug initiation and growth in horizontal ducts. Transient two-fluid equations are solved numerically using a class of high-resolution shock capturing methods. The advantage of this method is that slug formation and growth in a stratified regime can be calculated directly from the solutions to the flow field differential equations. In addition, by using high-resolution shock capturing methods that do not contain numerical diffusion, the discontinuity generated by slugging in the flow field can be modeled with good accuracy. The two-fluid model is shown to be well-posed mathematically only under certain conditions. Under these circumstances, the two-fluid model is capable of correctly predicting and modeling the flow physics. When ill-posed, an unbounded instability occurs in the flow field ...

Energy Technology Data Exchange (ETDEWEB)

Numerical simulation of dynamic fracture of concrete slabs, impacted by steel projectiles, was carried out in this study. The concrete response was described by a continuum damage model. This continuum damage model was originally developed to study rock fragmentation and was modified in the present study with an emphasis on the post-limit structural response. The model was implemented into a transient dynamic explicit finite element code LS-DYNA2D and the code was then used for the numerical simulations. The specific impact configuration of this study follows the experiment series conducted by Hanchak et al. Comparisons between calculated results and measured data were made. Good agreements were found.

Energy Technology Data Exchange (ETDEWEB)

Stylized numerical models of the tooth are usually employed in qualification procedure related to Electronic Paramagnetic Resonance in long-term accidental contamination dose reconstruction. In this work a voxel model was developed from the microCT image data set of a human non contaminated molar tooth. A stylized model, reproducing the characteristics of the voxel model, was also created in order to investigate the level of accuracy that can be obtained in this kind of study.

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)

Energy Technology Data Exchange (ETDEWEB)

The following paper is focused on the numerical modelling of the transient radon diffusion through radon-proof membranes during the measurement of their radon diffusion coefficient. The major aim of such numerical modelling is to increase the accuracy of radon diffusion coefficients derived from the measured data sets. The developed complex ''transient'' numerical model is able to calculate the radon diffusion coefficient with sufficient accuracy from almost any data set - even from a short-time measurement with a non-linear course of results. This numerical model can also be used for various analyses of transient radon transfer processes (e.g. for the calculation of radon distribution curves within the membrane). The following paper presents governing equations for the simulation ...

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

Energy Technology Data Exchange (ETDEWEB)

Results from a numerical model of axisymmetric solar cavity receivers are compared with experimental data for tests of a novel test bed receiver in the Saudi National Laboratories solar furnace. The computed energy transfer rates and temperatures are compared with the experimental data for different receiver geometries, aperture sizes, and operating conditions. In general, the agreement between the numerical model and the experimental data is better for the small-to-midsized apertures than for the large apertures. The analysis indicates that for the larger apertures, the convective heat losses are overpredicted. It also suggests that these losses could be better characterized. Sensitivity analyses show that both the total solar energy input rate and the convective heat-loss coefficient significantly affect the receiver thermal performance and that the distribution of the input solar flux significantly ...

Energy Technology Data Exchange (ETDEWEB)

The aim of this work is the evaluation of the quality of a chemistry-transport model, not by a classical comparison with observations, but by the estimation of its uncertainties due to the input data, to the model formulation and to the numerical approximations. The study of these 3 sources of uncertainty is carried out with Monte Carlo simulations, with multi-model simulations and with comparisons between numerical schemes, respectively. A high uncertainty is shown for ozone concentrations. To overcome the uncertainty-related limitations, a strategy consists in using the overall forecasting. By combining several models (up to 48) on the basis of past observations, forecasts can be significantly improved. This work has been also the occasion of developing an innovative modeling system, named Polyphemus. (J.S.)

Energy Technology Data Exchange (ETDEWEB)

The erosion phenomena of coastal regions in Italy and in the world is remarkable from the environmental point of view. Are showed activities of mechanical defense and numerical models to solve complex problem.

UK PubMed Central (United Kingdom)

Numerous epidemiological studies have consistently demonstrated that individuals who eat more fruits and vegetables (which are rich in carotenoids) and who have higher serum β-carotene levels...Full Text Available

British Library Electronic Table of Contents (United Kingdom)

A mechanical model with local resonators is proposed as an acoustic metamaterial that exhibits an unusual frequency-dependent effective stiffness. If treated as an equivalent elastic solid, its effective Young's modulus can become unbounded or vanishingly small at two respective frequencies. Moreover, in a certain frequency range, the effective Young's modulus would become negative, resulting in a band gap that coincides with this frequency range. The wave attenuation behavior and mechanism are studied through numerical simulations on the acoustic metamaterial model. The capability of the metamaterial to selectively block or filter unwanted waves is demonstrated by a numerical simulation example.

Science.gov (United States)

A general model for the dissolution of particles in multi-component alloys is proposed and analyzed. The model is based on diffusion equations with cross-terms for the several species, combined with a Stefan condition as the equation of motion of the interface between the particle and diffusive phase. Several numerical schemes for the solution of the Stefan problem are proposed and compared. It turns out that diagonalization is useful for numerical purposes. However, for the case of position-dependent diffusion coefficients one has to use a different scheme. Here, we analyze stability and workload of several time integration methods.

International Nuclear Information System (INIS)

The formation of the ions charge and energy distributions of the ions slowed down or randomly changing their charge in the collisions with the medium particles are studied. The effect of the ions dispersion by the charge on the Bragg curve form is investigated. The proposed diffusion approximation for the heavy ions kinetic equation makes it possible to determine simply the parameters of the ions distribution by charge and energy on the whole way of the ions motion. The relation between the ions charge distribution characteristics and the cross sections of the ionization-recombination processes is indicated. The ions distributions, calculated in the proposed analytical model, are compared with the results of the numerical calculations. Good agreement between the analytical, numerical and experimental results is obtained

British Library Electronic Table of Contents (United Kingdom)

A combined numerical and theoretical study is conducted herein on the penetration of semi-infinite targets by jacketed rods with different r"j"0/r"c"0 ratios where r"j"0 and r"c"0 are the radii of the jacket and the core, respectively. The numerical results show that for smaller r"j"0/r"c"0 ratios the u-v relationship changes only a little compared to that of unitary long rod penetrator of the same core material, hence, the u-v relationship of unitary (homogeneous) long rod penetration is also applicable for jacketed rod penetration. Model for cratering in semi-infinite targets by jacketed rods is then suggested by using the laws of conversation of mass, momentum and energy, together with the u-v relationship of unitary (homogeneous) long rod penetration and an analytical model for predict...

Energy Technology Data Exchange (ETDEWEB)

The numerical modelling of gravity changes resulting from the simultaneous effects of mass relocation and rock deformation is described. The method is supported by FEMMA (Finite Element Method for Multipurpose Applications) software, and has been used for modelling expected gravity changes in a large open pit coal mine in Poland. The results are in good agreement with conventional calculations of gravity changes for a simple geometry of relocated mass. 10 refs., 6 figs., 2 tabs.

International Nuclear Information System (INIS)

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

International Nuclear Information System (INIS)

The problem in perturbation theory in general relativity is finding a natural gauge which simplifies the equations considerably. We solved the first and second order perturbation equations in the ''Lorentz gauge'' using the two-timing method. The equations were deduced from the equations of Einstein in the presence of a matter energy-momentum tensor of a perfect fluid with equation of state p = (#gamma#-1) rho. It is found that the 0--0 component of the first order equations together with the gauge conditions are consistent with the Bianchi type V model, without axial symmetry. It represents a ''ground state'' situation of the Bianchi type IX model (Mixmaster model.) The second order equations, describing the reaction back of the gravitational waves on the slowly varying background, give rise to rotation and will disturb the isotropization mechanism. The numerical solutions show all the features which ...

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

Science.gov (United States)

A three-dimensional numerical simulation of the loss-of-feed water transient at the horizontal steam generator of the Kola nuclear power plant is performed. Presented numerical results show transient change of integral steam generator parameters, such as steam generation rate, water mass inventory, outlet reactor coolant temperature, as well as detailed distribution of shell side thermal-hydraulic parameters: swell and collapsed levels, void fraction distributions, mass flux vectors, etc. Numerical results are compared with measurements at the Kola NPP. The agreement is satisfactory, while differences are close to or below the measurement uncertainties. Obtained numerical results are the first ones that give complete insight into the three-dimensional and transient horizontal steam generator thermal-hydraulics. Also, the presented results serve as benchmark tests for the assessment and further ...

Energy Technology Data Exchange (ETDEWEB)

A rigorous mathematical model is developed to describe chemical absorption of carbon dioxide in a falling-film micro-contactor. Contrary to previous models described in literature, hydrodynamics in both phases is accurately captured by the Navier-Stokes equations. By this means, backmixing effects can be accounted for. Furthermore, the interface position is determined directly using the level-set method. The developed model is validated with experimental data obtained from literature and a good agreement is established. With the validated model, detailed studies are performed to investigate the impact of process parameters.

Energy Technology Data Exchange (ETDEWEB)

This paper proposes a high accuracy numerical method to model oil spill trajectories using a particle-tracking algorithm. The Euler method, used to calculate oil trajectories, can give adequate solutions in most open ocean applications. However, this method may not predict accurate particle trajectories in certain highly non-uniform velocity fields near coastal zones or in river problems. Simple numerical experiments show that the Euler method may also introduce artificial numerical dispersion that could lead to overestimation of spill areas. This article proposes a fourth-order Runge-Kutta method with fourth-order velocity interpolation to calculate oil trajectories that minimise these problems. The algorithm is implemented in the OilTrack model to predict oil trajectories following the 'Nissos Amorgos' oil spill accident that occurred in the Gulf of Venezuela in ...

Energy Technology Data Exchange (ETDEWEB)

We present a 1-D heat transfer, melting, vaporization and resolidification model describing the interaction of a scanning continuous-wave laser with a metal surface wherein the beam power is constant. A physical model based on the Stefan problem is developed with appropriate boundary conditions. The effects of processing parameters on process variables are investigated numerically by varying beam diameters, scan speeds and substrate temperatures for Nickel. Relations are derived for the times to initiate melting, to initiate vaporization, to reach maximum melting depth, for melting-resolidification, and for maximum melting and vaporization depths. Surface temperatures are compared with approximate closed form solutions. (Author)

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)

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)

British Library Electronic Table of Contents (United Kingdom)

The use of numerical modeling in oil spill incidents is a well established technique that has proven to provide cost-effective and reasonable estimates of oil surface drift. Good predictability of such models depends highly on the quality of the input data of the incident and on the model calibration effort. This paper presents the results of simulating oil spillage trajectory in the Arabian (Persian) Gulf. The study employed a 3-D rectilinear hydrodynamic model combined with oil spill model. Typical representative environmental conditions of the Arabian Gulf were first setup into a hydrodynamic circulation model using data from various sources. The performance of the hydrodynamic model was then tested against measurements of tidal fluctuation and sea currents at selected locations. The sp...

Science.gov (United States)

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.

International Nuclear Information System (INIS)

The objective of this study is to establish evaluation and verification guideline for the APR 1400 and to investigate the thermal-hydraulic characteristics for fluidic device is analyzed using FLUENT. The scope and major results of research are flow characteristics for fluidic device. In this study, three-dimensional numerical model for fluidic device is developed adequately for, and results are compared with experimental data performed by VAPER (VAlve Performance Evaluation test Rig) in KAERI with an aim to verify numerical simulation. In addition, the parametric study has also carried out to investigate the effect of major parameters such as velocity and pressure inside FD chamber.

CERN Document Server

We derive the expression for spontaneous emission rate in finite one-dimensional photonic crystal with arbitrary defects using the effective resonator model to describe electromagnetic field distributions in the structure. We obtain explicit formulas for contributions of different types of modes, i.e. radiation, substrate and guided modes. Formal calculations are illustrated with a few numerical examples, which demonstrate that the application of effective resonator model simplifies interpretation of results.

International Nuclear Information System (INIS)

A numerical study of a natural convection in a rectangular cavity with the low-Reynolds-number differential stress and flux model is presented. The primary emphasis of the study is placed on the investigation of the accuracy and numerical stability of the low-Reynolds-number differential stress and flux model for a natural convection problem. The turbulence model considered in the study is that developed by Peeters and Henkes (1992) and further refined by Dol and Hanjalic (2001), and this model is applied to the prediction of a natural convection in a rectangular cavity together with the two-layer model, the shear stress transport model and the time-scale bound #upsilon#"2-f model, all with an algebraic heat flux model. The computed results are compared with the experimental data ...

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

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 higher angle of inclinations (60 {<=} {beta} {<=} 90 ) of the receiver. The present 3-D ...

International Nuclear Information System (INIS)

The development of light-weight vehicle is in great demand for enhancement of fuel efficiency and dynamic performance. The vehicle weight can be reduced effectively by using lightweight materials such as magnesium alloys. However, the use of magnesium alloys in sheet forming processes is still limited because of their low formability at room temperature and the lack of understanding of the forming process of magnesium alloys at elevated temperatures. In this study, uniaxial tensile tests of the magnesium alloy AZ31B-O at various temperatures were performed to evaluate the mechanical properties of this alloy relevant for forming of magnesium sheets. To construct a FLD (forming limit diagram), a forming limit test were conducted at temperature of 100 and 200 deg. C. For the evaluation of the effects of the punch temperature on the formability of a rectangular cup drawing with AZ31B-O, numerical modelling was conducted. The experiment results ...

CERN Document Server

We obtain an elegant and useful description of the dynamics of Szekeres dust models (in their full generality) by means of "quasi--local" scalar variables constructed by suitable integral distributions that can be interpreted as weighed proper volume averages of the local covariant scalars. In terms of these variables, the field equations and basic physical and geometric quantities are formally identical to their corresponding expressions in the spherically symmetric LTB dust models. Since we can map every Szekeres model to a unique LTB model, rigorous results valid for the latter models can be readily generalized to a non--spherical Szekeres geometry. The new variables lead naturally to an initial value formulation in which all scalars are expressed as scaling laws in terms of their values at an arbitrary initial space slice. These variables also yield a significant simplification ...

Energy Technology Data Exchange (ETDEWEB)

Thermal striping is identified as one of the causes of thermal fatigue failure in nuclear power plants. Numerical studies of thermal striping require three-dimensional, unsteady turbulent modeling that resolves both large and small-scale turbulent motions. Benchmark studies were carried out using the LES turbulence model solved by the commercial CFD code FLUENT. Two types of mixing tee configurations were modeled to evaluate the performance of the CFD code. The simulation results presented in normalized average temperature and normalized fluctuating temperatures are in good agreement with measurements.

CERN Document Server

Heat exchange during laminar flow is studied at the fracture scale on the basis of the Stokes equation. We used a synthetic aperture model (a self-affine model) that has been shown to be a realistic geometrical description of the fracture morphology. We developed a numerical modelling using a finite difference scheme of the hydrodynamic flow and its coupling with an advection/conduction description of the fluid heat. As a first step, temperature within the surrounding rock is supposed to be constant. Influence of the fracture roughness on the heat flux through the wall, is estimated and a thermalization length is shown to emerge. Implications for the Soultz-sous-For\\^{e}ts geothermal project are discussed.

Energy Technology Data Exchange (ETDEWEB)

This executive summary presents the motivation, structure, objectives, methodologies and results of the first stage of the international DECOVALEX project - DECOVALEX I (1992-1995). The acronym stands for Development of Coupled Models and their Validation against Experiment in Nuclear Waste Isolation, and the project is an international effort to develop mathematical models, numerical methods and computer codes for coupled thermo-hydro-mechanical processes in fractured rocks and buffer materials for geological isolation of spent nuclear fuel and other radioactive wastes, and validate them against laboratory and field experiments. 24 refs.

British Library Electronic Table of Contents (United Kingdom)

This paper studies the pricing of variance swap derivatives with stochastic volatility by the control variate method. A closed form solution is derived for the approximate model with deterministic volatility, which plays the key role in the paper, and an efficient control variate technique is therefore proposed when the volatility obeys the log-normal process. By the analysis of moments for the underlying processes, the optimal volatility function in the approximate model is constructed. The numerical results show the high efficiency of our method; the results coincide with the theoretical results. The idea in the paper is also applicable for the valuation of other types of variance swap, options with stochastic volatility and other financial derivatives with multi-factor models.

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)

Analytical and numerical methods have been used in this investigation to model the behavior of geothermal systems under exploitation. The work is divided into three parts: (1) development of a numerical code, (2) theoretical studies of geothermal systems, and (3) field applications. A new single-phase three-dimensional simulator, capable of solving heat and mass flow problems in a saturated, heterogeneous porous or fractured medium has been developed. The simulator uses the integrated finite difference method for formulating the governing equations and an efficient sparse solver for the solution of the linearized equations. In the theoretical studies, various reservoir engineering problems have been examined. These include (a) well-test analysis, (b) exploitation strategies, (c) injection into fractured rocks, and (d) fault-charged geothermal reservoirs.

Energy Technology Data Exchange (ETDEWEB)

In this paper, the moisture buffering capacity of spruce plywood is measured by recording the change in mass of a test specimen when the air relative humidity (RH) is changed between 33% RH and 75% RH. The aim is to represent diurnal cycles in indoor humidity with 33% RH maintained for 16 h and 75% RH maintained for 8 h. Measurements are taken using two different apparatuses, which provide different convective transfer coefficients between the air and the plywood, and the results are compared to a numerical model for validation. The validated numerical model is then used to investigate the effect of initial conditions, boundary conditions and thickness on the moisture buffering capacity of plywood. The results show that the buffering capacity of plywood depends on the initial conditions and thickness of the plywood as well as the surface film coefficient and humidity cycle. (author)

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)

A numerical model of multiphase air-water flow and contaminant transport in the unsaturated zone is presented. The multiphase flow equations are solved using the two-pressure, mixed form of the equations with a modified Picard linearization of the equations and a finite element spatial approximation. A volatile contaminant is assumed to be transported in either phase, or in both phases simultaneously. The contaminant partitions between phases with an equilibrium distribution given by Henry`s Law or via kinetic mass transfer. The transport equations are solved using a Galerkin finite element method with reduced integration to lump the resultant matrices. The numerical model is applied to published experimental studies to examine the behavior of the air phase and associated contaminant movement under water infiltration. The model is also used to evaluate a hypothetical design for a ...

International Nuclear Information System (INIS)

The strong coupling of subsystem elements, such as the accelerator, wiggler, and optics, greatly complicates the understanding and design of a free electron laser (FEL), even at the conceptual level. To address the strong coupling character of the FEL the concept of an Integrated Numerical Experiment (INEX) was proposed. Unique features of the INEX approach are consistency and numerical equivalence of experimental diagnostics. The equivalent numerical diagnostics mitigates the major problem of misinterpretation that often occurs when theoretical and experimental data are compared. The INEX approach has been applied to a large number of accelerator and FEL experiments. Overall, the agreement between INEX and the experiments is very good. Despite the success of INEX, the approach is difficult to apply to trade-off and initial design studies because of the significant manpower and computational requirements. On the other hand, ...

British Library Electronic Table of Contents (United Kingdom)

For optimal performances, proton exchange membrane fuel cells require fine water and thermal management. Accurate modelling of the physical phenomena occurring in the fuel cell is a key issue to improve fuel cell technology. Here, an analytic steady state diphasic 2D model of heat and mass transfer is presented. Through this model, the aim of this work is to study the influence of local events on the global performances of a fuel cell. A part of the complete model is a microscopic representation of the coupling between water transport and charge transfers in the electrodes. The thickness of the liquid layer around the reactive agglomerates is deduced from the saturation. The evolution of the quantity of water within the catalyst layer is monitored and its influence on the global performanc...

British Library Electronic Table of Contents (United Kingdom)

Regional Ocean Modeling System (ROMS v 3.0), a three-dimensional numerical ocean model, was previously enhanced for shallow water applications by including wave-induced radiation stress forcing provided through coupling to wave propagation models (SWAN, REF/DIF). This enhancement made it suitable for surf zone applications as demonstrated using examples of obliquely incident waves on a planar beach and rip current formation in longshore bar trough morphology (Haas and Warner, 2009). In this contribution, we present an update to the coupled model which implements a wave roller model and also a modified method of the radiation stress term based on Mellor (2008, 2011a,b,in press) that includes a vertical distribution which better simulates non-conservative (i.e., wave breaking) processes and ...

Energy Technology Data Exchange (ETDEWEB)

We suggest a procedure for estimating uncertainties in neutron cross sections calculated with a nuclear model descriptive of a specific mass region. It applies standard error propagation techniques, using a model-parameter covariance matrix. Generally, available codes do not generate covariance information in conjunction with their fitting algorithms. Therefore, we resort to estimating a relative covariance matrix a posteriori from a statistical examination of the scatter of elemental parameter values about the regional representation. We numerically demonstrate our method by considering an optical-statistical model analysis of a body of total and elastic scattering data for the light fission-fragment mass region. In this example, strong uncertainty correlations emerge and they conspire to reduce estimated errors to some 50% of those obtained from a naive uncorrelated summation in quadrature. 37 ...

Science.gov (United States)

The Predictive Capability Maturity Model (PCMM) is a new model that can be used to assess the level of maturity of computational modeling and simulation (M&S) efforts. The development of the model is based on both the authors experience and their analysis of similar investigations in the past. The perspective taken in this report is one of judging the usefulness of a predictive capability that relies on the numerical solution to partial differential equations to better inform and improve decision making. The review of past investigations, such as the Software Engineering Institute's Capability Maturity Model Integration and the National Aeronautics and Space Administration and Department of Defense Technology Readiness Levels, indicates that a more restricted, more interpretable method is needed to assess the maturity of an M&S effort. The PCMM ...

International Nuclear Information System (INIS)

A simple model is used to explain the decrease in superconducting critical temperature with damage observed for irradiated A-15 compounds. A truncated t-matrix approximation is used to describe the disorder along the one-dimensional transition metal chains. Three dimensionality is introduced by the inclusion of interaction between transition metal atoms on different chains. Numerical fits to experiment are discussed in the conclusion. (author).

Energy Technology Data Exchange (ETDEWEB)

The topics covered in this session include: slimhole testing and data acquisition, theoretical and numerical models for slimholes, and an overview of the analysis of slimhole data acquired by the Japanese. The fundamental issues discussed are concerned with assessing the efficacy of slimhole testing for the evaluation of geothermal reservoirs. the term reservoir evaluation is here taken to mean the assessment of the potential of the geothermal reservoir for the profitable production of electrical power. As an introduction to the subsequent presentations and discussions, a brief summary of the more important aspects of the use of slimholes in reservoir evaluation is given.

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.

International Nuclear Information System (INIS)

A numerical method for studying the flow field of an under-expanded axis-symmetric jet created by a converging-diverging nozzle and impinging on a flat plate is presented. The calculation domain of interest is a region that contains all the features of the jet that leaves the nozzle and impinges onto the plate. Axial symmetry of the domain allows a reduced 2-dimensional model to be used. FLUENT software is utilized to solve the continuity, momentum and energy equations using a coupled implicit scheme. The ideal-gas law is used to determine the gas density along with a k-#epsilon# turbulence model with a special modification to account for compressibility effects. In addition, the temperature dependency of viscosity has been taken into consideration. A number of different modeling techniques are investigated including different approximations to account for the flow inside the nozzle. Performance of ...

International Nuclear Information System (INIS)

The effect of two-body nature of the nuclear shell model potential on the recent numerical calculations of the nuclear level density has been examined. For the two most widely used single particle energy level schemes based on harmonic oscillator and Woods-Saxon potential, this effect is shown to significantly modify the excitation energy dependence of the level densities. (author).

Energy Technology Data Exchange (ETDEWEB)

The concept of a multi-cavity volumetric solar receiver is very attractive for the profitability of certain of its characteristics such as high efficiency and economy. The absorber is based on a pack of small ceramic cavities which intercept and absorb the inherent high solar flux reflected from an array of mirrors. Atmospheric air acts as a coolant medium when it is drawn through the pack. A model for an overall heat transfer performance of the receiver is given and numerically solved.

UK PubMed Central (United Kingdom)

Cordycepin has been shown to interfere with a myriad of molecular processes from RNA elongation to kinase activity, and prevents numerous inflammatory processes in animal models. Here we show in a mouse...Full Text Available

Energy Technology Data Exchange (ETDEWEB)

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)

The apex dimensions and length are calculated as a function of current for the ion-emitting jet in liquid-metal ion sources (LMIS). The results agree well with observations. Since the final expressions are analytical, they give more insight into the fundamental mechanisms involved than do numerical calculations. Some implications of the model are discussed concerning focused ion beam (FIB) systems employing LMIS. (author).

Energy Technology Data Exchange (ETDEWEB)

The main objective of this dissertation is to build a numerical 3D geological model of a turbidite reservoir using well and seismic data. For the, the first step was to elaborate a depositional model using a high-resolution genetic stratigraphy approach. Simultaneously, the stratigraphic inversion over the seismic data resulted in an optimal acoustic impedance model. Both models, were used to build, by geostatistical methods, a 3-D geological model, which represents the spatial distribution of the reservoir facies. The studied case is the late Cretaceous Namorado turbidite system in the Campos basin, which is the reservoir of the Namorado oil field, located 80 km from the Brazilian coast, under a bathymetry between 110 and 250 m. The reservoir is a turbidite succession with 90 to 180 m thickness. It is composed by sand-supported conglomerates, matrix-supported ...

Energy Technology Data Exchange (ETDEWEB)

Discussed in this report is a wavefield simulation in the 3-dimensional seismic survey. With the level of the object of exploration growing deeper and the object more complicated in structure, the survey method is now turning 3-dimensional. There are several modelling methods for numerical calculation of 3-dimensional wavefields, such as the difference method, pseudospectral method, and the like, all of which demand an exorbitantly large memory and long calculation time, and are costly. Such methods have of late become feasible, however, thanks to the advent of the parallel computer. As compared with the difference method, the pseudospectral method requires a smaller computer memory and shorter computation time, and is more flexible in accepting models. It outputs the result in fullwave just like the difference method, and does not cause wavefield numerical variance. As the computation platform, the ...

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

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)

CERN Document Server

The study of randomness in low-dimensional quantum antiferromagnets is at the forefront of research in the field of strongly correlated electron systems, yet there have been relatively few experimental model systems. Complementary neutron scattering and numerical experiments demonstrate that the spin-diluted Heisenberg antiferromagnet La2Cu(1-z)(Zn,Mg)zO4 is an excellent model material for square-lattice site percolation in the extreme quantum limit of spin one-half. Measurements of the ordered moment and spin correlations provide important quantitative information for tests of theories for this complex quantum-impurity problem.

Energy Technology Data Exchange (ETDEWEB)

A flow field and particle trajectory in a HDD (Hard Disk Drive) between two rotating disks in axisymmetric enclosures is investigated using CFD code FLUENT/UNS. The RNG k-{epsilon} model is used as a turbulent model. In this study, the flow field between two disks are symmetric, and the flow field near the enclosure is very complex. Cross stream vectors are shown both for blowing and no blowing from the hub. The larger a particle, the more fast the particle deposits at the walls. In the case of blowing from the hub, the more fast the particle deposits at the walls. (author). 9 refs., 12 figs., 1 tab.

CERN Document Server

We briefly report on zero-temperature photoassociation of a Bose-Einstein condensate, focusing on the many-body rate limit for atom-molecule conversion. An upgraded model that explicitly includes spontaneous radiative decay leads to an unanticipated shift in the position of the photoassociation resonance, which affects whether the rate (constant) maximizes or saturates, as well as the limiting value itself. A simple analytical model agrees with numerical experiments, but only for high density. Finally, an explicit comparison with the two-body unitary limit, set by the size of the condensate, finds that the many-body rate limit is generally more strict.

Energy Technology Data Exchange (ETDEWEB)

The effect of the electric field on an exciton confined in a pair of vertically coupled quantum dots is studied. We use a single-band approximation and a parabolic model potential. As a result of these idealizations, we obtain a numerically solvable model, which is used to describe the influence of the electron-hole interaction on the Stark effect for the lowest-energy photoluminescence lines. We show that for intermediate tunnel coupling between the dots this interaction leads to an anomalous Stark effect with an essential deviation of the recombination energy from the usual quadratic dependence on the electric field.

Science.gov (United States)

A computer model of an electrothermal accelerator has been developed which applies to the plasma generator as well as the gun barrel, and it for the first time allows study of the interdependent physical processes in both components simultaneously. The code comprises a 2D transient description of the plasma arc, a 2D, two-phase interior ballistics module, and a preprocessor delivering the required plasma physical data. Some results on the flow field, pressure and temperature in a typical accelerator are presented, and probable limitations on the performance of such a system are discussed. The code is to be applied to different concepts of electrothermal chemical guns.

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

Energy Technology Data Exchange (ETDEWEB)

Reactive chemical transport models developed over the past decade have generally relied on the assumption that local thermodynamic equilibrium is achieved at all times between aqueous species in a given system. Consequently, homogeneous aqueous systems characterized by a number of kinetically slow reactions, particularly problems involving organic species, cannot be satisfactorily modeled. In this study, we present a prototype computer model, KINETRAN, which is designed to handle kinetically-controlled homogeneous reactions in the aqueous phase, along with the transport of the various species involved, through geologic media. 31 refs., 53 figs., 10 tabs.

British Library Electronic Table of Contents (United Kingdom)

The rolls and stand deflections induced by rolling force lead to a significant decrease in the radius prediction accuracy of the ring product in the in-plane roll-bending of strip. To precisely predict the deflections and control the radius, a new analytical model is developed. Numerical implementation is presented to solve the theoretical rolling force, the deflections as well as the inherent force-deflection relationship considering the interactions of the strip workpiece, the rolls and the stand. A series of profile measurements of the formed ring parts are made to assess the accuracy of the predicted deflections. Different control approaches are used in the experiments to validate the radius control model. The experimental results show that the model is reliable to control the radii wi...

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

British Library Electronic Table of Contents (United Kingdom)

We used a 13-year time series of abundance estimates of breeding Northern Saw-whet Owls (Aegolius acadicus), and of small mammals from central Ontario, Canada, to assess the numerical response of the owls to small-mammal prey species. We found that the finite rate of increase of breeding owls was directly related to estimates of red-backed vole (Myodes gapperi) abundance. Thus, it appeared that the owls were nomadic, and made decisions about where to breed based on vole supply. The owls showed a much weaker response to deer mouse (Peromyscus maniculatus) abundance. Across all years, 55% of variation in owl rate of increase could be uniquely attributed to vole abundance, whereas only 3% could be attributed to mouse abundance. Consistent with the model of nomadism, there was only a weak rela...

British Library Electronic Table of Contents (United Kingdom)

Purpose - The purpose of this paper is to investigate the accuracy of different force calculation methods and their impact on mechanical deformations. For this purpose, a micrometer scaled actuator is considered, which consists of a micro-coil and of a permanent magnet (PM) embedded in a deformable elastomeric layer. Design/methodology/approach - For the magnetic field evaluation a hybrid numerical approach (finite element method/boundary element method (FEM/BEM) coupling and a FEM/BEM/Biot-Savart approach) is used, whereas FEM is implemented for the mechanical deformation analysis. Furthermore, for the magneto-mechanical coupling several force calculation methods, namely the Maxwell stress tensor, the virtual work approach and the equivalent magnetic sources methods, are considered and co...

Energy Technology Data Exchange (ETDEWEB)

RODbaffle heat exchanger is a kind of shell-and-tube heat exchanger with eminent performance. Because of the characteristics of tube bundle support structure, fluid flow on shellside is longitudinal and periodical, and fluid flow is of symmetry in landscape orientation. According to the fluid flow characteristics on shellside and ignoring the impact of inlet, outlet and shell wall, a periodic flow unit duct was taken as the simplified model of the shellside to perform numerical simulation by using CFD code, FLUENT. It was found that the both errors in magnitude of the main fluid velocities and heat transfer coefficients between results of simulation and that of experiment or correlations are less than 10%, and the errors of pressure drop between simulation and correlation are no more than 20%. The detailed characteristics and relation of fluid flow and heat transfer on shellside of the RODbaffle heat exchanger were analyzed using the simulation ...

International Nuclear Information System (INIS)

The interaction between a lower hybrid wave and a fusion alpha particle displaces the alpha particle simultaneously in space and energy. This results in coupled diffusion. Diffusion of alphas down the density gradient could lead to their transferring energy to the wave. This could, in turn, put energy into current drive. Here we calculate numerical solutions for the alpha energy transfer and study a range of conditions that are favourable for wave amplification from alpha energy. We find that it is possible for fusion alpha particles to transfer a large fraction of their energy to the lower hybrid wave. The numerical calculation shows that the net energy transfer is not sensitive to the value of the diffusion coefficient over a wide range of practical values. An extension of this idea, the use of a lossy boundary to enhance the energy transfer, is investigated. This technique is shown to offer a large potential benefit. (Author).

Energy Technology Data Exchange (ETDEWEB)

The interaction between a lower hybrid wave and a fusion alpha particle displaces the alpha particle simultaneously in space and energy. This results in coupled diffusion. Diffusion of alphas down the density gradient could lead to their transferring energy to the wave. This could, in turn, put energy into current drive. An initial analytic study was done by Fisch and Rax. Here the authors calculate numerical solutions for the alpha energy transfer and study a range of conditions that are favorable for wave amplification from alpha energy. They find that it is possible for fusion alpha particles to transfer a large fraction of their energy to the lower hybrid wave. The numerical calculation shows that the net energy transfer is not sensitive to the value of the diffusion coefficient over a wide range of practical values. An extension of this idea, the use of a lossy boundary to enhance the energy transfer, is investigated. This technique is ...

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.

British Library Electronic Table of Contents (United Kingdom)

In energy industries, it is always of difficulty to produce high heat value(HHV) gas continuously using pulverized coal. In this paper, a new type furnace for partitioned alternative gasification using pulverized coal is developed, in which the oxidation and reduction reaction occur alternatively with the dropping of pulverized coal and finally HHV gas could be continuously obtained at the reduction zone exit and low heat value(LHV) gas at the oxidation zone exit. Furthermore, the gasification characteristics and their factors in furnace are numerically simulated under two dimensional model with a self-coded program, based on heat, mass and energy transfer as well as reaction principles. It is found that the producing rate of HHV gas is 1.10Nm3/kg with heat value of 11.72MJ/Nm3, however, t...

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 interface developed with ...

Energy Technology Data Exchange (ETDEWEB)

The relationship between the gas to melt ratio (GMR) and the surface temperature of an evolving billet surface in spray forming is investigated numerically. The basis for the analysis is an integrated approach for modelling the entire spray forming process. This model includes the droplet atomization taking thermal coupling into consideration and the deposition of material at the surface of the billet taking geometrical aspects such as shading into account. The coupling between these two models is accomplished by ensuring that the total droplet size distribution of the spray is the summation of ''local'' droplet size distributions along the r-axis of the spray cone. The criterion for a successful process has been a predefined process window characterised by a desired fraction solid range at a certain distance from the atomizer. Inside this process window, the gas and ...

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

Energy Technology Data Exchange (ETDEWEB)

We present a description of the ground state and low-lying excited states of two holes in the 4x4 cluster t-J model in terms of a simple model for the motion of a single bipolaron. The existence of short-range antiferromagnetic correlations has been assumed. According to the suggested scenario, the formation of the bipolaron is mediated by the reduction of the magnetic energy in the case of two holes occupying nearest neighbor sites. The relevant part of the Hilbert space consists of wave functions corresponding to holes oscillating around pairs of nearest neighbor sites and trapped in a potential well due to strings of spin defects. Virtual processes which connect these states involve both the kinetic term and the transverse part of the Heisenberg Hamiltonian. Many properties of energy level schemes obtained by numerical diagonalizations such as the sequence of the lowest states for each irreducible representation of the k ...

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

CERN Document Server

We model the transport of cosmic ray nuclei in the Galaxy by means of a new numerical code. Differently from previous numerical models we account for a generic spatial distribution of the diffusion coefficient. We found that in the case of radially uniform diffusion, the main secondary/primary ratios (B/C, N/O and sub-Fe/Fe) and the modulated antiproton spectrum match consistently the available observations. Convection and re-acceleration do not seem to be required in the energy range we consider: $1 \\le E \\le 10^3$ GeV/nucleon. We generalize these results accounting for radial dependence of the diffusion coefficient, which is assumed to trace that of supernova remnants. While this does not affect the prediction of secondary/primary ratios, the simulated longitude profile of the diffuse $\\gamma$-ray emission is significantly different from the uniform case and may agree with EGRET measurements ...

CERN Document Server

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

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

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

Science.gov (United States)

A butterfly valve is a type of flow control device, typically used to regulate a fluid flowing. Currently, FEA is often used to predict the safety in the design of valve disc. Also, the study about the affection of butterfly valve's disc to the valve flow characteristics by using CFD has been done by many researchers. Along with the development of computer technique, design and analysis of computer experiments has becoming more and more important in engineering design and optimization. Hereinto Kriging model is one popular analysis approach for the purpose of creating a cheap ``meta-model'' as a surrogate to a computationally expensive simulation model. In this paper, the numerical analysis considered the strength, pressure loss coefficient and weight of valve disc simultaneously is investigated to improve the shape of a traditional butterfly valve disc. Firstly, an initial model of ...

Energy Technology Data Exchange (ETDEWEB)

The Water Frame work Directive (December 2000) proposes integrated water management in regard to both quality and quality. A mathematical model has therefore been produced at the CEIT and EPTISA research centre, in collaboration with the Ebro Hydrographic. Confederation (CHE), to predict water quality and thus begin integrated management of it in the river basins. A description is given of the characteristics of the simulation tool. The hydraulic model is based on the numerical (weighted, four-point implicit) resolution of the complete Saint Venant equations. The quality model is based on the IWA River Water Quality Model Number 1 which has clear advantages (consistency, mass balance and easy integration with the biological models of the waste water treatment plants and the collector networks compared to traditional models. The experimental ...

Energy Technology Data Exchange (ETDEWEB)

In the last two decades the performance of numerical wind-wave models has improved considerably. Several models have been routinely producing good quality wave estimates globally since the mid-1980s. The verifications of wind-wave models have mainly focused on the evaluation of the error of the significant wave height H{sub s} estimates. However, for wave energy purposes, the main parameters to be assessed are the wave power P{sub w} and the mean (energy) period T{sub e}. Since P{sub w} is proportional to H{sub s}{sup 2}T{sub e}, its expected error is much larger than for the single-wave parameters. This paper summarizes the intercomparison of two wind-wave models against buoy data in the North Atlantic and the Mediterranean Sea to select the most suitable one for the construction of an Atlas of the wave energy resource in European waters. A full verification in the two basins of ...

Energy Technology Data Exchange (ETDEWEB)

Results are presented from a numerical model of the steady-state energy transfer in molten-salt-in-tube solar cavity receivers that includes convective energy transfer at a local (spatially resolved) level. Molten salt energy absorption and gray radiative transfer between all cavity surfaces are also included. This model is applied to the Molten Salt Subsystem Component Test Experiment (MSS/CTE) cavity receiver. Results for this receiver indicate the global (entire cavity) receiver thermal efficiency is invariant within a few percent to most parameters investigated, although front surface temperatures of the nonabsorbing walls vary considerably, and are particularly sensitive to the type of convective submodel used. Absorption efficiencies indicate the effects of the cavity enclosure environment. For all conditions investigated, tube inner wall temperatures remain under 855 K, ensuring that the salt remains chemically ...

Energy Technology Data Exchange (ETDEWEB)

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

Energy Technology Data Exchange (ETDEWEB)

In recent years sufficient new spectra of actinides in their numerous valence states have been measured to encourage a broader scale analysis effort than was attempted in the past. Theoretical modelling in terms of effective operators has also undergone development. Well established electronic structure parameters for the trivalent actinides are being used as a basis for estimating parameters in other valence states and relationships to atomic spectra are being extended. Recent contributions to our understanding of the spectra of 4+ actinides have been particularly revealing and supportive of a developing general effort to progress beyond a preoccupation with modelling structure to consideration of the much broader area of structure-bonding relationships. We summarize here both the developments in modelling electronic structure and the interpretation of apparent trends in bonding. 60 refs., 9 figs., 1 ...

CERN Document Server

We discuss the foundations of factor or regression models in the light of the self-consistency condition that the market portfolio (and more generally the risk factors) is (are) constituted of the assets whose returns it is (they are) supposed to explain. As already reported in several articles, self-consistency implies correlations between the return disturbances. As a consequence, the alpha's and beta's of the factor model are unobservable. Self-consistency leads to renormalized beta's with zero effective alpha's, which are observable with standard OLS regressions. Analytical derivations and numerical simulations show that, for arbitrary choices of the proxy which are different from the true market portfolio, a modified linear regression holds with a non-zero value $\\alpha_i$ at the origin between an asset $i$'s return and the proxy's return. Self-consistency also introduces ``orthogonality'' and ``normality'' conditions ...

Science.gov (United States)

We present a numerical model of vapor-dominated reservoirs which is based on the well-known conceptual model of White, Muffler, and Truesdell. Computer simulations show that upon heat recharge at the base, a single phase liquid-dominated geothermal reservoir in fractured rock with low matrix permeability will evolve into a two-phase reservoir with B.P.D. (boiling point-for-depth) pressure and temperature profiles. A rather limited discharge event through cracks in the caprock, involving loss of only a few percent of fluids in place, is sufficient to set the system off to evolve a vapor-dominated state. The attributes of this state are discussed, and some features requiring further clarification are identified. 26 refs., 5 figs.

Science.gov (United States)

A full-scale model of a rigid sectorial ring and a set of bellows which was covered with a temperature control layer and a coolant supply equipment was constructed to verify the adequacy of the heating and cooling system designed for JT-60. To cool and heat it effectively and to decrease the temperature differences among the various parts in the vacuum vessel, heater units and cooling pipes were located on the surface of the ring. The temperature control layer is to heat the vessel to 500/degree/sub //C within 70 hours and to maintain it at the temperature for 48 hours for simulating the state of bakeout. Subsequently the vessel is cooled down within about the same time as in heating. Prior to the series of tests, numerical analyses were performed to predict the cooling and heating efficiencies on the model and to examine the method of the temperature regulation.

British Library Electronic Table of Contents (United Kingdom)

A model is developed to represent mechanical strain, stress-enhanced dissolution, and shear dilation as innately hysteretic and interlinked processes in rough contacting fractures. The model is incorporated into a numerical simulator designed to examine permeability change and thermal exchange in chemically active and deformable fractured reservoirs. A candidate engineered geothermal reservoir system (EGS) is targeted. The mechanistic model is able to distinguish differences between the evolution of fluid transmission characteristics of (1) small scale, closely spaced fractures, and (2) large-scale, more widely spaced fractures. Alternate realizations of fracture frequency and scale, exhibiting identical initial bulk permeability, lead to significantly different conclusions regarding perme...

CERN Document Server

The next decade promises an observational revolution which will change cosmology forever. The precise measurement of the angular anisotropy of the cosmic microwave background should specify to a few percent all of the parameters of the cosmological model which effect astrophysics. The growth of structure will then be determined (but not yet observed) until gravitational collapse becomes highly non-linear and stars, galaxies and active galactic nuclei (AGN) form. These processes are hard to model with basic physics because they are complex and allow a rich variety of expression. Instead observations will determine when the first stars and quasars formed, and how and when galaxies assembled. If we can reconcile the numerous contradictions which characterize the subject today, cosmology will become a mature subject, founded on the agreement between detailed, inclusive and realistic models, which make ...

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

Energy Technology Data Exchange (ETDEWEB)

A new stem taper function was established for individual trees of two poplar hybrid clones grown on a short-rotation coppice. The model could be easily fitted and required three parameters to be estimated. It can be used to estimate both diameter at a given height and height for a given top diameter. Two of the three parameters controlled the conical and the neiloid parts of the stem. Significant differences in these parameters were observed between the two clones even if no differences were observed for diameter at breast height or total height of the stem. The model could not be integrated to calculate volumes (total volume, merchantable volume), which were estimated by numerical integration. However, use of this new model allows the optimal length of billets to be determined and thus maximizes the merchantable biomass of poplar in short-rotation coppice by minimizing the biomass of residues.

Energy Technology Data Exchange (ETDEWEB)

Nuclear fuel cycle activities of the former Soviet Union (FSU) have resulted in massive contamination of the environment in western Siberia. We are developing three-dimensional numerical models of the hydrogeology and potential contaminant migration in the West Siberian Basin. Our long-term goal at Pacific Northwest Laboratory is to help determine future environmental and human impacts given the releases that have occurred to date and the current waste management practices. In FY 1993, our objectives were to (1) refine and implement the hydrogeologic conceptual models of the regional hydrogeology of western Siberia developed in FY 1992 and develop the detailed, spatially registered digital geologic and hydrologic databases to test them, (2) calibrate the computer implementation of the conceptual models developed in FY 1992, and (3) develop general geologic and hydrologic information and preliminary ...

International Nuclear Information System (INIS)

A new three-dimensional treatment planning system (TPS) based on convolution/superposition algorithms (TMS-Radix from HELAX AB, Uppsala, Sweden) was recently installed at the University Hospital in Lund. The purpose of the present study was to design a quality assurance and acceptance testing programme to meet the specific characteristics of this convolution model. The model is based on parametrization of a non-measurable quantity - the polyenergetic pencil beam. However, the verification of the treatment planning model is still dependent on numerous comparisons of measured depth-doses and dose profiles. The test programme was divided in two basic parts: (i) model implementation and beam data consistency and (ii) model performance and limitations in special situations. The first part was scheduled for all photon beam qualities available before they could be used ...

CERN Document Server

A simple analytical/numerical model has been developed for computing the evolution, over periods of up to a few hours, of the current and temperature profile in the upper layer of the ocean. The model is based upon conservation laws for heat and momentum, and employs an eddy diffusion parameterisation which is dependent on both the wind speed and the wind stress applied at the sea surface. Other parameters such as the bulk-skin surface temperature difference and CO$_2$ flux are determined by application of the Molecular Oceanic Boundary Layer Model (MOBLAM) of Schluessel and Soloviev. A similar model, for the current profile only, predicts a temporary increase in wave breaking intensity and decrease in wave height under conditions where the wind speed increases suddenly, such as, for example, during gusts and squalls. The model results are compared with ...

Energy Technology Data Exchange (ETDEWEB)

Numerical simulations are used to predict the migration of radionuclides from the disposal units at Material Disposal Area G through the vadose zone and into the main aquifer in support of a radiological performance assessment and composite analysis for the site. The calculations are performed with the finite element code, FEHM. The transport of nuclides through the vadose zone is computed using a three-dimensional model that describes the complex mesa top geology of the site. The model incorporates the positions and inventories of thirty-four disposal pits and four shaft fields located at Area G as well as those of proposed future pits and shafts. Only three nuclides, C-14, Tc-99, and I-129, proved to be of concern for the groundwater pathway over a 10,000-year period. The spatial and temporal flux of these three nuclides from the vadose zone is applied as a source term for the three-dimensional saturated zone ...

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)

With an objective to discuss applicability of resistivity imaging using electromagnetic migration to detection of underground cavity, an applicability test was carried out on calculation of a numerical model and measurement data. By using the numerical model, a calculation was performed on a hypothetical case that a cubic cavity with sides each at 40 m exists in a homogenous medium of 200 ohm-m, with the cavity top located 20 m below the ground surface. As a result, it was possible to structure in a very short calculation time an image of the cavity which cannot be identified by one-dimensional inverse analysis. In the case of this cavity, a center of the image was structured on its lower face. It was shown that a location to indicate the image must be changed according to difference in the measurement locations. In a test on data measured in an underground cavity in the city of Utsunomiya, Tochigi ...

CERN Document Server

To detect the gravitational-wave signal from binary neutron stars and extract information about the equation of state of matter at nuclear density, it is necessary to match the signal with a bank of accurate templates. We have performed the longest (to date) general-relativistic simulations of binary neutron stars with different compactnesses and used them to constrain a tidal extension of the effective-one-body model so that it reproduces the numerical waveforms accurately and essentially up to the merger. The typical errors in the phase over the $\\simeq 22$ gravitational-wave cycles are $\\Delta \\phi\\simeq \\pm 0.24$ rad, thus with relative phase errors $\\Delta \\phi/\\phi \\simeq 0.2%$. We also show that with a single choice of parameters, the effective-one-body approach is able to reproduce all of the numerically-computed phase evolutions, in contrast with what found when adopting a tidally corrected post-Newtonian ...

International Nuclear Information System (INIS)

Nuclear physics parameters of the Wolsung CANDU-PHW reactor are computed by use of the PHWCELL computer code that is an improved version of LATREP. The PHWCELL code mainly computes cell parameters of heavy water moderated reactors, and modeling scheme of heavy water reactor cell calculations has been developed with the PHWCELL computer code. The reactor operating conditions considered in the study are cold zero power (CZP) and hot full power (HFP) with equilibrium poison. The cell parameters are also computed as a function of fuel burnup and the numerical results are compared with the results in PSR of the Wolsung unit and in the previous study. (author).

Energy Technology Data Exchange (ETDEWEB)

This paper presents a reverse time migration (RTM) method for the migration of shot records in tilted transversely isotropic (TTI) media. It is based on the tilted TI acoustic wave equation that was derived from the dispersion relation. The RTM is a full depth migration allowing for velocity to vary laterally as well as vertically and has no dip limitations. The wave equation is solved by a tenth-order finite difference scheme. Using 2D numerical models, we demonstrate that ignoring the tilt angle will introduce both lateral and vertical shifts in imaging. The shifts can be larger than 0.5 wavelength in the vertical direction and 1.5 wavelength in the lateral direction.

International Nuclear Information System (INIS)

The numerical analysis of the experimental data on viscosity change by time in the period after fast heating and crystal-liquid phase transition for the copper nonequilibrium melts and the Fe_7_0Cr_1_0P_1_3C_7 glass-forming melt near their melting temperatures is accomplished through the method of correlation functions and the Rytov-Dimentberg criterion. The observed dependences have nonmonotonic oscillating character but they do not relate to auto-oscillations. The proposed theoretical model of the properties oscillations origination by relaxation in the nonequilibrium liquids explains them as manifestation of noise-induced transitions in the metastable area near the critical temperatures (melting, structural transitions)

British Library Electronic Table of Contents (United Kingdom)

Due to the low thermal conductivity of ceramics large temperature gradients are present through the powder particles during plasma spray deposition. As a result the particles often impinge at the substrate in a semi-molten form; which in turn substantially affects the final characteristics of the coating being formed. This study is dedicated to the novel modelling development and simulation of a semi-molten droplet impingement. The study examines the impingement process during impact, spreading and solidification of semi-molten zirconia. The simulation provides an insight to the heat transfer process during impact and solidification of a semi-molten powder particle and illustrates the freezing-induced break-up mechanism at the splat periphery.

Energy Technology Data Exchange (ETDEWEB)

We consider energy investment, when a choice has to be made between fossil fuel and biomass fired production technologies. A dynamic model is presented to illustrate the effect of the different degrees of input price uncertainty on the choice of technology and the timing of the investment. It is shown that when the choice of technology is irreversible, it may be optimal to postpone the investment even if it would otherwise be optimal to invest in one or both of the plant types. We provide a numerical example based on cost, estimates of two different power plant types. (author)

Energy Technology Data Exchange (ETDEWEB)

Dynamic response of pipelines buried in a back-filled rectangular trench in a semi-infinite medium has been investigated. The pipelines are modeled as long cylindrical shells of small thickness. By using the boundary integral representation and finite element method, we have studied the three-dimensional response to account for either pane P or SV wave incident at an arbitrary angle to the pipe-axis. In this paper numerical results are presented for the normal displacements, displacements along pipe-axis, and the hoop stresses in the pipe wall. It is shown that the response depends critically on the back-filled material as well as on the directions of propagation of the incident waves.

Energy Technology Data Exchange (ETDEWEB)

Numerical analysis of solar dish modified cavity receiver with Cone, CPC and Trumpet reflectors is presented. Three-dimensional modeling is carried out to estimate the convective and radiative heat loss from the receiver for different angles of inclination and operating temperatures. Incorporating reflectors in the modified cavity receiver for second stage concentration, the natural convection heat losses are reduced by 29.23, 19.81 and 19.16%, respectively. The receiver with the trumpet reflector has shown better performance as compared to other configurations. (orig.)

CERN Document Server

We provide a both qualitative and quantitative comparison among different approaches aimed to solve the problem of non-linear diffusive acceleration of particles at shocks. In particular, we show that state-of-the-art models (numerical, Monte Carlo and semi-analytical), even if based on different physical assumptions and implementations, for typical environmental parameters lead to very consistent results in terms of shock hydrodynamics, cosmic ray spectrum and also escaping flux spectrum and anisotropy. Strong points and limits of each approach are also discussed, as a function of the problem one wants to study.

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.

Science.gov (United States)

In the past few decades the need for improved nuclear reactor safety analyses has led to a rapid development of advanced methods for multidimensional thermal-hydraulic analyses. These methods have become progressively more complex in order to account for the many physical phenomena anticipated during steady state and transient Light Water Reactor (LWR) conditions. The advanced thermal-hydraulic subchannel code COBRA-TF (Thurgood, M. J. et al., 1983) is used worldwide for best-estimate evaluations of the nuclear reactor safety margins. In the framework of a joint research project between the Pennsylvania State University (PSU) and AREVA NP GmbH, the theoretical models and numerics of COBRA-TF have been improved. Under the name F-COBRA-TF, the code has been subjected to an extensive verification and validation program and has been applied to variety of LWR steady state and transient simulations. To enable F-COBRA-TF for industrial applications, ...

Energy Technology Data Exchange (ETDEWEB)

No abstract prepared.

Energy Technology Data Exchange (ETDEWEB)

MAGMA is a FORTRAN computer code designed to viscous flow in in situ vitrification melt pools. It models three-dimensional, incompressible, viscous flow and heat transfer. The momentum equation is coupled to the temperature field through the buoyancy force terms arising from the Boussinesq approximation. All fluid properties, except density, are assumed variable. Density is assumed constant except in the buoyancy force terms in the momentum equation. A simple melting model based on the enthalpy method allows the study of the melt front progression and latent heat effects. An indirect addressing scheme used in the numerical solution of the momentum equation voids unnecessary calculations in cells devoid of liquid. Two-dimensional calculations can be performed using either rectangular or cylindrical coordinates, while three-dimensional calculations use rectangular coordinates. All derivatives are approximated by finite ...

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

Energy Technology Data Exchange (ETDEWEB)

Numerical models based on finite-difference or finite-element methods are used by various research groups in studies of radon-222 transport through soil and building materials. Applications range from design of radon remediation systems to more fundamental studies of radon transport. To ascertain that results obtained with these models are of good quality, it is necessary that such models are tested. This document reports on a benchmark test organized by the EU project ERRICCA: European Research into Radon in Construction Concerted Action. The test comprises the following cases: 1) Steady-state diffusive radon profiles in dry and wet soils, 2) steady-state entry of soil gas and radon into a house, 3) time-dependent radon exhalation from a building-material sample. These cases cover features such as: soil heterogeneity, anisotropy, 3D-effects, time dependency, combined advective and diffusive transport ...

Energy Technology Data Exchange (ETDEWEB)

The 3-D model program was developed to apply a mise-a-la-masse method to linear current sources with optional shapes. Mise-a-la-masse method is a bipolar mapping survey method using current sources installed in wells, and recently has been widely used for direct exploration of geothermal reservoirs. This method is also used for monitoring underground fluid as electrode arrangement of fluid flow tomography for surveying underground fluid (geothermal fluid, underground water, petroleum). In the geothermal reservoir exploration, the casing pipes of wells are used as linear current sources, and measured data are processed as those based on vertical current source. In the largely inclined well, the inclination of current sources should be considered. The 3-D modeling program was developed by difference calculus using the theoretical potential equation and apparent resistivity based on linear current sources with optional 3-D shapes. The ...

International Nuclear Information System (INIS)

Full text.This paper describes a software tool for generation expansion planning based on dynamic programming, probabilistic production simulation, and environmental assessment. The aim is to determine the units needed to expand a given system in order to minimize either the cost or the environmental impact or some weighed function of the two. The problem of generation expansion planning is solved in stages using Tunnel Dynamic Programming (TDP) with Heuristic rules to limit the number of options analyzed. The production costing methodology is based on combining a probabilistic generation model known as the capacity outage table (COT) with the load duration curve (LDC) of the system to deduce a risk model from which the expected energy not supplied (EENS) is estimated. The generation model is built recursively using an efficient numerical convolution procedure and is combined at each step with the LDC ...

International Nuclear Information System (INIS)

A numerical investigation of the mixing of gaseous uranium and hydrogen inside an open-cycle gas core nuclear rocket engine (spherical geometry) is presented. The gaseous uranium fuel is injected near the centerline of the spherical engine cavity at a constant mass flow rate, and the hydrogen propellant is injected around the periphery of the engine at a five degree angle to the wall, at a constant mass flow rate. The main objective is to seek ways to minimize the mixing of uranium and hydrogen by choosing a suitable injector geometry for the mixing of light and heavy gas streams. Three different uranium inlet areas are presented, and also three different turbulent models (k-var-epsilon model, RNG k-var-epsilon model, and RSM model) are investigated. The commercial CFD code, FLUENT, is used to model the flow field. Uranium mole fraction, axial mass flux, and ...

Energy Technology Data Exchange (ETDEWEB)

In this work the development of a third order scheme of finite differences centered in mesh is presented and it is applied in the numerical solution of those diffusion equations in multi groups in stationary state and X Y geometry. Originally this scheme was developed by Hennart and del Valle for the monoenergetic diffusion equation with a well-known source and they show that the one scheme is of third order when comparing the numerical solution with the analytical solution of a model problem using several mesh refinements and boundary conditions. The scheme by them developed it also introduces the application of numeric quadratures to evaluate the rigidity matrices and of mass that its appear when making use of the finite elements method of Galerkin. One of the used quadratures is the open quadrature of 4 points, no-standard, of Newton-Cotes to evaluate in approximate form the elements of the rigidity ...

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

International Nuclear Information System (INIS)

Momentum transfer through an air-water interface was studied in a horizontal rectangular channel through experimental and numerical analysis. The liquid flow Reynolds numbers, based on hydraulic diameters, were 11500 and 17000. The air flow Reynolds number varied from 0 to 3260. In experiments, velocity distribution was measured by hot-film and hot-wire probes. The effect of the liquid flow Reynolds number and the direction of interfacial shear stress was investigated by performing the cocurrent and countercurrent experiments. The modified k-#epsilon# model was able to reproduce these experimental results qualitatively. The range of applicability of the modified k-#epsilon# model was also discussed. (author).

Science.gov (United States)

The integrity and safety of beam-like structures are dependent in part on their boundary conditions which can vary with time due to damage or aging. Structural health monitoring of such structures should therefore include attention to boundary conditions. Where the boundary conditions can be represented by a lumped spring then the identification of associated stiffness parameter values may be a means to quantifying the integrity of the support. This paper investigates such a method for identifying the equivalent translational and rotational stiffness of a constrained tapered beam-like structure. An analytical model of a beam of tapered width and thickness is adopted as a simplified representation of a tower-like structure. The model is used to explore in what scenarios natural frequencies and/or nodal points might be sufficiently sensitive to changes in support conditions to be measurable indicators of damage. The method is evaluated by Monte ...

British Library Electronic Table of Contents (United Kingdom)

The condensation reaction of furfural (F) on acetone (Ac) gives a high added value product, the 4-(2-furyl)-3-buten-2-one (FAc), used as aroma in alcohol free drinks, ice, candies, gelatines and other products of current life. This synthesis valorises the residues of sugar cane treatment since furfural is obtained by hydrolysis of sugar cane bagasse followed by vapor training extraction. In the face of numerous and complex reactions involved in this synthesis, it is very complicated to define the kinetic laws from exact stoichiometry. A solution allowing to cope the problem consists in identifying an appropriate stoichiometric model. It does not attempt to represent exactly all the reaction mechanisms, but proposes a mathematical support to integrate available knowledge on the transformati...

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.

British Library Electronic Table of Contents (United Kingdom)

Remotely sensed ocean color data and numerical modeling have been used to study the phenology of both spring and fall phytoplankton blooms (FPBs) in the Nova Scotian Shelf (NSS)-Gulf of Maine (GoM) region. The ocean color data reveal a general pattern of westward progression of the spring phytoplankton bloom (SPB), and an eastward progression of the FPB in the NSS-GoM region. The spatial pattern of mean chlorophyll concentration in spring is similar to that in fall, with a lower concentration in the NSS and higher in the GoM. Interannually, there is a weak but significant tendency for years with earlier (delayed) SPBs to be followed by delayed (earlier) FPBs, but the mean chlorophyll concentrations during SPBs are not correlated with those during FPBs. The interannual variability of SPB ti...

British Library Electronic Table of Contents (United Kingdom)

The offshore and inshore wave energy resource in Asturias (N Spain) is studied using wave buoy data and a hindcast dataset spanning 44 years (1958-2001). Offshore average wave power and annual wave energy values are found to exceed 30 kW/m and 250 MWh/m, respectively, at 7 of the 11 study sites. This substantial resource is characterised in terms of the sea states involved. Most of the energy is provided by IV quadrant waves with significant wave heights between 2 m and 5 m and energy periods between 11 s and 13 s. After analysing the offshore resource, numerical modelling is used to investigate the inshore wave patterns. A coastal wave model is validated with wave buoy data and applied to three case studies representative of storm, winter and summer conditions. Inshore wave energy concent...

Energy Technology Data Exchange (ETDEWEB)

The authors consider the dynamics of interacting elastic disks in the plane. This is an experimentally realizable two-dimensional model of dry granular flow where the stresses can be visualized using the photoelastic effect. As the elastic disks move in a vacuum, they interact through collisions with each other and with the surrounding geometry. Because of the finite propagation speed of deformations inside each grain it can be difficult to capture computationally even simple experiments involving just a few interacting grains. The goal of this project is to improve our ability to simulate dense granular flow in complex geometry. They begin this process by reviewing some past work, how they can improve upon previous work. the focus of this project is on capturing the elastic dynamics of each grain in an approximate, computationally tractable, model that can be coupled to a molecular dynamics scheme.

British Library Electronic Table of Contents (United Kingdom)

On the basis of a numerical model, the temperature and liquid fraction of spray-formed H13 tool steel are calculated as a function of time. Results show that a preheated substrate at the appropriate temperature can lead to very low porosity by increasing the liquid fraction in the deposited steel. The calculated cooling rate can lead to a microstructure consisting of martensite, lower bainite, retained austenite, and proeutectoid carbides in as-spray-formed material. In the temperature range between the solidus and liquidus temperatures, the calculated temperature of the spray-formed material increases with increasing substrate preheat temperature, resulting in a very low porosity by increasing the liquid fraction of the deposited steel. In the temperature region where austenite decomposit...

British Library Electronic Table of Contents (United Kingdom)

Glycyrrhizin is a triterpene glycoside, a major active constituent of licorice (Glycyrrhiza glabra) root and numerous pharmacological effects like anti-inflammatory, anti-viral, anti-tumour and hepatoprotective activities has been attributed to it. In this study we evaluated the anti-inflammatory activities of glycyrrhizin in mice model of acute inflammation, carrageenan-induced pleurisy. We report here that glycyrrhizin (given at 10mg/kg i.p. 5min prior to carrageenan) exerts potent anti-inflammatory effects in this model. Injection of carrageenan into the pleural cavity of mice elicited an acute inflammatory response characterized by fluid accumulation in the pleural cavity which contained a large number of neutrophils (PMNs) as well as an infiltration of PMNs in lung tissues and subsequ...

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

International Nuclear Information System (INIS)

The energy levels, wavelengths, oscillator strengths, Auger rates and level-to-level dielectronic recombination rate coefficients describing dielectronic recombination into excited levels of Ne-like titanium from F-like low-lying states are calculated. Our calculations are based on Dr. R.D. Cowan's semi-relativistic mass-velocity and Darwin corrections are included in the Hamiltonian, and the distorted-wave model is used for the calculation of free electron wavefunctions. In order to set the recombination rate coefficients on a level by level basis, in a manner compatible with detailed level population kinetics modelling of highly-stripped ions in plasma, the dielectronic recombination rate coefficients as a function of free electron temperatures are given in an analytical form, which is not only very convenient in practice, but also hopefully accurate compared with the exactly calculated numerical results. (orig.).

International Nuclear Information System (INIS)

This paper develops a fluid dynamic model using finite difference methods characterizing flow in phantoms simulating in vivo conditions and to compare those results with velocity encoded MR images. The phantom consisted of 1-inch (2.5-cm) tubing with semicircular insert and fluid with viscosity, T1, and T2 comparable to blood. Numeric solutions to Navier-Stokes equations for this system were obtained using finite difference methods, with velocity input function of zero at walls and parabolic at both ends. In resulting color raster (CR) images, color temperature represented velocity value. In velocity-encoded MR images acquired under the same flow conditions, phase is proportional to average velocity during application of flow-encoding gradients. Because these gradients are applied along one direction per acquisition, magnitude and direction of velocity are obtained.

Energy Technology Data Exchange (ETDEWEB)

An analytical model is proposed for the design and simulation of in-situ fracture of deep oil sand formations. This model is based on the finite element variational principle in conjunction with special empirical modules to characterize in-situ oil sands behavior. A computer code by the name of SANFRAC was developed to handle the dynamic fracture of formations induced by explosives. Simulation of hydraulic fracture processes can be treated by the same code as special cases using the quasi-static analysis. Numerical case studies by the SANFRAC code indicate that extensive horizontal fracture can be achieved by dynamic loads with proper fracture starters configured at the injection well. The unique advantage of the dynamic fracturing technique over the hydraulic fracture methods is also demonstrated by these studies.

Science.gov (United States)

Atmospheric and ionospheric perturbations associated with the acoustic-gravity waves (AGW) with typical frequencies of a few hertz -millihertz are considered. These events may be caused by the influence from space and atmosphere as well as by oscillations of the Earth surface and other near-surface phenomena. The surface sources include long-period oscillations of the Earth's surface, earthquakes, explosions, thermal heating, seisches and tsunami waves. The wavelike phenomena manifest themself as travelling disturbances of air (in the atmosphere) and of electron density (in the ionosphere). Travelling ionospheric disturbances (TIDs) are well detected by radio physical methods. AGW generation by near-surface sources is modeled by the numerical solution of the equation of geophysical fluid dynamics for different sources in two-dimensional non-linear dissipative compressible atmosphere. The numerical calculations are based on ...

International Nuclear Information System (INIS)

It is of utmost importance to have a computer code in order to analyze how different parameters (like test duration time) affect the unavailability of safety systems of nuclear. In this context, a study was performed in order to evaluate the model employed by the FRANTIC computer code, which performs detailed calculations on the contribution to the system unavailability originated by hardware failures, component tests and repairs, aiming at considering the influence of different test schemes on the system unavailability. It was shown, by means of the results attained that the numerical model used by the FRANTIC code and the analytical model proposed by APOSTOLAKIS and CHU (4) give unavailability values much similar when the component tests are supposed to be perfect. When a test is supposed to be imperfect (that is, when it may induce a test is supposed to be imperfect (that is, when it may induce a ...

Energy Technology Data Exchange (ETDEWEB)

In a scenario of aqueous transport from a high-level radioactive waste repository, the concentration of radionuclides in water in contact with the waste constitutes the source term for transport models, and as such represents a fundamental component of all performance assessment models. Many laboratory experiments have been done to characterize release rates and understand processes influencing radionuclide release rates from irradiated nuclear fuel. Natural analogues of these waste forms have been studied to obtain information regarding the long-term stability of potential waste forms in complex natural systems. This information from diverse sources must be brought together to develop and defend methods used to define source terms for performance assessment models. In this manuscript examples of measures of radionuclide release rates from spent nuclear fuel or analogues of nuclear fuel are presented. Each example ...

Energy Technology Data Exchange (ETDEWEB)

To address the need for integrated models of advanced power generation systems and for improved analysis techniques, the US Department of Energy (DOE), through the Morgantown Energy Technology Center (METC), has supported the development of performance, emissions, and cost models of several advanced power generation systems, including integrated gasification combined cycle (IGCC), externally-fired combined cycle (EFCC) and pressurized fluidized bed combustion (PFBC) concepts. Simultaneously, numerical methods for simulation of uncertainties and optimization of process flowsheets have been implemented in the DOE`s public version of ASPEN. The new process models and computational capabilities have been applied to demonstrate the benefits of quantitative approaches to dealing with uncertainty and for optimizing technologies in the face of uncertainty. This paper focuses on modeling and ...

International Nuclear Information System (INIS)

Modeling free-electron laser (FEL) oscillators requires calculation of both the light-beam interaction within the undulator and the light propagation outside the undulator. We have developed a paraxial optical propagation code that can be combined with various existing models of gain media, for example, Genesis 1.3 for FELs, to model oscillators with full paraxial wave propagation within the resonator. A flexible scripting interface is used both to describe the optical resonator and to control the codes for propagation and amplification. To illustrate its capabilities, we numerically investigate two significantly different FEL oscillators: the free-electron laser for infrared experiments (FELIX) system and the vacuum-ultraviolet (VUV)-FEL oscillator of the proposed high-gain fourth generation light source. For the FELIX system, we find that diffraction losses are a considerable part of the single-pass ...

CERN Document Server

The present work concerns the numeric modeling of the sea-wave instability under the effect of the logarithmic-wind profile at hurricane conditions. Non-linear effects, such as wave breaking, foam production, etc. Powell et al. (2003), Shtemler et al. (2010) are ignored. The central point of the study is the calculation of the wave growth rate, which is proportional to the fractional input energy from the wind to the wave exponentially varied with time. The present modeling demonstrates that the Miles-type model applying Charnock's formula for roughness to the hurricane-wind parameters underestimates the growth rate from 5 to 40 times as compared with the model employing the roughness and friction velocity adopted from experimental data for hurricane winds.1 This occurs due to Charnock's formula fails at large wind speeds. The stability characteristics found on the base of the ...

CERN Document Server

We study a model for a protein searching for a target, using facilitated diffusion, on a DNA molecule confined in a finite volume. The model includes three distinct pathways for facilitated diffusion: (a) sliding - in which the protein diffuses along the contour of the DNA (b) jumping - where the protein travels between two sites along the DNA by three-dimensional diffusion, and finally (c) intersegmental transfer - which allows the protein to move from one site to another by transiently binding both at the same time. The typical search time is calculated using scaling arguments which are verified numerically. Our results suggest that the inclusion of intersegmental transfer (i) decreases the search time considerably (ii) makes the search time much more robust to variations in the parameters of the model and (iii) that the optimal search time occurs in a regime very different than that found for ...

Science.gov (United States)

We compared experimental wave aberrations in pseudophakic eyes with aspheric intraocular lenses (IOLs) to simulate aberrations from numerical ray tracing on customized computer eye models using corneal topography, angle ?, ocular biometry, IOL geometry, and IOL tilt and decentration measured on the same eyes. We found high correlations between real and simulated aberrations even for the eye with only the cornea, and these increased on average when the IOL geometry and position were included. Relevant individual aberrations were well predicted by the complete eye model. Corneal spherical aberration and horizontal coma were compensated by the IOL, and in 58.3% of the cases IOL tilt and decentration contributed to compensation of horizontal coma. We conclude that customized computer eye models are a good representation of real eyes with IOLs and allow understanding of the relative contribution of optical, ...

Energy Technology Data Exchange (ETDEWEB)

Mixed-convection phenomena can occur within liquid-dominated geothermal reservoirs due to interactions of injected flows, or ground-water flows, with the buoyancy-induced fluid motion. This problem was studied experimentally and numerically for the case of opposing flows about a vertical heat source in a liquid-saturated porous medium. The ratio of the Rayleigh number (Ra) to the Peclet number (Pe) was identified as the nondimensional parameter which characterizes the relative influence of buoyancy-driven to pressure-gradient-driven fluid motion. The transition from mixed to forced convection was numerically determined to be (Ra/Pe) approx. = -0.5, where the minus sign denotes superimposed downflow. Agreement between measured and predicted thermal-field results showed that the finite-element code of Gartling and Hickox (1982 a,b) can be used to model low-temperature (single-phase) geothermal reservoirs throughout the ...

International Nuclear Information System (INIS)

The numerical calculation of oscillatory wavefunctions occurring in the theory of electron scattering from positive ions is considered, the aim being to develop methods for evaluating the oscillatory solutions using a logarithmic grid. Thus, the permitted step-size for the solution depends only upon the potential, enabling calculations at different energies to be performed using the same discretisation. Two similar generalisations of the W.K.B. transformation for second order ordinary differential equations are made and applied to the radial Dirac equations. Iterative schemes for the solution of the resulting non-linear amplitude equations are discussed, one for use in the asymptotic region and the other a Newton scheme for continuing the solution closer to the origin. Using these techniques a reliable and efficient program for calculating relativistic scattering phase shifts from ionic potentials has been written. Relativistic Optical Model ...

Energy Technology Data Exchange (ETDEWEB)

Squeeze film dampers (SFD) are devices utilized to control vibrations of the shafts of high-speed rotating machinery. The SFD - squirrel cage combination is probably the most used system for tuning the stiffness and damping of the supports for rotors installed on ball bearings. Squeeze film dampers are essentially hydrodynamic bearings which contain the ball bearings housings of ball-bearings supported shafts. Consequently, the oil film within the SFD are influenced only by the precession and nutation of the shaft, that is the flow of the oil within the damper is not directly influenced by the spin of the rotor. However, in the classical theory, the flow in the thin film is also governed by the Reynolds equation. In this paper, some of the limits of the classical theory of the SFD are discussed and theoretical and experimental studies, which illustrate the ideas presented herein, are presented as well. The orbits of an unbalanced rotor that is supported by a ...

Energy Technology Data Exchange (ETDEWEB)

We intend to solve equations governing turbulent plane-vertical isotherm and non isotherm jets by taking into account inflow conditions at the exit of the nozzle. The analysis is focused on the influence of these conditions on this type of flow. Two cases are considered (uniform and parabolic velocity and temperature profiles). A finite difference scheme is developed to solve the governing equations. This numeric model allows us to show that the region of fully developed regime begins much nearer the nozzle for the turbulent case than for the laminar flow case. Indeed, the turbulence increases the mixing between the incoming gas from the nozzle and the ambient fluid, and consequently the size of the potential core zone decreases. The results are compared to other works introducing mathematical variables based on the energy conservation for the case of the mixed convection and the momentum conservation for the forced convection, which allows the ...

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)

Uranium enrichment by gas centrifugation relies on the pressure diffusion of the isotopes when their gaseous mixture is submitted to the pressure diffusion in the strong radial pressure gradient due to the centrifugal force in a very rapidly rotating cylinder. To compute the isotope separation in a gas centrifuge, the viscous compressible Navier-Stokes equations are solved for the mixture, then the diffusion equations are solved for the isotopes. As far as the separation performance is concerned, the models deal with the steady state of the axi-symmetric gas motion in the rotating cylinder. The numerical tools are briefly presented: first, the MOLINA code which numerically solves the axi-symmetric gas motion in the specific conditions of rapid rotation by use of a change of variables suited to the high density gradient and by use of a Marker And Cell method; then, the diffusion equations solver which uses a Finite Volume ...

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 compared with Les ones. Large Eddy Simulations have been ...

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 core numerical benchmark and the 3D model of the Argonauta research reactor, built in ...

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 good agreement with ...

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

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 experiments in laboratory. Focus is given in this paper on the ...

International Nuclear Information System (INIS)

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

Energy Technology Data Exchange (ETDEWEB)

Gas turbines are used where large amounts of energy is needed, typically as engines in aircraft, ferries and power plants. From an efficiency point of view it is desirable to increase the service temperature as much as possible. One of the limiting factors is then the maximum allowable metal temperatures in the turbine stages, primarily in the blades of the first stage, that are exposed to the highest gas temperatures. Specially designed materials are used to cope with these severe conditions, such as the nickel base superalloy IN792. In order to be able to design the components for higher temperatures and tighter tolerances, a detailed understanding and computationel models of the material behaviour is needed. The models presented in this work have been developed with the objective of being physically well motivated, and with the intention of avoiding excessive numbers of parameters. The influence of the parameters should also be as easy as ...

CERN Document Server

(abridged) The accretion disk in AGN is expected to produce strong outflows, in particular a UV-line driven wind. Despite providing a good fit to the data, current spectral models of the X-ray spectrum of AGN observed through an accretion disk wind are ad-hoc in their treatment of the properties of the wind material. In order to address these limitations we adopt a numerical computation method that links a series of radiative transfer calculations, incorporating the effect of a global velocity field in a self-consistent manner (XSCORT). We present a series of example spectra from the XSCORT code that allow us to examine the shape of AGN X-ray spectra seen through a wind, for a range of velocity and density distributions, total column densities and initial ionization parameters. These detailed spectral models clearly show considerable complexity and structure that is strongly affected by all these factors. The presence of ...

Energy Technology Data Exchange (ETDEWEB)

The suitability of the turbulent combustion flamelets model in order to predict the index of NO{sub x} production in turbulent flames of hydrogen diffusion is analyzed. In the flamelet approach, the turbulent flame is equivalent to a group of laminar flames submitted to a mechanical stretching which generates a chemical disequilibrium. This effect can be described by the stretching or by the scalar dissipation ratio. A numerical modeling is performed in order to evaluate the advantages of both approaches and to compare the behaviour of the NO{sub x} emission index with the experiments of Chen and Driscoll. This study shows that predictions of NO{sub x} emission indexes have a correct behaviour with respect to the Damkoehler number only when the scalar dissipation ratio is used as a parameter to describe the chemical state outside equilibrium. Predictions of the flamelet models are improving when the ...

Energy Technology Data Exchange (ETDEWEB)

Combined energy plants operating on the Brayton-Rankine cycle have gained wide acceptance as the most efficient energy conversion system. The technology is well proven and numerous research activities are underway to develop it further past its present plateau of performance. In developed and emerging economies alike, the current situation in the global economic market has caused the liberalisation of the energy industry in some countries, and in others it is causing a rethink of economic and energy policies. A result is that electric utility companies that have enjoyed a monopoly in power generation and distribution are seeing this privilege being eroded. This is opening up opportunities for independent power production and, consequently, for gas/steam turbine plants because of the technical merits of the technology and its ability to fill the vacuum being created in the energy industry. An economic model has been developed for carrying out ...