The authors are developing MP-TOUGH2 for exploiting massively parallel computers. The goals of this effort are to (1) create a data-parallel subsurface transport code for solving larger problems than currently practical on workstations, (2) write portable code that can take advantage of scalability to run on machines with more processors, and (3) minimize the necessity for additional validation and verification of the resulting code. The initial strategy they have followed is to focus on optimizing the generic and time-consuming task of linear equation solution while leaving the bulk of TOUGH2 unmodified. In so doing, they have implemented a massively parallel direct solver (MPDS) that takes advantage of the banded structure of TOUGH2 Jacobian matrices. The authors have compared timings of the iterative conjugate gradient solvers DSLUBC, DSLUCS, and DSLUGM written in Fortran77 for the front end with the MPDS which uses the ...
British Nuclear Fuels plc owns and operates the near-surface Drigg disposal facility for low level radioactive waste. The long-term performance of the site is modelled by a suite of computer codes called MONDRIAN. One of the modules of MONDRIAN deals with the transport of radionuclides through the environment, and this paper reports on the current status of this module (BIOS). We derive the basic set of working equations from first principles and show clearly how the approximate nature of the final equations is arrived at. This is done by an averaging process leading to compartments, in and out of which radionuclides, solids and water can flow. The equations allow radioactive decay chains and an arbitrary number of compartments. There is also the facility to deal with changes in the rate coefficients, thereby simulating different environmental states. It is also possible to include the creation of new ...
In this dissertation, we design and analyze efficient numerical methods for obtaining accurate solutions to model problems arising in fuel cells. A basic fuel cell model consists of five principles of conservation, namely, mass, momentum, species, charges (electrons and ions), and thermal energy. Overall, transport equations couple with electrochemical processes through source terms to describe reaction kinetics and electro-osmotic drag in the polymer electrolyte. To model multiphase species transport in the porous media and the gas channel of fuel cells, we consider a multiphase mixture model framework. The diffusivity of the two-phase mixture water conservation equation in this model is nonlinear, discontinuous, and degenerate. To handle this difficulty, we developed efficient and fast nonlinear iterative solvers based on the Kirchhoff transformation and nonlinear Dirichlet-Neumann domain ...
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 aim of this paper is to illustrate the validity and efficiency of iterative methods for solving large linear systems arising from the finite element discretization of the equation governing conduction-controlled solidification processes. Starting from the basic enthalpy equation, two alternative formulations are obtained and fixed-grid finite element discretizations are developed. These discretizations yield a set of nonlinear equations that are linearized using the Newton-Raphson scheme. The linearized equations are used as a basis for evaluating different iterative methods of the conjugated gradient type. Symmetric scaling and incomplete factorization preconditioning of the linear equations are used to improve the convergence properties of the iterative methods. Vectorization and parallelization are also employed to make full use of the CRAY-2 supercomputer. The results ...
A stable, accurate, and efficient implementation of MacCormack's explicit algorithm for the Parabolized Navier-Stokes equations is demonstrated. The familiar problem of decoding the conservative axial flux vector is solved, resulting in accurate, smooth dependent variable profiles through the viscous-layer sonic line. Source terms due to transformation of the parabolized governing equations into the computational plane and the equations into the computational plane and the resulting metric differencing have been identified and eliminated through inclusion of appropriate geometric conservation law terms. Test cases computed include two- and three-dimensional supersonic and hypersonic flow at laminar and turbulent Reynolds numbers. The computed results demonstrate very good agreement with experiment and with solutions of the full Navier-Stokes equations. Computational times required for the ...
Production decline-curve analysis is accepted for prediction of future performance of oil and gas wells and fields. One type of decline-curve analysis involves the hyperbolic-decline-curve equation. The challenge with this equation is to determine simultaneously three parameters. The new method creates a computer spreadsheet containing a set of production data. Then, the specific arrangement of cells and columns containing data and/or formulas allows for statistical fits of the data. Finally, setting-up the nonlinear solver allows for optimization of the three parameters by maximization of the square of the regression coefficient. Eight field cases demonstrated the repeatable and rapid curve-fitting of the method.
We show that the eigenvalue equations describing a cylindrical ideal magnetophydrodynamicsw (MHD) plasma interacting with a thin resistive wall can be put into the standard mathematical form: ??? = ??? ?. This is accomplished by using a finite element basis for the plasma, and by adding an extra degree of freedom corresponding to the electrical current in the thin wall. The standard form allows the use of linear eigenvalue solvers, without additional interations, to compute the complete spectrum of plasma modes in the presence of a surrounding restrictive wall at arbitrary separation. We show that our method recovers standard results in the limits of (1) an infinitely resistive wall (no wall), and (2) a zero resistance wall (ideal wall).
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.
The Kadomtsev-Petviashvili equation describes nonlinear dispersive waves which travel mainly in one direction, generalizing the Korteweg-de Vries equation for purely uni-directional waves. In this Letter we derive an improved KP-equation that has exact dispersion in the main propagation direction and that is accurate in second order of the wave height. Moreover, different from the KP-equation, this new equation is also valid for waves on deep water. These properties are inherited from the AB-equation (E. van Groesen, Andonowati, 2007 ) which is the unidirectional improvement of the KdV equation. The derivation of the equation uses the variational formulation of surface water waves, and inherits the basic Hamiltonian structure.
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 equationssolver which ...
We study the well-posedness of the initial value problem for a wide class of singular evolution equations. We prove a general well-posedness theorem under three assumptions easy to check: the first controls the singular part of the equation, the second the behavior of the nonlinearities, and the third one assumes that an energy estimate can be found for the linearized system. We allow losses of derivatives in this energy estimate and therefore construct a solution by a Nash-Moser iterative scheme. As an application to this general theorem, we prove the well-posedness of the Serre and Green-Naghdi equation and discuss the problem of their validity as asymptotic models for the water-waves equations.
In order to reduce the air concentration of (sup 3)H in the reactor buiIding of Wolsung Heavy Water Reactor, a computer code for estimation of adsorption behavior was programmed based on an equation derived for analysis of water vapor adsorption, and a ba...
A new multigrid algorithm based on the method of self-correction for the solution of elliptic problems is described. The method exploits information contained in the residual to dynamically modify the source term (right-hand side) of the elliptic problem. It is shown that the self-correcting solver is more efficient at damping the short wavelength modes of the algebraic error than its standard equivalent. When used in conjunction with a multigrid method, the resulting solver displays an improved convergence rate with no additional computational work.
An essential part of fire risk assessment is the analysis of fire hazards and fire propagation. In this work, models and tools for two different aspects of numerical fire simulation have been developed. The primary objectives have been firstly to investigate the possibility of exploiting state-of-the-art fire models within probabilistic fire risk assessments and secondly to develop a computationally efficient solver of thermal radiation for the Fire Dynamics Simulator (FDS) code. In the first part of the work, an engineering tool for probabilistic fire risk assessment has been developed. The tool can be used to perform Monte Carlo simulations of fires and is called the Probabilistic Fire Simulator (PFS). In Monte Carlo simulation, the simulations are repeated multiple times, covering the whole range of variability of the input parameters and thus resulting in a distribution of results covering what can be expected in reality. In practical applications, advanced ...
This paper describes the development of a computational multiphase fluid dynamics (CMFD) model of the Fischer Tropsch (FT) process in a Slurry Bubble Column Reactor (SBCR). The CMFD model is fundamentally based which allows it to be applied to different industrial processes and reactor geometries. The NPHASE CMFD solver [1] is used as the robust computational platform. Results from the CMFD model include gas distribution, species concentration profiles, and local temperatures within the SBCR. This type of model can provide valuable information for process design, operations and troubleshooting of FT plants. An ensemble-averaged, turbulent, multi-fluid solution algorithm for the multiphase, reacting flow with heat transfer was employed. Mechanistic models applicable to churn turbulent flow have been developed to provide a fundamentally based closure set for the equations. In this four-field model formulation, two of the fields are used to track ...
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 ...
The aim of this study concerns the use of numerical methods for the resolution of the Reynolds Averaged Navier Stokes equations adapted to the simulation of the cooling of the trailing edge of a stator in a high pressure turbine. These methods, based on the elsA solver developed at ONERA, use a four steps Runge Kutta time discretization scheme and a Jameson centered space discretization scheme. The scheme is applied through a finite volume approach on control volume centered on the cells of a multi-block structured mesh. Turbulence is simulated either through the algebraic Michel model, or through the one-transport-equation Spalart-Allmaras model, or through the two-transport-equations k 1, k {omega} and k {epsilon} models, and through ASM model. A simulation of the flow in a bidimensional stator, without cooling, is carried out. The cooling, which is realized with trailing edge slots, is then simulated ...
The interfacial area transport equation dynamically models the two-phase flow regime transitions and predicts continuous change of the interfacial area concentration along the flow field. Hence, when employed in the numerical thermal-hydraulic system analysis codes, it eliminates artificial bifurcations stemming from the use of the static flow regime transition criteria. Accounting for the substantial differences in the transport phenomena of various sizes of bubbles, the two-group interfacial area transport equations have been developed. The group 1 equation describes the transport of small-dispersed bubbles that are either distorted or spherical in shapes, and the group 2 equation describes the transport of large cap, slug or churn-turbulent bubbles. The source and sink terms in the right-hand-side of the transport equations have been established by mechanistically modeling the ...
A series of studies conducted to identify the factors that block and unblock problem solving is described. Through the construction of an isomorph of the classic "water jar" problems developed by A. S. Luchins (1942) as a dynamic graphic micro-world, several factors involved in producing blocked states were identified. Subjects included 10 individuals and 10 pairs of subjects. By comparing the behavior of individuals tackling the "missionaries and cannibals" problem to that of pairs of subjects solving this problem, a study identified means by which problem solvers operating in a social context are able to overcome blocks that discourage individuals. These studies point to the importance of "reflection" (evaluation of problem-solving results) for flexible problem solving. This research suggests that teaching students to analyze what they have done will help them develop flexibility in using a new approach when blocked. The results may also ...
... of the cycle shows that the separation area ea. ... with suitable flow solvers based on the Reynolds ... AGARD 75th Fluid Dynamics Panel Meeting and ...
What will be discussed in this report represents a framework upon which multiphase and other real physical effects can be built. Chemical models of increasing complexity are envisioned and this methodology can provide a tool for evaluating new ideas against known experimental data. The recent work to be reported here addresses the multiphase issue of temperature deviation between phases undergoing chemical and heat transport processes. Modeling of the LLNL ODTX experiment will be performed with FLUENT, a commercially available computational fluid dynamics (CFD) code. FLUENT solves flows in 2D or 3D in Cartesian, cylindrical, or general curvilinear coordinates, with steady-state of fully time-dependent analysis. Multiphase flows in which two or more continuous phases are present can be solved with arbitrary volumetric sources of heat, mass, momentum, and chemical species applied through user-defined FORTRAN subroutines. FLUENT models these of phenomena by solving the conservation ...
An equation is described to estimate the value of the liquid-phase mass-transfer coefficient of aqueous solutions of surfactants from its value for pure water. The form of the equation is obtained from Levich's theory, and it is adjusted to experiments of absorption of CO/sub 2/ in aqueous solutions of five pure surfactants. The reduction of the absorption rate is correlated with the interfacial pressure of the solution.
The interfacial area transport equation dynamically models the two-phase flow regime transitions and predicts continuous change of the interfacial area concentration along the flow field. Hence, when employed in the numerical thermal-hydraulic system analysis codes, it eliminates artificial bifurcations stemming from the use of the static flow regime transition criteria. Accounting for the substantial differences in the transport phenomena of various sizes of bubbles, the two-group interfacial area transport equations have been developed. The group 1 equation describes the transport of small-dispersed bubbles that are either distorted or spherical in shapes, and the group 2 equation describes the transport of large cap, slug or churn-turbulent bubbles. The source and sink terms in the right-hand-side of the transport equations have been established by mechanistically modeling the ...
The performance of five well-known photosynthesis-light equations is compared by presenting a wide range of solutions in the form of dimensionless nomographs for the case where photosynthetically available radiation (PAR) reduces exponentially down the water column and is distributed sinusoidally through the photoperiod. These provide a simple means of calculating daily photosynthesis at any depth (e.g. at a benthic layer), or through a water column, avoiding the need to perform complex integrations. An examination is made of the accuracy of common approximate methods for calculating daily photosynthesis, assuming constant PAR. For optically deep water a modification is proposed to Talling's planimetric solution, to enable daily photosynthesis to be calculated more accurately, yet simply, over the whole range of possible PAR values. The errors induced by approximating the daily PAR distribution ...
Dippold, V., S. Mohler , Validation of the Wind-US Unstructured Flow Solver .... Hamed, A. and A. Mohamed, Assessment of Shock Induced Flow Separation and ...... Nozzles for Hypersonic Propulsion, NASA CR 185197, AGARD Symposium ...
We assessed the age-related changes in magnetization transfer ratio (MTR) value and the ratio of water content in rat brain (frontal lobe and pons) during normal development. All studies were done on 42 normal Wistar rats. MTR was calculated by using the equation (M_o_f_f-M_o_n)/M_o_f_f, where M_o_f_f and M_o_n were the measured signal intensities on the conventional and MT-prepared SPGR images, respectively. The ratios of water content were measured by comparing the wet weight with the dry weight of the frontal lobe and the pons, which were extracted from a rat brain soon after MTR measurements. The ratio of water content was calculated by using the equation 100 x (wet weight-dry weight)/wet weight (%). The ratio of water content decreased markedly during the first 30 days, but the decrease was very slight there after both in the frontal lobe and the pons. The ...
A new energy concept known as centrifugal potential energy was discussed. This new energy concept is capable of increasing the pressure, temperature and enthalpy of a fluid, without having to apply work or heat transfer to the fluid. It occurs through a change in the centrifugal potential energy of the flowing fluid in a rotating frame of reference or a centrifugal force field, where work is performed internally by the centrifugal weight of the fluid. This energy concept has resulted in new energy equations, such as the Rotational Frame Bernoulli's Equation for liquids and the Rotational Frame Steady-Flow Energy Equation for gases. Applications of these equations have been incorporated into the design of centrifugal field pumps and compressors. Rather than compressing a fluid with a physical load transfer, these devices can compress a fluid via the effect of centrifugal force applied to the ...
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 ...
We present the runtime comparison of the two versions of Super LU{_}DIST, using up to 128 processors of the IBM SP at NERSC. One version provides the global input interface, and another provides the distributed input interface. The comparison includes the total runtime of the solver with both 32-bit and 64-bit addressing modes, the time breakdown for different phases of the solver. We also present an in-depth comparison off our sparse matrix-vector multiplication methods in the context of iterative refinement. Finally, we describe our Fortran 90 interface that enhances the usability of the software.
Full text of publication follows:In the current thermal-hydraulic system analysis codes using the two-fluid model, the empirical correlations that are based on the two-phase flow regimes and regime transition criteria are being employed as closure relations for the interfacial transfer terms. Due to its inherent shortcomings, however, such static correlations are inaccurate and present serious problems in the numerical analysis. In view of this, a new dynamic approach employing the interfacial area transport equation has been studied. The interfacial area transport equation dynamically models the two-phase flow regime transitions and predicts continuous change of the interfacial area concentration along the flow field. Hence, when employed in the thermal-hydraulic system analysis codes, it eliminates artificial bifurcations stemming from the use of the static flow regime transition criteria. Therefore, the interfacial area transport ...
Abstract The available surface equation of state for pure pulmonary surfactant monolayers is generalised to binary mixture monolayers by introducing a group of parameters, i(i-=-1-4) in the form of $x_{r}^{\\beta _{i} } $ to express the influence of the components in new mixing rules and a new factor, I $\\left[ { = \\left( {\\prod\\limits_{i = 1}^{4} {\\beta _{i} } } \\right)^{{\\raise0.5ex\\hbox{$\\scriptstyle 1$}\\kern-0.1em/\\kern-0.15em\\lower0.25ex\\hbox{$\\scriptstyle {4}$}}} } \\right]$, is defined to represent the interaction intensity between two different components. The ---A isotherms getting by the surface equations of state agree with the experimental data for protein-lipid binary monolayers, and the average deviation is about 11.41%. The result shows the order of the interaction intensity b...
Real-time neutron radiography was demonstrated to be a viable and useful tool in studying water flow in unsaturated media. The experiment was designed to test the method for following water flow. The wetting front advance, water-flow patterns, and neutron attenuation density were observed in a 7.7-cm diameter sand column. The advance was fit with a simple i = A t/sup B/ equation yielding an excellent correlation. Theoretical calculations presented estimate the minimum and maximum detectable water content in various sizes of columns. The ability of real-time neutron radiography to allow both qualitative and quantitative spatial and time measurements to be made in an nonintrusive manner was demonstrated.
Real-time neutron radiography was demonstrated to be a viable and useful tool in studying water flow in unsaturated media. The experiment was designed to test the method for following water flow. The wetting front advance, water-flow patterns, and neutron attenuation density were observed in a 7.7-cm diameter sand column. The advance was fit with a simple i = A t/sup B/ equation yielding an excellent correlation. Theoretical calculations presented estimate the minimum and maximum detectable water content in various sizes of columns. The ability of real-time neutron radiography to allow both qualitative and quantitative spatial and time measurements to be made in an nonintrusive manner was demonstrated.
For the layout of a pressure suppression system of a light water cooled reactor (boiling water reactor) it is important to know the time dependent behavior of the vent clearing transient after a loss-of-coolant accident for two main reasons: time of the end of the vent clearing transient influences strongly the pressure and temperature maxima in the drywell and wetwell. Time-dependent behavior of the vent clearing transient influences pressure loads in the condensation pool of the wetwell and therefore pressure induced stresses to the structure. The time-dependent behavior of the water masses in the vent pipes and wetwell are described by the basic equations for a nonstationary incompressible friction flow: momentum equation, continuity equation and a correlation for the variation of the state of the gas volume in the wetwell above the water ...
A method of modeling freeze-thaw cycles of naturally deposited snowpacks is presented. The model involves the Stefan condition as an independent governing equation on the exterior moving boundary to calculate snowpack thinning, flow of water through a variably saturated layered porous medium as described by the Richards equation, and heat conduction with a phase change. The heat conduction problem was treated in two ways. Local heat conduction between a snow grain and its surrounding water film was treated by using a simple energy balance. Global heat conduction with a phase change (the Stefan problem) was introduced to calculate the space-time temperature distribution. In order to handle multiple interior moving boundaries, a specific form of the enthalpy formulation was used for heat conduction with a phase change. Changing material properties were considered according to the calculated meltwater ...
The lakes formed in the Amazon floodplain play a key role in sustaining the highest populational density in the Amazon Basin because of its high primary productivity. The water volume in these lakes vary greatly during the year and possible water sources are the feeding river during high water state of the hydrographs and local inputs from its drainage basin. As these sources are chemically distinct, to understand the nutrient dynamics in these lakes, it is very important to know the relative water contribution of different water sources. In this paper we use "1"8O as a conservative to know the water sources for the lakes and our hypothesis is that: 1) In the high water stage lake water is derived predominantly from its feeding river. Therefore their #delta#"1"8O should be similar. 2) In the low ...
The objective of this study is to develop an artificial neural network (ANN) model to predict the thermal conductivity of ethylene glycol-water solutions based on experimentally measured variables. The thermal conductivity of solutions at different concentrations and various temperatures was measured using the cylindrical cell method that physical properties of the solution are being determined fills the annular space between two concentric cylinders. During the experiment, heat flows in the radial direction outwards through the test liquid filled in the annual gap to cooling water. In the steady state, conduction inside the cell was described by the Fourier equation in cylindrical coordinates, with boundary conditions corresponding to heat transfer between the solution and cooling water. ...
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 waterequations (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...
An electrochemistry model was developed to analyse the J-V characteristics of a Proton Exchange Membrane (PEM) water electrolyzer for hydrogen production. The Butler-Volmer equation and water transport characteristics through electrolyte membrane were employed to simulate the electrode activation over-potential and membrane ohmic over-potential, respectively. The modeling results are found to agree reasonably well with experimental data published in the literature. The parametric simulations show that the ohmic over-potential is relatively small with typical water content in the membrane. Compared with the cathode over-potential, the anode over-potential is more significant and constitutes the major source of voltage loss. The high anode over-potential is due to the relatively slow oxidation kinetics, which is related to anode material property and microstructure. This model can be integrated with a ...
This paper describes a comprehensive model of wastewater treatment in secondary facultative ponds, which combines 3D hydrodynamics with a mechanistic water quality model. The hydrodynamics are based on the Navier-Stokes equation for incompressible fluids under shallow water and Boussinesq assumptions capturing the flow dynamics along length, breadth and depth of the pond. The water quality sub model is based on the Activated Sludge Model (ASM) concept, describing COD and nutrient removal as function of bacterial growth following Monod kinetics, except for Escherichia coli removal, which was modelled as first order decay. The model was implemented in the Delft3D software and was used to evaluate the effect of wind and the addition of baffles on the water flow pattern, temperature profiles i...
In this paper the thermal design of a horizontal steam generator for the Indian PHWR nuclear power plant is described. The main attraction is absence of tube sheet and use of stainless steel 'U' tubes. It is emphasised that with appropriate water chemistry it is possible to use stainless steel tubes, which is many times cheaper than the Incoloy tubes used elsewhere. The design approach, applicable equation for the design and the results of computation in the form of heat transfer area and some important dimensions of the steam generator are presented.
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.
The tritium in transpired water of plant leaves from four locations within site area (2=0.92) with that of measured values indicating that at Tarapur, the contribution to TFWT of plant is only from air activity and pick up through soil route is negligible during the period of study. (author)
The Water Frame work Directive (December 2000) proposes integrated water management in regard to both quality and quality. A mathematical model has therefore been produced at the CEIT and EPTISA research centre, in collaboration with the Ebro Hydrographic. Confederation (CHE), to predict water quality and thus begin integrated management of it in the river basins. A description is given of the characteristics of the simulation tool. The hydraulic model is based on the numerical (weighted, four-point implicit) resolution of the complete Saint Venant equations. The quality model is based on the IWA River Water Quality Model Number 1 which has clear advantages (consistency, mass balance and easy integration with the biological models of the waste water treatment plants and the collector networks compared to traditional models. The experimental validation, which ...
Beyond a great number of applications in mineral process, the hydrocyclone separation has been used recently in an increasing number of fields, such as: environmental engineering, petrochemical engineering, food engineering, electrochemical engineering, pulp bioengineering, processes and paper industries, among others. Nowadays, the produced water is treated and then reinjected in the surface and also in return to the formation. In the offshore platforms, the produced water can be directly discarded in the ocean, since the dispersed phase concentration (oil) is according to Brazilian regulating agency limits. The basic principle used in hydrocyclone separation is the centrifugal sedimentation machine: the suspended particles are submitted to centrifugal acceleration, which separates them from the fluid. The basic geometry of a hydrocyclone consists of four parts: a cylindrical section for feeding, a conical section for swirling, an overflow and ...
Level set method based segmentation provides an efficient tool for topological and geometrical shape handling. Conventional level set surfaces are only $C^0$ continuous since the level set evolution involves linear interpolation to compute derivatives. Bajaj et al. present a higher order method to evaluate level set surfaces that are $C^2$ continuous, but are slow due to high computational burden. In this paper, we provide a higher order GPU based solver for fast and efficient segmentation of large volumetric images. We also extend the higher order method to multi-domain segmentation. Our streaming solver is efficient in memory usage.
A dense Pd-Ag membrane reactor (MR) with 100% hydrogen selectivity packed with either Rh/La2O3 or Rh/La2O3-SiO2 as catalysts was used to carry out the dry reforming of methane. The membrane reactor simulation was performed using a well-known reactor model. For this purpose, we employed the equations derived from complete kinetic studies of the dry reforming of methane reaction in connection with both catalysts. In addition, we developed the kinetic equation for the reverse water gas shift reaction (RWGS). The combination of detailed kinetic studies with the measured permeation flux for the Pd-Ag membrane allowed a complete comparison between experimental and simulated operation variables. The variables studied for both catalysts were methane conversion and hydrogen permeation as a function...
The heating of water layer using microwave oven with a rectangular waveguide has been studied both numerically and experimentally. The mathematical model is validated with the experimental data. The transient Maxwell's equations are solved by using the Finite Difference Time Domain (FDTD) method to describe the electromagnetic field inside the waveguide and sample. The temperature profile and velocity field within sample are determined by the solutions of the momentum, energy and Maxwell's equations. In this study, the effects of physical parameters, e.g. microwave power level, placement of sample inside the waveguide, volume of sample, are studied. The distribution of electric field, temperature profile and velocity field are presented in details. The results show good agreement between s...
A one-dimensional, steady-state, two-phase direct methanol fuel cell (DMFC) model is developed to precisely investigate complex physiochemical phenomena inside DMFCs. In this model, two-phase species transport through the porous components of a DMFC is formulated based on Maxwell-Stefan multi-component diffusion equations, while capillary-induced liquid flow in the porous media is described by Darcy's equation. In addition, the model fully accounts for water and methanol crossover through the membrane, which is driven by the effects of electro-osmotic drag, diffusion, and the hydraulic pressure gradient. The developed model is validated against readily available experimental data in the literature. Then, a parametric study is carried out to investigate the effects of the operating temperat...
The shortage of fresh water is increasing at alarming levels, even in temperate zones. Desalination of brackish and seawater is being considered as a solution to growing water demands. Seawater desalination can be accomplished through reverse osmosis, vapor compression, vapor absorption and thermally driven systems. The most suitable process for large-scale desalination plants is successive brine evaporation-condensation with intermediate heat recovery. However, higher process efficiency could be achieved. The use of low temperature heat sources in seawater desalination requires an efficient method of evaporation and condensation at relatively low vapor pressures. This could be achieved in a humidification-dehumidification process using circulated air to enhance evaporation of water. Simulation models, however, neglect either the heat of evaporation or the changes in air and water properties within ...
In this study, adjoint sensitivity theory is developed for equations of two-dimensional steady-state flow in a confined aquifer. Both the primary flow equation and the adjoint sensitivity equation are solved using the Galerkin finite element method. The developed computer code is used to investigate the regional flow parameters of the Leadville Formation of the Paradox Basin in Utah and the Wolcamp carbonate/sandstone aquifer of the Palo Duro Basin in the Texas Panhandle. Two performance measures are evaluated, local heads and velocity in the vicinity of potential high-level nuclear waste repositories. The results illustrate the sensitivity of calculated local heads to the boundary conditions. Local velocity-related performance measures are more sensitive to hydraulic conductivities. The uncertainty in the performance measure is a function of the parameter sensitivity, parameter variance and the correlation between ...
The objective of this study is to develop an artificial neural network (ANN) model to predict the thermal conductivity of ethylene glycol-water solutions based on experimentally measured variables. The thermal conductivity of solutions at different concentrations and various temperatures was measured using the cylindrical cell method that physical properties of the solution are being determined fills the annular space between two concentric cylinders. During the experiment, heat flows in the radial direction outwards through the test liquid filled in the annual gap to cooling water. In the steady state, conduction inside the cell was described by the Fourier equation in cylindrical coordinates, with boundary conditions corresponding to heat transfer between the solution and cooling water. The performance of ANN was evaluated by a regression analysis between the predicted and the experimental values. The ...
A theory has been developed for the mechano-catalytic water-splitting, which is the system of simultaneous H{sub 2} and O{sub 2} evolution by stirring the powder of an oxide semiconductor in pure water under the condition that the stirring rod must be kept in contact with the surface of the glass vessel. The kinetic equations and the coupling strength of the frictional energy conversion between mechanical and electrical systems are calculated . The total system composed of the formation of the dangling bonds on the glass surface, the trapping of the semiconductor particles at the microcrevice of the glass surface, the strong field inside the fine particles due to the frictional electricity, the mechanism of charge transfer from the semiconductor to the stirring rod, the hopping conduction of positive hole, the electric current density injected into water from the semiconductors, and the tunnel chemical ...
Ependymal cells line the cerebral ventricles forming the interface that separates the cerebrospinal (CSF) and interstitial fluids (ISF). Extracellular molecules move between ependymal cells, whereas lipid soluble molecules pass both between and through cells. We measured the transfer of tritiated water (TOH) from CSF to blood across the ependymal and capillary interfaces by ventriculocisternal (VC) steady-state tissue clearance. Adult cats anesthetized with pentobarbital sodium underwent VC perfusion with the extracellular marker (/sup 14/C)sucrose and TOH added to the artificial CSF. Brain tissue was analyzed for depth of penetration of the isotopes into periventricular gray matter. We found that TOH distribution space was lower than expected from water content measurements, whereas sucrose space was normal. Using VC steady-state equations we calculated an ependymal permeability that was similar to the permeability of the ...
This document is a reference guide to the Xyce Parallel Electronic Simulator, and is a companion document to the Xyce Users Guide. The focus of this document is (to the extent possible) exhaustively list device parameters, solver options, parser options, and other usage details of Xyce. This document is not intended to be a tutorial. Users who are new to circuit simulation are better served by the Xyce Users Guide.
This document is a reference guide to the Xyce Parallel Electronic Simulator, and is a companion document to the Xyce Users' Guide. The focus of this document is (to the extent possible) exhaustively list device parameters, solver options, parser options, and other usage details of Xyce. This document is not intended to be a tutorial. Users who are new to circuit simulation are better served by the Xyce Users' Guide.
An update of the ODEtools Maple package, for the analytical solving of 1st and 2nd order ODEs using Lie group symmetry methods, is presented. The set of routines includes an ODE-solver and user-level commands realizing most of the relevant steps of the symmetry scheme. The package also includes commands for testing the returned results, and for classifying 1st and 2nd order ODEs.
The 3-dimensional (3D) extension C5G7 MOX benchmark problems were solved by CHAPLET-3D code which is based on the idea of dynamic linkage of the multi-plane method of characteristics solutions. The benchmark results are quite accurate in comparison with the reference solutions, independently of the axial solver incorporated in the CHAPLET-3D code. (author)
The UNSAT-H model was developed at Pacific Northwest National Laboratory (PNNL) to assess the water dynamics of arid sites and, in particular, estimate recharge fluxes for scenarios pertinent to waste disposal facilities. During the last 4 years, the UNSAT-H model received support from the Immobilized Waste Program (IWP) of the Hanford Site's River Protection Project. This program is designing and assessing the performance of on-site disposal facilities to receive radioactive wastes that are currently stored in single- and double-shell tanks at the Hanford Site (LMHC 1999). The IWP is interested in estimates of recharge rates for current conditions and long-term scenarios involving the vadose zone disposal of tank wastes. Simulation modeling with UNSAT-H is one of the methods being used to provide those estimates (e.g., Rockhold et al. 1995; Fayer et al. 1999). To achieve the above goals for assessing water dynamics and estimating ...
Condensing two-phase channel flow occurs in many industrial applications, including heating and refrigeration systems. It can also occur in certain nuclear reactor accidents. For example, during a small-break loss-of-coolant accident in a pressurized water reactor, following the partial depletion of the primary coolant, condensation of steam on the primary side of the steam generator tubes can provide a heat sink for disposal of the decay heat generated in the reactor core. Condensing two-phase flow can also play an important role in the operation of the passive emergency cooling system in the advanced simplified boiling water reactor. Here, steady-state condensation in the presence of a noncondensable in a concurrent two-phase channel flow is analyzed using a two-fluid model. The effect of noncondensables on the combined heat transfer at the liquid-gas mixture interphase is accounted for by using the stagnant film model, and closure relations ...
In this paper, a mathematical model for a counterflow wet cooling tower is derived, which is based on one-dimensional heat and mass balance equations using the measured heat transfer coefficient. The balance equations are solved numerically to predict the temperature change of air and water, as well as the humidity as a function of the cooling tower high. Experimental measurements on two pilot-scale cooling towers were carried out in order to analyze the performance of different cooling tower filling materials. Also, the performance of other cooling tower elements, such as droplet separators and water spray nozzles, was investigated in the pilot experiments. The flow distribution, i.e. the velocity field, upstream to the filling material was predicted using the three-dimensional version of the computational fluid dynamics (CFD) code Fluent/UNS, version 4.2. The calculated flow fields are presented for ...
A safety analysis code to design a pressurized water reactor and to obtain the licenses including entire proprietary rights is under development in domestic R and D project. The tasks of KAERI is to develop the constitutive relations including models for defining flow regimes and flow regime related models for inter-phase friction, wall frictions, wall heat transfer, and interphase heat and mass transfer in the two-phase three-field equations. In this paper, the process will be presented for choosing the best flow regime maps which occur in gas-liquid two-phase flow in horizontal and vertical tubes.
A safety analysis code to design a pressurized water reactor and to obtain the licenses including entire proprietary rights is under development in domestic R and D project. The tasks of KAERI is to develop the constitutive relations including models for defining flow regimes and flow regime related models for inter-phase friction, wall frictions, wall heat transfer, and interphase heat and mass transfer in the two-phase three-field equations. In this paper, the process will be presented for choosing the best flow regime maps which occur in gas-liquid two-phase flow in horizontal and vertical tubes.
Green River oil shale was air oxidized at subretorting temperatures. Off gases consisting of nitrogen, oxygen, carbon monoxide, carbon dioxide, and water were monitored and quantitatively determined. A mathematical model of the oxidation reactions based on a shrinking core model has been developed. This model incorporates the chemical reaction of oxygen and the organic material in the oil shale as well as the diffusivity of the oxygen into the shale particle. Diffusivity appears to be rate limiting for the oxidation. Arrhenius type equations, which include a term for oil shale grade, have been derived for both the chemical reaction and the diffusivity.
In aqueous solutions, dissolved ions interact strongly with the surrounding water, thereby modifying the solution properties in an ion-specific manner. These ion-hydration interactions can be accounted for theoretically on a mean-field level by including phenomenological terms in the free energy that correspond to the most dominant ion-specific interactions. Minimizing this free energy leads to modified Poisson-Boltzmann equations with appropriate boundary conditions. Here, we review how this strategy has been used to predict some of the ways ion-specific effects can modify the forces acting within and between charged interfaces immersed in salt solutions.
A safety analysis code, named as SPACE, for a pressurized water reactor is under development to obtain a licensing to be used for the PWR design and to hold entire proprietary rights. The task of KAERI is to develop the physical models and correlations which are required to solve the field equations. It can be divided into four parts; i) flow regime determination, ii) wall heat transfer, iii) wall and interfacial friction, iv) interfacial heat and mass transfer. This paper will describe the process to develop the models for the two-phase flow regime maps in the horizontal and vertical pipes.
A safety analysis code, named as SPACE, for a pressurized water reactor is under development to obtain a licensing to be used for the PWR design and to hold entire proprietary rights. The task of KAERI is to develop the physical models and correlations which are required to solve the field equations. It can be divided into four parts; i) flow regime determination, ii) wall heat transfer, iii) wall and interfacial friction, iv) interfacial heat and mass transfer. This paper will describe the process to develop the models for the two-phase flow regime maps in the horizontal and vertical pipes.
An optimal deployment pattern of flux mapping detectors for a Canada uranium-deuterium (CANDU)-600 pressurized heavy water reactor (PHWR) is determined by obtaining an optimal feedback relationship between flux measurements and zone controllers. The reactor core is modeled with a time-dependent two-group, two-dimensional diffusion equation, and flux perturbation are expressed by model expansions. The modal expansion coefficients are used as elements of the state vector representing the system dynamics. An optimal feedback matrix connecting the flux measurement vector to the control vector is derived by minimizing a quadratic performance index involving both the state and control vectors. We obtain the detector effectiveness in terms of the optimal feedback matrix and determine optimal detector locations for the Wolsung Unit 1 reactor in Korea. We have tested the methodology through evaluation of flux maps generated through the CANDU flux ...
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
Land snails move via adhesive locomotion. Through muscular contraction and expansion of their foot, they transmit waves of shear stress through a thin layer of mucus onto a solid substrate. Since a free surface cannot support shear stress, adhesive locomotion is not a viable propulsion mechanism for water snails that travel inverted beneath the free surface. Nevertheless, the motion of the freshwater snail, Sorbeoconcha physidae, is reminiscent of that of its terrestrial counterparts, being generated by the undulation of the snail foot that is separated from the free surface by a thin layer of mucus. Here, a lubrication model is used to describe the mucus flow in the limit of small amplitude interfacial deformations. By assuming the shape of the snail foot to be a traveling sine wave and the mucus to be Newtonian, an evolution equation for the interface shape is obtained and the resulting propulsive force on the snail is calculated. This ...
Purpose Nonpoint sources (NPS) pollution has been an important cause for water quality impairment worldwide. To take the temporal variations of both NPS pollution and in-stream attenuation into consideration, an inverse modeling approach and the load duration curve (LDC) method were combined for variable nutrient total maximum daily load (TMDL) development. Methods Water quality and hydrological parameters were monitored monthly along the ChangLe River system in 2004?2008. The catchment NPS export load (EL) and TMDL for total nitrogen (TN) were estimated by the inverse format of an existing stream nutrient transport equation. The LDC method was used to describe the variability of EL, TMDL, requiring load (RLR) and percent (the ratio between the RLR and the EL, RPR) reduction, and then to s...
Relative densities and relative massic heat capacities have been measured for aqueous solutions of Y(ClO_4)_3, Yb(ClO_4)_3, Dy(ClO_4)_3, and Sm(ClO_4)_3 at T=(288.15,298.15,313.15, and 328.15) K and p=0.1 MPa. These measurements were made in the concentration range 0.01624#<=#m/(mol#centre dot#kg"-"1)#<=#0.41822 using a Sodev 02D Vibrating Tube Densimeter and a Picker Microflow Calorimeter, respectively. To counter the potential effects of hydrolysis, aqueous solutions of the investigated salts were acidified with perchloric acid. After correcting for the presence of the acid, the measured properties were used to calculate apparent molar volumes and apparent molar heat capacities for solutions of the perchlorate salts in water. The calculated apparent molar properties were modeled at each investigated temperature using Pitzer ion interaction equations to produce estimates of apparent molar volumes and heat capacities at infinite dilution. ...
This thesis surveys the concepts of exergy and extends the exergy method of analysis from the standpoint of its applications to the power plant systems. After a brief historical review of exergy concepts, the general exergy equation is derived from the combined equation of First and Second Law, and it is shown that any special case of exergy equation is a simplified form of the general exergy equation. The mathematical method for the exergy analysis of a steady-state, steady-flow system, analogous to that of the First Law, is given. The exergy losses in a power plant are discussed. Then in order to examine these losses, the Second Law performance of major processes of combustion, compression, heat transfer, mixing and throttling have been analyzed analytically, and the exergy efficiencies are defined that accurately assess the thermodynamic performance of the corresponding processes. The methods for ...
A new N-body and hydrodynamical code, called RAMSES, is presented. It has been designed to study structure formation in the universe with high spatial resolution. The code is based on Adaptive Mesh Refinement (AMR) technique, with a tree based data structure allowing recursive grid refinements on a cell-by-cell basis. The N-body solver is very similar to the one developed for the ART code (Kravtsov et al. 97), with minor differences in the exact implementation. The hydrodynamical solver is based on a second-order Godunov method, a modern shock-capturing scheme known to compute accurately the thermal history of the fluid component. The accuracy of the code is carefully estimated using various test cases, from pure gas dynamical tests to cosmological ones. The specific refinement strategy used in cosmological simulations is described, and potential spurious effects associated to shock waves propagation in the resulting AMR grid are discussed and ...
Notched and smooth cylindrical plasma nitrided (PN) and quench and tempered (Q and T) steel specimens made of a B-Mn SS2131 ({approx}AISI 15B21H) steel have been exposed to constant amplitude plane reversed bending corrosion fatigue tests (R = -1) at 47 Hz in sea water. S - N curves show that sea water suppresses the fatigue limit and reduces fatigue strength (especially at long lives) of smooth and notched Q and T and PN specimens. Plasma nitriding improves the corrosion fatigue resistance of Q and T specimens; this is associated with the good corrosion resistance of {epsilon} and {gamma}`-phases, the enhancement of corrosion and fatigue by compressive residual stresses, and the consumption of H{sup +} ions during reduction of nitrogen. This improvement is more significant for smooth specimens and for long lives. Notch sensitivity of Q and T and PN specimens decreases with fatigue life. Pitting corrosion, cyclic applied stress and residual ...
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 ...
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 ...
The quaternionic formalism for subluminal field equations (Maxwell's equations) and its interrelationship with complex superluminal Lorentz transformations have been given and it has been shown that the quaternionic forms of relativistic equations describe tachyons.
In this paper, the authors first describe a fourth order accurate finite difference discretization for both the Laplace equation and the heat equation with Dirichlet boundary conditions on irregular domains. In the case of the heat equation, they use an i...
The analytical results of blowdown characteristics and its thrust force were compared with the experiment, which were performed as pipe whip tests under the PWR LOCA conditions on the hypothetical accident of guillotine break of pipes. The blowdown thrust force was obtained by the integral momentum equation about single-phase flow, homogeneous and separated two-phase flow, assuming critical pressure at the exit if critical flow condition was satisfied. The following results are obtained: (1) The node-junction method is useful for the analysis of water hammer phenomena and of the blowdown thrust force. (2) The Henry-Fauske model for subcooled critical flow is effective for the analysis of the maximum thrust force under the PWR LOCA conditions. The jet thrust parameter of analysis and experiment is 1.08. (3) The thrust parameter of saturated blowdown has the same one with the value under pressurized condition when the stagnant pressure is chosen ...
The analytical results of blowdown characteristics and thrust forces were compared with the experiments, which were performed as pipe whip and jet discharge tests under the PWR LOCA conditions. The blowdown thrust forces were obtained by Navier-Stokes momentum equation for a single-phase, homogeneous and separated two-phase flow, assuming critical pressure at the exit if a critical flow condition was satisfied. The following results are obtained: (1) The node-junction method is useful for both the analyses of the blowdown thrust force and of the water hammer phenomena. (2) The Henry-Fauske model for subcooled critical flow is effective for the analysis of the maximum thrust force under the PWR LOCA conditions. The jet thrust parameter of the analysis and experiment is equal to 1.08. (3) The thrust parameter of saturated blowdown has the same one with the value under pressurized condition when the stagnant pressure is chosen as the saturated ...
The 3-D model program was developed to apply a mise-a-la-masse method to linear current sources with optional shapes. Mise-a-la-masse method is a bipolar mapping survey method using current sources installed in wells, and recently has been widely used for direct exploration of geothermal reservoirs. This method is also used for monitoring underground fluid as electrode arrangement of fluid flow tomography for surveying underground fluid (geothermal fluid, underground water, petroleum). In the geothermal reservoir exploration, the casing pipes of wells are used as linear current sources, and measured data are processed as those based on vertical current source. In the largely inclined well, the inclination of current sources should be considered. The 3-D modeling program was developed by difference calculus using the theoretical potential equation and apparent resistivity based on linear current sources with optional 3-D shapes. The numerical ...
The main goal of this paper is to present a methodology for calculating the radioactivity in the moderator and heat transport systems of Cernavoda NPP Unit 1, with the intention to improve the knowledge on the radionuclides inventories in the operational waste streams, and to aid the licensing process of new near surface repository. In the present paper we describe our methodology for estimating H-3 and C-14 production rates in the heavy-water moderator and heat transport systems using the capacity factors from 1997 to 2007 years. The radioactivity of the difficult-to-measure nuclides is predicted by scaling method using measured concentration in reference CANDU 6 reactor Gentilly-2. The difficult-to-measure radionuclides of primary interest in this study were those with long half-lives which have a significant role for post-closure safety assessment. The equation used to scale fission products (parents and daughters) is based on the ...
Korea Advanced Energy Research Institute (KAERI) conceived and started the development of MARS code with the main objective of producing a state-of-the-art realistic thermal hydraulic systems analysis code with multi-dimensional analysis capability. MARS achieves this objective by very tightly integrating the one dimensional RELAP5/MOD3 with the multi-dimensional COBRA-TF codes. The method of integration of the two codes is based on the dynamic link library techniques, and the system pressure equation matrices of both codes are implicitly integrated and solved simultaneously. In addition, the Equation-Of-State (EOS) for the light water was unified by replacing the EOS of COBRA-TF by that of the RELAP5. This models and correlations manual provides a complete list of detailed information of the thermal-hydraulic models used in MARS, so that this report would be very useful for the code users. The overall structure of the ...
Recent modifications to the EQ3/6 geochemical modeling software package provide for the use of Pitzer's equations to calculate the activity coefficients of aqueous species and the activity of water. These changes extend the range of solute concentrations over which the codes can be used to dependably calculate equilibria in geochemical systems, and permit the inclusion of ion pairs, complexes, and undissociated acids and bases as explicit component species in the Pitzer model. Comparisons of calculations made by the EQ3NR and EQ6 compuer codes with experimental data confirm that the modifications not only allow the codes to accurately evaluate activity coefficients in concentrated solutions, but also permit prediction of solubility limits of evaporite minerals in brines at 25/sup 0/C and elevated temperatures. Calculations for a few salts can be made at temperatures up to approx. 300/sup 0/C, but the temperature range for most ...
The super simulator 'SAMPSON' has been developed to show that there exist certain safety margins for light water reactors under hypothetical severe accidents and to investigate realistic measures of accident management by simulating accidents with a parallel computer. Heat-up of fuel rods and release of fission products from fuels are important factors to evaluate source terms. Models for fuel rod heat-up, hydrogen production due to cladding oxidation and cladding deformation and failure in the core region have been developed in the fuel rod heat-up analysis module. Fuel temperatures were calculated by solving the heat conduction equation. The calculated results for fuel temperature and hydrogen production were compared with CORA-13 experiment results. The comparisons showed prediction capability for the heat-up of fuel rods. The fission product release analysis module incorporates with models for fission product transport within fuel pellets, ...
The present paper is concerned with development and application of a so-called Effective Convection Model (ECM), which aims to provide a detailed, mechanistic description of heat transfer processes in a BWR lower plenum. The ECM is a Computational Fluid Dynamics (CFD)-like tool which employs a simpler and more effective approach to compute heat transfer by solving only energy conservation equation instead of solving the full set of Navier-Stokes and energy equations by a CFD code. We implement the ECM in a CFD code (Fluent), with detailed description of the ECM development, implementation and validation. A dual approach is used to validate the ECM, namely validation against experimental data and against heat transfer results obtained by CFD predictions in the same geometries and conditions. Insights gained from CFD simulations are also used to improve ECM. The ECM capability as an effective tool to simulate heat transfer of an internally heated ...
The present paper is concerned with development and application of a so-called Effective Convection Model (ECM), which aims to provide a detailed, mechanistic description of heat transfer processes in a BWR lower plenum. The ECM is a Computational Fluid Dynamics (CFD)-like tool which employs a simpler and more effective approach to compute heat transfer by solving only energy conservation equation instead of solving the full set of Navier-Stokes and energy equations by a CFD code. We implement the ECM in a CFD code (Fluent), with detailed description of the ECM development, implementation and validation. A dual approach is used to validate the ECM, namely validation against experimental data and against heat transfer results obtained by CFD predictions in the same geometries and conditions. Insights gained from CFD simulations are also used to improve ECM. The ECM capability as an effective tool to simulate heat transfer of an internally heated ...
A new scheme to accurately retain kinetic electron effects in particle-in-cell (PIC) simulations for the case of electrostatic drift waves is presented. The splitting scheme, which is based on exact separation between adiabatic and on adiabatic electron responses, is shown to yield more accurate linear growth rates than the standard df scheme. The linear and nonlinear elliptic problems that arise in the splitting scheme are solved using a multi-grid solver. The multi-grid particle-in-cell approach offers an attractive path, both from the physics and numerical points of view, to simulate kinetic electron dynamics in global toroidal plasmas.
A set of Maple V R.3/4 computer algebra routines for the analytical solving of 1st. order ODEs, using Lie group symmetry methods, is presented. The set of commands includes a 1st. order ODE-solver and routines for, among other things: the explicit determination of the coefficients of the infinitesimal symmetry generator; the construction of the most general invariant 1st. order ODE under given symmetries; the determination of the canonical coordinates of the underlying invariant group; and the testing of the returned results.
The chemical thermodynamic concepts used in the calculation of solubility data of silica (quartz) are presented taking into account the PVT characteristics of water. The temperature-dependence trends between the thermodynamically calculated and the experimental quartz solubility data are very similar, but the values are widely different at high temperatures. The experimental solubility, especially along the saturation curve at high temperature and thermodynamic data for silica need to be reevaluated in order to use silica chemistry to understand geological processes. There could exist a wide range of values for silica solubility at a specified temperature, depending upon the amount of water in the reaction vessel. Thus the silica contents in geothermal fluid, in general, cannot be used as a geothermometer to estimate the reservoir temperature. The derivation of a silica geothermometer needs an extra assumption about the total amount of ...
In summary, a scaling analysis of a water-cooled Reactor Cavity Cooling System (RCCS) system was performed based on generic information on the RCCS design of PBMR. The analysis demonstrates that the water-cooled RCCS can be simulated at the ANL NSTF facility at a prototypic scale in the lateral direction and about half scale in the vertical direction. Because, by necessity, the scaling is based on a number of approximations, and because no analytical information is available on the performance of a reference water-cooled RCCS, the scaling analysis presented here needs to be 'validated' by analysis of the steady state and transient performance of a reference water-cooled RCCS design. The analysis of the RCCS performance by CFD and system codes presents a number of challenges including: strong 3-D effects in the cavity and the RCCS tubes; simulation of turbulence in flows characterized by natural ...
This thesis was carried out at Forsmark Nuclear Power Plant. The power plant in Forsmark consists of three boiling water reactors (BWR) which produce about 17% of Swedish electricity. In a BWR the nuclear reactions are used to boil water inside the reactor vessel. The water works both as a coolant and as a moderator and the resulting steam is used directly to run the turbines. A problem when running a BWR at low flow conditions is the density wave oscillations that might occur to the water flow inside the fuel assemblies. These oscillations arise due to the connection between power and flow rate in a heated channel with two-phase flow. In order to improve the stability performance of the channel an orifice plate is placed at the inlet of each fuel assembly. Today these orifice plates have sharp edges and a constant resistance coefficient. Experimental work has been done with progressive orifices, the ...
Soil erodibility (the K factor in the Universal Soil Loss Equation, USLE) is an important index to measure soil susceptibility to water erosion, and an essential parameter needed for soil erosion prediction. To evaluate the appropriateness of the nomograph and other methods for estimating the K factor for the USLE and to develop a relationship for soil erodibility estimation for Chinese soils, a set of soil erodibility values was calculated using soil loss data from natural runoff plots at 13 sites in eastern China. The definition of soil erodibility in relation to the USLE was strictly followed. Comparing these measured values to those estimated using the nomograph method, the method adopted for the EPIC model and the formula of Shirazi and Boersma, we found that all these estimated soil ...
Aqueous radioactive high-level waste slurries are combined during processing steps that ultimately produce a stable borosilicate glass waste form. Chemically treated waste slurries are combined with each other and with glass frit-water slurries to produce the melter feed. Understanding the evolution of the rheological properties of the slurries is an important aspect of removing and treating the stored waste. To a first approximation, combinations of colloidal waste slurry with {approx}0.1-mm mean diameter glass frit or glass beads act in an analogous matter to slurries of spherical beads in Newtonian liquids. The non-Newtonian rheological properties of the waste slurries without frit, however, add complexity to the hydrodynamic analysis. The use of shear rate dependent apparent viscosities with the modified Einstein equation was used to model the rheological properties of aqueous frit-waste slurries.
This report describes an investigation of the performance of air-augmented waterjet thrusters. The investigation is divided into two parts. Part one describes a study of the injection of large diameter air bubbles to reduce heat transfer between the air and water and, thus, increase the thrust. This study is done analytically by solving the appropriate, governing equations. The results indicate that the approach is very worth-while since thrust can be increased ten percent as compared to air-augmentation with small bubbles. Several practical systems are discussed for injecting large air bubbles. Part two is an experimental investigation of the thrust produced under a wind range of conditions. The results show that a thrust increase of more than 20 percent can be achieved by air-augmentation. Further, the performance improves for larger thrust nozzles.
By assuming that not only counter-ions but DNA molecules as well are thermally distributed according to a Boltzmann law, we propose a modified Poisson-Boltzmann equation, at the classical level, as a starting point to compute the effects of quantum fluctuations of the electric field on the interaction among DNA-cation complexes. The latter are modeled here as infinite one-dimensional wires (?-functions). Our goal is to single out such quantum-vacuum-driven interaction from the counterion-induced and water-related interactions. We obtain a universal, frustration-free Casimir-like (codimension 2) interaction that extensive numerical analysis show to be a good candidate to explain the formation and stability of DNA aggregates. Such Casimir energy is computed for a variety of configurations of...
Despite the proven utility of the Universal Soil Loss Equation (USLE) to predict soil rainfall erosion, it has limited precision in predicting soil loss in some soil series. Also, runoff is a factor in rainfall erosion and some models, such as the modified USLE (USLE-M), include runoff in predicting soil erosion. The soil erodibility factor is an essential component in both models and this may be affected by application of soil conditioners, such as polyacrylamide (PAM). The objectives of this study were to determine the effects of different rates of PAM (0, 1, 2, 4, and 6kgha-1) applied with sprinkler irrigation water, followed by two sprinkler irrigations with no PAM, on the erodibility of a loam soil under a rainfall simulator in a laboratory. The range of erodibility factor Ku obtained...
A one-dimensional tachyon Klein-Gordon equation is reduced to a nonrelativistic-tachyon equation of motion. The interpretation of this reduced equation leads to the following conclusions: 1) tachyons can be localized in time instead of in space as compared with bradyons, 2) space representation and momentum representation of bradyonic quantum equation of motion are replaced by time representation and energy representation in tachyon quantum equation of motion and 3) with the aid of these results, it has been found that the solutions of the tachyon Klein-Gordon equation of motion form a complete set. (author).
The densities and refractive indices of two binary systems (water + MEA and water + MDEA) and three ternary systems (water + MEA + CO{sub 2}, water + MDEA + CO{sub 2}, and water + MEA + MDEA) used for carbon dioxide (CO{sub 2}) capture were measured over the range of compositions of the aqueous alkanolamine(s) used for CO{sub 2} absorption at temperatures from 295 to 338 K. Experimental densities were modeled empirically, while the experimental refractive indices were modeled using well-established models from the known values of their pure-component densities and refractive indices. The density and Gladstone-Dale refractive index models were then used to obtain the compositions of unknown samples of the binary and ternary systems by simultaneous solution of the density and refractive index equations. The results from this technique have been compared with HPLC ...
The present situation of the researches that have been made for predicting the process and range of diffusion of warm drain is reviewed. This review is divided into eight sections. The first section deals with the present situation of warm drain from power plants. For the establishment of drainage standard, there are many difficult problems to be solved because water temperature differs in its nature from other regulation items. In the second section, the process of diffusion and cooling of warm drain is explained. The third section deals with the diffusion characteristics of warm drain in Japanese coastal sea due to water temperature. Two types of diffusion are known. One is dominant irregular current, and the other is periodical reciprocating stream. The fourth section deals with the methods of prediction of diffusion. Research methods and simulation models are described. The fifth section deals with the prediction with single layer model. ...
In order to study the boundary conditions for the occurrence of ratcheting, a testing rig was erected. In this rig, small-scale mock-ups of the first wall were tested under mechanical conditions similar to those expected during disruptions. These mock-ups were made of stainless steel AISI 316. They were heated on the surface and water cooled from inside. Disruption forces were simulated by a high-speed hydraulic bending system by impact loads. The mock-ups were heated by a high-frequency generator and cooled by water through two cooling channels.In addition to the experimental work, simplified analytical analysis and numerical calculations were carried out using the finite element program ABAQUS. In these calculations, the temperature and stress fields in the mock-ups were investigated and the influence of the different constitutive equations was studied.The experimental results confirm the ratcheting behaviour predicted by ...
The design of high level radioactive waste (HLW) repositories in deep geological media in which bentonite clay is proposed as a sealing material leads to the need of further studying the behaviour of highly compacted expansive soils when subjected to mechanical, hydraulic and thermal changes. Laboratory tests may help to understand the processes that take place in the clay barrier under simple and controlled conditions and to develop the governing equations. The laboratory tests enable to isolate the different processes, making their interpretation easier, and provide with fundamental data concerning the parameters to be used in the models. The extremely low permeability of these materials, their avidity for water (high suction) and their high swelling capacity make necessary the modification of the conventional laboratory techniques and procedures to determine basic physical parameters. The main hydraulic properties of the barrier to be ...
Full text of publication follows: A porous media solution PORFLO has been developed for the 3-dimensional two-phase flow by describing the process facility in Cartesian or cylindrical coordinates. The local porosity fraction is applied for distinguishing the fluid filled volumes from the solid structures. The solid structure contribute the two-phase flow through the wall friction, flow area and heat transfer. Optionally the solid structure may contain primary liquid of steam generators, steam in the higher temperature and pressure to be condensed or electrical heating power. By using these optional boundary conditions three different process facilities have been analysed. The thermohydraulic solution based on 5-equation approach, where the conservation equations are solved for the liquid and gas (vapour) mass, mixture momentum (giving the velocity only for the mixture), liquid and gas energy, is described shortly. In addition to that the ...
In order to practice design-by-analysis of thermohydraulics design of BWR fuel rod bundles, the subchannel analysis would play a major role. There, one of the immediate concerns is improvement in its predictive capability of boiling transition phenomena on the fuel rod surface. This capability strongly depends on the modeling of thermohydraulics phenomena of interests: 1) vapor-liquid redistribution by inter-subchannel exchanges due to the diversion cross flow, turbulent mixing and void drift, 2) liquid film behaviors, 3) transition of two-phase flow regimes, 4) droplet entrainment and deposition and 5) spacer-droplet interactions. These are considered to be five key factors in understanding the BT in BWR fuel rod bundles. This paper describes a progress and current status in the second year of the three year project on developing generalized boiling transition models with the above five key factors being focused on. A combined approach of experiment and computation is described in ...
The purpose of this study is to investigate the feasibility of visualization and void fraction measurement of air-water two-phase flow in a small diameter tube (I.D.: 4.08 mm) by using the real-time neutron radiography and image processing techniques. Video images of two-phase flow were taken by using the real-time neutron radiography system (thermal neutron radiography facility No.2) installed at the Japan Research Reactor 3M of the Japan Atomic Energy Research Institute. The shape of bubbles and its moving behavior were clearly observed from the video images. The image corrections for dark current, shading, field intensity fluctuation and electrical system drift were examined in order to measure the void fraction from the video images. Though, generally speaking, the effect of the scattered neutron could not be ignored for quantification of the images taken by the neutron radiography, the scattered neutron could not affect the final results of void fraction in ...
The purpose of this study is to investigate the feasibility of visualization and void fraction measurement of air-water two-phase flow in a small diameter tube (inner diameter; 4.08mm) by using the real-time neutron radiography and image processing techniques. Video images of two-phase flow were taken by using the real-time neutron radiography system (thermal neutron radiography facility No.2) installed at the Japan Research Reactor 3 M of the Japan Atomic Energy Research Institute. The shape of bubbles and its moving behavior were clearly observed from the video images. The image corrections for dark current, shading, field intensity fluctuation and electrical system drift were examined in order to measure the void fraction from the video images. Though, generally speaking, the effect of the scattered neutron could not be ignored for quantification of the images taken by the neutron radiography, the scattered neutron could not affect the final results of void ...
The purpose of this study is to investigate the feasibility of visualization and void fraction measurement of air-water two-phase flow in a small diameter tube (I.D.: 4.08 mm) by using the real-time neutron radiography and image processing techniques. Video images of two-phase flow were taken by using the real-time neutron radiography system (thermal neutron radiography facility No.2) installed at the Japan Research Reactor 3M of the Japan Atomic Energy Research Institute. The shape of bubbles and its moving behavior were clearly observed from the video images. The image corrections for dark current, shading, field intensity fluctuation and electrical system drift were examined in order to measure the void fraction from the video images. Though, generally speaking, the effect of the scattered neutron could not be ignored for quantification of the images taken by the neutron radiography, the scattered neutron could not affect the final results of void fraction in ...
According to the present concept, the low and intermediate level wastes generated during the Cernavoda NPP operation will be disposed in a near surface repository. The Saligny site, placed in the NPP protected area, has been proposed for their disposal. Geologically, the main components of this site are the quaternary loess, the Precambrian and Pre-quaternary clays, the Eocene and Barremian limestone. Hydrologically, the site can be divided into a vadose zone down to 45-50 m and three distinct aquifers, two of them in the limestone beds and the third in the lenses of sand and limestone existing in the pre-quaternary clay layer. A large research program for site characterization was initiated in 1996. At present, the site characteristics requested for safety analysis have been experimentally measured on soil samples or calculated by different computer programs. Hundreds of experimental values of the density, porosity, hydraulic conductivity, soil-water retention, ...
Research highlights: ? We model power oscillations in boiling water reactors using a lumped parameter model. ? The nature and amplitudes of oscillations is obtained using a nonlinear analysis. ? The method of multiple scales has been used for the analytical treatment. ? Fuel temperature coefficient of reactivity determines the nature of oscillations. ? The presented systematic method of analysis useful for reduced order reactor models. - Abstract: In this paper, we perform a parametric study of the nonlinear dynamics of a reduced order model for boiling water reactors (BWR) near the Hopf bifurcation point using the method of multiple scales (MMS). Analysis has been performed for general values of the parameters, but the results are demonstrated for parameter values of the model corresponding to the advanced heavy water reactor (AHWR). The neutronics of the AHWR is modeled using point reactor kinetic ...
Several steam generator tube rupture accidents have occurred at plants in the past. In this paper the Computational Multi-Fluid Dynamics (CMFD) investigation of the horizontal steam generator thermal-hydraulics during the tube rupture accident is performed. A guillotine of a steam generator U-tube is assumed with choked flow from the primary to the secondary side of the steam generator. We have computed water and steam velocity fields, steam volume fraction distribution on the steam generator secondary (shell) side, as well as the swell level increase. The simulation results are a support to the safety analyses of the steam generator tube rupture accident. Numerical simulation is performed with the multidimensional multi-fluid modelling approach. The two-phase flow around steam generator tubes in the bundle is modelled by the porous media approach. Interfacial mass, momentum and energy transfer are modelled with the closure laws. The governing ...
We obtain the radial Dirac equations for the Plebanski metric and transform these equations into one-dimensional wave equations. Finally we conclude with the result that the electron and the neutrino fields are not superradiant.
In the procedure of the hydraulic control system analysis, a linearized approximate equation described by the first order term of Taylor's series has been widely used. Such a linearized equation is effective just near the operating point. In this study, the authors estimate computational errors in the process of applying the existing linearized equation stated above. For evaluating the computational accuracy in practical applications of the linearized equations, dynamic behaviors of hydraulic control systems are investigated through simulations with several kinds of representative hydraulic systems and the linearized equations suggested in this study.
This paper presents an innovative way in which university education can help pre-service teachers become better problem-solvers. The central idea is to use the ?Technology Fair?? as a means for promoting pre-service teachers pedagogical content knowledge about technological problem solving skills. This innovation is supported with results from a study carried out in autumn 2004. The purpose of the study was to investigate the influence of a procedure of working with primary school children to complete and present a technology fair project, on the educational value and meanings attached to problem solving skills by pre-service teachers. Pre-tests, mid-test and post-tests were administered to the pre-service teachers before, during, and after the preparation of the technology fair, respectiv...
Lean premixed combustion (LPC) of natural gas is of considerable interest in land-based gas turbines for power generation. However, modeling such combustors and adequately addressing the concerns of LPC, which include emissions of nitrogen oxides, carbon monoxide and unburned hydrocarbons, remains a significant challenge. In this paper, characteristics of published simulations of gas turbine combustion are summarised and methods of modeling turbulent combustion are reviewed. The velocity-composition PDF method is selected for implementation in a new comprehensive model that uses an unstructured-grid flow solver. Reduced mechanisms for methane combustion are evaluated in a partially stirred reactor model. Comprehensive model predictions of swirl-stabilised LPC of natural gas are compared with detailed measurements obtained in a laboratory-scale combustor. The model is also applied to industrial combustor geometries. (Author)
Ramjet burners are known to produce highly unsteady operating conditions with strong couplings between combustion, acoustics and flow dynamics. Predicting such operating limit-cycles still remains a difficult task for Computational Fluid Dynamics (CFD) although recent use of Large Eddy Simulation (LES) clearly opens new possibilities. The main difficulties for LES are to properly address numerically specific flow features at the same time. For example, a proper representation of the acoustic ramjet eigenmodes necessitates for the solver to be able to treat shocks often present at the inflow conditions without interfering with the low Mach number flow in the region of combustion. Chemistry modelling is another difficulty and it is still not clear what level of description is sufficient to r...
This is the first of three planned papers describing ZAP, a satisfiability engine that substantially generalizes existing tools while retaining the performance characteristics of modern high-performance solvers. The fundamental idea underlying ZAP is that many problems passed to such engines contain rich internal structure that is obscured by the Boolean representation used; our goal is to define a representation in which this structure is apparent and can easily be exploited to improve computational performance. This paper is a survey of the work underlying ZAP, and discusses previous attempts to improve the performance of the Davis-Putnam-Logemann-Loveland algorithm by exploiting the structure of the problem being solved. We examine existing ideas including extensions of the Boolean language to allow cardinality constraints, pseudo-Boolean representations, symmetry, and a limited form of quantification. While this paper is intended as a survey, our research ...
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.
To investigate the 2D/1D fusion core analysis method, a 1D neutron transport problem solver, PEACH-ID, is developed. It is a code of method of characteristics (MOC), both the usual fiat-source step characteristics (SC) scheme and linear source (LS) approximation scheme are adopted for tracking calculation along the neutron flying trajectory. Exponential function interpolation table and fission source extrapolation are adopted as two major methods to accelerate the computational process. Numerical results demonstrate that PEACH-1D is accurate and efficient, and the proposed LS scheme is able to handle quite larger mesh division and deserves much more application in the MOC codes. (authors)
A method for calculating flow through an airfoil cascade drawn on a surface of revolution is discussed. The three-dimensional flow was assumed to be represented by part-channels of varying width. The basic equations are the equation of continuity and the ...
The Yang-Mills equations are formulated in the form of generalized Maurer-Cartan equations, such that the corresponding algebraic operations are shown to satisfy the defining relations of homotopy Lie superalgebra.
In order to practice a design-by-analysis of thermohydraulics design of BWR fuel rod bundles, the subchannel analysis would play a major role. There, the immediate concern is improvement in its predictive capability of CHF due in particular to the film dryout (boiling transition phenomena: BT) on the fuel rod surface. Constitutive equations in the subchannel analysis formulation are responsible for the quality of calculated results. The constitutive equations are a result of integration of the local and instantaneous description of two-phase flows over the subchannel control volume. In general, they are expressed in terms of subchannel-control-volume- as well as area-averaged two-phase flow state variables. In principle the information on local and instantaneous physical phenomena taking place inside subchannels must be counted for in the algebraic form of the equations on the basis of a more mechanistic modeling approach. ...
Full text of publication follows:In order to practice a design-by-analysis of thermohydraulics design of BWR fuel rod bundles, the subchannel analysis would play a major role. There, the immediate concern is improvement in its predictive capability of CHF due in particular to the film dryout (boiling transition phenomena: BT) on the fuel rod surface. Constitutive equations in the subchannel analysis formulation are responsible for the quality of calculated results. The constitutive equations are a result of integration of the local and instantaneous description of two-phase flows over the subchannel control volume. In general, they are expressed in terms of subchannel-control-volume- as well as area-averaged two-phase flow state variables. In principle the information on local and instantaneous physical phenomena taking place inside subchannels must be counted for in the algebraic form of the equations on the basis of a ...
In order to practice a design-by-analysis of thermohydraulics design of BWR fuel rod bundles, the subchannel analysis would play a major role. There, the immediate concern is improvement in its predictive capability of CHF due in particular to the film dryout (boiling transition phenomena: BT) on the fuel rod surface. Constitutive equations in the subchannel analysis formulation are responsible for the quality of calculated results. The constitutive equations are a result of integration of the local and instantaneous description of two-phase flows over the subchannel control volume. In general, they are expressed in terms of subchannel-control-volume- as well as area-averaged two-phase flow state variables. In principle the information on local and instantaneous physical phenomena taking place inside subchannels must be counted for in the algebraic form of the equations on the basis of a more mechanistic modeling approach. ...
Clean Water Act You are here: Water Laws & Regulations Policy & Guidance Clean Water Act Water Quality and 401 Certification Water Quality and 401 Certification Water Quality...
This investigation presents scale-up of a detailed heterogeneous geostatistical model to a full field reservoir simulation model, considering both single and two-phase flow properties. The model represents a typical low permeability Danish North Sea chalk reservoir and includes capillary pressure and saturation end-point variations. Two new up-scaling methods has been investigated, all based on fine scale simulation on a cross section of the geomodel. The first methods assumes piston style behaviour and a coupled viscosity is introduced into the basic Darcy`s equations. The second method is a modification of the JBN method traditionally applied in analysing results from core flooding experiments, which emerged as the most successful and therefore also the recommended method. 1. In addition to the up scaling work we review the Equivalent Radius Method for capillary pressure normalisation with explicit derivation of type functions for Maastrichtian and Danian chalk ...
Carbon dioxide is receiving renewed interest as an efficient and environmentally safe refrigerant in a number of applications, including mobile air conditioning and heat pump systems, and hot water heat pumps. Compact heat exchangers for CO{sub 2} systems are designed with small-diameter tubing. The purpose of this study is therefore to provide a better basis for understanding and predicting heat transfer and pressure drop during flow vaporization of CO{sub 2} in microchannels. The ''unusual'' properties of carbon dioxide give heat transfer and two-phase flow characteristics that are very different from those of conventional refrigerants. Examples of these differences are the much higher pressure, the resulting high vapour density, a very low surface tension, and a low liquid viscosity. High pressure and low surface tension has a major effect on nucleate boiling characteristics, and earlier test data have shown a clear dominance ...
High-energy ionizing radiation is a prominent modality for the treatment of many cancers. The approaches to electron dose calculation can be categorized into semi-empirical models (e.g. Fermi-Eyges, convolution-superposition) and probabilistic methods (e.g. Monte Carlo). A third approach to dose calculation has only recently attracted attention in the medical physics community. This approach is based on the deterministic kinetic equations of radiative transfer. We derive a macroscopic partial differential equation model for electron transport in tissue. This model involves an angular closure in the phase space. It is exact for the free streaming and the isotropic regime. We solve it numerically by a newly developed HLLC scheme based on Berthon et al (2007 J. Sci. Comput. 31 347-89) that exactly preserves the key properties of the analytical solution on the discrete level. We discuss several test cases taken from the medical physics literature. ...
This manual describes the use of the Xyce Parallel Electronic Simulator. Xyce has been designed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: (1) Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). Note that this includes support for most popular parallel and serial computers. (2) Improved performance for all numerical kernels (e.g., time integrator, nonlinear and linear solvers) through state-of-the-art algorithms and novel techniques. (3) Device models which are specifically tailored to meet Sandia's needs, including some radiation-aware devices (for Sandia users only). (4) Object-oriented code design and implementation using ...
This manual describes the use of the Xyce Parallel Electronic Simulator. Xyce has been designed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: (1) Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). Note that this includes support for most popular parallel and serial computers; (2) Improved performance for all numerical kernels (e.g., time integrator, nonlinear and linear solvers) through state-of-the-art algorithms and novel techniques. (3) Device models which are specifically tailored to meet Sandia's needs, including some radiation-aware devices (for Sandia users only); and (4) Object-oriented code design and implementation using ...
In order to predict or to evaluate the structural reliability of a variety of nuclear components under actual changing operating conditions, it is necessary to take the detailed geometrical features of components into consideration on one hand, and thermo-mechanical loading conditions with space and time variations during the lifetime, on the other hand. This consideration leads the authors to develop a super large-scale structural analysis program based on a parallel computational FEM with a variable scale of PC cluster architecture, whose prototype system is presented in a companion paper in this Symposium. Development of the PC cluster technology is one of the most important subjects in large-scale structural simulation, reducing computer costs and thus making simulation more acceptable to the engineering community. This paper presents the developmental status of the parallel FEM solver, demonstrating computational capabilities with three dimensional effects ...
A general theorem on the GBDT version of the B\\"acklund-Darboux transformation for systems rationally depending on the spectral parameter is treated and its applications to nonlinear equations are given. Explicit solutions of direct and inverse problems for Dirac-type systems, including systems with singularities, and for the system auxiliary to the $N$-wave equation are reviewed. New results on explicit construction of the wave functions for radial Dirac equation are obtained.
In this paper, we study the existence and uniqueness of periodic solutions of the differential equation of the form . Here, we obtain some sufficient conditions which guarantee the existence of periodic solutions. This equation is a quite general third- order nonlinear vector differential equation, and one example is given for illustration of the subject.
The hierarchy of integrable nonlinear equations associated with the quadratic bundle is considered. The expressions for the solution of linearization of these equations and their conservation law in the terms of solutions of corresponding Lax pairs are found. It is shown for the first member of the hierarchy that the conservation law is connected with the solution of linearized equation due to the Noether's theorem. The local hierarchy and three nonlocal ones of the infinitesimal symmetries and conservation laws explicitly expressed through the variables of the nonlinear equations are derived. (author)
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.
The time-dependent probability distribution of neutrons, precursors and detectrons in a space-independent, low-power, multiplying assembly with a source is developed in the multigroup approximation as forward and backward Kolmogorov equations. The relationship between these as adjoint equations is made explicit in a tensor notation and the equations developed in the generating-function formalism. (author).
As a result of the large limestone deposits available in Poland, the low cost of reagent acquisition for the largescale technological use and relatively well-documented processes of flue gas desulfurization (FGD) technologies based on limestone sorbent slurry, wet scrubbing desulfurization is a method of choice in Poland for flue gas treatment in energy production facilities, including power plants and industrial systems. The efficiency of FGD using the above method depends on several technological and kinetic parameters, particularly on the pH value of the sorbent (i.e., ground limestone suspended in water). Consequently, many studies in Poland and abroad address the impact of various parameters on the pH value of the sorbent suspension, such as the average diameter of sorbent particles (related to the limestone pulverization degree), sorbent quality (in terms of pure calcium carbonate [CaCO3] content of the sorbent material), stoichiometric surfeit of CaCO3 in ...
This document is a user`s manual for the Rn3D finite element code. Rn3D was developed to simulate gas flow and radon transport in variably saturated, nonisothermal porous media. The Rn3D model is applicable to a wide range of problems involving radon transport in soil because it can simulate either steady-state or transient flow and transport in one-, two- or three-dimensions (including radially symmetric two-dimensional problems). The porous materials may be heterogeneous and anisotropic. This manual describes all pertinent mathematics related to the governing, boundary, and constitutive equations of the model, as well as the development of the finite element equations used in the code. Instructions are given for constructing Rn3D input files and executing the code, as well as a description of all output files generated by the code. Five verification problems are given that test various aspects of code operation, complete with example input ...
Donnan dialysis with ion-exchange membranes was studied under various kinds of experimental conditions using ions of different valences. The diffusion coefficients (D{sub d}) of various kinds of ions in the ion-exchange membrane were obtained by curve fitting an equation derived from the mass balance to three kinds of Donnan dialytic experiments. It was found that the value of D{sub d}/D{sub s} using D{sub d} of monovalent ions in Donnan dialysis with a set of monovalent feed ions and bivalent driving ions was 1/175, where D{sub s} represents a diffusion coefficient in solution. D{sub s} was calculated from the Nernst-Einstein equation substituted by the ionic conductance of ions at infinite dilution in water. Using D{sub d} of bivalent ions in Donnan dialysis with the same set led to a D{sub d}/D{sub s} value of 1/438. Moreover, using D{sub d} in Donnan dialysis with the same set, the value of D{sub d}/D{sub e} was kept ...
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 ...
Precondition for the low-NO operation of gas turbine burners is the generation of as homogeneous as possible a gaseous fuel/air mixture. Such a burner behavior can be achieved by selective adaptation of the flow and injection. A computing model is introduced that calculates the dynamically balanced GASs/droplet flow through numerical solution of the transport equations in a curved orthogonal coordinate system. The feedback reaction of the fuel droplets to the gas phase is considered by source terms and/or sink terms. Test computations were carried out for comparison with analytical solutions from the flow mechanics and validated by comparison with measurement results. Droplet motion and evaporation were checked by means of published experimental results on single droplets. For the validation of the developed computing technique, model tests were carried out with water in place of fuel. It is pointed out that the choice of different parameters ...
The preparation of neem oil microemulsion and its acaricidal activity in vitro was developed in this study. In these systems, the mixture of Tween-80 and the sodium dodecyl benzene sulfonate (SDBS) (4:1, by weight) was used as compound surfactant; the mixture of compound surfactant and hexyl alcohol (4:1, by weight) was used as emulsifier system; the mixture of neem oil, emulsifier system and water (1:3.5:5.5, by weight) was used as neem oil microemulsion. All the mixtures were stired in 800rpm for 15min at 40^oC. The acaricidal activity was measured by the speed of kill. The whole lethal time value of 10% neem oil microemulsion was 192.50min against Sarcoptes scabiei var. cuniculi larvae in vitro. The median lethal time value was 81.7463min with the toxicity regression equations of Y=-6.0...
This article discussed a project formed to develop innovative seismic imaging technologies. The Kaleidoscope project aims to accelerate the processing of seismic sound waves by several orders of magnitude using advanced computer-based techniques to reveal oil and gas deposits buried deep in the earth in a manner that is both cost-effective and efficient. It is expected that the technology will be used in deep water applications with massive untapped reserves. The seismic technology will be used to locate hydrocarbons and oil reserves buried 20,000 feet beneath the seabed. It is estimated that the Gulf of Mexico contains 56 billion barrels of oil equivalent worth nearly $6 trillion. However, the reserves are difficult to locate due to the interbedded salt bodies in the subsurface. The project is simultaneously working on both hardware and software applications, and are currently writing the first petascale set of seismic imaging applications. The technology will ...
The pozzolanic behaviour of compound-activated red mud-coal gangue has been investigated through TG, DTA, XRD, FTIR and {sup 27}Al MAS NMR. From viewpoint of reaction kinetics, it is found that the pozzolanic reaction mechanism of the compound-activated red mud-coal gangue-lime system is clearly consistent with diffusion control up to 14 days, and the reaction rate constant calculated from Jander equation decreases with the increase of CaO addition in the system. The hydration products formed in the red mud coal gangue-lime systems at ambient temperature are essentially aluminous C-S-H and Ca{sub 3}Al{sub 2}O{sub 6} center dot xH{sub 2}O. From TG analysis results, it is thought that the high amount of Ca(OH){sub 2} in the pastes of studied system is not conducive to the continual increase of non-evaporable water content of the hydration products. Of particular interest. {sup 27}Al MAS NMR proved to be an effective technique to obtain valuable ...
The pozzolanic behaviour of compound-activated red mud-coal gangue has been investigated through TG, DTA, XRD, FTIR and 27Al MAS NMR. From viewpoint of reaction kinetics, it is found that the pozzolanic reaction mechanism of the compound-activated red mud-coal gangue - lime system is clearly consistent with diffusion control up to 14 days, and the reaction rate constant calculated from Jander equation decreases with the increase of CaO addition in the system. The hydration products formed in the red mud-coal gangue - lime systems at ambient temperature are essentially aluminous C-S-H and Ca3Al2O6.xH2O. From TG analysis results, it is thought that the high amount of Ca(OH)2 in the pastes of studied system is not conducive to the continual increase of non-evaporable water content of the hydration products. Of particular interest, 27Al MAS NMR proved to be an effective technique to obtain valuable information of Al[4] in C-S-H and Al[6] in ...
Isotope separation of Ar, C, /sup 3/He, Kr, Ne, O, and Xe isotopes is reported. TiFeH/sub x/, TiCoH/sub x/, TiCuH/sub x/, and VH/sub x/ were studied using NMR (proton relaxation times). VD/sub x/ and VT/sub x/ were synthesized. The problem of calculating the valence state of Pu is discussed. A series solution to the plutonium (N,H) characteristic equation is suggested. Shipments of /sup 231/Pa, /sup 230/Th, and /sup 229/Th are reported. Separation and processing of /sup 234/U are also reported. Theoretical methods were developed to calculate temperature distributions as functions of water flow rate in liquid thermal diffusion columns. Diffusion coefficients were measured from 300 to 1200/sup 0/K for Kr-Xe and Kr-Ar. New thermal diffusion factors are submitted for Ne-Ar.
Flow structures and heat transfer at the liquid-liquid interface are investigated experimentally and analytically in the entrance region of a liquid metal-water stratified flow in a horizontal rectangular channel. The two-equation model and the mixing length mode, including a damping turbulent region near the interface, are used to obtain analytical results for pressure drop, position of the interface, and heat transmission coefficients across the interface. Analytical results are compared with experimental results. Conclusions are presented and include the following: (1) vertical position of the interface increases or decreases along the stream according to the fluid flow rate of both fluids, and then approaches a constant value in the fully developed region; (2) pressure drop increases along the stream from the minimum value at the inlet; and (3) heat transmission coefficients across the interface decrease in the entrance region to a minimum ...
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 release rate of a nuclide from a reactor or a radioactive waste disposal plant at the accident is not steady, but varies with time. The various parameters of a nuclide migration into environment vary also day after day, or with the seasons. In such cases, dynamic behavior of the nuclide in the environment must be taken into consideration. It is difficult for a mathematical model to involve all of mechanisms for the nuclide migration. The environment for evaluation of doses are usually divided into some of compartments in which a nuclide concentration is uniform. Time variations of the nuclide concentration in the compartment are described in simultaneous differential equations. The nuclide concentration can be solved as a time function, and the radiation doses, therefore, can be estimated as a time function. Generic analysis code for dynamic compartment model (GACOM) is developed for the nuclide migration and the evaluation of doses in terrestrial biosphere. ...
The CATHENA (formerly ATHENA) has been used to simulate the thermalhydraulic behaviour of the WOLSUNG-1 CANDU-600 reactor during the D_20 spill incident of 1984 November 25. A 4-inch (nominal) Liquid Relief Valve inadvertently opened in the reactor auxiliary system during normal reactor operation, resulting in a discharge of heavy water from the primary heat transport system. The valve remained open for approximately 29 minutes. CATHENA is an advanced thermalhydraulic computer code for analysis of postulated loss-of-coolant accidents (LOCA) and transient faults in CANDU nuclear reactors. A full two-fluid (six-equation) representation of the two-phase flow is used. Component models are used to represent pumps, valves, critical discharge, etc., which are necessary to describe the behaviour of the CANDU system under upset conditions. Heat transfer between the fluid and piping walls (or fuel) is modelled using applicable correlations for boiling, ...
A safety analysis code to design a pressurized water reactor and to obtain the licences including entire proprietary rights is under development in domestic research and development project. The purpose and scope of this report is to develop the flow regimes related models for inter-phase friction, wall frictions, wall heat transfer, and inter-phase heat and mass transfer in two-phase three-field equations. In order to choose choose the flow regime criteria, we have investigated various exiting best-estimate T/H codes in this chapter 2. They are the RELAP5-3D, TRAC-M, CATHARE, MARS codes. Around 500 references used in these codes have been collected and reviewed. Also we have investigated eleven papers in detail. In chapter 3, based on the selected flow regimes, the flow regime maps for a gas-liquid flow in horizontal and vertical tubes have decided including the mechanisms of flow regime transition regions. Conclusively, the process will be ...
Programs & Services Technical Resources Land Use Soils Water Snow Survey & Water Supply Water Management Water Quality Watersheds Wetlands Air Plants & Animals Energy Climate...
Motivated by the numerical investigations of Laval, Dubrulle & Nazarenko (1999), we develop a quasilinear theory of the 2D Euler equation and derive an integro-differential equation for the evolution of the coarse-grained vorticity. This equation respects all the invariance properties of the Euler equation and conserves angular momentum in a circular domain and linear impulse in a channel (as well as in an infinite domain). The explicit energy is not rigorously conserved as it is partly transfered into fine-grained fluctuations but the total energy is conserved. We prove a H-theorem for the Fermi-Dirac entropy and make the connection with statistical theories of 2D turbulence.
We carry out a Lie group analysis of the Sachs equations for a time-dependent axisymmetric non-rotating space-time in which the Ricci tensor vanishes. These equations, which are the first two members of the set of Newman-Penrose equations, define the characteristic initial-value problem for the space-time. We find a particular form for the initial data such that these equations admit a Lie symmetry, and so defines a geometrically special class of such spacetimes. These should additionally be of particular physical interest because of this special geometric feature.
The top quark antiquark production system in the dilepton decay channel is described by a set of equations which is nonlinear in the unknown neutrino momenta. Its most precise and least time consuming solution is of major importance for measurements of top quark properties like the top quark mass and tt spin correlations. The initial system of equations can be transformed into two polynomial equations with two unknowns by means of elementary algebraic operations. These two polynomials of multidegree two can be reduced to one univariate polynomial of degree four by means of resultants. The obtained quartic equation is solved analytically.
n an early approach, we proposed a kinetic model with multiple translational temperature [K. Xu, H. Liu and J. Jiang, Phys. Fluids {\\bf 19}, 016101 (2007)], to simulate non-equilibrium flows. In this paper, instead of using three temperatures in $x-$, $y-$, and $z$-directions, we are going to further define the translational temperature as a second-order symmetric tensor. Based on a multiple stage BGK-type collision model and the Chapman-Enskog expansion, the corresponding macroscopic gas dynamics equations in three-dimensional space will be derived. The zeroth-order expansion gives the 10 moment closure equations of Levermore [C.D. Levermore, J. Stat. Phys {\\bf 83}, pp.1021 (1996)]. To the 1st-order expansion, the derived gas dynamic equations can be considered as a regularization of Levermore's 10 moments equations. The new gas dynamic equations have the same structure as the ...
The interfacial area transport equation dynamically models two-phase flow regime transitions and predicts continuous changes of the interfacial area concentration along the flow field. It replaces the flow regime-dependent correlations for the interfacial area concentration in thermal-hydraulic system analysis. In the present study, detailed formulation of the interfacial area transport equation is presented along with its evaluation results based on the detailed benchmark experiments. In view of model evaluation, the equation is simplified into one-dimensional steady state one-group interfacial area transport equation. The prediction made by model agrees well with the experimental data obtained in round pipes of various diameters. The framework for the two-group transport equation and the necessary constitutive relations are also presented in view of bubble transport of various ...
To develop the advanced subchannel analysis code, the dominant factors that influence the boiling transitional process must be taken into account in the mechanistic constitutive equations based on the flow geometries and the fluid properties. The dominant factors that influence the boiling transitional processes are (1) the gas-liquid re-distribution by cross flow, (2) the liquid film dryout, (3) the two-phase flow regime transition, (4) the droplet deposition, and (5) the spacer-droplet interaction. At first, we indicated the strategy for the development of the constitutive equations for the five dominant factors based on the experimental database by the latest measurement technique and the latest computational fluid dynamics method. Then, the problems of the present constitutive equations and the improvement plan of the constitutive equations were indicated. Finally, the layered structure for the ...
A sodium-water reaction drove from the single tube break in steam generator of FBR might overheat labor tubes rapidly under internal pressure loadings. If the temperature of tube wall becomes too high, it has to be evaluated that the stress of tube does not exceed the material strength limit to prevent the propagation of tube rupture. This study clarified the tensile and creep properties of Mod.9Cr-1Mo steel at ultra-high temperature which will be used in evaluation of the tube burst by sodium-water reaction. The strain rates for tensile test are from 10%/min to 10%/sec, and creep-rupture time is maximum 277sec. The range of test temperature is 700degC to 1300degC. The main results obtained were as follows; (1) The evaluation data on the relationship between tensile strength and strain rate and creep-rupture strength in shorter time on Mod.9Cr-1Mo steel were acquired. (2) Short-term mechanical properties of Mod.9Cr-1Mo steel were evaluated ...
The objective of the study was the kinetic modeling of the effect of storage temperature on the quality and shelf life of chilled fish, modified atmosphere-packed (MAP), and osmotically pretreated with the addition of nisin as antimicrobial agent. Fresh gilthead seabream (Sparus aurata) fillets were osmotically treated with 50% high dextrose equivalent maltodextrin (DE 47) plus 5% NaCl. Water loss, solid gain, salt content, and water activity were monitored throughout treatment and treatment conditions were selected for the shelf life study. Untreated and osmotically pretreated slices with and without nisin (2 x 10(4) IU/100 g osmotic solution), packed in air or modified atmosphere (50% CO(2)-50% air), and stored at controlled isothermal conditions (0, 5, 10, and 15 degrees C) were studied. Quality assessment and modeling were based on growth of several microbial indices, total volatile nitrogen, trimethylamine nitrogen, lipid oxidation ...
for? Learn the Issues Science & Technology Laws & Regulations About EPA Contact Us Water: Water Quality Criteria You are here: Water Science & Technology Surface Water Standards &...
This document is a reference guide to the Xyce Parallel Electronic Simulator, and is a companion document to the Xyce Users Guide. The focus of this document is (to the extent possible) exhaustively list device parameters, solver options, parser options, and other usage details of Xyce. This document is not intended to be a tutorial. Users who are new to circuit simulation are better served by the Xyce Users Guide. The Xyce Parallel Electronic Simulator has been written to support, in a rigorous manner, the simulation needs of the Sandia National Laboratories electrical designers. It is targeted specifically to run on large-scale parallel computing platforms but also runs well on a variety of architectures including single processor workstations. It also aims to support a variety of devices and models specific to Sandia needs. This document is intended to complement the Xyce Users Guide. It contains comprehensive, detailed information about a number of topics ...
In this paper we apply stochastic programming modelling and solution techniques to planning problems for a consortium of oil companies. A multiperiod supply, transformation and distribution scheduling problem - the Depot and Refinery Optimization Problem (DROP) - is formulated for strategic or tactical level planning of the consortium's activities. This deterministic model is used as a basis for implementing a stochastic programming formulation with uncertainty in the product demands and spot supply costs (DROPS), whose solution process utilizes the deterministic equivalent linear programming problem. We employ our STOCHGEN general purpose stochastic problem generator to 'recreate' the decision (scenario) tree for the unfolding future as this deterministic equivalent. To project random demands for oil products at different spatial locations into the future and to generate random fluctuations in their future prices/costs a stochastic input data simulator is developed and calibrated to ...
Staggered arrays of dimples printed on opposite surfaces of a cooling channel is formulated numerically and optimized with hybrid multi-objective evolutionary algorithm and Pareto optimal front. As Pareto optimal front produces a set of optimal solutions, the trends of objective functions with design variables are predicted by hybrid multi-objective evolutionary algorithm. The problem is defined by three non-dimensional geometric design variables composed of dimpled channel height, dimple print diameter, dimple spacing, and dimple depth, to maximize heat transfer rate compromising with pressure drop. Twenty designs generated by Latin hypercube sampling were evaluated by Reynolds-averaged Navier-Stokes solver and the evaluated objectives were used to construct Pareto optimal front through hybrid multi-objective evolutionary algorithm. The optimum designs were grouped by k-means clustering technique and some of the clustered points were evaluated by flow analysis. ...
Staggered arrays of dimples printed on opposite surfaces of a cooling channel is formulated numerically and optimized with hybrid multi-objective evolutionary algorithm and Pareto optimal front. As Pareto optimal front produces a set of optimal solutions, the trends of objective functions with design variables are predicted by hybrid multi-objective evolutionary algorithm. The problem is defined by three non-dimensional geometric design variables composed of dimpled channel height, dimple print diameter, dimple spacing, and dimple depth, to maximize heat transfer rate compromising with pressure drop. Twenty designs generated by Latin hypercube sampling were evaluated by Reynolds-averaged Navier-Stokes solver and the evaluated objectives were used to construct Pareto optimal front through hybrid multi-objective evolutionary algorithm. The optimum designs were grouped by k-means clustering technique and some of the clustered points were evaluated by flow analysis. ...
Aliasing, or sharing, is prominent in many domains, denoting that two differently-named objects are in fact identical: a change in one object (memory cell, circuit terminal, disk block) is instantly reflected in the other. Languages for modelling such domains should let the programmer explicitly define the sharing among objects or expressions. A DSL compiler may find other identical expressions and share them, implicitly. Such common subexpression elimination is crucial to the efficient implementation of DSLs. Sharing is tricky in embedded DSL, since host aliasing may correspond to copying of the underlying objects rather than their sharing. This tutorial summarizes discussions of implementing sharing in Haskell DSLs for automotive embedded systems and hardware description languages. The technique has since been used in a Haskell SAT solver and the DSL for music synthesis. We demonstrate the embedding in pure Haskell of a simple DSL with a language form for ...
An important issue towards a broader acceptance of answer-set programming (ASP) is the deployment of tools which support the programmer during the coding phase. In particular, methods for debugging an answer-set program are recognised as a crucial step in this regard. Initial work on debugging in ASP mainly focused on propositional programs, yet practical debuggers need to handle programs with variables as well. In this paper, we discuss a debugging technique that is directly geared towards non-ground programs. Following previous work, we address the central debugging question why some interpretation is not an answer set. The explanations provided by our method are computed by means of a meta-programming technique, using a uniform encoding of a debugging request in terms of ASP itself. Our method also permits programs containing comparison predicates and integer arithmetics, thus covering a relevant language class commonly supported by all state-of-the-art ASP ...
A computational fluid dynamics (CFD) is coupled with a computational structural dynamics (CSD) to simulate the unsteady rotor flow with aeroelasticity effects. An unstructured upwind Navier-Stokes solver was developed for this simulation, with 2nd order time-accurate dual-time stepping method for temporal discretization and low Mach number preconditioning method. For turbulent flows, both the Spalart-Allmaras and Menter's SST model are available. Mesh deformation is achieved through a fast dynamic grid method called Delaunay graph map method for unsteady flow simulation. The rotor blades are modeled as Hodges & Dowell's nonlinear beams coupled flap-lag-torsion. The rotorcraft computational structural dynamics code employs the 15-dof beam finite element formulation for modeling. The structure code was validated by comparing the natural frequencies of a rotor model with UMARC. The flow and structure codes are coupled tightly with information exchange several ...
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.
The static thermophysical properties of the molten salt system like LiF-NaF-BeF_2 influence the design and construction of the fuel salt and coolant in the Molten Salt Reactor for the new generation. In this paper, the equation of state of the ternary system 0.15LiF-0.58NaF-0.27BeF_2, over the temperature range from 873.15K to 1073.15K at one atmosphere pressure, is described by using modified Peng-Robinson equation. The density of the ternary system is evaluated by this equation directly, and compared with the experimental data. Base on the equation of state, the other static thermophysical properties such as the enthalpy, entropy and heat capacity at constant pressure are evaluated by the fugacity coefficient and residual function methods respectively. The density calculated by Peng-Robinson equation is in highly agreement with the experimental data, and the enthalpy, entropy and ...
The static thermodynamic properties of the molten salt system like LiF-NaF-BeF_2 influence the design and construction of the fuel salt and coolant in the Molten Salt Reactor for the new generation. In this paper, the equation of state of the ternary system 15%LiF-58%NaF-27%BeF_2, over the temperature range of 873.15K to 1073.15K at one atmosphere pressure, is described using Peng-Robinson equation modified by us. And the density of the ternary system is evaluated by this equation directly, and compared with the experimental data. Base on the equation of state, the other static thermodynamic properties such as the enthalpy, entropy and heat capacity at constant pressure are estimated by the residual function method and the fugacity coefficient method respectively. The density calculated by Peng-Robinson equation is in highly agreement with the experimental data, and the enthalpy, ...
The static thermophysical properties of the molten salt system like LiF-NaF-BeF2 influence the design and construction of the fuel salt and coolant in the Molten Salt Reactor for the new generation. In this paper, the equation of state of the ternary system 0.15LiF-0.58NaF-0.27BeF2, over the temperature range from 873.15K to 1073.15K at one atmosphere pressure, is described by using modified Peng-Robinson equation. The density of the ternary system is evaluated by this equation directly, and compared with the experimental data. Base on the equation of state, the other static thermophysical properties such as the enthalpy, entropy and heat capacity at constant pressure are evaluated by the fugacity coefficient and residual function methods respectively. The density calculated by Peng-Robinson equation is in highly agreement with the experimental data, and the enthalpy, entropy and ...
We extend Winger's work on the wave equations for integer-spin particles to the spinorial case. A recent suggestion that the neutrino might be a fermionic tachyon is examined. We point out that a four-component Dirac equation cannot describe fermionic tachyon. (orig.).
An elliptic equation in a rectangle with coefficients depending on a fast variable and with its period being a small parameter is considered. An asymptotic expansion of the solution up to an arbitrary degree of the small parameter is constructed and substantiated by applying the two-scale expansion method.
Allometric equations are often used to extrapolate traits in animals for which only body mass estimates are known, such as dinosaurs. One important decision can be whether these equations should be...Full Text Available
A generalized Fisher equation (GFE) relates the time derivative of the average of the intrinsic rate of growth to its variance. The GFE is an exact mathematical result that has been widely used in population...Full Text Available
The objective of this project was to develop new photocatalytic or other innovative process chemistry for the treatment of pink water and related contaminated water.
We apply the group theory to Kadomtsev-Petviashvili-Burgers (KPBII) equation which is a natural model for the propagation of the two-dimensional damped waves. In correspondence with the generators of the symmetry group allowed by the equation, new types of symmetry reductions are performed. Some new exact solutions are obtained, which can be in the form of solitary waves and periodic waves. Specially, our solutions indicate that the equation may have time-dependent nonlinear shears. Such exact explicit solutions and symmetry reductions are important in both applications and the theory of nonlinear science.
The wave equation for spin;1/2 tachyons is derived from the Dirac equation and the principle of relativity extended to superluminal Lorentz frames. From this wave equation and the Dirac equation infinite velocity spinor transformations are obtained. They yield bispinors of the plane-wave states of the tachyon, their interpretation and covariant orthogonality relations satisfied by them. The transformation properties of the bispinors under Lorentz transformation are discussed. The boundary conditions for the free propagator of wave functions of tachyons are obtained and the propagator is constructed. Then the covariant S-matrix for scattering from an electromagnetic field is derived. It is applied to the scattering of electron-tachyons from the Coulomb field.
P. Galenko et al. proposed a modified Cahn-Hilliard equation to model rapid spinodal decomposition in non-equilibrium phase separation processes. This equation contains an inertial term which causes the loss of any regularizing effect on the solutions. Here we consider an initial and boundary value problem for this equation in a two-dimensional bounded domain. We prove a number of results related to well-posedness and large time behavior of solutions. In particular, we analyze the existence of bounded absorbing sets in two different phase spaces and, correspondingly, we establish the existence of the global attractor. We also demonstrate the existence of an exponential attractor.
... Abstract : The purpose of this study was to evaluate the tooling cost estimating equation of the DAPCA III model and determine if more accurate ...
the progressive simplification of composite sandwich panels. The equations ... The thermal and structural behavior of composite sandwich panels with a ...
Heat transfer processes use fluids which are generally not pure and can react with transfer surfaces. These surfaces are subject to deposits which can be sediments harmful to heat transfer and to integrity of materials. For nuclear plant steam generators, sludge build-up accelerates secondary side corrosion by concentrating chemical species. A major safety problem involved with such a corrosion is the growing of circumferential cracks which are very difficult to detect and size with eddy current probes. With a view to understand and control this problem, it is necessary to develop a mathematical model for the prediction of sludge behavior in PWR steam generators. Based on fundamental principles, this work intends to use different models available in literature for the prediction of the phenomenon leading to the accumulation of sludge particles at the bottom (the tubesheet) of a PWR. For that, a three-dimensional simulation of magnetite particulate fouling with the finite elements code ...
In this Letter the bound states of (2+1) Dirac equation with the cylindrically symmetric $\\delta (r-r_{0})$-potential are discussed. It is surprisingly found that the relation between the radial functions at two sides of $r_{0}$ can be established by an SO(2) transformation. We obtain a transcendental equation for calculating the energy of the bound state from the matching condition in the configuration space. The condition for existence of bound states is determined by the Sturm-Liouville theorem.
One normally identifies the prediction of tachyons with a higher spin problem. We show that the same phenomenon also happens with certain spin-1/2 equations, in exactly the same circumstances. Furthermore, we prove that these equations also have the defect of possessing solutions with complex energies, even with the smallest external fields. The consequences of these results are analyzed.
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.
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.
We report on the calculation of multi-loop Feynman integrals for single-scale problems by means of difference equations in Mellin space. The solution to these difference equations in terms of harmonic sums can be constructed algorithmically over difference fields, the so-called {pi}{sigma}{sup *}-fields. We test the implementation of the Mathematica package Sigma on examples from recent higher order perturbative calculations in Quantum Chromodynamics. (orig.)
In this paper we present a new efficient sixth-order scheme for nonlinear equations. The method is compared to several members of the family of methods developed by Neta (1979) [B. Neta, A sixth-order family of methods for nonlinear equations, Int. J. Comput. Math. 7 (1979) 157-161]. It is shown that the new method is an improvement over this well known scheme.
Volume 2 of the water resources data for the 1992 water year for Oklahoma consists of record of stage, discharge, and water quality of streams; stage, contents, and water quality of lakes or reservoirs; and water levels of ground-water well. The report contains discharge records for 136 gaging stations; stage and contents for 31 lakes or reservoirs; water quality for 50 gaging stations. Also included are 61 partial-record or miscellaneous streamflow stations and 28 ground-water sites.
The linearized collective Schroedinger equation for nuclear quadrupole surface vibrations incorporates a new spin degree of freeedom with a spin value of 3/2. We use this equation to describe the low energy spectrum of certain even-odd Ir nuclei which have a spin 3/2 in their ground state. For that purpose we explicitly introduce collective spin-dependent potentials which simulate the interaction of the valence nucleon with the core. The linearized Schroedinger equation is transformed into an effective Schroedinger equation with collective spin-dependent potentials. Already collective spin-orbit couplings of SO(3) and SO(5) type are sufficient to reproduce the lowest excited states of even-odd Ir nuclei. (orig.).
The linearized collective Schroedinger equation for nuclear quadrupole surface vibrations incorporates a new spin degree of freeedom with a spin value of 3/2. We use this equation to describe the low energy spectrum of certain even-odd Ir nuclei which have a spin 3/2 in their ground state. For that purpose we explicitly introduce collective spin-dependent potentials which simulate the interaction of the valence nucleon with the core. The linearized Schroedinger equation is transformed into an effective Schroedinger equation with collective spin-dependent potentials. Already collective spin-orbit couplings of SO(3) and SO(5) type are sufficient to reproduce the lowest excited states of even-odd Ir nuclei. (orig.).
Abstract The importance of evaluating greenhouse gas (GHG) emissions from dairy cows within the whole farm setting is being realized as more important than evaluating these emissions in isolation. Current whole farm models aimed at evaluating GHG emissions make use of simple regression equations to predict enteric methane (CH4) production. The objective of the current paper is to evaluate the performance of nine CH4 prediction equations that are currently being used in whole farm GHG models. Data used to evaluate the prediction equations came from a collection of individual (IND) and treatment averaged (TRT) data. Equations were compared based on mean square prediction error (MSPE) and concordance correlation coefficient (CCC) analysis. In general, predictions were poor, with root MSPE (as...
The molten salt reactor (MSR), which is one of the generation IV reactors, can meet the demand of transmutation and breeding. The thermodynamic properties of the molten salt system like LiF-NaF-BeF2 influence the design and construction of the fuel salt and coolant in the MSR for the new generation. In this paper, the equation of state of the ternary system 15%LiF-58%NaF-27%BeF2, over the temperature range from 873.15 to 1 073.15 K at one atmosphere pressure, is described using a modified Peng-Robinson (PR) equation. The densities of the ternary system and its components are estimated by this equation directly, and compared with the experimental data. Based on the equation of state, the other thermodynamic properties such as the enthalpy, entropy and heat capacity at constant pressure are ...
In this paper, a complete solution of the steady-state Newton-Euler equations describing the motion of the balls in high speed thrust bill bearings is presented. Kinematic equations are derived and used as ancillary equations to obtain the numerical solution of the Newton-Euler equations. This differs from previously reported work where extraneous assumptions (instead of kinematic analysis) are used to allow solutions of these equations. It is shown that the no-slip condition for the ball motion reported in the literature is based on an unnecessary assumption. The correct form of the condition is developed and used in the analysis of the ball motion. The analysis shows the importance of contact angle on both ball motion and contact force. An approximate theory is then developed for determining the optimum contact angle that would result in minimum contact force in high speed ball ...
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 ...
Real-time mission-oriented embedded systems are much more difficult to design than ordinary software systems. They require highly reliable and efficient implementations to satisfy mission and time constraints imposed by the applications. The Ada language was designed to facilitate real-time-system software development. However, for many programmers the size and complexity of Ada itself are of concern. In the assertive programming paradigm, computations are specified as sets of assertions about properties of the solution, and not as a sequence of procedural steps. Solving procedures are automatically generated from the assertive description. Real-time programming for mission-oriented systems is supported by equational languages in which assertions are expressed as algebraic equations. Programs written in equational languages are concise, free from implementation details, and easily amenable to verification and parallel ...
This paper describes a semi-empirical calculation of the air-broadened half-widths and the air pressure-induced frequency shifts for the H_2"1"6O isotopologue. This semi-empirical calculation is based on fits of several recent high-quality measurements and theoretical calculations to the first-order terms in the expansion of the complex Robert-Bonamy (CRB) equations, which yields a second- and first-order polynomial function of the differences in the upper- and lower-state vibrational quantum numbers for the half-width and line shift, respectively. The aim of this work was to obtain a complete set of air-broadened half-widths and air pressure-induced frequency shifts for transitions of H_2"1"6O present in the HITRAN database from microwave to the visible in order to supplement the observed and calculated values. For around 700 sets of rotational quantum numbers (J"'K_a"'K_c"'<-J"''K_a"''K_c"''), semi-empirical coefficients describing the vibrational dependence ...
This report documents a one-dimensional numerical model, SAMFT1D, developed to simulate single-phase and multiphase fluid flow and solute transport in variably saturated porous media. The formulation of the governing equations and the numerical procedures used in the code for single-phase and multiphase flow and transport are presented. The code is constructed to handle single-phase as well as two or three-phase flow conditions using two integrated sets of computational modules. The fully implicit scheme is used in the code for both single-phase and multiphase flow simulations. Either the Crank-Nicholson scheme or the fully implicit scheme may be used in the transport simulation. The single-phase modules employ the Galerkin and upstream weighted residual finite element techniques to model flow and transport of water (aqueous phase) containing dissolved single-species contaminants concurrently or sequentially, and include the treatment of ...
We study the surface and interface roughness of a two-liquid film system using neutron reflectometry. The middle phase and bottom phase are extracted from a three-phase microemulsion containing 2 wt% non-ionic surfactant C{sub 10}E{sub 4}, 81 wt% D{sub 2}O and 17 wt% octane (C{sub 8}H{sub 18}). The system is formed by spreading a few drops of the middle phase of the microemulsion on the bottom phase of the microemulsion. The middle phase does not wet the bottom phase, but demixes to form a top octane-water mixture layer with a thickness of {proportional_to}580 A. The mean surface roughness of the system is found to be 7 A but the data are insensitive to interface roughness. An extended capillary wave model for a two-liquid-film system is derived to account for the coupling between the surface and interface. The dispersion relation yields two fundamental modes. In a calculation of mean-square roughness (left angle {zeta}{sup 2} right angle {sup 1/2}) a high-k-cutoff ...
The aim of this study was to assess the prevalence and correlates of disclosure to network members of being hepatitis C virus (HCV)- or human immunodeficiency virus (HIV)-infected among injecting dyads of infected injection drug users (IDUs) in Budapest, Hungary and Vilnius, Lithuania,. Multivariate generalized estimating equations (GEE) were used to assess associations. Very strong infection disclosure norms exist in Hungary, and HCV disclosure was associated with using drugs and having sex within the dyad. Non-ethnic Russian IDUs in Lithuania were more likely to disclose HCV infection to non-Roma, emotionally close and HCV-infected network members, and to those with whom they shared cookers, filters, drug solutions or rinse water or got used syringes from, and if they had fewer non-IDU or IDU network members. Ethnic Russian Lithuanian IDUs were more likely to disclose HCV if they had higher disclosure attitude and knowledge scores, 'trusted' ...
The Hyperion project was developed to determine an algorithm for assessing the risk of hydrate plug formation in the pipeline transport oil-water-gas mixtures at low temperatures. The project is a collaboration between physicists, chemists and engineers within the petroleum industry. This paper provided an overview of the project and outlined results obtained as the project entered its third and final year. The main objective of the project has been to understand the inherent mitigation effects of some oils on gas hydrate formation as well as to develop methods of predicting the risk of hydrate plugging. To date, the project has extracted and studied natural inhibiting components (NICs) in oils. Molecular modelling techniques have been used to study hydrate and fluid interfaces in order to estimate the driving force of agglomeration and growth through mechanical surface stress and measurement of surface wave fluctuations. A scheme is also being developed measure ...
Isothermal reaction experiments were performed in the temperature range of 1000 - 1300 C in order to determine the chemical interactions between Zircaloy-4 fuel rod cladding and Inconel-718 spacer grids of Pressurized Water Reactors (PWR) under severe accident conditions. It was not possible to apply even higher temperatures since fast and complete liquefaction of the components occurred as a result of eutectic interactions during heatup. The liquid reaction products formed enhance and accelerate the degradation of the material couples and the fuel elements, respectively. Only small amounts of Inconel are necessary to liquefy large amounts of Zircaloy. Thin oxide layers on the Zircaloy surface delay the beginning of the chemical interactions with Inconel but cannot prevent them. In this work the reaction kinetics have been determined for the system: as-received and pre-oxidized Zircaloy-4/Inconel 718. The interactions can be described by parabolic rate laws; the ...
Combining the approaches given by Baldwin [Baldwin D et al. Symbolic computation of exact solutions expressible in hyperbolic and elliptic functions for nonlinear PDEs. J Symbol Comput 2004;37:669-705], Peng [Peng YZ. A polynomial expansion method and new general solitary wave solutions to KS equation. Comm Theor Phys 2003;39:641-2] and by Schuermann [Schuermann HW, Serov VS. Weierstrass' solutions to certain nonlinear wave and evolution equations. Proc progress electromagnetics research symposium, 28-31 March 2004, Pisa. p. 651-4; Schuermann HW. Traveling-wave solutions to the cubic-quintic nonlinear Schroedinger equation. Phys Rev E 1996;54:4312-20] leads to a method for finding exact travelling wave solutions of nonlinear wave and evolution equations (NLWEE). The first idea is to generalize ansaetze given by Baldwin and Peng to find elliptic solutions of NLWEEs. Secondly, conditions used by ...
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 ...
The biosorption of Cu(II) from aqueous solutions by valonia tannin resin was investigated as a function of particle size, initial pH, contact time and initial metal ion concentration. The aim of this study was to understand the mechanisms that govern copper removal and find a suitable equilibrium isotherm and kinetic model for the copper removal in a batch reactor. The experimental isotherm data were analysed using the Langmuir, Freundlich and Temkin equations. The equilibrium data fit well in the Langmuir isotherm. The experimental data were analysed using four sorption kinetic models - the pseudo-first- and second-order equations, the Elovich and the intraparticle diffusion model equation - to determine the best fit equation for the biosorption of copper ions onto valonia tannin resin. Results show that the pseudo-second-order equation provides the best correlation for the ...
Formulation is introduced for discretizing a boundary integral equation into an indirect boundary element method for the solution of 3-dimensional topographic problems. Yokoi and Takenaka propose an analytical solution-capable reference solution (solution for the half space elastic body with flat free surface) to problems of topographic response to seismic motion in a 2-dimensional in-plane field. That is to say, they propose a boundary integral equation capable of effectively suppressing the non-physical waves that emerge in the result of computation in the wake of the truncation of the discretized ground surface making use of the wave field in a semi-infinite elastic body with flat free surface. They apply the proposed boundary integral equation discretized into the indirect boundary element method to solve some examples, and succeed in proving its validity. In this report, the equation is expanded to ...
Water-resources data for the 1993 water year for Minnesota consists of records of stage, discharge and water quality of streams; stage, contents, and water quality of lakes and reservoirs; and water levels and water quality in wells and springs. This volume contains discharge records for 58 gaging stations; stage and contents for 9 lakes and reservoirs; water quality 15 stream stations, 1 lakes station, 1 precipitation station, and water levels for 15 observation wells. Also included are 61 high-flow partial-record stations. Additional water data were collected at various sites, not part of the systematic data collection program and are pulished as miscellaneous measurements.
Water-resources data for the 1991 water year for Minnesota consist of records of stage, discharge and water quality of streams; stage, contents, and water quality of lakes and reservoirs; and water levels and water quality in wells and springs. The volume contains discharge records for 60 gaging stations; stage and contents for 9 lakes and reservoirs; water quality for 18 stream stations, 1 lake station, 22 partial-record sites, 1 precipitation station, 108 wells; and water levels for 119 observation wells. Also included are 59 high-flow partial-record stations. Additional water data were collected at various sites, not part of the systematic data collection program and are published as miscellaneous measurements or low-flow investigations.
In the National Water Plan it is described which measures must be taken to keep the Netherlands safe and livable for future generations and to make use of the chances offered by water.
swimmers in or around the water. Designate a responsible adult who can swim and knows CPR to watch swimmers in or around water. The supervisor should not be involved in any...
Solving a coupled set of differential equations formed by a Fokker-Planck equation and a set of master equations, we can calculate the number of 'prefission' particles (n,p,a), which are emitted from a hot, rotating compound nucleus. Furthermore we get their energy spectrum and their angular momentum distribution. Our results depend sensitively on the friction parameter #gamma#, the inertia M and the fission potential. Comparing the results with experimental data, we may learn something about these quantities. The investigations have been made for "1"8"7Ir and "1"8"5Os. (orig.).
Solving a coupled set of differential equations formed by a Fokker-Planck equation and a set of master equations, we can calculate the number of 'prefission' particles (n,p,a), which are emitted from a hot, rotating compound nucleus. Furthermore we get their energy spectrum and their angular momentum distribution. Our results depend sensitively on the friction parameter {gamma}, the inertia M and the fission potential. Comparing the results with experimental data, we may learn something about these quantities. The investigations have been made for {sup 187}Ir and {sup 185}Os. (orig.).
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.
Constrained quantum dynamics is used to propose a nonlinear dynamical equation for pure states of a generalized coarse-grained system. The relevant constraint is given either by the generalized purity or by the generalized invariant fluctuation, and the coarse-grained pure states correspond to the generalized coherent i.e. generalized nonentangled states. Open system model of the coarse-graining is discussed. It is shown that in this model and in the weak coupling limit the constrained dynamical equations coincide with an equation for pointer states, based on Hilbert-Schmidt distance, that was previously suggested in the context of the decoherence theory.
Exact solutions are found for the Chandrasekhar Page angular equation which results when the Dirac equation in a Kerr Newman space time is separated into its radial and angular parts. The solutions turn out to be remarkably simple in form while satisfying the asymptotic conditions deduced earlier. The eigenvalues are found to be the square root of the total angular momentum as first found by Dirac for flat space; supplemented by a term which is the product of the mass of the Dirac particle times the specific angular momentum of the black hole. The additional contribution is what is expected from frame dragging.
Big Bang nucleosynthesis requires a fine balance between equations of state for photons and relativistic fermions. Several corrections to equation of state parameters arise from classical and quantum physics, which are derived here from a canonical perspective. In particular, loop quantum gravity allows one to compute quantum gravity corrections for Maxwell and Dirac fields. Although the classical actions are very different, quantum corrections to the equation of state are remarkably similar. To lowest order, these corrections take the form of an overall expansion-dependent multiplicative factor in the total density. We use these results, along with the predictions of Big Bang nucleosynthesis, to place bounds on these corrections.
The spherical harmonic analysis of cosmic ray neutron data from the worldwide network neutron monitor stations during the years 1966 to 1969 was carried out. The second zonal harmonic component obtained from the analysis corresponds to the Pole-Equator anisotropy of the cosmic ray neutron intensity. Such an anisotropy makes a semiannual variation. In addition to this, it is shown that the Pole-Equator anisotropy makes a variation depending on the interplanetary magnetic field (IMF) sector polarities around the passages of the IMF sector boundary. A mechanism to interpret these results is also discussed.
Pippenger (2011) recently proposed a solution to the longstanding forward-bias puzzle. He argues that the puzzling estimates obtained using the standard equation for the efficient markets hypothesis are due to omitted variable bias. He identifies the missing variables as the future change in the forward exchange rate and the future interest differential. When these are added to the standard equation, he finds a one-to-one relationship between the future change in the spot rate and the forward premium. However, we argue that his equation can only test covered interest parity and offers no insight into the forward-bias puzzle.
pletely new PW fuel cell model had to be developed for CINDA. Several ...... following linear equation, which is used in the SINDA fuel cell model: ...
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.).
Within the framework of a modified Balitsky-Kovchegov equation, we calculated and provide estimates of non-linear saturation effects expected in the LHC range of energies.
In this paper, research on pass schedule of tandem cold mill (TCM) is carried out. According to load (reduction, rolling force, motor power) balance, non-linear equations set with variables of inter-stand thickness is constructed. The pass schedule optimization is carried out by solving the non-linear equations set. As the traditional method, the Newton-Raphson method is used for solving the non-linear equations set. In this paper a new simple method is brought up. On basis of the monotone relations between thickness and load, the inter-stands thickness is adjusted dynamically. The solution of non-linear equations set can be converged by iterative calculation. This method can avoid the derivative calculation used by traditional method. So, this method is simple and calculation speed is high. It is suitable for on-line control. (author)
In this paper, research on pass schedule of tandem cold mill (TCM) is carried out. According to load (reduction, rolling force, motor power) balance, non-linear equations set with variables of inter-stand thickness is constructed. The pass schedule optimization is carried out by solving the non-linear equations set. As the traditional method, the Newton-Raphson method is used for solving the non-linear equations set. In this paper a new simple method is brought up. On basis of the monotone relations between thickness and load, the inter-stands thickness is adjusted dynamically. The solution of non-linear equations set can be converged by iterative calculation. This method can avoid the derivative calculation used by traditional method. So, this method is simple and calculation speed is high. It is suitable for on-line control. (author)
The geographical arrangement of people on the surface of the earth is approximated by a mathematical equation of 361 terms. This is a convenient form for comparison with other distributions and for...Full Text Available
The validity is given to the newly proposed two weighting {delta}f scheme (Wang et al., Research Report of National Institute for Fusion Science NIFS-588, 1999) for collisional or neoclassical transport calculations, which can solve the drift kinetic equation taking account of effects of steep plasma gradients, large radial electric field, finite banana width, and the non-standard orbit topology near the axis. The marker density functions in weight equations are successively solved by using the idea of {delta}f method and a hierarchy of equations for weight and marker density functions is obtained. These hierarchy equations are solved by choosing an appropriate source function for each marker density. Thus the validity of the two weighting {delta}f scheme is mathematically proved. (author)
Electrodynamic properties of processes involving charged tachyons are considered in connection with the real superluminal Lorentz transformations. These transformations are interpreted without making use of a tachon corridor and so that the interaction of tachyons from subluminal sources with ordinary matter is Lorentz invariant. Transformations of the electromagnetic fields are deduced from the extended principle of relativity. They are necessary in order to obtain the field equations that are satisfied by fields obeying Maxwell equations relative to superluminal inertial frames with respect to subluminal inertial frames. The field equations allow one to get the dependence of the field energy density and its current on the field strengths. The resulting equations are applied in a discussion of plane electromegnetic waves, sent out by a superluminal source. The obtained electrodynamics allows the ...
a large number of polycrystalline materials, has the fot_ of the. Manson-Coffin- Basquin equation, wherein the exponents of the two power law terms are given ...
Most gas dynamic computations in industrial ducts are done in one dimension with cross-section-averaged Euler equations. This poses a fundamental difficulty as soon as geometrical discontinuities are present. The momentum equation contains a non-conservative term involving a surface pressure integral, responsible for momentum loss. Definition of this integral is very difficult from a mathematical standpoint as the flow may contain other discontinuities (shocks, contact discontinuities). From a physical standpoint, geometrical discontinuities induce multidimensional vortices that modify the surface pressure integral. In the present paper, an improved 1D flow model is proposed. An extra energy (or entropy) equation is added to the Euler equations expressing the energy and turbulent pressure ...
The solution of the continuity equation in practical applications is examined, and the values needed for approximate solutions are indicated. Models are adequate for investigating what could happen,...Full Text Available
In this report we present a suite of subroutines for the solution of sparse unsymmetric sets of linear equations using a variant of Gaussian elimination. The subroutines are divided into three distinct phases. The first phase optionally preorders the matrix to block triangular form and then performs a sparsity oriented factorization, the second factorizes a matrix of a similar sparsity pattern, while the third uses these decompositions to solve the equations. In this revised edition, the actual Fortran listings have been replaced by a reference to their availability in machine readable form. Other changes to the first edition are minor. This revision is essentially a reprint of the 1979 revision, the only changes being that the MA28A specification sheet has been typeset and has a few minor changes and the existence of a version of the package for complex equations is announced. 11 refs., 11 tabs.
The objective of this research is to replace the current flow regime dependent interfacial area correlations in the thermal-hydraulic system analysis codes with an interfacial area transport equation, which dynamically models two-phase flow regime transitions. The interfacial area transport equation applicable to bubbly flows has been developed based on mechanistic modeling of bubble interaction. Detailed experiments have been performed for bubbly conditions in 2.54, 5.08 and 10.16-cm ID pipes to benchmark the model. This new approach predicts the continuous changes of the interfacial area and will eliminate artificial bifurcations stemming from the use of the static flow regime transition criteria. This paper presents the preliminary results of incorporating the one-group transport equation into the USNRC consolidated code, currently under development. The framework for the two-group transport equation ...
No Kerr-like exact solution has yet been found in Chern-Simons modified gravity. Intrigued by this absence, we study stationary and axisymmetric metrics that could represent the exterior field of spinning black holes. For the standard choice of the background scalar, the modified field equations decouple into the Einstein equations and additional constraints. These constraints eliminate essentially all solutions except for Schwarzschild. For non-canonical choices of the background scalar, we find several exact solutions of the modified field equations, including mathematical black holes and pp-waves. We show that the ultrarelativistically boosted Kerr metric can satisfy the modified field equations, and we argue that physical spinning black holes may exist in Chern-Simons modified gravity only if the metric breaks stationarity, axisymmetry or energy-momentum conservation.
A well-known diffuse interface model consists of the Navier-Stokes equations nonlinearly coupled with a convective Cahn-Hilliard type equation. This system describes the evolution of an incompressible isothermal mixture of binary-fluids and it has been investigated by many authors. Here we consider a variant of this model where the standard Cahn-Hilliard equation is replaced by its nonlocal version. More precisely, the gradient term in the free energy functional is replaced by a spatial con- volution operator acting on the order parameter phi. Therefore the coupling with the Navier-Stokes equations is difficult to handle even in two spatial dimensions because of the lack of regularity of phi. We establish the global existence of a weak solution.
major radius of the toroidal containment volume. The containment time was calculated using equation (4) as derived in references 16 and 17. The solid ...
The possible physical linkage between galactic cosmic rays intensity and the Earth's cloud cover is discussed using the analysis of the first indirect aerosol effect (Twomey effect) and its experimental representation as the dependence of average cloud droplet effective radius on aerosol index characterizing the aerosol concentration in the atmospheric air column of unit section. It is shown that the basic kinetic equation of the Earth's climate energy-balance model is described by the bifurcation equation (with respect to the temperature of the Earth's surface) in the form of fold catastrophe with two governing parameters defining the variations of insolation and Earth's magnetic field (or galactic cosmic rays intensity in the atmosphere), respectively. The principle of hierarchical climatic models construction, which consists in the structural invariance of balance equations of these models evolving on the different time ...
The DN-tracking method is used to prove the existence of a closed trajectory in a quadratic system of ordinary differential equations in three dimensions.
... Abstract : The purpose of this paper is to introduce strapdown inertial navigation error equations based on a quaternion relation between body-fixed ...
A new mathematical framework is formulated to derive the effective equations of motion for the constrained quantum system which possesses an internal clock. In the realm close to classical behavior, the quantum evolution is approximated by a finite system of coupled but ordinary differential equations adhered to the weakly imposed Hamiltonian constraint. For the simplified version of loop quantum cosmology in the Bianchi I model with a free massless scalar filed, the resulting effective equations of motion affirm the bouncing scenario predicted by the previous studies: The big bang singularity is resolved and replaced by the big bounces, which take place up to three times, once in each diagonal direction, whenever the directional density approaches the critical value in the regime of Planckian density. It is also revealed that back-reaction arises from the quantum corrections and modifies the precise value of the ...
Diffusion in the extracellular space (ECS) of the brain is constrained by the volume fraction and the tortuosity and a modified diffusion equation represents the transport behavior of many molecules...Full Text Available
In Randall-Sundrum-type brane-world cosmologies, the dynamical equations on the three-brane differ from the general relativity equations by terms that carry the effects of imbedding and of the free gravitational field in the five-dimensional bulk. Instead of starting from an ansatz for the metric, we derive the covariant nonlinear dynamical equations for the gravitational and matter fields on the brane. The local energy-momentum corrections are significant only at very high energies, and in this regime we show that fluid world-lines have a non-gravitational acceleration off the brane. The imprint on the brane of the nonlocal gravitational field in the bulk is more subtle, and we provide a careful decomposition of this effect. The nonlocal energy density determines the tidal acceleration in the off-brane direction, and can oppose singularity formation via the generalized Raychaudhuri equation. Unlike the ...
The interaction of thermal radiation with conduction and convection in thermally developing absorbing, emitting, nongray gas-particulate turbulent suspension flow through a circular tube is investigated. The contribution of thermal radiation is obtained through evaluation of the total hemispherical emittance of the particulate cloud and through evaluation of single band absorptances for molecular gases, modified to account for the interaction with the particles. The governing differential equation is derived as a (nonlinear) energy equation, coupled with integral equations to find the thermal radiation contributions. The energy equation is solved numerically by an implicit finite difference method with an iterative procedure. Qualitative results for Nusselt numbers are shown for a variety and range of parameters, such as optical thickness of particulates and single molecular gas bands, relative gas band ...
The interaction of thermal radiation with conduction and convection in thermally developing absorbing, emitting, non-gray gas particulate turbulent suspension flow through a circular tube is investigated. The contribution of thermal radiation is obtained through evaluation of the total hemispherical emittance of the particulate cloud and through evaluation of single band absorptances for molecular gases, modified to account for the interaction with the particles. The governing differential equation is derived as a (nonlinear) energy equation, coupled with integral equations to find the thermal radiation contributions. The energy equation is solved numerically by an implicit finite difference with its iterative procedure. Qualitative results for Nusselt numbers are shown for a variety and range of parameters, such as optical thickness of particulates and single molecular gas bands, relative gas band ...
The presence of an extraneous light gas must be taken into account in consideration of centrifugal separation of uranium isotopes, when there is inevitable leakage of such a gas through the gas sealant and/or leaking in from the atmosphere. Consideration is first given to the influence of the presence of the light gas on the maximum separative power. Then the basic equation for isotope separation containing a light gas is derived from Hirschfelder's diffusion equations. This equation is solved and the separative performance is expressed in terms of the shape factor and reflux parameter. The formulas for expressing the flow configurations of the gases are obtained for a simple model in which inflow and outflow prevail throughout the centrifuge. The corresponding equation for a model in which the gases flow in two concentric thin streams is also derived. It is concluded that the influence provided by the ...
The solution of the governing steady transport equations for momentum, heat and mass transfer in flowing fluids can be very difficult. These difficulties arise from the nonlinear, coupled, nonsymmetric nature of the system of algebraic equations that results from spatial discretization of the PDEs. In this manuscript the authors focus on evaluating a proposed nonlinear solution method based on an inexact Newton method with backtracking. In this context they use a particular spatial discretization based on a pressure stabilized Petrov-Galerkin finite element formulation of the low Mach number Navier-Stokes equations with heat and mass transport. The discussion considers computational efficiency, robustness and some implementation issues related to the proposed nonlinear solution scheme. Computational results are presented for several challenging CFD benchmark problems as well as two large scale 3D flow simulations.
The set of nonlinear equations describing the standard model kinematics of the top quark antiquark production system in the dilepton decay channel has at most a fourfold ambiguity due to two not fully reconstructed neutrinos. Its most precise solution is of major importance for measurements of top quark properties like the top quark mass and tt spin correlations. Simple algebraic operations allow one to transform the nonlinear equations into a system of two polynomial equations with two unknowns. These two polynomials of multidegree eight can in turn be analytically reduced to one polynomial with one unknown by means of resultants. The obtained univariate polynomial is of degree 16. The number of its real solutions is determined analytically by means of Sturm's theorem, which is as well used to isolate each real solution into a unique pairwise disjoint interval. The solutions are polished by seeking the sign change of the ...
We consider particle evaporation and fission of an ensemble of hot, rotating compound nuclei as a stochastic process. We derive a set of coupled differential equations formed by a Fokker-Planck equation describing fission, and master equations for calculating particle evaporation. From these equations, we are able to determine multiplicities of prefission neutrons, protons and {alpha}-particles, their energy spectra and their angular momentum distributions. A comparison of our results with experimental data provides us with information regarding the reduced friction coefficient {beta}, the fission barrier height and the level density parameter. For different iridium isotopes, ({sup 181,185,187}Ir), {sup 185}Os and {sup 158}Er, we obtain as an upper limit {beta}{<=}8.0x10{sup 21} s{sup -1}. (orig.).
We consider particle evaporation and fission of an ensemble of hot, rotating compound nuclei as a stochastic process. We derive a set of coupled differential equations formed by a Fokker-Planck equation describing fission, and master equations for calculating particle evaporation. From these equations, we are able to determine multiplicities of prefission neutrons, protons and #alpha#-particles, their energy spectra and their angular momentum distributions. A comparison of our results with experimental data provides us with information regarding the reduced friction coefficient #beta#, the fission barrier height and the level density parameter. For different iridium isotopes, ("1"8"1","1"8"5","1"8"7Ir), "1"8"5Os and "1"5"8Er, we obtain as an upper limit #beta##<=#8.0x10"2"1 s"-"1. (orig.).
A set of mode-coupled excitation equations for the slowly-growing amplitudes of dielectric waveguide eigenmodes is derived as a description of the electromagnetic signal field of a high-gain free-electron laser, or FEL, including the effects of longitudinal space-charge. This approach to describing the field basis set has notable advantages for FEL analysis in providing an efficient characterization of such eigenmodes, and in allowing a clear connection to free-space propagation of the input (seeding) and output radiation. A simple transformation converts the coupled differential excitation equations into a set of coupled algebraic equations and yields a matrix determinant equation for the FEL eigenmodes. A quadratic index medium is used as a model dielectric waveguide to obtain an expression for the predicted spot size of the dominant eigenmode, in the approximation that it consists of a single ...
... 3. Critical constants, and triple point 4. Compressibility isothermE S. Density, molar volume (of a condensed phase) 6. Equations of state, general ...
An improved process is disclosed for reducing the water content of coal containing bound water by releasing at least a portion of the bound water by maintaining the coal at a temperature from about 220/sup 0/ to about 500/sup 0/ F. in the presence of water at a pressure sufficient to maintain at least portion of the water in a liquid phase for a time sufficient to release at least a portion of the bound water wherein the improvement comprises contacting the coal during such treatment with an active material selected from the group consisting of carboxylic organic acids containing up to about 6 carbon atoms, phenol, phenolic acids and inorganic acids.
In the United States, economic growth increasingly requires that greater volumes of freshwater be made available for new users, yet supplies of freshwater are already allocated to existing users. Currently, water for new users is made available through re-allocation of xisting water supplies-for example, by cities purchasing agricultural water rights. Water may also be made available through conservation efforts and, in some locales, through the development of ''new'' water from non-traditional sources such as the oceans, deep aquifer rackish groundwater, and water reuse.
Jet Propulsion Laboratory. California Institute of Technology. Pasadena, California. Energy and Water Conservation. Steve Rigdon. Facilities Energy Manager ...
The water and steam circuits of steam generators in pressurized-water nuclear power plants are described together with the mechanism of denting, and the corrosion of spacer plates that leads to cracks in tubes by constriction. The different chemical specifications applicable to the water of the secondary circuit of the generators in normal operation and on first commissioning are listed. The results obtained and the measurements of chemical values taken in operation on the water in the secondary circuits of steam generators at Fessenheim and Bugey are presented.
The water and steam circuits of steam generators in pressurized-water nuclear power plants are described together with the mechanism of denting, and the corrosion of spacer plates that leads to cracks in tubes by constriction. The different chemical specifications applicable to the water of the secondary circuit of the generators in normal operation and on first commissioning are listed. The results obtained and the measurements of chemical values taken in operation on the water in the secondary circuits of steam generators at Fessenheim and Bugey are presented.
The International Atomic Energy Agency has initiated a co-ordinated research programme on implementation of base-isolation for nuclear structures. This paper discusses two areas relevant to modelling elastomeric base-isolators. These are the use of simplified models to predict the response of isolated structures to earthquake inputs and finite element analysis for calculating the stress distributions within the isolators. In the former, a curvilinear hysteretic model of the high damping natural rubber able to accommodate the stiffening of the rubber at large shear deflections is presented. Its predictions of structural accelerations and bearing displacement produced by design earthquakes and those above the design level are compared with those using a linear spring and dashpot model. A comparison has been made between two finite element analyses using MARC and ABAQUS of the force-deformation behaviour of a single disc of rubber bonded on both sides. The disc was loaded both in ...
Novice problem solvers are rather sensitive to surface problem features, and they often resort to trial and error formula matching rather than identifying an appropriate solution approach. These observations have been interpreted to imply that novices structure their knowledge according to surface features rather than according to problem type categories. However, it may also be the case that novices do know problem types, but cannot map the problem at hand to a known type, because they fail to create a sufficiently well-elaborated problem representation. This study aims to distinguish between these explanations. In this study novice physics students at high and low levels of proficiency completed two problem-sorting tasks from the domain of electricity and magnetism, one with and one without elaboration support. Results confirm that these students do distinguish problem types in accordance with their required solution approaches, and that their problem-sorting ...
This book examines some of the basic principles behind chemical reactions kinetics. Topics considered include ab initio determination of potential energy for chemical reactions; semi-empirical potential energy surfaces; general theory of reactive scattering: different equation approach; and integral equation approach to reactive scattering.
An operator of FE-closure is introduced on the set of functions of a multivalued logic based on the systems of functional equations. It is proved that, for every k ? 2, the FE-closure operator generates a finite classification on the set P k of functions of k-valued logic. The least class in this classification is shown to be the class H k of all homogeneous functions. Also a series of corollaries are obtained concerning the finite FE-generating sets in the FE-closed classes.
The pressure of a 6:12 fluid is obtained from the Percus-Yevick theory at the temperature kT/? = 1.35. The energy equation results are in excellent agreement with machine calculation results in contrast to the pressure and compressibility equation results. Present address: Diffusion Research Unit, The Australian National University, Canberra, ACT, Australia.
The spin-up problem of He II in a cylindrical container of infinite and finite height is considered on the basis of the Hall-Vinen-Bekarevich-Khalatnikov equations. The boundary conditions include the assumption of the generation of superfluid vorticity at the walls of the container. In the lowest order approximation with respect to viscosity, the problem is reduced to a Stefan problem for a single diffusion-type equation. Examples of solutions are given which show that the method used gives a good insight into the physical development of the spin-up flow.
An explicite PN solution of the multi-dimensional homogeneous neutron transport equation is given by expanding the angular flux into a series of geometry-independent spherical harmonics operators. An algorithm is developed for representing the spherical harmonic operators in orthogonal curvilinear coordinates. The general formulae are applied to two-dimensional spherical geometry; detailed P3 formulae are given. (orig.).
For a system of noninteracting electrons scattered by resonant levels of neutral impurities, we show that virial and quasiparticle corrections have nearly equal magnitudes. We propose a modification of the Boltzmann equation that includes quasiparticle and virial corrections and discuss their interplay on a dielectric function. copyright 1997 The American Physical Society.
The wave equation of a spinless tachyon is studied in Schwarzschild space-time. In contrast to earlier approaches to the problem, it is shown that tachyonic static solutions satisfy a simple second-order linear differential equation regardless of the mass of the black hole and the mass parameter of the tachyon. Physical implication of the present approach is discussed. Using Langer modification of the WKB (Wentzel-Kramers-Brillouin) boundary condition an expression similar to the Bohr-Sommerfeld quantization condition is derived.
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.
We study the dynamics of the scalar tensor cosmological model in the presence of tachyon field. In an alternative approach, in two exponential and power law form of the scalar field functions in the model, field equations are solved by simultaneously best fitting the model parameters with the most recent observational data. This approach gives us an observationally verified interpretation of the dynamics of the universe. We then discuss the best fitted of equation of state parameter, the statefinder parameters and the reconstructed scalar field in the model.
Diffraction by an impermeable scatterer in ?3 with a Dirichlet boundary condition is considered. It is shown that the null field equation is equivalent to the original diffraction problem. The general theoretical constructions are illustrated by analyzing the excitation of a sphere by a point source. The evolution of the residual of the null field equation as the scatterer surface is approached by an auxiliary surface is estimated.
A method of auxiliary spectrum is modified so that matrix elements of the reaction matrix in final nuclei are determined by means of algebraic operations alone (inversion of matrices). No differential equations need to be solved; Pauli's exclusion principle is accurately taken into account. A single-particle potential may be of any kind, but a two-particle interaction must have no solid core.
1 - Description of program or function: This program solves an eigenvalue problem zBx=Ax where A and B are real block tri-diagonal matrices. This eigenvalue problem is derived from a reduced set of linear resistive MHD equations which is often employed to study tokamak plasma stability problem. 2 - Method of solution: Both the determinant and inverse iteration methods are employed. 3 - Restrictions on the complexity of the problem: The eigenvalue z must be real
The loop quantum cosmology 'improved dynamics' of the Bianchi type IX model are studied. The action of the Hamiltonian constraint operator is obtained via techniques developed for the Bianchi type I and type II models, no new input is required. It is shown that the big bang and big crunch singularities are resolved by quantum gravity effects. We also present effective equations which provide quantum geometry corrections to the classical equations of motion.
By a variant version of Fountain Theorem due to Zou [W. Zou, Variant fountain theorems and their applications, Manuscripta Math. 104 (2001) 343-358], the existence of infinitely many solutions is obtained for a class of biharmonic equations where the nonlinearity involves a combination of superlinear and asymptotically linear terms.
The optimization of ethylene furnace operations using cfd-based simulations has been addressed. The optimization problems have been cast into various formulations: the Multidisciplinary feasible (MDF) approach, the All-At-Once (AAO) approach and the Individual discipline feasible (IDF) approach. These approaches mainly differ in their handling of the state equations as constraints, and hence some of the formulations place restrictions on the methods used to solve the state equations.
Zanker has recently presented nomographs for determining particle settling velocities in solid-liquid systems. These nomographs were based on the general correlations developed by Barnea and Mizrahi and Barnea and Mednick. This work presents an equation directly computing particle settling velocities, eliminating the uncertainty associated with nomographs.
In this paper we first establish global pointwise time-space estimates for a class of oscillatory integrals. Then, we use them to establish $L^p-L^q$ estimates for a class of higher order wave-type equations of the form $\\partial_{tt}u+P(D_{x})u=0$, where the symbol $P(\\xi)$ is a real non-degenerate elliptic polynomial of ${\\bf R}^n$.
This demonstration project uses modern simulation techniques to illustrate the important technologies and design variables that an auto-designer would consider in production a high efficiency, low emissions vehicle. Simulation and modeling techniques use the idea of capturing the relationships between real components of the systems with mathematical equations. These equations are then solved on a computer to simulate the behavior or performance of the system under various conditions. In the current demonstration project, we focus on many variations of a hydrogen-powered vehicle.
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.
Four solutions of the Cauchy problem for Mathieu?s equation away from parametric resonance domains are analytically constructed using an asymptotic averaging method in the fourth approximation. Three solutions occur near fractional parameter values at which slow combination phases exist. The fourth solution occurs in the absence of slow phases away from parametric resonance domains and the fractional parameter values.
Equations are compiled for thermal balance in which for simplification, no consideration is made for heat conductivity along the axis of the cable and dependence of losses, heat capacitance and heat conductivity on temperature. Equations are modeled on a transistor analog calculator 42 TA. The solution to the task on the computer produced values of maximum temperature on the cable and coordinates of the point of maximum overheating. Using the analog model, one can study other parameters of the thermal mode.
To calculate the energy band structures in semiconductors using the relativistic augmented plane wave method, atomic potential and charge density are needed, which are calculated by self-consistent method. Wave function for one electron is determined by solving the Dirac equation with the Hartree-Fock equation based on the slater's exchange potential. The results of calculation for Cu"+"1 are given. (Author).
AimsTwo-dimensional echocardiographic (2DE) continuity-equation derived aortic valve area (AVA) in aortic stenosis (AS) relies on non-simultaneous measurement of left ventricular...Full Text Available
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)
A formula for calculating the Dirac S-matrix for central Lorentz scalar and vector potentials is derived by use of a new amplitude-phase method. The derivation also makes use of certain invariants of the reduced 2-spinor radial Dirac equations.
A formula for calculating the Dirac S-matrix for central Lorentz scalar and vector potentials is derived by use of a new amplitude-phase method. The derivation also makes use of certain invariants of the reduced 2-spinor radial Dirac equations.
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.
We consider the problem of rotation in a homogeneous Bianchi type IX cosmological model. Studying the Wheeler-De Witt equation corresponding to this minisuperspace, and adopting a particular choice of the factor ordering, we are able to find a particular solution which is strongly peaked about isotropy. This result confirms all the previous investigations in this field, and suggests the conclusion that the machian nature of our universe is the natural, i.e. most likely, outcome of the Planck epoch. (orig.).
We consider the problem of rotation in a homogeneous Bianchi type IX cosmological model. Studying the Wheeler-De Witt equation corresponding to this minisuperspace, and adopting a particular choice of the factor ordering, we are able to find a particular solution which is strongly peaked about isotropy. This result confirms all the previous investigations in this field, and suggests the conclusion that the machian nature of our universe is the natural, i.e. most likely, outcome of the Planck epoch. (orig.).
We treat symmetric semi-infinite nuclear matter in the relativistic mean-field approximation for the scalar-vector field theoretical model. Using special-type Dirac spinors the nucleonic Dirac equation is decoupled into two sets of differential equations for the spin-orientation dependent orbital nucleon Dirac spinors. We also rewrite the Dirac equation in terms of second-order differential equations with the spin-orbit interaction appearing explicitly. These equations can be solved if the spin-orbit part is left out. The spin-orbit effects thus can be isolated, and are shown to reduce the surface energy coefficient a/sub s/ as well as the surface thickness t in such a way that their ratio remains practically unchanged. For realistic lagrangians - in linear as well as non-linear forms - consistent with the empirical spin-orbit single-partial level splittings, a/sub s/ is reduced by ...
The Schroedinger equation for nuclear quadruple surface vibrations is linearized with the consequence that a new spin degree of freedom appears in the wave function of the linearized equation. This spin is called collective spin and has a value of 3/2. The linearized Schroedinger equation for quadrupole vibrations is used for the description of certain collective aspects of even-odd {sup 187,189,191}Ir nuclei which have a spin 3/2 in their ground state. As a potential we use the {gamma}-soft collective potential of the neighboring even-even nuclei, which is inserted into the linearized Schroedinger equation via a scalar coupling. This leads to a collective spin-dependent fine structure splitting of the energy levels governed by a collective SO(5) spin-orbit coupling and a correction to the kinetic energy. Further, we consider explicitly spin-dependent potentials which effectively describe the ...
The Schroedinger equation for nuclear quadruple surface vibrations is linearized with the consequence that a new spin degree of freedom appears in the wave function of the linearized equation. This spin is called collective spin and has a value of 3/2. The linearized Schroedinger equation for quadrupole vibrations is used for the description of certain collective aspects of even-odd "1"8"7","1"8"9","1"9"1Ir nuclei which have a spin 3/2 in their ground state. As a potential we use the #gamma#-soft collective potential of the neighboring even-even nuclei, which is inserted into the linearized Schroedinger equation via a scalar coupling. This leads to a collective spin-dependent fine structure splitting of the energy levels governed by a collective SO(5) spin-orbit coupling and a correction to the kinetic energy. Further, we consider explicitly spin-dependent potentials which effectively describe the ...
We analyze the stability of the structure equations of the vacuum in the brane world models, by using both the linear (Lyapunov) stability analysis, and the Jacobi stability analysis, the Kosambi-Cartan-Chern theory. In the brane world models the four-dimensional effective Einstein equations acquire extra terms, called dark radiation and dark pressure, respectively, which arise from the embedding of the three-brane in the bulk. Generally, the spherically symmetric vacuum solutions of the brane gravitational field equations have properties quite distinct as compared to the standard black hole solutions of general relativity. We close the structure equations by assuming a simple linear equation of state for the dark pressure. In this case the vacuum is Jacobi stable only for a small range of values of the proportionality constant relating the dark pressure and the dark radiation. The ...
This paper develops a clear procedure for solving the nonlinear Vlasov-Maxwell equations for a one-component intense charged particle beam or finite-length charge bunch propagating through a cylindrical conducting pipe (radius r = r(subscript)w = const.), and confined by an applied focusing force. In particular, the nonlinear Vlasov-Maxwell equations are Lorentz-transformed to the beam frame ('primed' variables) moving with axial velocity relative to the laboratory. In the beam frame, the particle motions are nonrelativistic for the applications of practical interest, already a major simplification. Then, in the beam frame, we make the electrostatic approximation which fully incorporates beam space-charge effects, but neglects any fast electromagnetic processes with transverse polarization (e.g., light waves). The resulting Vlasov-Maxwell equations are then Lorentz-transformed back to the laboratory frame, ...
A linear structural equation model relates random variables of interest and corresponding Gaussian noise terms via a linear equation system. Each such model can be represented by a mixed graph in which directed edges encode the linear equations, and bidirected edges indicate possible correlations among noise terms. We study parameter identifiability in these models, that is, we ask for conditions that ensure that the edge coefficients and correlations appearing in a linear structural equation model can be uniquely recovered from the covariance matrix of the associated normal distribution. We treat the case of generic identifiability, where unique recovery is possible for almost every choice of parameters. We give a new graphical criterion that is sufficient for generic identifiability. It improves criteria from prior work and does not require the directed part of the graph to be acyclic. We also develop ...
For good image quality using ultrasound inverse scattering, one alternately solves the well-posed forward scattering equation for an estimated total field and the ill-posed inverse scattering equation for the desired object property function. In estimating the total field, error or noise contaminates the coefficients of both matrix and data of the inverse scattering equation. Previous work on ill-posed inverse ultrasonic scattering commonly used Tikhonov regularization, which considers error only in the data. The solution so obtained is not precise enough to reconstruct the quantitative internal structure of a large or high-contrast object. This paper adopts the truncated total least squares method, simultaneously considering error and noise on both sides of the inverse scattering equation, and compares it with the classical Tikhonov regularization. We show that it can substantially improve ...
An efficient method of hydrogen storage in nano-porous carbons is its reversible sorption by electrochemical decomposition of a KOH water solution according to the following equation: C + xH{sub 2}O + xe{sup -} {yields} (CH{sub x}) + xOH{sup -} where (CH{sub x}) stands for the hydrogen inserted into the nano-porous carbon during charging and oxidized during discharging. Although various carbon materials have been investigated as hydrogen adsorbents, the information about the storage mechanism as well as the nature of the hydrogen/carbon interaction is still not sufficient. In order to extend the understanding of the process, carbon samples charged electrochemically were investigated by temperature programmed desorption (TPD). The nature of the hydrogen/carbon interaction was studied by electrochemical analysis at different temperatures. The TPD experiments consist of heating the samples from room temperature to 950 C and of quantitative ...
As nuclear research has become more prevalent, environmental contamination from the disposal of radioactive waste has become a prominent issue. At Los Alamos National Laboratory (LANL) in northern New Mexico, radioactive contamination from disposal operations has raised some very specific concerns. Material Disposal Area G (Area G) is the primary low-level radioactive waste disposal site at LANL and occupies an area adjacent to land belonging to the Native American community of the Pueblo of San Ildefonso. Analyses of soil and vegetation collected from the perimeter of Area G have shown concentrations of radionuclides greater than background concentrations established for northern New Mexico. As a result, Pueblo residents had become concerned that contaminants from Area G could enter tribal lands through various ecological pathways. The residents specifically questioned the safety of consuming meat from elk and deer that forage near Area G and then migrate onto tribal lands. ...
Cooling water is sprayed on the outer surface of an upper portion of a container, and a pool is formed by the cooling water flowing down while cooling the container. Further, the cooling water stored in the cooling water pool is recycled by a pump for spraying the cooling water on the outer surface of the upper portion of the container. Sufficient amount of cooling water is supplied for spraying the cooling water to the outer surface of the upper portion of the container so that the outer surface of the container is free from drying and a liquid membrane is formed on the entire surface. The amount of the cooling water is made greater than that of the cooling water evaporated when the entire amount of the heat generate in the reactor core of the reactor is transferred to the cooling ...
Hydrothermal oxidation (HTO) is a promising technology for the treatment of aqueous-fluid hazardous and mixed waste streams. Waste streams identified as likely candidates for treatment by this technology are primarily aqueous fluids containing hazardous organic compounds, and often containing inorganic compounds including radioisotopes (mixed wastes). These wastes are difficult and expensive to treat by conventional technologies (e.g. incineration) due to their high water content; in addition, incineration can lead to concerns related to stack releases. An especially attractive potential advantage of HTO over conventional treatment methods is the total containment of all reaction products within the overall system. The potential application of hydrothermal oxidation (HTO) technology for the treatment of DOE hazardous or mixed wastes has been uncertain due to concerns about safe and efficient operation of the technology. In principle, aqueous DOE wastes, including ...
Water-resources data for the 1993 water year for Texas are presented in four volumes, and consist of records of stage, discharge, and water quality of streams and canals; and stage, contents, and water quality of lakes and reservoirs; and water levels and water quality of ground-water wells. Volume 3 contains records for water discharge at 134 gaging stations; stage only at 1 gaging stations; stage and contents at 13 lakes and reservoirs; water quality at 81 gaging stations; and data for 30 partial-record and 4 flood-hydrograph partial-record stations. Also included are lists of discontinued surface-water discharge or stage-only stations and discontinued surface-water-quality stations; crest-stage and flood-hydrograph partial-record stations, reconnaissance ...
Singlet oxygen generators are multiphase flow chemical reactors used to generate energetic oxygen to be used as a fuel for chemical oxygen iodine lasers. In this paper, a theoretical model of the generator is presented along with its solutions over ranges of parameter space and oxygen maximizing optimizations. The singlet oxygen generator (SOG) is a low-pressure, multiphase flow chemical reactor that is used to produce molecular oxygen in an electronically excited state, i.e. singlet delta oxygen. The primary product of the reactor, the energetic oxygen, is used in a stage immediately succeeding the SOG to dissociate and energize iodine. The gas mixture including the iodine is accelerated to a supersonic speed and lased. Thus the SOG is the fuel generator for the chemical oxygen iodine laser (COIL). The COIL has important application for both military purposes--it was developed by the US Air Force in the 1970s--and, as the infrared beam is readily absorbed by metals, industrial cutting ...
The chemistries of reactants, plasticizers, solvents and additives in an epoxy paint are discussed. Polyamide additives may play an important role in the absorption of molecular iodine by epoxy paints. It is recommended that the unsaturation of the polyamide additive in the epoxy cure be determined. Experimental studies of water absorption by epoxy resins are discussed. These studies show that absorption can disrupt hydrogen bonds among segments of the polymers and cause swelling of the polymer. The water absorption increases the diffusion coefficient of water within the polymer. Permanent damage to the polymer can result if water causes hydrolysis of ether linkages. Water desorption studies are recommended to ascertain how water absorption affects epoxy paint.
...Species - see CITES Trade of Dangerous Chemicals Transport and Environment U Urban Environment Use of natural resources V Volatile Organic Compounds W Waste Landfills Waste electrical and electronic equipment - see WEEE Waste incineration Waste legislation - reporting Waste management planning Waste oils Waste shipments Waste water - see Urban waste waterWater Bathing Water Drinking water Floods Marine environment Urban waste waterWater Framework Directive WEEE - Waste Electric and Electronic Equipment Wildlife ...
A copyrighted computational fluid dynamics computer code, ICOMFLO, has been developed for the simulation of multiphase reacting flows. The code solves conservation equations for gaseous species and droplets (or solid particles) of various sizes. General conservation laws, expressed by elliptic type partial differential equations, are used in conjunction with rate equations governing the mass, momentum, enthalpy, species, turbulent kinetic energy, and turbulent dissipation. Associated phenomenological submodels of the code include integral combustion, two parameter turbulence, particle evaporation, and interfacial submodels. A newly developed integral combustion submodel replacing an Arrhenius type differential reaction submodel has been implemented to improve numerical convergence and enhance numerical stability. A two parameter turbulence submodel is modified for both gas and solid phases. An evaporation submodel treats ...
This report is a summary of the contribution of nitrogen and phosphorus in the eutrophication process of inland and coastal waters. Special attention was paid to the mechanisms of these nutrients in regulating biological processes and to the methods available in estimating their effects in the eutrophication of water bodies. The report includes five chapters which are entitled: Introduction, which is a general background to the subject with special attention to the requirements of the Finnish Water Act. Phosphorus and nitrogen as factors regulating biological processes. The topics included are: definition of eutrophication, forms of phosphorus and nitrogen and their sources to inland and coastal waters, effects of these nutrients as growth factors of phytoplankton and macrophytes and consequences of eutrophication. Estimation of the effects of phosphorus and nitrogen. The topics discussed from the point ...
This procedure describes the formulation and make-up of Simulated Dilute Water (SOW), a low-ionic-content water to be used for Activity E-20-50, Long-Term Corrosion Studies. This water has an ionic content which is nominally a factor of ten higher than that of representative waters at or near Yucca Mountain. Representative waters were chosen as J-13 well water [Harrar, 1990] and perched water at Yucca Mountain [Glassley, 1996]. J-13 well water is obtained from ground water that is in contact with the Topopah Spring tuff, which is the repository horizon rock. The perched water is located in the Topopah Spring tuff, but below the repository horizon and above the water table. A nominal times ten higher ionic content was chosen to simulate the effect of ionic ...
The study investigates the possibility of enhancing crop water productivity in the parts of Northwest India where groundwater quality is marginal and canal water supply is severely scarce. Soil, Water, Atmosphere and Plant (SWAP) model was calibrated and validated in three farmers' fields with varying canal water availability and groundwater quality in the Kaithal Irrigation Circle of the Bhakra Canal system, Haryana. On the basis of predicted and observed soil water content, pressure heads, salt concentration at 2 week intervals and crop yields, the model was found suitable for use in the region. A few nomographs were prepared to provide a graphical method to predict the effect of different combinations of water quality and depth of water application on crop yield and soil salinity and to...
The research committee of the Atomic Energy Society of Japan on water chemistry standard aims at establishing the private standard of water chemistry of nuclear power plants. The committee gathers up 'BWR water chemistry management manual', 'PWR primary system water chemistry management manual' and 'PWR water chemical analysis standard method', and furthermore aims at the standardization of those in future. Looking back on the committee's activities for the past four years, latest results of research of water chemistry mainly contributing to safe and reliable nuclear power plants were described with the future perspective of water chemistry and a demanded break-through. (T.T.)
Using the Laser-Fluorescence technique, an investigation was made, adopting the standard mix method, on trace uranium concentrations in hot spring water and geothermal water from Lishan region, and in tap water from some major cities in Shanxi province. Totally 40 samples from 27 sites were investigated. Measurement showed that the tap water contains around 10"-"6 g/L of uranium, whose concentrations in both hot spring water and geothermal water are 10"-"5 g/L. Most of samples are at normal radioactive background level, some higher contents were determined in a few samples
In equipment industries, the equipments handling industrial water and pure water are numerous. In power generation including nuclear power generation, water serves as a working medium. Review is made on the experiences in the corrosion of iron and nickel base alloys in high temperature, high pressure water and the results of researches derived from them. Under high temperature and high pressure, carbon steel, low alloy steel, stainless steel and high nickel alloy cause corrosion even in pure water. But in the case of serious corrosion, chlorine, oxygen, alkali and others in water take part. The following matters are described: corrosion by steam; stress corrosion cracking in pure water; corrosion by impurities in high temperature, high pressure pure water, i.e. chlorine ions, dissolved oxygen, and alkali; corrosion under ...
...AquaMagna Water Conditioner and Descaler OCETA Environmental Technology & Business Profiles AQUAMAGNA WATER CONDITIONER AND DESCALER CAT. #03-018 CONTACT: Magna-Tek ...removal high-strength ceramic magnet no energy input The AquaMagna Water Conditioner and Descaler is a non-chemical method of water treatment. The technology uses ... Figure 1: The AquaMagna Water Conditioner and Descaler PROCESS/ PRODUCT APPLICATION: hard water domestic water process water reverse osmosis heat exchangers cooling ... For once-through or recirculating heating or cooling equipment, the AquaMagna descaler improves thermal efficiencies and reduces downtime. Figure 2 illustrates a cross-section ...
This paper presents a novel domestic hot water heater model to be used in a multi-objective demand side management program. The model incorporates both the thermal losses and the water usage to determine the temperature of the water in the tank. Water heater loads are extracted from household load data and then used to determine the household water usage patterns. The benefits of the model are: (1) the on/off state of the water heater and temperature of the water in the tank can be accurately predicted, and (2) it enables the development of water usage profiles so that users can be classified based on usage behaviour. As a result, the amount of ancillary services and peak shaving that can be achieved are accurately predictable and can be maximized without adversely affecting users. (author)
...Tota-MaharajE-Mail: Interests: water and wastewater treatment, environmental engineering and sustainable systems; sustainable water management; sustainable urban drainage systems (SUDS); combined renewable energy applications with reverse osmosis desalination; microbial fuel cells for bioenergy production and treatment of urban wastewater; solar photocatalytic treatment and disinfection of water/wastewater Dr. Simon Toze CSIRO Land and Water, Queensland Bioscience Precinct - St Lucia, 306 Carmody Road, St Lucia QLD 4067,...
Indicators of source water use in different economic branches, water consumption in recycling and recycling-successive water supply in 1980?2006 are analyzed. Trends were identified in variations in specific water consumption in industry, agriculture, and municipal economy. Territorial differentiation of wetness indicators of regional product over federal districts and constituent territories of the RF is shown. Wetness values of gross domestic product in different economies are compared.
A wavelet transformation is performed over each of the spatial coordinates of the scalar wave equation. This transformed equation is solved directly with a finite-difference scheme for both homogeneous and smooth inhomogeneous media. Wavefield extrapolation is performed completely in the spatial wavelet domain without transforming back into the space domain at each time step. The wavelet coefficients are extrapolated, rather than the wavefield itself. The numerical solution of the scalar wave equation in the spatial wavelet domain is closely related to the finite-difference method because of the compact support of the wavelet bases. Poststack reverse-time migration is implemented as an application. The resolution spaces of the wavelet transform provide a natural framework for multigrid analysis. Migrated images are constructed from various resolution spaces.
VSP reverse-time migration is a well adaptable wave equation migration method. Its control equation not only describes all-direction propagation of seismic wave but also removes interbed multiples. Clearbout's image principle is generalized to determine image conditions, real VSP data are used to determine boundary condition, and two way reflection-free wave equation is solved by making reverse-time extrapolation. In each step of extrapolation, the migration value at relevant image point is obtained by using the image condition. The complete migration of a seismic section is achieved when reverse-time extrapolation reaches the minimum image time. In this paper it is proved theoretically and practically that this method is applicable to any velocity variation and makes the migrated section have both good resolution and high S/N ratio. Besides, this method results in high processing efficiency.
Heat capacities and densities of aqueous MgCl/sub 2/, CaCl/sub 2/, and SrCl/sub 2/ from the accompanying paper are combined with literature data up to 473 K to yield temperature-dependent equations by using the ion-interaction model of Pitzer. These heat capacity equations have been integrated to yield the enthalpy and the Gibbs energy. The enthalpy parameters for 298 K are evaluated in separate calculations using published high-temperature osmotic data as well as heats of dilution, while the Gibbs energy parameters for 298 K are taken from the literature. The range of validity of the final equations is described.
A new formulation of the theory of tachyons using the same two postulates as in special relativity is applied to the electrodynamics of material media. A discussion of Lagrange`s equations and Hamilton`s equations for `classical` charged tachyons shows that such a formalism is a viable approach. An essay is included on why tachyons can be considered to be localised particles for the purpose of calculations. Tachyonic transformations of the electromagnetic fields D, P, H and M are shown to be the same as for bradyonic transformations. Examples discussed include the electric dipole moment of a tachyonic current loop, constitutive equations, polarisation in tachyonic dielectric materials and the velocity of light in tachyonic dielectric media. This is followed by discussions of the collision energy loss for charged tachyons interacting with a material medium and a mathematical proof that tachyons cannot emit Cherenkov ...
A new formulation of the theory of tachyons using the same two postulates as in special relativity is applied to the electrodynamics of material media. A discussion of Lagrange's equations and Hamilton's equations for 'classical' charged tachyons shows that such a formalism is a viable approach. An essay is included on why tachyons can be considered to be localised particles for the purpose of calculations. Tachyonic transformations of the electromagnetic fields D, P, H and M are shown to be the same as for bradyonic transformations. Examples discussed include the electric dipole moment of a tachyonic current loop, constitutive equations, polarisation in tachyonic dielectric materials and the velocity of light in tachyonic dielectric media. This is followed by discussions of the collision energy loss for charged tachyons interacting with a material medium and a mathematical proof that tachyons cannot emit Cherenkov ...
The space derivation term of a wave equation is accurately calculated using Fourier transform method, and the wave equation only relating to time derivative in time-space domain is derived. Solving this equation with the aid of central difference method brings the numerical recursion formula for forward modeling or reverse time migration of seismic wave field in the medium in which seismic velocity varies. The key to this method lies in introducing two auxiliary wave fields relating to velocity and spacial frequency of Fourier transform respectively after multi-dimensional Fourier transform respectively after multi-dimensional Fourier transform of space vector is made. Theoretically, this method is suitable to the forward modeling and migration of seismic wave field in complex area where seismic velocity and structure shape vary arbitrarily. Theoretical results proved this method satisfactory.
We investigate the existence of time-periodic solutions of the Dirac equation in the Kerr-Newman background metric. To this end, the solutions are expanded in a Fourier series with respect to the time variable $t$ and the Chandrasekhar separation ansatz is applied so that the question of existence of a time-periodic solution is reduced to the solvability of a certain coupled system of ordinary differential equations. First, we prove the already known result that there are no time-periodic solutions in the non-extreme case. Then it is shown that in the extreme case for fixed black hole data there is a sequence of particle masses $(m_N)_{N\\in\\mathbb N}$ for which a time-periodic solution of the Dirac equation does exist. The period of the solution depends only on the data of the black hole described by the Kerr-Newman metric.
We study perturbations of a Schwarzschild black hole in Chern-Simons modified gravity. We begin by showing that Birkhoff's theorem holds for a wide family of Chern-Simons coupling functions, a scalar field present in the theory that controls the strength of the Chern-Simons correction to the Einstein-Hilbert action. After decomposing the perturbations in spherical harmonics, we study the linearized modified field equations and find that axial and polar modes are coupled, in contrast to general relativity. The divergence of the modified equations leads to the Pontryagin constraint, which forces the vanishing of the Cunningham-Price-Moncrief master function associated with axial modes. We analyze the structure of these equations and find that the appearance of the Pontryagin constraint yields an overconstrained system that does not allow for generic black hole oscillations. We illustrate this situation by studying the case ...
A macroscopic realization of the strange virtual particles is presented. The classical Helmholtz and the quantum mechanical Schr\\"odinger equations are analogous differential equations. Their imaginary solutions are called evanescent modes in the case of elastic and electromagnetic fields. In the case of non-relativistic quantum mechanical fields they are called tunneling solutions. The imaginary solutions of this differential equation point to strange consequences: They are non local, they are not observable, and they described as virtual particles. During the last two decades QED calculations of the imaginary solutions have been experimentally confirmed for phonons, photons, and for electrons. The experimental proofs of the predictions of the non-relativistic quantum mechanics and of the Wigner phase time approach for the elastic, the electromagnetic and the Schr\\"odinger fields will be presented in this article. The ...
Equilibrium and energy equations have been developed in describing the solid-fluid transition zone during the melting and solidification of the binary alloys. Due to the existence of the transition region exhibited both solid and fluid phases at the same material point from continuum point of view, mixture theory was utilized to analyze the region. Unlike the Stefan problem, the latent heat due to the phase change appears as a source term in the heat equation. The molten fluid is treated as a thermoviscous and incompressible fluid, whereas the solid is thermoviscoplastic described by the Bodner-Partom/Walker type of constitutive equations. Thermal mechanical behaviors of the solid and the fluid phases are determined separately because of insignificant mechanical interactions between them. Volume fractions of the phases are obtained according to the equilibrium phase diagram. The simulation process of the transition zone and ...
We obtain x-ray absorption near-edge structures (XANES) by solving the equation of motion for the two-particle Green's function for the electron-hole pair, the Bethe-Salpeter equation (BSE), within the all-electron full-potential linearized augmented plane wave method (FPLAPW). The excited states are calculated for the Li K-edge in the insulating solids LiF, Li_2O and Li_2S, and absorption spectra are compared with independent particle results using the random phase approximation (RPA), as well as supercell calculations using the core-hole approximation within density functional theory (DFT). The binding energies of strongly bound excitations are determined in the materials, and core-exciton wavefunctions are demonstrated for LiF.
Lateral plume dispersion in deep river valleys during neutral and stable lapse rate conditions can be exceptionally high due to the intense horizontal turbulence generated by prominent variations in the width, orientation and surface roughness of the valley walls. Use of the standard Gaussian dispersion equation to represent this enhanced lateral plume spread in a narrow valley may be inappropriate since consideration should also be given to the restriction of the horizontal plume spread due to impingement against the valley walls. The basic concept employed in the modification of the Gaussian dispersion equation was to assume that multiple eddy reflections occur between the valley walls in a manner similar to the vertical eddy reflections between the ground plane and an inversion layer aloft. The expressions are developed, but no attempt has been made to validate them with actual measured field data. 4 references, 3 figures, 1 table.
This paper presents research work concerning the mathematical modelling of air pollution in the Athens basin. A mathematical model based on aerodynamics and conservation equations together with turbulence modelling in body-fitted coordinates is used to provide the necessary input information to a similar but more detailed (i.e. small scale) model around Athens. Both models consist of solving the full Navier-Stokes equations and take into account the precise topography of the region. Different roughness heights are used for sea, coastal regions and land. The smaller scale model makes use of a two-equation turbulence model ({kappa}-{epsilon}) and the large scale model assumes a constant turbulence viscosity. The pollutant sources used in these models are based on actual emissions in the Athens area and are distributed in space. The results of this work appear promising and demonstrate the capability of the model as a tool for ...
Within the framework of the Bianchi type-IX homogeneous space, we set up a system of coupled equations for the cosmic scale factors, scalar field, and radiative energy density. At the tree level, the equations are written in a self-consistent, Hartree-Fock form. For phi/sup 4/ theory, the system of nine first-order differential equations is solved numerically for a varying ratio of the energy of anisotropy to the vacuum energy. As the vacuum energy increases, there appears to be less reheating, since the energy of anisotropy is more efficiently converted into isotropic expansion. If the energy of anisotropy is large enough, the inflationary phase is prevented. In this case, a series of cosmological phase transitions will take place each time the square of the effective mass changes its sign.
Within the framework of the Bianchi type-IX homogeneous space, we set up a system of coupled equations for the cosmic scale factors, scalar field, and radiative energy density. At the tree level, the equations are written in a self-consistent, Hartree-Fock form. For phi"4 theory, the system of nine first-order differential equations is solved numerically for a varying ratio of the energy of anisotropy to the vacuum energy. As the vacuum energy increases, there appears to be less reheating, since the energy of anisotropy is more efficiently converted into isotropic expansion. If the energy of anisotropy is large enough, the inflationary phase is prevented. In this case, a series of cosmological phase transitions will take place each time the square of the effective mass changes its sign.
Many methods have been proposed for the determination of the corrosion current density from polarization curves; the two most commonly used ones are the three-point method and the curve-fitting method. The basis of these methods is a polarization equation that was derived for the complete absence of mass-transport effect or considers the cathodic partial reaction to be under complete mass-transport control. The mixed-control case has been much less investigated, even though this case may quite frequently reflect practical situations. A completely generalized polarization equation is proposed that can be used as a basis for curve-fitting data evaluation under any conditions. This equation was used to determine the error caused by the neglect of mass-transport effect in conventional data evaluation.
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)
The present investigation is concerned with the effect of Hall currents on boundary layer flow, and heat and mass transfer of an electrically conducting fluid over an unsteady stretching sheet in the presence of a strong magnetic field. The electron-atom collision frequency is assumed to be relatively high, so that the Hall effect is assumed to exist, while the induced magnetic field is neglected. The governing time-dependent boundary layer equations for momentum, thermal energy, and concentration are reduced using a similarity transformation to a set of coupled ordinary differential equations. The similarity ordinary differential equations are then solved numerically by the successive linearization method together with the Chebyshev pseudo-spectral collocation method. Effects of the Prand...
We study the thermodynamics of a one-dimensional attractive Fermi gas (the Gaudin-Yang model) with spin imbalance. The exact solution has been known from the thermodynamic Bethe ansatz for decades, but it involves an infinite number of coupled nonlinear integral equations whose physics is difficult to extract. Here the solution is analytically reduced to a simple, powerful set of four algebraic equations. The simplified equations become universal and exact in the experimental regime of strong interaction and relatively low temperature. Using the new formulation, we discuss the qualitative features of finite-temperature crossover and make quantitative predictions on the density profiles in traps. We propose a practical two-stage scheme to achieve accurate thermometry for a trapped spin-imbalanced Fermi gas.
The equilibrium composition and temperature resulting from the combustion of hydrocarbon fuels in combination with hydrogen have been analytically evaluated. Rich, lean and stoichiometric mixtures have been considered for this purpose. The flame temperature data obtained from computer calculations has been fitted with an equation expressed in terms of calorific value, C/H ratio, air-fuel ratio and heat of dissociation of the fuel mixture. This equation can be used to calculate the flame temperature directly without resorting to the elaborate calculation of species concentration. The results obtained using the equation have been found to have an accuracy of 0.05 to 2.5% for different fuel combinations. (author).
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...
We obtain an elegant and useful description of the dynamics of Szekeres dust models (in their full generality) by means of "quasi--local" scalar variables constructed by suitable integral distributions that can be interpreted as weighed proper volume averages of the local covariant scalars. In terms of these variables, the field equations and basic physical and geometric quantities are formally identical to their corresponding expressions in the spherically symmetric LTB dust models. Since we can map every Szekeres model to a unique LTB model, rigorous results valid for the latter models can be readily generalized to a non--spherical Szekeres geometry. The new variables lead naturally to an initial value formulation in which all scalars are expressed as scaling laws in terms of their values at an arbitrary initial space slice. These variables also yield a significant simplification of numerical work, since the fluid flow evolution equations ...
A nodal integral method is derived for the monoenergetic, steady-state, fixed source neutron diffusion equation in hexagonal geometry based on a coordinate transformation that maps a parallelogram into a rectangle. The new hexagonal nodal diffusion method is implemented in the computer code HND where the discrete-variable equations are solved via an iterative scheme. Because the new method's equations are derived for a rhombus, they can be solved on a sequence of embedded meshes to study the method's error order. Indeed a preliminary numerical error analysis reveals a second-order error in the mesh size, and comparison with finite difference results obtained with the finite difference based BOLD-VENTURE code indicate the superior accuracy of our new nodal method.
Self-excited oscillations of a distributed medium can be investigated in terms of discrete physical models described by multidimensional systems of nonlinear equations which transform into equations of self-oscillating distributed media as the number of dimensions goes to infinity. Differential equations of such a self-oscillating system are set up and investigated. Coupling between elements of the medium, conditions under which the system constitutes an open system, internal impedance between the energy source and the system elements, and boundary conditions defining open and closed systems are discussed. The Gunn effect, eolian noise, self-oscillations of DNA molecules and of cardiac muscle, and wave processes in the cerebral cortex are mentioned as examples of self-oscillatory processes in question.
This contribution deals with identification of fractional-order dynamical systems. System identification, which refers to estimation of process parameters, is a necessity in control theory. Real processes are usually of fractional order as opposed to the ideal integral order models. A simple and elegant scheme of estimating the parameters for such a fractional order process is proposed. This method employs fractional calculus theory to find equations relating the parameters that are to be estimated, and then estimates the process parameters after solving the simultaneous equations. The said simultaneous equations are generated and updated using particle swarm optimization (PSO) technique, the fitness function being the sum of squared deviations from the actual set of observations. The data used for the calculations are intentionally corrupted to simulate real-life conditions. Results show that the proposed scheme offers a ...
number Fax (501) 228-3600 (501) 228-3601 Surface Water Data Information Requests Bill Baldwin (501) 228-3602 Surface Water General Information Jaysson Funkhouser (501) 228-3663...
Background Microbial water-quality indicators, in high concentrations in sewage, are used to determine whether water is safe for recreational purposes. Recently, the use of these indicators...Full Text Available
Before the implementation Water Framework directive, it was usual to forget that a good environment protection of the receiving waters needs a correct and coordinated operation of the subsystems of the water cycle, specially sewerage system, WWTP and receiving waters. This explains that most of the countries have focused their efforts in the treatment of dry weather flows forgetting the management of wet weather flows. Actually the idea that a sewerage system or a WWTP can not be planned or managed independently without considering the effects on the receiving waters is commonly accepted because not only each one of these systems must work correctly but also it is required a minimum impact in the receiving waters of the sewerage and WWTP overflows in dry and wet weather. All these links will affect the management strategy of the sewerage system (storm ...
The effect of the water-chemical regime (WCR) on damage sustained by heating surfaces of steam generators at NPP is analyzed. It is indicated that phosphate treatment with minimal excesses of phosphates in the steamgenerator water is the most optimal method of managing the WCR regime of horizontal steam generators.
The effect of water stress on patterns of nitrate reductase activity in the leaves and nodules and on nitrogen fixation were investigated in Medicago sativa L. plants watered 1 week...Full Text Available
The reference medium for brachytherapy dose measurements is water. Accuracy of dose measurements of brachytherapy sources is critically dependent on precise measurement of the source–detector...Full Text Available
A model was developed to compute concentration of impurities in the WWER 440 steam generator (SG) secondary water along the tube bundle. Calculated values were verified by concentration values obtained from secondary water sample chemical analysis. (orig.). 2 refs.
A method of measuring the water potential of stored potato tubers (Solanum tuberosum L.) was needed to investigate the relationship of bacterial soft rot in tubers to water potential....Full Text Available
By the time of primary 21st century, water hyacinth had become a serious environmental problem in China. Water hyacinth contributes to the major part of ecological hazards from the invasion of foreign...Full Text Available
Background and AimsXylem vessels containing gases (embolized) must be refilled with water if they are to resume transport of water through the plant, so refilling is of great importance...Full Text Available
average bulk water temper ature rise, OF bulk water temperature at elevation z, OF bulk water temperature in channels 0 and 1, O F film temperature, OF ...
Object: To smoothly control automatic water supply for realizing stable operation of a nuclear reactor by providing a flow rate limiting signal selection circuit and a preferential circuit in a water supply control device for a nuclear reactor wherein the speed of a recirculation pump may be changed in two-steps. Structure: Opening angle signals for a water supply regulating valve are controlled by a nuclear reactor water level signal, a vapor flow rate signal and a supplied water flow rate signal through an adder and an adjuster in response to a predetermined water level setting signal. When the water in the reactor is maintained at a predetermined level, a selection circuit receives a water pump condition signal for selecting one of the signals from a supplied water rate limiting signal generator ...
Inadequate measures of water quality have been used in many studies of the health effects associated with water supplies in developing countries. The present 1-year epidemiological-microbiological study...Full Text Available
THMs are disinfection by-products (DBPs) generated during water chlorination. Concentration of individual and total THMs, depends on treatment process and THMs precursors level. ATLL water utility has two DWTP (Llobregar and Ter) that produce and supply drinking water to Barcelona and regional area. This work studies the levels of THMs along the ATLL distribution system (450 km). Although, no differences were observed along water pipes system, changes of water resource and mix procedures were related. (Author) 12 refs.
In a supercritical water-cooled reactor, property of water changes significantly around the critical point. It is expected that irradiation and change of water property will affect the chemistry and material corrosion. Deep understanding of interactions between supercritical water and materials under irradiation is important. However, comprehensive data on radiolysis, kinetics, corrosion and thermodynamics have not been obtained due to the severe experimental condition. To get such data by experiments and computer simulations, a national program funded by Ministry of Education, Culture, Sports, Science and Technology (MEXT) has been started since December 2002. (author)
In this paper the diffusion of tritiated water vapor into concrete walls is studied to evaluate tritiated water retention capacity of a fusion reactor concrete building. Using a model of the tritiated water diffusion determined form experimental results, depth profiles of tritiated water in concrete are calculated in the case of being exposed to air containing tritiated water vapor during the normal operational condition of a fusion reactor. A 0.5-m-thick concrete is sufficient for reactor hall walls from a viewpoint of the tritium containment.
In this paper the diffusion of tritiated water vapor into concrete walls is studied to evaluate tritiated water retention capacity of a fusion reactor concrete building. Using a model of the tritiated water diffusion determined form experimental results, depth profiles of tritiated water in concrete are calculated in the case of being exposed to air containing tritiated water vapor during the normal operational condition of a fusion reactor. A 0.5-m-thick concrete is sufficient for reactor hall walls from a viewpoint of the tritium containment.
A non-linear analysis of a Cerenkov maser is presented. The system consists of a ring configuration of a cylindrical waveguide filled with a dielectric material. A single transverse-magnetic mode is assumed to propagate in the system. A low-density pencil electron beam travels in part of the ring, confined by a strong axial magnetic field. Using the single-particle description for the beam and the wave equation for the field, we obtain a set of two coupled non-linear differential equations describing the slowly varying amplitude and phase of the electromagnetic mode. The gain per path is assumed to be small and the spatial growth of the field is neglected. The resulting time dependent amplitude includes the exponential gain of the linear stage and the saturation to its maximum value. The time dependent frequency is also calculated. The two equations are combined to a single Van Der Pol equation with a ...
VSP data are usually recorded in common-shot gather and contain less informations. It is difficult, therefore, to obtain good migration result by using conventional migration methods. Both inner interface reflections and interformational multiples can be removed by reverse-time migration with the use of two-way nonreflection wave equation: besides, boundary reflection in VSP can be eliminated using Clayton's absorption boundary condition. It is proved that reverse-time migration is able to improve migration quality. The imaging condition is still obtained by generalizing Claerbout's imaging principle and solving eikonal equation. This is a very practical wave equation migration method. This method is applicable to prestack VSP data, conventional multiple-offset shot gather and the simultaneous migration of VSP data and surface seismic data. Besides, it is adaptable to various velocity variations. Before ...
Determination of the specific surface of uranium oxides: UO_2 and U_3O_8 using a simplified equation derived from the Brunauer, Emmett and Teller theory. The method is suitable for samples having a surface between 6 to 50 m"2.
This paper describes a long term course of research into the fundamentals of two phase flow regime dynamics that had lead to a differential description of regime properties. The development has included experiment, microscale simulation, and mathematical averaging to produce a dynamic equation for interfacial area.
We solve the loop equations of the $\\beta$-ensemble model analogously to the solution found for the Hermitian matrices $\\beta=1$. For \\beta=1$, the solution was expressed using the algebraic spectral curve of equation $y^2=U(x)$. For arbitrary $\\beta$, the spectral curve converts into a Schr\\"odinger equation $((\\hbar\\partial)^2-U(x))\\psi(x)=0$ with $\\hbar\\propto (\\sqrt\\beta-1/\\sqrt\\beta)/N$. This paper is similar to the sister paper~I, in particular, all the main ingredients specific for the algebraic solution of the problem remain the same, but here we present the second approach to finding a solution of loop equations using sectorwise definition of resolvents. Being technically more involved, it allows defining consistently the B-cycle structure of the obtained quantum algebraic curve (a D-module of the form $y^2-U(x)$, where $[y,x]=\\hbar$) and to construct explicitly the correlation ...
Three new classes of finite, homogeneous and nonsingular solutions of Einstein's equations which have time-dependent expansion, shear and rotation are presented. The t-constant sections are of Bianchi type IX. The source of these geometries is a fluid which has not been thermalized. (Author).
An exact cosmological solution of Einstein's equations which has time-dependent rotation is presented. The t-constant sections are of Bianchi type II. The source of this geometry is a fluid which has not been thermalized. (Author).
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.
We show that, in addition to radiation travelling at the speed of light, QED theory predicts a second type of radiation with an infinite velocity. We also show that charge, as it appears in the Dirac equation, may have a triune nature.
This paper considers the time dependent Stefan problem with convection in the fluid phase governed by the Stokes equation, and with adherence of the fluid on the lateral boundaries. The existence of a weak solution is obtained via the introduction of a te...
Test firings of 120-mm rounds were performed to characterize effects due to the combustible cartridge cases used. Comparisons were made between measured ballistic parameters, especially gas pressures, and those predicted by computerized simulation calcula...
A certain model of one-dimensional detonation waves leads to a Stefan problem: the unknown f satisfies Burgers equations on the two sides of a moving discontinuity at which it is given (f, say) and the jump in it derivative (corresponding to the exothermi...
A theoretical model for spray evaporation in a direct-injection diesel engine is developed taking into consideration the unsteady evaporation period. For this purpose, the liquid droplets temperatures, masses and velocities are studied. Solving the conservation of energy and mass equations together with the equation of state, a set of ordinary differential equations is achieved. These equations are solved using a special numerical technique. This technique, based on a fourth-order Runge-Kutta integration method is described. Calculated results are compared with experimental data and were found to be in good agreement. Analysis of the obtained results shows that for the combustion-chamber constituents, temperature is a dominant factor in the spray-evaporation process and in shaping the temperature history of the evaporating droplets. The combustion-chamber pressure and the spray volume play a secondary ...
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. ...
5 days ago... Sun sensors, 2 secondary Sun sensors, and an Earth sensor. ... This continued intermittently until the rocket burn, at which time ... with incoming asymptotic direction at an angle of -13.6 degrees from the lunar equator. ...
Background. Detection of subjects with early chronic kidney disease (CKD) is important because some will progress up to stage 5 CKD, and most are at high risk of cardiovascular morbidity...Full Text Available
Critical properties of a fluid or fluid mixtures are important for describing fluid phase behavior, predicting physical properties, developing equations of state, and designing supercritical-fluid extraction processes, and compression and refrigeration units. The predictive capability of the Peng-Robinson-Styjek-Vera (PRSV-2) equation of state (1986) for critical properties of binary mixtures was investigated. The procedure adopted by Heidemann and Khalil (1980) and discussed by Abu-Eishah et al. (1998) was followed. An optimized value for the binary interaction parameter based on minimization of error between experimental and predicted critical temperatures was used. The standard and the average of the absolute relative deviations in critical properties are included. The predicted critical temperature and pressure for several nonpolar and polar systems agree well with experimental data and are always better than those predicted by the ...
We study a rapidly convergent modification of a two-sided method for the approximate integration of a boundary-value problem with parameters in boundary conditions for a system of quasilinear second-order differential equations.
A (n +1)-dimensional Einstein-Gauss-Bonnet (EGB) model is considered. For diagonal cosmological type metrics the equations of motion are reduced to a set of Lagrange equations. The effective Lagrangian contains two minisuperspace metrics on R^n. The first one is the well-known 2-metric of pseudo-Euclidean signature and the second one is the Finslerian 4-metric that is proportional to n-dimensional Berwald-Moor 4-metric. When a synchronous-like time gauge is considered the equations of motion are reduced to an autonomous system of first order differential equations. For the case of the pure Gauss-Bonnet model two exact solutions with power-law and exponential dependence of scale factors (w.r.t. synchronous-like variable) are obtained. In cosmological case the first solution was obtained earlier by A. Toporensky and P. Tretyakov for n = 5 (and verified for n = 6,7). A generalization of the effective ...
The B test series from the course of ISP Nr. 43 was analysed. The boundary conditions measured include the initial temperature of the primary system, the front/slug injection flowrate and temperature, and the pressure drop across the core. Temperature data were collected at 185 thermocouple positions in the downcomer and 38 positions in the lower plenum. The frequency of data acquisition for code predictions was set to 2 Hz, which corresponds to the acquisition frequency of the two experimental setups. Calculations were performed using the FLUENT Computational Fluid Dynamics (CFD) code. This program is based on the finite volume method, The computational mesh was generated by the pre-processor - GAMBIT program. FLUENT uses a control-volume-based technique to convert the governing equations to algebraic equations, which can be solved numerically. This control volume technique consists of integrating the governing equations ...
Using the Klein-Gordon equation with a box potential, a mass formula describing the family of nonflavored meson states with I = 1 is derived. The energy levels calculated agree with those observed within an accuracy of approx. 5%. In the model discussed quarks behave like tachyons.
The methods of statistical dynamics are applied to a fluid with 5 conserved fields (the mass, the energy, and the three components of momentum) moving in a given external potential. When the potential is zero, we recover a previously derived system of parabolic differential equations, called "corrections to fluid dynamics".
We present several new results on the classic problem of estimating Gaussian profile parameters from a set of noisy data, showing that an exact solution of the maximum likelihood equations exists...Full Text Available
We show the existence of global-in-time weak solutions to a general class of coupled FENE-type bead-spring chain models that arise from the kinetic theory of dilute solutions of polymeric liquids with noninteracting polymer chains. The class of models involves the unsteady incompressible Navier-Stokes equations in a bounded domain in two or three space dimensions for the velocity and the pressure of the fluid, with an elastic extra-stress tensor appearing on the right-hand side in the momentum equation. The extra-stress tensor stems from the random movement of the polymer chains and is defined by the Kramers expression through the associated probability density function that satisfies a Fokker-Planck-type parabolic equation, a crucial feature of which is the presence of a center-of-mass diffusion term. We require no structural assumptions on the drag term in the Fokker-Planck equation; in particular, ...
Main features of star evolution are recalled. Then the general structure of white dwarf stars is examined. From the equation of state of an electron gas completely degenerated are deduced: mechanical equilibrium, Viriel theorem, mass-radius relationship and Chandrasekhar limit. These results are applied to neutron stars.
We introduce in this paper a new method for reducing neurodynamical data to an effective diffusion equation, either experimentally or using simulations of biophysically detailed models. The dimensionality...Full Text Available
We investigate the default time of a firm when a stochastic discount factor is used so that both diffusion and regime switching risks are priced. We establish the relationship between the probability distribution of the default time and the solution of a system of coupled partial differential equations.
It is urged that the lesson of gauge invariance in quantum electrodynamics implies the irrelevance of `Schwinger term` difficulties in current algebra. The divergence equations of Veltman form the basis of a gauge-variation formalism in which these questions are avoided. (author). 9 refs.
The scattering of a charged scalar field on Coulomb potential on de Sitter space-time is studied using the solution of the free Klein-Gordon equation. We find that the scattering amplitude is independent of the choice of the picture and in addition the total energy is conserved in the scattering process.
A method for calculating plume rise above stack level is presented. The equations set forth by Briggs, which are presently the most popular for such calculations, are discussed. A method using 2 nomographs, simplifying the calculations is given. (JMT)
This paper describes the development of an interactive internet-based cost-estimating tool for commonly used urban storm runoff best management practices (BMP), including: retention and detention ponds, grassed swales, and constructed wetlands. The paper presents the cost data, c...
We consider a class of nonlinear Schrodinger / Gross-Pitaevskii (NLS/GP) equations with periodic potentials, having an even symmetry. We construct "solitons", centered about any point of symmetry of the potential. For focusing (attractive) nonlinearities, these solutions bifurcate from the zero state at the lowest band edge frequency, into the semi-infinite spectral gap. Our results extend to bifurcations into finite spectral gaps, for focusing or defocusing (repulsive) nonlinearities under more restrictive hypotheses. Soliton nonlinear bound states with frequencies near a band edge are well-approximated by a slowly decaying solution of a homogenized NLS/GP equation, with constant homogenized effective mass tensor and effective nonlinear coupling coefficient, modulated by a Bloch state. For the critical NLS equation with a periodic potential, e.g. the cubic two dimensional NLS/GP with a periodic potential, our results imply ...
In this paper, an amount of confining reinforcement required for the sufficient flexural ductility of reinforced concrete column was investigated. An equation to calculate the required confining reinforcement ratio was deduced from the analytical investigation and examination on the previous experimental data as to the effects of various factors on the required confining reinforcement ratio. The rotation angle required in plastic hinges, the axial force ratio, the material properties such as concrete strength, and the ratio of core area to the whole area of column sections were employed as the factors. From the examination of the equation, it was shown that the ultimate rotation in plastic hinges calculated by the equation was more conservative than the ultimate rotaion angle of members obtained by test column specimens. It was also indicated that the ultimate rotation angle of members of columns failed in shear was smaller ...
An ascending node occurs when the subsatellite point crosses the earth's equator going from .... 11dB sigma-naught. This error budget includes the altimeter noise, ...... it is also recommended to filter the data as follows to ...
Wood, R. A.; and Ogden, H. R.: The All-Beta Titanium Alloy (Ti-13V-llCr-3Al). DMIC Rept. 110 (ASTIA AD 214002), Battelle Mem. Inst., Apr. 17, 1959. 22. White ...
...Colormetric System, which lie within the region bounded by the spectrum locus and lines defined by the following equations: X...intensity of a light in candela as defined by the Illuminating Engineering Society's Guide for Calculating the Effective...
Our work has focused on the development and analysis of domain decomposition algorithms for a variety of problems arising in continuum mechanics modeling. In particular, we have extended and analyzed FETI-DP and BDDC algorithms; these iterative solvers were first introduced and studied by Charbel Farhat and his collaborators, see [11, 45, 12], and by Clark Dohrmann of SANDIA, Albuquerque, see [43, 2, 1], respectively. These two closely related families of methods are of particular interest since they are used more extensively than other iterative substructuring methods to solve very large and difficult problems. Thus, the FETI algorithms are part of the SALINAS system developed by the SANDIA National Laboratories for very large scale computations, and as already noted, BDDC was first developed by a SANDIA scientist, Dr. Clark Dohrmann. The FETI algorithms are also making inroads in commercial engineering software systems. We also note that the analysis of these ...
The transport of liquid water and gaseous reactants through a gas diffusion layer (GDL) is one of the most important water management issues in a proton exchange membrane fuel cell (PEMFC). In this work, the liquid water breakthrough dynamics, characterized by the capillary pressure and water saturation, across GDLs with and without a microporous layer (MPL) are studied in an ex-situ setup which closely simulates a real fuel cell configuration and operating conditions. The results reveal that recurrent breakthroughs are observed for all of the GDL samples tested, indicating the presence of an intermittent water drainage mechanism in the GDL. This is accounted for by the breakdown and redevelopment of the continuous water paths during water drainage as demonstrated by Haines jumps. For GDL ...
Water is one of our most precious resources. However, for many in the United States, having fresh, safe drinking water is taken for granted, and due to this perceived lack of relevance, students may not fully appreciate the luxury of having safe running water--in the home. One approach to resolving water-quality issues in the United States may reside in providing education that presents accurate information in a meaningful way. Accordingly, this article describes a unit designed to emphasize the importance of water-quality testing and purification and to introduce students to local water-quality issues. The engineering-based module of this eighth-grade science activity is particularly important due to the design-build-test component. (Contains 5 figures.)
Most plans for the disposition of surplus nuclear materials involve storage in sealed containers where the evolution of gases from reactions of adsorbed water could present both pressure and flammability hazards[l] . Despite efforts such as calcining the material to minimize the water content prior to packaging, both residual moisture and readsorbed water may be present in the final containers . Given the anticipated temperature excursions during transportation and storage, this water may thermally desorb, increasing the pressure, and/or thermally dissociate to produce H2 gas, increasing flammability hazards . In addition, the radiation from the nuclear material may induce radiolysis of the water with the likely products being water vapor, H2, 02 and H2O2. In order to better understand the relative importance of the thermal- and radiationinduced chemistry, we ...
Ground water is treated in the Czech Republic so that small water tanks are built above each water source to serve the primary ground water treatment; water so pretreated is then concentrated in large basins for subsequent treatment. Some water tanks where the first contact of the ground water with air takes place were selected as sites predisposed to radon accumulation. The examination was carried out near the town of Jihlava, where the bedrock contains slightly elevated radium concentrations. The average radon concentrations lay within the region of 2-4 kBq/m"3; the instantaneous values, however, exhibited appreciable periodical variations during the day. The relatively high radon concentrations will not pose a marked hazard for the personnel because the employees only reside at the sites in question for 10 to 15 minutes within 2 days, not ...
A single cylinder diesel engine study of water-in-diesel emulsions was conducted to investigate the effect of water emulsification on the engine performance and gases exhaust temperature. Emulsified Diesel fuels of 0, 5, 10, 15 and 20 water/diesel ratios by volume were used in a single cylinder, direct injection diesel engine, operating at 1200-3300 rpm. The results indicate that the addition of water in the form of emulsion improves combustion efficiency. The engine torque, power and brake thermal efficiency increase as the water percentage in the emulsion increases. The average increase in the brake thermal efficiency for 20% water emulsion is approximately 3.5% over the use of diesel for the engine speed range studied. The proper brake specific fuel consumption and gases exhaust temperature decrease as the percentage of water in the ...
... Marine transport and oil disposal by boats and Matadi pose threats in the region. When the water ... Marine transport and oil disposal by boats and Matadi pose threats in the region. When the water ... ...
This site explains how temperature, pressure, and salinity work together to determine the density of ocean water. The three density layers of the ocean are described by means of text description and a graphic illustration.
Using air as a carrier gas is common in humidification?dehumidification desalination processes. A computer program was written using mass and energy balances for modeling the process behavior. The parameters considered in this work were inlet air and fresh water recycle temperatures, inlet air flow rate, saline water and fresh water recycle flow rates, and saline water to air flow ratio. Results of simulation showed that increasing inlet air and fresh water recycle flow rate increases fresh water production. It was also found that heating the inlet air to humidification column or cooling the inlet water to dehumidification column increases the production rate but increasing water to air flow ratio in a humidifier leads to a lower production rate. The predicted effects of the parameters on ...
The largest outbreak of cryptosporidiosis reported in the United Kingdom, involving 575 confirmed cases (of which 474 met an agreed case definition), occurred in the county of Devon during August and September of 1995. The descriptive epidemiology supports the hypothesis that the outbreak was associated with the consumption of cold tap water in the area served by a particular water treatment works. Cryptosporidium oocysts were detected in treated water samples at the time of the outbreak. Although the epidemiological analysis provided strong circumstantial evidence of a waterborne outbreak, the data were not recorded in a manner that made them admissible in criminal proceedings taken by the Drinking Water Inspectorate against the water company involved. The need to carry out an analytical study in conjunction with the identification and characterisation of the pathogen in the ...
metric temperatures has been prepared showing the surface water temper- ... occasion, gray levels show different densities for similar water temper- ...
... should give full con- As for the effects of the Three Gorges on the mouth of sideration to the rising water process and water temper- the river, the ...
ObjectiveTo evaluate the effect of moderate aerobic physical activity in water on fetal cardiotocography patterns in sedentary pregnant women.MethodIn...Full Text Available
The very low-pressure expansion tank of the title invention is connected to the water in the central heating installation via a connecting pipe with a pump and valves on one side, and on the other side the tank is connected via a connecting pipe with valve to the tap water mains, so that the supply of water can be regulated automatically. Within the expansion tank contact with the outside air is not possible because of an air/water separating floater. By means of recording and control (also remote) of the contents of the expansion tank, the installation pressure and the quantity of supplied water from the expansion tank and the tap water mains, failures and water damage are prevented. 4 figs.
An environmental study was done to examine the prevalence of Vibrio mimicus in some aquatic environments of Dhaka, Bangladesh, and of Okayama, Japan. Water samples from Dhaka environments and water...Full Text Available
density is slightly higher, its electrical conductivity is much higher, and it is slightly ... Sea water has characteristic properties (e.g. density) that are independent of ...
In the Mediterranean arc, vines for wine production are mainly grown without the support of irrigation. Under such conditions, site variables affecting the extent and seasonal timing of water deficits are the dominant environmental constraints for grape production. Moreover, water availability and vine water status are the factors most comprehensively determining fruit composition and, thus, wine quality. Therefore, monitoring the extent of water stress in vines might be a valuable tool for the optimisation of grape yield and quality. The objective of this study was to evaluate the feasibility of using the reflectance based Water Index (WI) to estimate vine water status at the leaf and canopy levels. The study was conducted on Vitis vinifera cv. Chardonnay potted plants submitted to contra...
Halothane, chloroform, and carbon tetrachloride, in the vapor and liquid phases, stimulate the water receptor of the blowfly Phormia regina. There are three successive phases of response to long-lasting...Full Text Available
An exact solution method for the free vibration problem of thin circular cylindrical shells is presented. The differential equations of motion are solved directly with the use of simple Fourier series as the modal displacement functions. Stokes' transformation is exploited to obtain correct series expressions for the derivatives of the Fourier series. From this method an explicit expression of the exact frequency equation can be obtained for any kind of boundary conditions. The accuracy of the present method is checked against available data. The proposed method is then used to find the modal characteristics of the thermal liner of the Fast Test Reactor (FTR). The numerical results obtained are compared with finite element method solutions. (auth)
A new class of solutions of the euclidean Einstein equations with #LAMBDA#-term ( A-class ) is found by solving the complex two-dimensional Kaehler Einstein equations with the following realization of complex metrics. The A-Class includes two gravitational instantons already known: the CP"2 and the Eguchi-Hanson metric, and allows a U(1)-generalized spin structure. It is shown that all Einstein euclidean two-axial Bianchi type IX metrics are exhausted by the Taub-NUT-de Sitter family and the A-class. (orig.).
In terms of Berezins's theory of symbols of operators, the integral formulation is suggested for the free differential algebra which gives rise to consistent equations of motion of interacting massless fields of all spins 0#<=#s<#infinity# in the frameworks of gravity. In the first nontrivial order of the expansion in powers of curvatures, Frobenius consistency conditions for higher-spin equations of motion are shown to reduce to the simple geometrical fast that there are two ways for splitting any quadrangle in two triangles. To clarify our construction, we illustrate how it works in the simplest case of pure gravity. (orig.).
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.
The Rabi Hamiltonian, describing the coupling of a two-level system to a single quantized boson mode, is studied in the Bargmann-Fock representation. The corresponding system of differential equations is transformed into a canonical form in which all regular singularities between zero and infinity have been removed. The canonical or Birkhoff-transformed equations give rise to a two-dimensional eigenvalue problem, involving the energy and a transformational parameter which affects the coupling strength. The known isolated exact solutions of the Rabi Hamiltonian are found to correspond to the uncoupled form of the canonical system.
Muonic-atom energy levels and wave functions are calculated. The results are corrected for nearly all important static effects, including finite nuclear size as produced by a user-specified nuclear charge distribution. Fourth-order Runge-Kutta integration of the radial Dirac equations, with the secant method determining the eigenvalues. Corrections are either included as potentials in the Dirac equations or computed as first-order perturbations. The code assumes spherical symmetry, so that effects of nuclear deformation may be treated only as static first-order perturbations. Dynamic muon-electron and muon-nuclear effects are neglected (the latter are calculated by RURP). (Auth.).
We present a dynamical spectral model for Large Eddy Simulation of the incompressible magnetohydrodynamic (MHD) equations based on the Eddy Damped Quasi Normal Markovian approximation. This model extends classical spectral Large Eddy Simulations for the Navier-Stokes equations to incorporate general (non Kolmogorovian) spectra as well as eddy noise. We derive the model for MHD and show that introducing a new eddy-damping time for the dynamics of spectral tensors in the absence of equipartition between the velocity and magnetic fields leads to better agreement with direct numerical simulations, an important point for dynamo computations.
In symmetrical nuclear matter the solutions of pion dispersion equation are investigated in the complex plane of the pion frequency $\\omega$. There are three well-known branches of solutions on the physical sheet : sound, pion and isobar, at the matter density less than the critical one $\\rho <\\rho_c$. At the condition $\\omega^2_c\\leq0$ (in general case Re$ \\omega^2_c\\leq0$) takes place. This points out the instability of the ground state which is possibly related to the pion condensation.
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.
This paper presents a reverse time migration (RTM) method for the migration of shot records in tilted transversely isotropic (TTI) media. It is based on the tilted TI acoustic wave equation that was derived from the dispersion relation. The RTM is a full depth migration allowing for velocity to vary laterally as well as vertically and has no dip limitations. The wave equation is solved by a tenth-order finite difference scheme. Using 2D numerical models, we demonstrate that ignoring the tilt angle will introduce both lateral and vertical shifts in imaging. The shifts can be larger than 0.5 wavelength in the vertical direction and 1.5 wavelength in the lateral direction.
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 paper develops a new class of financial market models. These models are based on generalized telegraph processes: Markov random flows with alternating velocities and jumps occurring when the velocities are switching. While such markets may admit an arbitrage opportunity, the model under consideration is arbitrage-free and complete if directions of jumps in stock prices are in a certain correspondence with their velocity and interest rate behaviour. An analog of the Black-Scholes fundamental differential equation is derived, but, in contrast with the Black-Scholes model, this equation is hyperbolic. Explicit formulas for prices of European options are obtained using perfect and quantile hedging.
A simplified approach has been proposed to model the drying of capillary-porous materials with small geometric parameters for practical purposes. Luikov`s theory of mass and heat transfer for capillary-porous bodies was applied to drying of biological materials and it was shown that the thermo diffusion effect can be neglected if the thermal gradient remains small during the whole period of drying. The resulting simplified equations can then be solved using available mathematical software. The simplified equations were applied to drying of wood and a peanut pod
Isosteric heats and adsorption isotherms of methane on nanometer active carbon were measured at supercritical temperature (273-373 K) and pressure from 0 to 10 MPa. The measured data agreed well with Dubinin-Astakhov (DA) model at lower pressure but failed when pressure exceeded a special range. General Freundlich (GF) equation was used to modify the DA equation at high pressure and thus formed a combined bisection model GFDA. The adsorption mechanism of methane on nanometer active carbon was raised according to GFDA model and the heterogeneous energy distribution of the adsorbent was analyzed.
We investigate stability properties of indirectly damped systems of evolution equations in Hilbert spaces, under new compatibility assumptions. We prove polynomial decay for the energy of solutions and optimize our results by interpolation techniques, obtaining a full range of power-like decay rates. In particular, we give explicit estimates with respect to the initial data. We discuss several applications to hyperbolic systems with {\\em hybrid} boundary conditions, including the coupling of two wave equations subject to Dirichlet and Robin type boundary conditions, respectively.
The entrainment flow rate distribution, the gas velocity profile, and the concentration profile of droplets across the channel cross section in fully developed region of a horizontal rectangular channel of 150mm width and 50mm height were measured. The concentration profile of droplets was expressed by a simple equation based on a constant diffusion coefficient model. From this equation the effects of gravity and turbulent diffusion of droplets on the concentration profile were evaluated. The characteristic mean settling velocity of a group of droplets with various diameters was derived, and using this value the mean diffusion coefficient of the group of droplets was obtained
The entrainment flow rate distribution, the gas velocity profile, and the concentration profile of droplets across the channel cross section in fully developed region of a horizontal rectangular channel of 150 mm width and 50 mm height were measured. The concentration profile of droplets was expressed by a simple equation based on a constant diffusion coefficient model. From this equation the effects of gravity and turbulent diffusion of droplets on the concentration profile were evaluated. The characteristic mean settling velocity of a group of droplets with various diameters was derived, and using this value the mean diffusion coefficient of the group of droplets was obtained. (author).
Bases of the algorithm for assessing the reinforcement status in reinforced concrete products by gamma-absorption method are presented. Analytical equations are obtained for estimation of error of the parameter characterizing the degree of reinforcement destruction. It is recommended to use high-energy Bremsstrahlung sources-betatrons of 4-10 MeV maximum energy for testing products of 500-600 mm thick. Linear radiation attenuation factor (LAF) of concrete in estimated equation is replaced by effective LAF of concrete, and LAF of reinforcement and corrosion materials - by differential LAF. Corresponding LAF of nonevident form in the algorithm is assessed by the results of processing of direct-shadow radiographs of defectometers
A simple closed-form equation to calculate the envelope correlation between any two receiver or transmitter antennas in a multi-input multi-output (MIMO) system of an arbitrary number of elements is derived. The equation uses the scattering parameters obtained at the antenna feed point to calculate the envelope correlation coefficient. This approach has the advantage that it does not require knowledge of the antenna radiation pattern. Numerical data that include conductor and permittivity loss are shown to validate the approach. Using the scattering parameters for calculating the envelope correlation is less laborious and relatively cheaper, as compared to the radiation pattern approach.
An energy control strategy (10) for a hybrid electric vehicle that controls an electric motor during bleed and charge modes of operation. The control strategy (10) establishes (12) a value of the power level at which the battery is to be charged. The power level is used to calculate (14) the torque to be commanded to the electric motor. The strategy (10) of the present invention identifies a transition region (22) for the electric motor's operation that is bounded by upper and lower speed limits. According to the present invention, the desired torque is calculated by applying equations to the regions before, during and after the transition region (22), the equations being a function of the power level and the predetermined limits and boundaries.
The present contribution deals with the thermophoresis particle deposition and thermal radiation effects on the flow, heat and mass transfer characteristics in a viscous fluid over a semi-infinite vertical porous plate. The governing boundary layer equations are written into a dimensionless form by similarity transformations. The transformed coupled nonlinear ordinary differential equations are solved numerically by means of the fourth-order Runge?Kutta method with a shooting technique. The effects of different parameters on the dimensionless velocity, temperature, and concentration profiles are shown graphically. In addition, results for the local skin-friction coefficient, the local Nusselt number, and the local Sherwood number are tabulated and discussed.
The present contribution deals with the thermophoresis particle deposition and thermal radiation effects on the flow, heat and mass transfer characteristics in a viscous fluid over a semi-infinite vertical porous plate. The governing boundary layer equations are written into a dimensionless form by similarity transformations. The transformed coupled nonlinear ordinary differential equations are solved numerically by means of the fourth-order Runge-Kutta method with a shooting technique. The effects of different parameters on the dimensionless velocity, temperature, and concentration profiles are shown graphically. In addition, results for the local skin-friction coefficient, the local Nusselt number, and the local Sherwood number are tabulated and discussed.
We find a large family of solutions to the Dirac equation on a manifold of $G_2$ holonomy asymptotic to a cone over $S^3 \\times S^3$, including all radial solutions. The behaviour of these solutions is studied as the manifold developes a conical singularity. None of the solutions found are both localised and square integrable at the origin. This result suggests that M-theory on the conifold over $S^3\\times S^3$ does not produce a chiral theory. The approach here is complementary to previous analyses using dualities and anomaly cancellation which found chiral fermions on all known $G_2$ conifolds except that over $S^3\\times S^3$.
A generalization of Faddeev's group cohomology applicable to diffeomorphism groups is presented. This cohomology is used to calculate the two cocycles associated with a projective representation of the diffeomorphism group on the circle. The group version of the n"3 term descends from a three dimensional Chern-Simons action based on the diffeomorphism group. The group version of the n term arises from an ambiguity in the descent equations of adding closed but not exact forms and is trivial only if its appropriately normalized coefficient is quantized to be an integer. Finally, a hamiltonian interpretation of global anomalies is suggested in the language of group cohomology. (orig.).
The chemical reaction effect on an unsteady magnetohydrodynamic (MHD) flow past a semi-infinite vertical porous plate with viscous dissipation is analyzed. The governing equations of motion, energy, and species are transformed into ordinary differential equations (ODEs) using the time dependent similarity parameter. The resultant ODEs are then solved numerically by a finite element method. The effects of various parameters on the velocity, temperature, and concentration profiles are presented graphically, and the values of the skin-friction, Nusselt number, and Sherwood number for various values of physical parameters are presented through tables.
Burnout data with low pressure Freon-113 for even- and odd- numbered peripheral rod clusters with relatively large spacings were used to derive equations in terms of dimensionless parameters suggested by Barnett. The equations which are for three different flow regimes for each rod geometry (even or odd) were found to predict burnout data with maximum RMS deviation being 3.8%. (author). 11 figs., 3 tabs., 15 refs.
This article proposes a model for estimating the effect of the overall thermal resistance of the absorber and the coolant on a solar thermal energy system utilizing advanced heat transfer mediums. Two equations can be used to model a dynamic system: one for the maximum conversion efficiency and one to specify the constraint that, in quasi-steady state, the net energy collected must be transferred to the coolant. These two equations define a constrained optimization problem in two variables, the fluid and the absorber temperatures. The results present interesting implications for liquid metal, direct absorption, and solar evaporating systems.
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.
The entire Virasoro, Ramond and Neveu-Schwarz algebras can each be constructed from a finite number of well-chosen generators satisfying a small number of conditions. Our most economical sets consist of just two starting generators in all cases, subject to no more than six conditions for the Virasoro case, five conditions for the Ramond case, and nine conditions for the Neveu-Schwarz case. Consequently, the Virasoro algebra simply amounts to 6 equations in two operator unknowns, and correspondingly 5 to 9 equations for the foregoing superalgebras. 2 refs.
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.
MAGMA is a FORTRAN computer code designed to viscous flow in in situ vitrification melt pools. It models three-dimensional, incompressible, viscous flow and heat transfer. The momentum equation is coupled to the temperature field through the buoyancy force terms arising from the Boussinesq approximation. All fluid properties, except density, are assumed variable. Density is assumed constant except in the buoyancy force terms in the momentum equation. A simple melting model based on the enthalpy method allows the study of the melt front progression and latent heat effects. An indirect addressing scheme used in the numerical solution of the momentum equation voids unnecessary calculations in cells devoid of liquid. Two-dimensional calculations can be performed using either rectangular or cylindrical coordinates, while three-dimensional calculations use rectangular coordinates. All derivatives are approximated by finite ...
This paper, firstly, introduces the background in developing a high quality water supply system, in which the drinking water system is isolated inside a building. Results of questionnaire on the high quality service water are illustrated. The results of questionnaire have revealed that the high quality service water is extremely interested. Then, are described the target quality of high quality water, the constitution of high quality water supply system, the treatment process, the measures to secure safety and sanitary, and the method of maintenance. The high quality water is produced through the activated charcoal absorption treatment, membrane treatment, ozonation, cooling, mineral addition, and disinfection of city water. Furthermore, application examples of this system adopted to a department store and a bachelor ...
In an attempt to obtain basic data for evaluation of exposure doses in Niigata Prefecture, the concentrations of tritium in atmospheric water, precipitation, river water, and tap water were measured. Samples of atmospheric water, river water, and tap water were collected once for 2 weeks; and precipitation collected for one week was used as sample. The concentration of atmospheric tritium depended on the concentration of tritium in moisture and the content of water in atmosphere. Tritium levels were high in May, November and March, and low in September. Regarding tritium concentrations, there was a good correlation between atmospheric water and precipitation. Tritium concentrations in both of them varied from sample to sample. The concentrations of tritium in river and tap water tended to be high ...
An efficient method of hydrogen storage in nano-porous carbons is its reversible sorption by electrochemical decomposition of a KOH water solution [1-3] according to the following equation: C + xH{sub 2}O + xe{sup -} {yields} (CH{sub x}) + xOH{sup -} where (CH{sub x}) stands for the hydrogen inserted into the nano-porous carbon during charging and oxidized during discharging. Although various carbon materials have been investigated as hydrogen adsorbents, the information about the storage mechanism as well as the nature of the hydrogen/carbon interaction is still not sufficient. In order to extend the understanding of the process, carbon samples charged electrochemically were investigated by temperature programmed desorpt(TPD). The nature of the hydrogen/carbon interaction was studied by electrochemical analysis at different temperatures. The TPD experiments consist of heating the samples from room temperature to 950 C and of quantitative ...
The transport of radon in concrete takes place through the complicated network of interconnected pores that is, at any time, the result of the process of hydration of cement and of moisture distribution and transport. Initially the microstructure of concrete depends on the mix proportions and curing conditions, its time-evolution being conditioned by its surrounding environment. Radon transport will be consequently a function of time, as it is influenced by the changing microstructure (total porosity and its distribution) and by the amount and distribution of the moisture contained in the pore system. A selection of information from the large amount of literature available on concrete is presented in chapter 2. A model that describes the process of hydration, of microstructure development and of moisture transport is presented in chapter 3. The physics of radon diffusion in homogeneous porous materials is outlined in chapter 4. The coupling of the numerical implementation of the ...
The transport of radon in concrete takes place through the complicated network of interconnected pores that is, at any time, the result of the process of hydration of cement and of moisture distribution and transport. Initially the microstructure of concrete depends on the mix proportions and curing conditions, its time-evolution being conditioned by its surrounding environment. Radon transport will be consequently a function of time, as it is influenced by the changing microstructure (total porosity and its distribution) and by the amount and distribution of the moisture contained in the pore system. A selection of information from the large amount of literature available on concrete is presented in chapter 2. A model that describes the process of hydration, of microstructure development and of moisture transport is presented in chapter 3. The physics of radon diffusion in homogeneous porous materials is outlined in chapter 4. The coupling of the numerical implementation of the ...
This document summarizes the degree to which State water quality standards include wetlands in their definitions of State waters. The detailed review of other elements of standards, such as uses, criteria and antidegradation, is beyond the scope of this review. The purpose of this review is to provide a baseline of information for the development of EPA guidance and policies related to water quality standards for wetlands.
Water transport in the ionomeric membrane, typically Nafion{reg_sign}, has profound influence on the performance of the polymer electrolyte fuel cell, in terms of internal resistance and overall water balance. In this work, high resolution neutron imaging of the Nafion{reg_sign} membrane is presented in order to measure water content and through-plane gradients in situ under disparate temperature and humidification conditions.