Sur l'origine des chiffres arabes A. Boucenna 1 From the pagination of an Algerian Arabic manuscript of the beginning of the 19th century,we rediscover the original shape that the Arabicnumerals had before passing in Europe and underwent the transformation that gave the modern Arabicnumerals. This original shape,whose use disappeared completely, proves that these numerals have their origin in the Arabic letters. Contrary to what some hypotheses pretend, particularly those that present them as drifting of Indian characters, the 10 Arabicnumerals that we use are, nothing else, 10 Arabic letters more or less modified and taken in the "Abjadi" order. The hypothesis of the Indian origin of the Arabicnumerals is ...
A large body of evidence indicates clear relationships between number and space processing in healthy and brain-damaged adults, as well as in children. The present paper addressed this issue regarding atypical math development. Adults with a diagnosis of dyscalculia (DYS) during childhood were compared to adults with average or high abilities in mathematics across two bisection tasks. Participants were presented with Arabic number triplets and had to judge either the number magnitude or the spatial location of the middle number relative to the two outer numbers. For the numerical judgment, adults with DYS were slower than both groups of control peers. They were also more strongly affected by the factors related to number magnitude such as the range of the triplets or the distance between t...
Few data exist about the impact of computerized prescribing systems on outpatient medication errors (MEs) and adverse drug events (ADEs). We compared the rates of MEs and ADEs in handwritten sites versus...Full Text Available
Artificial Neural Networks (ANNs) are parallel distributed processing machines. The unique characteristics of ANNs are: Fault tolerance, robustness, plasticity and generalization. These offer great potential in many AI applications such as character recognition. Handwritten character recognition is an intrinsically interesting problem, but the difficulties of this task are the many variations in the characters. A robust new incremental learning method, which combines supervised and unsupervised learning paradigms implemented by the Functional Link Net, is illustrated with experimental results. Clustering, based on unsupervised learning, classifies the input data into several categories. The supervised learning paradigm then further classifies the data in the clustered categories.
The purpose of this study was to assess the utility of the Arabic translation of the Alcohol Use Disorders Identification Test (AUDIT) among a sample of male Muslim prisoners (N = 107) in Dubai, United Arab Emirates (UAE). Confirmatory factor analysis indicated a two-factor model to be the best fit of the data. Good internal reliability (a =.91) and predictive validity were also observed. While the limitations of the study are observed, these promising results suggest that the Arabic version of the AUDIT is a reliable and sound measure of alcohol use disorders among prisoners in the UAE.
Various physics applications involve the computation of the standard hypergeometric function {sub 2}F{sub 1} and its derivatives. Because it is not an intrinsic in the common programming languages, automatic differentiation tools will either differentiate through the code that computes {sub 2}F{sub 1} if that code is available or require the user to provide a hand-written derivative code. We present options for the derivative computation in the context of an ionization problem and compare the approach implemented in the Diamant library to standard methods.
An observational and modeling study was performed over the Arabian Gulf region to investigate the coastal circulations and aerosol transport in the area. Climatological data and observations from the United Arab Emirates' Unified Aerosol Experiment were used to develop a better understanding of the complex meteorological processes in the Arabian Gulf region. Climatological data suggests that sea breezes occur on more than 77 percent of days in all months of the year and land breezes occur on more than 70 percent of the days. The occurrence of the sea and land breeze circulations are higher (90-99 percent) during the summer months when large-scale weather patterns are quiescent. Measurements of aerosol concentrations taken during the Unified Aerosol Experiment(2004) experiment are used to investigate aerosols, namely dust, transported in the Arabian Gulf region. Vertical profiles of dust concentration along with vertical profiles of potential temperature and wind ...
BACKGROUND AND OBJECTIVE:Streptococcus pneumoniae is a major cause of mortality and morbidity in both developing and industrialized countries, especially...Full Text Available
... t.., tt:,tifv! th - cui,nt !shor tfal 1 in breaching capAbiii tN wi h i- 1!,-ia . ~ ... Before the 1967 war, the Israeli Training ... trained Arab opponents. ...
QUESTION: You've already had one trip where you went to deal with fallout from WikiLeaks. You are visiting a couple of places where the leaders were quoted saying some very...
OBJECTIVE: To determine the common types of injuries among children (0-14 years) in Al-Ain, United Arab Emirates (UAE). DESIGN: A retrospective descriptive hospital based study. SETTING: Al-Ain Medical...Full Text Available
Due to shortage of fresh water resources, the vegetation of the eastern region of the United Arab Emirates (UAE) has experienced a series of declines resulting from salinization of groundwater, which is the major source of irrigation. To assess these changes, field measurements combined with Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) based Soil Adjusted Vegetation Index (SAVI) were analysed. TM and ETM+ images from two dates, 1987 and 2000 were acquired to enable the computation of the greenness anomalies for three sites in the eastern region, Fujairah, Kalba and Hatta. The results show an overall increase in agricultural area, associated with a severe decrease in vegetation greenness and health conditions, particularly in the Kalba study area. The SAVI values decreased ...
HAMILTON is a computer code performing all algebraic operations necessary for an analytic determination of the power series of the Hamiltonian equations of motion in the electromagnetic fields with at least one plane of symmetry. It is written entirely in FORTRAN in order to achieve fast machine performance, a requirement which is essential due to the complexity of the equations of motion in higher orders. HAMILTON is considerably faster than common more versatile formula manipulators and uses noticeably less storage. Besides the mere solution of the equations of motion, HAMILTON also produces FORTRAN code compatible with the program COSY 5.0 allowing the computation of matrix elements of individual optical elements and their concatenation. The produced FORTRAN code is highly optimized and on average requires only 30% of the execution time of a handwritten comparable code. (orig.).
Clastic sediments of Albian age provide hydrocarbon source rocks, reservoirs and seals. Major reservoirs occur in sandstones to the north, and have a limited reservoir potential in offshore Qatar. Distal shales are a major seal in offshore Qatar, the United Arab Emirates and Oman. Local source-rock potential exists in the northern part of the Arabian Gulf and in the extreme SE of Saudi Arabia (in the Rub Al Khali Basin). (author)
...Events Archive Worldwide Africa Programme Arab Programme Asia and Pacific Programme Europe and NIS Programme Latin American and Caribbean Programme Least-Developed Countries South-South Cooperation Human Security Directory of Offices Research & Statistics Publications Research and Statistics Subsites Annual Report Evaluation Group Investment and Technology Promotion Policymaking Organs Legal Resources Goodwill Ambassadors Industrial Competitiveness and Trade Private Sector Development Environmental Management Energy and Climate Change Research and Statistics Financial Institutions Partnerships Thematic ...Events Archive Worldwide Africa Programme Arab Programme Asia and Pacific Programme Europe and NIS Programme Latin American and Caribbean Programme Least-Developed Countries South-South Cooperation Human Security Directory of Offices Research & Statistics Publications Research and Statistics Subsites Annual Report Evaluation Group ...
The papers presented at the symposium had covered three general areas in which isotopes could have been beneficially used. these areas are: -Water use and water use efficiency studies. -Ground water investigations -Water problems in the arab countries. The individual papers had dealt with these subjects: -Hydrological research in the arab countries by use of radioisotopes. -The perspectives of use of radioisotopes in hydrological studies in Syria. -Water use efficiency and sub-soil water studies. -Sea water inclusion in a coast el aquifers in Lebanon. -Irrigation requirements of crops in Lebanon as determined by a Neutron probe with reference to other methods. -The use of the neutron moisture meter and other methods of the determination of the evapotranspiration of maize. -Ground water investigations, dating and nuclear methods applied to hydrology. -Ground water investigation in Wa di El-Nat run, U.A.R. -Velocity distribution along the pumped ...
A number of herbal drugs are used in the Unani (Greco-Arab) System of Medicine for cardiovascular diseases. The herbs were analyzed by flame AAS and ICP-AES to determine if their therapeutic actions can be associated with the elements present in them. Cadmium, cobalt, chromium, copper, iron, potassium, magnesium, manganese, sodium, nickel, phosphorus, lead and zinc were some of the elements which play various roles in cardiovascular affections. An effort was made to correlate the role of these elements in cardiac diseases. (Auth.). 2 tabs., 32 refs.
This Publications Catalogue lists all sales publications of the IAEA published in 2004, 2005 and forthcoming in 2006. Most IAEA publications are issued in English, some are also available in Arabic, Chinese, French, Russian or Spanish. This is indicated at the bottom of the book entry. A complete listing of all IAEA priced publications is available on the IAEA's web site: http://www.iaea.org/books
These classroom guides, designed to accompany the daily CNN (Cable News Network) Newsroom broadcasts for the month of October 2000, provide program rundowns, suggestions for class activities and discussion, student handouts, and a list of related news terms. Top stories include: Chinese authorities detain Falun Gong protesters on Tiananmen Square and Pope John Paul II angers China by canonizing Chinese martyrs, U.S. presidential candidates prepare for the first presidential debate, U.S. presidential candidates debate education issues, Danville, Kentucky prepares for the first U.S. vice presidential debate, and Yugoslav protesters seize parliament (October 2-6); Yugoslavian president Vojislav Kostunica takes office, recent clashes between Arabs and Israelis concern neighboring Arab states, North Korea observes 55 years of communism with celebrations, Al Gore and George W. Bush meet in the second presidential debate, and Violence in the Mideast ...
Possible factors regulating phytoplankton variability in the Arabian (Persian) Gulf were analyzed on the basis of satellite observations and meteorological data (1997-2009), including remotely-sensed chlorophyll a concentration (CHL), sea surface temperature, wind, solar radiation, precipitation, and aerosols. Shallow waters of northwestern Gulf influenced by Shatt Al-Arab River discharge were more productive than open Gulf waters, although seasonal CHL patterns in this and other shallow regions looked unrealistic likely because the CHL signal was obscured by bottom reflection. Therefore our further analyses focused on the open Gulf waters, which show a subtropical seasonal CHL cycle with maximum in winter and minimum in spring-summer. This cycle, however, was decoupled from the seasonal e...
This paper summarizes the IHY and BSS activities in West Asia and their importance in many Arab countries, such as Algeria, Egypt, Iraq, Jordan, Kuwait, Qatar, Saudi Arabia, UAE, etc. BSS future plans for some of these countries are as follows: It is proposed by the astronomers from the Arabian Gulf Region to build the Gulf Observatory on top of Jabal Shams (2980 msl) which will have a 2-3 m optical telescope. Libya signed a contract with a French company for building an observatory which will have a 2-m optical robotic telescope. It is also proposed to rebuild the Iraqi National Astronomical Observatory (INAO) which was destroyed during the two wars. It is planned to build a 5-6 m optical telescope and a small solar telescope on the top of Korek mountain, which has excellent observing conditions.
The effects of time of seed maturation and dry seed storage and of light and temperature requirements during seed incubation on final germination percentage and germination rate were assessed for the invasive shrub Prosopis juliflora (Sw.) D.C., grown under desert environmental conditions of the United Arab Emirates (UAE). Seeds were collected from Fujira on the northern coast of the UAE at different times during the growing seasons (autumn, winter and spring) and were germinated immediately and after 8 months of dry storage under room temperature (20+-3degreeC). Seeds were germinated at three temperatures (15, 25 and 40degreeC) in both continuous light and darkness. The results showed significant effects for time of seed collection, seed storage, light and temperature of seed incubation a...
Author.Field, laboratory and market studies are a must for proper evaluation of natural resources for the mineral industry of lebanon. Sites selectively convenient to the existing major Lebanese ceramics industry centered in the Beqa'a region, were investigated as to their geology, geography and economic character. The raw materials are shales, mud stones, siltstones and other argillaceous rocks of Jurassic to Cretaceous age coming from selective sites in south and central Lebanon. The finished products include wall, floor, roofing tiles, pipes, sanitary ware, pottery and brick specimens. Differential thermal analysis, scanning electron microscopy, firing and physico-chemical tests and analyses characterized the raw materials into two major groups: the suitable are siliceous argillaceous rocks and unsuitable calcareous argillaceous rocks. The suitable group is divided into two varieties. The first is dominantly a disordered Kaolinite with low drying and firing values, low plasticity ...
The Upper Jurassic Hith Anhydrite is a major hydrocarbon seal in the Arabian Gulf region. Outcrops, core samples from the subsurface, and the literature indicate that the Hith Formation is composed mainly of anhydrite. In most locations where a section of the Hith Formation has been measured, this unit contains less than 20% carbonate much of which is in the form of thin laminations. This lack of carbonate, locally thick layers of salt, and the predominance of anhydrite favor a playa for the setting in which this sediment was accumulated. In fact, much of the Hith has the sedimentary characteristics of the Holocene Lake MacLeod playa of Western Australia, which is dominated by layers of gypsum and halite (what little carbonate that occurs is found in layers at the base of the section). Locally the Hith appears to have accumulated in a sabkha setting, particularly toward central Abu Dhabi where it pinches out into shallow-water, and peritidal carbonate. This sabkha setting is indicated ...
Research into methods for reasoning under uncertainty is currently one of the most exciting areas of artificial intelligence, largely because it has recently become possible to record, store, and process large amounts of data. While impressive achievements have been made in pattern classification problems such as handwritten character recognition, face detection, speaker identification, and prediction of gene function, it is even more exciting that researchers are on the verge of introducing systems that can perform large-scale combinatorial analyses of data, decomposing the data into interacting components. For example, computational methods for automatic scene analysis are now emerging in the computer vision community. These methods decompose an input image into its constituent objects, lighting conditions, motion patterns, etc. Two of the main challenges are finding effective representations and models in specific applications and finding efficient algorithms ...
The application of foam in dealing with underground fire is well known due to its smothering action by cutting off air feed to burning fuel as well as acting as coolant. Besides plugging air feed to fire, water could be virtually reached to the fire affected areas much beyond the jet range as underground galleries with low roof restrict jet range of water. This method also enables a closer approach of a fire fighting team by isolating the toxic gases and smoke with a foam plug. The paper describes the development of high expansion foam composition and its application technology in order that foam plug method can be suitably utilized for combating mine fires in India. Three compositions were recommended for generation of high expansion foam: (a) 0.5% sodium/ammonium lauryl sulphate, 0.15 to 0.2% sodium carboxy methyl cellulose, 0.1% booster; (b) 0.5% sodium/ammonium lauryl sulfate, 0.12 to 0.15% alkaline solution of gum arabic, 0.1 to 0.2% ferrous gluconate; and (c) ...
Diesel engines operating on gaseous fuels are commonly known as dual-fuel engines. In the present work, a single-cylinder, compression ignition, indirect injection research (Ricardo E6) engine has been installed at United Arab Emirates University for investigation of the exhaust emissions when the engine is operating as a dual-fuel engine. The influence of changes in major operating and design parameters, such as the concentration of gaseous fuel in the cylinder charge, pilot fuel quantity, injection timing and intake temperature, on the production of exhaust emissions was investigated. Diesel fuel was used as the pilot fuel, while methane or propane was used as the main fuel which was inducted in the intake manifold and mixed with the intake air. The experimental investigations showed that the poor emissions at light loads can be improved significantly by increasing the concentration of gaseous fuel (total equivalence ratio), employing a large pilot fuel quantity, ...
... Title : The Numerical Weather Prediction System at the Italian Air Force Weather Service: Impact of Non-Conventional Observations and Increased ...
The report gives a summary of important developments in numerical weather prediction and discusses their relationship to, and effect on, the practice of synoptic meteorology.
... been possible to prove the existence and uniqueness of a classical solution which tends to the well-known solution of the Stefan problem as the ...
Numerical calculations of static, switched, and radiofrequency (RF) electromagnetic (EM) fields considering the geometry and EM properties of the human body are used increasingly in MRI to explain...Full Text Available
This report presents a summary of highlights of the year 1978. It then discusses: Research and development in numerical analysis and prediction systems; Techniques, development and application of new products; Analysis-forecast systems in operation use du...
Pfister, Robert M. (Lamont Geological Observatory, Palisades, N.Y.), and Paul R. Burkholder. Numerical taxonomy of some bacteria isolated from Antarctic and tropical seawaters. J....Full Text Available
BackgroundControversy exists about how to effectively measure patient reported outcomes in IBS clinical trials. Pain numeric rating scales (NRS) are widely used in...Full Text Available
Animals in a variety of taxa discriminate between a greater quantity and a lesser quantity of the same object, an ability that is referred to as relative numerousness judgment. For example,...Full Text Available
Coke-fired cupola furnaces were improved and made suitable for the production of high-quality casting melts by numerous additional devices. Moreover, they were equipped with environmental protection systems with numerous dust separation devices and afterb...
The vortex-excited dynamics of a uniform pivoted cylinder in uniform and sheared flow was investigated experimentally. The experiments were numerically simulated using a diffusive Van der Pol oscillator model developed by Balasubramanian & Skop recently. Salient features of the experimental investigations and the numerical simulations are presented here. Comparisons between the experimentally recorded and numerically predicted structural response to vortex-excited vibrations, power spectral density measurements of near-wake velocity fluctuations and lock-in ranges are made. A comparison of the numerical predictions and the experimental data reveals good agreement.
The radiative transport problem is cast in integral form using a transport kernel. The transport kernel has an explicit representation in terms of a Feynman Path Integral over all paths between selected points in a volume. This representation is setup in detail. Numerical evaluation of this Path Integral is formulated with a Frenet-Serret based procedure for generating valid random paths, and with a numerical evaluation of the weight for each valid path. Very early sanity checks of a numerical implementation are reported. Approaches to optimization are identified. (authors)
A physical model of free-electron laser (FEL) amplifier with variable-parameter wiggler magnets for one-dimensional numerical simulation is presented and a numerical example is given. The wiggler parameters, efficiency of energy conversion between electron beam and laser field, laser intensity, phase-space distributions and energy spectrum of electrons are computed. The period of synchronous oscillation and saturation value of laser intensity agree with estimated one.
A boundary value problem for Laplace?s equation in a bounded two-dimensional domain filled with a piecewise homogeneous medium is considered. The boundary of the inhomogeneity is assumed to be unknown. The inverse problem of determining the inhomogeneity boundary and the solution of the equation given the solution and its normal derivative on the boundary of the domain is discussed. Numerical methods are proposed for solving the inverse problem, and the results of numerical experiments are presented.
The enhancement of longitudinal heat transfer by means of fluid pulsation in a pipe has been investigated analytically and numerically, including the transient state. The effects of pulsation amplitude, frequency, and pipe length on thermal properties such as effective thermal diffusivity and delay time are clarified. Their effects on numerical calculations are also presented and suggestions for efficient numerical calculations of this problem are made concerning the combination of parameters.
An numerical method was developed for measuring radiant temperature using a set of absorption-foils coupled to soft x ray streak camera, SXRSC, in laser plasma experiments. An numerical code, SCC, was given. An x-ray intensity ratio vs temperature was calculated for Parylene, C_3H_6, Mylar and Aluminium. A suitable range of the measuring temperature was discussed to some absorber.
This paper presents a case study in the design and implementation of a numerical weather prediction model on a supercomputer (CRAY-1). Following a historical introduction to the evolution of the model, the governing equations of the model are presented and the numerical solution of these forecast equations is described. A brief tutorial on the architecture of the CRAY-1 is presented with a discussion of how it affects the choice of algorithms and code design of the model. A summary of the advantages gained by use of the vector aspects of the CRAY-1 is included.
Numerical simulations of compressible flows are commonly based on the Euler equations when effects of viscosity are thought to be negligible. These equations admit singular solutions, even in cases where the initial and boundary conditions are smooth. So-called ''Euler solvers'' rely on numerical dissipation, explicitly or implicitly present in the scheme, to regularize the problem, such that physical solutions are selected.
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.
Processing of symbolic sequences represented by mapping of symbolic data into numerical signals is commonly used in various applications. It is a particularly popular approach in genomic and proteomic sequence analysis. Numerous mappings of symbolic sequences have been proposed for various applications. It is unclear however whether the processing of symbolic data provides an artifact of the numerical mapping or is an inherent property of the symbolic data. This issue has been long ignored in the engineering and scientific literature. It is possible that many of the results obtained in symbolic signal processing could be a byproduct of the mapping and might not shed any light on the underlying properties embedded in the data. Moreover, in many applications, conflicting conclusions may arise due to the choice of the mapping used for numerical representation of symbolic data. In this paper, we present a ...
Full text.Changes in land use and urbanization yield more pressure put on limited soil and water resources, including the risk of pollution with toxic heavy metals. The study area lies in the Bekaa valley totaling about 12753 ha. The valley receives from the west torrential fan deposits and a mixture of colluvial and alluvial material. The principal soil classes are Fluvisols, Cambisols, Regosols, Vertisols and Luvisols. The area is populated and also the most important agricultural part of the plain. Agriculture in the plain is being practiced mainly with cash, field crops and vegetables. The western surrounding area is being used mainly for terraced fruit trees. This Arab-German Technical Cooperation Project (ACSAD-BGR) aimed, following the ISO standards and Eikman-Klocke recommendations, at investigating the nature of the extends of soil pollution by heavy metals in two pilot areas: The central Bekaa-Lebanon and Ghouta-Damascus. Different institutions cooperate ...
The MISR instrument aboard NASA's Terra satellite participated in the UAE-2 campaign, August-October 2004. This campaign represented a unique opportunity to study the complex aerosol situation in the Arabian Gulf region, in the context of a first-rate collection of aircraft and surface-based instruments, giving us the opportunity to do some groundbreaking satellite aerosol validation work. We aimed (1) to validate MISR aerosol retrieval results for dust and pollution particles over dark and light surfaces, and (2) to contribute regional maps of aerosol optical thickness (AOT) and particle micro-physical properties, giving spatial context to the field-instrument measurements, and moving toward a satellite-based regional aerosol climatology. The validation effort benefited from the combination of an instrumented aircraft and a regional network of surface-based sun photometers. We obtained high-quality sub-orbital data coincident with MISR overpasses on three days: September 01, 10, and ...
The present paper is concerned with the prediction method for the motions and tensions of mooring lines on slack-moored floating oceanic architectural building in coastal zone. The motions and the tensions of mooring lines were estimated by the time domain numerical simulation in the present method. The present method includes the nonlinearity of mooring system, slow varying wave drift force. The authors performed a time domain numerical simulation for motions and the tensions of mooring lines on floating barge in two component waves. Also, the results obtained by numerical simulation were compared with the experimental results. The validity of the present method was confirmed.
This paper presents a numerical simulation with the EUROPLEXUS fast dynamics software of a pipe whipping phenomenon occurring in the thermal hydraulic conditions of a loss of coolant accident in a PWR primary circuit. Different physical phenomena take place simultaneously during the rupture and the whipping of the pipe such as plasticity, contact, large displacements, two-phase flow regime and fluid structure interaction. Two kinds of numerical models - a simplified pipeline model and a mixed 1D/3D model - are considered and compared throughout modelling and computation. Numerical results are compared with experimental data belonging to the Aquitaine II test campaign.
We review an analysis of a consistent renormalization of the top and bottom quark/squark sector of the MSSM with complex parameters (cMSSM). Various renormalization schemes are defined, analyzed analytically and tested numerically in the decays Stop_2 -> Sbottom_i H^+/W^+ (i = 1,2). No scheme is found that produces numerically acceptable results over all the cMSSM parameter space, where problems occur mostly already for real parameters. Some numerical examples for Gamma(Stop_2 -> Sbottom_1 H^+) in our preferred scheme, "m_b, A_b DRbar" are shown.
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 representation of the system. (authors)
drical Enclosures Using a Numerical Ray Tracing Technique. NASA. TM-I02527, 1990 . Buckley, H.: Radiation from the Interior of a Reflecting Cylinder. Philos. ...
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
This bibliography lists unclassified publications of work done at the Los Alamos National Laboratory for 1988. The entries, which are subdivided by broad subject categories, are cross-referenced with an author index and a numeric index.
terms of 11) numerical models, 2D SPH simulations, and analytic models. The 1D model was useful for determining the evolution over very long timescales. ...
Brush fires, especially during the summer months, are hazards that occasionally threaten the leeward (Waianae) coast of Oahu in Hawaii. The prediction of future ... ...
SummaryA meshless numerical model is proposed to investigate shallow-water dam break flows in 1D open channels. The numerical model is to solve the shallow water equations (SWE) based on smoothed particle hydrodynamics (SPH). The concept of slice water particles (SWP) is adopted in the SPH-SWE formulation. The numerical sensitivity analysis is first performed to study the appropriate SWP number and variable smoothing length through dam break flows in an idealized 1D channel with dry/wet beds. Extensive validation by comparison with laboratory and field data is next conducted for four benchmark problems, including dam break flows through a rough flat channel, a rough bumpy channel with various downstream boundary conditions, a nonprismatic channel, and a realistic scale model of the Toce ri...
A numerical analysis method for melting/solidification phenomena has been developed to evaluate a feasibility of several candidate techniques in the nuclear fuel cycle. Our method is based on the eXtended Finite Element Method (X-FEM) which has been used for moving boundary problems. Key technique of the X-FEM is to incorporate signed distance function into finite element interpolation to represent a discontinuous gradient of the temperature at a moving solid-liquid interface. Construction of the finite element equation, the technique of quadrature and the method to solve the equation are reported here. The numerical solutions of the one-dimensional Stefan problem, solidification in a two-dimensional square corner and melting of pure gallium are compared to the exact solutions or to the experimental data. Through these analyses, validity of the newly developed numerical analysis method has been demonstrated. (author)
... that range from relatively low frequency strong ground motion modeling [Vidale et al., 1985] to high frequency reverse time migration of seismic ...
... being pulled out of water, the movement is restricted ... This heating depended on the particle size ... hydrometeor injection into the upper tropopause. ...
The quest to understand why we age has given rise to numerous lines of investigation that have gradually converged to include metabolic control by mitochondrial activity as a major player. That...Full Text Available
Numerous reports of chronic renal disease in patients who habitually use phenacetin-containing compounds to excess have aroused considerable controversy over the possible relationship between phenacetin,...Full Text Available
... these features and yields results which are two orders of magnitude faster than those available using conventional numerical ray tracing techniques ...
ground strike hazards 3) Advancements in the initialization of numerical weather prediction models through better identification of deep convection 4) Improved routing of...
We develop and investigate numerically a thermodynamically consistent model of two-dimensional multicomponent vesicles in an incompressible viscous fluid. The model is derived using an energy...Full Text Available
Numerous bacteria from a salt-marsh environment are capable of degrading malathion, an organophosphate insecticide, when supplied with additional nutrients as energy and carbon sources. Seven isolates...Full Text Available
Jan 7, 2010 ... Deforestation. The environmental implications of the word are as numerous as the syllables. And scientists like Jim Collatz have the job of ...
(Lactuca scariole),morning glory (Ipomea sp.), broomsedge (Andropogon virginicus ) and numerous composites. Because this area is leveed -and managed it ...
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.
Carboxylesterases (CE) are ubiquitous enzymes that hydrolyze numerous ester-containing xenobiotics, including complex molecules, such as the anticancer drugs irinotecan (CPT-11) and capecitabine...Full Text Available
Numerical results on aeroelastic standard configurations are presented. The methods used for two dimensional configurations include the small perturbations approach for inviscid flow, coupling methods for unseparated flow, coupling methods for unseparated or separated, laminar or turbulent boundary layers, and a numerical solution of the Euler equations for inviscid flow. The three dimensional configurations are studied by the transonic small disturbance approach. The detailed results are given.
Parameters of the vortex air flow in an isothermal model of the cyclone chamber are studied. The effect of the air velocity in nozzles on the values and character of the distribution of pressure and tangential and longitudinal velocities of air in the chamber is determined. Numerical modeling of the flow is performed, and it is shown that numerical calculation by the k-w turbulence model is in agreement with experimental data.
The structure of polydispersed populations of unilamellar vesicles is studied by small-angle neutron scattering for three types of lipid systems, namely, single-, two-and four-component vesicular systems. Results of the numerical analysis based on the separated-form-factor model are reported.
In this paper the stress due to torsion and the torsional buckling behavior of bellows are investigated through both theoretical and numerical analyses. The numerical solutions agree with the theoretical results very well in a wide range of dimensions and the number of convolutions. Based on these results, the framework of the design rule related to the torsion of bellows is proposed. (J.P.N.).
We study by means of Quantum Monte Carlo simulations based on the Worm Algorithm the low temperature (down to T = 0.05 K) properties of parahydrogen clusters comprising up to 40 molecules. Three different intermolecular interactions are employed: the Silvera-Goldman, the Buck and the Lennard-Jones potential. Despite important discrepancies observed in the numerical estimates of energy and superfluid fraction, the mechanism by which clusters melt at low T is independent of the particular choice of the potential, whose only effect is to alter the temperature scale.
Results are given of numerical simulation of three-dimensional pattern of flow of a two-phase steam-water mixture in the house of a PGV-1000 horizontal steam generator obtained using the BAGIRA best-estimate thermohydrodynamic computer codes. The space distributions of velocities and local void fractions in the steam generator housing for different modes of operation of power-generating unit are calculated and compared with available experimental data.
Results are given of numerical simulation of three-dimensional pattern of flow of a two-phase steam-water mixture in the house of a PGV-1000 horizontal steam generator obtained using the BAGIRA best-estimate thermohydrodynamic computer codes. The space distributions of velocities and local void fractions in the steam generator housing for different modes of operation of power-generating unit are calculated and compared with available experimental data.
Numerical errors of PALLAS calculation due to spatial mesh sizes are examined for a typical deep penetration shielding problem of isotropic incident fission neutrons penetrating a 200-cm-thick water slab. The exponential approximation for the source spatial distribution to solve the transport equation based on the direct integration method is verified to be effective for radiation transport in attenuating medium, while the linear approximation employed in the old PALLAS code is found to produce large errors for calculations with large mesh sizes.
A numerical calculation procedure applicable to cavity-type receiver configurations and flow conditions was developed. Flow visualization experiments were performed, and experimental measurements of quantities valuable for the development of the numerical calculation procedure were made. The investigation is focussed on a configuration which is strongly two-dimensional in the mean flow structure (but turbulent in a truly three-dimensional sense). (LEW)
Tthe acceleration of a liquid from a cylindrical container by the products of instantaneous detonation of a superimposed explosive charge (implosion method) is studies numerically in a two-dimensional formulation. Analytical formulas for the asymptotic velocities of liquid expansion and acceleration of the container body were obtained using numerical calculations in a one-dimensional approximation. The effect of the fracture of the body on the velocity of liquid expansion was studied in a two-dimensional formulation.
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.
A three-dimensional numerical simulation of mammatus-like clouds is presented. A portion of a cirrus outflow anvil cloud is simulated including cloud ice and snow microphysical representations. The simulated mammatus clouds appear in a cellular pattern and are compared with the few available previously published physical observations of mammatus. Copyright Copyright 2006 Royal Meteorological Society
In this paper, finite-difference time-domain method is used to research the transient electromagnetic pulse (EMP) coupling to the cylinder object which is shielded in another shielding cavity with a hole numerically. The different coupling effectiveness of three transient EMP is simulated. The distribution of the electric field in the shielding cavity and the current on the cavity are researched. The shielding effectiveness is evaluated. (authors)
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 ...
Finite-difference acoustic-wave modeling and reverse-time depth migration based on the full wave equation are general approaches that can take into account arbitrary variations in velocity and density and can handle turning waves as well. However, conventional finite-difference methods for solving the acoustic- or elastic-wave equation suffer from numerical dispersion when too few samples per wavelength are used. The flux-corrected transport (FCT) algorithm, adapted from hydrodynamics, reduces the numerical dispersion in finite-difference wavefield continuation. The flux-correction procedure endeavors to incorporate diffusion into the wavefield continuation process only where needed to suppress the numerical dispersion. Incorporating the flux-correction procedure in conventional finite-difference modeling or reverse-time migration can provide finite-difference solutions with no numerical dispersion even ...
A numerical optimization technique is used to obtain low-energy momentum transfer, j = 0 [yields] 2 rotational and v = 0 [yields] vibrational sections from measured electron swarm data for parahydrogen. The downhill simplex algorithm is used to find cross sections that represent the best numerical fit to the measured electron drift velocity and characteristic energy over a range of E/N. These results, which are in excellent agreement with published cross sections derived using traditional swarm analysis techniques, demonstrates the feasibility of using automated computational algorithms for swarm analysis involving the estimation of multiple cross sections. (Author).
We study analytically and numerically the stability of the pressure-less, viscously spreading accretion ring. We show that the ring is unstable to small non-axisymmetric perturbations. To perform the perturbation analysis of the ring we use a stretching transformation of the time coordinate. We find that to 1st order, one-armed spiral structures, and to 2nd order additionally two-armed spiral features may appear. Furthermore, we identify a dispersion relation determining the instability of the ring. The theoretical results are confirmed in several simulations, using two different numerical methods. These computations prove independently the existence of a secular spiral instability driven by viscosity, which evolves into persisting leading and trailing spiral waves. Our results settle the question whether the spiral structures found in earlier simulations of the spreading ring are numerical artifacts or genuine ...
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)
In this paper, numerical simulation of a concentric tube heat exchanger is presented to determine the convective heat transfer coefficient and friction factor in a smooth tube. Increasing the convective heat transfer coefficient can increase heat transfer rate in a concentric tube heat exchanger from a given tubular surface area. This can be achieved by using heat transfer augmentation devices. This work constitutes the initial phase of the numerical simulation of heat transfer from tubes employing augmentation devices, such as twisted tapes, wire-coil inserts, for heat transfer enhancement. A computational fluid dynamics (CFD) simulation tool was developed with CFX software and the results obtained from the simulations are validated with the empirical correlations for a smooth tube heat exchanger. The difficulties associated with the simulation of a heat exchanger augmented with wire-coil inserts are discussed. (author)
This paper presents a numerical study on the formation history of coalbed methane (CBM) reservoir in the southeast edge of Ordos Basin, China. The coal seams studied belong to the Late Palaeozoic coal-bearing series. These coal seams have a burial history and experienced the process of subsidence, rapid subsidence alternated with uplift and then uplift, sequentially, and underwent the geothermal actions at normal, extremely high, and then normal temperatures, respectively. Coal organic matter of the coal seams matured in the Triassic Period and in the Late Jurassic to Early Cretaceous Period. The results from numerical simulation reveal that CBM reservoir evolution history can be classified into five stages, namely primary, initial, stagnant, active and dissipative stages. In the first (pr...
This paper presents results, numeric and experimental, due to installation operation of a Risers Uphold Sub-Surface Buoy, (BSR). This kind of installation developed by PETROBRAS is unique in the world. The work of BSR installation was based on a numeric pre analysis to verify the system and determine the main parameters to be experimentally verified. The second phase of the work was the experimental analysis in a deep water ocean basin. s. The work describes the BSR and their main accessories, the experimental environment and the model constructed in aluminum in a 1:12 scale and the main results. (author)
Turbulent premixed methane-air combustion in a cylindrical chamber is numerically simulated considering radiation effect. Reaction rates are considered as minimum rates between Arrhenius rates and eddy break up rates. A five step reduced mechanism is used. Turbulent modeling is done via standard k-{epsilon} model imposed by empirical inlet boundary conditions. Source terms of energy equation consist of reaction rates and radiation effects. The discrete ordinate method (DOM) is employed to solve the radiative transfer equation (RTE) and the weighted sum of gray gas model (WSGGM) is imposed to consider radiation effect of non-gray gases. The results indicate that in the case of turbulent combusting flows, the effect of radiation of gases can affect the temperature and species concentrations. The numerical results obtained considering radiation effect are closer to the experimental data than that of the case without radiation effect. (author)
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-? model imposed by empirical inlet boundary conditions. Source terms of energy equation consist of reaction rates and radiation effects. The discrete ordinate method (DOM) is employed to solve the radiative transfer equation (RTE) and the weighted sum of gray gas model (WSGGM) is imposed to consider radiation effect of non-gray gases. The results indicate that in the case of turbulent combusting flows, the effect of radiation of gases can affect the temperature and species concentrations. The numerical results obtained considering radiation effect are closer to the experimental data than that of the case without radiation effect.
Baxter (Australas J Philos 79:449?464, 2001) proposes an ingenious solution to the problem of instantiation based on his theory of cross-count identity. His idea is that where a particular instantiates a universal it shares an aspect with that universal. Both the particular and the universal are numerically identical with the shared aspect in different counts. Although Baxter does not say exactly what a count is, it appears that he takes ways of counting as mysterious primitives against which different numerical identities are defined. In contrast, I defend the idea?suggested, though not quite endorsed, by Baxter himself?that counts are independent dimensions of numerical identity. Different ways of counting are explained by the existence of these different sorts of identity (i.e., counts)...
We present a new set of 3.5 Post-Newtonian equations in which Newtonian hydrodynamics is coupled to the nonconservative effects of gravitational radiation emission. Our formalism differs in two significant ways from a similar 3.5 Post-Newtonian approach proposed by Blanchet (1993, 1997). Firstly we concentrate only on the radiation-reaction effects produced by a time-varying mass-current quadrupole $S_{ij}$. Secondly, we adopt a gauge in which the radiation-reaction force densities depend on the fourth time derivative of $S_{ij}$, rather than on the fifth, as in Blanchet's approach. This difference makes our formalism particularly well-suited to numerical implementation and could prove useful in performing fully numerical simulations of the recently discovered $r$-mode instability for rotating neutron stars subject to axial perturbations.
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 affects the temperature distribution in the ...
The control volume finite element method (CVFEM) was developed to combine the local numerical conservation property of control volume methods with the unstructured grid and generality of finite element methods (FEMs). Most implementations of CVFEM include mass-lumping and upwinding techniques typical of control volume schemes. In this work we compare, via numerical error analysis, CVFEM and FEM utilizing consistent and lumped mass implementations, and stabilized Petrov-Galerkin streamline upwind schemes in the context of advection-diffusion processes. For this type of problem, we find no apparent advantage to the local numerical conservation aspect of CVFEM as compared to FEM. The stabilized schemes improve accuracy and degree of positivity on coarse grids, and also reduce iteration counts for advection-dominated problems.
Heterogeneous anisotropic diffusion problems arise in the various areas of science and engineering including plasma physics, petroleum engineering, and image processing. Standard numerical methods can produce spurious oscillations when they are used to solve those problems. A common approach to avoid this difficulty is to design a proper numerical scheme and/or a proper mesh so that the numerical solution validates the discrete counterpart (DMP) of the maximum principle satisfied by the continuous solution. A well known mesh condition for the DMP satisfaction by the linear finite element solution of isotropic diffusion problems is the non-obtuse angle condition that requires the dihedral angles of mesh elements to be non-obtuse. In this paper, a generalization of the condition, the so-called anisotropic non-obtuse angle condition, is developed for the finite element solution of heterogeneous anisotropic diffusion problems. ...
Proteinaceous aggregation occurs through self-assembly-- a process not entirely understood. In a recent article [1], an analytical theory for amyloid fibril growth via secondary rather than primary nucleation was presented. Remarkably, with only a single kinetic parameter, the authors were able to unify growth characteristics for a variety of experimental data. In essence, they seem to have uncovered the underlying allometric laws governing the evolution of filament elongation simply from two coupled non-linear ordinary differential equations (ODEs) stemming from a master equation. While this work adds significantly to our understanding of filament self-assembly, it required an approximate analytical solution representation. Here, we show that the same results are found by purely numerical means once a straightforward and reliable numerical solution to the set of ODEs has been established.
Ultra-short-pulse reflectometry is studied by means of the numerical integration of a one-dimensional full-wave equation for ordinary modes propagating in a plasma. The numerical calculations illustrate the potential of using the reflection of ultra-short-pulse, microwaves as an effective probe of the density profile even in the presence of significant density fluctuations. The difference in time delays of differing frequency components of the microwaves can be used to deduce the density profile. The modification of the reflected pulses in the presence of density fluctuations is examined and can be understood based on considerations of Bragg resonance. A simple and effective profile-reconstruction algorithm using the zero-crossings of the reflected pulse and subsequent Abel inversion is demonstrated. The robustness of the profile reconstruction algorithm in the presence of a sufficiently small amplitude density perturbation is assessed.
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...
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)
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)
The phosphatidylinositol pathway is implicated in the regulation of numerous cellular functions and responses to extracellular signals. An important branching point in the pathway is the phosphorylation...Full Text Available
The side effects of cancer therapy on normal tissues limit the success of therapy. Generation of reactive oxygen species (ROS) has been implicated for numerous chemotherapeutic agents including doxorubicin...Full Text Available
BackgroundChanges in ionic concentration have a fundamental effect on numerous physiological processes. For example, IP3-gated thapsigargin sensitive intracellular calcium...Full Text Available
Mating signals are often directed at numerous senses and provide information about species identity, gender, receptiveness, individual identity and mate quality. Given the diversity of colourful body...Full Text Available
BackgroundThe epidemiology of avian influenza viruses (AIVs) in gulls is only partially known. The role of the world's most numerous gull species, the black-legged kittiwake (Rissa...Full Text Available
Man-made endocrine-disrupting chemicals (EDCs) present a threat to biodiversity, even in remote areas. To date, numerous wildlife species have been affected by EDCs in the environment, but it is likely...Full Text Available
An ultrasmall spot size scanning laser ophthalmoscope has been developed that employs an annular aberration-corrected incident beam to increase the effective numerical aperture of the eye thereby reducing...Full Text Available
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 model, together with a brief description of algorithms incorporated in the newly ...
Lithium isotope separation by displacement chromatography is studied using fundamental principles; the equations are derived assuming theoretical stages in lithium adsorption bands. The concentration profiles in the band are calculated numerically under unsteady state.
The electric explosion of a conductor being part of an electric cable is studied. It is ascertained that the electric properties of a polymeric dielectric depend nonlinearly on the maintenance of the ultradisperse metal filler. In the course of the numerous heating-cooling processes of an insulator, a streamlining of the structure of a polymeric composite occurs.
The intracellular concentration of protein may be as high as 400 mg per ml; thus it seems inevitable that within the cell, numerous protein-protein contacts are constantly occurring. A basic biochemical...Full Text Available
BackgroundSurface electromyography (sEMG) signals have been used in numerous studies for the classification of hand gestures and movements and successfully implemented in the position...Full Text Available
Although there have been numerous studies investigating subjective spatial impression in rooms, only a few of those studies have addressed the influence of visual cues on the judgment of auditory measures....Full Text Available
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
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.
Adipocyte differentiation is a complex developmental process that involves the coordinated interplay of numerous transcription factors. PPARγ has emerged as a master regulator of adipogenesis...Full Text Available
The Neurofibromatosis-2 (NF2) tumor suppressor merlin negatively regulates cell proliferation in numerous cell types. We have previously shown that the NF2...Full Text Available
The impacts of climate change on crop productivity are often assessed using simulations from a numerical climate model as an input to a crop simulation model. The precision of these predictions reflects...Full Text Available
... His work has been recognised with numerous honours, including Fellow of the Australian Academy of Science; International Research Scholarships from the Howard Hughes Medical Institute ; a Wellcome Trust Fellowship; the Glaxo-Wellcome Australia Medal; the Australian Academy of Science Gottschalk Medal; the Boehringer Mannheim Medal;...
SummaryThe numerous applications of hydrogelators have led to rapid expansion of this field. In the present work we report the facile synthesis of amphiphilic hydrogelators having a quaternary...Full Text Available
This work presents an experimental and numerical study of pulsated Dean flow, three-dimensional pulsatile flow in a curved pipe. The numerical study is performed by CFD code (Fluent 6) in which a pulsated velocity field is imposed as an inlet condition. The experimental setup involves principally a 'Scotch-yoke' pulsatile generator and a 90o bend. Laser Doppler Velocimetry (LDV) measurements have shown that the Scotch-yoke generator produces pure sinusoidal instantaneous mean velocities with a mean deviation of 3%. Visualizations by laser-induced fluorescence (LIF) and velocity measurements, coupled with the numerical results, have permitted analysis of the evolution of the swirling secondary flow structures that develop along the bend during the pulsation phase. These measurements were made for a range of stationary Reynolds number (300?Rest?1200), frequency parameter (1??=r0(?/?)1/2max,osc/Ust). We observe satisfactory ...
Allergic asthma is characterized by airway hyperresponsiveness, inflammation, and a cellular infiltrate dominated by eosinophils. Numerous epidemiological studies have related the exacerbation of allergic...Full Text Available
The continuing rapid increase in available computing power has not reduced the importance of efficient methods of optical system assessment for automatic lens design. On the contrary, the new capabilities simply show that truly automatic optical design wi...
This paper reports on Monte Carlo path integral simulations of a Bianchi Type IX cosmology which yield a wave function which reflects the known chaotic dynamics of the classical system and predicts a significant probability for a correlation between large universe volume and large anisotropy.
Operation of an X-ray spectrometer based on a spherical variable-line-spacing (VLS) grating is analyzed using dedicated ray-tracing software allowing fast optimization of the grating parameters and...Full Text Available
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 solution, and the instability amplitude increases exponentially with decreasing mesh sizes. This work shows ...
Stable and accurate numerical analytical method even at high Rayleigh numbers is desired, and in addition flexibility and economical efficiency are very important for numerical analysis. Unsteady flows of natural convention in a square cavity are investigated using the GSMAC (generalized and simplified marker and cell) finite-element method at high Rayleigh numbers from 10 {sup 6} to 10 {sup 8}. Validities of the multi-pass algorithm and BTD (balancing tensor diffusivity), which are highly accurate solutions for finite-element method, are investigated by introducing them to the GSMAC finite-element method. As the result, it was found that the multi-pass algorithm has little effect for the GSMAC finite-element method. The steady solutions at the Rayleigh numbers of 10 {sup 5} and 10 {sup 7} agreed well with the bench mark solutions and numerical solutions. Complex transient phenomena at Rayleigh numbers higher than 10 {sup ...
BackgroundNumerous microRNAs (miRNAs) have heterogeneous ends resulting from imprecise cleavages by processing nucleases and from various non-templated nucleotide additions. The...Full Text Available
The estimation of numerical values of the mean distance of closest approach of ions, a, of alkaline-earth metal ion salts in aqueous solutions, determined from activity coefficients, as well as from different theoretical approaches, is presented and discussed.
BackgroundGrapevine is subjected to numerous pests and diseases resulting in the use of phytochemicals in large quantities. The will to decrease the use of phytochemicals leads to...Full Text Available
BackgroundThe epidermal growth factor receptor (Egfr) with its numerous ligands has fundamental roles in development, cell differentiation and physiology. Dysfunction of the receptor-ligand...Full Text Available
provides global ocean-surface wind fields with more structures than numerical weather prediction (Geophys. Res. Lett., Vol 25, No. 6, 761-764). The NSCAT data sets was produced...
Lepidoptera, butterflies and moths, is the second largest animal order and includes numerous agricultural pests. To facilitate comparative genomics in Lepidoptera, we isolated BAC clones containing...Full Text Available
Although cross-resistance in houseflies to the organophosphates has eliminated numerous potentially useful compounds from field use, the ”subgroup” specificity of this phenomenon has...Full Text Available
Information circulars are published from time to time under the symbol INFCIRC/... for the purpose of bringing matters of general interest to the attention of all Members of the Agency. The present revision contains INFCIRCs published up to mid-August 1994. A complete numerical list of information circulars is reproduced with their titles in the Annex.
Skeletal formation is an essential and intricately regulated part of vertebrate development. Humans and mice deficient in Growth and Differentiation Factor 6 (Gdf6) have numerous...Full Text Available
An analysis is made of the manner in which the process of primordial black-hole formation and the subsequent accretion of gas depend on the equation of state. On the assumption that the process is spherically symmetric, the problem is solved numerically.
Numerous studies have noted that the evolution of new enzymatic specificities is accompanied by loss of the protein's thermodynamic stability (ΔΔG), thus suggesting...Full Text Available
BackgroundDespite numerous treatments available for deteriorated cutaneous wound healing such as a diabetic foot, there is still the need for more effective therapy. Adipose-derived...Full Text Available
A numerical technique is developed for estimating water quality violation frequencies due to pollutant discharges from urban areas during combined sewer overflow events. The first four moments of in-stream pollutant concentration are found by integrating a pollutant loading - wat...
IntroductionTyrosine phosphorylated signal transducer and activator of transcription 3 (pStat3) is expressed in numerous cancers and is required for mediating tumorigenesis. Autocrine...Full Text Available
In recent years, the National Toxicology Program (NTP) has selected numerous glycol ethers for testing in routine laboratory mammals to ascertain the magnitude of their ability to injure the conceptus....Full Text Available
BackgroundRed ripe tomatoes are the result of numerous physiological changes controlled by hormonal and developmental signals, causing maturation or differentiation of various fruit...Full Text Available
Motivation: Cell tracking is an important method to quantitatively analyze time-lapse microscopy data. While numerous methods and tools exist for tracking cells in 2D time-lapse images,...Full Text Available
Clostridium perfringens produces numerous toxins, which are responsible for severe diseases in man and animals. Delta toxin is one of the three hemolysins released by a number of C....Full Text Available
In this college-level introductory text the authors outline the fundamental principles of meteorology. The text is divided into three sections: meteorology, regionalization of climate, and climate change. Numerous charts, graphs, and photographs are displayed with each subject. A glossary of terms is also provided at the end of the book.
Organisms belonging to the genus Streptomyces produce numerous important secondary metabolites and undergo a sophisticated morphological differentiation program. In many instances these...Full Text Available
Stem cell therapies for neurodegenerative disorders require accurate delivery of the transplanted cells to the sites of damage. Numerous studies have established that fluid injections to the hippocampus...Full Text Available
The comprehensive survey explains the indications and the specific value of the modern tomographic modalities with respect to the various diagnostic issues. Numerous reproductions of CT, MRI, US scans illustrate the aspects discussed. (orig./CB)
Piping in nuclear power plants is vital to the proper operation and safety of these facilities. To assure safety in the unlikely event of a pipe break, it is necessary to evaluate the consequences from the resulting whipping pipe on neighboring components and structures. Numerical simulations allow for rapid evaluation of these consequences. Before simulations can be accepted, however, the methodology and computer codes must be validated against experimental results. This paper uses a probabilistic approach to validate pipe whip simulations against limited experimental results. Probabilistic analysis software was developed and coupled to existing deterministic finite element software. An example of a whipping pipe impacting against a reinforced concrete slab was simulated. The described probabilistic approach was used to validate the numerical simulations. The conclusions obtained were that the numerical simulations of ...
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
Numerous xenobiotics, including therapeutics agents, are substrates for bioactivation to electrophilic reactive intermediates that may covalently modify biomolecules. Selective estrogen receptor...Full Text Available
Despite numerous investigations, the correlation between alcohol consumption and prostate cancer risk remains uncertain. This report investigated the association between alcohol use and prostate cancer...Full Text Available
MRI plays an increasingly important role for assessment of patients with chronic liver disease. MRI has numerous advantages, including lack of ionizing radiation and the possibility of performing...Full Text Available
Neurons are born throughout adulthood in the hippocampus and show enhanced plasticity compared to mature neurons. However, there are conflicting reports on whether or not young neurons contribute...Full Text Available
The state of chromatin (the packaging of DNA in eukaryotes) has long been recognized to have major effects on levels of gene expression, and numerous chromatin-altering strategies—including...Full Text Available
Aberrant transglutaminase (TG) activity has been implicated in the pathology of numerous diseases including Huntington disease and Alzheimer disease. To fully characterize the role of TGs in...Full Text Available
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.
A modelling approach is presented that appears to have the potential to provide reliable assessments of pollution concentration in coastal zones and complex terrain, where the assumptions behind current operational air quality models often are inadequate. With the use of a numerical mesoscale model (NMM), physically consistent flow fields can be predicted, thereby providing higher spatial and temporal resolution in the meteorological fields than would be available from a limited number of observation points. These predictions are used to calculate mean trajectories of pollutant parcels, as well as to provide quantitative estimates of pollution concentration using two techniques.
An equation of state of neutron star matter with strange baryons has been obtained. The effects of the strength of hyperon-hyperon interactions on the equations of state constructed for the chosen parameter sets have been analysed. Numerous neutron star models show that the appearance of hyperons is connected with the increasing density in neutron star interiors. The performed calculations have indicated that a change of the hyperon-hyperon coupling constants affects the chemical composition of a neutron star. The obtained numerical hyperon star models exclude a large population of strange baryons in the star interior.
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.
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.
Asymptotic and numerical analysis of relaxation self-oscillations in a three-dimensional system of Volterra ordinary differential equations that models the well-known Belousov reaction is carried out. A numerical study of the corresponding distributed model-the parabolic system obtained from the original system of ordinary differential equations with the diffusive terms taken into account subject to the zero Neumann boundary conditions at the endpoints of a finite interval is attempted. It is shown that, when the diffusion coefficients are proportionally decreased while the other parameters remain intact, the distributed model exhibits the diffusion chaos phenomenon; that is, chaotic attractors of arbitrarily high dimension emerge.
The two-dimensional WKB model has been the basis for several investigations of the quasi-Landau hydrogenic spectrum. Whereas other authors have used numerical integration, it is shown that the results can be obtained analytically in terms of elliptic integrals. The later are easily generated by even programmable pocket calculators, from which numerical results - which are in good agreement with experiments - are easily obtained. A further advantage of using elliptic integrals is that several strong-field problems can be shown to be special cases of a general potential whose energy and spacing are expressible in terms of them.
We present and explain numerical results illustrating the mechanism of a type of discontinuous bifurcation of a chaotic set that occurs in typical dynamical systems. After the bifurcation, the chaotic set acquires new pieces located at a finite distance from its location just before the bifurcation, and these new pieces were not part of a previously existing chaotic set. A scaling law is given describing the creation of unstable periodic orbits following such a bifurcation. We also provide numerical evidence of such a bifurcation for a nonattracting chaotic set of the H{acute e}non map. {copyright} {ital 1998} {ital The American Physical Society}
Necessary and sufficient conditions for a minimax solution to the Cauchy problem for the Hamilton-Jacobi-Bellman equation are obtained as viability conditions for classical characteristics inside the graph of this solution. Using this property, a representative formula for a one-dimensional conservation law in terms of classical characteristics is derived. An estimate of the numerical integration of the characteristic system is presented and errors of numerical realizations of representative formulas are determined for the conservation law and its potential equal to the minimax solution of the Hamilton-Jacobi-Bellman equation.
We study diffusion of Cosmic Rays (CRs) in turbulent magnetic fields using test particle simulations. Electromagnetic fields are produced in direct numerical MHD simulations of turbulence and used as an input for particle tracing, particle feedback on turbulence being ignored. Statistical transport coefficients from the test particle runs are compared with earlier analytical predictions. We find qualitative correspondence between them in various aspects of CR diffusion. In the incompressible case, that we consider in this paper, the dominant scattering mechanism occurs to be the non-resonant mirror interactions with the slow-mode perturbations. Perpendicular transport roughly agrees with being produced by magnetic field wandering.
Three-dimensional magnetohydrodynamic simulations are executed in a full toroidal geometry to clarify the physical mechanisms of the Internal Reconnection Event (IRE), which is observed in the spherical tokamak experiments. The simulation results reproduce several main properties of IRE. Comparison between the numerical results and experimental observation indicates fairly good agreements regarding nonlinear behavior, such as appearance of localized helical distortion, appearance of characteristic conical shape in the pressure profile during thermal quench, and subsequent appearance of the m=2/n=1 type helical distortion of the torus. (author)
Approximate methods are presented for solving nonlinear boundary value problems with an unknown (free) boundary, such as the Stefan problem in heat conductivity. Examples of this kind of problems in thermal physics, hydrodynamics, elasticity, and plasma physics are examined. Attention is given to the principal computational methods for solving stationary problems for second- and fourth-order elliptic equations, such as the domain transformation and penalty methods. To illustrate the capabilities of the methods presented here, examples of numerical solutions are presented for various specific applied problems. 255 references.
We present a numerical method for computing diffusive transport on a surface derived from image data. Our underlying discretization method uses a Cartesian grid embedded boundary method for computing the volume transport in region consisting of all points a small distance from the surface. We obtain a representation of this region from image data using a front propagation computation based on level set methods for solving the Hamilton-Jacobi and eikonal equations. We demonstrate that the method is second-order accurate in space and time, and is capable of computing solutions on complex surface geometries obtained from image data of cells.
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.
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.
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.
This paper focuses on a theoretical model that predicts the temperature increase of Magneto-rheological (MR) fluid dampers experiencing a sinusoidal input motion. A theoretical model is developed to estimate the temperature rise based on the non-linear behavior of the MR fluid damper. This model is solved numerically, and the numerical solution is compared with a known linear solution and experimental result in order to validate the accuracy of the model. Also, a non-dimensional form of the governing equations are developed to examine the key parameters. The non-dimensional terms show the effect of external and internal parameters on the trends of heat dissipation as well as heat generation within the MR fluid damper.
Numerical solutions for fully developed laminar flow in internally finned tubes with trapezoidal and triangular fin profiles were given with Finite Element Method (FEM): The heat transfer characteristics were obtained and compared under the boundary conditions of uniform heat flux, uniform wall temperature, and the third boundary condition with finite wall thermal conductivity considered. The numerical results show that boundary conditions have pronounced effects on the temperature field. Furthermore, a new mechanism on the heat transfer augmentation of internally finned tubes is proposed.
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
The goal of the book is to summarize those methods for evaluating Feynman integrals that have been developed over a span of more than fifty years. The book characterizes the most powerful methods and illustrates them with numerous examples starting from very simple ones and progressing to nontrivial examples. The book demonstrates how to choose adequate methods and combine evaluation methods in a non-trivial way. The most powerful methods are characterized and then illustrated through numerous examples. This is an updated textbook version of the previous book (Evaluating Feynman integrals, STMP 211) of the author. (orig.)
We present results for a complete set of polarization observables for jet production in lepton proton collision, where the final state lepton is not observed. The calculations are carried out in collinear factorization at the level of Born diagrams. For all the observables we also provide numerical estimates for typical kinematics of a potential future Electron Ion Collider. On the basis of this numerical study, the prospects for the transverse single target spin asymmetry are particularly promising. This observable is given by a certain quark-gluon correlation function, which has a direct relation to the transverse momentum dependent Sivers parton distribution.
Stabilized formation temperatures were estimated at different depths in 40 wells from the Los Humeros geothermal field, Mexico, using the Horner and the spherical radial flow (SRF) methods. The results showed that the Horner method underestimates formation temperatures, while the SRF method gives temperatures that are closer to the true formation temperatures. This was supported by numerical simulation of a combined circulation and shut-in period in several wells, and results for well H-26 are presented. Numerical reproduction of logged temperature is more feasible if an initial temperature profile based on the SRF method is employed instead of using an initial temperature profile based on the Horner method. (Author)
A new aiming point guidance (APG) law in 3-dimensional space is derived analytically. The guidance law is capable of guiding the missile (interceptor) to a straight collision course with respect to the target. The missile turning rate command proportional to the error angle between the current missile heading and the predicted aim point is calculated based on the predicted target position. Performance of the APG law is numerically evaluated and compared with the proportional navigation guidance (PNG) law. The result of numerical simulation shows that the APG law is superior to the PNG law in terms of miss distance, time of interception and final stage missile turning rate demand.
This work deals with the development of wind maps in the Venezuelan coastal zones (South America). The Navier-Stokes equations for real viscous incompressible flow are discretized via the finite element method and they are applied to a full 3D numerical model. Two regions were analyzed: the Margarita island and the Paraguana peninsula, both of them located in the Caribbean Sea at the North of Venezuela. These regions were modeled by using a three dimensional control volume, discretized through hexaedric eight-noded finite elements. Some wind maps obtained in this research are presented and discussed. (author)
This paper is concerned with the determination of a unique scaling parameter in complex scaling analysis and with accurate calculation of dynamics resonances. In the preceding paper we have presented a theoretical analysis and provided a formalism for dynamical resonance calculations. In this paper we present accurate numerical results for two non-trivial dynamical processes, namely, models of diatomic molecular predissociation and of barrier potential scattering for resonances. The results presented in this paper confirm our theoretical analysis, remove a theoretical ambiguity on determination of the complex scaling parameter, and provide an improved understanding for dynamical resonance calculations in rigged Hilbert space.
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.)
The formation of a quasiequilibrium beam distribution matched to an alternating-gradient quadrupole focusing lattice by means of the adiabatic turn-on of the oscillating focusing field is studied numerically using particle-in-cell simulations. Quiescent beam propagation over several hundred lattice periods is demonstrated for a broad range of beam intensities and vacuum phase advances describing the strength of the oscillating focusing field. Properties of the matched-beam distribution are investigated. In particular, self-similar evolution of the beam density profile is observed over a wide range of system parameters. The numerical simulations are performed using the WARP particle-in-cell code.
We propose a numerical method for resummation of perturbative series, which is based on the stochastic perturbative solution of Schwinger-Dyson equations. The method stochastically estimates the coefficients of perturbative series, and incorporates Borel resummation in a natural way. Similarly to the "worm" algorithm, the method samples open Feynman diagrams, but with an arbitrary number of external legs. As a test of our numerical algorithm, we study the scale dependence of the renormalized coupling constant in a theory of one-component scalar field with quartic interaction. We confirm the triviality of this theory in four and five space-time dimensions, and the instability of the trivial fixed point in three dimensions.
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 ...
DECOVALEX II project started in November 1995 as a continuation of the DECOVALEX I project, which was completed at the end of 1994. The project was initiated by recognizing the fact that a proper evaluation of the current capacities of numerical modelling of the coupled T-H-M processes in fractured media is needed not only for small scale, well controlled laboratory test cases such as those studied in DECOVALEX I, but also for less characterised, more complex and realistic in-situ experiments. This will contribute to validation and confidence building in the current mathematical models, numerical methods and computer codes. Four tasks were defined in the DECOVALEX II project: TASK 1 - numerical study of the RCF3 pumping test and shaft excavation at Sellafield by Nirex, UK; TASK 2 - numerical study of the in-situ T-H-M experiments at Kamaishi Mine by PNC, Japan; TASK 3 - review of current ...
Hyogo-ken Nanbu earthquake with a focus in the Akashi straits has given huge earthquake damages in and around Awaji Island and Kobe City in 1995. It is clear that the basement structure, which is steeply deepened at Kobe City from Rokko Mountains towards the coast, and the focus under this related closely to the local generation of strong ground motion. Generation process of the strong ground motion was discussed using 2D and 3D numerical simulation methods. The 3D pseudospectral method was used for the calculation. Space of 51.2km{times}25.6km{times}25.6km was selected for the calculation. This space was discretized with the lattice interval of 200m. Consequently, it was found that the basement structure with a steeply deepened basement, soft and weak geological structure thickly deposited on the basement, and earthquake faults running under the boundary of base rock and sediments related greatly to the generation of strong ground motion. ...
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.)
An outline is given of time-dependent wavepacket methods as applied to calculations of molecular collisions with solid surfaces. The methods reviewed include numerical integration algorithms for the time-dependent Schroedinger equation, semiclassical wavepacket treatments, and approximations that treat some of the degrees-of-freedom quantum-mechanically and others classically. The computational and numerical characteristics of these methods are discussed, with emphasis on their particular advantages and relevance in the context of certain molecule/surface scattering problems. For the semiclassical and mixed quantal-classical treatments, the approximation errors and their physical origins are discussed. For the quantum wavepacket techniques a numerical error analysis is presented. The computational efficiency of the various algorithms is considered and examined in the context of several applications. The main focus is on ...
In this paper, we calculated the radioactive concentration distribution of radioactive waste water, the temperature distribution of drained cooling water and the effect of implement from the Daya Bay Nuclear Power Plant on nearby waters range, discussed and analysed some problems of computational results and computation with Alternating Direction Implicit Method (ADI). The contents of the article included: the establishment of two-dimension tidal current equation, radioactive waste water pollutant dispersion equation and cooling water heat convection diffusion equation, the numerical difference calculation model of tidal current field, concentration field as well as temperature field, effect impingement with ADI method, numerical calculation results. The result of research showed that: when the Daya Bay Nuclear Power Plant is on normal operation and after the low level radioactive waste water and low temperature cooling water have been drained ...
The stability of a forced-flow cooled superconducting coil is investigated by use of the numerical simulation. The numerical code to integrate the simultaneous partial differential system composed of the 1 D hydrodynamic equations and the 1 D thermal conduction equation has been developed and stability margins are evaluated as functions of coolant mass flow rate, operation current and imposed magnetic field. The results of computations show that the stability margin is multi-valued with respect to these operation parameters, as expected from the experimental results. It is also shown that the appearance of the first unstable regime is closely related to the existance of the stagnant region located at the upstream side of the heated zone and that the second stable regime appears because the heat transfer is appreciably enhanced by the induced backflow due to the thermal expansion of coolant. 13 refs., 13 figs., 1 tab.
This dissertation reports on the application of numerical optimization techniques as applied to fuel cell simulation and design. Due to the "multi-physics" inherent in a fuel cell, which results in a highly coupled and non-linear behavior, an experimental program to analyze and improve the performance of fuel cells is extremely difficult. This program applies new optimization techniques with computational methods from the field of aerospace engineering to the fuel cell design problem. After an overview of fuel cell history, importance, and classification, a mathematical model of solid oxide fuel cells (SOFC) is presented. The governing equations are discretized and solved with computational fluid dynamics (CFD) techniques including unstructured meshes, non-linear solution methods, numerical derivatives with complex variables, and sensitivity analysis with adjoint methods. Following the validation of the fuel cell model in 2-D and 3-D, the ...
In the conventional visualization system for the computational results, only Japanese (Nihongo) Line Printer (NLP) was available to print two dimensional cross sectional plots of vector and scalar fields. To evaluate the phenomena, an analyst had to print many plots on the NLP. This task makes difficult to check the computational results immediately after the calculation. Recently, as the visualization tools, we introduced Micro AVS and Field View which are utilized widely in the scientific and the industrial fields. In order to show the numerical results on the visualization software, we constructed a post processing system which convert the results of the numerical code to 'intermediate files' which can be read by the visualization tools. As using this system, the examination of the numerical results can be executed on the display of the personal computer. Furthermore, the persuasive report and paper with high quality can ...
The paper deals with the modeling and the development of a numerical procedure for the analysis of shape-memory alloy (SMA) elements in order to predict the main features of SMA devices. A 3D SMA model in the framework of small strain theory is developed starting from the thermo-mechanical model proposed by Souza et al. (Eur J Mech A/Solids 17:789-806, 1998) and modified by Auricchio and Petrini (Int J Numer Methods Eng 55:1255-1284, 2002). The aim of this paper is to propose some more modifications to the original model, to derive its consistent 1D formulation, to clarify the mechanical meaning of the material parameters governing the constitutive model. A robust time integration algorithm is developed in the framework of the finite element method and a new beam finite element is proposed. Some numerical applications and a comparison with experimental data available in literature are carried out in order to assess the ...
The paper is concerned with the economic operation of power systems comprising pumped storage and thermal plants. In power systems, pumped storage plants may receive the energy required for pumping water either from the system in which they are integrated or from a neighbouring system. Therefore, two mathematical formulations of the problem are presented depending on the source of energy during pumping. In both cases, the problem is reduced to a non linear optimization problem. Two numerical techniques are suggested for the solution of this optimization problem, namely the Lagrangian Multipliers technique and the Rosen Gradient Projection technique. Numerical examples are presented to demonstrate the techniques. Although the Lagrangian Multipliers technique has the advantage of being a well known technique, however, computer wise, the Rosen Gradient Projection technique proves to be more efficient. The results of the ...
The improvement in the functions of the viscous flow calculation method VEGA-SHIP around a ship and the expansion of application range were described as the numerical flow simulation in ship and ocean engineering and at the same time application examples to the ocean engineering by the general-purpose flow simulation code FLOW-3D handling the non-steady flow with a free surface were introduced as the numerical simulation regarding such products as a water gate and a dam. In the VEGA-SHIP, water surface was handled as a fixed wall so that wave could not be calculated. Therefore, an algorithm for calculating wave on the water surface was added to the VEGA-SHIP and a calculation method simultaneously considering the creation of wave around the ship and viscosity was developed. The FLOW-3D was used to calculate the phenomenon where inside liquid moved greatly due to the oscillation of a tank and hit against and damaged the tank ceiling in the tank, ...
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 indicate that the ...
A fractographic and numerical approach is presented to analyze hydrogen-plasticity interactions in pearlitic steel and to elucidate the main hydrogen transport mechanism in this material under triaxial stress states produced by notches. Fractographic analysis showed that the microdamage produced by the hydrogen was clearly detectable by scanning electron microscopy (SEM), through a specific microscopic topography associated with hydrogen effects (tearing topography surface or TTS). Numerical computations obtained by using an elastic-plastic finite element program gave the progressive spreading of the plastic zone, closely associated with the movement of dislocations. In the majority of cases, the plastic zone (PZ) clearly exceeds the hydrogen affected region (TTS) and has no relation with it. In some tests, however, the hydrogen-induced micro-damage surpasses the only region in which there is dislocation movement, and in this case the net ...
A numerical study of laminar flow and heat transfer in an array of stacked rectangular plates is presented. The array is placed in a uniform stream, and the plates are subjected to a constant surface heat flux. This flow configuration is relevant to a number of practical heat transfer devices with finned surfaces. The computations were performed using a finite volume solution of the steady, two-dimensional Navier-Stokes equations and energy equation. A numerical scheme that reduces numerical diffusion is used to discretize the equations. The dominant feature of the flow is the separation, and subsequent reattachment of, the boundary layer, which takes place at Reynolds numbers greater than about 75. The separation first occurs downstream of the leading edge of the plate; then as Re increases, the separation point moves upstream and remains fixed at the leading edge, and the reattachment length increases linearly with Re. ...
Turbulent heat transfer performance of a fuel rod with three-dimensional trapezoidal spacer ribs for high temperature gas-cooled reactors was studied for various Reynolds numbers using an annular channel at the same coolant condition as the reactor operation, maximum outlet temperature of 1000 C and pressure of 4 MPa, and analytically by a numerical simulation using the k-#epsilon# turbulence model. The turbulent heat transfer coefficients of the fuel rod were 18-80% higher than those of a concentric smooth annulus at a region of Reynolds number exceeding 2000. On the other hand, the predicted average Nusselt number of the fuel rod agreed well with the empirical correlation obtained from the experimental data within a relative error of 10% with Reynolds number of more than 5000. It was verified that the numerical analysis results had sufficient accuracy. Furthermore, the numerical prediction could clarify quantitatively the ...
Lining shafts and placing them in a given strata needs extensive computation. Therefore, it is natural that the use of electronic computing and of numerical methods in this mining area started early. For shafts being sunk in non-stable water-bearing strata the deformations and stresses of the freezing wall are analysed. The complex system of an outer lining in the shaft sinking stage is designed with regard to the elastic support of the strata. In cases of strata-embedded systems the tubbing liner can be shaped as the loads require. Here analysis of safety against buckling is very important. If a sliding shaft is built as a composite construction, there are advantages in determining the distribution of the material by numerical methods. The analysis of stability is possible only by numerical methods and could not be solved otherwise. The Finite Element Method proves itself as a very general tool to solve the problem of ...
A Langmuir film is a molecularly thin film on the surface of a fluid; we study the evolution of a Langmuir film with two co-existing fluid phases driven by an inter-phase line tension and damped by the viscous drag of the underlying subfluid. Experimentally, we study an 8CB Langmuir film via digitally-imaged Brewster Angle Microscopy (BAM) in a four-roll mill setup which applies a transient strain and images the response. When a compact domain is stretched by the imposed strain, it first assumes a bola shape with two tear-drop shaped reservoirs connected by a thin tether which then slowly relaxes to a circular domain which minimizes the interfacial energy of the system. We process the digital images of the experiment to extract the domain shapes. We then use one of these shapes as an initial condition for the numerical solution of a boundary-integral model of the underlying hydrodynamics and compare the subsequent images of the experiment to the ...
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 1997. Despite lack of adequate field ...
We examine the effect of a threshold bias on the power spectrum and the bispectrum in an ensemble of numerical simulations (Gaussian initial perturbations with power law spectra P(k) \\sim k^n, n=+1, 0, -1, -2) and compare our results with theoretical predictions. Our simulations are evolved sufficiently that on the scale where we apply the threshold the rms fluctuation has developed significantly into the nonlinear regime. Thus, predictions based on perturbation theory do not necessarily apply. Nevertheless, we find our results for the power spectrum, biased power simply amplified by a numerical factor, follow predicted trends, far beyond the regime where perturbation theory is expected to be valid. We find that the biased bispectrum continues to follow the so-called hierarchical form, with reduced three-point amplitude Q \\approx 1 in the strongly nonlinear regime, independent of initial spectrum. In the quasi-linear perturbative regime the ...
We consider the supercooled Stefan problem with a general anisotropic curvature- and velocity-dependent boundary condition on the moving interface. We present numerical methods, based on an integral equation formulation and including a new algorithm for moving curves with curvature-dependent velocity. These methods compute a periodic interface with {ital O}({Delta}{ital t}) accuracy, where {Delta}{ital t} is the time step. Previous work has been limited to short time spans and achieved slightly less than {ital O}({Delta}{ital t}{sup 1/2}) accuracy. Accurate numerical results are seen to agree with the predictions of linear stability theory. This agreement has eluded previous authors, because their numerical methods suffered from grid effects and their linear stability theory was incorrect. We study the long-time evolution of an unstable interface. Our computations exhibit the beginnings of a sidebranching instability when ...
Results of three-dimensional laminar and standard K-#epsilon# turbulent numerical simulations of natural convection cooling of ten cubic aluminum blocks mounted on an insulated plate, facing a shrouding wall, are presented. This geometry is chosen so that comparison with experimental results is possible. The considered problem is of great practical importance because it simulates the case of heated electronic chips, mounted on printed board assemblies, which are frequently encountered in electronic industry applications. The problem is mathematically modeled by the three-dimensional conservation differential equations of mass, momentum, energy and turbulent kinetic energy and dissipation (for the turbulent flow model). IN this paper, these equations are numerically solved by a finite volume method and the laminar and turbulent results are compared to the experimental results obtained with similar parameters.
This article presents numerical study of 3-D thermosolutale mixed convection (TSMC) in horizontal rectangular channels. The contribution of this work is to characterize the travelling wave's appearance and to generalize the behavior of Poiseuille-Rayleigh-Benard (PRB) systems for a broad range of dimensionless parameters, which control the double diffusive mixed convection. The numerical results consist of analyzing the flow regimes of the steady longitudinal thermoconvectives rolls for the case of purely thermal mixed convection (TMC) and for both thermal and mass transfer (TSMC). The transition from opposed volume forces to cooperating ones at fixed Rayleigh (Ra), Reynolds (Re), and Lewis (Le) numbers, considerably affects the birth and the development of the longitudinal rolls (noted R/...
The direct injection of steam into a water pool is a method of heat transfer used in many process industries. The amount of research in this area however is limited to the nuclear industry, with applications relating to reactor cooling systems. Electrical resistance tomography (ERT), a low cost, non-invasive and which has high temporal resolution characteristics, can be used as a visualization tool for the resistivity distribution for the steam injection into water pool such as IRWST. In this paper, three dimensional resistivity distribution of the process is obtained through ERT using iterative Gauss-Newton method. Numerical experiments are performed by assuming different resistive objects in the water pool. Numerical results show that ERT is successful in estimating the resistivity distribution for the injection of steam in the water pool
Thermal-fluid flow analysis and demonstration test were performed for a spent fuel storage system. The commercial computational fluid dynamics (CFD) code, FLUENT was used for the numerical analysis. Effective thermal conductivities of a spent fuel assembly and a fuel basket were derived to optimize a thermal analysis model. Also, a porous model, which can simplify a complex configuration of a fuel assembly, was used in the thermal analysis. Demonstration test were performed to verify the thermal analysis method and procedure using a half scaled-down model and an electrically heated dummy fuel. The numerical analysis results were compared with the experimental data. Thermal analyses of the storage system were carried out for normal and off-normal conditions by using the verified analysis method.
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...
Abstract in english Numerical solutions for the effects of radiation on a MHD convective heat transfer past a semi-infinite porous plate with a magnetic field are obtained. It is assumed that the porous plate moves with a constant velocity in the direction of fluid flow, and the free stream velocity follows the exponentially increasing small perturbation law. The magnetic field acts perpendicular to the porous surface which absorbs the fluid with a suction velocity varying with time. The gov (more) erning equations for the flow are transformed into a system of nonlinear ordinary differential equations by perturbation technique and then are solved numerically by using the shooting method. The effects of the various parameters on the velocity, temperature profiles as well as the surface skin-friction and surface heat transfer are illustrated graphically.
A free-electron laser (FEL) two-beam accelerator (TBA) is proposed, in which the FEL interaction takes place in a series of drive cavities, rather than in a waveguide. Each drive cavity is 'beat-coupled' to a section of the accelerating structure. This standing-wave TBA is investigated theoretically and numerically, with analyses included of microwave extraction, growth of the FEL signal through saturation, equilibrium longitudinal beam dynamics following saturation, and sensitivity of the microwave amplitude and phase to errors in current and energy. It is found that phase errors due to current jitter are substantially reduced from previous versions of the TBA. Analytic scalings and numerical simulations are used to obtain an illustrative TBA parameter set.
A detailed treatment is introduced to measure the dynamic stability of the relativistic electrons in a self-amplified spontaneous emission free-electron laser (FEL) system, which includes the numerical approach of the Kolmogorov entropy (entropy-like quantity), the general equations of motion for a charged particle and the method of monitoring the simulation accuracy. Numerical experiments reveal a new phenomenon that there exists the possibility of the transition from chaotic to non-chaotic phase-space trajectories of the strongly relativistic electrons due to the effect of their self-fields. The adiabatic magnetic field of a one-dimensional wiggler may have a slight influence on the electron transportation in the absence of the FEL fields, but substantially affects the dynamic stability of the electrons in the process of the FEL interaction. Moreover, the laser fields diminish the dynamic stability of the electrons as the FEL interaction ...
The design and construction of a Browne-Buechner type round pole face magnetic electron spectrograph is described. The design is based on the first order transfer properties of the bending magnet including the focal surface, dispersion, magnification and resolution. Second order aberrations are compared wit numerical ray tracing results using measured field data to simulate the magnet. The numerical results are also compared with floating wire measurements on the system. Agreement is within one part in 10"4. Data extraction from the instrument is done with an array of copper detectors placed along the focal surface. The spectrograph has been successfully employed in free-electron laser (FEL) experiments for over a year.
The understanding of transport processes that occur in naturally fractured geothermal systems is far from being complete. Often, evaluation and numerical simulations of fractured geothermal reservoirs, are carried out by assuming equivalent porous media and homogeneous petrophysical properties within big matrix blocks. The purpose of this paper, is to present a comparison between results obtained from numerical studies of a naturally fractured reservoir treated as a simple porous medium and the simulation of some real aspects of the fractured reservoir. A general conclusion outlines the great practical importance of considering even approximately, the true nature of such systems. Our results show that the homogeneous simplified evaluation of the energy resource in a fractured system, could result in unrealistic estimates of the reservoir capacity to generate electricity.
The understanding of transport processes that occur in naturally fractured geothermal systems is far from being complete. Often, evaluation and numerical simulations of fractured geothermal reservoirs, are carried out by assuming equivalent porous media and homogeneous petrophysical properties within big matrix blocks. The purpose of this paper, is to present a comparison between results obtained from numerical studies of a naturally fractured reservoir treated as a simple porous medium and the simulation of some real aspects of the fractured reservoir. A general conclusion outlines the great practical importance of considering even approximately, the true nature of such systems. Our results show that the homogeneous simplified evaluation of the energy resource in a fractured system, could result in unrealistic estimates of the reservoir capacity to generate electricity.
Extraskeletal Ewing's sarcoma (ES) and primitive neuroectodermal tumor (PNET) are widely regarded as clinically and histologically identical tumors which consist of small blue round cells. Extraskeletal ESs/PNETs usually occur in the deep soft tissues of the paraspinal region, chest wall, or lower extremities. However, superficially located cases, so-called cutaneous ESs/PNETs, are exceedingly rare, and the vast majority of the reported cases present as a single small mass. We present magnetic resonance imaging (MRI) findings and clinical course of a unique case of primary cutaneous ES/PNET presenting as numerous huge masses with severe ulceration on them. (orig.)
The Fourier transform of cosmological density perturbations can be represented in terms of amplitudes and phases for each Fourier mode. We investigate the phase evolution of these modes using a mixture of analytical and numerical techniques. Using a toy model of one-dimensional perturbations evolving under the Zel'dovich approximation as an initial motivation, we develop a statistic that quantifies the information content of the distribution of phases. Using numerical simulations beginning with more realistic Gaussian random-phase initial conditions, we show that the information content of the phases grows from zero in the initial conditions, first slowly and then rapidly when structures become non-linear. This growth of phase information can be expressed in terms of an effective entropy: Gaussian initial conditions are a maximum entropy realisation of the initial power spectrum, gravitational evolution decreases the phase entropy. We show that ...
The detection and estimation of gravitational wave (GW) signals belonging to a parameterized family of waveforms requires, in general, the numerical maximization of a data-dependent function of the signal parameters. Due to noise in the data, the function to be maximized is often highly multi-modal with numerous local maxima. Searching for the global maximum then becomes computationally expensive, which in turn can limit the scientific scope of the search. Stochastic optimization is one possible approach to reducing computational costs in such applications. We report results from a first investigation of the Particle Swarm Optimization (PSO) method in this context. The method is applied to a testbed motivated by the problem of detection and estimation of a binary inspiral signal. Our results show that PSO works well in the presence of high multi-modality, making it a viable candidate method for further applications in GW data analysis.
The effect of the wind on the sustain of extreme water waves is investigated experimentally and numerically. A series of experiments conducted in the Large Air-Sea Interactions Facility (LASIF) showed that a wind blowing over a strongly nonlinear short wave group due to the linear focusing of a modulated wave train may increase the life time of the extreme wave event. The expriments suggested that the air flow separation that occurs on the leeward side of the steep crests may sustain longer the maximum of modulation of the focusing-defocusing cycle. Based on a Boundary-Integral Equation Method and a pressure distribution over the steep crests given by the Jeffreys'sheltering theory, similar numerical simulations have confirmed the experimental results
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 are expected to appear. Moreover, it is found that the frequency ...
The present paper describes the application of the commercial CFD-code FLUENT 6.1 to the numerical simulation of the flow field in a high specific speed radial fan. Two-dimensional as well as three-dimensional computations are performed using the steady 'Frozen-Rotor method' as well as the unsteady 'Sliding-Mesh method'. The predicted performance curves of the radial fan are compared with the results obtained from measurements on a test stand. Due to the relatively low computation times, the 'Frozen-Rotor method' can be used routinely for design purposes. The results of this approach can be improved if the leakage flow through the gap between the inlet nozzle and the impeller is taken into account. (orig.)
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.
Numerical results on the translocation of long biopolymers through mid-sized and wide pores are presented. The simulations are based on a novel methodology which couples molecular motion to a mesoscopic fluid solvent. Thousands of events of long polymers (up to 8000 monomers) are monitored as they pass through nanopores. Comparison between the different pore sizes shows that wide pores can host a larger number of multiple biopolymer segments, as compared to smaller pores. The simulations provide clear evidence of folding quantization in the translocation process as the biopolymers undertake multi-folded configurations, characterized by a well-defined integer number of folds. Accordingly, the translocation time is no longer represented by a single-exponent power law dependence on the length, as it is the case for single-file translocation through narrow pores. The folding quantization increases with the biopolymer length, while the rate of translocated beads at each ...
In this paper we study theoretically a pulsed corona discharge plasma in air at atmospheric pressure. A self consistent one-dimensional model, using a flux-corrected transport numerical scheme, is presented to achieve it. The charged particle kinetics in this model are described by one-dimensional continuity equations coupled with Poisson's equations. The spatio-temporal local field, charge density variations and the velocity profile of an ionizing front are calculated to describe the dynamical behavior of corona discharge plasmas. The simulation results of a wire-in-cylinder corona discharge plasma explain the physical mechanisms of discharge processes. These results may also be apply to obtain the optimizing parameters for designing the plasma reactor.
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)
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...
Designs of nuclear reactor plants aim for high performance under safety consideration. Because of large scale and high pressure/temperature conditions, data from costly mockup tests have been required to verify simulation codes of systems and components. Establishment of design by analysis (DBA) in nuclear engineering is required for development of next generation nuclear reactors. Recent powerful computers and simulation technique enable numerical analyses to predict realistic behaviors of thermo-fluid flow, structure and do on. The present report describes resent simulation results of complex gas-liquid two-phase flow, large scale structure dynamics and fluid-structure interaction. (author)
Theoretical investigations of fatty acids as a phase change material (PCM) for energy storage system have been conducted in this study. The selected fatty acids were capric acid, lauric acid, myristic acid, palmitic acid and stearic acid. For the two-dimensional simulation model based on the enthalpy approach, calculations have been made for the melt fraction with conduction only. Glass, stainless steel, tin, aluminium mixed, aluminium and copper were used as heat exchanger materials in the numerical calculations. Theoretical results show that capric acid was found good compatibility with latent heat storage system. The large value of thermal conductivity of heat exchanger materials did not make significant contribution on the melt fraction. (author)
To represent a heterogeneous unsaturated fractured rock by its homogeneous equivalent, Monte Carlo simulations are used to obtain upscaled (effective) flow properties. In this study, we present a numerical procedure for upscaling the van Genuchten parameters of unsaturated fractured rocks by conducting Monte Carlo simulations of the unsaturated flow in a domain under gravity-dominated regime. The simulation domain can be chosen as the scale of block size in the field-scale modeling. The effective conductivity is computed from the steady-state flux at the lower boundary and plotted as a function of the averaging pressure head or saturation over the domain. The scatter plot is then fitted using van Genuchten model and three parameters, i.e., the saturated conductivity K{sub s}, the air-entry parameter {alpha}, the pore-size distribution parameter n, corresponding to this model are considered as the effective K{sub s}, effective {alpha}, and effective n, respectively.
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 results. (author)
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).
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).
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
The scope of the present article is two-fold. Firstly, to conduct an experiment to provide the temperature-time history of the cooling of a hot ball bearing in quiescent ambient air. Secondly, to predict the temporal variation of the bearing under the hypothesis of natural convection, radiation or natural convection coexists with radiation for a non-vanishing total hemispherical emissivity of the surface of the bearing. Numerical solutions of the three governing nonlinear lumped heat equations were carried out with a Runge-Kutta-Fehlberg (RKF45) algorithm accounting for automatic step size control. The experimental data was obtained with chrome steel ball bearings of diameter 0.953 cm (7/16 in) heated in an electric oven to a pre-set temperature. The heated bearing was exposed later to ambient air at atmospheric temperature and pressure. (orig.)
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 ...
Mixed convection in a rectangular channel (width/height = 2) with bottom-heated and top-cooled sections is studied by laser Doppler anemometry in nitrogen at Ra = 22,200 and Re = 18.75, 36, and 54. At the lower Re values, symmetry breaking is observed in steady but spatially oscillating flows that prevail over a certain distance from the leading edge of the differentially heated section. Further downstream, unsteady flows are found even for Re = 18.75. Numerical models are used to investigate the effects of adiabatic, conducting (with a conductive-convective heat transfer coefficient), and perfectly conducting side walls; channel tilts and Prandtl number dependence. Good agreement between calculations and experiment is obtained for longitudinal convective roll velocities. The transverse velocities are found to be independent of Re.
Dark matter particles form halos that contribute the major part of the mass of galaxy clusters. The formation of these cosmological structures have been investigated both observationally and in numerical simulations, which have confirmed the existence of a universal mass profile. However, the dynamic behaviour of dark matter in halos is not as well understood. We have used observations of 16 equilibrated galaxy clusters to show that the random velocities of dark matter particles are larger on average along the radial direction than along the tangential, and that the magnitude of this velocity anisotropy is radially varying. Our measurement implies that the collective behaviour of dark matter particles is fundamentally different from that of normal particles and the radial variation of the anisotropy velocity agrees with the predictions of numerical simulation.
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.
Eddy current testing is one of the most widely used methods in non-destructive testing for the inspection of conductive materials. Numerical modelling of eddy current testing has emerged as an important approach alongside experimental studies. This paper investigates an application of numerical modelling and experimental study as a means of the quantitative non-destructive evaluation (QNDE) of defects in conductive samples. There are two methods of measuring eddy current response, more commonly by measuring the change in impedance of the eddy current probe coil, or as used in this work, by measuring the change in magnetic field directly using magnetic field sensors such as superconducting quantum interference devices, giant magneto resistance, or as in this case Hall sensors. Specifically,...
An analysis is presented for the effects of chemical reaction and thermal radiation on hydromagnetic free convection heat and mass transfer for a micropolar fluid via a porous medium bounded by a semi-infinite vertical porous plate in the presence of heat generation. The plate moves with a constant velocity in the longitudinal direction and the free stream velocity follows an exponentially small perturbation law. A uniform magnetic field acts perpendicularly to the porous surface in which absorbs the micropolar fluid with a suction velocity varying with time. Analytical expressions are computed numerically. Numerical calculations are carried out the purpose of the discussion of the results which are shown on graphs and the effects of the various dimensionless parameters entering into the p...
Determination of the maximum temperature and temperature distribution along the rake face of the cutting tool is of particular importance because of its controlling influence on tool life, as well as, the quality of the machined part. Numerous attempts have been made to approach the problem with different methods including experimental, analytical and numerical analysis. Although considerable research effort has been made on the thermal problem in metal cutting, there is hardly a consensus on the basics principles. The unique tribological contact phenomenon, which occur in metal cutting is highly localized and non-linear, and occurs at high temperatures, high pressures and high strains. This has made it extremely difficult to predict in a precise manner or even assess the performance of va...
The abstract is additional with repect to the paper published in JGCD. Ordinary Earth satellites are usually placed into five categories of special orbits: sun-synchronous orbits, orbits at the critical inclination, frozen orbits, repeating ground track orbits, and geostationary orbits. This paper investigates their counterparts around Mars and examines the basic nature of these orbits, which are of special interest for missions conducted around Mars, including Mars reconnaissance. Mars' gravity field is much more complicated, with relatively smaller J2, compared to Earth's, which makes the behaviors of these Martian orbits different from those of Earth. Analytical formulations and numerical simulations are used to analyze these Martian orbits and compare them with their Earth counterparts. First, mean element theory is employed to describe variations of orbital elements and give the constraint conditions for achieving these special orbits. Then, ...
In this research, the interactive effect of grain and specimen sizes on the flow stress of sheet metal in microforming is investigated via the tensile test of pure copper and numerical modeling. Models based on different assumptions are proposed to analyze the size effect phenomenon. It is found that the flow stress decreases linearly with the decrease of the ratio of specimen to grain sizes. The grain boundary thickness decreases and its volume fraction increases with the decrease of grain size. The variation of grain boundary thickness is not proportional to the variation of grain size. Furthermore, the fraction of grain boundary increases with the strain and the ratio of specimen to grain sizes. Based on the FE simulation, it is found that the simulated flow stress, which is modeled bas...
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.
The present invention concerns a method of evaluating corrosion resistance of a zirconium alloy easily and in a short period of time. An anode polarization curve of the zirconium alloy is measured to obtain an anode polarization curve expressed by logarithm. The curve is converted to a potential-current density curve expressed by absolute values. The peak area in the curve of the converted potential-current density is indicated by numerical values. Further, the corrosion resistance of the zirconium alloy is evaluated based on the peak area converted into the numerical values as a reference. This method is based on the finding that the peak area has a close relation with nodular corrosion resistance, and the corrosion resistance can be judged with respect to a specific zircaloy-2. (T.M.).
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)
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.
The code is being developed starting from the steady-state finite element code FENDER for the solution of the diffusion equation by extending it to become time-dependent. The numerical solution of the time-dependent multigroup diffusion equations within TRANSFUSION is performed at the present stage of development by using a backward difference scheme for the time variable, leading to a rearrangement of FENDER by adding a new loop over time steps. The code retains the multigroup coupled neutron-gamma features of FENDER, and provides a consistent two-, and quasi three-dimensional numerical solution of both static and time-dependent multigroup diffusion equations. (orig./DG)
Design of pipe whip restraints requires a knowledge of the maximum allowable unsupported pipe length. This paper presents a numerical method for calculating this critical length of the pipe. Salient features of the method are: (1) as a flow rounds an elbow, it exerts a transverse kick and an axial thrust to the pipe, both the axial thrust and the bending moment are considered; (2) the jet force is applied in an abrupt manner, the dynamic amplification factor (DAF) is determined from the load-deflection (H-{Delta}) curve of the pipe, by taking into consideration large strain, large deformation and the nonlinear stress-strain relationship of the piping material; (3) the ultimate capacity of the pipe under the combined action of an axial force and a bending moment is governed by an interaction formula. The maximum unsupported pipe length is calculated from this formula. Only the formulas and derivation are presented. Numerical results and design ...
A review of the deep-level defects observed in both electron- and proton-irradiated GaAs solar cells is presented. Studies of the effects of periodic and continuous thermal annealing on the radiation-induced electron and hole traps and the recombination parameters in GaAs solar cells were made for a wide range of electron and proton energies, fluence, annealing temperature and annealing time. A refined model for numerical simulations of the displacement damage was developed for computing the defect density and the cell parameters in the electron- and proton-irradiated GaAs solar cells. Excellent agreement was obtained between the calculated values and the experimental data for the proton-irradiated GaAs solar cells. (orig.).
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...
This paper describes the study of particles' dispersion in an isotropic turbulent flow. The particle's motion and the turbulent flow characteristics are calculated independently. While the particles' displacement is computed by the author's code, the flow is simulated with a commercial code : PowerFLOW. The particles and the flow are coupled through the relative velocity component of the aerodynamic force. When the simulated flow is turbulent, a turbulence regeneration model is used in order to get the flow instantaneous velocity. Validation of the method is done by comparing the particles' dispersion obtained with experimental results from literature and with the results calculated by FLUENT. Good accordance is found between numerical studies and experimental results. However, comparison between results of PowerFLOW coupled to the author's code and results from FLUENT shows differences when the particle's path goes through the turbulent boundary layer. (author)
A procedure for code Verification by the Method of Manufactured Solutions (MMS) is presented. Although the procedure requires a certain amount of creativity and skill, we show that MMS can be applied to a variety of engineering codes which numerically solve partial differential equations. This is illustrated by detailed examples from computational fluid dynamics. The strength of the MMS procedure is that it can identify any coding mistake that affects the order-of-accuracy of the numerical method. A set of examples which use a blind-test protocol demonstrates the kinds of coding mistakes that can (and cannot) be exposed via the MMS code Verification procedure. The principle advantage of the MMS procedure over traditional methods of code Verification is that code capabilities are tested in full generality. The procedure thus results in a high degree of confidence that all coding mistakes which prevent the equations from being solved correctly ...
Electromagnetic forming (EMF) is a forming process that relies on the inductive electromagnetic force to deform metallic workpiece at high speed. Calculation of the electromagnetic force is essential to understand the EMF process. However, accurate calculation requires complex numerical solution, in which the coupling between the electromagnetic process and the deformation of workpiece needs be considered. In this paper, an appropriate formula has been developed to calculate the electromagnetic force in metal work-piece in the sheet EMF process. The effects of the geometric size of coil, the material properties, and the parameters of discharge circuit on electromagnetic force are taken into consideration. Through the formula, the electromagnetic force at different time and in different positions of the workpiece can be predicted. The calculated electromagnetic force and magnetic field are in good agreement with the numerical and experimental ...
We present a theoretical study of the resonant interaction between dynamical localized states (discrete breathers) and linear electromagnetic excitations (EE's) in Josephson junction ladders. By making use of direct numerical simulations we find that such an interaction manifests itself by resonant steps and various sharp switchings (voltage jumps) in the current-voltage characteristics. Moreover, the power of ac oscillations away from the breather center (the breather tail) displays singularities as the externally applied dc bias decreases. All these features may be mapped to the spectrum of EE's that has been derived analytically and numerically. Using an improved analysis of the breather tail, a spectroscopy of the EE's is developed, The nature of breather instability driven by localized EE's is established.
A numerical analysis is performed to clarify the heat transfer characteristics of combined convection around the tubes which form a single row perpendicularly arranged to a forced flow in the direction of gravity. A boundary-fitted coordinate transformation technique is adopted to solve the governing equations numerically. It is found in the case of liquid sodium that gravity-induced flow gives no substantial effect to heat transfer augmentation in a low Reynolds number region while it can be sufficiently expected in the ordinary fluids in the same Reynolds number region. This means that such heat transfer augmentation as expected in ordinary fluids can not be realized when crossflow-type heat exchangers installed for decay-heat removal operates in a low Reynolds number region.
In order to reduce the size and cost of heat exchangers, an air-side wavy fin-and-tube heat transfer surface with three-row tubes needs to be replaced by two-row tubes with some appropriate enhancing techniques. The major purpose of the present paper is to search for such new structure by numerical simulation. First, longitudinal vortex generators of Delta-winglet type are tried. The influence of number and of arrangement of the winglets on the performance of the heat transfer surface is studied in detail. The numerical results show that the fin with two winglets aligned spanwise in the front and rear of each tube (Fin W6) has higher heat transfer capability than other enhanced structures with vortex generators, but it still unable to meet the heat transfer requirement. Then a combination ...
The types of investigations undertaken to test for possible biological effects of extremely low frequency electric fields have been numerous. However, neither animal and plant experimentation nor clinical studies nor experience with operating extremely high voltage transmission lines have to date provided convincing evidence of a harmful effect from exposure to electric fields associated with transmission lines in spite of numerous attempts to find such effects. Analysis of internal fields and currents supports these observations as the levels appear to be too low to affect mammalian cells. Thus, while one can never prove the negative (i.e. that there is no effect), the overwhelming body of evidence indicates that the electric fields associated with high-voltage lines have no deleterious biological effects. (author).
The bulk removal of CO/sub 2/ from industrial gases is a frequently applied process in the field of gas-treating. Usually chemical solvents are used in order to increase both the solvent-loading capacity and the absorption rate. These solvents often are aqueous solutions of alkanolamines. The loaded-liquid is regenerated at high temperatures with the aid steam. In the present study absorption rates are experimentally determined in various mixtures of amines over a wide range of concentrations. Simultaneously, a model was developed for the description of mass transfer accompanied by parallel reversible chemical reactions. The equations of this mass transfer model could not be solved analytically and therefore numerical techniques were used to obtain an exact description of this process. This model also includes the asymptotic situations of homogeneous catalysts and shuttle mechanism for which in the literature several approximate solutions have been presented. The ...
Biostimulation by nutrient application is a viable technology for restoring oil-contaminated beaches. Maximizing the nutrient residence time is key for achieving a rapid cost-effective cleanup. We considered the nutrient injection strategy through a perforated pipe at the high tide line and we simulated numerically beach hydraulics, which allowed us to estimate the optimal injection flow rate of nutrient solution. Our results indicate that the optimal application is one that starts following the falling high tide and lasts for half tidal cycle. The saturated wet-front of the nutrient solution on the beach surface would move seaward with the same speed of the falling tide keeping a constant distance with the tide line. The numerical results were generalized to beaches of wide ranges of hydr...
The FVM-LES-acoustic analogy method (FVM-LES-AAM), which is a hybrid prediction technique for the acoustical property computation, is presented and performed in this paper. The FVM-LES-AAM was developed by combining the finite volume method (FVM), the large eddy simulation (LES), and the Ffowcs Williams-Hawkings analogy algorithm (FWH-AA). To predict the acoustical properties of induction cookers, the FVM is used for discretizing the calculation field and building numerical equations, and the LES and FWH-AA are performed for computing the sound sources and predicting the far-field sound, respectively. Using the FVM with the unstructured grids method to discretize the control equation of Navier-Stokes was introduced for illuminating the above numerical simulation procedure. To prove the FVM...
In this paper, the unsteady MHD free convection heat and mass transfer of viscous fluid flowing through a Darcian porous regime adjacent to a moving vertical semi-infinite plate under Soret and Dufour effect have been examined. Viscous dissipation effects are included in the energy equation. A uniform magnetic field is applied transversely to the direction of the flow. The differential equations governing the problem have been transformed by a similarity transformation into a system of non-dimensional differential equations which are solved numerically by element free Galerkin method. The influence of Grashof number (Gr), magnetic parameter (M), heat absorption parameter (Q), permeability parameter (K), Schmidt number (Sc), Soret number (Sr), and Dufour number (Du) on the velocity, tempera...
A numerical analysis is performed to clarify the heat transfer characteristics of combined convection around the tubes which form a single row perpendicularly arranged to a forced flow in the direction of gravity. A boundary-fitted coordinate transformation technique is adopted to solve the governing equations numerically. It is found in the case of liquid sodium that gravity-induced flow gives no substantial effect to heat transfer augmentation in a low Reynolds number region while it can be sufficiently expected in the ordinary fluids in the same Reynolds number region. This means that such heat transfer augmentation as expected in ordinary fluids can not be realized when crossflow-type heat exchangers installed for decay-heat removal operates in a low Reynolds number region. (author).
In this paper, we have designed power system stabilizer (PSS) using object-oriented method. There are several types of power system stabilizer, A proportional-integral (PI) controller is very simple for practical implementation. Therefore it has been widely employed by the industry. The methods of obtaining the gains (Ki,Kp) of PI controller are root-locus method and sub-optimal regulator approach. But these methods are cannot be applicable to nonlinear system and faulted power system. So we proposed a new method which can be applied to nonlinear system by numerical analysis method. The method of dynamic system simulation by numerical method is very difficult and complex. We proposed flexible simulation method for complex power system analysis using object-oriented programming (OOP) and applied to PI controller design. (author). 7 refs., 14 figs., 2 tabs.
A new semi-submersible floating structure is proposed on which three wind turbine towers are installed. This paper presents a basic characteristic of the wave-induced motion of this semi-submersible floating structure via. numerical computations and 1/150 scaled rigid model experiments in a wave tank. In the numerical computations, nonlinear damping effect due to drag forces modeled by the Morison's formula is considered in the equation of motion, where the linear hydrodynamic forces are obtained from the Green's function model. As a result, the response characteristics around the resonant frequency region were successfully improved. In addition to such basic examination, major results of feasibility studies, including the structural stability for severe wave conditions and the long-term fatigue limit state, are presented for a realistic situation.
The WBURN (2-D, 2-group, coarse mesh) code is developed to analyze the equilibrium core characteristics of CANDU-PHWR. The equilibrium characteristics of Wolsung reactor computed by using WBURN are compared with the values given in the Wolsung FSR. The changes of equilibrium core characteristics caused by the variation of design parameters for operating conditions are also investigated. The numerical results indicate that the average discharge irradiation in the Wolsung reactor can be increased up to about 5%.
The WBURN (2-D, 2-group, coarse mesh) code is developed to analyze the equilibrium core characteristics of CANDU-PHWR. The equilibrium characteristics of Wolsung reactor computed by using WBURN are compared with the values given in the Wolsung FSR. The changes of equilibrium core characteristics caused by the variation of design parameters for operating conditions are also investigated. The numerical results indicate that the average discharge irradiation in the Wolsung reactor can be increased up to about 5%. (Author).
This paper describes a dissipativity-based framework for the study of low-frequency oscillations in power systems and for power system stabilizer design. This framework leads to a robust controller design formulation, amenable to both H{sub {infinity}} and QFT tools. An illustrating numerical example presents QFT based design for a widely used benchmark two area, four machine power system.
A tissue culture process, organogenesis, is described which involves the stepwise production of a shoot apex, which grows to have leaves and stem and finally roots, culminating in the production of a complete plant. The production and subsequent development of numerous adventitious shoots, preferably from adult tree tissues offers the greatest potential for commercial cloned plantlet production. The process describes these steps for Pinus taeda, the loblolly pine. (Refs. 25).
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).
The mechanism proposed by Kormendy (1984) for the formation of counterrotating cores in elliptical galaxies is investigated using self-consistent numerical simulations of mergers between a high- and a low-luminosity elliptical galaxies. The conditions for a counterrotation to appear are determined, observational properties of the remnants are described, and the evolution of the structural and kinematic parameters of the larger galaxy is analyzed. It is shown that a counterrotation results only when the merging orbits are retrograde, due to a large change in the secondary spin during the merger. 36 refs.
We demonstrate a two-dimensional device simulator for MOSFET structures that incorporates models for defect distributions and show predicted effects on device switching performance for various spatial distributions of defects in amorphous and polycrystalline silicon.
Ultrasound tomography is an attractive imaging method for the detection of breast cancer. The complex anatomy of the breast with its different spatial scales and material property contrasts make accurate reconstructions very challenging. This paper proposes a hybrid approach whereby Travel-of-Flight and Diffraction Tomography are combined together to achieve high-resolution and high-accuracy sound-speed reconstructions. The method is validated with several numerical phantoms.
The temperature coefficient has been investigated on the Wolsung nuclear power reactor, in which fuel is natural uranium dioxide and moderator heavy water. The numerical computations are carried out in terms of changes of the effective neutron multiplication factor with respect to fuel, moderator, and coolant temperatures. Those results are compared with the computed values of temperature coefficient based on the LATREP computer code. (author).
Surface deposition probe techniques have been shown to be an effective means of determining the characteristics of the plasma edge region of magnetic confinement devices. Experimental results on ISX-B and EBT-S have led to a consistent picture of the edge plasma and the dominant impurity introduction mechanisms. This picture is supported by numerous other measurements and emphasizes the importance of concentrating several techniques on the complex problems of plasma fusion.
This paper is concerned with the application of a new fuzzy inference method to the design of fuzzy power system stabilizer. Based on Truth-Value-Flow-Inference (TVFI), a fuzzy logic controlled power system stabilizer is proposed to enhance the stability of power system with weak or negative damping. The design of the proposed fuzzy power system stabilizer is presented. Numerical simulation results are presented and compared with two other kinds of fuzzy power system stabilizers. (author)
The structure of the P-2, S-1 accelerator mode with one step, S-1, on double period, P-2 for systems, described by the Harper mapping as well as by standard mapping, is analyzed. Detailed analysis of stability criteria is given and comparison with the results of numerical analysis is performed. Three-period compression of the P-2, S-1 island is revealed in the standard mapping. Refs. 7, figs. 7.
A numerically feasible method, based on the use of deformed phonons, is developed for the diagonalization of the collective quadrupole Hamiltonian for a system with an odd particle coupled to an anharmonic even core. Examples: the transition from prolate to oblate via #gamma#-unstable shapes and furthermore the hsub(9/2) spectra of the nuclei "1"8"7Ir and "1"9"7Tl. (Auth.).
In this work, the thermal performance of a conventional collector is improved by inserting porous substrates at the inner walls of the collector tubes. The porous substrates improve the convective heat transfer coefficient between the tube wall and the fluid. This improvement is investigated numerically and its effects on the efficiency and the useful gain of the collector are evaluated. It is found that inserting the porous substrate may raise the collector efficiency considerably, especially at high values of the overall heat loss coefficient.
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.
The peculiarities of pulsed laser melting and evaporation of the superconducting ceramics are analyzed by means of numerical simulation. The appearance of the overheated metastable states in solid and liquid phases is shown as a result of the phase front dynamics and volume nature of laser energy release. A method of dynamic adaptation for the multifront Stefan problem is proposed.
The model of transient enhanced diffusion of ion-implanted As is formulated and the finite-difference method for numerical solution of the system of equations obtained is developed. The nonuniform distribution of point defects near the interface and more accurate description of arsenic clustering are simultaneously taken into account. Simulation of As diffusion during rapid annealing gives a reasonable agreement with the experimental data. (authors)
The article explains the method of selecting the optimum energy for correcting the energy selectivity of tomographic projection and the optimum length of the crystals corresponding to it when the dual energy method in computerized tomography is used. Results of numerical calculations for actual materials and count results are presented.
Sandia National Laboratories maintains several libraries of equation of state tables, in a modified Sesame format, for use in hydrocode calculations and other applications. This report discusses one of those libraries, the seslan file, which contains 78 tables from the Los Alamos equation of state library. Minor changes have been made to these tables, making them more convenient for code users and reducing numerical difficulties that occasionally arise in hydrocode calculations.
Vibration of machinery plays an important role in preventive maintenance programs of machinery. It can be applied in different phases ranging from simple to computerized programs. The advantages of applying machine condition monitoring programs are numerous since vibration measurement and analysis leads to the early detection of component deterioration. The cost, savings and some details of the vibration preventive maintenance program are discussed.
The relativistic motion of an electron is numerically analyzed in a tokamak having macroscopic magnetic turbulence. Stochasticity induced by the relativistic motion overwhelms the phase averaging effect, which provides a tokamak with an effective loss mechanism for the avoidance/suppression of runaway electron generation at a major disruption. On the other hand, electrons in the KAM (Kolmogorov-Arnold-Moser) region will be observed as a runaway snake. (author)
The proton therapy gives an irradiation adapted to the shape and the depth of tumors ( eyes, head and neck) with a millimetric, even sub millimetric precision for the ophthalmological tumors. It is equally an excellent tool for the pediatrics tumors. For numerous treatments in France the proton therapy is associated to the photons treatments (x radiation). however, others physical treatments could modulate the medical offer, such thermal therapy or treatments by ultrasonic waves. (N.C.)
An introduction is presented on the radiation hazards connected with biomedical radiography and nuclear medicine. The frequency of radiodiagnostic efforts was rather high in the Netherlands. This was reduced considerably by abolishing the thorax screening of the population. About diagnostic nuclear medicine less can be said because far fewer numerical data are available. An exposition of genetically and somatically significant doses and how to compute them is given. The drawing up of a profit versus risk evaluation for medical applications of ionizing radiations is recommended. (Auth.).
Pumped storage plants (PSP) supply up to several GW of power and several ten GWh of capacity. Short start-up times and low start-up costs predestine PSP for the control energy market. Grid charges for PSP endanger the further integration of renewable energy and lead to higher electricity supply costs. Therefore, an amendment of legislation is urgently required. Requirements for building large PSP are met at numerous sites in the low mountain ranges of Germany. (orig.)
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.
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.
Plant roots release about 5% to 20% of all photosynthetically-fixed carbon, and as a result create a carbon-rich environment for numerous rhizosphere organisms, including plant pathogens and symbiotic...Full Text Available
After reviewing some of the mathematical foundations and numerical difficulties facing lattice QCD, I review the status of several calculations relevant to experimental high-energy physics. The topics considered are moments of structure functions, which may prove relevant to search for new phenomena at the LHC, and several aspects of flavor physics, which are relevant to understanding CP and flavor violation.
This paper suggests how nonlinear periodic optimal control of a pumped storage plant can be realized. The control problem consists in maximization of the plant benefits over an operational period. An optimal control law is proposed, yielding a bang-off-bang mode of operation. An algorithm for numerical solution of the problem was developed, and its effectiveness was demonstrated by simulation experiments.
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 non-equilibrium wall functions for near-wall turbulence treatment ...
This paper presents the results of numerical investigations on electron trajectories in different 3D fringing fields of spherical defectors. We start with a chosen set of geometrical values for the electrodes such as radius of curvature, electrode gap, angle between entrance and exit face. We do a first numerical calculation of the spatial potential distribution using a 3D electron optics program based on the finite difference method (FDM). Ray tracing through this potential distribution is then performed using another program, which also gives the potential and the electric field strength along the curved beam path. In the next step we estimate the position and the geometry of the field clamp using Herzog's formulae assuming the entrance of a parallel plate condenser to be a good approximation for the entrance of the spherical deflector in first order. For the consecutive numerical recalculations of the real potential ...
The basic objective of this project was to consider a large class of matrix computations with particular emphasis on algorithms that can be implemented on arrays of processors. In particular, methods useful for sparse matrix computations were investigated. These computations arise in a variety of applications such as the solution of partial differential equations by multigrid methods and in the fitting of geodetic data. Some of the methods developed have already found their use on some of the newly developed architectures.
To determine the transient response of linear and time-invariant transmission systems which are only described for discret frequencies, it is possible to find a good functional approximation by means of the generalized impulse method, combined with the most important natural spline-interpolation functions of the first and the third degree. This procedure has been applied to calculate the step response of the transient behaviour of the zero-sequence impedance system of the high voltage transmission line.
The study of the melting of combustible rods is of great importance for the safety analysis of nuclear reactors. Due to the special characteristics of the problem, a sharp interface between the solid and liquid region does not exist, but appears a 'mushy' region in which the material is partially melted. The Finite Element Method is employed here, together with a regularized enthalpy formulation. Finally, the results obtained are presented and discussed. (Author).
This paper describes the use of numerical analysis for studying natural convection cooling systems for long term storage of heat producing radioactive materials, including special nuclear materials and nuclear waste. The paper explains the major design philosophy, and shares the experiences of numerical modeling. The strategy of storing radioactive material is to immobilize nuclear high-level waste by a vitrification process, convertion it into borosilicate glass, and cast the glass into stainless steel canisters. These canisters are seal welded, decontaminated, inspected, and temporarily stored in an underground vault until they can be sent to a geologic repository for permanent storage. These canisters generate heat by nuclear decay of radioactive isotopes. The function of the storage facility ventilation system is to ensure that the glass centerline temperature does not exceed the glass transition temperature during storage and the vault ...
The qualitative technique of the analysis of the efficiency of the RF plasma production is presented in which the solution of boundary problem for Maxwell's equations is only necessary. The analysis of the character of the plasma production process with the crankshaft antenna in Uragan-2M stellarator is carried out. The discussion of the calculations results is presented.
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).
The primary of this thesis is the derivation of nonisothermal equations for isotropic viscoelastic materials with the help of the thermodynamics. Because the description of all phenomena would be too ambitious for one thesis, the area of investigation has...
A mathematical model describing the behavior of a gas of any density released into a turbulent atmosphere was developed. A numerical treatment was established for two dimensional flow. An understanding of the effects of accidental or continuous release wa...
Standard methods in non-linear analysis are used to show that there exists a parabolic branching of solutions of the Lichnerowicz-York equation with an unscaled source. We also apply these methods to the extended conformal thin sandwich formulation and show that if the linearised system develops a kernel solution for sufficiently large initial data then we obtain parabolic solution curves for the conformal factor, lapse and shift identical to those found numerically by Pfeiffer and York. The implications of these results for constrained evolutions are discussed.
An inner-city development with pile foundations as well as the proximity to the river Neckar with a groundwater flow offered ideal conditions for this new building in Mannheim in order to heat and cool geothermally. Already with the planning phase, detailed numerical flow computations were necessary in order to prove the functionality of the energy concept.
A generalized method for the design of excitation control or a power system stabilizer (PSS) based on complex frequency is described. The method selects PSS parameters such that exact assignment of eigenvalues associated with the mechanical modes of oscillation to desired locations is achieved. Numerical examples are used to illustrate the concepts of the proposed new technique.
High pressure experiments in a jet-stirred reactor have been performed to study the NO{sub x} formation in lean premixed combustion of methane/air mixtures. The experimental results are compared with numerical predictions using four well known reaction mechanisms and a model which consists of a series of two perfectly stirred reactors and a plug flow reactor. (author) 2 figs., 7 refs.
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.
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.
This paper presents a numerical investigation of laminar flows and heat transfer in a horizontal rectangular channel whose top and bottom plates have been punched out in the form of a delta wing. The flow structure with respect to the generation, transport, and stability of vortices in laminar horizontal channel flows with combined forced and free convection are reported. To include free convection, Boussinesq approximation of the buoyancy is used and the flow medium is treated as incompressible.
We study the possibility of utilizing the superfluid to Mott-insulator quantum phase transition in an array of quantum well exciton-polariton traps to generate indistinguishable single photons in a massive parallel fashion. By means of analytical and numerical methods, the device operations and system properties are examined using realistic experimental parameters. Such a deterministic, massive parallel generation may find new applications in photonic quantum information processing.
This project is designed to establish one basis for empirically linking aspects of management and organizational design to the safe operation of nuclear power plants. Current work is focusing on (a) reviewing existing literature relevant to this linkage, and (b) isolating and incorporating numerous safety relevant indicators. Later stages of this project will involve hypothesis development, data collection and analysis.
It is proposed to measure the spot sizes and polarization of #gamma# beams of future #gamma#e and #gamma##gamma# colliders detecting e"+e"- pairs produced as a result of interaction of high energy #gamma# quanta with density modulated and not modulated laser photons. The quantum electrodynamics cross sections, necessary numerical results as well as a short comparison of the proposed method with some other methods are given. (orig.).
Methods of solving Navier-Stokes equations for turbulent channel flow are given. These numerical solutions utilize either Neumann or Dirichlet boundary conditions. Computer codes were written and tested and are discussed. Digital image processing of flow visualization video sequences, taken simultaneously with vorticity probe measurements in a turbulent boundary layer, were carried out. 4 figs. (GHH)
This report is a brief review of biological and medical applications of ion-track membranes. The review aims at informing nuclear physicists about alternative (i.e. non-fundamental-science) use of heavy ion accelerators such as production of micro- and nano-porous materials. The ion-track membranes produced this way are employed in life sciences and numerous technological applications. The author focuses on recent results from the Flerov laboratory in co-operation with other scientific institutions and industrial partners.
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, ...
This paper deals with the efficient simulation of the dynamical behaviour of molten carbonate fuel cells (MCFCs). MCFCs allow an efficient and environmentally friendly energy production via electrochemical reactions. Their dynamics can be described by large scale systems of up to currently 22 nonlinear partial differential algebraic equations (PDAE). The paper also serves as a basis for later parameter identification and optimal control purposes. Therefore, the numerical simulations are particularly based on hierarchically embedded systems of PDAE, first of all in one space dimension. The PDAE are of mixed parabolic-hyperbolic type and are completed by nonlinear initial and boundary conditions of mixed type. For a series of embedded models in one space dimension, the vertical method of lines (MOL) is used throughout this paper. For the semi-discretization in space appropriate difference schemes are applied depending on the type of equations. The resulting system of ...
We study the dynamics of scroll vortices in excitable reaction-diffusion systems analytically and numerically. We demonstrate that intrinsic three-dimensional instability of a straight scroll leads to the formation of helicoidal structures. This behavior originates from the competition between the scroll curvature and unstable core dynamics. We show that the obtained instability persists even beyond the meander core instability of the two-dimensional spiral wave. copyright 1998 The American Physical Society.
Heavy metals (Pb, Cd, Mn, Cu, Zn, Fe) were determined in surface waters in the surroundings of the depositories of the mining shrubs in the region of Banska Stiavnica (SR) by radionuclide X-ray fluorescence analysis. Allowed concentrations of the determined metals are exceeded numerously (multiplicity in some cases was 10"3-10"6). (author).
Heat transfer has emerged as a leading-edge technology supporting the rapid pace of developments in science and engineering ranging from energy production to electronic cooling. Subjects covered at the 6th UK National Conference on heat transfer held at Edinburgh from 15-16 September 1999 include: convection; boiling and evaporation; condensation; numerical techniques; heat exchangers; radiation and combustion; conduction; two-phase flow and heat transfer augmentation. Fifteen of the papers are abstracted here.
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.
The periodic mode is analyzed together with two conventional boundary handling modes for particle swarm. By providing an infinite space that comprises periodic copies of original search space, it avoids possible disorganizing of particle swarm that is induced by the undesired mutations at the boundary. The results on benchmark functions show that particle swarm with periodic mode is capable of improving the search performance significantly, by compared with that of conventional modes and other algorithms.
An analytic solution is obtained for the equations of resonance coherent SRS by neglecting the population of the final level of the Raman transition for the systems with the active-medium length that is smaller than the wavelength of the incident light. For the extended systems, a numerical solution is obtained. The energy distribution of the Stokes pulses is found. The large-scale (about 100%) fluctuations of the Stokes radiation energy were observed in the case of unsaturated amplified spontaneous emission. (nonlinear optical phenomena)
Nature-inspired algorithms such as Particle Swarm Optimization and Firefly Algorithm are among the most powerful algorithms for optimization. In this paper, we intend to formulate a new metaheuristic algorithm by combining Levy flights with the search strategy via the Firefly Algorithm. Numerical studies and results suggest that the proposed Levy-flight firefly algorithm is superior to existing metaheuristic algorithms. Finally implications for further research and wider applications will be discussed.
Solar cookers offer a viable alternative to conventional cooking methods in many areas, and can be an effective tool in the fight against the deforestation and desertification that plague many developing countries. However, there have been numerous obstacles to the successful dissemination of solar cookers in the past. The purpose of this paper is to identify opportunities, review constraints and develop a marketing strategy to disseminate the Sundyne Solar Cooker (SSC) in developing countries.
In this letter we report a clear and unambiguous observation of the out-of-plane quadrupole magnetic field suggested by numerical simulations in the reconnecting current sheet in the Magnetic Reconnection Experiment (MRX). Measurements show that the Hall effect is large in collisionless regime and becomes small as the collisionality increases, indicating that the Hall effect plays an important role in collisionless reconnection.
The estimation of numerical values of the mean distance of closest approach of ions, a, of heavy metal ion salts in aqueous solutions, determined from activity coefficients, as well as from different theoretical approaches, is presented and discussed.
The interaction between gravity convection and Marangoni convection in a horizontal rectangular channel filled with a liquid containing a surfactant and a drop of another liquid is numerically investigated. For large Schmidt numbers the occurring oscillatory regime of solutal convection is analyzed. In the model with a surface phase the effect of the adsorption and desorption processes on the convective flow structure is determined. The corresponding initial and boundary value problem is solved using a difference method.
Numerous studies have documented the general toxicity of the polychlorinated biphenyls (PCBs) and a sizeable literature has accumulated concerning the toxic effects of PCBs on a variety of homeostatic systems. However, current information regarding the effects of PCBs on rhythmic endocrine function is quite limited. Thus, the present study was undertaken to determine whether circadian periodicity in pituitary-adrenal function is affected by exposure to PCBs.
The mainly nonstationary operation of a solar-heated receiver can be simulated with sufficient accuracy only if data about the dynamic behavior are available. For this reason, the dynamic behavior of a solar cavity receiver with parabolic dish collector is investigated. The development of a mathematical simulation considering heat transfer and storage processes is presented and the procedure for a numerical solution is illustrated. The performance of the calculation method is finally demonstrated by simulating the passage of a cloud.
The mainly instationary operation of a solar heated receiver can be simulated with sufficient accuracy only if data about the dynamic behavior are available. For this reason, the dynamic behavior of a solar cavity receiver with parabolic dish collector is investigated. The development of a mathematical simulation considering heat transfer and storage processes is presented and the procedure for a numerical solution is illustrated. The performance of the calculation method is finally demonstrated by simulating the passage of a cloud.
A method for transforming the Euler and Navier-Stokes equations and a complete system of fluid dynamics equations in three dimensions to a closed system on any moving surface is proposed. As a result, for an arbitrary geometric configuration, the dimension of the equations is reduced by one, which makes them convenient for numerical simulation. The general principles of the method are described, and verifying examples are presented.
A distinguishing design feature of CANDU nuclear reactors is the use of horizontal fuel channels housed in a horizontal vessel called a calandria, which is made of stainless steel 304L. Each channel consists of a Zr-2.5%Nb alloy pressure tube and an externally concentric Zr-2 calandria tube. The calandria tubes are joined to the end plates (tubesheets) of the calandria vessel by joints formed by roller expansion. The bores in the tubesheets are grooved. Roller expanded joints provide a cost effective means of joining dissimilar materials, require minimal space and no maintenance. The quality of these roller expanded joints is important from a sealing, strength and stress corrosion point of view. The roller expansion process consists of expanding the calandria tubes to deform them plastically against the bores and into the grooves of the tubesheets. Therefore, understanding the effect of the number, geometry and the pitch of the grooves on the quality of a roller expanded joint is very ...
In this work the development of a third order scheme of finite differences centered in mesh is presented and it is applied in the numerical solution of those diffusion equations in multi groups in stationary state and X Y geometry. Originally this scheme was developed by Hennart and del Valle for the monoenergetic diffusion equation with a well-known source and they show that the one scheme is of third order when comparing the numerical solution with the analytical solution of a model problem using several mesh refinements and boundary conditions. The scheme by them developed it also introduces the application of numeric quadratures to evaluate the rigidity matrices and of mass that its appear when making use of the finite elements method of Galerkin. One of the used quadratures is the open quadrature of 4 points, no-standard, of Newton-Cotes to evaluate in approximate form the elements of the rigidity matrices. The other ...
Protein bodies within the endosperm of castor bean (Ricinus communis L. cv. Hale) seeds arise from numerous small vacuoles which progressively become filled with storage protein, of...Full Text Available
Extending the usual endpoint and midpoint interactions, we introduce numerous kinds of interactions, labelled by a parameter lambda and obtain a non-commutative and associative string field algebra by adding up all interactions. With this algebra we develop a covariant open bosonic string field theory, which reduces to Witten's open bosonic string field theory under a special string length choice.
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
An extremely simple and convenient method is presented for computing eigenvalues in quantum mechanics by representing position and momentum operators in matrix form. The simplicity and success of the method is illustrated by numerical results concerning eigenvalues of bound systems and resonances for Hermitian and non-Hermitian Hamiltonians as well as driven quantum systems. Various MATLAB program codes are listed. (author)
The accelerating flow of a lighter continuous phase through a heavier one is considered. Small nonuniformities grow into large ones due to the Rayleigh-Taylor instability. An experiment exemplifying the large bubble formation due to Rayleigh-Taylor instability was performed and simulated using the PHOENICS 84 computer code. The same numerical procedure was applied to the two-phase flow in a gun barrel. It shows that the acceleration provided by the movement of the projectile can cause initial nonuniformities to grow with time.
We consider the numerical discretization of the time-domain Maxwell's equations with an energy-conserving discontinuous Galerkin finite element formulation. This particular formulation allows for higher order approximations of the electric and magnetic field. Special emphasis is placed on an efficient implementation which is achieved by taking advantage of recurrence properties and the tensor-product structure of the chosen shape functions. These recurrences have been derived symbolically with computer algebra methods reminiscent of the holonomic systems approach.
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 ...
''Climates of the Oceans'' deals with the atmosphere over the world oceans and provides a treatment not only of the climatic elements such as temperature, pressure, wind, precipitation etc., but also of the circulation of the atmosphere and its changes throughout the year. The connection between sea and overlying air is examined through studies of the water and heat balance. Numerous maps and tables illustrate the elements and their fluctuations over the water and on islands and coasts.
Chlorination of water can lead to the formation of chloropicrin. The numerous potential precursors (of various reactivities) observed during this study, confirm this hypothesis. Combination of ozonation and chlorination can also lead to the formation of this compound, dangerous to health; however, the conditions of the formation and particularly the impact of a nitration reaction in the gas phase are still not clearly defined.
With the help of the symmetry at the P point, we have determined the relations between the force constants and the moduli of elasticity for A-15 compounds. Based on these relations, a semi-experimental method for the study of the lattice dynamic behaviors of A-15 compounds is developed. Numerical calculations of the phonon spectrum of Nb/sub 3/Sn and its temperature dependence have been made, and they are in good agreement with the experimental data.
The paper contains a numerical evaluation of the expressions for the absorptive potential in heavy ion reactions given earlier. With a standard folding expression for the real part of the ion-ion potential general good agreement is found with experimental data for the angular distributions of elastic and inelastic scattering. Special interest is attached to the case of /sup 16/O + /sup 28/Si where the calculated imaginary potential is very small at low bombarding energies.
The paper contains a numerical evaluation of the expressions for the absorptive potential in heavy ion reactions given earlier. With a standard folding expression for the real part of the ion-ion potential general good agreement is found with experimental data for the angular distributions of elastic and inelastic scattering. Special interest is attached to the case of "1"6O + "2"8Si where the calculated imaginary potential is very small at low bombarding energies. (orig.).
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.
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.
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.
A numerical method for the calculation of unsteady, inviscid turbomachine flow in the transonic velocity range is described. The calculation of complete stages is possible. If the flow layer geometry is given, the flow field can be determined by independe...
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.
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.).
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).
The effect of finite gyroradius on the current diffusive ballooning mode is examined. Starting from the reduced MHD equations including turbulent transports, coupling with drift motion and finite gyroradius effect of ions, we derive a ballooning mode equation with complex transport coefficients. The eigenfrequency, saturation level and thermal diffusivity are evaluated numerically from the marginal stability condition. Preliminary results of their parameter dependence is presented. (author)
This report documents the investigation regarding the failure of CPVC piping that was used to connect a solar hot water system to standard plumbing in a home. Details of the failure are described along with numerous pictures and diagrams. A potential failure mechanism is described and recommendations are outlined to prevent such a failure.
The numerous types of plastic materials and their varying structures of use call for different recycling techniques. This situation is illustrated by the survey of the most important recycling methods. Depending on origin and type of plastic waste, not only different recycling methods have to be adopted, but the resulting recycling materials are also used for different purposes. Disintegration of thermoplastic particles and waste into granular matter and grains which can be directly reused is particularly important.
In the paper is studied the physics of a free electron laser (FEL) based on a two-frequency undulator (TFU) which induces large non linear effects, especially on the spectral dynamics. These effects are analyzed in an extended formalism where the spontaneous emission, the low-gain regime and the strong-field saturation regime are studied. Numerical simulations show that the optimized TFU generates a laser field having both a large extraction efficiency and a narrow spectrum.
An interpretive theory of trace impurity transport in tokamaks is presented. It is based on an earlier phenomenological approach developed by the authors for describing tokamak plasma transport experiments. The present model is used to explain disparate phenomena observed in ALCATOR C and DITE in a unified phenomenological framework. Predictive numerical calculations are presented for DITE conditions. If experimentally confirmed, the model could be used to correlate plasma properties with trace impurity behaviour. (orig.)
A new nonlinear S{sub n} transport differencing scheme for slab geometry is presented that is fourth order accurate for small meshes and is strictly positive. The new scheme has been coded into the existing ONELD code and tested. Numerical results to demonstrate the accuracy and positivity of this new scheme are presented.
We investigate an application of the method of fundamental solutions (MFS) to the one-dimensional inverse Stefan problem for the heat equation by extending the MFS proposed in [5] for the one-dimensional direct Stefan problem. The sources are placed outside the space domain of interest and in the time interval (-T,T). Theoretical properties of the method, as well as numerical investigations, are included, showing that accurate and stable results can be obtained efficiently with small computational cost.
Ethylene (C2H4) is a unique plant-signaling molecule that regulates numerous developmental processes. The key enzyme in the two-step biosynthetic pathway of ethylene is 1-aminocyclopropane-1-carboxylate...Full Text Available
A conference on geotechnics presented papers on soil properties, soil liquefaction, slope stability, pile foundations, numerical modelling, engineering geology, rock mechanics, cold regions engineering, oil sands, irrigation, environmental engineering, geosynthetics, and geotechnical case histories. Separate abstracts have been prepared for 29 papers from the conference.
This paper is concerned with the validation of the 3D deterministic neutral-particle transport theory code EVENT for shielding applications. The code is based on the finite element-spherical harmonics (FE-P{sub N}) method which has been extensively developed over the last decade. A general multi-group, anisotropic scattering formalism enables the code to address realistic steady state and time dependent, multi-dimensional coupled neutron/gamma radiation transport problems involving high scattering and deep penetration alike. The powerful geometrical flexibility and competitive computational effort makes the code an attractive tool for shielding applications. In recognition of this, EVENT is currently in the process of being adopted by the UK nuclear industry. The theory behind EVENT is described and its numerical implementation is outlined. Numerical results obtained by the code are compared with predictions of the Monte Carlo code MCBEND and ...
The microbial transport simulator (MTS) is a three-dimensional, three-phase, multiple-component numerical model that permits the study of the transport of microorganisms and nutrients in porous media. Microbial parameters incorporated into MTS include: Microbial growth and decay, microbial deposition, chemotaxis, diffusion, convective dispersion, tumbling, and nutrient consumption. Governing equations for microbial and nutrient transport are coupled with continuity and flow equations under conditions appropriate for a black oil reservoir. The model`s mathematical formulations and preparation procedures of data files for conducting simulations using MTS are described. A general background of microbial transport simulation is given in Section I and the governing equations, mechanisms, and numerical solutions of MTS are given in Section II. Explanations for preparing an input data file with reservoir and microbial transport data are described in ...
The wavenumber-frequency spectral densities of turbulent wall pressure fluctuations are investigated over a rigid flat plate. Nonlinear Reynolds stress terms of the inhomogeneous Orr-Sommerfeld equation are regarded as a known forcing function. The forcing function is modeled after Bark{close_quote}s hydrodynamic bursting formulation. The inhomogeneous Orr-Sommerfeld equation is solved by the method of Eckhaus in terms of discrete homogeneous solutions. The method of Eckhaus is then extended and proved for the continuous Orr-Sommerfeld eigenfunctions. Turbulent wall pressure fluctuations in terms of wavenumber-frequency spectral densities are numerically computed and compared to the experimental results of Martin as well as to his transformation of Blake{close_quote}s data fitted to a modified Corcos model. The wavenumber-frequency spectral densities numerically computed from the discrete eigenfunctions compared well with Martin{close_quote}s ...
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 natural, mixed, and forced ...
In the present study, applications of the SPH method to industrial related issues are considered by starting from an existing open source 2D SPH code, namely the SPHYSICS code, which offers an effective ground for numerical developments, which are performed in order to bring an answer to industrial problems, such as simulations of solid/fluid coupling in a free surface flow context. The purpose of the present paper is therefore to expose the numerical developments which yield an enhanced version (referred to as "SPHYSIC2") of the initial code. Firstly, the different features added to obtain the operational code needed for engineering applications are described, and so are the problems raised on this way, offering a kind of review of SPH methods for engineers. Secondly, the validation of the proposed code is partially presented with two well known but difficult test cases, namely the classical "dam break" and "wedge entry" problems. Thirdly, ...
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 parallel computer nCUBE-2S is used. The ...
In the heat exchangers of power plants, scale deposition may occur, especially at the leading edge of contraction. The growth of scale can lead to an increase in pressure and cause oscillation of the water level. In our previous study, the dependence of flow on contraction geometry was analyzed numerically and empirically. It was shown that the contraction ratio of hydraulic diameter of the flow path contributed greatly to the difference in pressure drop and turbulent kinetic energy at the leading edge of the flow path. In this study, the effect of contraction geometry and flow rate on scale adhesion was studied in AVT chemistry at 270degC using a high-velocity scale-adhesion test loop. The differential pressure and the amount of scale deposited due to scale adhesion at the leading edge of the flow path increased more in a quatrefoil-type flow path than in drill-type flow path. The differential pressure and the amount of scale deposited also increased with ...
Experimental and numerical investigations were performed for the laminar burning velocity and the flame structure of laminar premixed CH{sub 4}/O{sub 2}/CO{sub 2} flames. Measurements of the laminar burning velocity were conducted by using a flame cone angle method for a circular nozzle burner. Numerical simulation was performed using one-dimensional plane flame code including radiation heat loss with an optically thin model. It was shown that the laminar burning velocity decrease with CO{sub 2} addition even though the adiabatic flame temperature is the same as that for CH{sub 4}/Air flames. The radiation heat loss is significant for the CH{sub 4}/O{sub 2}/CO{sub 2}, flames, and the flame temperature and laminar burning velocity decreases when the radiation heat loss is considered. Effects of thermal properties, radiation, and chemical reaction on the determination of the laminar burning velocity of CH{sub 4}/O{sub 2}/CO{sub 2} flames were ...
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 melt flows have been varied and the influence of the gas and ...
Observations show the ubiquitous presence of propagating magnetohydrodynamic (MHD) kink waves in the solar atmosphere. Waves and flows are often observed simultaneously. Due to plasma inhomogeneity in the perpendicular direction to the magnetic field, kink waves are spatially damped by resonant absorption. The presence of flow may affect the wave spatial damping. Here, we investigate the effect of longitudinal background flow on the propagation and spatial damping of resonant kink waves in transversely nonuniform magnetic flux tubes. We combine approximate analytical theory with numerical investigation. The analytical theory uses the thin tube (TT) and thin boundary (TB) approximations to obtain expressions for the wavelength and the damping length. Numerically, we verify the previously obtained analytical expressions by means of the full solution of the resistive MHD eigenvalue problem beyond the TT and TB approximations. We find that the ...
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 problems have been ...
A numerical investigation of a humidification dehumidification desalination (HDD) process using solar energy is presented. The HDD system consists mainly of a concentrating solar water heating collector, flat plate solar air heating collector, humidifying tower and dehumidifying exchanger. Two separate circulating loops constitute the HDD system, the first for heating the feed water and the second for heating air. A mathematical model is developed, simulating the HDD system, to study the influence of the different system configurations, weather and operating conditions on the system productivity. The model validity is examined by comparing the theoretical and experimental results of the same authors. It is found that the results of the developed mathematical model are in good agreement with the experimental results and other published works. The results show also that the productivity of the unit is strongly influenced by the air flow rate, cooling water flow rate ...
A numerical investigation of a humidification dehumidification desalination (HDD) process using solar energy is presented. The HDD system consists mainly of a concentrating solar water heating collector, flat plate solar air heating collector, humidifying tower and dehumidifying exchanger. Two separate circulating loops constitute the HDD system, the first for heating the feed water and the second for heating air. A mathematical model is developed, simulating the HDD system, to study the influence of the different system configurations, weather and operating conditions on the system productivity. The model validity is examined by comparing the theoretical and experimental results of the same authors. It is found that the results of the developed mathematical model are in good agreement with the experimental results and other published works. The results show also that the productivity of the unit is strongly influenced by the air flow rate, cooling water flow rate ...
Starting from the Molten Salt Breeder Reactor project of Oak-Ridge, we have performed parametric studies in terms of safety coefficients, reprocessing requirements and breeding capabilities. In the frame of this major re-evaluation of the molten salt reactor (MSR), we have developed a new concept called Molten Salt Fast Reactor or MSFR, based on the Thorium fuel cycle and a fast neutron spectrum. This concept has been selected for further studies by the MSR steering committee of the Generation IV International Forum in 2009. Our reactor's studies of the MSFR concept rely on numerical simulations making use of the MCNP neutron transport code coupled with a code for materials evolution which resolves the Bateman's equations giving the population of each nucleus inside each part of the reactor at each moment. Because of MSR's fundamental characteristics compared to classical solid-fuelled reactors, the classical Bateman equations have to be modified by adding two ...
The RELAP5 code has been developed for best-estimate transient simulation of light water reactor coolant systems during postulated accidents. The code models the coupled behavior of the reactor coolant system and the core for loss-of-coolant accidents and operational transients such as anticipated transient without scram, loss of offsite power, loss of feedwater, and loss of flow. A generic modeling approach is used that permits simulating a variety of thermal hydraulic systems. Control system and secondary system components are included to permit modeling of plant controls, turbines, condensers, and secondary feedwater systems. RELAP5/MOD3 code documentation is divided into seven volumes: Volume I presents modeling theory and associated numerical schemes; Volume II details instructions for code application and input data preparation; Volume III presents the results of developmental assessment cases that demonstrate and verify the models used in the code; Volume IV ...
Increase usage of magnesium components in the portable electronics industries has propelled considerable research interest on the development of magnesium forging. There have been some studies done in the area of sheet metal forming and forging of axis-symmetrical parts of magnesium alloys but there are very little work done on the plane strain backward extrusion of magnesium alloy. This paper deals with the analysis of plane strain extrusion of magnesium alloy, AZ31, using numerical modeling at different extrusion conditions of temperatures and friction coefficients. A commercially available finite element package, ANSYS/LSDYNA, was used in the finite element analysis. FE simulated results, material deformation, punch force and stress/strain distributions, were compared with actual experimental data. From the numerical analysis, a decrease in forging temperature from 300 C to 250 C will lead to a moderate increase in forming resistance. The ...
This manuscript contains the lecture notes for a course taught from July 7th through July 11th at the 1997 Numerical Analysis Summer School sponsored by C.E.A., I.N.R.I.A., and E.D.F. The subject area was chosen to support the general theme of that year`s school which is ``Multiscale Methods and Wavelets in Numerical Simulation.`` The first topic covered in these notes is a description of the problem domain. This coverage is limited to classical PDEs with a heavier emphasis on hyperbolic systems and constrained hyperbolic systems. The next topic is difference schemes. These schemes are the foundation for the adaptive methods. After the background material is covered, attention is focused on a simple patched based adaptive algorithm and its associated data structures for square grids and hyperbolic conservation laws. Embellishments include curvilinear meshes, embedded boundary and overset meshes. Next, several strategies for parallel ...
The study of natural convection flow and heat transfer within a cylindrical annulus has received considerable attention because of its numerous applications, such as in nuclear reactor design, electronic component cooling, thermal storage systems, energy conservation, energy storage, and energy transmission. Here, the effects of multiple geometric perturbations on the inner and outer cylinders of an annulus with impermeable end walls are investigated in this work. A three-dimensional study was done using a numerical scheme based on a Galerkin method of finite element formulation. The nature of the buoyancy-induced flow field has been analyzed in detail. The flow fields for the cases considered were found to be qualitatively similar, and the introduction of each additional perturbation altered the flow field in a regular and recurring manner. The introduction of each perturbation on the outer cylinder causes clockwise and counterclock-wise ...
The study of natural convection flow and heat transfer within a cylindrical annulus has received considerable attention because of its numerous applications, such as in nuclear reactor design, electronic component cooling, thermal storage systems, energy conservation, energy storage, and energy transmission. Here, the effects of multiple geometric perturbations on the inner and outer cylinders of an annulus with impermeable end walls are investigated in this work. A three-dimensional study was done using a numerical scheme based on a Galerkin method of finite element formulation. The nature of the buoyancy-induced flow field has been analyzed in detail. The flow fields for the cases considered were found to be qualitatively similar, and the introduction of each additional perturbation altered the flow field in a regular and recurring manner. The introduction of each perturbation on the outer cylinder causes clockwise and counterclock-wise ...
We present an optimal control approach for the solidification process of a melt in a container. The process is described by a two phase Stefan problem with the free boundary (interface between the two phases) modelled as a graph. We control the evolution of the free boundary using the temperature on the container wall. The control goal consists in tracking a prescribed evolution of the free boundary. We achieve this goal by minimizing a appropriate cost functional. The resulting minimization problem is solved numerically by a steepest descent method with step size control, where the gradient of the cost functional is expressed in terms of the adjoint variables. Several numerical examples are presented which illustrate the performance of the method. The novelty of the approach presented consists in using a sharp interface model for the control of the free boundary. This guarantees direct access to the free boundary as optimization variable in ...
The PLANCK SURVEYOR is a European Space Agency satellite mission to image the very faint anisotropies in the temperature of the Cosmic Microwave Background (CMB) radiation. Maynooth is actively participating in an international collaboration of scientists involved with the optical design of the High Frequency Instrument (HFI). This paper outlines research which has been undertaken in Maynooth concerned with numerical modelling of the optical characteristics of the multi-frequency array of detectors making up the HFI instrument. In the study the commercial software package ZEMAX was used to model the coupling of the focal plane HFI detectors to the PLANCK telescope. This package is particularly useful in the optical design of PLANCK because of the powerful optimisation features of the software. It is thus possible to readily determine the optimum positioning of the detectors in the focal plane of the telescope. Although the package is based on a ...
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 ...
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 vector point groups can be reproduced ...
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 vector point groups can be reproduced ...
ObjectivesThe overall scientific aim of this project is, through development of an appropriate numerical modelling tool, to identify and investigate the mechanisms that control the current observed rapid thinning and acceleration of Greenland outlet glaciers and to investigate their likely future behaviour in the context of climate warming. This aim is addressed through the four major specific objectives below. These objectives build on each and therefore the order of priority is driven by the order in wh [continued...]DescriptionCurrently, the Greenland ice sheet is undergoing rapid changes in the coastal regions which have been attributed to a general warming trend to the regions climate over the last decade. Several of the narrow and fast flowing outlet glaciers that drain the ice sheet into the sea are observed to have accelerated their flow and thinned. For instance, Jakobshavn Isbrae on the West coast of Greenland has almost doubled its speed since 2000 from ...
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 validation of our ...
Numerical studies are performed to identify an optimal range of geometric parameters in connection with the design of an impingement-conditioned small direct-injection diesel engine. The parametric studies were made in a pressure chamber with a circular protrusion on to which sprays are directed to impinge. Computational runs were made for 24 cases carefully chosen to allow examination of the effects on post-impingement spray characteristics of the following parameters: (a) the nozzle-plate distance and (b) the size and geometry pattern of the protrusion. An optimal range of these parameters has been identified through an analysis of the calculated results in terms of vapour concentration, the Sauter mean diameter of the wall spray and the dispersed wall spray volume. The results obtained also reveals that, by just changing the size of the impingement surface and the impingement distance, local fuel vaporization efficiency can be increased as high as 40 per cent ...
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 ...
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 ...
In this dissertation, a two-phase, air-droplets, dilute, turbulent, and steady state flow in a horizontal rectangular channel, is modeled and numerically simulated using a modified KIVA-3V code. The deposition of different sizes of droplets on the walls of the channel is also studied. In this model, the interaction effects between the phases (two-way coupling) are considered by source terms in the momentum and energy equations for the continuous phase and by the instantaneous local velocity of the air in the droplet equation of motion, which includes the aerodynamic and gravitational forces. The turbulence is modeled by a k-? model. The interaction effects between the turbulence and the dispersed droplets are also taken into account. The effects of the turbulence on the droplets are modeled by a fluctuating component added to the local air velocity in the droplet equation of motion. The effects of the droplets on the turbulence are modeled by two extra terms in the ...
Theoretical investigations on the phase change materials (PCMs) used as the heat storage media for box-type solar cookers have been conducted in this study. The selected PCMs are magnesium nitrate hexahydrate, stearic acid, acetamide, acetanilide and erythritol. For a two-dimensional simulation model based on the enthalpy approach, calculations have been made for the melt fraction with conduction only. Different materials such as glass, stainless steel, tin, aluminum mixed, aluminum and copper are used as the heat exchanger container materials in the numerical calculations. The large value of thermal conductivity of heat exchanger container material did not make a significant contribution on the melt fraction except for at very low thermal conductivities. Based on the theoretical results, stearic acid and acetamide are found to be good compatibility with latent heat storage system. It is also found that the initial temperature of PCM does not have very important ...
The scalar mixing time scale, a key quantity in many turbulent combustion models, is investigated for reactive scalars in premixed combustion. Direct numerical simulations (DNS) of three-dimensional, turbulent Bunsen flames with reduced methane-air chemistry have been analyzed in the thin reaction zones regime. Previous conclusions from single step chemistry DNS studies are confirmed regarding the role of dilatation and turbulence-chemistry interactions on the progress variable dissipation rate. Compared to the progress variable, the mixing rates of intermediate species is found to be several times greater. The variation of species mixing rates are explained with reference to the structure of one-dimensional premixed laminar flames. According to this analysis, mixing rates are governed by the strong gradients which are imposed by flamelet structures at high Damkoehler numbers. This suggests a modeling approach to estimate the mixing rate of individual species which ...
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, resulting in a reduction ...
We review results from general relativistic axisymmetric magnetohydrodynamic simulations of accretion in Sgr A*. We use general relativistic radiative transfer methods and to produce a broad band (from millimeter to gamma-rays) spectrum. Using a ray tracing scheme we also model images of Sgr A* and compare the size of image to the VLBI observations at 230 GHz. We perform a parameter survey and study radiative properties of the flow models for various black hole spins, ion to electron temperature ratios, and inclinations. We scale our models to reconstruct the flux and the spectral slope around 230 GHz. The combination of Monte Carlo spectral energy distribution calculations and 230 GHz image modeling constrains the parameter space of the numerical models. Our models suggest rather high black hole spin ($a_*\\approx 0.9$), electron temperatures close to the ion temperature ($T_i/T_e \\sim 3$) and high inclination angles ($i \\approx 90 \\deg$).
This is a numerical verification of two groups of pipe whip tests sponsored or cosponsored by EPRI. Experimental data of the two pipe whip tests, one by Tennessee Valley Authority (TVA) and by FRAMATOME/CEA, were provided by EPRI. A nonlinear finite element code, ABAQUS-EPGEN, developed under partial sponsorship by EPRI was used for modeling the pipe whip tests. Beam elements together with an equivalent nonlinear spring element or a partial shell mesh were used to model pipes and elbow in the pipe whip tests. Material nonlinearity due to plasticity, strain rate effects, and temperature, as well as geometric nonlinearity due to large rotation and boundary conditions were included in the study. Effects of strain rate and modeling techniques were assessed. Results by current industry approach were also included as a reference solution. This report can be used as a guideline for numerical simulation of pipe whip phenomena. 74 figs.
The availability of new, high-intensity, cold and thermal neutron sources has opened the possibility of performing high-precision fundamental neutron physics experiments, including measurements that study the hadronic weak interaction and standard model test measurements, using neutron decay. The observables in these experiments are usually correlated with the direction of neutron polarization and are often very small (10 8 10 6). Mott-Schwinger scattering of polarized neutrons can produce spin-dependent shifts in beam centroids, which has the potential to produce significant systematic effects for these types of experiments. An accurate calculation of this process for neutral atoms and basic molecules has not been carried out for low neutron energies. In this work, we derive a general expression for the electromagnetic (Mott-Schwinger) contributions to the analyzing power for low-energy neutron scattering. We obtain numerical results for 11 nuclei in the range of ...
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 method with a mesh ...
A numerical model was developed to simulate the combined effects of heat transfer, magnetite oxidation, and NO{sub x} formation in a grate-kiln furnace for taconite pellet induration. Heat transfer from the flame in the kiln was described by the net radiation method. The shrinking core model was used to account for magnetite oxidation on the grate. A novel approach to oxidation of tumbling pellets in a kiln was derived. The Zeldovich mechanism was used to predict thermal NO generation. Temperature fluctuations in the gas streams were estimated with a clipped Gaussian probability density function. The thermal energy balances and mass balances resulted in coupled systems of first-order differential equations, which were solved numerically. The model is capable of predicting NO production and magnetite oxidation in agreement with observation of plant performance. Although the design of the grate-kiln system is for efficient heat and mass transfer, ...
Sudden changes in the temperature of flowing liquids can result in transient buoyancy forces which strongly impact the flow hydrodynamics via flow stratification. These effects have been studied for the case of potential flow of stratified liquids to line sinks, but not for moderator flow in SRS reactors. Standard codes, such as TRAC and COMMIX, do not have the capability to capture the stratification effect, due to strong numerical diffusion which smears away the hot/cold fluid interface. A related problem with standard codes is the inability to track plumes injected into the liquid flow, again due to numerical diffusion. The combined effects of buoyant stratification and plume dispersion have been identified as being important in operation the Supplementary Safety System which injects neutron-poison ink into SRS reactors to provide safe shutdown in the event of safety rod failure. The MODFLOW code discussed here provides transient moderator ...
Sudden changes in the temperature of flowing liquids can result in transient buoyancy forces which strongly impact the flow hydrodynamics via flow stratification. These effects have been studied for the case of potential flow of stratified liquids to line sinks, but not for moderator flow in SRS reactors. Standard codes, such as TRAC and COMMIX, do not have the capability to capture the stratification effect, due to strong numerical diffusion which smears away the hot/cold fluid interface. A related problem with standard codes is the inability to track plumes injected into the liquid flow, again due to numerical diffusion. The combined effects of buoyant stratification and plume dispersion have been identified as being important in operation the Supplementary Safety System which injects neutron-poison ink into SRS reactors to provide safe shutdown in the event of safety rod failure. The MODFLOW code discussed here provides transient moderator ...
In this article we derive longitudinal impedance and wake function for an undulator setup with arbitrary undulator parameter, taking into account a finite transverse size of the electron bunch. Earlier studies considered a line density-distribution of electrons instead. We focus our attention on the long-wavelength asymptote (compared with resonance wavelength), at large distance of the electron bunch from the undulator entrance compared to the overtaking length, and for large vacuumchamber size compared to the typical transverse size of the field. These restrictions define a parameter region of interest for practical applications. We calculate a closed expression for impedance and wake function that may be evaluated numerically in the most general case. Such expression allows us to derive an analytical solution for a Gaussian transverse and longitudinal bunch shape. Finally, we study the feasibility of current-enhanced SASE schemes (ESASE) recently proposed for ...
In this article we derive longitudinal impedance and wake function for an undulator setup with arbitrary undulator parameter, taking into account a finite transverse size of the electron bunch. Earlier studies considered a line density-distribution of electrons instead. We focus our attention on the long-wavelength asymptote (compared with resonance wavelength), at large distance of the electron bunch from the undulator entrance compared to the overtaking length, and for large vacuum-chamber size compared to the typical transverse size of the field. These restrictions define a parameter region of interest for practical applications. We calculate a closed expression for impedance and wake function that may be evaluated numerically in the most general case. Such expression allows us to derive an analytical solution for a Gaussian transverse and longitudinal bunch shape. Finally, we study the feasibility of current-enhanced SASE schemes (ESASE) recently proposed for ...
The behavior of a small power system consisting of two interconnected generators is simulated in real-time by a prototype Digital Transient Network Analyzer (TNA). The prototype Digital TNA consists of two Computational Modules and one I/O Module. The Modules communicate with each other through ribbon cables. Each Computational Module simulates one Turbo-Generator, its Transformer, its Governor, Exciter, and Power System Stabilizer Systems. The numerical integration is shared by two TMS320C30 DSPs at a step-size of 100 microseconds in real-time. The I/O module post-processes the state variables and presents selected information for analog display. The paper presents oscillograms from a test program which includes symmetry checks and behavioral checks against well known waveforms of hunting oscillations, synchronization out-of-phase torques, and subsynchronous resonance phenomena. The success of the Digital TNA depends on: (a) the theoretical method of Decoupled ...
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)
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)
The influence of the velocity profile at the inlet boundary on the simulation of air velocity distribution inside an electrostatic precipitator is presented in this study. Measurements and simulations were performed in a duct and an electrostatic precipitator (ESP). A four-hole cobra probe was used for the measurement of velocity distribution. The flow simulation was performed by using the computational fluid dynamics (CFD) code FLUENT. Numerical calculations for the air flow were carried out by solving the Reynolds-averaged Navier-Stokes equations coupled with the realizable k-{epsilon} turbulence model equations. Simulations were performed with two different velocity profiles at the inlet boundary - one with a uniform (ideal) velocity profile and the other with a non-uniform (real) velocity profile to demonstrate the effect of velocity inlet boundary condition on the flow simulation results inside an ESP. The real velocity profile was obtained from the velocity ...
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 ...
We propose a method to determine the cosmic mass density Omega from redshift-space distortions induced by large-scale flows in the presence of nonlinear clustering. Nonlinear structures in redshift space such as fingers of God can contaminate distortions from linear flows on scales as large as several times the small-scale pairwise velocity dispersion sigma_v. Following Peacock & Dodds (1994), we work in the Fourier domain and propose a model to describe the anisotropy in the redshift-space power spectrum; tests with high-resolution numerical data demonstrate that the model is robust for both mass and biased galaxy halos on translinear scales and above. On the basis of this model, we propose an estimator of the linear growth parameter beta = Omega^0.6/b, where b measures bias, derived from sampling functions which are tuned to eliminate distortions from nonlinear clustering. The measure is tested on the numerical data and found to recover ...
Through DNA methylation, histone modifications, and small regulatory RNAs the epigenome systematically controls gene expression during development, both in utero and throughout life. The epigenome is also a very reactive system; its labile nature allows it to sense and respond to environmental perturbations to ensure survival during fetal growth. This pliability can lead to aberrant epigenetic modifications that persist into later life and induce numerous disease states. Endocrine-disrupting compounds (EDCs) are ubiquitous chemicals that interfere with growth and development. Several EDCs also interfere with epigenetic programming. The investigation of the epigenotoxic effects of bisphenol A (BPA), an EDC used in the production of plastics and resins, has further raised concern over the impact of EDCs on the epigenome. Using the Agouti viable yellow (A(vy)) mouse model, dietary BPA exposure was shown to hypomethylate both the A(vy) and the Cabp(IAP) metastable ...
As a three-year joint university-industry effort, development of a generalized boiling transition analysis method has been started in 2002 aiming at enhanced capabilities of subchannel analysis for a wide variety of BWR-type fuel bundle geometry from ordinary BWR to tight lattice fuel bundles. For this purpose, five dominant factors affecting boiling transition phenomena have been identified on which our efforts of experimentation and numerical analyses are focused. In this report, as the first-year achievement, we will describe a master plan of the development and contents for experimental approaches to construct thermal-hydraulic databases. The databases will be utilized for the developments of constitutive equations to describe the basic characteristics of the elementary processes. The planned experiments are divided into two groups. One is air-water experiments at atmospheric pressure, and the other is steam-water experiments up to 1 MPa. The former group of ...
As a three-year joint university-industry effort, development of a generalized boiling transition analysis method has been started in 2002 aiming at enhanced capabilities of subchannel analysis for a wide variety of BWR-type fuel bundle geometry from ordinary BWR to tight lattice fuel bundles. For this purpose, five dominant factors affecting boiling transition phenomena have been identified on which our efforts of experimentation and numerical analyses are focused. In this report, as the first-year achievement, we will describe a master plan of the development and contents for experimental approaches to construct thermal-hydraulic databases. The databases will be utilized for the developments of constitutive equations to describe the basic characteristics of the elementary processes. The planned experiments are divided into two groups. One is air-water experiments at atmospheric pressure, and the other is steam-water experiments up to 1 MPa. The former group of ...
We present ultraviolet spectra of Mira's companion star from the Space Telescope Imaging Spectrograph (STIS) instrument on board the Hubble Space Telescope (HST). The companion is generally assumed to be a white dwarf surrounded by an accretion disk fed by Mira's wind, which dominates the UV emission from the system. The STIS UV spectrum is dominated by numerous, narrow H2 lines fluoresced by H I Ly-alpha, which were not detected in any of the numerous observations of Mira B by the International Ultraviolet Explorer (IUE). The high temperature lines detected by IUE (e.g., C IV 1550) still exist in the STIS spectrum but with dramatically lower fluxes. The continuum fluxes in the STIS spectra are also much lower, being more than an order of magnitude lower than ever observed by IUE, and also an order of magnitude lower than fluxes observed in more recent HST Faint Object Camera objective prism spectra from 1995. Thus, the accretion rate onto Mira ...
DF (decontamination factor) of two types of air filter was investigated against /sup 198/Au aerosol particles with AMAD (activity median aerodynamic diameter) = 1.2 ..mu..m and sigmag (geometric standard deviation) = 2.0. One of the filters, SO filter, is HEPA (high efficiency particulate air) filter media with 46 mmphi, and the other is HE-40T (cellulose-glass fiber filter for dust sampling) filter media with 40 mmphi. DF of the SO filter was found to be over 5 x 10/sup 4/ at air flow rate of 1 lmin. And DF of HE-40 T filter was determined to be 94+-4 at the same airflow rate. In the multistage HE-40 T filter system, DF decreased at each succeeding stage. It was confirmed by numerical simulation that the decrease of DF was not inconsistent with the decrease of collection performance in each successive filter. It was also confirmed that DF could be predicted by the numerical simulaiton, if particle penetration function of filter and size ...
DF (decontamination factor) of two types of air filter was investigated against "1"9"8Au aerosol particles with AMAD (activity median aerodynamic diameter) = 1.2 #mu#m and #sigma#g (geometric standard deviation) = 2.0. One of the filters, SO filter, is HEPA (high efficiency particulate air) filter media with 46 mm#phi#, and the other is HE-40T (cellulose-glass fiber filter for dust sampling) filter media with 40 mm#phi#. DF of the SO filter was found to be over 5 x 10"4 at air flow rate of 1 l/min. And DF of HE-40 T filter was determined to be 94#+-#4 at the same airflow rate. In the multistage HE-40 T filter system, DF decreased at each succeeding stage. It was confirmed by numerical simulation that the decrease of DF was not inconsistent with the decrease of collection performance in each successive filter. It was also confirmed that DF could be predicted by the numerical simulaiton, if particle penetration function of filter and size ...
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 without invoking {\\it ad hoc} assumptions on ...
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 Brasil. (author)
The convection heat loss from cavity receiver in parabolic dish solar thermal power system can significantly reduce the efficiency and consequently the cost effectiveness of the system. It is important to assess this heat loss and subsequently improve the thermal performance of the receiver. This paper aims to present a comprehensive review and systematic summarization of the state of the art in the research and progress in this area. The efforts include the convection heat loss mechanism, experimental and numerical investigations on the cavity receivers with varied shapes that have been considered up to date, and the Nusselt number correlations developed for convection heat loss prediction as well as the wind effect. One of the most important features of this paper is that it has covered numerous cavity literatures encountered in various other engineering systems, such as those in electronic cooling devices and buildings. The studies related ...
The SONACO experiments are conducted on an electrically heated 37-pin rod bundle, immersed in liquid sodium and contained within a hexagonal wrapper. The rig was designed to investigate natural convection cooling for a geometry representative of fast reactor fuel assemblies. Heat can be removed from the test section in several ways, but in this paper only the axial cooling mode is examined. Above the heated bundle is a plenum, at the top of which is a cooling coil containing a separate, forced sodium flow. Heat transfer from the bundle to this cooling coil is effected by means of buoyancy driven circulatory flow in the sodium, and in the axial cooling mode almost all the heat is removed by the coil. This mode is intended to simulate the natural convection cooling of a blocked fuel assembly by way of thermosyphon coupling to the inner pool. In this paper experimental results are presented, for the temperatures measured under such conditions, and these are compared with results from a ...
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 Prefecture, it was found that the result may ...
During the course of a fiscal year, Oak Ridge National Laboratory`s Carbon Dioxide Information Analysis Center (CDIAC) distributes thousands of specialty publications-numeric data packages (NDPs), computer model packages (CMPs), technical reports, public communication publications, newsletters, article reprints, and reference books-in response to requests for information related to global environmental issues, primarily those pertaining to climate change. CDIACs staff also provides technical responses to specific inquiries related to carbon dioxide (CO{sub 2}), other trace gases, and climate. Hundreds of referrals to other researchers, policy analysts, information specialists, or organizations are also facilitated by CDIAC`s staff. This report provides an account of the activities accomplished by CDIAC during the period October 1, 1991 to September 30, 1992. An organizational overview of CDIAC and its staff is supplemented by a detailed description of inquiries ...
Abstract We present the first numerical,-N-body, hydrodynamical, chemical simulations of cosmic structure formation in the framework of non-Gaussian models. We study the impact of primordial non-Gaussianities on early chemistry (e-, H, H+, H-, He, He+, He++, H2, H+2, D, D+, HD, HeH+), molecular and atomic gas cooling, star formation, metal (C, O, Si, Fe, Mg, S) enrichment, Population-III (popIII) and Population-II-I (popII) transition and on the evolution of -visible- objects. We find that non-Gaussianities can have some consequences on baryonic structure formation at very early epochs, but the subsequent evolution at later times washes out any difference among the various models. When assuming reasonable values for primordial non-Gaussian perturbations, it turns out that they are responsi...
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 low-level radioactive waste disposal facility. The model has been ...
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 Taylor-T4 expansion.
The purpose of this paper is to present an analysis of lightning fault detection, location and protection using numeric distance relays applied in high voltage transmission lines, more specifically in the 500 kV transmission lines of CEMIG (Brazilian Energy Utility) between the Vespasiano 2 - Neves 1 (short line - 23.9 km) and Vespasiano 2 - Mesquita (long line - 148.6 km) substations. The analysis was based on the simulations results of numeric distance protective relays on power transmission lines, realized in September 02 to 06, 2002, at Siemens AG's facilities (Erlangen - Germany), using Real Time Digital Simulator (RTDS{sup TM}). Several lightning faults simulations were accomplished, in several conditions of the electrical power system where the protective relays would be installed. The results are presented not only with the times of lightning faults elimination, but also all the functionality of a protection system, including ...
The present work addresses an off-axis Cassegrain optical concentration system. The specific primary collector analyzed, a short focal length parabolic concentrator, is at the University of Florida`s Energy Park. A secondary hyperbolic reflective element was designed to redirect the solar radiation from the primary focal plane to an off-axis target on the polar axis of the primary concentrator. This ground level target will be required for planned experimental work. The analysis was performed using a numerical ray tracing procedure that incorporates both random and systematic errors due to slope and surface irregularities. The optimization process varied secondary element size, curvature, and offset angle, and yielded information required for optimum design. As a single secondary element was found impractical, three elements were designed for use at various time of the year. The numerical analysis predicts that typically 70 to 75 percent of the ...
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 numerical model is compared with other well ...
The results of an experimental investigation of the performance of a small bore engine with a swirl chamber when operating as a dual fuel engine with commercial methane as the gaseous fuel were presented in this paper. The experiment involved using a 3-dimensional computational fluid dynamics model to predict the performance of the engine. A detailed chemical kinetics for the gaseous fuel component, consisting primarily of methane and a reduced detailed chemical kinetics for the diesel fuel while considering the turbulent combustion processes an associated performance of a dual fuel engine with a swirl chamber were incorporated in the simulation. The study experimentally and numerically investigated the effects of changes in the quantities of the liquid fuel pilot and gaseous fuels on the combustion processes, engine performance, cyclic variations, and emissions. The paper discussed the experimental approach and results. It also discussed the simulation of the dual ...
A force measurement technique has been developed for large-scale aerodynamic models with a short test time. The technique is based on direct acceleration measurements, with miniature accelerometers mounted on a test model suspended by wires. Measuring acceleration at two different locations, the technique can eliminate oscillations from natural vibration of the model. The technique was used for drag force measurements on a 3 m long supersonic combustor model in the HIEST free-piston driven shock tunnel. A time resolution of 350 #mu#s is guaranteed during measurements, whose resolution is enough for ms order test time in HIEST. To evaluate measurement reliability and accuracy, measured values were compared with results from a three-dimensional Navier-Stokes numerical simulation. The difference between measured values and numerical simulation values was less than 5%. We conclude that this measurement technique is sufficiently reliable for ...
The choice of vacuum state for a quantum scalar field propagating in a de Sitter spacetime (massive and arbitrarily coupled to the gravitational field) is discussed. The problem of finite-time initial conditions for the mode functions is analyzed, as well as how these determine the vacuum state of the quantum system. The principle guiding the choice of vacuum state is the following: one wants the vacuum contribution to the energy-momentum tensor to contain all the ultraviolet divergent terms, so that the particle creation terms are finite, and covariantly conserved. There is a suitable set of modes (instantaneous adiabatic basis) in which this splitting of the expectation value of the energy-momentum tensor can be carried out. Numerical results are presented for different finite-time initial conditions (m = 0.6, {zeta} = 1/6). The nature of the particle creation effect is described and its relationship to the concept of a horizon crossing time is shown. These ...
The evolution of FRI jets has been long studied in the framework of the FRI-FRII dichotomy. In this paper, we test the present theoretical and observational models via a relativistic numerical simulation of the jets in the radio galaxy 3C 31. We use the parameters derived from the modelling presented by \\cite{lb02a,lb02b} as input parameters for the simulation of the evolution of the source, thus assuming that they have not varied over the lifetime of the source. We simulate about 10 % of the total lifetime of the jets in 3C 31. Realistic density and pressure gradients for the atmosphere are used. The simulation includes an equation of state for a two-component relativistic gas that allows a separate treatment of leptonic and baryonic matter. We compare our results with the modelling of the observational data of the source. Our results show that the bow shock evolves self-similarly at a quasi-constant speed, with slight deceleration by the end of the simulation, ...
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 experiments. The ...
Abstract in english In this work the air flow in a furnace was computationally investigated. The furnace, for which experimental test data are available, is composed of a movable block burner connected to a cylindrical combustion chamber by a conical quarl. The apertures between the movable and the fixed blocks of the burner determine the ratio of the tangential to the radial air streams supplied to the furnace. Three different positions of the movable blocks were studied at this time. A thr (more) ee-dimensional investigation was performed by means of the finite volume method. The numerical grid was developed by the multiblock technique. The turbulence phenomenon was addressed by the RNG k-epsilon model. Profiles of the axial, tangential and radial velocities in the combustion chamber were outlined. The map of the predicted axial velocity in the combustion chamber was compared with a map of the experimental axial velocity. The internal space of the furnace was ...
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 introduced to find ...
We studied two-photon exchange for elastic electron-proton scattering at low Q{sup 2}. Compact approximate formulas for the amplitudes were obtained. Numerical calculations were done for Q{sup 2}{<=}0.1 GeV{sup 2} with several realistic form factor parametrizations, yielding similar results. They indicate that the corrections to the magnetic form factor can visibly affect the cross-section and proton radii. For low-Q{sup 2} electron-neutron scattering two-photon exchange corrections were shown to be negligibly small00.
Different calculated models are presented for turbulent break-off flows and their classification, reflecting the sequence of historical development. The study was done based on equations of viscous liquid of the Navier-Stokes type with development of special phenomenological models of turbulence which take into consideration real properties of the break-off flows based on simpler models of flow presented in the work. In order to calculate two-dimensional turbulent flows, a method of viscous-nonviscous interaction is used. It employs numerical solutions for nonviscous flow and integrated methods of calculating the dissipative region. This method can be extended for calculating the transonic break-off flows, and also break-off in an incompressible fluid when there is cavitation.
The phase space beam analyzer is a measurement instrument that is applied in laser technology to perform analyses of the spatial and angular distribution of rays. We are interested in this instrument as a means to characterize non-coherent light sources. In this context, a closer look at the tolerances of this optical instrument was considered useful. Having a so-called quadrupole lens as a key element, the phase space beam analyzer is a device that features anamorphic optical properties. To describe these anamorphic properties, recurrence was made to a description by extended ray-transfer matrices. This formalism allows for an analysis of the alignment tolerances of the phase space beam analyzer and facilitates a study of the sensitivities of the instrument. The analysis is complemented using numerical ray tracing.
The theoretical hardness against total dose of the six-transistor SRAM cell is investigated in detail. An explicit analytical expression of the maximum tolerable threshold voltage shift is derived for two cross-coupled inverters. A numerical method is used to explore the hardness of the read and write operations. Both N- and P-channel access transistors designs are considered and their respective advantages are compared. The study points out that the radiation hardness mainly relies on the technology. Results obtained with the very robust Gate-All-Around process are finally presented.
The load-carrying capacity and ultimate deformation capacity of timber-concrete composite systems can be significantly influenced by the ductility of the connection between the two materials. The use of more ductile connections can increase the load-carrying capacity of the composite system as well as its ultimate deformation capacity. In this work, the potential increase that might be expected for these two parameters due to the use of ductile connections will be assessed through numerical simulations, taking the non-linear behaviour of the connections into account. Furthermore, the connection ductility required to achieve the maximum load capacity depends on the mechanical properties of the connection as well as on the geometric and mechanical properties of the composite system. There ar...
We numerically demonstrate selective near-field localization determined by the polarization state of a single emitter coupled to plasmonic nano-cluster. Seven gold nanospheres are carefully arranged such that up to ten polarization states of the single emitter, including linear, circular, and elliptical polarizations, can be distinguished via the distinct field localization in four gaps. The ability to transform polarization state into field spatial localization may find application for single emitter polarization analysis.
Axial offset anomaly (AOA) is caused by the deposition of crud on the fuel cladding of a PWR. When significant levels of crud build up on the cladding, boron can accumulate in the pores of the crud as a concentrated solution or solid phase, and cause the flux depression. Numerous studies have been conducted on the primary water chemistry to reduce the amount of crud in the primary circuit to avoid radioactivity buildup and unexpected power transition in the plant. However, experiments on the crud are restricted in the laboratory because the crud is a highly radioactive material. The objective of this study is to develop a test method for simulating the deposition of crud in a nuclear power plant
We consider non-relativistic systems in quantum mechanics interacting through the Coulomb potential, and discuss the existence of bound states which are stable against spontaneous dissociation into smaller atoms or ions. We review the studies that have been made of specific mass configurations and also the properties of the domain of stability in the space of masses or inverse masses. These rigorous results are supplemented by numerical investigations using accurate variational methods. A section is devoted to systems of three arbitrary charges and another to molecules in a world with two space-dimensions.
This report details an investigation into the efficacy of two approaches to solving the radiation diffusion equation within a radiation hydrodynamic simulation. Because leading-edge scientific computing platforms have evolved from large single-node vector processors to parallel aggregates containing tens to thousands of individual CPU's, the ability of an algorithm to maintain high compute efficiency when distributed over a large array of nodes is critically important. The viability of an algorithm thus hinges upon the tripartite question of numerical accuracy, total time to solution, and parallel efficiency.
The Idaho National Engineering Laboratory (INEL) has been conducting geothermal reservoir research and testing sponsored by the US Department of Energy (DOE) since 1983. The INEL research program is primarily aimed at the development of reservoir engineering techniques for fractured geothermal reservoirs. Numerical methods have been developed which allow the simulation of fluid flow and heat transfer in complex fractured reservoirs. Sensitivity studies have illustrated the importance of incorporating the influence of fractures in reservoir simulations. Related efforts include fracture characterization, geochemical reaction kinetics and field testing.
Nonlinear distortion added by the loudspeaker in a hearing aid lowers the signal-to-noise ratio and may degrade the hearing aid user's ability to understand speech. The balancedarmature- type loudspeakers, predominantly used in hearing aids, are inherently nonlinear devices, as any displacement of the loudspeaker diaphragm inevitably changes the magnetic and electrical characteristics of the loudspeaker. A numerical time-domain model capable of describing these nonlinearities is presented. By simulation it is demonstrated how the output distortion could potentially be reduced significantly through careful design of the mechanical properties of the armature.
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.
The author briefly discusses definition of terms, gives an introduction to measurement techniques and describes the characteristics of various low-frequency fields and their causes using typical examples: natural electric fields (thunderstroms), natural magnetic fields, technical electric constant fields (urban transportation, households), static magnetic fields (urban transportation, nuclear magnetic resonance imaging), technical electric alternating fields (high-voltage transmission lines, households), and magnetic alternating fields (high-voltage transmission lines). The author discusses both occupational exposure and that of the general public while underpinning his statements by numerous tables, measurement diagrams and charts. (Uhe).
We have developed a novel laser ray-tracing method to measure aberrations in optical systems. It consists of delivering narrow laser pencils (by a laser scanner), recording the spots that are formed on the image plane (with a CCD camera), and computing the position of each centroid. This approach could be considered an experimental (approximate) implementation of standard numerical ray tracing. Several tests and experiments, including a direct comparison with a Hartmann-Shack wave-front sensor, provided highly satisfactory results that confirmed the validity of the method and revealed potential advantages. PMID:18073906
Trevor Keel speaks to Cara Sutton, Assistant Commissioning Editor Within the World Gold Council's technology sector, Trevor Keel is responsible for a range of activities in relation to the development of new gold-based technologies. He has particular interest in the use of gold in the fields of medicine and diagnostics, catalysis and nanotechnology. Previously, Trevor was a principal chemist at GlaxoSmithKline where he was involved in numerous New Chemical Entity projects. Dr Keel has an extensive knowledge of nanotechnology, having completed a PhD in Pharmaceutical Nanotechnology at the University of Nottingham, Nottingham, UK.
The measurement of electron densities using interferometry and refraction measurements in plasmas with cross-sections where the electron density contours are concentric ellipses is examined. Transforms are found for both interferometrically deduced optical path-length differences and refraction angle data obtained from elliptical cross-section plasmas. The transformed data can be inverted to give electron densities using standard techniques developed for circular cross-section plasmas. Refraction of light in elliptical cross-section plasmas is examined using numerical ray tracing.
We discuss the inclusive production of D{sup *{+-}} mesons in {gamma}p collisions at DESY HERA, based on a calculation at next-to-leading order in the general-mass variable-flavor-number scheme. In this approach, MS subtraction is applied in such a way that large logarithmic corrections are resummed in universal parton distribution and fragmentation functions and finite mass terms are taken into account. We present detailed numerical results for a comparison with data obtained at HERA and discuss various sources of theoretical uncertainties. (orig.)
We study the production of gravitational waves from cosmic domain walls created during phase transition in the early universe. We investigate the process of formation and evolution of domain walls by running three dimensional lattice simulations. If we introduce an approximate discrete symmetry, walls become metastable and finally disappear. We calculate the spectrum of gravitational waves produced by collapsing metastable domain walls. Extrapolating the numerical results, we find the signal of gravitational waves produced by domain walls whose energy scale is around 10^10-10^12GeV will be observable in the next generation gravitational wave interferometers.
The primary objective is to perform a fundamental study of filtration of nanoparticles, and to obtain filtration knowledge necessary to design particle collection devices/systems for nanoparticle processing and for preventing nanoparticle emissions into the environment. The research covered a wide area relevant to nanoparticle filtration, under these main topics: (1) nanoparticle filtration and molecular dynamics simulation, (2) nanoparticle virtual impactor, (3) particle transport under low pressure, and (4) development of a high-throughput nanoparticle generator. A number of novel tools and numerical models have been developed under the DOE support.
Recent surveys of supergiants stars in the LMC indicate that the post-main-sequence region of the colour-magnitude diagram is well populated, although numerical evolution of massive stars with normal surface hydrogen indicates to the contrary. Supergiant stars having surface enrichment of helium acquired for example from a previous phase of accretion from a binary companion, however, evolve in a way so that the evolved models and observed data are consistent. We compare the available data with computed evolutionary tracks of massive stars of metallicity relevant to the LMC with and without helium-enriched envelopes and conclude that a large fraction of supergiant stars may occur in binaries. (author).
A methodology of tsunami PSA was developed in this study. A tsunami PSA consists of tsunami hazard analysis, tsunami fragility analysis and system analysis. In the case of tsunami hazard analysis, evaluation of tsunami return period is major task. For the evaluation of tsunami return period, numerical analysis and empirical method can be applied. The application of this method was applied to a nuclear power plant, Ulchin 56 NPP, which is located in the east coast of Korean peninsula. Through this study, whole tsunami PSA working procedure was established and example calculation was performed for one of real nuclear power plant in Korea
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.
Long cylindrical mirror sections with an 'Alliptical' cross- section and an interfoci distance of approximately 1 m can be produced by means of diamond fly cutting. However, because of the finite tool radius, the generated profile is basically not an ellipse, but an equidistant ellipse or 'Allipse'. By numerical ray tracing analysis two unique conjugated 'A-points' can be found that are optically nearly equivalent with pure elliptical focal points. A reversed modeling procedure is used to predict the optimum machine configuration for producing just the required ellipse sections. Optical figure quality, focal line position and straightness are tested by a simple deflectometer set-up.
A brief review is presented of transport in the boundary region of magnetized plasmas by blob-like filaments. Such structures have enhanced levels of particles and heat, are elongated along the magnetic field lines and are localized in the drift plane across the field. The motion of an isolated blob structure is described in some detail and the contribution of such filaments to turbulence-driven transport are discussed. Results are presented from numerical simulations and probe measurements in tokamak plasmas. An interpretation is given of the measured dependence of particle density and transport on experimental control parameters in the scrape-off layer. (author)
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.
This work contains 27 chapters which cover recent advances in fracture mechanics applied to cement-based composites: concrete, reinforced concrete, and fiber reinforced concrete. The book, based on a NATO Advanced Study Institute, is a comprehensive survey of treating nonlinearity associated with crack growth in brittle materials such as cement composites, rocks and ice. The book is divided into 8 major themes: advances in nonlinear fracture mechanics; fracture processes; numerical modelling of fracture; experimental methods of determining fracture parameters; damage and continuum modelling; strain-rate and dynamic effects of crack propagation; stress-corrosion, time and temperature effects on fracture; and implications for concrete structures.
Diffusion of radon in a bed of activated charcoal is described by diffusion equations. An analytical solution of these equations is presented for the case of constant radon concentration in the atmosphere. The solutions are given separately for short term and long term exposure. An analytical form of the calibration constant f for long term exposure and constant radon concentration in air, was found. A numerical method and computer code based on the method of finite elements is developed for the case of variable radon concentration in air. This program simulates radon adsorption by the activated charcoal bed, enabling determination of sensitivity. The dependence of sensitivity on different parameters, such as temperature, thickness of the charcoal, etc. was studied using this program. (orig.) 10 refs.
This paper presents a supernetwork equilibrium model integrating supply chain networks with a transport network, namely, a supply chain-transport supernetwork equilibrium model. The model takes into account the behaviour of freight carriers and transport network users to endogenously determine the transport costs generated in the supply chain networks. The interaction between transport network and supply chain networks can also be examined. Results of the numerical tests reveal that the improvement of transport network could enhance the efficiency of supply chain networks. The paper makes contributions to modelling of supply chain networks as well as to that of transport networks.