Development and analysis of finite volume methods

International Nuclear Information System (INIS)

Omnes, P.

2010-05-01

This document is a synthesis of a set of works concerning the development and the analysis of finite volume methods used for the numerical approximation of partial differential equations (PDEs) stemming from physics. In the first part, the document deals with co-localized Godunov type schemes for the Maxwell and wave equations, with a study on the loss of precision of this scheme at low Mach number. In the second part, discrete differential operators are built on fairly general, in particular very distorted or nonconforming, bidimensional meshes. These operators are used to approach the solutions of PDEs modelling diffusion, electro and magneto-statics and electromagnetism by the discrete duality finite volume method (DDFV) on staggered meshes. The third part presents the numerical analysis and some a priori as well as a posteriori error estimations for the discretization of the Laplace equation by the DDFV scheme. The last part is devoted to the order of convergence in the L2 norm of the finite volume approximation of the solution of the Laplace equation in one dimension and on meshes with orthogonality properties in two dimensions. Necessary and sufficient conditions, relatively to the mesh geometry and to the regularity of the data, are provided that ensure the second-order convergence of the method. (author)

Numerical methods of mathematical optimization with Algol and Fortran programs

CERN Document Server

Künzi, Hans P; Zehnder, C A; Rheinboldt, Werner

1971-01-01

Numerical Methods of Mathematical Optimization: With ALGOL and FORTRAN Programs reviews the theory and the practical application of the numerical methods of mathematical optimization. An ALGOL and a FORTRAN program was developed for each one of the algorithms described in the theoretical section. This should result in easy access to the application of the different optimization methods.Comprised of four chapters, this volume begins with a discussion on the theory of linear and nonlinear optimization, with the main stress on an easily understood, mathematically precise presentation. In addition

Account of volume heat capacity on interface in numerical solution of the Stephen problem using the strained coordinates method

International Nuclear Information System (INIS)

Latynin, V.A.; Reshetov, V.A.; Karaseva, L.N.

1988-01-01

Numerical solution of the Stephen problem by the strained coordinate method is presented for an one-dimensional sphere. Differential formulae of heat fluxes from moving interfaces do not take into account volume heat capacities of the front nodes. Calculations, carried out according to these balanced formulae, as well as according to those usually used, have shown that the balanced formulae permit to reduce approximately by an order the number of nodes on the sphere radius, if similar accuracy of heat balance of the whole process of melting or crystallization is observed. 2 refs.; 1 fig

Numerical simulation of homogenization time measurement by probes with different volume size

International Nuclear Information System (INIS)

Thyn, J.; Novy, M.; Zitny, R.; Mostek, M.; Jahoda, M.

2004-01-01

Results of continuous homogenization time measurement of liquid in a stirred tank depend on the scale of scrutiny. Experimental techniques use a probe, which is situated inside as a conductivity method, or outside of the tank as in the case of gamma-radiotracer methods. Expected value of homogenization time evaluated for a given degree of homogenization is higher when using the conductivity method because the conductivity probe measures relatively small volume in contrast to application of radiotracer, when the volume is much greater. Measurement through the wall of tank is a great advantage of radiotracer application but a comparison of the results with another method supposes a determination of measured volume, which is not easy. Simulation of measurement by CFD code can help to solve the problem. Methodology for CFD simulation of radiotracer experiments was suggested. Commercial software was used for simulation of liquid homogenization in mixed vessel with Rushton turbine. Numerical simulation of liquid homogenization time by CFD for different values of detected volume was confronted with measurement of homogenization time with conductivity probe and with different radioisotopes 198 Au, 82 Br and 24 Na. Detected size of the tank volume was affected by different energy of radioisotope used. (author)

The volume of fluid method in spherical coordinates

NARCIS (Netherlands)

Janse, A.M.C.; Janse, A.M.C.; Dijk, P.E.; Kuipers, J.A.M.

2000-01-01

The volume of fluid (VOF) method is a numerical technique to track the developing free surfaces of liquids in motion. This method can, for example, be applied to compute the liquid flow patterns in a rotating cone reactor. For this application a spherical coordinate system is most suited. The novel

Development op finite volume methods for fluid dynamics; Developpement de methodes de volumes finis pour la mecanique des fluides

Energy Technology Data Exchange (ETDEWEB)

Delcourte, S

2007-09-15

We aim to develop a finite volume method which applies to a greater class of meshes than other finite volume methods, restricted by orthogonality constraints. We build discrete differential operators over the three staggered tessellations needed for the construction of the method. These operators verify some analogous properties to those of the continuous operators. At first, the method is applied to the Div-Curl problem, which can be viewed as a building block of the Stokes problem. Then, the Stokes problem is dealt with with various boundary conditions. It is well known that when the computational domain is polygonal and non-convex, the order of convergence of numerical methods is deteriorated. Consequently, we have studied how an appropriate local refinement is able to restore the optimal order of convergence for the Laplacian problem. At last, we have discretized the non-linear Navier-Stokes problem, using the rotational formulation of the convection term, associated to the Bernoulli pressure. With an iterative algorithm, we are led to solve a saddle-point problem at each iteration. We give a particular interest to this linear problem by testing some pre-conditioners issued from finite elements, which we adapt to our method. Each problem is illustrated by numerical results on arbitrary meshes, such as strongly non-conforming meshes. (author)

The spectral volume method as applied to transport problems

International Nuclear Information System (INIS)

McClarren, Ryan G.

2011-01-01

We present a new spatial discretization for transport problems: the spectral volume method. This method, rst developed by Wang for computational fluid dynamics, divides each computational cell into several sub-cells and enforces particle balance on each of these sub-cells. Also, these sub-cells are used to build a polynomial reconstruction in the cell. The idea of dividing cells into many cells is a generalization of the simple corner balance and other similar schemes. The spectral volume method preserves particle conservation and preserves the asymptotic diffusion limit. We present results from the method on two transport problems in slab geometry using discrete ordinates and second through sixth order spectral volume schemes. The numerical results demonstrate the accuracy and preservation of the diffusion limit of the spectral volume method. Future work will explore possible bene ts of the scheme for high-performance computing and for resolving diffusive boundary layers. (author)

Numerical evaluation of acoustic characteristics and their damping of a thrust chamber using a constant-volume bomb model

OpenAIRE

Jianxiu QIN; Huiqiang ZHANG; Bing WANG

2018-01-01

In order to numerically evaluate the acoustic characteristics of liquid rocket engine thrust chambers by means of a computational fluid dynamics method, a mathematical model of an artificial constant-volume bomb is proposed in this paper. A localized pressure pulse with a very high amplitude can be imposed on specified regions in a combustion chamber, the numerical procedure of which is described. Pressure oscillations actuated by the released constant-volume bomb can then be analyzed via Fas...

Advanced numerical methods for three dimensional two-phase flow calculations

Energy Technology Data Exchange (ETDEWEB)

Toumi, I. [Laboratoire d`Etudes Thermiques des Reacteurs, Gif sur Yvette (France); Caruge, D. [Institut de Protection et de Surete Nucleaire, Fontenay aux Roses (France)

1997-07-01

This paper is devoted to new numerical methods developed for both one and three dimensional two-phase flow calculations. These methods are finite volume numerical methods and are based on the use of Approximate Riemann Solvers concepts to define convective fluxes versus mean cell quantities. The first part of the paper presents the numerical method for a one dimensional hyperbolic two-fluid model including differential terms as added mass and interface pressure. This numerical solution scheme makes use of the Riemann problem solution to define backward and forward differencing to approximate spatial derivatives. The construction of this approximate Riemann solver uses an extension of Roe`s method that has been successfully used to solve gas dynamic equations. As far as the two-fluid model is hyperbolic, this numerical method seems very efficient for the numerical solution of two-phase flow problems. The scheme was applied both to shock tube problems and to standard tests for two-fluid computer codes. The second part describes the numerical method in the three dimensional case. The authors discuss also some improvements performed to obtain a fully implicit solution method that provides fast running steady state calculations. Such a scheme is not implemented in a thermal-hydraulic computer code devoted to 3-D steady-state and transient computations. Some results obtained for Pressurised Water Reactors concerning upper plenum calculations and a steady state flow in the core with rod bow effect evaluation are presented. In practice these new numerical methods have proved to be stable on non staggered grids and capable of generating accurate non oscillating solutions for two-phase flow calculations.

Advanced numerical methods for three dimensional two-phase flow calculations

International Nuclear Information System (INIS)

Toumi, I.; Caruge, D.

1997-01-01

This paper is devoted to new numerical methods developed for both one and three dimensional two-phase flow calculations. These methods are finite volume numerical methods and are based on the use of Approximate Riemann Solvers concepts to define convective fluxes versus mean cell quantities. The first part of the paper presents the numerical method for a one dimensional hyperbolic two-fluid model including differential terms as added mass and interface pressure. This numerical solution scheme makes use of the Riemann problem solution to define backward and forward differencing to approximate spatial derivatives. The construction of this approximate Riemann solver uses an extension of Roe's method that has been successfully used to solve gas dynamic equations. As far as the two-fluid model is hyperbolic, this numerical method seems very efficient for the numerical solution of two-phase flow problems. The scheme was applied both to shock tube problems and to standard tests for two-fluid computer codes. The second part describes the numerical method in the three dimensional case. The authors discuss also some improvements performed to obtain a fully implicit solution method that provides fast running steady state calculations. Such a scheme is not implemented in a thermal-hydraulic computer code devoted to 3-D steady-state and transient computations. Some results obtained for Pressurised Water Reactors concerning upper plenum calculations and a steady state flow in the core with rod bow effect evaluation are presented. In practice these new numerical methods have proved to be stable on non staggered grids and capable of generating accurate non oscillating solutions for two-phase flow calculations

Development op finite volume methods for fluid dynamics

International Nuclear Information System (INIS)

Delcourte, S.

2007-09-01

We aim to develop a finite volume method which applies to a greater class of meshes than other finite volume methods, restricted by orthogonality constraints. We build discrete differential operators over the three staggered tessellations needed for the construction of the method. These operators verify some analogous properties to those of the continuous operators. At first, the method is applied to the Div-Curl problem, which can be viewed as a building block of the Stokes problem. Then, the Stokes problem is dealt with with various boundary conditions. It is well known that when the computational domain is polygonal and non-convex, the order of convergence of numerical methods is deteriorated. Consequently, we have studied how an appropriate local refinement is able to restore the optimal order of convergence for the Laplacian problem. At last, we have discretized the non-linear Navier-Stokes problem, using the rotational formulation of the convection term, associated to the Bernoulli pressure. With an iterative algorithm, we are led to solve a saddle-point problem at each iteration. We give a particular interest to this linear problem by testing some pre-conditioners issued from finite elements, which we adapt to our method. Each problem is illustrated by numerical results on arbitrary meshes, such as strongly non-conforming meshes. (author)

Methods of numerical relativity

International Nuclear Information System (INIS)

Piran, T.

1983-01-01

Numerical Relativity is an alternative to analytical methods for obtaining solutions for Einstein equations. Numerical methods are particularly useful for studying generation of gravitational radiation by potential strong sources. The author reviews the analytical background, the numerical analysis aspects and techniques and some of the difficulties involved in numerical relativity. (Auth.)

Simulation of Jetting in Injection Molding Using a Finite Volume Method

Directory of Open Access Journals (Sweden)

Shaozhen Hua

2016-05-01

Full Text Available In order to predict the jetting and the subsequent buckling flow more accurately, a three dimensional melt flow model was established on a viscous, incompressible, and non-isothermal fluid, and a control volume-based finite volume method was employed to discretize the governing equations. A two-fold iterative method was proposed to decouple the dependence among pressure, velocity, and temperature so as to reduce the computation and improve the numerical stability. Based on the proposed theoretical model and numerical method, a program code was developed to simulate melt front progress and flow fields. The numerical simulations for different injection speeds, melt temperatures, and gate locations were carried out to explore the jetting mechanism. The results indicate the filling pattern depends on the competition between inertial and viscous forces. When inertial force exceeds the viscous force jetting occurs, then it changes to a buckling flow as the viscous force competes over the inertial force. Once the melt contacts with the mold wall, the melt filling switches to conventional sequential filling mode. Numerical results also indicate jetting length increases with injection speed but changes little with melt temperature. The reasonable agreements between simulated and experimental jetting length and buckling frequency imply the proposed method is valid for jetting simulation.

Numerical methods using Matlab

CERN Document Server

Lindfield, George

2012-01-01

Numerical Methods using MATLAB, 3e, is an extensive reference offering hundreds of useful and important numerical algorithms that can be implemented into MATLAB for a graphical interpretation to help researchers analyze a particular outcome. Many worked examples are given together with exercises and solutions to illustrate how numerical methods can be used to study problems that have applications in the biosciences, chaos, optimization, engineering and science across the board. Numerical Methods using MATLAB, 3e, is an extensive reference offering hundreds of use

Numerical method for estimating the size of chaotic regions of phase space

International Nuclear Information System (INIS)

Henyey, F.S.; Pomphrey, N.

1987-10-01

A numerical method for estimating irregular volumes of phase space is derived. The estimate weights the irregular area on a surface of section with the average return time to the section. We illustrate the method by application to the stadium and oval billiard systems and also apply the method to the continuous Henon-Heiles system. 15 refs., 10 figs

Finite Volume Method for Pricing European Call Option with Regime-switching Volatility

Science.gov (United States)

Lista Tauryawati, Mey; Imron, Chairul; Putri, Endah RM

2018-03-01

In this paper, we present a finite volume method for pricing European call option using Black-Scholes equation with regime-switching volatility. In the first step, we formulate the Black-Scholes equations with regime-switching volatility. we use a finite volume method based on fitted finite volume with spatial discretization and an implicit time stepping technique for the case. We show that the regime-switching scheme can revert to the non-switching Black Scholes equation, both in theoretical evidence and numerical simulations.

Advanced numerical methods for three dimensional two-phase flow calculations in PWR

International Nuclear Information System (INIS)

Toumi, I.; Gallo, D.; Royer, E.

1997-01-01

This paper is devoted to new numerical methods developed for three dimensional two-phase flow calculations. These methods are finite volume numerical methods. They are based on an extension of Roe's approximate Riemann solver to define convective fluxes versus mean cell quantities. To go forward in time, a linearized conservative implicit integrating step is used, together with a Newton iterative method. We also present here some improvements performed to obtain a fully implicit solution method that provides fast running steady state calculations. This kind of numerical method, which is widely used for fluid dynamic calculations, is proved to be very efficient for the numerical solution to two-phase flow problems. This numerical method has been implemented for the three dimensional thermal-hydraulic code FLICA-4 which is mainly dedicated to core thermal-hydraulic transient and steady-state analysis. Hereafter, we will also find some results obtained for the EPR reactor running in a steady-state at 60% of nominal power with 3 pumps out of 4, and a thermal-hydraulic core analysis for a 1300 MW PWR at low flow steam-line-break conditions. (author)

Ion diode simulation with a finite-volume PIC approach for the numerical solution of the Maxwell-Lorentz system

Energy Technology Data Exchange (ETDEWEB)

Munz, C D; Schneider, R; Stein, E; Voss, U [Forschungszentrum Karlsruhe (Germany). Institut fuer Neutronenphysik und Reaktortechnik; Westermann, T [FH Karlsruhe (Germany). Fachbereich Naturwissenschaften; Krauss, M [Forschungszentrum Karlsruhe (Germany). Hauptabteilung Informations- und Kommunikationstechik

1997-12-31

The numerical concept realized in the the Karlsruhe Diode Code KADI2D is briefly reviewed. Several new aspects concerning the Maxwell field solver based on high resolution finite-volume methods are presented. A new approach maintaining charge conservation numerically for the Maxwell-Lorentz equations is shortly summarized. (author). 2 figs., 12 refs.

Ion diode simulation with a finite-volume PIC approach for the numerical solution of the Maxwell-Lorentz system

International Nuclear Information System (INIS)

Munz, C.D.; Schneider, R.; Stein, E.; Voss, U.; Westermann, T.; Krauss, M.

1996-01-01

The numerical concept realized in the the Karlsruhe Diode Code KADI2D is briefly reviewed. Several new aspects concerning the Maxwell field solver based on high resolution finite-volume methods are presented. A new approach maintaining charge conservation numerically for the Maxwell-Lorentz equations is shortly summarized. (author). 2 figs., 12 refs

Flow simulation of a Pelton bucket using finite volume particle method

International Nuclear Information System (INIS)

Vessaz, C; Jahanbakhsh, E; Avellan, F

2014-01-01

The objective of the present paper is to perform an accurate numerical simulation of the high-speed water jet impinging on a Pelton bucket. To reach this goal, the Finite Volume Particle Method (FVPM) is used to discretize the governing equations. FVPM is an arbitrary Lagrangian-Eulerian method, which combines attractive features of Smoothed Particle Hydrodynamics and conventional mesh-based Finite Volume Method. This method is able to satisfy free surface and no-slip wall boundary conditions precisely. The fluid flow is assumed weakly compressible and the wall boundary is represented by one layer of particles located on the bucket surface. In the present study, the simulations of the flow in a stationary bucket are investigated for three different impinging angles: 72°, 90° and 108°. The particles resolution is first validated by a convergence study. Then, the FVPM results are validated with available experimental data and conventional grid-based Volume Of Fluid simulations. It is shown that the wall pressure field is in good agreement with the experimental and numerical data. Finally, the torque evolution and water sheet location are presented for a simulation of five rotating Pelton buckets

3rd International Conference on Numerical Analysis and Optimization : Theory, Methods, Applications and Technology Transfer

CERN Document Server

Grandinetti, Lucio; Purnama, Anton

2015-01-01

Presenting the latest findings in the field of numerical analysis and optimization, this volume balances pure research with practical applications of the subject. Accompanied by detailed tables, figures, and examinations of useful software tools, this volume will equip the reader to perform detailed and layered analysis of complex datasets. Many real-world complex problems can be formulated as optimization tasks. Such problems can be characterized as large scale, unconstrained, constrained, non-convex, non-differentiable, and discontinuous, and therefore require adequate computational methods, algorithms, and software tools. These same tools are often employed by researchers working in current IT hot topics such as big data, optimization and other complex numerical algorithms on the cloud, devising special techniques for supercomputing systems. The list of topics covered include, but are not limited to: numerical analysis, numerical optimization, numerical linear algebra, numerical differential equations, opt...

A spatial discretization of the MHD equations based on the finite volume - spectral method

International Nuclear Information System (INIS)

Miyoshi, Takahiro

2000-05-01

Based on the finite volume - spectral method, we present new discretization formulae for the spatial differential operators in the full system of the compressible MHD equations. In this approach, the cell-centered finite volume method is adopted in a bounded plane (poloidal plane), while the spectral method is applied to the differential with respect to the periodic direction perpendicular to the poloidal plane (toroidal direction). Here, an unstructured grid system composed of the arbitrary triangular elements is utilized for constructing the cell-centered finite volume method. In order to maintain the divergence free constraint of the magnetic field numerically, only the poloidal component of the rotation is defined at three edges of the triangular element. This poloidal component is evaluated under the assumption that the toroidal component of the operated vector times the radius, RA φ , is linearly distributed in the element. The present method will be applied to the nonlinear MHD dynamics in an realistic torus geometry without the numerical singularities. (author)

Numerical Methods for Free Boundary Problems

CERN Document Server

1991-01-01

About 80 participants from 16 countries attended the Conference on Numerical Methods for Free Boundary Problems, held at the University of Jyviiskylii, Finland, July 23-27, 1990. The main purpose of this conference was to provide up-to-date information on important directions of research in the field of free boundary problems and their numerical solutions. The contributions contained in this volume cover the lectures given in the conference. The invited lectures were given by H.W. Alt, V. Barbu, K-H. Hoffmann, H. Mittelmann and V. Rivkind. In his lecture H.W. Alt considered a mathematical model and existence theory for non-isothermal phase separations in binary systems. The lecture of V. Barbu was on the approximate solvability of the inverse one phase Stefan problem. K-H. Hoff­ mann gave an up-to-date survey of several directions in free boundary problems and listed several applications, but the material of his lecture is not included in this proceedings. H.D. Mittelmann handled the stability of thermo capi...

Numerical modeling of contaminant transport in fractured porous media using mixed finite-element and finitevolume methods

KAUST Repository

Dong, Chen; Sun, Shuyu; Taylor, Glenn A.

2011-01-01

A mathematical model for contaminant species passing through fractured porous media is presented. In the numerical model, we combine two locally conservative methods; i.e., the mixed finite-element (MFE) method and the finite-volume method. Adaptive

A New Class of Non-Linear, Finite-Volume Methods for Vlasov Simulation

International Nuclear Information System (INIS)

Banks, J.W.; Hittinger, J.A.

2010-01-01

Methods for the numerical discretization of the Vlasov equation should efficiently use the phase space discretization and should introduce only enough numerical dissipation to promote stability and control oscillations. A new high-order, non-linear, finite-volume algorithm for the Vlasov equation that discretely conserves particle number and controls oscillations is presented. The method is fourth-order in space and time in well-resolved regions, but smoothly reduces to a third-order upwind scheme as features become poorly resolved. The new scheme is applied to several standard problems for the Vlasov-Poisson system, and the results are compared with those from other finite-volume approaches, including an artificial viscosity scheme and the Piecewise Parabolic Method. It is shown that the new scheme is able to control oscillations while preserving a higher degree of fidelity of the solution than the other approaches.

Development of parallel implementation of adaptive numerical methods with industrial applications in fluid mechanics

International Nuclear Information System (INIS)

Laucoin, E.

2008-10-01

Numerical resolution of partial differential equations can be made reliable and efficient through the use of adaptive numerical methods.We present here the work we have done for the design, the implementation and the validation of such a method within an industrial software platform with applications in thermohydraulics. From the geometric point of view, this method can deal both with mesh refinement and mesh coarsening, while ensuring the quality of the mesh cells. Numerically, we use the mortar elements formalism in order to extend the Finite Volumes-Elements method implemented in the Trio-U platform and to deal with the non-conforming meshes arising from the adaptation procedure. Finally, we present an implementation of this method using concepts from domain decomposition methods for ensuring its efficiency while running in a parallel execution context. (author)

Numerical evaluation of acoustic characteristics and their damping of a thrust chamber using a constant-volume bomb model

Directory of Open Access Journals (Sweden)

Jianxiu QIN

2018-03-01

Full Text Available In order to numerically evaluate the acoustic characteristics of liquid rocket engine thrust chambers by means of a computational fluid dynamics method, a mathematical model of an artificial constant-volume bomb is proposed in this paper. A localized pressure pulse with a very high amplitude can be imposed on specified regions in a combustion chamber, the numerical procedure of which is described. Pressure oscillations actuated by the released constant-volume bomb can then be analyzed via Fast Fourier Transformation (FFT, and their modes can be identified according to the theoretical acoustic eigenfrequencies of the thrust chamber. The damping performances of the corresponding acoustic modes are evaluated by the half-power bandwidth method. The predicted acoustic characteristics and their damping for a special engine combustor agree well with the experimental data, validating the mathematical model and its numerical procedures. A small-thrust liquid rocket engine chamber is then analyzed by the present model. The First Longitudinal (1L acoustic mode can be excited easily and is hard to be damped. The axial position of the central constant-volume bomb has little influence on the amplitude and damping capacity of the First Radial (1R and 1L acoustic modes. Tangential acoustic modes can only be triggered by an off-centered constant-volume bomb, among which the First Tangential (1T mode is the strongest and regarded as the most harmful one. The amplitude of the 1L acoustic mode is smaller, but its damping factor is larger, as a constant-volume bomb is imposed approaching the injector face. These results are contributed to evaluate the acoustic characteristics and their damping of the combustion chamber. Keywords: Acoustic mode, Constant-volume bomb, Damping characteristics, Damping factor, Half-power bandwidth, Pressure oscillation

Introduction to precise numerical methods

CERN Document Server

Aberth, Oliver

2007-01-01

Precise numerical analysis may be defined as the study of computer methods for solving mathematical problems either exactly or to prescribed accuracy. This book explains how precise numerical analysis is constructed. The book also provides exercises which illustrate points from the text and references for the methods presented. All disc-based content for this title is now available on the Web. · Clearer, simpler descriptions and explanations ofthe various numerical methods· Two new types of numerical problems; accurately solving partial differential equations with the included software and computing line integrals in the complex plane.

New finite volume methods for approximating partial differential equations on arbitrary meshes

International Nuclear Information System (INIS)

Hermeline, F.

2008-12-01

This dissertation presents some new methods of finite volume type for approximating partial differential equations on arbitrary meshes. The main idea lies in solving twice the problem to be dealt with. One addresses the elliptic equations with variable (anisotropic, antisymmetric, discontinuous) coefficients, the parabolic linear or non linear equations (heat equation, radiative diffusion, magnetic diffusion with Hall effect), the wave type equations (Maxwell, acoustics), the elasticity and Stokes'equations. Numerous numerical experiments show the good behaviour of this type of method. (author)

A finite volume method for cylindrical heat conduction problems based on local analytical solution

KAUST Repository

Li, Wang

2012-10-01

A new finite volume method for cylindrical heat conduction problems based on local analytical solution is proposed in this paper with detailed derivation. The calculation results of this new method are compared with the traditional second-order finite volume method. The newly proposed method is more accurate than conventional ones, even though the discretized expression of this proposed method is slightly more complex than the second-order central finite volume method, making it cost more calculation time on the same grids. Numerical result shows that the total CPU time of the new method is significantly less than conventional methods for achieving the same level of accuracy. © 2012 Elsevier Ltd. All rights reserved.

A finite volume method for cylindrical heat conduction problems based on local analytical solution

KAUST Repository

Li, Wang; Yu, Bo; Wang, Xinran; Wang, Peng; Sun, Shuyu

2012-01-01

A new finite volume method for cylindrical heat conduction problems based on local analytical solution is proposed in this paper with detailed derivation. The calculation results of this new method are compared with the traditional second-order finite volume method. The newly proposed method is more accurate than conventional ones, even though the discretized expression of this proposed method is slightly more complex than the second-order central finite volume method, making it cost more calculation time on the same grids. Numerical result shows that the total CPU time of the new method is significantly less than conventional methods for achieving the same level of accuracy. © 2012 Elsevier Ltd. All rights reserved.

Numerical methods for the simulation of particle generated electromagnetic fields in acclerator physics

International Nuclear Information System (INIS)

Lau, T.

2006-01-01

In this work modifications of the classical Particle-In-Cell method for the solution of the Maxwell-Vlasov equations are investigated with respect to their application in particle accelerator physics. The aim of the work is to find modifications of the method which minimize and under certain conditions even eliminate the numerical dispersion effect along the beam axis in the numerical solution of Maxwell's equations. This is achieved by the development of dedicated time-integration methods for the Finite Integration Technique and two Finite Volume Methods. The methods are theoretically investigated regarding the conservation of a discrete energy and the existence of a discrete continuity equation. Finally, some of the methods are applied to the simulation of a high frequency rf-gun. (orig.)

Numerical solution of viscous and viscoelastic fluids flow through the branching channel by finite volume scheme

Science.gov (United States)

Keslerová, Radka; Trdlička, David

2015-09-01

This work deals with the numerical modelling of steady flows of incompressible viscous and viscoelastic fluids through the three dimensional channel with T-junction. The fundamental system of equations is the system of generalized Navier-Stokes equations for incompressible fluids. This system is based on the system of balance laws of mass and momentum for incompressible fluids. Two different mathematical models for the stress tensor are used for simulation of Newtonian and Oldroyd-B fluids flow. Numerical solution of the described models is based on cetral finite volume method using explicit Runge-Kutta time integration.

Numerical modeling of contaminant transport in fractured porous media using mixed finite-element and finitevolume methods

KAUST Repository

Dong, Chen

2011-01-01

A mathematical model for contaminant species passing through fractured porous media is presented. In the numerical model, we combine two locally conservative methods; i.e., the mixed finite-element (MFE) method and the finite-volume method. Adaptive triangle mesh is used for effective treatment of the fractures. A hybrid MFE method is employed to provide an accurate approximation of velocity fields for both the fractures and matrix, which are crucial to the convection part of the transport equation. The finite-volume method and the standard MFE method are used to approximate the convection and dispersion terms, respectively. The temporary evolution for the pressure distributions, streamline fields, and concentration profiles are obtained for six different arrangements of fractures. The results clearly show the distorted concentration effects caused by the ordered and disordered (random) patterns of the fractures and illustrate the robustness and efficiency of the proposed numerical model. © 2011 by Begell House Inc.

Error estimates for a numerical method for the compressible Navier-Stokes system on sufficiently smooth domains

Czech Academy of Sciences Publication Activity Database

Feireisl, Eduard; Hošek, Radim; Maltese, D.; Novotný, A.

2017-01-01

Roč. 51, č. 1 (2017), s. 279-319 ISSN 0764-583X EU Projects: European Commission(XE) 320078 - MATHEF Institutional support: RVO:67985840 Keywords : Navier-Stokes system * finite element numerical method * finite volume numerical method Subject RIV: BA - General Mathematics OBOR OECD: Pure mathematics Impact factor: 1.727, year: 2016 http://www.esaim-m2an.org/ articles /m2an/abs/2017/01/m2an150157/m2an150157.html

Strongly correlated systems numerical methods

CERN Document Server

Mancini, Ferdinando

2013-01-01

This volume presents, for the very first time, an exhaustive collection of those modern numerical methods specifically tailored for the analysis of Strongly Correlated Systems. Many novel materials, with functional properties emerging from macroscopic quantum behaviors at the frontier of modern research in physics, chemistry and material science, belong to this class of systems. Any technique is presented in great detail by its own inventor or by one of the world-wide recognized main contributors. The exposition has a clear pedagogical cut and fully reports on the most relevant case study where the specific technique showed to be very successful in describing and enlightening the puzzling physics of a particular strongly correlated system. The book is intended for advanced graduate students and post-docs in the field as textbook and/or main reference, but also for other researchers in the field who appreciate consulting a single, but comprehensive, source or wishes to get acquainted, in a as painless as possi...

Continuous modelling study of numerical volumes - Applications to the visualization of anatomical structures

International Nuclear Information System (INIS)

Goret, C.

1990-12-01

Several technics of imaging (IRM, image scanners, tomoscintigraphy, echography) give numerical informations presented by means of a stack of parallel cross-sectional images. Since many years, 3-D mathematical tools have been developed and allow the 3 D images synthesis of surfaces. In first part, we give the technics of numerical volume exploitation and their medical applications to diagnosis and therapy. The second part is about a continuous modelling of the volume with a tensor product of cubic splines. We study the characteristics of this representation and its clinical validation. Finally, we treat of the problem of surface visualization of objects contained in the volume. The results show the interest of this model and allow to propose specifications for 3-D workstation realization [fr

Study on numerical methods for transient flow induced by speed-changing impeller of fluid machinery

International Nuclear Information System (INIS)

Wu, Dazhuan; Chen, Tao; Wang, Leqin; Cheng, Wentao; Sun, Youbo

2013-01-01

In order to establish a reliable numerical method for solving the transient rotating flow induced by a speed-changing impeller, two numerical methods based on finite volume method (FVM) were presented and analyzed in this study. Two-dimensional numerical simulations of incompressible transient unsteady flow induced by an impeller during starting process were carried out respectively by using DM and DSR methods. The accuracy and adaptability of the two methods were evaluated by comprehensively comparing the calculation results. Moreover, an intensive study on the application of DSR method was conducted subsequently. The results showed that transient flow structure evolution and transient characteristics of the starting impeller are obviously affected by the starting process. The transient flow can be captured by both two methods, and the DSR method shows a higher computational efficiency. As an application example, the starting process of a mixed-flow pump was simulated by using DSR method. The calculation results were analyzed by comparing with the experiment data.

A lattice Boltzmann coupled to finite volumes method for solving phase change problems

Directory of Open Access Journals (Sweden)

El Ganaoui Mohammed

2009-01-01

Full Text Available A numerical scheme coupling lattice Boltzmann and finite volumes approaches has been developed and qualified for test cases of phase change problems. In this work, the coupled partial differential equations of momentum conservation equations are solved with a non uniform lattice Boltzmann method. The energy equation is discretized by using a finite volume method. Simulations show the ability of this developed hybrid method to model the effects of convection, and to predict transfers. Benchmarking is operated both for conductive and convective situation dominating solid/liquid transition. Comparisons are achieved with respect to available analytical solutions and experimental results.

Numerical methods in software and analysis

CERN Document Server

Rice, John R

1992-01-01

Numerical Methods, Software, and Analysis, Second Edition introduces science and engineering students to the methods, tools, and ideas of numerical computation. Introductory courses in numerical methods face a fundamental problem-there is too little time to learn too much. This text solves that problem by using high-quality mathematical software. In fact, the objective of the text is to present scientific problem solving using standard mathematical software. This book discusses numerous programs and software packages focusing on the IMSL library (including the PROTRAN system) and ACM Algorithm

Excel spreadsheet in teaching numerical methods

Science.gov (United States)

Djamila, Harimi

2017-09-01

One of the important objectives in teaching numerical methods for undergraduates’ students is to bring into the comprehension of numerical methods algorithms. Although, manual calculation is important in understanding the procedure, it is time consuming and prone to error. This is specifically the case when considering the iteration procedure used in many numerical methods. Currently, many commercial programs are useful in teaching numerical methods such as Matlab, Maple, and Mathematica. These are usually not user-friendly by the uninitiated. Excel spreadsheet offers an initial level of programming, which it can be used either in or off campus. The students will not be distracted with writing codes. It must be emphasized that general commercial software is required to be introduced later to more elaborated questions. This article aims to report on a teaching numerical methods strategy for undergraduates engineering programs. It is directed to students, lecturers and researchers in engineering field.

Design of a micro-irrigation system based on the control volume method

Directory of Open Access Journals (Sweden)

Chasseriaux G.

2006-01-01

Full Text Available A micro-irrigation system design based on control volume method using the back step procedure is presented in this study. The proposed numerical method is simple and consists of delimiting an elementary volume of the lateral equipped with an emitter, called « control volume » on which the conservation equations of the fl uid hydrodynamicʼs are applied. Control volume method is an iterative method to calculate velocity and pressure step by step throughout the micro-irrigation network based on an assumed pressure at the end of the line. A simple microcomputer program was used for the calculation and the convergence was very fast. When the average water requirement of plants was estimated, it is easy to choose the sum of the average emitter discharge as the total average fl ow rate of the network. The design consists of exploring an economical and effi cient network to deliver uniformly the input fl ow rate for all emitters. This program permitted the design of a large complex network of thousands of emitters very quickly. Three subroutine programs calculate velocity and pressure at a lateral pipe and submain pipe. The control volume method has already been tested for lateral design, the results from which were validated by other methods as fi nite element method, so it permits to determine the optimal design for such micro-irrigation network

Numerical methods for hydrodynamic stability problems

International Nuclear Information System (INIS)

Fujimura, Kaoru

1985-11-01

Numerical methods for solving the Orr-Sommerfeld equation, which is the fundamental equation of the hydrodynamic stability theory for various shear flows, are reviewed and typical numerical results are presented. The methods of asymptotic solution, finite difference methods, initial value methods and expansions in orthogonal functions are compared. (author)

Numerical methods for Eulerian and Lagrangian conservation laws

CERN Document Server

Després, Bruno

2017-01-01

This book focuses on the interplay between Eulerian and Lagrangian conservation laws for systems that admit physical motivation and originate from continuum mechanics. Ultimately, it highlights what is specific to and beneficial in the Lagrangian approach and its numerical methods. The two first chapters present a selection of well-known features of conservation laws and prepare readers for the subsequent chapters, which are dedicated to the analysis and discretization of Lagrangian systems. The text is at the frontier of applied mathematics and scientific computing and appeals to students and researchers interested in Lagrangian-based computational fluid dynamics. It also serves as an introduction to the recent corner-based Lagrangian finite volume techniques.

Three-Dimensional Phase Field Simulations of Hysteresis and Butterfly Loops by the Finite Volume Method

International Nuclear Information System (INIS)

Xi Li-Ying; Chen Huan-Ming; Zheng Fu; Gao Hua; Tong Yang; Ma Zhi

2015-01-01

Three-dimensional simulations of ferroelectric hysteresis and butterfly loops are carried out based on solving the time dependent Ginzburg–Landau equations using a finite volume method. The influence of externally mechanical loadings with a tensile strain and a compressive strain on the hysteresis and butterfly loops is studied numerically. Different from the traditional finite element and finite difference methods, the finite volume method is applicable to simulate the ferroelectric phase transitions and properties of ferroelectric materials even for more realistic and physical problems. (paper)

HYDRA-II: A hydrothermal analysis computer code: Volume 1, Equations and numerics

International Nuclear Information System (INIS)

McCann, R.A.

1987-04-01

HYDRA-II is a hydrothermal computer code capable of three-dimensional analysis of coupled conduction, convection, and thermal radiation problems. This code is especially appropriate for simulating the steady-state performance of spent fuel storage systems. The code has been evaluated for this application for the US Department of Energy's Commercial Spent Fuel Management Program. HYDRA-II provides a finite difference solution in Cartesian coordinates to the equations governing the conservation of mass, momentum, and energy. A cylindrical coordinate system may also be used to enclose the Cartesian coordinate system. This exterior coordinate system is useful for modeling cylindrical cask bodies. The difference equations for conservation of momentum are enhanced by the incorporation of directional porosities and permeabilities that aid in modeling solid structures whose dimensions may be smaller than the computational mesh. The equation for conservation of energy permits of modeling of orthotropic physical properties and film resistances. Several automated procedures are available to model radiation transfer within enclosures and from fuel rod to fuel rod. The documentation of HYDRA-II is presented in three separate volumes. This volume, Volume I - Equations and Numerics, describes the basic differential equations, illustrates how the difference equations are formulated, and gives the solution procedures employed. Volume II - User's Manual contains code flow charts, discusses the code structure, provides detailed instructions for preparing an input file, and illustrates the operation of the code by means of a model problem. The final volume, Volume III - Verification/Validation Assessments, presents results of numerical simulations of single- and multiassembly storage systems and comparisons with experimental data. 4 refs

Advances in Numerical Methods

CERN Document Server

Mastorakis, Nikos E

2009-01-01

Features contributions that are focused on significant aspects of current numerical methods and computational mathematics. This book carries chapters that advanced methods and various variations on known techniques that can solve difficult scientific problems efficiently.

Some numerical methods for two-fluid two-phase flows in oil pipes; Quelques methodes numeriques pour les ecoulements diphasiques bi-fluide en conduites petrolieres

Energy Technology Data Exchange (ETDEWEB)

Masella, J.M.

1997-05-29

This thesis is devoted to the numerical simulation of some two-fluid models describing gas-liquid two-phase flow in pipes. The numerical models developed here can be more generally used in the modelling of a wide class of physical models which can be put under an hyperbolic form. We introduce first two isothermal two-fluid models, composed of a mass balance equation and a momentum equation written in each phase, describing respectively a stratified two-phase flow and a dispersed two-phase flow. These models are hyperbolic under some physical assumptions and can be written under a nonconservative vectorial system. We define and analyse a new numerical finite volume scheme (v{integral}Roe) founded on a linearized Riemann solver. This scheme does not need any analytical calculation and gives good results in the tracking of shocks. We compare this new scheme with the classical Roe scheme. Then we propose and study some numerical models, with and without flux splitting method, which are adapted to the discretization of the two-fluid models. This numerical models are given by a finite volume integration of the equations, and lean on the v{integral} scheme. In order to reducing cpu time, due to the low Mach number of two-phase flows, acoustic waves are implicit. Afterwards we proposed a discretization of boundary conditions, which allows the generation of transient flows in pipe. Some numerical academic and more physical tests show the good behaviour of the numerical methods. (author) 77 refs.

Extrusion Process by Finite Volume Method Using OpenFoam Software

International Nuclear Information System (INIS)

Matos Martins, Marcelo; Tonini Button, Sergio; Divo Bressan, Jose; Ivankovic, Alojz

2011-01-01

The computational codes are very important tools to solve engineering problems. In the analysis of metal forming process, such as extrusion, this is not different because the computational codes allow analyzing the process with reduced cost. Traditionally, the Finite Element Method is used to solve solid mechanic problems, however, the Finite Volume Method (FVM) have been gaining force in this field of applications. This paper presents the velocity field and friction coefficient variation results, obtained by numerical simulation using the OpenFoam Software and the FVM to solve an aluminum direct cold extrusion process.

Numerical Analysis of Indoor Sound Quality Evaluation Using Finite Element Method

Directory of Open Access Journals (Sweden)

Yu-Tuan Chou

2013-01-01

Full Text Available Indoors sound field distribution is important to Room Acoustics, but the field suffers numerous problems, for example, multipath propagation and scattering owing to sound absorption by furniture and other aspects of décor. Generally, an ideal interior space must have a sound field with clear quality. This provides both the speaker and the listener with a pleasant conversational environment. This investigation uses the Finite Element Method to assess the acoustic distribution based on the indoor space and chamber volume. In this situation, a fixed sound source at different frequencies is used to simulate the acoustic characteristics of the indoor space. This method considers the furniture and decoration sound absorbing material and thus different sound absorption coefficients and configurations. The preliminary numerical simulation provides a method that can forecast the distribution of sound in an indoor room in complex situations. Consequently, it is possible to arrange interior furnishings and appliances to optimize acoustic distribution and environmental friendliness. Additionally, the analytical results can also be used to calculate the Reverberation Time and speech intelligibility for specified indoor space.

Numerical investigation on compressible flow characteristics in axial compressors using a multi block finite-volume scheme

International Nuclear Information System (INIS)

Farhanieh, B.; Amanifard, N.; Ghorbanian, K.

2002-01-01

An unsteady two-dimensional numerical investigation was performed on the viscous flow passing through a multi-blade cascade. A Cartesian finite-volume approach was linked to Van-Leer's and Roe's flux splitting schemes to evaluate inviscid flux terms. To prevent the oscillatory behavior of numerical results and to increase the accuracy, Mon tonic Upstream Scheme for Conservation Laws was added to flux splitting schemes. The Baldwin-Lo max (B L) turbulence model was implemented to solve the turbulent case studies. Implicit solution was also provided using Lower and Upper (L U) decomposition technique to compare with explicit solutions. To validate the numerical procedure, two test cases are prepared and flow over a Na Ca 0012 airfoil was investigated and the pressure coefficients were compared to the reference data. The numerical solver was implemented to study the flow passing over a compressor cascade. The results of various combinations of splitting schemes and the Mon tonic Upstream Scheme for Conventional Laws limiter were compared with each other to find the suitable methods in cascade problems. Finally the convergence histories of implemented schemes were compared to each other to show the behavior of the solver in using various methods before implementation of them in flow instability studies

Simple method for the generation of multiple homogeneous field volumes inside the bore of superconducting magnets.

Science.gov (United States)

Chou, Ching-Yu; Ferrage, Fabien; Aubert, Guy; Sakellariou, Dimitris

2015-07-17

Standard Magnetic Resonance magnets produce a single homogeneous field volume, where the analysis is performed. Nonetheless, several modern applications could benefit from the generation of multiple homogeneous field volumes along the axis and inside the bore of the magnet. In this communication, we propose a straightforward method using a combination of ring structures of permanent magnets in order to cancel the gradient of the stray field in a series of distinct volumes. These concepts were demonstrated numerically on an experimentally measured magnetic field profile. We discuss advantages and limitations of our method and present the key steps required for an experimental validation.

Numerical simulation of seismic wave propagation from land-excited large volume air-gun source

Science.gov (United States)

Cao, W.; Zhang, W.

2017-12-01

The land-excited large volume air-gun source can be used to study regional underground structures and to detect temporal velocity changes. The air-gun source is characterized by rich low frequency energy (from bubble oscillation, 2-8Hz) and high repeatability. It can be excited in rivers, reservoirs or man-made pool. Numerical simulation of the seismic wave propagation from the air-gun source helps to understand the energy partitioning and characteristics of the waveform records at stations. However, the effective energy recorded at a distance station is from the process of bubble oscillation, which can not be approximated by a single point source. We propose a method to simulate the seismic wave propagation from the land-excited large volume air-gun source by finite difference method. The process can be divided into three parts: bubble oscillation and source coupling, solid-fluid coupling and the propagation in the solid medium. For the first part, the wavelet of the bubble oscillation can be simulated by bubble model. We use wave injection method combining the bubble wavelet with elastic wave equation to achieve the source coupling. Then, the solid-fluid boundary condition is implemented along the water bottom. And the last part is the seismic wave propagation in the solid medium, which can be readily implemented by the finite difference method. Our method can get accuracy waveform of land-excited large volume air-gun source. Based on the above forward modeling technology, we analysis the effect of the excited P wave and the energy of converted S wave due to different water shapes. We study two land-excited large volume air-gun fields, one is Binchuan in Yunnan, and the other is Hutubi in Xinjiang. The station in Binchuan, Yunnan is located in a large irregular reservoir, the waveform records have a clear S wave. Nevertheless, the station in Hutubi, Xinjiang is located in a small man-made pool, the waveform records have very weak S wave. Better understanding of

Numerical Methods for Stochastic Computations A Spectral Method Approach

CERN Document Server

Xiu, Dongbin

2010-01-01

The first graduate-level textbook to focus on fundamental aspects of numerical methods for stochastic computations, this book describes the class of numerical methods based on generalized polynomial chaos (gPC). These fast, efficient, and accurate methods are an extension of the classical spectral methods of high-dimensional random spaces. Designed to simulate complex systems subject to random inputs, these methods are widely used in many areas of computer science and engineering. The book introduces polynomial approximation theory and probability theory; describes the basic theory of gPC meth

A novel finite volume discretization method for advection-diffusion systems on stretched meshes

Science.gov (United States)

Merrick, D. G.; Malan, A. G.; van Rooyen, J. A.

2018-06-01

This work is concerned with spatial advection and diffusion discretization technology within the field of Computational Fluid Dynamics (CFD). In this context, a novel method is proposed, which is dubbed the Enhanced Taylor Advection-Diffusion (ETAD) scheme. The model equation employed for design of the scheme is the scalar advection-diffusion equation, the industrial application being incompressible laminar and turbulent flow. Developed to be implementable into finite volume codes, ETAD places specific emphasis on improving accuracy on stretched structured and unstructured meshes while considering both advection and diffusion aspects in a holistic manner. A vertex-centered structured and unstructured finite volume scheme is used, and only data available on either side of the volume face is employed. This includes the addition of a so-called mesh stretching metric. Additionally, non-linear blending with the existing NVSF scheme was performed in the interest of robustness and stability, particularly on equispaced meshes. The developed scheme is assessed in terms of accuracy - this is done analytically and numerically, via comparison to upwind methods which include the popular QUICK and CUI techniques. Numerical tests involved the 1D scalar advection-diffusion equation, a 2D lid driven cavity and turbulent flow case. Significant improvements in accuracy were achieved, with L2 error reductions of up to 75%.

Numerical computations with GPUs

CERN Document Server

Kindratenko, Volodymyr

2014-01-01

This book brings together research on numerical methods adapted for Graphics Processing Units (GPUs). It explains recent efforts to adapt classic numerical methods, including solution of linear equations and FFT, for massively parallel GPU architectures. This volume consolidates recent research and adaptations, covering widely used methods that are at the core of many scientific and engineering computations. Each chapter is written by authors working on a specific group of methods; these leading experts provide mathematical background, parallel algorithms and implementation details leading to

Constitutive Modelling in Thermomechanical Processes, Using The Control Volume Method on Staggered Grid

DEFF Research Database (Denmark)

Thorborg, Jesper

, however, is constituted by the implementation of the $J_2$ flow theory in the control volume method. To apply the control volume formulation on the process of hardening concrete viscoelastic stress-strain models has been examined in terms of various rheological models. The generalized 3D models are based...... on two different suggestions in the literature, that is compressible or incompressible behaviour of the viscos response in the dashpot element. Numerical implementation of the models has shown very good agreement with corresponding analytical solutions. The viscoelastic solid mechanical model is used...

An outline review of numerical transport methods

International Nuclear Information System (INIS)

Budd, C.

1981-01-01

A brief review is presented of numerical methods for solving the neutron transport equation in the context of reactor physics. First the various forms of transport equation are given. Second, the various ways of classifying numerical transport methods are discussed. Finally each method (or class of methods) is outlined in turn. (U.K.)

Coupled numerical approach combining finite volume and lattice Boltzmann methods for multi-scale multi-physicochemical processes

Energy Technology Data Exchange (ETDEWEB)

Chen, Li; He, Ya-Ling [Key Laboratory of Thermo-Fluid Science and Engineering of MOE, School of Energy and Power Engineering, Xi' an Jiaotong University, Xi' an, Shaanxi 710049 (China); Kang, Qinjun [Computational Earth Science Group (EES-16), Los Alamos National Laboratory, Los Alamos, NM (United States); Tao, Wen-Quan, E-mail: wqtao@mail.xjtu.edu.cn [Key Laboratory of Thermo-Fluid Science and Engineering of MOE, School of Energy and Power Engineering, Xi' an Jiaotong University, Xi' an, Shaanxi 710049 (China)

2013-12-15

A coupled (hybrid) simulation strategy spatially combining the finite volume method (FVM) and the lattice Boltzmann method (LBM), called CFVLBM, is developed to simulate coupled multi-scale multi-physicochemical processes. In the CFVLBM, computational domain of multi-scale problems is divided into two sub-domains, i.e., an open, free fluid region and a region filled with porous materials. The FVM and LBM are used for these two regions, respectively, with information exchanged at the interface between the two sub-domains. A general reconstruction operator (RO) is proposed to derive the distribution functions in the LBM from the corresponding macro scalar, the governing equation of which obeys the convection–diffusion equation. The CFVLBM and the RO are validated in several typical physicochemical problems and then are applied to simulate complex multi-scale coupled fluid flow, heat transfer, mass transport, and chemical reaction in a wall-coated micro reactor. The maximum ratio of the grid size between the FVM and LBM regions is explored and discussed. -- Highlights: •A coupled simulation strategy for simulating multi-scale phenomena is developed. •Finite volume method and lattice Boltzmann method are coupled. •A reconstruction operator is derived to transfer information at the sub-domains interface. •Coupled multi-scale multiple physicochemical processes in micro reactor are simulated. •Techniques to save computational resources and improve the efficiency are discussed.

Numerical Modelling of Three-Fluid Flow Using The Level-set Method

Science.gov (United States)

Li, Hongying; Lou, Jing; Shang, Zhi

2014-11-01

This work presents a numerical model for simulation of three-fluid flow involving two different moving interfaces. These interfaces are captured using the level-set method via two different level-set functions. A combined formulation with only one set of conservation equations for the whole physical domain, consisting of the three different immiscible fluids, is employed. Numerical solution is performed on a fixed mesh using the finite volume method. Surface tension effect is incorporated using the Continuum Surface Force model. Validation of the present model is made against available results for stratified flow and rising bubble in a container with a free surface. Applications of the present model are demonstrated by a variety of three-fluid flow systems including (1) three-fluid stratified flow, (2) two-fluid stratified flow carrying the third fluid in the form of drops and (3) simultaneous rising and settling of two drops in a stationary third fluid. The work is supported by a Thematic and Strategic Research from A*STAR, Singapore (Ref. #: 1021640075).

The finite volume method in computational fluid dynamics an advanced introduction with OpenFOAM and Matlab

CERN Document Server

Moukalled, F; Darwish, M

2016-01-01

This textbook explores both the theoretical foundation of the Finite Volume Method (FVM) and its applications in Computational Fluid Dynamics (CFD). Readers will discover a thorough explanation of the FVM numerics and algorithms used for the simulation of incompressible and compressible fluid flows, along with a detailed examination of the components needed for the development of a collocated unstructured pressure-based CFD solver. Two particular CFD codes are explored. The first is uFVM, a three-dimensional unstructured pressure-based finite volume academic CFD code, implemented within Matlab. The second is OpenFOAM®, an open source framework used in the development of a range of CFD programs for the simulation of industrial scale flow problems. With over 220 figures, numerous examples and more than one hundred exercise on FVM numerics, programming, and applications, this textbook is suitable for use in an introductory course on the FVM, in an advanced course on numerics, and as a reference for CFD programm...

On nitrogen condensation in hypersonic nozzle flows: Numerical method and parametric study

KAUST Repository

Lin, Longyuan

2013-12-17

A numerical method for calculating two-dimensional planar and axisymmetric hypersonic nozzle flows with nitrogen condensation is developed. The classical nucleation theory with an empirical correction function and the modified Gyarmathy model are used to describe the nucleation rate and the droplet growth, respectively. The conservation of the liquid phase is described by a finite number of moments of the size distribution function. The moment equations are then combined with the Euler equations and are solved by the finite-volume method. The numerical method is first validated by comparing its prediction with experimental results from the literature. The effects of nitrogen condensation on hypersonic nozzle flows are then numerically examined. The parameters at the nozzle exit under the conditions of condensation and no-condensation are evaluated. For the condensation case, the static pressure, the static temperature, and the amount of condensed fluid at the nozzle exit decrease with the increase of the total temperature. Compared with the no-condensation case, both the static pressure and temperature at the nozzle exit increase, and the Mach number decreases due to the nitrogen condensation. It is also indicated that preheating the nitrogen gas is necessary to avoid the nitrogen condensation even for a hypersonic nozzle with a Mach number of 5 operating at room temperatures. © 2013 Springer-Verlag Berlin Heidelberg.

Numerical methods and modelling for engineering

CERN Document Server

Khoury, Richard

2016-01-01

This textbook provides a step-by-step approach to numerical methods in engineering modelling. The authors provide a consistent treatment of the topic, from the ground up, to reinforce for students that numerical methods are a set of mathematical modelling tools which allow engineers to represent real-world systems and compute features of these systems with a predictable error rate. Each method presented addresses a specific type of problem, namely root-finding, optimization, integral, derivative, initial value problem, or boundary value problem, and each one encompasses a set of algorithms to solve the problem given some information and to a known error bound. The authors demonstrate that after developing a proper model and understanding of the engineering situation they are working on, engineers can break down a model into a set of specific mathematical problems, and then implement the appropriate numerical methods to solve these problems. Uses a “building-block” approach, starting with simpler mathemati...

Cálculo do volume na equação de van der Waals pelo método de cardano Volume calculation in van der Waals equation by the cardano method

Directory of Open Access Journals (Sweden)

Nelson H. T. Lemes

2010-01-01

Full Text Available Analytical solutions of a cubic equation with real coefficients are established using the Cardano method. The method is first applied to simple third order equation. Calculation of volume in the van der Waals equation of state is afterwards established. These results are exemplified to calculate the volumes below and above critical temperatures. Analytical and numerical values for the compressibility factor are presented as a function of the pressure. As a final example, coexistence volumes in the liquid-vapor equilibrium are calculated. The Cardano approach is very simple to apply, requiring only elementary operations, indicating an attractive method to be used in teaching elementary thermodynamics.

Numerical analysis

CERN Document Server

Brezinski, C

2012-01-01

Numerical analysis has witnessed many significant developments in the 20th century. This book brings together 16 papers dealing with historical developments, survey papers and papers on recent trends in selected areas of numerical analysis, such as: approximation and interpolation, solution of linear systems and eigenvalue problems, iterative methods, quadrature rules, solution of ordinary-, partial- and integral equations. The papers are reprinted from the 7-volume project of the Journal of Computational and Applied Mathematics on '/homepage/sac/cam/na2000/index.html<

An introduction to numerical methods and analysis

CERN Document Server

Epperson, James F

2013-01-01

Praise for the First Edition "". . . outstandingly appealing with regard to its style, contents, considerations of requirements of practice, choice of examples, and exercises.""-Zentralblatt MATH "". . . carefully structured with many detailed worked examples.""-The Mathematical Gazette The Second Edition of the highly regarded An Introduction to Numerical Methods and Analysis provides a fully revised guide to numerical approximation. The book continues to be accessible and expertly guides readers through the many available techniques of numerical methods and analysis. An Introduction to

Numerical Methods for Radiation Magnetohydrodynamics in Astrophysics

Energy Technology Data Exchange (ETDEWEB)

Klein, R I; Stone, J M

2007-11-20

We describe numerical methods for solving the equations of radiation magnetohydrodynamics (MHD) for astrophysical fluid flow. Such methods are essential for the investigation of the time-dependent and multidimensional dynamics of a variety of astrophysical systems, although our particular interest is motivated by problems in star formation. Over the past few years, the authors have been members of two parallel code development efforts, and this review reflects that organization. In particular, we discuss numerical methods for MHD as implemented in the Athena code, and numerical methods for radiation hydrodynamics as implemented in the Orion code. We discuss the challenges introduced by the use of adaptive mesh refinement in both codes, as well as the most promising directions for future developments.

Numerical Methods for Radiation Magnetohydrodynamics in Astrophysics

International Nuclear Information System (INIS)

Klein, R I; Stone, J M

2007-01-01

We describe numerical methods for solving the equations of radiation magnetohydrodynamics (MHD) for astrophysical fluid flow. Such methods are essential for the investigation of the time-dependent and multidimensional dynamics of a variety of astrophysical systems, although our particular interest is motivated by problems in star formation. Over the past few years, the authors have been members of two parallel code development efforts, and this review reflects that organization. In particular, we discuss numerical methods for MHD as implemented in the Athena code, and numerical methods for radiation hydrodynamics as implemented in the Orion code. We discuss the challenges introduced by the use of adaptive mesh refinement in both codes, as well as the most promising directions for future developments

A numerical method for PCM-based pin fin heat sinks optimization

International Nuclear Information System (INIS)

Pakrouh, R.; Hosseini, M.J.; Ranjbar, A.A.; Bahrampoury, R.

2015-01-01

Highlights: • Optimization of PCM-based heat sink by using the Taguchi method. • Derivation of optimal PCM percentage to reach the maximum critical time. • Optimization is performed for four different critical temperatures. • Effective design factors are fins’ height and fins’ number. • The optimum configuration depends on geometric properties and the critical temperature. - Abstract: This paper presents a numerical investigation on geometric optimization of PCM-based pin fin heat sinks. Paraffin RT44HC is used as PCM while the fins and heat sink base is made of aluminum. The fins act as thermal conductivity enhancers (TCEs). The main goal of the study is to obtain the configurations that maximize the heat sink operational time. An approach witch couples Taguchi method with numerical simulations is utilized for this purpose. Number of fins, fins height, fins thickness and the base thickness are parameters which are studied for optimization. In this study natural convection and PCM volume variation during melting process are considered in the simulations. Optimization is performed for different critical temperatures of 50 °C, 60 °C, 70 °C and 80 °C. Results show that a complex relation exists between PCM and TCE volume percentages. The optimal case strongly depends on the fins’ number, fins’ height and thickness and also the critical temperature. The optimum PCM percentages are found to be 60.61% (corresponds to 100 pin fin heat sink with 4 mm thick fins) for critical temperature of 50 °C and 82.65% (corresponds to 100 pin fin heat sink with 2 mm thick fins) for other critical temperatures

ALE finite volume method for free-surface Bingham plastic fluids with general curvilinear coordinates

International Nuclear Information System (INIS)

Nagai, Katsuaki; Ushijima, Satoru

2010-01-01

A numerical prediction method has been proposed to predict Bingham plastic fluids with free-surface in a two-dimensional container. Since the linear relationships between stress tensors and strain rate tensors are not assumed for non-Newtonian fluids, the liquid motions are described with Cauchy momentum equations rather than Navier-Stokes equations. The profile of a liquid surface is represented with the two-dimensional curvilinear coordinates which are represented in each computational step on the basis of the arbitrary Lagrangian-Eulerian (ALE) method. Since the volumes of the fluid cells are transiently changed in the physical space, the geometric conservation law is applied to the finite volume discretizations. As a result, it has been shown that the present method enables us to predict reasonably the Bingham plastic fluids with free-surface in a container.

ALE finite volume method for free-surface Bingham plastic fluids with general curvilinear coordinates

Science.gov (United States)

Nagai, Katsuaki; Ushijima, Satoru

2010-06-01

A numerical prediction method has been proposed to predict Bingham plastic fluids with free-surface in a two-dimensional container. Since the linear relationships between stress tensors and strain rate tensors are not assumed for non-Newtonian fluids, the liquid motions are described with Cauchy momentum equations rather than Navier-Stokes equations. The profile of a liquid surface is represented with the two-dimensional curvilinear coordinates which are represented in each computational step on the basis of the arbitrary Lagrangian-Eulerian (ALE) method. Since the volumes of the fluid cells are transiently changed in the physical space, the geometric conservation law is applied to the finite volume discretizations. As a result, it has been shown that the present method enables us to predict reasonably the Bingham plastic fluids with free-surface in a container.

Proceedings of the fifth international symposium on numerical methods in engineering. Vol. 1 and 2

Energy Technology Data Exchange (ETDEWEB)

Gruber, R [Ecole Polytechnique Federale, Lausanne (Switzerland); Periaux, J [Avions Marcel Dassault-Breguet Aviation, 92 - Saint-Cloud (France). Aerodynamique, Methodes Numeriques; Shaw, R P [Buffalo Univ., NY (USA). Dept. of Civil Engineering; eds.

1989-01-01

The present two volumes survey the state of the art in advanced scientific computing as applied to engineering science. Many fields ranging from modelization (partial differential equations, integral equations, boundary conditions, macroscopic models, cellular automata) to numerical methods (finite elements, optimization, software engineering tools, parallel and vector computing methods, domain decomposition, multigrid) and applications (nonlinear solid mechanics, fracture mechanics, composite materials, friction and contact, fluid mechanics, chemical flows, convection, free boundaries, combustion, electromagnetics) are covered by invited papers, minisymposia and contributed papers. (orig./HP).

A relaxation-projection method for compressible flows. Part I: The numerical equation of state for the Euler equations

International Nuclear Information System (INIS)

Saurel, Richard; Franquet, Erwin; Daniel, Eric; Le Metayer, Olivier

2007-01-01

A new projection method is developed for the Euler equations to determine the thermodynamic state in computational cells. It consists in the resolution of a mechanical relaxation problem between the various sub-volumes present in a computational cell. These sub-volumes correspond to the ones traveled by the various waves that produce states with different pressures, velocities, densities and temperatures. Contrarily to Godunov type schemes the relaxed state corresponds to mechanical equilibrium only and remains out of thermal equilibrium. The pressure computation with this relaxation process replaces the use of the conventional equation of state (EOS). A simplified relaxation method is also derived and provides a specific EOS (named the Numerical EOS). The use of the Numerical EOS gives a cure to spurious pressure oscillations that appear at contact discontinuities for fluids governed by real gas EOS. It is then extended to the computation of interface problems separating fluids with different EOS (liquid-gas interface for example) with the Euler equations. The resulting method is very robust, accurate, oscillation free and conservative. For the sake of simplicity and efficiency the method is developed in a Lagrange-projection context and is validated over exact solutions. In a companion paper [F. Petitpas, E. Franquet, R. Saurel, A relaxation-projection method for compressible flows. Part II: computation of interfaces and multiphase mixtures with stiff mechanical relaxation. J. Comput. Phys. (submitted for publication)], the method is extended to the numerical approximation of a non-conservative hyperbolic multiphase flow model for interface computation and shock propagation into mixtures

Analysis of numerical methods

CERN Document Server

Isaacson, Eugene

1994-01-01

This excellent text for advanced undergraduates and graduate students covers norms, numerical solution of linear systems and matrix factoring, iterative solutions of nonlinear equations, eigenvalues and eigenvectors, polynomial approximation, and other topics. It offers a careful analysis and stresses techniques for developing new methods, plus many examples and problems. 1966 edition.

Hybrid methods for airframe noise numerical prediction

Energy Technology Data Exchange (ETDEWEB)

Terracol, M.; Manoha, E.; Herrero, C.; Labourasse, E.; Redonnet, S. [ONERA, Department of CFD and Aeroacoustics, BP 72, Chatillon (France); Sagaut, P. [Laboratoire de Modelisation en Mecanique - UPMC/CNRS, Paris (France)

2005-07-01

This paper describes some significant steps made towards the numerical simulation of the noise radiated by the high-lift devices of a plane. Since the full numerical simulation of such configuration is still out of reach for present supercomputers, some hybrid strategies have been developed to reduce the overall cost of such simulations. The proposed strategy relies on the coupling of an unsteady nearfield CFD with an acoustic propagation solver based on the resolution of the Euler equations for midfield propagation in an inhomogeneous field, and the use of an integral solver for farfield acoustic predictions. In the first part of this paper, this CFD/CAA coupling strategy is presented. In particular, the numerical method used in the propagation solver is detailed, and two applications of this coupling method to the numerical prediction of the aerodynamic noise of an airfoil are presented. Then, a hybrid RANS/LES method is proposed in order to perform some unsteady simulations of complex noise sources. This method allows for significant reduction of the cost of such a simulation by considerably reducing the extent of the LES zone. This method is described and some results of the numerical simulation of the three-dimensional unsteady flow in the slat cove of a high-lift profile are presented. While these results remain very difficult to validate with experiments on similar configurations, they represent up to now the first 3D computations of this kind of flow. (orig.)

Stratal Control Volumes and Stratal Control Trajectories: A New Method to Constrain, Understand and Reconcile Results from Stratigraphic Outcrop Analysis, Subsurface Analysis and Analogue and Numerical Modelling

Science.gov (United States)

Burgess, P. M.; Steel, R. J.

2016-12-01

Decoding a history of Earth's surface dynamics from strata requires robust quantitative understanding of supply and accommodation controls. The concept of stratigraphic solution sets has proven useful in this decoding, but application and development of this approach has so far been surprisingly limited. Stratal control volumes, areas and trajectories are new approaches defined here, building on previous ideas about stratigraphic solution sets, to help analyse and understand the sedimentary record of Earth surface dynamics. They may have particular application reconciling results from outcrop and subsurface analysis with results from analogue and numerical experiments. Stratal control volumes are sets of points in a three-dimensional volume, with axes of subsidence, sediment supply and eustatic rates of change, populated with probabilities derived from analysis of subsidence, supply and eustasy timeseries (Figure 1). These empirical probabilities indicate the likelihood of occurrence of any particular combination of control rates defined by any point in the volume. The stratal control volume can then by analysed to determine which parts of the volume represent relative sea-level fall and rise, where in the volume particular stacking patterns will occur, and how probable those stacking patterns are. For outcrop and subsurface analysis, using a stratal control area with eustasy and subsidence combined on a relative sea-level axis allows similar analysis, and may be preferable. A stratal control trajectory is a history of supply and accommodation creation rates, interpreted from outcrop or subsurface data, or observed in analogue and numerical experiments, and plotted as a series of linked points forming a trajectory through the stratal control volume (Figure 1) or area. Three examples are presented, one from outcrop and two theoretical. Much work remains to be done to build a properly representative database of stratal controls, but careful comparison of stratal

A novel numerical method for the analysis of electron transport in the presence of pointlike magnetic scatterers

International Nuclear Information System (INIS)

Miyagawa, Yuu; Ueta, Tsuyoshi

2008-01-01

The boundary element method (BEM) is so extended as to treat two-dimensional (2D) electron systems in the presence of pointlike islands of magnetic moment. In the present paper, the pointlike magnetic scatterer is modeled by a cylindrical barrier. The radius of the cylindrical barrier is assumed to be so small, keeping the volume definite, that the pointlike magnetic scatterer is approximated by a Dirac δ function. Then, we make an approximation on the BEM formulation, wherefore we derive a novel numerical method for electron transport in the presence of pointlike magnetic scatterers. In a numerical implementation of the method extended here, the numerical errors of probability conservation are less than 1% for any cases and the computational costs, that is, the required memory amount and CPU time, are much reduced. As examples, the proposed method is applied to transport problems through a quantum wire with four pointlike magnetic scatterers. It is clearly shown that magnetic scatterers, even pointlike magnetic moments, lead to spin flip-flop, localization and resonance

Numerical analysis in electromagnetics the TLM method

CERN Document Server

Saguet, Pierre

2013-01-01

The aim of this book is to give a broad overview of the TLM (Transmission Line Matrix) method, which is one of the "time-domain numerical methods". These methods are reputed for their significant reliance on computer resources. However, they have the advantage of being highly general.The TLM method has acquired a reputation for being a powerful and effective tool by numerous teams and still benefits today from significant theoretical developments. In particular, in recent years, its ability to simulate various situations with excellent precision, including complex materials, has been

Numerical and adaptive grid methods for ideal magnetohydrodynamics

Science.gov (United States)

Loring, Burlen

2008-02-01

In this thesis numerical finite difference methods for ideal magnetohydrodynamics(MHD) are investigated. A review of the relevant physics, essential for interpreting the results of numerical solutions and constructing validation cases, is presented. This review includes a discusion of the propagation of small amplitude waves in the MHD system as well as a thorough discussion of MHD shocks, contacts and rarefactions and how they can be piece together to obtain a solutions to the MHD Riemann problem. Numerical issues relevant to the MHD system such as: the loss of nonlinear numerical stability in the presence of discontinuous solutions, the introduction of spurious forces due to the growth of the divergence of the magnetic flux density, the loss of pressure positivity, and the effects of non-conservative numerical methods are discussed, along with the practical approaches which can be used to remedy or minimize the negative consequences of each. The use of block structured adaptive mesh refinement is investigated in the context of a divergence free MHD code. A new method for conserving magnetic flux across AMR grid interfaces is developed and a detailed discussion of our implementation of this method using the CHOMBO AMR framework is given. A preliminary validation of the new method for conserving magnetic flux density across AMR grid interfaces illustrates that the method works. Finally a number of code validation cases are examined spurring a discussion of the strengths and weaknesses of the numerics employed.

Numerical simulation of overflow at vertical weirs using a hybrid level set/VOF method

Science.gov (United States)

Lv, Xin; Zou, Qingping; Reeve, Dominic

2011-10-01

This paper presents the applications of a newly developed free surface flow model to the practical, while challenging overflow problems for weirs. Since the model takes advantage of the strengths of both the level set and volume of fluid methods and solves the Navier-Stokes equations on an unstructured mesh, it is capable of resolving the time evolution of very complex vortical motions, air entrainment and pressure variations due to violent deformations following overflow of the weir crest. In the present study, two different types of vertical weir, namely broad-crested and sharp-crested, are considered for validation purposes. The calculated overflow parameters such as pressure head distributions, velocity distributions, and water surface profiles are compared against experimental data as well as numerical results available in literature. A very good quantitative agreement has been obtained. The numerical model, thus, offers a good alternative to traditional experimental methods in the study of weir problems.

A numerical method for resonance integral calculations

International Nuclear Information System (INIS)

Tanbay, Tayfun; Ozgener, Bilge

2013-01-01

A numerical method has been proposed for resonance integral calculations and a cubic fit based on least squares approximation to compute the optimum Bell factor is given. The numerical method is based on the discretization of the neutron slowing down equation. The scattering integral is approximated by taking into account the location of the upper limit in energy domain. The accuracy of the method has been tested by performing computations of resonance integrals for uranium dioxide isolated rods and comparing the results with empirical values. (orig.)

Design of heat exchangers by numerical methods

International Nuclear Information System (INIS)

Konuk, A.A.

1981-01-01

Differential equations describing the heat tranfer in shell - and tube heat exchangers are derived and solved numerically. The method of ΔT sub(lm) is compared with the proposed method in cases where the specific heat at constant pressure, Cp and the overall heat transfer coefficient, U, vary with temperature. The error of the method of ΔT sub (lm) for the computation of the exchanger lenght is less than + 10%. However, the numerical method, being more accurate and at the same time easy to use and economical, is recommended for the design of shell-and-tube heat exchangers. (Author) [pt

Numerical methods for semiconductor heterostructures with band nonparabolicity

International Nuclear Information System (INIS)

Wang Weichung; Hwang Tsungmin; Lin Wenwei; Liu Jinnliang

2003-01-01

This article presents numerical methods for computing bound state energies and associated wave functions of three-dimensional semiconductor heterostructures with special interest in the numerical treatment of the effect of band nonparabolicity. A nonuniform finite difference method is presented to approximate a model of a cylindrical-shaped semiconductor quantum dot embedded in another semiconductor matrix. A matrix reduction method is then proposed to dramatically reduce huge eigenvalue systems to relatively very small subsystems. Moreover, the nonparabolic band structure results in a cubic type of nonlinear eigenvalue problems for which a cubic Jacobi-Davidson method with an explicit nonequivalence deflation method are proposed to compute all the desired eigenpairs. Numerical results are given to illustrate the spectrum of energy levels and the corresponding wave functions in rather detail

3D modeling of forces between magnet and HTS in a levitation system using new approach of the control volume method based on an unstructured grid

Energy Technology Data Exchange (ETDEWEB)

Alloui, L., E-mail: lotfi.alloui@lgep.supelec.fr [Laboratoire de Genie Electrique de Paris - LGEP, CNRS UMR 8507, Supelec, Universite Pierre et Marie Curie-Paris 6, Universite Paris Sud-Paris 11, Plateau de Moulon, 11 rue Joliot Curie, 91192 Gif-Sur-Yvette Cedex (France); Laboratoire de modelisation des systemes energetiques (LMSE), Universite de Biskra, 07000 Biskra (Algeria); Bouillault, F., E-mail: bouillault@lgep.supelec.fr [Laboratoire de Genie Electrique de Paris - LGEP, CNRS UMR 8507, Supelec, Universite Pierre et Marie Curie-Paris 6, Universite Paris Sud-Paris 11, Plateau de Moulon, 11 rue Joliot Curie, 91192 Gif-Sur-Yvette Cedex (France); Bernard, L., E-mail: laurent.bernardl@lgep.supelc.fr [Laboratoire de Genie Electrique de Paris - LGEP, CNRS UMR 8507, Supelec, Universite Pierre et Marie Curie-Paris 6, Universite Paris Sud-Paris 11, Plateau de Moulon, 11 rue Joliot Curie, 91192 Gif-Sur-Yvette Cedex (France); Leveque, J., E-mail: jean.leveque@green.uhp-nancy.fr [Groupe de recherche en electronique et electrotechnique de Nancy, Universite Henry Poincare, BP 239, 54506 Vandoeuvre les Nancy (France)

2012-05-15

In this paper we present new 3D numerical model to calculate the vertical and the guidance forces in high temperature superconductors taking into account the influence of the flux creep phenomena. In the suggested numerical model, we adopt a new approach of the control volume method. This approach is based on the use of an unstructured grid which can be used to model more complex geometries. A comparison of the control volume method results with experiments verifies the validity of this approach and the proposed numerical model. Based on this model, the levitation force's relaxation at different temperatures was also studied.

Numerical analysis of partially molten splat during thermal spray process using the finite element method

Science.gov (United States)

Zirari, M.; Abdellah El-Hadj, A.; Bacha, N.

2010-03-01

A finite element method is used to simulate the deposition of the thermal spray coating process. A set of governing equations is solving by a volume of fluid method. For the solidification phenomenon, we use the specific heat method (SHM). We begin by comparing the present model with experimental and numerical model available in the literature. In this study, completely molten or semi-molten aluminum particle impacts a H13 tool steel substrate is considered. Next we investigate the effect of inclination of impact of a partially molten particle on flat substrate. It was found that the melting state of the particle has great effects on the morphologies of the splat.

Numerical implementation of the loop-tree duality method

Energy Technology Data Exchange (ETDEWEB)

Buchta, Sebastian; Rodrigo, German [Universitat de Valencia-Consejo Superior de Investigaciones Cientificas, Parc Cientific, Instituto de Fisica Corpuscular, Valencia (Spain); Chachamis, Grigorios [Universidad Autonoma de Madrid, Instituto de Fisica Teorica UAM/CSIC, Madrid (Spain); Draggiotis, Petros [Institute of Nuclear and Particle Physics, NCSR ' ' Demokritos' ' , Agia Paraskevi (Greece)

2017-05-15

We present a first numerical implementation of the loop-tree duality (LTD) method for the direct numerical computation of multi-leg one-loop Feynman integrals. We discuss in detail the singular structure of the dual integrands and define a suitable contour deformation in the loop three-momentum space to carry out the numerical integration. Then we apply the LTD method to the computation of ultraviolet and infrared finite integrals, and we present explicit results for scalar and tensor integrals with up to eight external legs (octagons). The LTD method features an excellent performance independently of the number of external legs. (orig.)

Numerical methods for metamaterial design

CERN Document Server

2013-01-01

This book describes a relatively new approach for the design of electromagnetic metamaterials.  Numerical optimization routines are combined with electromagnetic simulations to tailor the broadband optical properties of a metamaterial to have predetermined responses at predetermined wavelengths. After a review of both the major efforts within the field of metamaterials and the field of mathematical optimization, chapters covering both gradient-based and derivative-free design methods are considered.  Selected topics including surrogate-base optimization, adaptive mesh search, and genetic algorithms are shown to be effective, gradient-free optimization strategies.  Additionally, new techniques for representing dielectric distributions in two dimensions, including level sets, are demonstrated as effective methods for gradient-based optimization.  Each chapter begins with a rigorous review of the optimization strategy used, and is followed by numerous examples that combine the strategy with either electromag...

Lagrangian numerical methods for ocean biogeochemical simulations

Science.gov (United States)

Paparella, Francesco; Popolizio, Marina

2018-05-01

We propose two closely-related Lagrangian numerical methods for the simulation of physical processes involving advection, reaction and diffusion. The methods are intended to be used in settings where the flow is nearly incompressible and the Péclet numbers are so high that resolving all the scales of motion is unfeasible. This is commonplace in ocean flows. Our methods consist in augmenting the method of characteristics, which is suitable for advection-reaction problems, with couplings among nearby particles, producing fluxes that mimic diffusion, or unresolved small-scale transport. The methods conserve mass, obey the maximum principle, and allow to tune the strength of the diffusive terms down to zero, while avoiding unwanted numerical dissipation effects.

Mathematical modelling and numerical simulation of casting processes

DEFF Research Database (Denmark)

Hattel, Jesper Henri

1998-01-01

The control volume method applied to numerical modelling of castning. Analytical solutions based on the error function.Riemann-temperature. Modelling of release of latent heat with the enthalpy method....

Numerical investigation on the influence of surface tension and viscous force on the bubble dynamics with a CLSVOF method

Energy Technology Data Exchange (ETDEWEB)

Wang, Zhiying; Li, Yikai; Huang, Biao; Gao, Deming [Beijing Institute of Technology, Beijing (China)

2016-06-15

We numerically investigated the rising of bubbles in a quiescent liquid layer. The numerical simulation is performed by solving the incompressible, multiphase Navier-Stokes equations via computational code in axisymmetric coordinates using a Coupled level-set and volume-of-fluid (CLSVOF) method. The numerical results show that the CLSVOF method with a novel algebraic relation between F and f for axisymmetric two-phase flows not only can predict the bubble surface accurately, but also overcome the deficiency in preserving volume conservation. The effects of the Reynolds number Re and the Bond number Bo on the bubble deformation and its motion are investigated. The results show that with the increasing of Re (10 < Re < 150), the bubble shape transfers from oblate ellipsoidal cap to toroidal when Bo = 116. With the increasing of Bo (10 < Bo < 700), the bubble shape transfers from oblate ellipsoidal to toroidal when Re = 30. Although the toroidal bubble shapes are reached in these two cases, the transition modes are different. For the case Bo = 116, the bubble front is pierced by an upward jet from the rear of the bubble. While for the case Re = 30, the rear of the bubble is pierced by a downward jet from the front part.

Molecular dynamics with deterministic and stochastic numerical methods

CERN Document Server

Leimkuhler, Ben

2015-01-01

This book describes the mathematical underpinnings of algorithms used for molecular dynamics simulation, including both deterministic and stochastic numerical methods. Molecular dynamics is one of the most versatile and powerful methods of modern computational science and engineering and is used widely in chemistry, physics, materials science and biology. Understanding the foundations of numerical methods means knowing how to select the best one for a given problem (from the wide range of techniques on offer) and how to create new, efficient methods to address particular challenges as they arise in complex applications.  Aimed at a broad audience, this book presents the basic theory of Hamiltonian mechanics and stochastic differential equations, as well as topics including symplectic numerical methods, the handling of constraints and rigid bodies, the efficient treatment of Langevin dynamics, thermostats to control the molecular ensemble, multiple time-stepping, and the dissipative particle dynamics method...

Tensor viscosity method for convection in numerical fluid dynamics

International Nuclear Information System (INIS)

Dukowicz, J.K.; Ramshaw, J.D.

1979-01-01

A new method, called the tensor viscosity method, is described for differencing the convective terms in multidimensional numerical fluid dynamics. The method is the proper generalization to two or three dimensions of interpolated donor cell differencing in one dimension, and is designed to achieve numerical stability with minimal numerical damping. It is a single-step method that is distinguished by simplicity and case of implementation, even in the case of an arbitrary non-rectangular mesh. It should therefore be useful in finite-element as well as finite-difference formulations

Numerical simulation methods for phase-transitional flow

NARCIS (Netherlands)

Pecenko, A.

2010-01-01

The object of the present dissertation is a numerical study of multiphase flow of one fluid component. In particular, the research described in this thesis focuses on the development of numerical methods that are based on a diffuse-interface model (DIM). With this approach, the modeling problem

Assessing numerical methods used in nuclear aerosol transport models

International Nuclear Information System (INIS)

McDonald, B.H.

1987-01-01

Several computer codes are in use for predicting the behaviour of nuclear aerosols released into containment during postulated accidents in water-cooled reactors. Each of these codes uses numerical methods to discretize and integrate the equations that govern the aerosol transport process. Computers perform only algebraic operations and generate only numbers. It is in the numerical methods that sense can be made of these numbers and where they can be related to the actual solution of the equations. In this report, the numerical methods most commonly used in the aerosol transport codes are examined as special cases of a general solution procedure, the Method of Weighted Residuals. It would appear that the numerical methods used in the codes are all capable of producing reasonable answers to the mathematical problem when used with skill and care. 27 refs

Nonlinear ordinary differential equations analytical approximation and numerical methods

CERN Document Server

Hermann, Martin

2016-01-01

The book discusses the solutions to nonlinear ordinary differential equations (ODEs) using analytical and numerical approximation methods. Recently, analytical approximation methods have been largely used in solving linear and nonlinear lower-order ODEs. It also discusses using these methods to solve some strong nonlinear ODEs. There are two chapters devoted to solving nonlinear ODEs using numerical methods, as in practice high-dimensional systems of nonlinear ODEs that cannot be solved by analytical approximate methods are common. Moreover, it studies analytical and numerical techniques for the treatment of parameter-depending ODEs. The book explains various methods for solving nonlinear-oscillator and structural-system problems, including the energy balance method, harmonic balance method, amplitude frequency formulation, variational iteration method, homotopy perturbation method, iteration perturbation method, homotopy analysis method, simple and multiple shooting method, and the nonlinear stabilized march...

Numerical methods in multibody dynamics

CERN Document Server

Eich-Soellner, Edda

1998-01-01

Today computers play an important role in the development of complex mechanical systems, such as cars, railway vehicles or machines. Efficient simulation of these systems is only possible when based on methods that explore the strong link between numerics and computational mechanics. This book gives insight into modern techniques of numerical mathematics in the light of an interesting field of applications: multibody dynamics. The important interaction between modeling and solution techniques is demonstrated by using a simplified multibody model of a truck. Different versions of this mechanical model illustrate all key concepts in static and dynamic analysis as well as in parameter identification. The book focuses in particular on constrained mechanical systems. Their formulation in terms of differential-algebraic equations is the backbone of nearly all chapters. The book is written for students and teachers in numerical analysis and mechanical engineering as well as for engineers in industrial research labor...

A finite volume method for density driven flows in porous media

Directory of Open Access Journals (Sweden)

Hilhorst Danielle

2013-01-01

Full Text Available In this paper, we apply a semi-implicit finite volume method for the numerical simulation of density driven flows in porous media; this amounts to solving a nonlinear convection-diffusion parabolic equation for the concentration coupled with an elliptic equation for the pressure. We compute the solutions for two specific problems: a problem involving a rotating interface between salt and fresh water and the classical but difficult Henry’s problem. All solutions are compared to results obtained by running FEflow, a commercial software package for the simulation of groundwater flow, mass and heat transfer in porous media.

Numerical Methods for a Class of Differential Algebraic Equations

Directory of Open Access Journals (Sweden)

Lei Ren

2017-01-01

Full Text Available This paper is devoted to the study of some efficient numerical methods for the differential algebraic equations (DAEs. At first, we propose a finite algorithm to compute the Drazin inverse of the time varying DAEs. Numerical experiments are presented by Drazin inverse and Radau IIA method, which illustrate that the precision of the Drazin inverse method is higher than the Radau IIA method. Then, Drazin inverse, Radau IIA, and Padé approximation are applied to the constant coefficient DAEs, respectively. Numerical results demonstrate that the Padé approximation is powerful for solving constant coefficient DAEs.

Numerical computer methods part D

CERN Document Server

Johnson, Michael L

2004-01-01

The aim of this volume is to brief researchers of the importance of data analysis in enzymology, and of the modern methods that have developed concomitantly with computer hardware. It is also to validate researchers' computer programs with real and synthetic data to ascertain that the results produced are what they expected. Selected Contents: Prediction of protein structure; modeling and studying proteins with molecular dynamics; statistical error in isothermal titration calorimetry; analysis of circular dichroism data; model comparison methods.

Finite Volume Element Predictor-corrector Method for a Class of Nonlinear Parabolic Systems%一类非线性抛物型方程组的有限体积元预估-校正方法

Institute of Scientific and Technical Information of China (English)

高夫征

2005-01-01

A finite volume element predictor-correetor method for a class of nonlinear parabolic system of equations is presented and analyzed. Suboptimal L2 error estimate for the finite volume element predictor-corrector method is derived. A numerical experiment shows that the numerical results are consistent with theoretical analysis.

Mathematical analysis and numerical methods for science and technology

CERN Document Server

Dautray, Robert

These 6 volumes - the result of a 10 year collaboration between the authors, two of France's leading scientists and both distinguished international figures - compile the mathematical knowledge required by researchers in mechanics, physics, engineering, chemistry and other branches of application of mathematics for the theoretical and numerical resolution of physical models on computers. Since the publication in 1924 of the "Methoden der mathematischen Physik" by Courant and Hilbert, there has been no other comprehensive and up-to-date publication presenting the mathematical tools needed in applications of mathematics in directly implementable form. The advent of large computers has in the meantime revolutionised methods of computation and made this gap in the literature intolerable: the objective of the present work is to fill just this gap. Many phenomena in physical mathematics may be modeled by a system of partial differential equations in distributed systems: a model here means a set of equations, which ...

Numerical methods and analysis of multiscale problems

CERN Document Server

Madureira, Alexandre L

2017-01-01

This book is about numerical modeling of multiscale problems, and introduces several asymptotic analysis and numerical techniques which are necessary for a proper approximation of equations that depend on different physical scales. Aimed at advanced undergraduate and graduate students in mathematics, engineering and physics – or researchers seeking a no-nonsense approach –, it discusses examples in their simplest possible settings, removing mathematical hurdles that might hinder a clear understanding of the methods. The problems considered are given by singular perturbed reaction advection diffusion equations in one and two-dimensional domains, partial differential equations in domains with rough boundaries, and equations with oscillatory coefficients. This work shows how asymptotic analysis can be used to develop and analyze models and numerical methods that are robust and work well for a wide range of parameters.

Spectral Methods in Numerical Plasma Simulation

DEFF Research Database (Denmark)

Coutsias, E.A.; Hansen, F.R.; Huld, T.

1989-01-01

An introduction is given to the use of spectral methods in numerical plasma simulation. As examples of the use of spectral methods, solutions to the two-dimensional Euler equations in both a simple, doubly periodic region, and on an annulus will be shown. In the first case, the solution is expanded...

Evaluation of the Reference Numerical Parameters of the Monthly Method in ISO 13790 Considering S/V Ratio

Directory of Open Access Journals (Sweden)

Hee-Jeong Kwak

2015-01-01

Full Text Available Many studies have investigated the accuracy of the numerical parameters in the application of the quasi steady-state calculation method. The aim of this study is to derive the reference numerical parameters of the ISO 13790 monthly method by reflecting the surface-to-volume (S/V ratio and the characteristics of the structures. The calculation process was established, and the parameters necessary to derive the reference numerical parameters were calculated based on the input data prepared for the established calculation processes. The reference numerical parameters were then derived through regression analyses of the calculated parameters and the time constant. The parameters obtained from an apartment building and the parameters of the international standard were both applied to the Passive House Planning Package (PHPP and EnergyPlus programs, and the results were analyzed in order to evaluate the validity of the results. The analysis revealed that the calculation results based on the parameters derived from this study yielded lower error rates than those based on the default parameters in ISO 13790. However, the differences were shown to be negligible in the case of high heat capacity.

Colorado Conference on iterative methods. Volume 1

Energy Technology Data Exchange (ETDEWEB)

NONE

1994-12-31

The conference provided a forum on many aspects of iterative methods. Volume I topics were:Session: domain decomposition, nonlinear problems, integral equations and inverse problems, eigenvalue problems, iterative software kernels. Volume II presents nonsymmetric solvers, parallel computation, theory of iterative methods, software and programming environment, ODE solvers, multigrid and multilevel methods, applications, robust iterative methods, preconditioners, Toeplitz and circulation solvers, and saddle point problems. Individual papers are indexed separately on the EDB.

Numerical solutions of the N-body problem

International Nuclear Information System (INIS)

Marciniak, A.

1985-01-01

Devoted to the study of numerical methods for solving the general N-body problem and related problems, this volume starts with an overview of the conventional numerical methods for solving the initial value problem. The major part of the book contains original work and features a presentation of special numerical methods conserving the constants of motion in the general N-body problem and methods conserving the Jacobi constant in the problem of motion of N bodies in a rotating frame, as well as an analysis of the applications of both (conventional and special) kinds of methods for solving these problems. For all the methods considered, the author presents algorithms which are easily programmable in any computer language. Moreover, the author compares various methods and presents adequate numerical results. The appendix contains PL/I procedures for all the special methods conserving the constants of motion. 91 refs.; 35 figs.; 41 tabs

Numerical methods and optimization a consumer guide

CERN Document Server

Walter, Éric

2014-01-01

Initial training in pure and applied sciences tends to present problem-solving as the process of elaborating explicit closed-form solutions from basic principles, and then using these solutions in numerical applications. This approach is only applicable to very limited classes of problems that are simple enough for such closed-form solutions to exist. Unfortunately, most real-life problems are too complex to be amenable to this type of treatment. Numerical Methods and Optimization – A Consumer Guide presents methods for dealing with them. Shifting the paradigm from formal calculus to numerical computation, the text makes it possible for the reader to ·         discover how to escape the dictatorship of those particular cases that are simple enough to receive a closed-form solution, and thus gain the ability to solve complex, real-life problems; ·         understand the principles behind recognized algorithms used in state-of-the-art numerical software; ·         learn the advantag...

Operator theory and numerical methods

CERN Document Server

Fujita, H; Suzuki, T

2001-01-01

In accordance with the developments in computation, theoretical studies on numerical schemes are now fruitful and highly needed. In 1991 an article on the finite element method applied to evolutionary problems was published. Following the method, basically this book studies various schemes from operator theoretical points of view. Many parts are devoted to the finite element method, but other schemes and problems (charge simulation method, domain decomposition method, nonlinear problems, and so forth) are also discussed, motivated by the observation that practically useful schemes have fine mathematical structures and the converses are also true. This book has the following chapters: 1. Boundary Value Problems and FEM. 2. Semigroup Theory and FEM. 3. Evolution Equations and FEM. 4. Other Methods in Time Discretization. 5. Other Methods in Space Discretization. 6. Nonlinear Problems. 7. Domain Decomposition Method.

Numerical models for differential problems

CERN Document Server

Quarteroni, Alfio

2017-01-01

In this text, we introduce the basic concepts for the numerical modelling of partial differential equations. We consider the classical elliptic, parabolic and hyperbolic linear equations, but also the diffusion, transport, and Navier-Stokes equations, as well as equations representing conservation laws, saddle-point problems and optimal control problems. Furthermore, we provide numerous physical examples which underline such equations. We then analyze numerical solution methods based on finite elements, finite differences, finite volumes, spectral methods and domain decomposition methods, and reduced basis methods. In particular, we discuss the algorithmic and computer implementation aspects and provide a number of easy-to-use programs. The text does not require any previous advanced mathematical knowledge of partial differential equations: the absolutely essential concepts are reported in a preliminary chapter. It is therefore suitable for students of bachelor and master courses in scientific disciplines, an...

Numerical methods for scientists and engineers

CERN Document Server

Antia, H M

2012-01-01

This book presents an exhaustive and in-depth exposition of the various numerical methods used in scientific and engineering computations. It emphasises the practical aspects of numerical computation and discusses various techniques in sufficient detail to enable their implementation in solving a wide range of problems. The main addition in the third edition is a new Chapter on Statistical Inferences. There is also some addition and editing in the next chapter on Approximations. With this addition 12 new programs have also been added.

Explicit appropriate basis function method for numerical solution of stiff systems

International Nuclear Information System (INIS)

Chen, Wenzhen; Xiao, Hongguang; Li, Haofeng; Chen, Ling

2015-01-01

Highlights: • An explicit numerical method called the appropriate basis function method is presented. • The method differs from the power series method for obtaining approximate numerical solutions. • Two cases show the method is fit for linear and nonlinear stiff systems. • The method is very simple and effective for most of differential equation systems. - Abstract: In this paper, an explicit numerical method, called the appropriate basis function method, is presented. The explicit appropriate basis function method differs from the power series method because it employs an appropriate basis function such as the exponential function, or periodic function, other than a polynomial, to obtain approximate numerical solutions. The method is successful and effective for the numerical solution of the first order ordinary differential equations. Two examples are presented to show the ability of the method for dealing with linear and nonlinear systems of differential equations

Numerical methods for stochastic partial differential equations with white noise

CERN Document Server

Zhang, Zhongqiang

2017-01-01

This book covers numerical methods for stochastic partial differential equations with white noise using the framework of Wong-Zakai approximation. The book begins with some motivational and background material in the introductory chapters and is divided into three parts. Part I covers numerical stochastic ordinary differential equations. Here the authors start with numerical methods for SDEs with delay using the Wong-Zakai approximation and finite difference in time. Part II covers temporal white noise. Here the authors consider SPDEs as PDEs driven by white noise, where discretization of white noise (Brownian motion) leads to PDEs with smooth noise, which can then be treated by numerical methods for PDEs. In this part, recursive algorithms based on Wiener chaos expansion and stochastic collocation methods are presented for linear stochastic advection-diffusion-reaction equations. In addition, stochastic Euler equations are exploited as an application of stochastic collocation methods, where a numerical compa...

Numerical method of singular problems on singular integrals

International Nuclear Information System (INIS)

Zhao Huaiguo; Mou Zongze

1992-02-01

As first part on the numerical research of singular problems, a numerical method is proposed for singular integrals. It is shown that the procedure is quite powerful for solving physics calculation with singularity such as the plasma dispersion function. Useful quadrature formulas for some class of the singular integrals are derived. In general, integrals with more complex singularities can be dealt by this method easily

Analysis of one-dimensional nonequilibrium two-phase flow using control volume method

International Nuclear Information System (INIS)

Minato, Akihiko; Naitoh, Masanori

1987-01-01

A one-dimensional numerical analysis model was developed for prediction of rapid flow transient behavior involving boiling. This model was based on six conservation equations of time averaged parameters of gas and liquid behavior. These equations were solved by using a control volume method with an explicit time integration. This model did not use staggered mesh scheme, which had been commonly used in two-phase flow analysis. Because void fraction and velocity of each phase were defined at the same location in the present model, effects of void fraction on phase velocity calculation were treated directly without interpolation. Though non-staggered mesh scheme was liable to cause numerical instability with zigzag pressure field, stability was achieved by employing the Godunov method. In order to verify the present analytical model, Edwards' pipe blow down and Zaloudek's initially subcooled critical two-phase flow experiments were analyzed. Stable solutions were obtained for rarefaction wave propagation with boiling and transient two-phase flow behavior in a broken pipe by using this model. (author)

Numerical methods for characterization of synchrotron radiation based on the Wigner function method

Directory of Open Access Journals (Sweden)

Takashi Tanaka

2014-06-01

Full Text Available Numerical characterization of synchrotron radiation based on the Wigner function method is explored in order to accurately evaluate the light source performance. A number of numerical methods to compute the Wigner functions for typical synchrotron radiation sources such as bending magnets, undulators and wigglers, are presented, which significantly improve the computation efficiency and reduce the total computation time. As a practical example of the numerical characterization, optimization of betatron functions to maximize the brilliance of undulator radiation is discussed.

Coincidental match of numerical simulation and physics

Science.gov (United States)

Pierre, B.; Gudmundsson, J. S.

2010-08-01

Consequences of rapid pressure transients in pipelines range from increased fatigue to leakages and to complete ruptures of pipeline. Therefore, accurate predictions of rapid pressure transients in pipelines using numerical simulations are critical. State of the art modelling of pressure transient in general, and water hammer in particular include unsteady friction in addition to the steady frictional pressure drop, and numerical simulations rely on the method of characteristics. Comparison of rapid pressure transient calculations by the method of characteristics and a selected high resolution finite volume method highlights issues related to modelling of pressure waves and illustrates that matches between numerical simulations and physics are purely coincidental.

Numerical methods for engine-airframe integration

International Nuclear Information System (INIS)

Murthy, S.N.B.; Paynter, G.C.

1986-01-01

Various papers on numerical methods for engine-airframe integration are presented. The individual topics considered include: scientific computing environment for the 1980s, overview of prediction of complex turbulent flows, numerical solutions of the compressible Navier-Stokes equations, elements of computational engine/airframe integrations, computational requirements for efficient engine installation, application of CAE and CFD techniques to complete tactical missile design, CFD applications to engine/airframe integration, and application of a second-generation low-order panel methods to powerplant installation studies. Also addressed are: three-dimensional flow analysis of turboprop inlet and nacelle configurations, application of computational methods to the design of large turbofan engine nacelles, comparison of full potential and Euler solution algorithms for aeropropulsive flow field computations, subsonic/transonic, supersonic nozzle flows and nozzle integration, subsonic/transonic prediction capabilities for nozzle/afterbody configurations, three-dimensional viscous design methodology of supersonic inlet systems for advanced technology aircraft, and a user's technology assessment

NUMERICAL AND ANALYTIC METHODS OF ESTIMATION BRIDGES’ CONSTRUCTIONS

Directory of Open Access Journals (Sweden)

Y. Y. Luchko

2010-03-01

Full Text Available In this article the numerical and analytical methods of calculation of the stressed-and-strained state of bridge constructions are considered. The task on increasing of reliability and accuracy of the numerical method and its solution by means of calculations in two bases are formulated. The analytical solution of the differential equation of deformation of a ferro-concrete plate under the action of local loads is also obtained.

A critical analysis of some popular methods for the discretisation of the gradient operator in finite volume methods

Science.gov (United States)

Syrakos, Alexandros; Varchanis, Stylianos; Dimakopoulos, Yannis; Goulas, Apostolos; Tsamopoulos, John

2017-12-01

Finite volume methods (FVMs) constitute a popular class of methods for the numerical simulation of fluid flows. Among the various components of these methods, the discretisation of the gradient operator has received less attention despite its fundamental importance with regards to the accuracy of the FVM. The most popular gradient schemes are the divergence theorem (DT) (or Green-Gauss) scheme and the least-squares (LS) scheme. Both are widely believed to be second-order accurate, but the present study shows that in fact the common variant of the DT gradient is second-order accurate only on structured meshes whereas it is zeroth-order accurate on general unstructured meshes, and the LS gradient is second-order and first-order accurate, respectively. This is explained through a theoretical analysis and is confirmed by numerical tests. The schemes are then used within a FVM to solve a simple diffusion equation on unstructured grids generated by several methods; the results reveal that the zeroth-order accuracy of the DT gradient is inherited by the FVM as a whole, and the discretisation error does not decrease with grid refinement. On the other hand, use of the LS gradient leads to second-order accurate results, as does the use of alternative, consistent, DT gradient schemes, including a new iterative scheme that makes the common DT gradient consistent at almost no extra cost. The numerical tests are performed using both an in-house code and the popular public domain partial differential equation solver OpenFOAM.

Stability, accuracy and numerical diffusion analysis of nodal expansion method for steady convection diffusion equation

International Nuclear Information System (INIS)

Zhou, Xiafeng; Guo, Jiong; Li, Fu

2015-01-01

Highlights: • NEMs are innovatively applied to solve convection diffusion equation. • Stability, accuracy and numerical diffusion for NEM are analyzed for the first time. • Stability and numerical diffusion depend on the NEM expansion order and its parity. • NEMs have higher accuracy than both second order upwind and QUICK scheme. • NEMs with different expansion orders are integrated into a unified discrete form. - Abstract: The traditional finite difference method or finite volume method (FDM or FVM) is used for HTGR thermal-hydraulic calculation at present. However, both FDM and FVM require the fine mesh sizes to achieve the desired precision and thus result in a limited efficiency. Therefore, a more efficient and accurate numerical method needs to be developed. Nodal expansion method (NEM) can achieve high accuracy even on the coarse meshes in the reactor physics analysis so that the number of spatial meshes and computational cost can be largely decreased. Because of higher efficiency and accuracy, NEM can be innovatively applied to thermal-hydraulic calculation. In the paper, NEMs with different orders of basis functions are successfully developed and applied to multi-dimensional steady convection diffusion equation. Numerical results show that NEMs with three or higher order basis functions can track the reference solutions very well and are superior to second order upwind scheme and QUICK scheme. However, the false diffusion and unphysical oscillation behavior are discovered for NEMs. To explain the reasons for the above-mentioned behaviors, the stability, accuracy and numerical diffusion properties of NEM are analyzed by the Fourier analysis, and by comparing with exact solutions of difference and differential equation. The theoretical analysis results show that the accuracy of NEM increases with the expansion order. However, the stability and numerical diffusion properties depend not only on the order of basis functions but also on the parity of

Stability, accuracy and numerical diffusion analysis of nodal expansion method for steady convection diffusion equation

Energy Technology Data Exchange (ETDEWEB)

Zhou, Xiafeng, E-mail: zhou-xf11@mails.tsinghua.edu.cn; Guo, Jiong, E-mail: guojiong12@tsinghua.edu.cn; Li, Fu, E-mail: lifu@tsinghua.edu.cn

2015-12-15

Highlights: • NEMs are innovatively applied to solve convection diffusion equation. • Stability, accuracy and numerical diffusion for NEM are analyzed for the first time. • Stability and numerical diffusion depend on the NEM expansion order and its parity. • NEMs have higher accuracy than both second order upwind and QUICK scheme. • NEMs with different expansion orders are integrated into a unified discrete form. - Abstract: The traditional finite difference method or finite volume method (FDM or FVM) is used for HTGR thermal-hydraulic calculation at present. However, both FDM and FVM require the fine mesh sizes to achieve the desired precision and thus result in a limited efficiency. Therefore, a more efficient and accurate numerical method needs to be developed. Nodal expansion method (NEM) can achieve high accuracy even on the coarse meshes in the reactor physics analysis so that the number of spatial meshes and computational cost can be largely decreased. Because of higher efficiency and accuracy, NEM can be innovatively applied to thermal-hydraulic calculation. In the paper, NEMs with different orders of basis functions are successfully developed and applied to multi-dimensional steady convection diffusion equation. Numerical results show that NEMs with three or higher order basis functions can track the reference solutions very well and are superior to second order upwind scheme and QUICK scheme. However, the false diffusion and unphysical oscillation behavior are discovered for NEMs. To explain the reasons for the above-mentioned behaviors, the stability, accuracy and numerical diffusion properties of NEM are analyzed by the Fourier analysis, and by comparing with exact solutions of difference and differential equation. The theoretical analysis results show that the accuracy of NEM increases with the expansion order. However, the stability and numerical diffusion properties depend not only on the order of basis functions but also on the parity of

Theoretical and numerical method in aeroacoustics

Directory of Open Access Journals (Sweden)

Nicuşor ALEXANDRESCU

2010-06-01

Full Text Available The paper deals with the mathematical and numerical modeling of the aerodynamic noisegenerated by the fluid flow interaction with the solid structure of a rotor blade.Our analysis use Lighthill’s acoustic analogy. Lighthill idea was to express the fundamental equationsof motion into a wave equation for acoustic fluctuation with a source term on the right-hand side. Theobtained wave equation is solved numerically by the spatial discretization. The method is applied inthe case of monopole source placed in different points of blade surfaces to find this effect of noisepropagation.

Numerical methods in finance and economics a MATLAB-based introduction

CERN Document Server

Brandimarte, Paolo

2006-01-01

A state-of-the-art introduction to the powerful mathematical and statistical tools used in the field of financeThe use of mathematical models and numerical techniques is a practice employed by a growing number of applied mathematicians working on applications in finance. Reflecting this development, Numerical Methods in Finance and Economics: A MATLAB?-Based Introduction, Second Edition bridges the gap between financial theory and computational practice while showing readers how to utilize MATLAB?--the powerful numerical computing environment--for financial applications.The author provides an essential foundation in finance and numerical analysis in addition to background material for students from both engineering and economics perspectives. A wide range of topics is covered, including standard numerical analysis methods, Monte Carlo methods to simulate systems affected by significant uncertainty, and optimization methods to find an optimal set of decisions.Among this book''s most outstanding features is the...

Stochastic numerical methods an introduction for students and scientists

CERN Document Server

Toral, Raul

2014-01-01

Stochastic Numerical Methods introduces at Master level the numerical methods that use probability or stochastic concepts to analyze random processes. The book aims at being rather general and is addressed at students of natural sciences (Physics, Chemistry, Mathematics, Biology, etc.) and Engineering, but also social sciences (Economy, Sociology, etc.) where some of the techniques have been used recently to numerically simulate different agent-based models. Examples included in the book range from phase-transitions and critical phenomena, including details of data analysis (extraction of critical exponents, finite-size effects, etc.), to population dynamics, interfacial growth, chemical reactions, etc. Program listings are integrated in the discussion of numerical algorithms to facilitate their understanding. From the contents: Review of Probability ConceptsMonte Carlo IntegrationGeneration of Uniform and Non-uniformRandom Numbers: Non-correlated ValuesDynamical MethodsApplications to Statistical MechanicsIn...

Numerical methods design, analysis, and computer implementation of algorithms

CERN Document Server

Greenbaum, Anne

2012-01-01

Numerical Methods provides a clear and concise exploration of standard numerical analysis topics, as well as nontraditional ones, including mathematical modeling, Monte Carlo methods, Markov chains, and fractals. Filled with appealing examples that will motivate students, the textbook considers modern application areas, such as information retrieval and animation, and classical topics from physics and engineering. Exercises use MATLAB and promote understanding of computational results. The book gives instructors the flexibility to emphasize different aspects--design, analysis, or computer implementation--of numerical algorithms, depending on the background and interests of students. Designed for upper-division undergraduates in mathematics or computer science classes, the textbook assumes that students have prior knowledge of linear algebra and calculus, although these topics are reviewed in the text. Short discussions of the history of numerical methods are interspersed throughout the chapters. The book a...

Three-dimensional reconstruction volume: a novel method for volume measurement in kidney cancer.

Science.gov (United States)

Durso, Timothy A; Carnell, Jonathan; Turk, Thomas T; Gupta, Gopal N

2014-06-01

The role of volumetric estimation is becoming increasingly important in the staging, management, and prognostication of benign and cancerous conditions of the kidney. We evaluated the use of three-dimensional reconstruction volume (3DV) in determining renal parenchymal volumes (RPV) and renal tumor volumes (RTV). We compared 3DV with the currently available methods of volume assessment and determined its interuser reliability. RPV and RTV were assessed in 28 patients who underwent robot-assisted laparoscopic partial nephrectomy for kidney cancer. Patients with a preoperative creatinine level of kidney pre- and postsurgery overestimated 3D reconstruction volumes by 15% to 102% and 12% to 101%, respectively. In addition, volumes obtained from 3DV displayed high interuser reliability regardless of experience. 3DV provides a highly reliable way of assessing kidney volumes. Given that 3DV takes into account visible anatomy, the differences observed using previously published methods can be attributed to the failure of geometry to accurately approximate kidney or tumor shape. 3DV provides a more accurate, reproducible, and clinically useful tool for urologists looking to improve patient care using analysis related to volume.

Development of numerical methods for reactive transport

International Nuclear Information System (INIS)

Bouillard, N.

2006-12-01

When a radioactive waste is stored in deep geological disposals, it is expected that the waste package will be damaged under water action (concrete leaching, iron corrosion). Then, to understand these damaging processes, chemical reactions and solutes transport are modelled. Numerical simulations of reactive transport can be done sequentially by the coupling of several codes. This is the case of the software platform ALLIANCES which is developed jointly with CEA, ANDRA and EDF. Stiff reactions like precipitation-dissolution are crucial for the radioactive waste storage applications, but standard sequential iterative approaches like Picard's fail in solving rapidly reactive transport simulations with such stiff reactions. In the first part of this work, we focus on a simplified precipitation and dissolution process: a system made up with one solid species and two aqueous species moving by diffusion is studied mathematically. It is assumed that a precipitation dissolution reaction occurs in between them, and it is modelled by a discontinuous kinetics law of unknown sign. By using monotonicity properties, the convergence of a finite volume scheme on admissible mesh is proved. Existence of a weak solution is obtained as a by-product of the convergence of the scheme. The second part is dedicated to coupling algorithms which improve Picard's method and can be easily used in an existing coupling code. By extending previous works, we propose a general and adaptable framework to solve nonlinear systems. Indeed by selecting special options, we can either recover well known methods, like nonlinear conjugate gradient methods, or design specific method. This algorithm has two main steps, a preconditioning one and an acceleration one. This algorithm is tested on several examples, some of them being rather academical and others being more realistic. We test it on the 'three species model'' example. Other reactive transport simulations use an external chemical code CHESS. For a

A consistent method for finite volume discretization of body forces on collocated grids applied to flow through an actuator disk

DEFF Research Database (Denmark)

Troldborg, Niels; Sørensen, Niels N.; Réthoré, Pierre-Elouan

2015-01-01

This paper describes a consistent algorithm for eliminating the numerical wiggles appearing when solving the finite volume discretized Navier-Stokes equations with discrete body forces in a collocated grid arrangement. The proposed method is a modification of the Rhie-Chow algorithm where the for...

A student's guide to numerical methods

CERN Document Server

Hutchinson, Ian H

2015-01-01

This concise, plain-language guide for senior undergraduates and graduate students aims to develop intuition, practical skills and an understanding of the framework of numerical methods for the physical sciences and engineering. It provides accessible self-contained explanations of mathematical principles, avoiding intimidating formal proofs. Worked examples and targeted exercises enable the student to master the realities of using numerical techniques for common needs such as solution of ordinary and partial differential equations, fitting experimental data, and simulation using particle and Monte Carlo methods. Topics are carefully selected and structured to build understanding, and illustrate key principles such as: accuracy, stability, order of convergence, iterative refinement, and computational effort estimation. Enrichment sections and in-depth footnotes form a springboard to more advanced material and provide additional background. Whether used for self-study, or as the basis of an accelerated introdu...

Partial differential equations with numerical methods

CERN Document Server

Larsson, Stig

2003-01-01

The book is suitable for advanced undergraduate and beginning graduate students of applied mathematics and engineering. The main theme is the integration of the theory of linear PDEs and the numerical solution of such equations. For each type of PDE, elliptic, parabolic, and hyperbolic, the text contains one chapter on the mathematical theory of the differential equation, followed by one chapter on finite difference methods and one on finite element methods. As preparation, the two-point boundary value problem and the initial-value problem for ODEs are discussed in separate chapters. There is also one chapter on the elliptic eigenvalue problem and eigenfunction expansion. The presentation does not presume a deep knowledge of mathematical and functional analysis. Some background on linear functional analysis and Sobolev spaces, and also on numerical linear algebra, is reviewed in two appendices.

Hybrid RANS-LES using high order numerical methods

Science.gov (United States)

Henry de Frahan, Marc; Yellapantula, Shashank; Vijayakumar, Ganesh; Knaus, Robert; Sprague, Michael

2017-11-01

Understanding the impact of wind turbine wake dynamics on downstream turbines is particularly important for the design of efficient wind farms. Due to their tractable computational cost, hybrid RANS/LES models are an attractive framework for simulating separation flows such as the wake dynamics behind a wind turbine. High-order numerical methods can be computationally efficient and provide increased accuracy in simulating complex flows. In the context of LES, high-order numerical methods have shown some success in predictions of turbulent flows. However, the specifics of hybrid RANS-LES models, including the transition region between both modeling frameworks, pose unique challenges for high-order numerical methods. In this work, we study the effect of increasing the order of accuracy of the numerical scheme in simulations of canonical turbulent flows using RANS, LES, and hybrid RANS-LES models. We describe the interactions between filtering, model transition, and order of accuracy and their effect on turbulence quantities such as kinetic energy spectra, boundary layer evolution, and dissipation rate. This work was funded by the U.S. Department of Energy, Exascale Computing Project, under Contract No. DE-AC36-08-GO28308 with the National Renewable Energy Laboratory.

Biquartic Finite Volume Element Metho d Based on Lobatto-Guass Structure

Institute of Scientific and Technical Information of China (English)

Gao Yan-ni; Chen Yan-li

2015-01-01

In this paper, a biquartic finite volume element method based on Lobatto-Guass structure is presented for variable coefficient elliptic equation on rectangular partition. Not only the optimal H1 and L2 error estimates but also some super-convergent properties are available and could be proved for this method. The numer-ical results obtained by this finite volume element scheme confirm the validity of the theoretical analysis and the effectiveness of this method.

Numerical analysis of transport phenomena in Y-Ba-Cu-O melt during growth of superconducting crystal Y123 by Czochralski method

Science.gov (United States)

Szmyd, J. S.; Suzuki, K.

2003-10-01

In 1993, at the Superconductivity Research Laboratory (SRL), International Superconductivity Technology Centre (ISTEC), in Tokyo, continuous growth of large single crystals of YBa 2Cu 3O 7- x (Y123) was achieved by the application of a modified Czochralski method. This paper presents the numerical computations of the flow, thermal and Y concentration fields in the Ba-Cu-O melt for Y123 single crystal growth by this modified method. The finite volume method was used to calculate the fluid flow, heat transfer and yttrium distribution in the melt with staggered numerical grid. The flow in the melt was modelled as an incompressible Newtonian and Boussinesque fluid. Calculations are presented for a combined flow regime of buoyancy-driven natural convection and crystal-rotation-driven forced convection.

Application of the photoelastic experimental hybrid method with new numerical method to the high stress distribution

International Nuclear Information System (INIS)

Hawong, Jai Sug; Lee, Dong Hun; Lee, Dong Ha; Tche, Konstantin

2004-01-01

In this research, the photoelastic experimental hybrid method with Hook-Jeeves numerical method has been developed: This method is more precise and stable than the photoelastic experimental hybrid method with Newton-Rapson numerical method with Gaussian elimination method. Using the photoelastic experimental hybrid method with Hook-Jeeves numerical method, we can separate stress components from isochromatics only and stress intensity factors and stress concentration factors can be determined. The photoelastic experimental hybrid method with Hook-Jeeves had better be used in the full field experiment than the photoelastic experimental hybrid method with Newton-Rapson with Gaussian elimination method

Numerical simulation of GEW equation using RBF collocation method

Directory of Open Access Journals (Sweden)

Hamid Panahipour

2012-08-01

Full Text Available The generalized equal width (GEW equation is solved numerically by a meshless method based on a global collocation with standard types of radial basis functions (RBFs. Test problems including propagation of single solitons, interaction of two and three solitons, development of the Maxwellian initial condition pulses, wave undulation and wave generation are used to indicate the efficiency and accuracy of the method. Comparisons are made between the results of the proposed method and some other published numerical methods.

Monte Carlo method for critical systems in infinite volume: The planar Ising model.

Science.gov (United States)

Herdeiro, Victor; Doyon, Benjamin

2016-10-01

In this paper we propose a Monte Carlo method for generating finite-domain marginals of critical distributions of statistical models in infinite volume. The algorithm corrects the problem of the long-range effects of boundaries associated to generating critical distributions on finite lattices. It uses the advantage of scale invariance combined with ideas of the renormalization group in order to construct a type of "holographic" boundary condition that encodes the presence of an infinite volume beyond it. We check the quality of the distribution obtained in the case of the planar Ising model by comparing various observables with their infinite-plane prediction. We accurately reproduce planar two-, three-, and four-point of spin and energy operators. We also define a lattice stress-energy tensor, and numerically obtain the associated conformal Ward identities and the Ising central charge.

Implementation of a revised numerical integration technique into QAD

International Nuclear Information System (INIS)

De Gangi, N.L.

1983-01-01

A technique for numerical integration through a uniform volume source is developed. It is applied to gamma radiation transport shielding problems. The method is based on performing a numerical angular and ray point kernel integration and is incorporated into the QAD-CG computer code (i.e. QAD-UE). Several test problems are analyzed with this technique. Convergence properties of the method are analyzed. Gamma dose rates from a large tank and post LOCA dose rates inside a containment building are evaluated. Results are consistent with data from other methods. The new technique provides several advantages. User setup requirements for large volume source problems are reduced from standard point kernel requirements. Calculational efficiencies are improved. An order of magnitude improvement is seen with a test problem

Teaching numerical methods with IPython notebooks and inquiry-based learning

KAUST Repository

Ketcheson, David I.

2014-01-01

A course in numerical methods should teach both the mathematical theory of numerical analysis and the craft of implementing numerical algorithms. The IPython notebook provides a single medium in which mathematics, explanations, executable code, and visualizations can be combined, and with which the student can interact in order to learn both the theory and the craft of numerical methods. The use of notebooks also lends itself naturally to inquiry-based learning methods. I discuss the motivation and practice of teaching a course based on the use of IPython notebooks and inquiry-based learning, including some specific practical aspects. The discussion is based on my experience teaching a Masters-level course in numerical analysis at King Abdullah University of Science and Technology (KAUST), but is intended to be useful for those who teach at other levels or in industry.

Two numerical methods for mean-field games

KAUST Repository

Gomes, Diogo A.

2016-01-01

Here, we consider numerical methods for stationary mean-field games (MFG) and investigate two classes of algorithms. The first one is a gradient flow method based on the variational characterization of certain MFG. The second one uses monotonicity properties of MFG. We illustrate our methods with various examples, including one-dimensional periodic MFG, congestion problems, and higher-dimensional models.

Two numerical methods for mean-field games

KAUST Repository

Gomes, Diogo A.

2016-01-09

Here, we consider numerical methods for stationary mean-field games (MFG) and investigate two classes of algorithms. The first one is a gradient flow method based on the variational characterization of certain MFG. The second one uses monotonicity properties of MFG. We illustrate our methods with various examples, including one-dimensional periodic MFG, congestion problems, and higher-dimensional models.

Numerical methods in electron magnetic resonance

International Nuclear Information System (INIS)

Soernes, A.R.

1998-01-01

The focal point of the thesis is the development and use of numerical methods in the analysis, simulation and interpretation of Electron Magnetic Resonance experiments on free radicals in solids to uncover the structure, the dynamics and the environment of the system

Numerical methods in electron magnetic resonance

Energy Technology Data Exchange (ETDEWEB)

Soernes, A.R

1998-07-01

The focal point of the thesis is the development and use of numerical methods in the analysis, simulation and interpretation of Electron Magnetic Resonance experiments on free radicals in solids to uncover the structure, the dynamics and the environment of the system.

Numerical Modeling of Ablation Heat Transfer

Science.gov (United States)

Ewing, Mark E.; Laker, Travis S.; Walker, David T.

2013-01-01

A unique numerical method has been developed for solving one-dimensional ablation heat transfer problems. This paper provides a comprehensive description of the method, along with detailed derivations of the governing equations. This methodology supports solutions for traditional ablation modeling including such effects as heat transfer, material decomposition, pyrolysis gas permeation and heat exchange, and thermochemical surface erosion. The numerical scheme utilizes a control-volume approach with a variable grid to account for surface movement. This method directly supports implementation of nontraditional models such as material swelling and mechanical erosion, extending capabilities for modeling complex ablation phenomena. Verifications of the numerical implementation are provided using analytical solutions, code comparisons, and the method of manufactured solutions. These verifications are used to demonstrate solution accuracy and proper error convergence rates. A simple demonstration of a mechanical erosion (spallation) model is also provided to illustrate the unique capabilities of the method.

Viscous wing theory development. Volume 1: Analysis, method and results

Science.gov (United States)

Chow, R. R.; Melnik, R. E.; Marconi, F.; Steinhoff, J.

1986-01-01

Viscous transonic flows at large Reynolds numbers over 3-D wings were analyzed using a zonal viscid-inviscid interaction approach. A new numerical AFZ scheme was developed in conjunction with the finite volume formulation for the solution of the inviscid full-potential equation. A special far-field asymptotic boundary condition was developed and a second-order artificial viscosity included for an improved inviscid solution methodology. The integral method was used for the laminar/turbulent boundary layer and 3-D viscous wake calculation. The interaction calculation included the coupling conditions of the source flux due to the wing surface boundary layer, the flux jump due to the viscous wake, and the wake curvature effect. A method was also devised incorporating the 2-D trailing edge strong interaction solution for the normal pressure correction near the trailing edge region. A fully automated computer program was developed to perform the proposed method with one scalar version to be used on an IBM-3081 and two vectorized versions on Cray-1 and Cyber-205 computers.

Numerical computer methods part E

CERN Document Server

Johnson, Michael L

2004-01-01

The contributions in this volume emphasize analysis of experimental data and analytical biochemistry, with examples taken from biochemistry. They serve to inform biomedical researchers of the modern data analysis methods that have developed concomitantly with computer hardware. Selected Contents: A practical approach to interpretation of SVD results; modeling of oscillations in endocrine networks with feedback; quantifying asynchronous breathing; sample entropy; wavelet modeling and processing of nasal airflow traces.

Recent advances in two-phase flow numerics

Energy Technology Data Exchange (ETDEWEB)

Mahaffy, J.H.; Macian, R. [Pennsylvania State Univ., University Park, PA (United States)

1997-07-01

The authors review three topics in the broad field of numerical methods that may be of interest to individuals modeling two-phase flow in nuclear power plants. The first topic is iterative solution of linear equations created during the solution of finite volume equations. The second is numerical tracking of macroscopic liquid interfaces. The final area surveyed is the use of higher spatial difference techniques.

Recent advances in two-phase flow numerics

International Nuclear Information System (INIS)

Mahaffy, J.H.; Macian, R.

1997-01-01

The authors review three topics in the broad field of numerical methods that may be of interest to individuals modeling two-phase flow in nuclear power plants. The first topic is iterative solution of linear equations created during the solution of finite volume equations. The second is numerical tracking of macroscopic liquid interfaces. The final area surveyed is the use of higher spatial difference techniques

Numerical methods for the Lévy LIBOR model

DEFF Research Database (Denmark)

Papapantoleon, Antonis; Skovmand, David

2010-01-01

but the methods are generally slow. We propose an alternative approximation scheme based on Picard iterations. Our approach is similar in accuracy to the full numerical solution, but with the feature that each rate is, unlike the standard method, evolved independently of the other rates in the term structure....... This enables simultaneous calculation of derivative prices of different maturities using parallel computing. We include numerical illustrations of the accuracy and speed of our method pricing caplets.......The aim of this work is to provide fast and accurate approximation schemes for the Monte-Carlo pricing of derivatives in the L\\'evy LIBOR model of Eberlein and \\"Ozkan (2005). Standard methods can be applied to solve the stochastic differential equations of the successive LIBOR rates...

Numerical Methods for the Lévy LIBOR Model

DEFF Research Database (Denmark)

Papapantoleon, Antonis; Skovmand, David

are generally slow. We propose an alternative approximation scheme based on Picard iterations. Our approach is similar in accuracy to the full numerical solution, but with the feature that each rate is, unlike the standard method, evolved independently of the other rates in the term structure. This enables...... simultaneous calculation of derivative prices of different maturities using parallel computing. We include numerical illustrations of the accuracy and speed of our method pricing caplets.......The aim of this work is to provide fast and accurate approximation schemes for the Monte-Carlo pricing of derivatives in the Lévy LIBOR model of Eberlein and Özkan (2005). Standard methods can be applied to solve the stochastic differential equations of the successive LIBOR rates but the methods...

Volume Sculpting Using the Level-Set Method

DEFF Research Database (Denmark)

Bærentzen, Jakob Andreas; Christensen, Niels Jørgen

2002-01-01

In this paper, we propose the use of the Level--Set Method as the underlying technology of a volume sculpting system. The main motivation is that this leads to a very generic technique for deformation of volumetric solids. In addition, our method preserves a distance field volume representation....... A scaling window is used to adapt the Level--Set Method to local deformations and to allow the user to control the intensity of the tool. Level--Set based tools have been implemented in an interactive sculpting system, and we show sculptures created using the system....

Taylor bubbles at high viscosity ratios: experiments and numerical simulations

Science.gov (United States)

Hewakandamby, Buddhika; Hasan, Abbas; Azzopardi, Barry; Xie, Zhihua; Pain, Chris; Matar, Omar

2015-11-01

The Taylor bubble is a single long bubble which nearly fills the entire cross section of a liquid-filled circular tube, often occurring in gas-liquid slug flows in many industrial applications, particularly oil and gas production. The objective of this study is to investigate the fluid dynamics of three-dimensional Taylor bubble rising in highly viscous silicone oil in a vertical pipe. An adaptive unstructured mesh modelling framework is adopted here which can modify and adapt anisotropic unstructured meshes to better represent the underlying physics of bubble rising and reduce computational effort without sacrificing accuracy. The numerical framework consists of a mixed control volume and finite element formulation, a `volume of fluid'-type method for the interface-capturing based on a compressive control volume advection method, and a force-balanced algorithm for the surface tension implementation. Experimental results for the Taylor bubble shape and rise velocity are presented, together with numerical results for the dynamics of the bubbles. A comparison of the simulation predictions with experimental data available in the literature is also presented to demonstrate the capabilities of our numerical method. EPSRC Programme Grant, MEMPHIS, EP/K0039761/1.

On mesh refinement and accuracy of numerical solutions

NARCIS (Netherlands)

Zhou, Hong; Peters, Maria; van Oosterom, Adriaan

1993-01-01

This paper investigates mesh refinement and its relation with the accuracy of the boundary element method (BEM) and the finite element method (FEM). TO this end an isotropic homogeneous spherical volume conductor, for which the analytical solution is available, wag used. The numerical results

Classical and modern numerical analysis theory, methods and practice

CERN Document Server

Ackleh, Azmy S; Kearfott, R Baker; Seshaiyer, Padmanabhan

2009-01-01

Mathematical Review and Computer Arithmetic Mathematical Review Computer Arithmetic Interval ComputationsNumerical Solution of Nonlinear Equations of One Variable Introduction Bisection Method The Fixed Point Method Newton's Method (Newton-Raphson Method) The Univariate Interval Newton MethodSecant Method and Müller's Method Aitken Acceleration and Steffensen's Method Roots of Polynomials Additional Notes and SummaryNumerical Linear Algebra Basic Results from Linear Algebra Normed Linear Spaces Direct Methods for Solving Linear SystemsIterative Methods for Solving Linear SystemsThe Singular Value DecompositionApproximation TheoryIntroduction Norms, Projections, Inner Product Spaces, and Orthogonalization in Function SpacesPolynomial ApproximationPiecewise Polynomial ApproximationTrigonometric ApproximationRational ApproximationWavelet BasesLeast Squares Approximation on a Finite Point SetEigenvalue-Eigenvector Computation Basic Results from Linear Algebra The Power Method The Inverse Power Method Deflation T...

Numerical method for the nonlinear Fokker-Planck equation

International Nuclear Information System (INIS)

Zhang, D.S.; Wei, G.W.; Kouri, D.J.; Hoffman, D.K.

1997-01-01

A practical method based on distributed approximating functionals (DAFs) is proposed for numerically solving a general class of nonlinear time-dependent Fokker-Planck equations. The method relies on a numerical scheme that couples the usual path-integral concept to the DAF idea. The high accuracy and reliability of the method are illustrated by applying it to an exactly solvable nonlinear Fokker-Planck equation, and the method is compared with the accurate K-point Stirling interpolation formula finite-difference method. The approach is also used successfully to solve a nonlinear self-consistent dynamic mean-field problem for which both the cumulant expansion and scaling theory have been found by Drozdov and Morillo [Phys. Rev. E 54, 931 (1996)] to be inadequate to describe the occurrence of a long-lived transient bimodality. The standard interpretation of the transient bimodality in terms of the flat region in the kinetic potential fails for the present case. An alternative analysis based on the effective potential of the Schroedinger-like Fokker-Planck equation is suggested. Our analysis of the transient bimodality is strongly supported by two examples that are numerically much more challenging than other examples that have been previously reported for this problem. copyright 1997 The American Physical Society

FINITE VOLUME METHOD FOR SOLVING THREE-DIMENSIONAL ELECTRIC FIELD DISTRIBUTION

Directory of Open Access Journals (Sweden)

Paţiuc V.I.

2011-04-01

Full Text Available The paper examines a new approach to finite volume method which is used to calculate the electric field spatially homogeneous three-dimensional environment. It is formulated the problem Dirihle with building of the computational grid on base of space partition, which is known as Delone triangulation with the use of Voronoi cells. It is proposed numerical algorithm for calculating the potential and electric field strength in the space formed by a cylinder placed in the air. It is developed algorithm and software which were for the case, when the potential on the inner surface of the cylinder has been assigned and on the outer surface and the bottom of cylinder it was assigned zero potential. There are presented results of calculations of distribution in the potential space and electric field strength.

A Numerical Method for Lane-Emden Equations Using Hybrid Functions and the Collocation Method

Directory of Open Access Journals (Sweden)

Changqing Yang

2012-01-01

Full Text Available A numerical method to solve Lane-Emden equations as singular initial value problems is presented in this work. This method is based on the replacement of unknown functions through a truncated series of hybrid of block-pulse functions and Chebyshev polynomials. The collocation method transforms the differential equation into a system of algebraic equations. It also has application in a wide area of differential equations. Corresponding numerical examples are presented to demonstrate the accuracy of the proposed method.

Numerical methods for modeling photonic-crystal VCSELs

DEFF Research Database (Denmark)

Dems, Maciej; Chung, Il-Sug; Nyakas, Peter

2010-01-01

We show comparison of four different numerical methods for simulating Photonic-Crystal (PC) VCSELs. We present the theoretical basis behind each method and analyze the differences by studying a benchmark VCSEL structure, where the PC structure penetrates all VCSEL layers, the entire top-mirror DBR...... to the effective index method. The simulation results elucidate the strength and weaknesses of the analyzed methods; and outline the limits of applicability of the different models....

A parallel finite-volume finite-element method for transient compressible turbulent flows with heat transfer

International Nuclear Information System (INIS)

Masoud Ziaei-Rad

2010-01-01

In this paper, a two-dimensional numerical scheme is presented for the simulation of turbulent, viscous, transient compressible flows in the simultaneously developing hydraulic and thermal boundary layer region. The numerical procedure is a finite-volume-based finite-element method applied to unstructured grids. This combination together with a new method applied for the boundary conditions allows for accurate computation of the variables in the entrance region and for a wide range of flow fields from subsonic to transonic. The Roe-Riemann solver is used for the convective terms, whereas the standard Galerkin technique is applied for the viscous terms. A modified κ-ε model with a two-layer equation for the near-wall region combined with a compressibility correction is used to predict the turbulent viscosity. Parallel processing is also employed to divide the computational domain among the different processors to reduce the computational time. The method is applied to some test cases in order to verify the numerical accuracy. The results show significant differences between incompressible and compressible flows in the friction coefficient, Nusselt number, shear stress and the ratio of the compressible turbulent viscosity to the molecular viscosity along the developing region. A transient flow generated after an accidental rupture in a pipeline was also studied as a test case. The results show that the present numerical scheme is stable, accurate and efficient enough to solve the problem of transient wall-bounded flow.

Numerical Cerebrospinal System Modeling in Fluid-Structure Interaction.

Science.gov (United States)

Garnotel, Simon; Salmon, Stéphanie; Balédent, Olivier

2018-01-01

Cerebrospinal fluid (CSF) stroke volume in the aqueduct is widely used to evaluate CSF dynamics disorders. In a healthy population, aqueduct stroke volume represents around 10% of the spinal stroke volume while intracranial subarachnoid space stroke volume represents 90%. The amplitude of the CSF oscillations through the different compartments of the cerebrospinal system is a function of the geometry and the compliances of each compartment, but we suspect that it could also be impacted be the cardiac cycle frequency. To study this CSF distribution, we have developed a numerical model of the cerebrospinal system taking into account cerebral ventricles, intracranial subarachnoid spaces, spinal canal and brain tissue in fluid-structure interactions. A numerical fluid-structure interaction model is implemented using a finite-element method library to model the cerebrospinal system and its interaction with the brain based on fluid mechanics equations and linear elasticity equations coupled in a monolithic formulation. The model geometry, simplified in a first approach, is designed in accordance with realistic volume ratios of the different compartments: a thin tube is used to mimic the high flow resistance of the aqueduct. CSF velocity and pressure and brain displacements are obtained as simulation results, and CSF flow and stroke volume are calculated from these results. Simulation results show a significant variability of aqueduct stroke volume and intracranial subarachnoid space stroke volume in the physiological range of cardiac frequencies. Fluid-structure interactions are numerous in the cerebrospinal system and difficult to understand in the rigid skull. The presented model highlights significant variations of stroke volumes under cardiac frequency variations only.

Design of steel energy absorbing restrainers and their incorporation into nuclear power plants for enhanced safety. Volume 1C. Numerical method for dynamic substructure analysis

International Nuclear Information System (INIS)

Dickens, J.M.; Wilson, E.L.

1980-06-01

This report presents several new or improved numerical methods for the dynamic response analysis of large complex structural systems. The methods presented are primarily concerned with linear dynamic analysis. However, one of the significant results of the research is the development of efficient numerical methods for large systems with a small number of non-linear members. The technique of dynamic substructure analysis is used to accurately model the behavior of the linear part of the system. This approach allows the behavior of the linear structural members to be reduced to a small number of generalized dynamic coordinates. Those linear dynamic coordinates are then used in direct combination with the nonlinear members to approximate the solution of the large nonlinear problem. The methods presented can be applied to structures of arbitrary geometry. However, the use of techniques such as cyclic symmetry can further decrease the numerical effort

Coupled Finite Volume and Finite Element Method Analysis of a Complex Large-Span Roof Structure

Science.gov (United States)

Szafran, J.; Juszczyk, K.; Kamiński, M.

2017-12-01

The main goal of this paper is to present coupled Computational Fluid Dynamics and structural analysis for the precise determination of wind impact on internal forces and deformations of structural elements of a longspan roof structure. The Finite Volume Method (FVM) serves for a solution of the fluid flow problem to model the air flow around the structure, whose results are applied in turn as the boundary tractions in the Finite Element Method problem structural solution for the linear elastostatics with small deformations. The first part is carried out with the use of ANSYS 15.0 computer system, whereas the FEM system Robot supports stress analysis in particular roof members. A comparison of the wind pressure distribution throughout the roof surface shows some differences with respect to that available in the engineering designing codes like Eurocode, which deserves separate further numerical studies. Coupling of these two separate numerical techniques appears to be promising in view of future computational models of stochastic nature in large scale structural systems due to the stochastic perturbation method.

A numerical method for transient gas-liquid two-phase flow using a general curvilinear coordinate system. 1. Governing equations and numerical method

International Nuclear Information System (INIS)

Tomiyama, Akio; Matsuoka, Toshiyuki.

1995-01-01

A simple numerical method for solving a transient incompressible two-fluid model was proposed in the present study. A general curvilinear coordinate system was adopted in this method for predicting transient flows in practical engineering devices. The simplicity of the present method is due to the fact that the field equations and constitutive equations were expressed in a tensor form in the general curvilinear coordinate system. When a conventional rectangular mesh is adopted in a calculation, the method reduces to a numerical method for a Cartesian coordinate system. As an example, the present method was applied to transient air-water bubbly flow in a vertical U-tube. It was confirmed that the effects of centrifugal and gravitational forces on the phase distribution in the U-tube were reasonably predicted. (author)

High accuracy mantle convection simulation through modern numerical methods

KAUST Repository

Kronbichler, Martin

2012-08-21

Numerical simulation of the processes in the Earth\\'s mantle is a key piece in understanding its dynamics, composition, history and interaction with the lithosphere and the Earth\\'s core. However, doing so presents many practical difficulties related to the numerical methods that can accurately represent these processes at relevant scales. This paper presents an overview of the state of the art in algorithms for high-Rayleigh number flows such as those in the Earth\\'s mantle, and discusses their implementation in the Open Source code Aspect (Advanced Solver for Problems in Earth\\'s ConvecTion). Specifically, we show how an interconnected set of methods for adaptive mesh refinement (AMR), higher order spatial and temporal discretizations, advection stabilization and efficient linear solvers can provide high accuracy at a numerical cost unachievable with traditional methods, and how these methods can be designed in a way so that they scale to large numbers of processors on compute clusters. Aspect relies on the numerical software packages deal.II and Trilinos, enabling us to focus on high level code and keeping our implementation compact. We present results from validation tests using widely used benchmarks for our code, as well as scaling results from parallel runs. © 2012 The Authors Geophysical Journal International © 2012 RAS.

Computing the partial volume of pressure vessels

Energy Technology Data Exchange (ETDEWEB)

Wiencke, Bent [Nestle USA, Corporate Engineering, 800 N. Brand Blvd, Glendale, CA 91203 (United States)

2010-06-15

The computation of the partial and total volume of pressure vessels with various type of head profiles requires detailed knowledge of the head profile geometry. Depending on the type of head profile the derivation of the equations can become very complex and the calculation process cumbersome. Certain head profiles require numerical methods to obtain the partial volume, which for most application is beyond the scope of practicability. This paper suggests a unique method that simplifies the calculation procedure for the various types of head profiles by using one common set of equations without the need for numerical or complex computation methods. For ease of use, all equations presented in this paper are summarized in a single table format for horizontal and vertical vessels. (author)

Numerical methods: Analytical benchmarking in transport theory

International Nuclear Information System (INIS)

Ganapol, B.D.

1988-01-01

Numerical methods applied to reactor technology have reached a high degree of maturity. Certainly one- and two-dimensional neutron transport calculations have become routine, with several programs available on personal computer and the most widely used programs adapted to workstation and minicomputer computational environments. With the introduction of massive parallelism and as experience with multitasking increases, even more improvement in the development of transport algorithms can be expected. Benchmarking an algorithm is usually not a very pleasant experience for the code developer. Proper algorithmic verification by benchmarking involves the following considerations: (1) conservation of particles, (2) confirmation of intuitive physical behavior, and (3) reproduction of analytical benchmark results. By using today's computational advantages, new basic numerical methods have been developed that allow a wider class of benchmark problems to be considered

A numerical method for a transient two-fluid model

International Nuclear Information System (INIS)

Le Coq, G.; Libmann, M.

1978-01-01

The transient boiling two-phase flow is studied. In nuclear reactors, the driving conditions for the transient boiling are a pump power decay or/and an increase in heating power. The physical model adopted for the two-phase flow is the two fluid model with the assumption that the vapor remains at saturation. The numerical method for solving the thermohydraulics problems is a shooting method, this method is highly implicit. A particular problem exists at the boiling and condensation front. A computer code using this numerical method allow the calculation of a transient boiling initiated by a steady state for a PWR or for a LMFBR

Numerical Simulation of Two-Fluid Mingling Using the Particle Finite Element Method with Applications to Magmatic and Volcanic Processes

Science.gov (United States)

de Mier, M.; Costa, F.; Idelsohn, S.

2008-12-01

Many magmatic and volcanic processes (e.g., magma differentiation, mingling, transport in the volcanic conduit) are controlled by the physical properties and flow styles of high-temperature silicate melts. Such processes can be experimentally investigated using analog systems and scaling methods, but it is difficult to find the suitable material and it is generally not possible to quantitatively extrapolate the results to the natural system. An alternative means of studying fluid dynamics in volcanic systems is with numerical models. We have chosen the Particle Finite Element Method (PFEM), which is based on a Delaunay mesh that moves with the fluid velocity, the Navier-Stokes equations in Lagrangian formulation, and linear elements for velocity, pressure, and temperature. Remeshing is performed when the grid becomes too distorted [E. Oñate et al., 2004. The Particle Finite Element Method: An Overview. Int. J. Comput. Meth. 1, 267-307]. The method is ideal for tracking material interfaces between different fluids or media. Methods based on Eulerian reference frames need special techniques, such as level-set or volume-of-fluid, to capture the interface position, and these techniques add a significant numerical diffusion at the interface. We have performed a series of two-dimensional simulations of a classical problem of fluid dynamics in magmatic and volcanic systems: intrusion of a basaltic melt in a silica-rich magma reservoir. We have used realistic physical properties and equations of state for the silicate melts (e.g., temperature, viscosity, and density) and tracked the changes in the system for geologically relevant time scales (up to 100 years). The problem is modeled by the low-Mach-number equations derived from an asymptotic analysis of the compressible Navier-Stokes equations that removes shock waves from the flow but allows however large variations of density due to temperature variations. Non-constant viscosity and volume changes are taken into account

Fast numerical upscaling of heat equation for fibrous materials

KAUST Repository

Iliev, Oleg; Lazarov, Raytcho; Willems, Joerg

2010-01-01

We are interested in numerical methods for computing the effective heat conductivities of fibrous insulation materials, such as glass or mineral wool, characterized by low solid volume fractions and high contrasts, i.e., high ratios between the thermal conductivities of the fibers and the surrounding air. We consider a fast numerical method for solving some auxiliary cell problems appearing in this upscaling procedure. The auxiliary problems are boundary value problems of the steady-state heat equation in a representative elementary volume occupied by fibers and air. We make a simplification by replacing these problems with appropriate boundary value problems in the domain occupied by the fibers only. Finally, the obtained problems are further simplified by taking advantage of the slender shape of the fibers and assuming that they form a network. A discretization on the graph defined by the fibers is presented and error estimates are provided. The resulting algorithm is discussed and the accuracy and the performance of the method are illusrated on a number of numerical experiments. © Springer-Verlag 2010.

Fast numerical upscaling of heat equation for fibrous materials

KAUST Repository

Iliev, Oleg

2010-08-01

We are interested in numerical methods for computing the effective heat conductivities of fibrous insulation materials, such as glass or mineral wool, characterized by low solid volume fractions and high contrasts, i.e., high ratios between the thermal conductivities of the fibers and the surrounding air. We consider a fast numerical method for solving some auxiliary cell problems appearing in this upscaling procedure. The auxiliary problems are boundary value problems of the steady-state heat equation in a representative elementary volume occupied by fibers and air. We make a simplification by replacing these problems with appropriate boundary value problems in the domain occupied by the fibers only. Finally, the obtained problems are further simplified by taking advantage of the slender shape of the fibers and assuming that they form a network. A discretization on the graph defined by the fibers is presented and error estimates are provided. The resulting algorithm is discussed and the accuracy and the performance of the method are illusrated on a number of numerical experiments. © Springer-Verlag 2010.

Method for numerical simulation of two-term exponentially correlated colored noise

International Nuclear Information System (INIS)

Yilmaz, B.; Ayik, S.; Abe, Y.; Gokalp, A.; Yilmaz, O.

2006-01-01

A method for numerical simulation of two-term exponentially correlated colored noise is proposed. The method is an extension of traditional method for one-term exponentially correlated colored noise. The validity of the algorithm is tested by comparing numerical simulations with analytical results in two physical applications

Mathematica with a Numerical Methods Course

Science.gov (United States)

Varley, Rodney

2003-04-01

An interdisciplinary "Numerical Methods" course has been shared between physics, mathematics and computer science since 1992 at Hunter C. Recently, the lectures and workshops for this course have become formalized and placed on the internet at http://www.ph.hunter.cuny.edu (follow the links "Course Listings and Websites" >> "PHYS385 (Numerical Methods)". Mathematica notebooks for the lectures are available for automatic download (by "double clicking" the lecture icon) for student use in the classroom or at home. AOL (or Netscape/Explorer) can be used provided Mathematica (or the "free" MathReader) has been made a "helper application". Using Mathematica has the virtue that mathematical equations (no LaTex required) can easily be included with the text and Mathematica's graphing is easy to use. Computational cells can be included within the notebook and students may easily modify the calculation to see the result of "what if..." questions. Homework is sent as Mathematica notebooks to the instructor via the internet and the corrected workshops are returned in the same manner. Most exam questions require computational solutions.

Numerical perturbative methods in the quantum theory of physical systems

International Nuclear Information System (INIS)

Adam, G.

1980-01-01

During the last two decades, development of digital electronic computers has led to the deployment of new, distinct methods in theoretical physics. These methods, based on the advances of modern numerical analysis as well as on specific equations describing physical processes, enabled to perform precise calculations of high complexity which have completed and sometimes changed our image of many physical phenomena. Our efforts have concentrated on the development of numerical methods with such intrinsic performances as to allow a successful approach of some Key issues in present theoretical physics on smaller computation systems. The basic principle of such methods is to translate, in numerical analysis language, the theory of perturbations which is suited to numerical rather than to analytical computation. This idea has been illustrated by working out two problems which arise from the time independent Schroedinger equation in the non-relativistic approximation, within both quantum systems with a small number of particles and systems with a large number of particles, respectively. In the first case, we are led to the numerical solution of some quadratic ordinary differential equations (first section of the thesis) and in the second case, to the solution of some secular equations in the Brillouin area (second section). (author)

Well balanced finite volume methods for nearly hydrostatic flows

International Nuclear Information System (INIS)

Botta, N.; Klein, R.; Langenberg, S.; Luetzenkirchen, S.

2004-01-01

In numerical approximations of nearly hydrostatic flows, a proper representation of the dominant hydrostatic balance is of crucial importance: unbalanced truncation errors can induce unacceptable spurious motions, e.g., in dynamical cores of models for numerical weather prediction (NWP) in particular near steep topography. In this paper we develop a new strategy for the construction of discretizations that are 'well-balanced' with respect to dominant hydrostatics. The classical idea of formulating the momentum balance in terms of deviations of pressure from a balanced background distribution is realized here through local, time dependent hydrostatic reconstructions. Balanced discretizations of the pressure gradient and of the gravitation source term are achieved through a 'discrete Archimedes' buoyancy principle'. This strategy is applied to extend an explicit standard finite volume Godunov-type scheme for compressible flows with minimal modifications. The resulting method has the following features: (i) It inherits its conservation properties from the underlying base scheme. (ii) It is exactly balanced, even on curvilinear grids, for a large class of near-hydrostatic flows. (iii) It solves the full compressible flow equations without reference to a background state that is defined for an entire vertical column of air. (iv) It is robust with respect to details of the implementation, such as the choice of slope limiting functions, or the particularities of boundary condition discretizations

Method and apparatus for imaging volume data

International Nuclear Information System (INIS)

Drebin, R.; Carpenter, L.C.

1987-01-01

An imaging system projects a two dimensional representation of three dimensional volumes where surface boundaries and objects internal to the volumes are readily shown, and hidden surfaces and the surface boundaries themselves are accurately rendered by determining volume elements or voxels. An image volume representing a volume object or data structure is written into memory. A color and opacity is assigned to each voxel within the volume and stored as a red (R), green (G), blue (B), and opacity (A) component, three dimensional data volume. The RGBA assignment for each voxel is determined based on the percentage component composition of the materials represented in the volume, and thus, the percentage of color and transparency associated with those materials. The voxels in the RGBA volume are used as mathematical filters such that each successive voxel filter is overlayed over a prior background voxel filter. Through a linear interpolation, a new background filter is determined and generated. The interpolation is successively performed for all voxels up to the front most voxel for the plane of view. The method is repeated until all display voxels are determined for the plane of view. (author)

Numerical Calculation of Transport Based on the Drift-Kinetic Equation for Plasmas in General Toroidal Magnetic Geometry: Numerical Methods

International Nuclear Information System (INIS)

Reynolds, J. M.; Lopez-Bruna, D.

2009-01-01

In this report we continue with the description of a newly developed numerical method to solve the drift kinetic equation for ions and electrons in toroidal plasmas. Several numerical aspects, already outlined in a previous report [Informes Tecnicos Ciemat 1165, mayo 2009], will be treated now in more detail. Aside from discussing the method in the context of other existing codes, various aspects will be now explained from the viewpoint of numerical methods: the way to solve convection equations, the adopted boundary conditions, the real-space meshing procedures along with a new software developed to build them, and some additional questions related with the parallelization and the numerical integration. (Author) 16 refs

Towards numerical simulations of supersonic liquid jets using ghost fluid method

International Nuclear Information System (INIS)

Majidi, Sahand; Afshari, Asghar

2015-01-01

Highlights: • A ghost fluid method based solver is developed for numerical simulation of compressible multiphase flows. • The performance of the numerical tool is validated via several benchmark problems. • Emergence of supersonic liquid jets in quiescent gaseous environment is simulated using ghost fluid method for the first time. • Bow-shock formation ahead of the liquid jet is clearly observed in the obtained numerical results. • Radiation of mach waves from the phase-interface witnessed experimentally is evidently captured in our numerical simulations. - Abstract: A computational tool based on the ghost fluid method (GFM) is developed to study supersonic liquid jets involving strong shocks and contact discontinuities with high density ratios. The solver utilizes constrained reinitialization method and is capable of switching between the exact and approximate Riemann solvers to increase the robustness. The numerical methodology is validated through several benchmark test problems; these include one-dimensional multiphase shock tube problem, shock–bubble interaction, air cavity collapse in water, and underwater-explosion. A comparison between our results and numerical and experimental observations indicate that the developed solver performs well investigating these problems. The code is then used to simulate the emergence of a supersonic liquid jet into a quiescent gaseous medium, which is the very first time to be studied by a ghost fluid method. The results of simulations are in good agreement with the experimental investigations. Also some of the famous flow characteristics, like the propagation of pressure-waves from the liquid jet interface and dependence of the Mach cone structure on the inlet Mach number, are reproduced numerically. The numerical simulations conducted here suggest that the ghost fluid method is an affordable and reliable scheme to study complicated interfacial evolutions in complex multiphase systems such as supersonic liquid

A solution of two-dimensional magnetohydrodynamic flow using the finite volume method

Directory of Open Access Journals (Sweden)

Naceur Sonia

2014-01-01

Full Text Available This paper presents the two dimensional numerical modeling of the coupling electromagnetic-hydrodynamic phenomena in a conduction MHD pump using the Finite volume Method. Magnetohydrodynamic problems are, thus, interdisciplinary and coupled, since the effect of the velocity field appears in the magnetic transport equations, and the interaction between the electric current and the magnetic field appears in the momentum transport equations. The resolution of the Maxwell's and Navier Stokes equations is obtained by introducing the magnetic vector potential A, the vorticity z and the stream function y. The flux density, the electromagnetic force, and the velocity are graphically presented. Also, the simulation results agree with those obtained by Ansys Workbench Fluent software.

Development of a high-order finite volume method with multiblock partition techniques

Directory of Open Access Journals (Sweden)

E. M. Lemos

2012-03-01

Full Text Available This work deals with a new numerical methodology to solve the Navier-Stokes equations based on a finite volume method applied to structured meshes with co-located grids. High-order schemes used to approximate advective, diffusive and non-linear terms, connected with multiblock partition techniques, are the main contributions of this paper. Combination of these two techniques resulted in a computer code that involves high accuracy due the high-order schemes and great flexibility to generate locally refined meshes based on the multiblock approach. This computer code has been able to obtain results with higher or equal accuracy in comparison with results obtained using classical procedures, with considerably less computational effort.

Numerical methods in matrix computations

CERN Document Server

Björck, Åke

2015-01-01

Matrix algorithms are at the core of scientific computing and are indispensable tools in most applications in engineering. This book offers a comprehensive and up-to-date treatment of modern methods in matrix computation. It uses a unified approach to direct and iterative methods for linear systems, least squares and eigenvalue problems. A thorough analysis of the stability, accuracy, and complexity of the treated methods is given. Numerical Methods in Matrix Computations is suitable for use in courses on scientific computing and applied technical areas at advanced undergraduate and graduate level. A large bibliography is provided, which includes both historical and review papers as well as recent research papers. This makes the book useful also as a reference and guide to further study and research work. Åke Björck is a professor emeritus at the Department of Mathematics, Linköping University. He is a Fellow of the Society of Industrial and Applied Mathematics.

Numerical simulation of a bubble rising in an environment consisting of Xanthan gum

Science.gov (United States)

Aguirre, Víctor A.; Castillo, Byron A.; Narvaez, Christian P.

2017-09-01

An improved numerical algorithm for front tracking method is developed to simulate a bubble rising in viscous liquid. In the new numerical algorithm, volume correction is introduced to conserve the bubble volume while tracking the bubble's rising and deforming. Volume flux conservation is adopted to solve the Navier-Stokes equation for fluid flow using finite volume method. Non-Newtonian fluids are widely used in industry such as feed and energy industries. In this research we used Xanthan gum which is a microbiological polysaccharide. In order to obtain the properties of the Xanthan gum, such as viscosity, storage and loss modulus, shear rate, etc., it was necessary to do an amplitude sweep and steady flow test in a rheometer with a concentric cylinder as geometry. Based on the data given and using a numerical regression, the coefficients required by Giesekus model are obtained. With these coefficients, it is possible to simulate the comportment of the fluid by the use of the developed algorithm. Once the data given by OpenFOAM is acquired, it is compared with the experimental data.

Hybrid numerical calculation method for bend waveguides

OpenAIRE

Garnier , Lucas; Saavedra , C.; Castro-Beltran , Rigoberto; Lucio , José Luis; Bêche , Bruno

2017-01-01

National audience; The knowledge of how the light will behave in a waveguide with a radius of curvature becomes more and more important because of the development of integrated photonics, which include ring micro-resonators, phasars, and other devices with a radius of curvature. This work presents a numerical calculation method to determine the eigenvalues and eigenvectors of curved waveguides. This method is a hybrid method which uses at first conform transformation of the complex plane gene...

Intelligent numerical methods applications to fractional calculus

CERN Document Server

Anastassiou, George A

2016-01-01

In this monograph the authors present Newton-type, Newton-like and other numerical methods, which involve fractional derivatives and fractional integral operators, for the first time studied in the literature. All for the purpose to solve numerically equations whose associated functions can be also non-differentiable in the ordinary sense. That is among others extending the classical Newton method theory which requires usual differentiability of function. Chapters are self-contained and can be read independently and several advanced courses can be taught out of this book. An extensive list of references is given per chapter. The book’s results are expected to find applications in many areas of applied mathematics, stochastics, computer science and engineering. As such this monograph is suitable for researchers, graduate students, and seminars of the above subjects, also to be in all science and engineering libraries.

Numerical Continuation Methods for Intrusive Uncertainty Quantification Studies

Energy Technology Data Exchange (ETDEWEB)

Safta, Cosmin [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Najm, Habib N. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Phipps, Eric Todd [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

2014-09-01

Rigorous modeling of engineering systems relies on efficient propagation of uncertainty from input parameters to model outputs. In recent years, there has been substantial development of probabilistic polynomial chaos (PC) Uncertainty Quantification (UQ) methods, enabling studies in expensive computational models. One approach, termed ”intrusive”, involving reformulation of the governing equations, has been found to have superior computational performance compared to non-intrusive sampling-based methods in relevant large-scale problems, particularly in the context of emerging architectures. However, the utility of intrusive methods has been severely limited due to detrimental numerical instabilities associated with strong nonlinear physics. Previous methods for stabilizing these constructions tend to add unacceptably high computational costs, particularly in problems with many uncertain parameters. In order to address these challenges, we propose to adapt and improve numerical continuation methods for the robust time integration of intrusive PC system dynamics. We propose adaptive methods, starting with a small uncertainty for which the model has stable behavior and gradually moving to larger uncertainty where the instabilities are rampant, in a manner that provides a suitable solution.

Assembly Discontinuity Factors for the Neutron Diffusion Equation discretized with the Finite Volume Method. Application to BWR

International Nuclear Information System (INIS)

Bernal, A.; Roman, J.E.; Miró, R.; Verdú, G.

2016-01-01

Highlights: • A method is proposed to solve the eigenvalue problem of the Neutron Diffusion Equation in BWR. • The Neutron Diffusion Equation is discretized with the Finite Volume Method. • The currents are calculated by using a polynomial expansion of the neutron flux. • The current continuity and boundary conditions are defined implicitly to reduce the size of the matrices. • Different structured and unstructured meshes were used to discretize the BWR. - Abstract: The neutron flux spatial distribution in Boiling Water Reactors (BWRs) can be calculated by means of the Neutron Diffusion Equation (NDE), which is a space- and time-dependent differential equation. In steady state conditions, the time derivative terms are zero and this equation is rewritten as an eigenvalue problem. In addition, the spatial partial derivatives terms are transformed into algebraic terms by discretizing the geometry and using numerical methods. As regards the geometrical discretization, BWRs are complex systems containing different components of different geometries and materials, but they are usually modelled as parallelepiped nodes each one containing only one homogenized material to simplify the solution of the NDE. There are several techniques to correct the homogenization in the node, but the most commonly used in BWRs is that based on Assembly Discontinuity Factors (ADFs). As regards numerical methods, the Finite Volume Method (FVM) is feasible and suitable to be applied to the NDE. In this paper, a FVM based on a polynomial expansion method has been used to obtain the matrices of the eigenvalue problem, assuring the accomplishment of the ADFs for a BWR. This eigenvalue problem has been solved by means of the SLEPc library.

A Numerical Matrix-Based method in Harmonic Studies in Wind Power Plants

DEFF Research Database (Denmark)

Dowlatabadi, Mohammadkazem Bakhshizadeh; Hjerrild, Jesper; Kocewiak, Łukasz Hubert

2016-01-01

In the low frequency range, there are some couplings between the positive- and negative-sequence small-signal impedances of the power converter due to the nonlinear and low bandwidth control loops such as the synchronization loop. In this paper, a new numerical method which also considers...... these couplings will be presented. The numerical data are advantageous to the parametric differential equations, because analysing the high order and complex transfer functions is very difficult, and finally one uses the numerical evaluation methods. This paper proposes a numerical matrix-based method, which...

Interdisciplinary Study of Numerical Methods and Power Plants Engineering

Directory of Open Access Journals (Sweden)

Ioana OPRIS

2014-08-01

Full Text Available The development of technology, electronics and computing opened the way for a cross-disciplinary research that brings benefits by combining the achievements of different fields. To prepare the students for their future interdisciplinary approach,aninterdisciplinary teaching is adopted. This ensures their progress in knowledge, understanding and ability to navigate through different fields. Aiming these results, the Universities introduce new interdisciplinary courses which explore complex problems by studying subjects from different domains. The paper presents a problem encountered in designingpower plants. The method of solvingthe problem isused to explain the numerical methods and to exercise programming.The goal of understanding a numerical algorithm that solves a linear system of equations is achieved by using the knowledge of heat transfer to design the regenerative circuit of a thermal power plant. In this way, the outcomes from the prior courses (mathematics and physics are used to explain a new subject (numerical methods and to advance future ones (power plants.

Numerical Solution of Fuzzy Differential Equations by Runge-Kutta Verner Method

Directory of Open Access Journals (Sweden)

T. Jayakumar

2015-01-01

Full Text Available In this paper we study the numerical methods for Fuzzy Differential equations by an application of the Runge-Kutta Verner method for fuzzy differential equations. We prove a convergence result and give numerical examples to illustrate the theory.

Method of moments solution of volume integral equations using higher-order hierarchical Legendre basis functions

DEFF Research Database (Denmark)

Kim, Oleksiy S.; Jørgensen, Erik; Meincke, Peter

2004-01-01

An efficient higher-order method of moments (MoM) solution of volume integral equations is presented. The higher-order MoM solution is based on higher-order hierarchical Legendre basis functions and higher-order geometry modeling. An unstructured mesh composed of 8-node trilinear and/or curved 27...... of magnitude in comparison to existing higher-order hierarchical basis functions. Consequently, an iterative solver can be applied even for high expansion orders. Numerical results demonstrate excellent agreement with the analytical Mie series solution for a dielectric sphere as well as with results obtained...

Numerical methods for Bayesian inference in the face of aging

International Nuclear Information System (INIS)

Clarotti, C.A.; Villain, B.; Procaccia, H.

1996-01-01

In recent years, much attention has been paid to Bayesian methods for Risk Assessment. Until now, these methods have been studied from a theoretical point of view. Researchers have been mainly interested in: studying the effectiveness of Bayesian methods in handling rare events; debating about the problem of priors and other philosophical issues. An aspect central to the Bayesian approach is numerical computation because any safety/reliability problem, in a Bayesian frame, ends with a problem of numerical integration. This aspect has been neglected until now because most Risk studies assumed the Exponential model as the basic probabilistic model. The existence of conjugate priors makes numerical integration unnecessary in this case. If aging is to be taken into account, no conjugate family is available and the use of numerical integration becomes compulsory. EDF (National Board of Electricity, of France) and ENEA (National Committee for Energy, New Technologies and Environment, of Italy) jointly carried out a research program aimed at developing quadrature methods suitable for Bayesian Interference with underlying Weibull or gamma distributions. The paper will illustrate the main results achieved during the above research program and will discuss, via some sample cases, the performances of the numerical algorithms which on the appearance of stress corrosion cracking in the tubes of Steam Generators of PWR French power plants. (authors)

Numerical methods for image registration

CERN Document Server

Modersitzki, Jan

2003-01-01

Based on the author's lecture notes and research, this well-illustrated and comprehensive text is one of the first to provide an introduction to image registration with particular emphasis on numerical methods in medical imaging. Ideal for researchers in industry and academia, it is also a suitable study guide for graduate mathematicians, computer scientists, engineers, medical physicists, and radiologists.Image registration is utilised whenever information obtained from different viewpoints needs to be combined or compared and unwanted distortion needs to be eliminated. For example, CCTV imag

INTRA/Mod3.2. Manual and Code Description. Volume I - Physical Modelling

International Nuclear Information System (INIS)

Andersson, Jenny; Edlund, O.; Hermann, J.; Johansson, Lise-Lotte

1999-01-01

The INTRA Manual consists of two volumes. Volume I of the manual is a thorough description of the code INTRA, the Physical modelling of INTRA and the ruling numerical methods and volume II, the User's Manual is an input description. This document, the Physical modelling of INTRA, contains code characteristics, integration methods and applications

INTRA/Mod3.2. Manual and Code Description. Volume I - Physical Modelling

Energy Technology Data Exchange (ETDEWEB)

Andersson, Jenny; Edlund, O; Hermann, J; Johansson, Lise-Lotte

1999-01-01

The INTRA Manual consists of two volumes. Volume I of the manual is a thorough description of the code INTRA, the Physical modelling of INTRA and the ruling numerical methods and volume II, the User`s Manual is an input description. This document, the Physical modelling of INTRA, contains code characteristics, integration methods and applications

An integrated algorithm for hypersonic fluid-thermal-structural numerical simulation

Science.gov (United States)

Li, Jia-Wei; Wang, Jiang-Feng

2018-05-01

In this paper, a fluid-structural-thermal integrated method is presented based on finite volume method. A unified integral equations system is developed as the control equations for physical process of aero-heating and structural heat transfer. The whole physical field is discretized by using an up-wind finite volume method. To demonstrate its capability, the numerical simulation of Mach 6.47 flow over stainless steel cylinder shows a good agreement with measured values, and this method dynamically simulates the objective physical processes. Thus, the integrated algorithm proves to be efficient and reliable.

Monotone numerical methods for finite-state mean-field games

KAUST Repository

Gomes, Diogo A.; Saude, Joao

2017-01-01

Here, we develop numerical methods for finite-state mean-field games (MFGs) that satisfy a monotonicity condition. MFGs are determined by a system of differential equations with initial and terminal boundary conditions. These non-standard conditions are the main difficulty in the numerical approximation of solutions. Using the monotonicity condition, we build a flow that is a contraction and whose fixed points solve the MFG, both for stationary and time-dependent problems. We illustrate our methods in a MFG modeling the paradigm-shift problem.

Monotone numerical methods for finite-state mean-field games

KAUST Repository

Gomes, Diogo A.

2017-04-29

Here, we develop numerical methods for finite-state mean-field games (MFGs) that satisfy a monotonicity condition. MFGs are determined by a system of differential equations with initial and terminal boundary conditions. These non-standard conditions are the main difficulty in the numerical approximation of solutions. Using the monotonicity condition, we build a flow that is a contraction and whose fixed points solve the MFG, both for stationary and time-dependent problems. We illustrate our methods in a MFG modeling the paradigm-shift problem.

Developing Teaching Material Software Assisted for Numerical Methods

Science.gov (United States)

Handayani, A. D.; Herman, T.; Fatimah, S.

2017-09-01

The NCTM vision shows the importance of two things in school mathematics, which is knowing the mathematics of the 21st century and the need to continue to improve mathematics education to answer the challenges of a changing world. One of the competencies associated with the great challenges of the 21st century is the use of help and tools (including IT), such as: knowing the existence of various tools for mathematical activity. One of the significant challenges in mathematical learning is how to teach students about abstract concepts. In this case, technology in the form of mathematics learning software can be used more widely to embed the abstract concept in mathematics. In mathematics learning, the use of mathematical software can make high level math activity become easier accepted by student. Technology can strengthen student learning by delivering numerical, graphic, and symbolic content without spending the time to calculate complex computing problems manually. The purpose of this research is to design and develop teaching materials software assisted for numerical method. The process of developing the teaching material starts from the defining step, the process of designing the learning material developed based on information obtained from the step of early analysis, learners, materials, tasks that support then done the design step or design, then the last step is the development step. The development of teaching materials software assisted for numerical methods is valid in content. While validator assessment for teaching material in numerical methods is good and can be used with little revision.

Introduction to numerical methods for time dependent differential equations

CERN Document Server

Kreiss, Heinz-Otto

2014-01-01

Introduces both the fundamentals of time dependent differential equations and their numerical solutions Introduction to Numerical Methods for Time Dependent Differential Equations delves into the underlying mathematical theory needed to solve time dependent differential equations numerically. Written as a self-contained introduction, the book is divided into two parts to emphasize both ordinary differential equations (ODEs) and partial differential equations (PDEs). Beginning with ODEs and their approximations, the authors provide a crucial presentation of fundamental notions, such as the t

TEMPEST: A three-dimensional time-dependent computer program for hydrothermal analysis: Volume 1, Numerical methods and input instructions

International Nuclear Information System (INIS)

Trent, D.S.; Eyler, L.L.; Budden, M.J.

1983-09-01

This document describes the numerical methods, current capabilities, and the use of the TEMPEST (Version L, MOD 2) computer program. TEMPEST is a transient, three-dimensional, hydrothermal computer program that is designed to analyze a broad range of coupled fluid dynamic and heat transfer systems of particular interest to the Fast Breeder Reactor thermal-hydraulic design community. The full three-dimensional, time-dependent equations of motion, continuity, and heat transport are solved for either laminar or turbulent fluid flow, including heat diffusion and generation in both solid and liquid materials. 10 refs., 22 figs., 2 tabs

Geothermal-Related Thermo-Elastic Fracture Analysis by Numerical Manifold Method

Directory of Open Access Journals (Sweden)

Jun He

2018-05-01

Full Text Available One significant factor influencing geothermal energy exploitation is the variation of the mechanical properties of rock in high temperature environments. Since rock is typically a heterogeneous granular material, thermal fracturing frequently occurs in the rock when the ambient temperature changes, which can greatly influence the geothermal energy exploitation. A numerical method based on the numerical manifold method (NMM is developed in this study to simulate the thermo-elastic fracturing of rocklike granular materials. The Voronoi tessellation is incorporated into the pre-processor of NMM to represent the grain structure. A contact-based heat transfer model is developed to reflect heat interaction among grains. Based on the model, the transient thermal conduction algorithm for granular materials is established. To simulate the cohesion effects among grains and the fracturing process between grains, a damage-based contact fracture model is developed to improve the contact algorithm of NMM. In the developed numerical method, the heat interaction among grains as well as the heat transfer inside each solid grain are both simulated. Additionally, as damage evolution and fracturing at grain interfaces are also considered, the developed numerical method is applicable to simulate the geothermal-related thermal fracturing process.

Rigid inclusions-Comparison between analytical and numerical methods

International Nuclear Information System (INIS)

Gomez Perez, R.; Melentijevic, S.

2014-01-01

This paper compares different analytical methods for analysis of rigid inclusions with finite element modeling. First of all, the load transfer in the distribution layer is analyzed for its different thicknesses and different inclusion grids to define the range between results obtained by analytical and numerical methods. The interaction between the soft soil and the inclusion in the estimation of settlements is studied as well. Considering different stiffness of the soft soil, settlements obtained analytical and numerically are compared. The influence of the soft soil modulus of elasticity on the neutral point depth was also performed by finite elements. This depth has a great importance for the definition of the total length of rigid inclusion. (Author)

Efficient numerical methods for fluid- and electrodynamics on massively parallel systems

Energy Technology Data Exchange (ETDEWEB)

Zudrop, Jens

2016-07-01

In the last decade, computer technology has evolved rapidly. Modern high performance computing systems offer a tremendous amount of computing power in the range of a few peta floating point operations per second. In contrast, numerical software development is much slower and most existing simulation codes cannot exploit the full computing power of these systems. Partially, this is due to the numerical methods themselves and partially it is related to bottlenecks within the parallelization concept and its data structures. The goal of the thesis is the development of numerical algorithms and corresponding data structures to remedy both kinds of parallelization bottlenecks. The approach is based on a co-design of the numerical schemes (including numerical analysis) and their realizations in algorithms and software. Various kinds of applications, from multicomponent flows (Lattice Boltzmann Method) to electrodynamics (Discontinuous Galerkin Method) to embedded geometries (Octree), are considered and efficiency of the developed approaches is demonstrated for large scale simulations.

Direct Calculation of Permeability by High-Accurate Finite Difference and Numerical Integration Methods

KAUST Repository

Wang, Yi

2016-07-21

Velocity of fluid flow in underground porous media is 6~12 orders of magnitudes lower than that in pipelines. If numerical errors are not carefully controlled in this kind of simulations, high distortion of the final results may occur [1-4]. To fit the high accuracy demands of fluid flow simulations in porous media, traditional finite difference methods and numerical integration methods are discussed and corresponding high-accurate methods are developed. When applied to the direct calculation of full-tensor permeability for underground flow, the high-accurate finite difference method is confirmed to have numerical error as low as 10-5% while the high-accurate numerical integration method has numerical error around 0%. Thus, the approach combining the high-accurate finite difference and numerical integration methods is a reliable way to efficiently determine the characteristics of general full-tensor permeability such as maximum and minimum permeability components, principal direction and anisotropic ratio. Copyright © Global-Science Press 2016.

Implementation and assessment of high-resolution numerical methods in TRACE

Energy Technology Data Exchange (ETDEWEB)

Wang, Dean, E-mail: wangda@ornl.gov [Oak Ridge National Laboratory, 1 Bethel Valley RD 6167, Oak Ridge, TN 37831 (United States); Mahaffy, John H.; Staudenmeier, Joseph; Thurston, Carl G. [U.S. Nuclear Regulatory Commission, Washington, DC 20555 (United States)

2013-10-15

Highlights: • Study and implement high-resolution numerical methods for two-phase flow. • They can achieve better numerical accuracy than the 1st-order upwind scheme. • They are of great numerical robustness and efficiency. • Great application for BWR stability analysis and boron injection. -- Abstract: The 1st-order upwind differencing numerical scheme is widely employed to discretize the convective terms of the two-phase flow transport equations in reactor systems analysis codes such as TRACE and RELAP. While very robust and efficient, 1st-order upwinding leads to excessive numerical diffusion. Standard 2nd-order numerical methods (e.g., Lax–Wendroff and Beam–Warming) can effectively reduce numerical diffusion but often produce spurious oscillations for steep gradients. To overcome the difficulties with the standard higher-order schemes, high-resolution schemes such as nonlinear flux limiters have been developed and successfully applied in numerical simulation of fluid-flow problems in recent years. The present work contains a detailed study on the implementation and assessment of six nonlinear flux limiters in TRACE. These flux limiters selected are MUSCL, Van Leer (VL), OSPRE, Van Albada (VA), ENO, and Van Albada 2 (VA2). The assessment is focused on numerical stability, convergence, and accuracy of the flux limiters and their applicability for boiling water reactor (BWR) stability analysis. It is found that VA and MUSCL work best among of the six flux limiters. Both of them not only have better numerical accuracy than the 1st-order upwind scheme but also preserve great robustness and efficiency.

Implementation and assessment of high-resolution numerical methods in TRACE

International Nuclear Information System (INIS)

Wang, Dean; Mahaffy, John H.; Staudenmeier, Joseph; Thurston, Carl G.

2013-01-01

Highlights: • Study and implement high-resolution numerical methods for two-phase flow. • They can achieve better numerical accuracy than the 1st-order upwind scheme. • They are of great numerical robustness and efficiency. • Great application for BWR stability analysis and boron injection. -- Abstract: The 1st-order upwind differencing numerical scheme is widely employed to discretize the convective terms of the two-phase flow transport equations in reactor systems analysis codes such as TRACE and RELAP. While very robust and efficient, 1st-order upwinding leads to excessive numerical diffusion. Standard 2nd-order numerical methods (e.g., Lax–Wendroff and Beam–Warming) can effectively reduce numerical diffusion but often produce spurious oscillations for steep gradients. To overcome the difficulties with the standard higher-order schemes, high-resolution schemes such as nonlinear flux limiters have been developed and successfully applied in numerical simulation of fluid-flow problems in recent years. The present work contains a detailed study on the implementation and assessment of six nonlinear flux limiters in TRACE. These flux limiters selected are MUSCL, Van Leer (VL), OSPRE, Van Albada (VA), ENO, and Van Albada 2 (VA2). The assessment is focused on numerical stability, convergence, and accuracy of the flux limiters and their applicability for boiling water reactor (BWR) stability analysis. It is found that VA and MUSCL work best among of the six flux limiters. Both of them not only have better numerical accuracy than the 1st-order upwind scheme but also preserve great robustness and efficiency

Numerical study of turbulent heat transfer from confined impinging jets using a pseudo-compressibility method

Energy Technology Data Exchange (ETDEWEB)

Rahman, M.; Rautaheimo, P.; Siikonen, T.

1997-12-31

A numerical investigation is carried out to predict the turbulent fluid flow and heat transfer characteristics of two-dimensional single and three impinging slot jets. Two low-Reynolds-number {kappa}-{epsilon} models, namely the classical model of Chien and the explicit algebraic stress model of Gatski and Speziale, are considered in the simulation. A cell-centered finite-volume scheme combined with an artificial compressibility approach is employed to solve the flow equations, using a diagonally dominant alternating direction implicit (DDADI) time integration method. A fully upwinded second order spatial differencing is adopted to approximate the convective terms. Roe`s damping term is used to calculate the flux on the cell face. A multigrid method is utilized for the acceleration of convergence. On average, the heat transfer coefficients predicted by both models show good agreement with the experimental results. (orig.) 17 refs.

Numerical methods in nuclear engineering. Part 1

International Nuclear Information System (INIS)

Phillips, G.J.

1983-08-01

These proceedings, published in two parts contain the full text of 56 papers and summaries of six papers presented at the conference. They cover the use of numerical methods in thermal hydraulics, reactor physics, neutron diffusion, subchannel analysis, risk assessment, transport theory, and fuel behaviour

Numerical solution of the Navier-Stokes equations by discontinuous Galerkin method

Science.gov (United States)

Krasnov, M. M.; Kuchugov, P. A.; E Ladonkina, M.; E Lutsky, A.; Tishkin, V. F.

2017-02-01

Detailed unstructured grids and numerical methods of high accuracy are frequently used in the numerical simulation of gasdynamic flows in areas with complex geometry. Galerkin method with discontinuous basis functions or Discontinuous Galerkin Method (DGM) works well in dealing with such problems. This approach offers a number of advantages inherent to both finite-element and finite-difference approximations. Moreover, the present paper shows that DGM schemes can be viewed as Godunov method extension to piecewise-polynomial functions. As is known, DGM involves significant computational complexity, and this brings up the question of ensuring the most effective use of all the computational capacity available. In order to speed up the calculations, operator programming method has been applied while creating the computational module. This approach makes possible compact encoding of mathematical formulas and facilitates the porting of programs to parallel architectures, such as NVidia CUDA and Intel Xeon Phi. With the software package, based on DGM, numerical simulations of supersonic flow past solid bodies has been carried out. The numerical results are in good agreement with the experimental ones.

Spectral methods in numerical plasma simulation

International Nuclear Information System (INIS)

Coutsias, E.A.; Hansen, F.R.; Huld, T.; Knorr, G.; Lynov, J.P.

1989-01-01

An introduction is given to the use of spectral methods in numerical plasma simulation. As examples of the use of spectral methods, solutions to the two-dimensional Euler equations in both a simple, doubly periodic region, and on an annulus will be shown. In the first case, the solution is expanded in a two-dimensional Fourier series, while a Chebyshev-Fourier expansion is employed in the second case. A new, efficient algorithm for the solution of Poisson's equation on an annulus is introduced. Problems connected to aliasing and to short wavelength noise generated by gradient steepening are discussed. (orig.)

Numerical methods to solve the two-dimensional heat conduction equation

International Nuclear Information System (INIS)

Santos, R.S. dos.

1981-09-01

A class of numerical methods, called 'Hopscotch Algorithms', was used to solve the heat conduction equation in cylindrical geometry. Using a time dependent heat source, the temperature versus time behaviour of cylindric rod was analysed. Numerical simulation was used to study the stability and the convergence of each different method. Another test had the temperature specified on the outer surface as boundary condition. The various Hopscotch methods analysed exhibit differing degrees of accuracy, few of them being so accurate as the ADE method, but requiring more computational operations than the later, were observed. Finally, compared with the so called ODD-EVEN method, two other Hopscotch methods, are more time consuming. (Author) [pt

Numerical study of Taylor bubbles with adaptive unstructured meshes

Science.gov (United States)

Xie, Zhihua; Pavlidis, Dimitrios; Percival, James; Pain, Chris; Matar, Omar; Hasan, Abbas; Azzopardi, Barry

2014-11-01

The Taylor bubble is a single long bubble which nearly fills the entire cross section of a liquid-filled circular tube. This type of bubble flow regime often occurs in gas-liquid slug flows in many industrial applications, including oil-and-gas production, chemical and nuclear reactors, and heat exchangers. The objective of this study is to investigate the fluid dynamics of Taylor bubbles rising in a vertical pipe filled with oils of extremely high viscosity (mimicking the ``heavy oils'' found in the oil-and-gas industry). A modelling and simulation framework is presented here which can modify and adapt anisotropic unstructured meshes to better represent the underlying physics of bubble rise and reduce the computational effort without sacrificing accuracy. The numerical framework consists of a mixed control-volume and finite-element formulation, a ``volume of fluid''-type method for the interface capturing based on a compressive control volume advection method, and a force-balanced algorithm for the surface tension implementation. Numerical examples of some benchmark tests and the dynamics of Taylor bubbles are presented to show the capability of this method. EPSRC Programme Grant, MEMPHIS, EP/K0039761/1.

Geothermal-Related Thermo-Elastic Fracture Analysis by Numerical Manifold Method

OpenAIRE

Jun He; Quansheng Liu; Zhijun Wu; Yalong Jiang

2018-01-01

One significant factor influencing geothermal energy exploitation is the variation of the mechanical properties of rock in high temperature environments. Since rock is typically a heterogeneous granular material, thermal fracturing frequently occurs in the rock when the ambient temperature changes, which can greatly influence the geothermal energy exploitation. A numerical method based on the numerical manifold method (NMM) is developed in this study to simulate the thermo-elastic fracturing ...

A hybrid numerical method for orbit correction

International Nuclear Information System (INIS)

White, G.; Himel, T.; Shoaee, H.

1997-09-01

The authors describe a simple hybrid numerical method for beam orbit correction in particle accelerators. The method overcomes both degeneracy in the linear system being solved and respects boundaries on the solution. It uses the Singular Value Decomposition (SVD) to find and remove the null-space in the system, followed by a bounded Linear Least Squares analysis of the remaining recast problem. It was developed for correcting orbit and dispersion in the B-factory rings

Conservative numerical methods for solitary wave interactions

Energy Technology Data Exchange (ETDEWEB)

Duran, A; Lopez-Marcos, M A [Departamento de Matematica Aplicada y Computacion, Facultad de Ciencias, Universidad de Valladolid, Paseo del Prado de la Magdalena s/n, 47005 Valladolid (Spain)

2003-07-18

The purpose of this paper is to show the advantages that represent the use of numerical methods that preserve invariant quantities in the study of solitary wave interactions for the regularized long wave equation. It is shown that the so-called conservative methods are more appropriate to study the phenomenon and provide a dynamic point of view that allows us to estimate the changes in the parameters of the solitary waves after the collision.

FORECASTING PILE SETTLEMENT ON CLAYSTONE USING NUMERICAL AND ANALYTICAL METHODS

Directory of Open Access Journals (Sweden)

Ponomarev Andrey Budimirovich

2016-06-01

Full Text Available In the article the problem of designing pile foundations on claystones is reviewed. The purpose of this paper is comparative analysis of the analytical and numerical methods for forecasting the settlement of piles on claystones. The following tasks were solved during the study: 1 The existing researches of pile settlement are analyzed; 2 The characteristics of experimental studies and the parameters for numerical modeling are presented, methods of field research of single piles’ operation are described; 3 Calculation of single pile settlement is performed using numerical methods in the software package Plaxis 2D and analytical method according to the requirements SP 24.13330.2011; 4 Experimental data is compared with the results of analytical and numerical calculations; 5 Basing on these results recommendations for forecasting pile settlement on claystone are presented. Much attention is paid to the calculation of pile settlement considering the impacted areas in ground space beside pile and the comparison with the results of field experiments. Basing on the obtained results, for the prediction of settlement of single pile on claystone the authors recommend using the analytical method considered in SP 24.13330.2011 with account for the impacted areas in ground space beside driven pile. In the case of forecasting the settlement of single pile on claystone by numerical methods in Plaxis 2D the authors recommend using the Hardening Soil model considering the impacted areas in ground space beside the driven pile. The analyses of the results and calculations are presented for examination and verification; therefore it is necessary to continue the research work of deep foundation at another experimental sites to improve the reliability of the calculation of pile foundation settlement. The work is of great interest for geotechnical engineers engaged in research, design and construction of pile foundations.

Numerical dissipation and dispersion of the homogenenous and complete flux schemes

NARCIS (Netherlands)

Thije Boonkkamp, ten J.H.M.; Anthonissen, M.J.H.

2014-01-01

We analyse numerical dissipation and dispersion of the homogeneous ¿ux (HF) and complete ¿ux (CF) schemes, ¿nite volume methods introduced in [1]. To that purpose we derive the modi¿ed equation of both schemes. We show that the HF scheme suffers from numerical diffusion for dominant advection, which

Application of numerical analysis methods to thermoluminescence dosimetry

International Nuclear Information System (INIS)

Gomez Ros, J. M.; Delgado, A.

1989-01-01

This report presents the application of numerical methods to thermoluminescence dosimetry (TLD), showing the advantages obtained over conventional evaluation systems. Different configurations of the analysis method are presented to operate in specific dosimetric applications of TLD, such as environmental monitoring and mailed dosimetry systems for quality assurance in radiotherapy facilities. (Author) 10 refs

Numerical method for three dimensional steady-state two-phase flow calculations

International Nuclear Information System (INIS)

Raymond, P.; Toumi, I.

1992-01-01

This paper presents the numerical scheme which was developed for the FLICA-4 computer code to calculate three dimensional steady state two phase flows. This computer code is devoted to steady state and transient thermal hydraulics analysis of nuclear reactor cores 1,3 . The first section briefly describes the FLICA-4 flow modelling. Then in order to introduce the numerical method for steady state computations, some details are given about the implicit numerical scheme based upon an approximate Riemann solver which was developed for calculation of flow transients. The third section deals with the numerical method for steady state computations, which is derived from this previous general scheme and its optimization. We give some numerical results for steady state calculations and comparisons on required CPU time and memory for various meshing and linear system solvers

A numerical method for investigating crystal settling in convecting magma chambers

Science.gov (United States)

Verhoeven, J.; Schmalzl, J.

2009-12-01

Magma chambers can be considered as thermochemically driven convection systems. We present a new numerical method that describes the movement of crystallized minerals in terms of active spherical particles in a convecting magma that is represented by an infinite Prandtl number fluid. The main part focuses on the results we obtained. A finite volume thermochemical convection model for two and three dimensions and a discrete element method, which is used to model granular material, are combined. The new model is validated with floating experiments using particles of different densities and an investigation of single and multiparticle settling velocities. The resulting velocities are compared with theoretical predictions by Stokes's law and a hindered settling function for the multiparticle system. Two fundamental convection regimes are identified in the parameter space that is spanned by the Rayleigh number and the chemical Rayleigh number, which is a measure for the density of the particles. We define the T regime that is dominated by thermal convection. Here the thermal driving force is strong enough to keep all particles in suspension. As the particles get denser, they start settling to the ground, which results in a C regime. The C regime is characterized by the existence of a sediment layer with particle-rich material and a suspension layer with few particles. It is shown that the presence of particles can reduce the vigor of thermal convection. In the frame of a parameter study we discuss the change between the regimes that is systematically investigated. We show that the so-called TC transition fits a power law. Furthermore, we investigate the settling behavior of the particles in vigorous thermal convection, which can be linked to crystal settling in magma chambers. We develop an analytical settling law that describes the number of settled particles against time and show that the results fit the observations from numerical and laboratory experiments.

Application of volume of fluid method for simulation of a droplet impacting a fiber

Directory of Open Access Journals (Sweden)

M. Khalili

2016-06-01

Full Text Available In the present work, impact of a Newtonian drop on horizontal thin fibers with circular cross section is simulated in 2D views. The numerical simulations of the phenomena are carried out using volume of fluid (VOF method for tracking the free surface motion. Impacting of a Newtonian droplet on a circular thin fiber (350μm radius investigated numerically. The main focus of this simulation is to acquire threshold radius and velocity of a drop which is entirely captured by the fiber. The model agrees well with the experiments and demonstrates the threshold radius decreased generally with the increase of impact velocity. In other words, for velocity larger than threshold velocity of capture perhaps only a small portion of fluid is stuck on the solid and the rest of the drop is ejected for impact velocity smaller than critical velocity the drop is totally captured. This threshold velocity has been determined when the impact is centered.

Numerical study of cesium effects on negative ion production in volume sources

Energy Technology Data Exchange (ETDEWEB)

Fukumasa, Osamu; Niitani, Eiji [Yamaguchi Univ., Ube (Japan). Faculty of Engineering

1997-02-01

Effects of cesium vapor injection of H{sup -} production in a tandem negative ion source are studied numerically as a function of plasma parameters. Model calculation is done by solving a set of particle balance equations in a steady-state hydrogen discharge plasmas. Here, the results which focus on gas pressure and electron temperature dependences of H{sup -} volume production are presented and discussed. With including H{sup -} surface production processes caused by both H atoms and positive hydrogen ions, enhancement of H{sup -} production and pressure dependence of H{sup -} production observed experimentally are well reproduced in the model. To enhance H{sup -} production, however, so-called electron cooling is not so effective if plasma parameters are initially optimized with the use of magnetic filter. (author)

Prediction of the Individual Wave Overtopping Volumes of a Wave Energy Converter using Experimental Testing and First Numerical Model Results

DEFF Research Database (Denmark)

Victor, L.; Troch, P.; Kofoed, Jens Peter

2009-01-01

For overtopping wave energy converters (WECs) a more efficient energy conversion can be achieved when the volumes of water, wave by wave, that enter their reservoir are known and can be predicted. A numerical tool is being developed using a commercial CFD-solver to study and optimize...... nearshore 2Dstructure. First numerical model results are given for a specific test with regular waves, and are compared with the corresponding experimental results in this paper....

A numerical method for solving singular De`s

Energy Technology Data Exchange (ETDEWEB)

Mahaver, W.T.

1996-12-31

A numerical method is developed for solving singular differential equations using steepest descent based on weighted Sobolev gradients. The method is demonstrated on a variety of first and second order problems, including linear constrained, unconstrained, and partially constrained first order problems, a nonlinear first order problem with irregular singularity, and two second order variational problems.

Preface of "The Second Symposium on Border Zones Between Experimental and Numerical Application Including Solution Approaches By Extensions of Standard Numerical Methods"

Science.gov (United States)

Ortleb, Sigrun; Seidel, Christian

2017-07-01

In this second symposium at the limits of experimental and numerical methods, recent research is presented on practically relevant problems. Presentations discuss experimental investigation as well as numerical methods with a strong focus on application. In addition, problems are identified which require a hybrid experimental-numerical approach. Topics include fast explicit diffusion applied to a geothermal energy storage tank, noise in experimental measurements of electrical quantities, thermal fluid structure interaction, tensegrity structures, experimental and numerical methods for Chladni figures, optimized construction of hydroelectric power stations, experimental and numerical limits in the investigation of rain-wind induced vibrations as well as the application of exponential integrators in a domain-based IMEX setting.

New simple method for fast and accurate measurement of volumes

International Nuclear Information System (INIS)

Frattolillo, Antonio

2006-01-01

A new simple method is presented, which allows us to measure in just a few minutes but with reasonable accuracy (less than 1%) the volume confined inside a generic enclosure, regardless of the complexity of its shape. The technique proposed also allows us to measure the volume of any portion of a complex manifold, including, for instance, pipes and pipe fittings, valves, gauge heads, and so on, without disassembling the manifold at all. To this purpose an airtight variable volume is used, whose volume adjustment can be precisely measured; it has an overall capacity larger than that of the unknown volume. Such a variable volume is initially filled with a suitable test gas (for instance, air) at a known pressure, as carefully measured by means of a high precision capacitive gauge. By opening a valve, the test gas is allowed to expand into the previously evacuated unknown volume. A feedback control loop reacts to the resulting finite pressure drop, thus contracting the variable volume until the pressure exactly retrieves its initial value. The overall reduction of the variable volume achieved at the end of this process gives a direct measurement of the unknown volume, and definitively gets rid of the problem of dead spaces. The method proposed actually does not require the test gas to be rigorously held at a constant temperature, thus resulting in a huge simplification as compared to complex arrangements commonly used in metrology (gas expansion method), which can grant extremely accurate measurement but requires rather expensive equipments and results in time consuming methods, being therefore impractical in most applications. A simple theoretical analysis of the thermodynamic cycle and the results of experimental tests are described, which demonstrate that, in spite of its simplicity, the method provides a measurement accuracy within 0.5%. The system requires just a few minutes to complete a single measurement, and is ready immediately at the end of the process. The

Numerical method for two phase flow with a unstable interface

International Nuclear Information System (INIS)

Glimm, J.; Marchesin, D.; McBryan, O.

1981-01-01

The random choice method is used to compute the oil-water interface for two dimensional porous media equations. The equations used are a pair of coupled equations; the (elliptic) pressure equation and the (hyperbolic) saturation equation. The equations do not include the dispersive capillary pressure term and the computation does not introduce numerical diffusion. The method resolves saturation discontinuities sharply. The main conclusion of this paper is that the random choice is a correct numerical procedure for this problem even in the highly fingered case. Two methods of inducing fingers are considered: deterministically, through choice of Cauchy data and heterogeneity, through maximizing the randomness of the random choice method

An implicit second order numerical method for two-fluid models

International Nuclear Information System (INIS)

Toumi, I.

1995-01-01

We present an implicit upwind numerical method for a six equation two-fluid model based on a linearized Riemann solver. The construction of this approximate Riemann solver uses an extension of Roe's scheme. Extension to second order accurate method is achieved using a piecewise linear approximation of the solution and a slope limiter method. For advancing in time, a linearized implicit integrating step is used. In practice this new numerical method has proved to be stable and capable of generating accurate non-oscillating solutions for two-phase flow calculations. The scheme was applied both to shock tube problems and to standard tests for two-fluid codes. (author)

NUMERICAL WITHOUT ITERATION METHOD OF MODELING OF ELECTROMECHANICAL PROCESSES IN ASYNCHRONOUS ENGINES

Directory of Open Access Journals (Sweden)

D. G. Patalakh

2018-02-01

Full Text Available Purpose. Development of calculation of electromagnetic and electromechanic transients is in asynchronous engines without iterations. Methodology. Numeral methods of integration of usual differential equations, programming. Findings. As the system of equations, describing the dynamics of asynchronous engine, contents the products of rotor and stator currents and product of rotation frequency of rotor and currents, so this system is nonlinear one. The numeral solution of nonlinear differential equations supposes an iteration process on every step of integration. Time-continuing and badly converging iteration process may be the reason of calculation slowing. The improvement of numeral method by the way of an iteration process removing is offered. As result the modeling time is reduced. The improved numeral method is applied for integration of differential equations, describing the dynamics of asynchronous engine. Originality. The improvement of numeral method allowing to execute numeral integrations of differential equations containing product of functions is offered, that allows to avoid an iteration process on every step of integration and shorten modeling time. Practical value. On the basis of the offered methodology the universal program of modeling of electromechanics processes in asynchronous engines could be developed as taking advantage on fast-acting.

3rd GAMM/IFIP-Workshop on “Stochastic Optimization: Numerical Methods and Technical Applications” held at the Federal Armed Forces University Munich

CERN Document Server

Kall, Peter

1998-01-01

Optimization problems arising in practice usually contain several random parameters. Hence, in order to obtain optimal solutions being robust with respect to random parameter variations, the mostly available statistical information about the random parameters should be considered already at the planning phase. The original problem with random parameters must be replaced by an appropriate deterministic substitute problem, and efficient numerical solution or approximation techniques have to be developed for those problems. This proceedings volume contains a selection of papers on modelling techniques, approximation methods, numerical solution procedures for stochastic optimization problems and applications to the reliability-based optimization of concrete technical or economic systems.

Numerical methods for differential equations and applications

International Nuclear Information System (INIS)

Ixaru, L.G.

1984-01-01

This book is addressed to persons who, without being professionals in applied mathematics, are often faced with the problem of numerically solving differential equations. In each of the first three chapters a definite class of methods is discussed for the solution of the initial value problem for ordinary differential equations: multistep methods; one-step methods; and piecewise perturbation methods. The fourth chapter is mainly focussed on the boundary value problems for linear second-order equations, with a section devoted to the Schroedinger equation. In the fifth chapter the eigenvalue problem for the radial Schroedinger equation is solved in several ways, with computer programs included. (Auth.)

A numerical test of the collective coordinate method

International Nuclear Information System (INIS)

Dobrowolski, T.; Tatrocki, P.

2008-01-01

The purpose of this Letter is to compare the dynamics of the kink interacting with the imperfection which follows from the collective coordinate method with the numerical results obtained on the ground of the field theoretical model. We showed that for weekly interacting kinks the collective coordinate method works similarly well for low and extremely large speeds

Free surface modelling with two-fluid model and reduced numerical diffusion of the interface

International Nuclear Information System (INIS)

Strubelj, Luka; Tiselj, Izrok

2008-01-01

Full text of publication follows: The free surface flows are successfully modelled with one of existing free surface models, such as: level set method, volume of fluid method (with/without surface reconstruction), front tracking, two-fluid model (two momentum equations) with modified interphase force and others. The main disadvantage of two-fluid model used for simulations of free surface flows is numerical diffusion of the interface, which can be significantly reduced using the method presented in this paper. Several techniques for reduction of numerical diffusion of the interface have been implemented in the volume of fluid model and are based on modified numerical schemes for advection of volume fraction near the interface. The same approach could be used also for two-fluid method, but according to our experience more successful reduction of numerical diffusion of the interface can be achieved with conservative level set method. Within the conservative level set method, continuity equation for volume fraction is solved and after that the numerical diffusion of the interface is reduced in such a way that the thickness of the interface is kept constant during the simulation. Reduction of the interface diffusion can be also called interface sharpening. In present paper the two-fluid model with interface sharpening is validated on Rayleigh-Taylor instability. Under assumptions of isothermal and incompressible flow of two immiscible fluids, we simulated a system with the fluid of higher density located above the fluid of smaller density in two dimensions. Due to gravity in the system, fluid with higher density moves below the fluid with smaller density. Initial condition is not a flat interface between the fluids, but a sine wave with small amplitude, which develops into a mushroom-like structure. Mushroom-like structure in simulation of Rayleigh-Taylor instability later develops to small droplets as result of numerical dispersion of interface (interface sharpening

RELAP-7 Numerical Stabilization: Entropy Viscosity Method

Energy Technology Data Exchange (ETDEWEB)

R. A. Berry; M. O. Delchini; J. Ragusa

2014-06-01

The RELAP-7 code is the next generation nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL). The code is based on the INL's modern scientific software development framework, MOOSE (Multi-Physics Object Oriented Simulation Environment). The overall design goal of RELAP-7 is to take advantage of the previous thirty years of advancements in computer architecture, software design, numerical integration methods, and physical models. The end result will be a reactor systems analysis capability that retains and improves upon RELAP5's capability and extends the analysis capability for all reactor system simulation scenarios. RELAP-7 utilizes a single phase and a novel seven-equation two-phase flow models as described in the RELAP-7 Theory Manual (INL/EXT-14-31366). The basic equation systems are hyperbolic, which generally require some type of stabilization (or artificial viscosity) to capture nonlinear discontinuities and to suppress advection-caused oscillations. This report documents one of the available options for this stabilization in RELAP-7 -- a new and novel approach known as the entropy viscosity method. Because the code is an ongoing development effort in which the physical sub models, numerics, and coding are evolving, so too must the specific details of the entropy viscosity stabilization method. Here the fundamentals of the method in their current state are presented.

A new method to estimate the atomic volume of ternary intermetallic compounds

International Nuclear Information System (INIS)

Pani, M.; Merlo, F.

2011-01-01

The atomic volume of an A x B y C z ternary intermetallic compound can be calculated starting from volumes of some proper A-B, A-C and B-C binary phases. The three methods by Colinet, Muggianu and Kohler, originally used to estimate thermodynamic quantities, and a new method here proposed, were tested to derive volume data in eight systems containing 91 ternary phases with the known structure. The comparison between experimental and calculated volume values shows the best agreement both for the Kohler method and for the new proposed procedure. -- Graphical abstract: Synopsys: the volume of a ternary intermetallic compound can be calculated starting from volumes of some binary phases, selected by the methods of Colinet, Muggianu, Kohler and a new method proposed here. The so obtained values are compared with the experimental ones for eight ternary systems. Display Omitted Research highlights: → The application of some thermodinamic methods to a crystallochemical problem. → The prevision of the average atomic volume of ternary intermetallic phases. → The proposal of a new procedure to select the proper starting set of binary phases.

Advanced Numerical and Theoretical Methods for Photonic Crystals and Metamaterials

Science.gov (United States)

Felbacq, Didier

2016-11-01

This book provides a set of theoretical and numerical tools useful for the study of wave propagation in metamaterials and photonic crystals. While concentrating on electromagnetic waves, most of the material can be used for acoustic (or quantum) waves. For each presented numerical method, numerical code written in MATLAB® is presented. The codes are limited to 2D problems and can be easily translated in Python or Scilab, and used directly with Octave as well.

A asymptotic numerical method for the steady-state convection diffusion equation

International Nuclear Information System (INIS)

Wu Qiguang

1988-01-01

In this paper, A asymptotic numerical method for the steady-state Convection diffusion equation is proposed, which need not take very fine mesh size in the neighbourhood of the boundary layer. Numerical computation for model problem show that we can obtain the numerical solution in the boundary layer with moderate step size

Numerical methods and computers used in elastohydrodynamic lubrication

Science.gov (United States)

Hamrock, B. J.; Tripp, J. H.

1982-01-01

Some of the methods of obtaining approximate numerical solutions to boundary value problems that arise in elastohydrodynamic lubrication are reviewed. The highlights of four general approaches (direct, inverse, quasi-inverse, and Newton-Raphson) are sketched. Advantages and disadvantages of these approaches are presented along with a flow chart showing some of the details of each. The basic question of numerical stability of the elastohydrodynamic lubrication solutions, especially in the pressure spike region, is considered. Computers used to solve this important class of lubrication problems are briefly described, with emphasis on supercomputers.

Numerical modeling in materials science and engineering

CERN Document Server

Rappaz, Michel; Deville, Michel

2003-01-01

This book introduces the concepts and methodologies related to the modelling of the complex phenomena occurring in materials processing. After a short reminder of conservation laws and constitutive relationships, the authors introduce the main numerical methods: finite differences, finite volumes and finite elements. These techniques are developed in three main chapters of the book that tackle more specific problems: phase transformation, solid mechanics and fluid flow. The two last chapters treat inverse methods to obtain the boundary conditions or the material properties and stochastic methods for microstructural simulation. This book is intended for undergraduate and graduate students in materials science and engineering, mechanical engineering and physics and for engineering professionals or researchers who want to get acquainted with numerical simulation to model and compute materials processing.

A first course in ordinary differential equations analytical and numerical methods

CERN Document Server

Hermann, Martin

2014-01-01

This book presents a modern introduction to analytical and numerical techniques for solving ordinary differential equations (ODEs). Contrary to the traditional format—the theorem-and-proof format—the book is focusing on analytical and numerical methods. The book supplies a variety of problems and examples, ranging from the elementary to the advanced level, to introduce and study the mathematics of ODEs. The analytical part of the book deals with solution techniques for scalar first-order and second-order linear ODEs, and systems of linear ODEs—with a special focus on the Laplace transform, operator techniques and power series solutions. In the numerical part, theoretical and practical aspects of Runge-Kutta methods for solving initial-value problems and shooting methods for linear two-point boundary-value problems are considered. The book is intended as a primary text for courses on the theory of ODEs and numerical treatment of ODEs for advanced undergraduate and early graduate students. It is assumed t...

8th conference on Finite Volumes for Complex Applications

CERN Document Server

Omnes, Pascal

2017-01-01

This first volume of the proceedings of the 8th conference on "Finite Volumes for Complex Applications" (Lille, June 2017) covers various topics including convergence and stability analysis, as well as investigations of these methods from the point of view of compatibility with physical principles. It collects together the focused invited papers comparing advanced numerical methods for Stokes and Navier–Stokes equations on a benchmark, as well as reviewed contributions from internationally leading researchers in the field of analysis of finite volume and related methods, offering a comprehensive overview of the state of the art in the field. The finite volume method in its various forms is a space discretization technique for partial differential equations based on the fundamental physical principle of conservation, and recent decades have brought significant advances in the theoretical understanding of the method. Many finite volume methods preserve further qualitative or asymptotic properties, including m...

Direct numerical methods of mathematical modeling in mechanical structural design

International Nuclear Information System (INIS)

Sahili, Jihad; Verchery, Georges; Ghaddar, Ahmad; Zoaeter, Mohamed

2002-01-01

Full text.Structural design and numerical methods are generally interactive; requiring optimization procedures as the structure is analyzed. This analysis leads to define some mathematical terms, as the stiffness matrix, which are resulting from the modeling and then used in numerical techniques during the dimensioning procedure. These techniques and many others involve the calculation of the generalized inverse of the stiffness matrix, called also the 'compliance matrix'. The aim of this paper is to introduce first, some different existing mathematical procedures, used to calculate the compliance matrix from the stiffness matrix, then apply direct numerical methods to solve the obtained system with the lowest computational time, and to compare the obtained results. The results show a big difference of the computational time between the different procedures

A numerical simulation of wheel spray for simplified vehicle model based on discrete phase method

Directory of Open Access Journals (Sweden)

Xingjun Hu

2015-07-01

Full Text Available Road spray greatly affects vehicle body soiling and driving safety. The study of road spray has attracted increasing attention. In this article, computational fluid dynamics software with widely used finite volume method code was employed to investigate the numerical simulation of spray induced by a simplified wheel model and a modified square-back model proposed by the Motor Industry Research Association. Shear stress transport k-omega turbulence model, discrete phase model, and Eulerian wall-film model were selected. In the simulation process, the phenomenon of breakup and coalescence of drops were considered, and the continuous and discrete phases were treated as two-way coupled in momentum and turbulent motion. The relationship between the vehicle external flow structure and body soiling was also discussed.

Finite volume form factors in the presence of integrable defects

International Nuclear Information System (INIS)

Bajnok, Z.; Buccheri, F.; Hollo, L.; Konczer, J.; Takacs, G.

2014-01-01

We developed the theory of finite volume form factors in the presence of integrable defects. These finite volume form factors are expressed in terms of the infinite volume form factors and the finite volume density of states and incorporate all polynomial corrections in the inverse of the volume. We tested our results, in the defect Lee–Yang model, against numerical data obtained by truncated conformal space approach (TCSA), which we improved by renormalization group methods adopted to the defect case. To perform these checks we determined the infinite volume defect form factors in the Lee–Yang model exactly, including their vacuum expectation values. We used these data to calculate the two point functions, which we compared, at short distance, to defect CFT. We also derived explicit expressions for the exact finite volume one point functions, which we checked numerically. In all of these comparisons excellent agreement was found

Numerical Methods for the Design and Analysis of Photonic Crystal Fibres

DEFF Research Database (Denmark)

Roberts, John

2008-01-01

The numerical methods available for calculating the electromagnetic mode properties of photonic crystal fibres are reviewed. The preferred schemes for analyzing TIR guiding and band gap guiding fibres are contrasted.......The numerical methods available for calculating the electromagnetic mode properties of photonic crystal fibres are reviewed. The preferred schemes for analyzing TIR guiding and band gap guiding fibres are contrasted....

A New Method to Solve Numeric Solution of Nonlinear Dynamic System

Directory of Open Access Journals (Sweden)

Min Hu

2016-01-01

Full Text Available It is well known that the cubic spline function has advantages of simple forms, good convergence, approximation, and second-order smoothness. A particular class of cubic spline function is constructed and an effective method to solve the numerical solution of nonlinear dynamic system is proposed based on the cubic spline function. Compared with existing methods, this method not only has high approximation precision, but also avoids the Runge phenomenon. The error analysis of several methods is given via two numeric examples, which turned out that the proposed method is a much more feasible tool applied to the engineering practice.

NUMERICAL SIMULATION OF SHOCK WAVE REFRACTION ON INCLINED CONTACT DISCONTINUITY

Directory of Open Access Journals (Sweden)

P. V. Bulat

2016-05-01

Full Text Available We consider numerical simulation of shock wave refraction on plane contact discontinuity, separating two gases with different density. Discretization of Euler equations is based on finite volume method and WENO finite difference schemes, implemented on unstructured meshes. Integration over time is performed with the use of the third-order Runge–Kutta stepping procedure. The procedure of identification and classification of gas dynamic discontinuities based on conditions of dynamic consistency and image processing methods is applied to visualize and interpret the results of numerical calculations. The flow structure and its quantitative characteristics are defined. The results of numerical and experimental visualization (shadowgraphs, schlieren images, and interferograms are compared.

On numerical solution of Burgers' equation by homotopy analysis method

International Nuclear Information System (INIS)

Inc, Mustafa

2008-01-01

In this Letter, we present the Homotopy Analysis Method (shortly HAM) for obtaining the numerical solution of the one-dimensional nonlinear Burgers' equation. The initial approximation can be freely chosen with possible unknown constants which can be determined by imposing the boundary and initial conditions. Convergence of the solution and effects for the method is discussed. The comparison of the HAM results with the Homotopy Perturbation Method (HPM) and the results of [E.N. Aksan, Appl. Math. Comput. 174 (2006) 884; S. Kutluay, A. Esen, Int. J. Comput. Math. 81 (2004) 1433; S. Abbasbandy, M.T. Darvishi, Appl. Math. Comput. 163 (2005) 1265] are made. The results reveal that HAM is very simple and effective. The HAM contains the auxiliary parameter h, which provides us with a simple way to adjust and control the convergence region of solution series. The numerical solutions are compared with the known analytical and some numerical solutions

Using high-order polynomial basis in 3-D EM forward modeling based on volume integral equation method

Science.gov (United States)

Kruglyakov, Mikhail; Kuvshinov, Alexey

2018-05-01

3-D interpretation of electromagnetic (EM) data of different origin and scale becomes a common practice worldwide. However, 3-D EM numerical simulations (modeling)—a key part of any 3-D EM data analysis—with realistic levels of complexity, accuracy and spatial detail still remains challenging from the computational point of view. We present a novel, efficient 3-D numerical solver based on a volume integral equation (IE) method. The efficiency is achieved by using a high-order polynomial (HOP) basis instead of the zero-order (piecewise constant) basis that is invoked in all routinely used IE-based solvers. We demonstrate that usage of the HOP basis allows us to decrease substantially the number of unknowns (preserving the same accuracy), with corresponding speed increase and memory saving.

Theoretical and applied aerodynamics and related numerical methods

CERN Document Server

Chattot, J J

2015-01-01

This book covers classical and modern aerodynamics, theories and related numerical methods, for senior and first-year graduate engineering students, including: -The classical potential (incompressible) flow theories for low speed aerodynamics of thin airfoils and high and low aspect ratio wings. - The linearized theories for compressible subsonic and supersonic aerodynamics. - The nonlinear transonic small disturbance potential flow theory, including supercritical wing sections, the extended transonic area rule with lift effect, transonic lifting line and swept or oblique wings to minimize wave drag. Unsteady flow is also briefly discussed. Numerical simulations based on relaxation mixed-finite difference methods are presented and explained. - Boundary layer theory for all Mach number regimes and viscous/inviscid interaction procedures used in practical aerodynamics calculations. There are also four chapters covering special topics, including wind turbines and propellers, airplane design, flow analogies and h...

A numerical method to compute interior transmission eigenvalues

International Nuclear Information System (INIS)

Kleefeld, Andreas

2013-01-01

In this paper the numerical calculation of eigenvalues of the interior transmission problem arising in acoustic scattering for constant contrast in three dimensions is considered. From the computational point of view existing methods are very expensive, and are only able to show the existence of such transmission eigenvalues. Furthermore, they have trouble finding them if two or more eigenvalues are situated closely together. We present a new method based on complex-valued contour integrals and the boundary integral equation method which is able to calculate highly accurate transmission eigenvalues. So far, this is the first paper providing such accurate values for various surfaces different from a sphere in three dimensions. Additionally, the computational cost is even lower than those of existing methods. Furthermore, the algorithm is capable of finding complex-valued eigenvalues for which no numerical results have been reported yet. Until now, the proof of existence of such eigenvalues is still open. Finally, highly accurate eigenvalues of the interior Dirichlet problem are provided and might serve as test cases to check newly derived Faber–Krahn type inequalities for larger transmission eigenvalues that are not yet available. (paper)

EFFECTS OF DIFFERENT NUMERICAL INTERFACE METHODS ON HYDRODYNAMICS INSTABILITY

Energy Technology Data Exchange (ETDEWEB)

FRANCOIS, MARIANNE M. [Los Alamos National Laboratory; DENDY, EDWARD D. [Los Alamos National Laboratory; LOWRIE, ROBERT B. [Los Alamos National Laboratory; LIVESCU, DANIEL [Los Alamos National Laboratory; STEINKAMP, MICHAEL J. [Los Alamos National Laboratory

2007-01-11

The authors compare the effects of different numerical schemes for the advection and material interface treatments on the single-mode Rayleigh-Taylor instability, using the RAGE hydro-code. The interface growth and its surface density (interfacial area) versus time are investigated. The surface density metric shows to be better suited to characterize the difference in the flow, than the conventional interface growth metric. They have found that Van Leer's limiter combined to no interface treatment leads to the largest surface area. Finally, to quantify the difference between the numerical methods they have estimated the numerical viscosity in the linear-regime at different scales.

Conservation properties of numerical integration methods for systems of ordinary differential equations

Science.gov (United States)

Rosenbaum, J. S.

1976-01-01

If a system of ordinary differential equations represents a property conserving system that can be expressed linearly (e.g., conservation of mass), it is then desirable that the numerical integration method used conserve the same quantity. It is shown that both linear multistep methods and Runge-Kutta methods are 'conservative' and that Newton-type methods used to solve the implicit equations preserve the inherent conservation of the numerical method. It is further shown that a method used by several authors is not conservative.

Numerical simulation of spray coalescence in an Eulerian framework: Direct quadrature method of moments and multi-fluid method

International Nuclear Information System (INIS)

Fox, R.O.; Laurent, F.; Massot, M.

2008-01-01

The scope of the present study is Eulerian modeling and simulation of polydisperse liquid sprays undergoing droplet coalescence and evaporation. The fundamental mathematical description is the Williams spray equation governing the joint number density function f(v,u;x,t) of droplet volume and velocity. Eulerian multi-fluid models have already been rigorously derived from this equation in Laurent et al. [F. Laurent, M. Massot, P. Villedieu, Eulerian multi-fluid modeling for the numerical simulation of coalescence in polydisperse dense liquid sprays, J. Comput. Phys. 194 (2004) 505-543]. The first key feature of the paper is the application of direct quadrature method of moments (DQMOM) introduced by Marchisio and Fox [D.L. Marchisio, R.O. Fox, Solution of population balance equations using the direct quadrature method of moments, J. Aerosol Sci. 36 (2005) 43-73] to the Williams spray equation. Both the multi-fluid method and DQMOM yield systems of Eulerian conservation equations with complicated interaction terms representing coalescence. In order to focus on the difficulties associated with treating size-dependent coalescence and to avoid numerical uncertainty issues associated with two-way coupling, only one-way coupling between the droplets and a given gas velocity field is considered. In order to validate and compare these approaches, the chosen configuration is a self-similar 2D axisymmetrical decelerating nozzle with sprays having various size distributions, ranging from smooth ones up to Dirac delta functions. The second key feature of the paper is a thorough comparison of the two approaches for various test-cases to a reference solution obtained through a classical stochastic Lagrangian solver. Both Eulerian models prove to describe adequately spray coalescence and yield a very interesting alternative to the Lagrangian solver. The third key point of the study is a detailed description of the limitations associated with each method, thus giving criteria for

Comparison of the accuracy of three angiographic methods for calculating left ventricular volume measurement

International Nuclear Information System (INIS)

Hu Lin; Cui Wei; Shi Hanwen; Tian Yingping; Wang Weigang; Feng Yanguang; Huang Xueyan; Liu Zhisheng

2003-01-01

Objective: To compare the relative accuracy of three methods measuring left ventricular volume by X-ray ventriculography: single plane area-length method, biplane area-length method, and single-plane Simpson's method. Methods: Left ventricular casts were obtained within 24 hours after death from 12 persons who died from non-cardiac causes. The true left ventricular cast volume was measured by water displacement. The calculated volume of the casts was obtained with 3 angiographic methods, i.e., single-plane area-length method, biplane area-length method, and single-plane Simpson's method. Results: The actual average volume of left ventricular casts was (61.17±26.49) ml. The left ventricular volume was averagely (97.50±35.56) ml with single plane area-length method, (90.51±36.33) ml with biplane area-length method, and (65.00± 23.63) ml with single-plane Simpson's method. The left ventricular volumes calculated with single-plane and biplane area-length method were significantly larger than that the actual volumes (P 0.05). The left ventricular volumes calculated with single-plane and biplane area-length method were significantly larger than those calculated with single-plane Simpson's method (P 0.05). The over-estimation of left ventricular volume by single plane area-length method (36.34±17.98) ml and biplane area-length method (29.34±15.59) ml was more obvious than that calculated by single-plane Simpson's method (3.83±8.48) ml. Linear regression analysis showed that there was close correlations between left ventricular volumes calculated with single plane area-length method, biplane area-length method, Simpson's method and the true volume (all r>0.98). Conclusion: Single-plane Simpson's method is more accurate than single plane area-length method and biplane area-length method for left ventricular volume measurement; however, both the single-plane and biplane area-length methods could be used in clinical practice, especially in those imaging modality

Numerical simulation of compressible two-phase flow using a diffuse interface method

International Nuclear Information System (INIS)

Ansari, M.R.; Daramizadeh, A.

2013-01-01

Highlights: ► Compressible two-phase gas–gas and gas–liquid flows simulation are conducted. ► Interface conditions contain shock wave and cavitations. ► A high-resolution diffuse interface method is investigated. ► The numerical results exhibit very good agreement with experimental results. -- Abstract: In this article, a high-resolution diffuse interface method is investigated for simulation of compressible two-phase gas–gas and gas–liquid flows, both in the presence of shock wave and in flows with strong rarefaction waves similar to cavitations. A Godunov method and HLLC Riemann solver is used for discretization of the Kapila five-equation model and a modified Schmidt equation of state (EOS) is used to simulate the cavitation regions. This method is applied successfully to some one- and two-dimensional compressible two-phase flows with interface conditions that contain shock wave and cavitations. The numerical results obtained in this attempt exhibit very good agreement with experimental results, as well as previous numerical results presented by other researchers based on other numerical methods. In particular, the algorithm can capture the complex flow features of transient shocks, such as the material discontinuities and interfacial instabilities, without any oscillation and additional diffusion. Numerical examples show that the results of the method presented here compare well with other sophisticated modeling methods like adaptive mesh refinement (AMR) and local mesh refinement (LMR) for one- and two-dimensional problems

Simulation of 3D parachute fluid–structure interaction based on nonlinear finite element method and preconditioning finite volume method

Directory of Open Access Journals (Sweden)

Fan Yuxin

2014-12-01

Full Text Available A fluid–structure interaction method combining a nonlinear finite element algorithm with a preconditioning finite volume method is proposed in this paper to simulate parachute transient dynamics. This method uses a three-dimensional membrane–cable fabric model to represent a parachute system at a highly folded configuration. The large shape change during parachute inflation is computed by the nonlinear Newton–Raphson iteration and the linear system equation is solved by the generalized minimal residual (GMRES method. A membrane wrinkling algorithm is also utilized to evaluate the special uniaxial tension state of membrane elements on the parachute canopy. In order to avoid large time expenses during structural nonlinear iteration, the implicit Hilber–Hughes–Taylor (HHT time integration method is employed. For the fluid dynamic simulations, the Roe and HLLC (Harten–Lax–van Leer contact scheme has been modified and extended to compute flow problems at all speeds. The lower–upper symmetric Gauss–Seidel (LU-SGS approximate factorization is applied to accelerate the numerical convergence speed. Finally, the test model of a highly folded C-9 parachute is simulated at a prescribed speed and the results show similar characteristics compared with experimental results and previous literature.

A new electric method for non-invasive continuous monitoring of stroke volume and ventricular volume-time curves

Directory of Open Access Journals (Sweden)

Konings Maurits K

2012-08-01

Full Text Available Abstract Background In this paper a new non-invasive, operator-free, continuous ventricular stroke volume monitoring device (Hemodynamic Cardiac Profiler, HCP is presented, that measures the average stroke volume (SV for each period of 20 seconds, as well as ventricular volume-time curves for each cardiac cycle, using a new electric method (Ventricular Field Recognition with six independent electrode pairs distributed over the frontal thoracic skin. In contrast to existing non-invasive electric methods, our method does not use the algorithms of impedance or bioreactance cardiography. Instead, our method is based on specific 2D spatial patterns on the thoracic skin, representing the distribution, over the thorax, of changes in the applied current field caused by cardiac volume changes during the cardiac cycle. Since total heart volume variation during the cardiac cycle is a poor indicator for ventricular stroke volume, our HCP separates atrial filling effects from ventricular filling effects, and retrieves the volume changes of only the ventricles. Methods ex-vivo experiments on a post-mortem human heart have been performed to measure the effects of increasing the blood volume inside the ventricles in isolation, leaving the atrial volume invariant (which can not be done in-vivo. These effects have been measured as a specific 2D pattern of voltage changes on the thoracic skin. Furthermore, a working prototype of the HCP has been developed that uses these ex-vivo results in an algorithm to decompose voltage changes, that were measured in-vivo by the HCP on the thoracic skin of a human volunteer, into an atrial component and a ventricular component, in almost real-time (with a delay of maximally 39 seconds. The HCP prototype has been tested in-vivo on 7 human volunteers, using G-suit inflation and deflation to provoke stroke volume changes, and LVot Doppler as a reference technique. Results The ex-vivo measurements showed that ventricular filling

Deformation data modeling through numerical models: an efficient method for tracking magma transport

Science.gov (United States)

Charco, M.; Gonzalez, P. J.; Galán del Sastre, P.

2017-12-01

Nowadays, multivariate collected data and robust physical models at volcano observatories are becoming crucial for providing effective volcano monitoring. Nevertheless, the forecast of volcanic eruption is notoriously difficult. Wthin this frame one of the most promising methods to evaluate the volcano hazard is the use of surface ground deformation and in the last decades many developments in the field of deformation modeling has been achieved. In particular, numerical modeling allows realistic media features such as topography and crustal heterogeneities to be included, although it is still very time cosuming to solve the inverse problem for near-real time interpretations. Here, we present a method that can be efficiently used to estimate the location and evolution of magmatic sources base on real-time surface deformation data and Finite Element (FE) models. Generally, the search for the best-fitting magmatic (point) source(s) is conducted for an array of 3-D locations extending below a predefined volume region and the Green functions for all the array components have to be precomputed. We propose a FE model for the pre-computation of Green functions in a mechanically heterogeneous domain which eventually will lead to a better description of the status of the volcanic area. The number of Green functions is reduced here to the number of observational points by using their reciprocity relationship. We present and test this methodology with an optimization method base on a Genetic Algorithm. Following synthetic and sensitivity test to estimate the uncertainty of the model parameters, we apply the tool for magma tracking during 2007 Kilauea volcano intrusion and eruption. We show how data inversion with numerical models can speed up the source parameters estimations for a given volcano showing signs of unrest.

Second GAMM-conference on numerical methods in fluid mechanics

International Nuclear Information System (INIS)

Hirschel, E.H.; Geller, W.

1977-01-01

Proceedings of the Second GAMM-Conference on Numerical Methods in Fluid Mechanics held at the DFVLR, Koeln, October 11 to 13, 1977. The conference was attended by approximately 100 participants from 13 European countries representing quite different fields ranging from Aerodynamics to Nuclear Energy. At the meeting 34 papers were presented, many of them concerned with basic problems in the field. It was well demonstrated that Numerical Methods in Fluid Mechanics do not only serve as means for the computation of flow fields but also as tools in the analysis of fluid mechanical phenomena, a role of large future importance if one considers the complexity especially of three-dimensional flows. (orig./RW) [de

Advanced Topics in Computational Partial Differential Equations: Numerical Methods and Diffpack Programming

International Nuclear Information System (INIS)

Katsaounis, T D

2005-01-01

The scope of this book is to present well known simple and advanced numerical methods for solving partial differential equations (PDEs) and how to implement these methods using the programming environment of the software package Diffpack. A basic background in PDEs and numerical methods is required by the potential reader. Further, a basic knowledge of the finite element method and its implementation in one and two space dimensions is required. The authors claim that no prior knowledge of the package Diffpack is required, which is true, but the reader should be at least familiar with an object oriented programming language like C++ in order to better comprehend the programming environment of Diffpack. Certainly, a prior knowledge or usage of Diffpack would be a great advantage to the reader. The book consists of 15 chapters, each one written by one or more authors. Each chapter is basically divided into two parts: the first part is about mathematical models described by PDEs and numerical methods to solve these models and the second part describes how to implement the numerical methods using the programming environment of Diffpack. Each chapter closes with a list of references on its subject. The first nine chapters cover well known numerical methods for solving the basic types of PDEs. Further, programming techniques on the serial as well as on the parallel implementation of numerical methods are also included in these chapters. The last five chapters are dedicated to applications, modelled by PDEs, in a variety of fields. In summary, the book focuses on the computational and implementational issues involved in solving partial differential equations. The potential reader should have a basic knowledge of PDEs and the finite difference and finite element methods. The examples presented are solved within the programming framework of Diffpack and the reader should have prior experience with the particular software in order to take full advantage of the book. Overall

Uniqueness and numerical methods in inverse obstacle scattering

International Nuclear Information System (INIS)

Kress, Rainer

2007-01-01

The inverse problem we consider in this tutorial is to determine the shape of an obstacle from the knowledge of the far field pattern for scattering of time-harmonic plane waves. In the first part we will concentrate on the issue of uniqueness, i.e., we will investigate under what conditions an obstacle and its boundary condition can be identified from a knowledge of its far field pattern for incident plane waves. We will review some classical and some recent results and draw attention to open problems. In the second part we will survey on numerical methods for solving inverse obstacle scattering problems. Roughly speaking, these methods can be classified into three groups. Iterative methods interpret the inverse obstacle scattering problem as a nonlinear ill-posed operator equation and apply iterative schemes such as regularized Newton methods, Landweber iterations or conjugate gradient methods for its solution. Decomposition methods, in principle, separate the inverse scattering problem into an ill-posed linear problem to reconstruct the scattered wave from its far field and the subsequent determination of the boundary of the scatterer from the boundary condition. Finally, the third group consists of the more recently developed sampling methods. These are based on the numerical evaluation of criteria in terms of indicator functions that decide whether a point lies inside or outside the scatterer. The tutorial will give a survey by describing one or two representatives of each group including a discussion on the various advantages and disadvantages

Monte Carlo Method with Heuristic Adjustment for Irregularly Shaped Food Product Volume Measurement

Directory of Open Access Journals (Sweden)

Joko Siswantoro

2014-01-01

Full Text Available Volume measurement plays an important role in the production and processing of food products. Various methods have been proposed to measure the volume of food products with irregular shapes based on 3D reconstruction. However, 3D reconstruction comes with a high-priced computational cost. Furthermore, some of the volume measurement methods based on 3D reconstruction have a low accuracy. Another method for measuring volume of objects uses Monte Carlo method. Monte Carlo method performs volume measurements using random points. Monte Carlo method only requires information regarding whether random points fall inside or outside an object and does not require a 3D reconstruction. This paper proposes volume measurement using a computer vision system for irregularly shaped food products without 3D reconstruction based on Monte Carlo method with heuristic adjustment. Five images of food product were captured using five cameras and processed to produce binary images. Monte Carlo integration with heuristic adjustment was performed to measure the volume based on the information extracted from binary images. The experimental results show that the proposed method provided high accuracy and precision compared to the water displacement method. In addition, the proposed method is more accurate and faster than the space carving method.

Monte Carlo method with heuristic adjustment for irregularly shaped food product volume measurement.

Science.gov (United States)

Siswantoro, Joko; Prabuwono, Anton Satria; Abdullah, Azizi; Idrus, Bahari

2014-01-01

Volume measurement plays an important role in the production and processing of food products. Various methods have been proposed to measure the volume of food products with irregular shapes based on 3D reconstruction. However, 3D reconstruction comes with a high-priced computational cost. Furthermore, some of the volume measurement methods based on 3D reconstruction have a low accuracy. Another method for measuring volume of objects uses Monte Carlo method. Monte Carlo method performs volume measurements using random points. Monte Carlo method only requires information regarding whether random points fall inside or outside an object and does not require a 3D reconstruction. This paper proposes volume measurement using a computer vision system for irregularly shaped food products without 3D reconstruction based on Monte Carlo method with heuristic adjustment. Five images of food product were captured using five cameras and processed to produce binary images. Monte Carlo integration with heuristic adjustment was performed to measure the volume based on the information extracted from binary images. The experimental results show that the proposed method provided high accuracy and precision compared to the water displacement method. In addition, the proposed method is more accurate and faster than the space carving method.

An analytical-numerical comprehensive method for optimizing the fringing magnetic field

International Nuclear Information System (INIS)

Xiao Meiqin; Mao Naifeng

1991-01-01

The criterion of optimizing the fringing magnetic field is discussed, and an analytical-numerical comprehensive method for realizing the optimization is introduced. The method mentioned above consists of two parts, the analytical part calculates the field of the shims, which corrects the fringing magnetic field by using uniform magnetizing method; the numerical part fulfils the whole calculation of the field distribution by solving the equation of magnetic vector potential A within the region covered by arbitrary triangular meshes with the aid of finite difference method and successive over relaxation method. On the basis of the method, the optimization of the fringing magnetic field for a large-scale electromagnetic isotope separator is finished

Numerical methods for hyperbolic differential functional problems

Directory of Open Access Journals (Sweden)

Roman Ciarski

2008-01-01

Full Text Available The paper deals with the initial boundary value problem for quasilinear first order partial differential functional systems. A general class of difference methods for the problem is constructed. Theorems on the error estimate of approximate solutions for difference functional systems are presented. The convergence results are proved by means of consistency and stability arguments. A numerical example is given.

Numerical Study for a Large Volume Droplet on the Dual-rough Surface: Apparent Contact Angle, Contact Angle Hysteresis and Transition Barrier.

Science.gov (United States)

Dong, Jian; Jin, Yanli; Dong, He; Liu, Jiawei; Ye, Senbin

2018-06-14

The profile, apparent contact angle (ACA), contact angle hysteresis (CAH) and wetting state transmission energy barrier (WSTEB) are important static and dynamic properties of a large volume droplet on the hierarchical surface. Understanding them can provide us with important insights to functional surfaces and promote the application in corresponding areas. In this paper, we established three theoretical models (Model 1, Model 2 and Model 3) and corresponding numerical methods, which were obtained by the free energy minimization and the nonlinear optimization algorithm, to predict the profile, ACA, CAH and WSTEB of a large volume droplet on the horizontal regular dual-rough surface. In consideration of the gravity, the energy barrier on the contact circle, the dual heterogenous structures and their roughness on the surface, the models are more universal and accurate than previous models. It showed that the predictions of the models were in good agreement with the results from the experiment or literature. The models are promising to become novel design approaches of functional surfaces, which are frequently applied in microfluidic chips, water self-catchment system and dropwise condensation heat transfer system.

A virtual component method in numerical computation of cascades for isotope separation

International Nuclear Information System (INIS)

Zeng Shi; Cheng Lu

2014-01-01

The analysis, optimization, design and operation of cascades for isotope separation involve computations of cascades. In analytical analysis of cascades, using virtual components is a very useful analysis method. For complicated cases of cascades, numerical analysis has to be employed. However, bound up to the conventional idea that the concentration of a virtual component should be vanishingly small, virtual component is not yet applied to numerical computations. Here a method of introducing the method of using virtual components to numerical computations is elucidated, and its application to a few types of cascades is explained and tested by means of numerical experiments. The results show that the concentration of a virtual component is not restrained at all by the 'vanishingly small' idea. For the same requirements on cascades, the cascades obtained do not depend on the concentrations of virtual components. (authors)

Method of measuring a liquid pool volume

Science.gov (United States)

Garcia, G.V.; Carlson, N.M.; Donaldson, A.D.

1991-03-19

A method of measuring a molten metal liquid pool volume and in particular molten titanium liquid pools is disclosed, including the steps of (a) generating an ultrasonic wave at the surface of the molten metal liquid pool, (b) shining a light on the surface of a molten metal liquid pool, (c) detecting a change in the frequency of light, (d) detecting an ultrasonic wave echo at the surface of the molten metal liquid pool, and (e) computing the volume of the molten metal liquid. 3 figures.

Numerical methods for axisymmetric and 3D nonlinear beams

Science.gov (United States)

Pinton, Gianmarco F.; Trahey, Gregg E.

2005-04-01

Time domain algorithms that solve the Khokhlov--Zabolotzskaya--Kuznetsov (KZK) equation are described and implemented. This equation represents the propagation of finite amplitude sound beams in a homogenous thermoviscous fluid for axisymmetric and fully three dimensional geometries. In the numerical solution each of the terms is considered separately and the numerical methods are compared with known solutions. First and second order operator splitting are used to combine the separate terms in the KZK equation and their convergence is examined.

Methods for enhancing numerical integration

International Nuclear Information System (INIS)

Doncker, Elise de

2003-01-01

We give a survey of common strategies for numerical integration (adaptive, Monte-Carlo, Quasi-Monte Carlo), and attempt to delineate their realm of applicability. The inherent accuracy and error bounds for basic integration methods are given via such measures as the degree of precision of cubature rules, the index of a family of lattice rules, and the discrepancy of uniformly distributed point sets. Strategies incorporating these basic methods often use paradigms to reduce the error by, e.g., increasing the number of points in the domain or decreasing the mesh size, locally or uniformly. For these processes the order of convergence of the strategy is determined by the asymptotic behavior of the error, and may be too slow in practice for the type of problem at hand. For certain problem classes we may be able to improve the effectiveness of the method or strategy by such techniques as transformations, absorbing a difficult part of the integrand into a weight function, suitable partitioning of the domain, transformations and extrapolation or convergence acceleration. Situations warranting the use of these techniques (possibly in an 'automated' way) are described and illustrated by sample applications

Valve cam design using numerical step-by-step method

OpenAIRE

Vasilyev, Aleksandr; Bakhracheva, Yuliya; Kabore, Ousman; Zelenskiy, Yuriy

2014-01-01

This article studies the numerical step-by-step method of cam profile design. The results of the study are used for designing the internal combustion engine valve gear. This method allows to profile the peak efficiency of cams in view of many restrictions, connected with valve gear serviceability and reliability.

Investigating Convergence Patterns for Numerical Methods Using Data Analysis

Science.gov (United States)

Gordon, Sheldon P.

2013-01-01

The article investigates the patterns that arise in the convergence of numerical methods, particularly those in the errors involved in successive iterations, using data analysis and curve fitting methods. In particular, the results obtained are used to convey a deeper level of understanding of the concepts of linear, quadratic, and cubic…

Solutions manual to accompany An introduction to numerical methods and analysis

CERN Document Server

Epperson, James F

2014-01-01

A solutions manual to accompany An Introduction to Numerical Methods and Analysis, Second Edition An Introduction to Numerical Methods and Analysis, Second Edition reflects the latest trends in the field, includes new material and revised exercises, and offers a unique emphasis on applications. The author clearly explains how to both construct and evaluate approximations for accuracy and performance, which are key skills in a variety of fields. A wide range of higher-level methods and solutions, including new topics such as the roots of polynomials, sp

Finite volume multigrid method of the planar contraction flow of a viscoelastic fluid

Science.gov (United States)

Moatssime, H. Al; Esselaoui, D.; Hakim, A.; Raghay, S.

2001-08-01

This paper reports on a numerical algorithm for the steady flow of viscoelastic fluid. The conservative and constitutive equations are solved using the finite volume method (FVM) with a hybrid scheme for the velocities and first-order upwind approximation for the viscoelastic stress. A non-uniform staggered grid system is used. The iterative SIMPLE algorithm is employed to relax the coupled momentum and continuity equations. The non-linear algebraic equations over the flow domain are solved iteratively by the symmetrical coupled Gauss-Seidel (SCGS) method. In both, the full approximation storage (FAS) multigrid algorithm is used. An Oldroyd-B fluid model was selected for the calculation. Results are reported for planar 4:1 abrupt contraction at various Weissenberg numbers. The solutions are found to be stable and smooth. The solutions show that at high Weissenberg number the domain must be long enough. The convergence of the method has been verified with grid refinement. All the calculations have been performed on a PC equipped with a Pentium III processor at 550 MHz. Copyright

A detailed survey of numerical methods for unconstrained minimization. Pt. 1

International Nuclear Information System (INIS)

Mika, K.; Chaves, T.

1980-01-01

A detailed description of numerical methods for unconstrained minimization is presented. This first part surveys in particular conjugate direction and gradient methods, whereas variable metric methods will be the subject of the second part. Among the results of special interest we quote the following. The conjugate direction methods of Powell, Zangwill and Sutti can be best interpreted if the Smith approach is adopted. The conditions for quadratic termination of Powell's first procedure are analyzed. Numerical results based on nonlinear least squares problems are presented for the following conjugate direction codes: VA04AD from Harwell Subroutine Library and ZXPOW from IMSL, both implementations of Powell's second procedure, DFMND from IBM-SILMATH (Zangwill's method) and Brent's algorithm PRAXIS. VA04AD turns out to be superior in all cases, PRAXIS improves for high-dimensional problems. All codes clearly exhibit superlinear convergence. Akaike's result for the method of steepest descent is derived directly from a set of nonlinear recurrence relations. Numerical results obtained with the highly ill conditioned Hilbert function confirm the theoretical predictions. Several properties of the conjugate gradient method are presented and a new derivation of the equivalence of steepest descent partan and the CG method is given. A comparison of numerical results from the CG codes VA08AD (Fletcher-Reeves), DFMCG (the SSP version of the Fletcher-Reevens algorithm) and VA14AD (Powell's implementation of the Polak-Ribiere formula) reveals that VA14AD is clearly superior in all cases, but that the convergence rate of these codes is only weakly superlinear such that high accuracy solutions require extremely large numbers of function calls. (orig.)

A numerical method for two-dimensional anisotropic transport problem in cylindrical geometry

International Nuclear Information System (INIS)

Du Mingsheng; Feng Tiekai; Fu Lianxiang; Cao Changshu; Liu Yulan

1988-01-01

The authors deal with the triangular mesh-discontinuous finite element method for solving the time-dependent anisotropic neutron transport problem in two-dimensional cylindrical geometry. A prior estimate of the numerical solution is given. Stability is proved. The authors have computed a two dimensional anisotropic neutron transport problem and a Tungsten-Carbide critical assembly problem by using the numerical method. In comparision with DSN method and the experimental results obtained by others both at home and abroad, the method is satisfactory

MATH: A Scientific Tool for Numerical Methods Calculation and Visualization

Directory of Open Access Journals (Sweden)

Henrich Glaser-Opitz

2016-02-01

Full Text Available MATH is an easy to use application for various numerical methods calculations with graphical user interface and integrated plotting tool written in Qt with extensive use of Qwt library for plotting options and use of Gsl and MuParser libraries as a numerical and parser helping libraries. It can be found at http://sourceforge.net/projects/nummath. MATH is a convenient tool for use in education process because of its capability of showing every important step in solution process to better understand how it is done. MATH also enables fast comparison of similar method speed and precision.

Stress strain modelling of casting processes in the framework of the control volume method

DEFF Research Database (Denmark)

Hattel, Jesper; Andersen, Søren; Thorborg, Jesper

1998-01-01

Realistic computer simulations of casting processes call for the solution of both thermal, fluid-flow and stress/strain related problems. The multitude of the influencing parameters, and their non-linear, transient and temperature dependent nature, make the calculations complex. Therefore the nee......, the present model is based on the mainly decoupled representation of the thermal, mechanical and microstructural processes. Examples of industrial applications, such as predicting residual deformations in castings and stress levels in die casting dies, are presented...... for fast, flexible, multidimensional numerical methods is obvious. The basis of the deformation and stress/strain calculation is a transient heat transfer analysis including solidification. This paper presents an approach where the stress/strain and the heat transfer analysis uses the same computational...... domain, which is highly convenient. The basis of the method is the control volume finite difference approach on structured meshes. The basic assumptions of the method are shortly reviewed and discussed. As for other methods which aim at application oriented analysis of casting deformations and stresses...

Microstructure Optimization of Dual-Phase Steels Using a Representative Volume Element and a Response Surface Method: Parametric Study

Science.gov (United States)

Belgasam, Tarek M.; Zbib, Hussein M.

2017-12-01

Dual-phase (DP) steels have received widespread attention for their low density and high strength. This low density is of value to the automotive industry for the weight reduction it offers and the attendant fuel savings and emission reductions. Recent studies on developing DP steels showed that the combination of strength/ductility could be significantly improved when changing the volume fraction and grain size of phases in the microstructure depending on microstructure properties. Consequently, DP steel manufacturers are interested in predicting microstructure properties and in optimizing microstructure design. In this work, a microstructure-based approach using representative volume elements (RVEs) was developed. The approach examined the flow behavior of DP steels using virtual tension tests with an RVE to identify specific mechanical properties. Microstructures with varied martensite and ferrite grain sizes, martensite volume fractions, carbon content, and morphologies were studied in 3D RVE approaches. The effect of these microstructure parameters on a combination of strength/ductility of DP steels was examined numerically using the finite element method by implementing a dislocation density-based elastic-plastic constitutive model, and a Response surface methodology to determine the optimum conditions for a required combination of strength/ductility. The results from the numerical simulations are compared with experimental results found in the literature. The developed methodology proves to be a powerful tool for studying the effect and interaction of key microstructural parameters on strength and ductility and thus can be used to identify optimum microstructural conditions.

Parametric methods outperformed non-parametric methods in comparisons of discrete numerical variables

Directory of Open Access Journals (Sweden)

Sandvik Leiv

2011-04-01

Full Text Available Abstract Background The number of events per individual is a widely reported variable in medical research papers. Such variables are the most common representation of the general variable type called discrete numerical. There is currently no consensus on how to compare and present such variables, and recommendations are lacking. The objective of this paper is to present recommendations for analysis and presentation of results for discrete numerical variables. Methods Two simulation studies were used to investigate the performance of hypothesis tests and confidence interval methods for variables with outcomes {0, 1, 2}, {0, 1, 2, 3}, {0, 1, 2, 3, 4}, and {0, 1, 2, 3, 4, 5}, using the difference between the means as an effect measure. Results The Welch U test (the T test with adjustment for unequal variances and its associated confidence interval performed well for almost all situations considered. The Brunner-Munzel test also performed well, except for small sample sizes (10 in each group. The ordinary T test, the Wilcoxon-Mann-Whitney test, the percentile bootstrap interval, and the bootstrap-t interval did not perform satisfactorily. Conclusions The difference between the means is an appropriate effect measure for comparing two independent discrete numerical variables that has both lower and upper bounds. To analyze this problem, we encourage more frequent use of parametric hypothesis tests and confidence intervals.

Numerical method for partial equilibrium flow

International Nuclear Information System (INIS)

Ramshaw, J.D.; Cloutman, L.D.; Los Alamos, New Mexico 87545)

1981-01-01

A numerical method is presented for chemically reactive fluid flow in which equilibrium and nonequilibrium reactions occur simultaneously. The equilibrium constraints on the species concentrations are established by a quadratic iterative procedure. If the equilibrium reactions are uncoupled and of second or lower order, the procedure converges in a single step. In general, convergence is most rapid when the reactions are weakly coupled. This can frequently be achieved by a judicious choice of the independent reactions. In typical transient calculations, satisfactory accuracy has been achieved with about five iterations per time step

On the numerical stability analysis of pipelined Krylov subspace methods

Czech Academy of Sciences Publication Activity Database

Carson, E.T.; Rozložník, Miroslav; Strakoš, Z.; Tichý, P.; Tůma, M.

submitted 2017 (2018) R&D Projects: GA ČR GA13-06684S Grant - others:GA MŠk(CZ) LL1202 Institutional support: RVO:67985807 Keywords : Krylov subspace methods * the conjugate gradient method * numerical stability * inexact computations * delay of convergence * maximal attainable accuracy * pipelined Krylov subspace methods * exascale computations

Hydrothermal analysis in engineering using control volume finite element method

CERN Document Server

Sheikholeslami, Mohsen

2015-01-01

Control volume finite element methods (CVFEM) bridge the gap between finite difference and finite element methods, using the advantages of both methods for simulation of multi-physics problems in complex geometries. In Hydrothermal Analysis in Engineering Using Control Volume Finite Element Method, CVFEM is covered in detail and applied to key areas of thermal engineering. Examples, exercises, and extensive references are used to show the use of the technique to model key engineering problems such as heat transfer in nanofluids (to enhance performance and compactness of energy systems),

Numerical experiments using deflation with the HISQ action

Science.gov (United States)

Davies, Christine; DeTar, Carleton; McNeile, Craig; Vaquero, Alejandro

2018-03-01

We report on numerical experiments using deflation to compute quark propagators for the highly improved staggered quark (HISQ) action. The method is tested on HISQ gauge configurations, generated by the MILC collaboration, with lattice spacings of 0.15 fm, with a range of volumes, and sea quark masses down to the physical quark mass.

ABC/2 Method Does not Accurately Predict Cerebral Arteriovenous Malformation Volume.

Science.gov (United States)

Roark, Christopher; Vadlamudi, Venu; Chaudhary, Neeraj; Gemmete, Joseph J; Seinfeld, Joshua; Thompson, B Gregory; Pandey, Aditya S

2018-02-01

Stereotactic radiosurgery (SRS) is a treatment option for cerebral arteriovenous malformations (AVMs) to prevent intracranial hemorrhage. The decision to proceed with SRS is usually based on calculated nidal volume. Physicians commonly use the ABC/2 formula, based on digital subtraction angiography (DSA), when counseling patients for SRS. To determine whether AVM volume calculated using the ABC/2 method on DSA is accurate when compared to the exact volume calculated from thin-cut axial sections used for SRS planning. Retrospective search of neurovascular database to identify AVMs treated with SRS from 1995 to 2015. Maximum nidal diameters in orthogonal planes on DSA images were recorded to determine volume using ABC/2 formula. Nidal target volume was extracted from operative reports of SRS. Volumes were then compared using descriptive statistics and paired t-tests. Ninety intracranial AVMs were identified. Median volume was 4.96 cm3 [interquartile range (IQR) 1.79-8.85] with SRS planning methods and 6.07 cm3 (IQR 1.3-13.6) with ABC/2 methodology. Moderate correlation was seen between SRS and ABC/2 (r = 0.662; P ABC/2 (t = -3.2; P = .002). When AVMs were dichotomized based on ABC/2 volume, significant differences remained (t = 3.1, P = .003 for ABC/2 volume ABC/2 volume > 7 cm3). The ABC/2 method overestimates cerebral AVM volume when compared to volumetric analysis from SRS planning software. For AVMs > 7 cm3, the overestimation is even greater. SRS planning techniques were also significantly different than values derived from equations for cones and cylinders. Copyright © 2017 by the Congress of Neurological Surgeons

Energy conserving numerical methods for the computation of complex vortical flows

Science.gov (United States)

Allaneau, Yves

One of the original goals of this thesis was to develop numerical tools to help with the design of micro air vehicles. Micro Air Vehicles (MAVs) are small flying devices of only a few inches in wing span. Some people consider that as their size becomes smaller and smaller, it would be increasingly more difficult to keep all the classical control surfaces such as the rudders, the ailerons and the usual propellers. Over the years, scientists took inspiration from nature. Birds, by flapping and deforming their wings, are capable of accurate attitude control and are able to generate propulsion. However, the biomimicry design has its own limitations and it is difficult to place a hummingbird in a wind tunnel to study precisely the motion of its wings. Our approach was to use numerical methods to tackle this challenging problem. In order to precisely evaluate the lift and drag generated by the wings, one needs to be able to capture with high fidelity the extremely complex vortical flow produced in the wake. This requires a numerical method that is stable yet not too dissipative, so that the vortices do not get diffused in an unphysical way. We solved this problem by developing a new Discontinuous Galerkin scheme that, in addition to conserving mass, momentum and total energy locally, also preserves kinetic energy globally. This property greatly improves the stability of the simulations, especially in the special case p=0 when the approximation polynomials are taken to be piecewise constant (we recover a finite volume scheme). In addition to needing an adequate numerical scheme, a high fidelity solution requires many degrees of freedom in the computations to represent the flow field. The size of the smallest eddies in the flow is given by the Kolmogoroff scale. Capturing these eddies requires a mesh counting in the order of Re³ cells, where Re is the Reynolds number of the flow. We show that under-resolving the system, to a certain extent, is acceptable. However our

Workshop on Numerical Methods for Ordinary Differential Equations

CERN Document Server

Gear, Charles; Russo, Elvira

1989-01-01

Developments in numerical initial value ode methods were the focal topic of the meeting at L'Aquila which explord the connections between the classical background and new research areas such as differental-algebraic equations, delay integral and integro-differential equations, stability properties, continuous extensions (interpolants for Runge-Kutta methods and their applications, effective stepsize control, parallel algorithms for small- and large-scale parallel architectures). The resulting proceedings address many of these topics in both research and survey papers.

Supplier Portfolio Selection and Optimum Volume Allocation: A Knowledge Based Method

NARCIS (Netherlands)

Aziz, Romana; Aziz, R.; van Hillegersberg, Jos; Kersten, W.; Blecker, T.; Luthje, C.

2010-01-01

Selection of suppliers and allocation of optimum volumes to suppliers is a strategic business decision. This paper presents a decision support method for supplier selection and the optimal allocation of volumes in a supplier portfolio. The requirements for the method were gathered during a case

Numerical methods for the simulation of continuous sedimentation in ideal clarifier-thickener units

Energy Technology Data Exchange (ETDEWEB)

Buerger, R.; Karlsen, K.H.; Risebro, N.H.; Towers, J.D.

2001-10-01

We consider a model of continuous sedimentation. Under idealizing assumptions, the settling of the solid particles under the influence of gravity can be described by the initial value problem for a nonlinear hyperbolic partial differential equation with a flux function that depends discontinuously on height. The purpose of this contribution is to present and demonstrate two numerical methods for simulating continuous sedimentation: a front tracking method and a finite finite difference method. The basic building blocks in the front tracking method are the solutions of a finite number of certain Riemann problems and a procedure for tracking local collisions of shocks. The solutions of the Riemann problems are recalled herein and the front tracking algorithm is described. As an alternative to the front tracking method, a simple scalar finite difference algorithm is proposed. This method is based on discretizing the spatially varying flux parameters on a mesh that is staggered with respect to that of the conserved variable, resulting in a straightforward generalization of the well-known Engquist-Osher upwind finite difference method. The result is an easily implemented upwind shock capturing method. Numerical examples demonstrate that the front tracking and finite difference methods can be used as efficient and accurate simulation tools for continuous sedimentation. The numerical results for the finite difference method indicate that discontinuities in the local solids concentration are resolved sharply and agree with those produced by the front tracking method. The latter is free of numerical dissipation, which leads to sharply resolved concentration discontinuities, but is more complicated to implement than the former. Available mathematical results for the proposed numerical methods are also briefly reviewed. (author)

MRI definition of target volumes using fuzzy logic method for three-dimensional conformal radiation therapy

International Nuclear Information System (INIS)

Caudrelier, Jean-Michel; Vial, Stephane; Gibon, David; Kulik, Carine; Fournier, Charles; Castelain, Bernard; Coche-Dequeant, Bernard; Rousseau, Jean

2003-01-01

Purpose: Three-dimensional (3D) volume determination is one of the most important problems in conformal radiation therapy. Techniques of volume determination from tomographic medical imaging are usually based on two-dimensional (2D) contour definition with the result dependent on the segmentation method used, as well as on the user's manual procedure. The goal of this work is to describe and evaluate a new method that reduces the inaccuracies generally observed in the 2D contour definition and 3D volume reconstruction process. Methods and Materials: This new method has been developed by integrating the fuzziness in the 3D volume definition. It first defines semiautomatically a minimal 2D contour on each slice that definitely contains the volume and a maximal 2D contour that definitely does not contain the volume. The fuzziness region in between is processed using possibility functions in possibility theory. A volume of voxels, including the membership degree to the target volume, is then created on each slice axis, taking into account the slice position and slice profile. A resulting fuzzy volume is obtained after data fusion between multiorientation slices. Different studies have been designed to evaluate and compare this new method of target volume reconstruction and a classical reconstruction method. First, target definition accuracy and robustness were studied on phantom targets. Second, intra- and interobserver variations were studied on radiosurgery clinical cases. Results: The absolute volume errors are less than or equal to 1.5% for phantom volumes calculated by the fuzzy logic method, whereas the values obtained with the classical method are much larger than the actual volumes (absolute volume errors up to 72%). With increasing MRI slice thickness (1 mm to 8 mm), the phantom volumes calculated by the classical method are increasing exponentially with a maximum absolute error up to 300%. In contrast, the absolute volume errors are less than 12% for phantom

Analytic-numerical method of determining the freezing front location

Directory of Open Access Journals (Sweden)

R. Grzymkowski

2011-07-01

Full Text Available Mathematical modeling of thermal processes combined with the reversible phase transitions of type: solid phase – liquid phase leads to formulation of the parabolic boundary problems with the moving boundary. Solution of such defined problem requires, most often, to use sophisticated numerical techniques and far advanced mathematical tools. Excellent illustration of the complexity of considered problems, as well as of the variety of approaches used for finding their solutions, gives the papers [1-4]. In the current paper, the authors present the, especially attractive from the engineer point of view, analytic-numerical method for finding the approximate solution of selected class of problems which can be reduced to the one-phase solidification problem of a plate with the unknown a priori, varying in time boundary of the region in which the solution is sought. Proposed method is based on the known formalism of initial expansion of the sought function describing the temperature field into the power series, some coefficients of which are determined with the aid of boundary conditions, and on the approximation of the function defining the location of freezing front with the broken line, parameters of which are numerically determined.

Stability of numerical method for semi-linear stochastic pantograph differential equations

Directory of Open Access Journals (Sweden)

Yu Zhang

2016-01-01

Full Text Available Abstract As a particular expression of stochastic delay differential equations, stochastic pantograph differential equations have been widely used in nonlinear dynamics, quantum mechanics, and electrodynamics. In this paper, we mainly study the stability of analytical solutions and numerical solutions of semi-linear stochastic pantograph differential equations. Some suitable conditions for the mean-square stability of an analytical solution are obtained. Then we proved the general mean-square stability of the exponential Euler method for a numerical solution of semi-linear stochastic pantograph differential equations, that is, if an analytical solution is stable, then the exponential Euler method applied to the system is mean-square stable for arbitrary step-size h > 0 $h>0$ . Numerical examples further illustrate the obtained theoretical results.

A method for solving the KDV equation and some numerical experiments

International Nuclear Information System (INIS)

Chang Jinjiang.

1993-01-01

In this paper, by means of difference method for discretization of space partial derivatives of KDV equation, an initial value problem in ordinary differential equations of large dimensions is produced. By using this ordinary differential equations the existence and the uniqueness of the solution of the KDV equation and the conservation of scheme are proved. This ordinary differential equation can be solved by using implicit Runge-Kutta methods, so a new method for finding the numerical solution of the KDV equation is presented. Numerical experiments not only describe in detail the procedure of two solitons collision, soliton reflex and soliton produce, but also show that this method is very effective. (author). 7 refs, 3 figs

The finite volume element (FVE) and multigrid method for the incompressible Navier-Stokes equations

International Nuclear Information System (INIS)

Gu Lizhen; Bao Weizhu

1992-01-01

The authors apply FVE method to discrete INS equations with the original variable, in which the bilinear square finite element and the square finite volume are chosen. The discrete schemes of INS equations are presented. The FMV multigrid algorithm is applied to solve that discrete system, where DGS iteration is used as smoother, DGS distributive mode for the INS discrete system is also presented. The sample problems for the square cavity flow with Reynolds number Re≤100 are successfully calculated. The numerical solutions show that the results with 1 FMV is satisfactory and when Re is not large, The FVE discrete scheme of the conservative INS equations and that of non-conservative INS equations with linearization both can provide almost same accuracy

Fast multiview three-dimensional reconstruction method using cost volume filtering

Science.gov (United States)

Lee, Seung Joo; Park, Min Ki; Jang, In Yeop; Lee, Kwan H.

2014-03-01

As the number of customers who want to record three-dimensional (3-D) information using a mobile electronic device increases, it becomes more and more important to develop a method which quickly reconstructs a 3-D model from multiview images. A fast multiview-based 3-D reconstruction method is presented, which is suitable for the mobile environment by constructing a cost volume of the 3-D height field. This method consists of two steps: the construction of a reliable base surface and the recovery of shape details. In each step, the cost volume is constructed using photoconsistency and then it is filtered according to the multiscale. The multiscale-based cost volume filtering allows the 3-D reconstruction to maintain the overall shape and to preserve the shape details. We demonstrate the strength of the proposed method in terms of computation time, accuracy, and unconstrained acquisition environment.

An analytically based numerical method for computing view factors in real urban environments

Science.gov (United States)

Lee, Doo-Il; Woo, Ju-Wan; Lee, Sang-Hyun

2018-01-01

A view factor is an important morphological parameter used in parameterizing in-canyon radiative energy exchange process as well as in characterizing local climate over urban environments. For realistic representation of the in-canyon radiative processes, a complete set of view factors at the horizontal and vertical surfaces of urban facets is required. Various analytical and numerical methods have been suggested to determine the view factors for urban environments, but most of the methods provide only sky-view factor at the ground level of a specific location or assume simplified morphology of complex urban environments. In this study, a numerical method that can determine the sky-view factors ( ψ ga and ψ wa ) and wall-view factors ( ψ gw and ψ ww ) at the horizontal and vertical surfaces is presented for application to real urban morphology, which are derived from an analytical formulation of the view factor between two blackbody surfaces of arbitrary geometry. The established numerical method is validated against the analytical sky-view factor estimation for ideal street canyon geometries, showing a consolidate confidence in accuracy with errors of less than 0.2 %. Using a three-dimensional building database, the numerical method is also demonstrated to be applicable in determining the sky-view factors at the horizontal (roofs and roads) and vertical (walls) surfaces in real urban environments. The results suggest that the analytically based numerical method can be used for the radiative process parameterization of urban numerical models as well as for the characterization of local urban climate.

Safety Analysis in Large Volume Vacuum Systems Like Tokamak: Experiments and Numerical Simulation to Analyze Vacuum Ruptures Consequences

Directory of Open Access Journals (Sweden)

A. Malizia

2014-01-01

Full Text Available The large volume vacuum systems are used in many industrial operations and research laboratories. Accidents in these systems should have a relevant economical and safety impact. A loss of vacuum accident (LOVA due to a failure of the main vacuum vessel can result in a fast pressurization of the vessel and consequent mobilization dispersion of hazardous internal material through the braches. It is clear that the influence of flow fields, consequence of accidents like LOVA, on dust resuspension is a key safety issue. In order to develop this analysis an experimental facility is been developed: STARDUST. This last facility has been used to improve the knowledge about LOVA to replicate a condition more similar to appropriate operative condition like to kamaks. By the experimental data the boundary conditions have been extrapolated to give the proper input for the 2D thermofluid-dynamics numerical simulations, developed by the commercial CFD numerical code. The benchmark of numerical simulation results with the experimental ones has been used to validate and tune the 2D thermofluid-dynamics numerical model that has been developed by the authors to replicate the LOVA conditions inside STARDUST. In present work, the facility, materials, numerical model, and relevant results will be presented.

Protection Parameters against the Cracks by the Method of Volume Compensation Dam

Directory of Open Access Journals (Sweden)

Bulatov Georgiy

2016-01-01

Full Text Available This article provides estimates the parameters of protection from cracking dam due to volume compensation method. This article discusses the method of compensation dam volume. This method allows calculating the settings of security causing cracks the dam. Presents graphs of horizontal deformations of elongation calculated surface along the length of the construction and in time. Showing horizontal stress distribution diagram in the ground around the pile in plan and in section. Given all the necessary formulas for the method of compensation of the dam volume.

Evaluation of two-phase flow solvers using Level Set and Volume of Fluid methods

Science.gov (United States)

Bilger, C.; Aboukhedr, M.; Vogiatzaki, K.; Cant, R. S.

2017-09-01

Two principal methods have been used to simulate the evolution of two-phase immiscible flows of liquid and gas separated by an interface. These are the Level-Set (LS) method and the Volume of Fluid (VoF) method. Both methods attempt to represent the very sharp interface between the phases and to deal with the large jumps in physical properties associated with it. Both methods have their own strengths and weaknesses. For example, the VoF method is known to be prone to excessive numerical diffusion, while the basic LS method has some difficulty in conserving mass. Major progress has been made in remedying these deficiencies, and both methods have now reached a high level of physical accuracy. Nevertheless, there remains an issue, in that each of these methods has been developed by different research groups, using different codes and most importantly the implementations have been fine tuned to tackle different applications. Thus, it remains unclear what are the remaining advantages and drawbacks of each method relative to the other, and what might be the optimal way to unify them. In this paper, we address this gap by performing a direct comparison of two current state-of-the-art variations of these methods (LS: RCLSFoam and VoF: interPore) and implemented in the same code (OpenFoam). We subject both methods to a pair of benchmark test cases while using the same numerical meshes to examine a) the accuracy of curvature representation, b) the effect of tuning parameters, c) the ability to minimise spurious velocities and d) the ability to tackle fluids with very different densities. For each method, one of the test cases is chosen to be fairly benign while the other test case is expected to present a greater challenge. The results indicate that both methods can be made to work well on both test cases, while displaying different sensitivity to the relevant parameters.

A calculation method for RF couplers design based on numerical simulation by microwave studio

International Nuclear Information System (INIS)

Wang Rong; Pei Yuanji; Jin Kai

2006-01-01

A numerical simulation method for coupler design is proposed. It is based on the matching procedure for the 2π/3 structure given by Dr. R.L. Kyhl. Microwave Studio EigenMode Solver is used for such numerical simulation. the simulation for a coupler has been finished with this method and the simulation data are compared with experimental measurements. The results show that this numerical simulation method is feasible for coupler design. (authors)

Numerical simulation methods for wave propagation through optical waveguides

International Nuclear Information System (INIS)

Sharma, A.

1993-01-01

The simulation of the field propagation through waveguides requires numerical solutions of the Helmholtz equation. For this purpose a method based on the principle of orthogonal collocation was recently developed. The method is also applicable to nonlinear pulse propagation through optical fibers. Some of the salient features of this method and its application to both linear and nonlinear wave propagation through optical waveguides are discussed in this report. 51 refs, 8 figs, 2 tabs

Methods for determining enzymatic activity comprising heating and agitation of closed volumes

Science.gov (United States)

Thompson, David Neil; Henriksen, Emily DeCrescenzo; Reed, David William; Jensen, Jill Renee

2016-03-15

Methods for determining thermophilic enzymatic activity include heating a substrate solution in a plurality of closed volumes to a predetermined reaction temperature. Without opening the closed volumes, at least one enzyme is added, substantially simultaneously, to the closed volumes. At the predetermined reaction temperature, the closed volumes are agitated and then the activity of the at least one enzyme is determined. The methods are conducive for characterizing enzymes of high-temperature reactions, with insoluble substrates, with substrates and enzymes that do not readily intermix, and with low volumes of substrate and enzyme. Systems for characterizing the enzymes are also disclosed.

A higher order numerical method for time fractional partial differential equations with nonsmooth data

Science.gov (United States)

Xing, Yanyuan; Yan, Yubin

2018-03-01

Gao et al. [11] (2014) introduced a numerical scheme to approximate the Caputo fractional derivative with the convergence rate O (k 3 - α), 0 equation is sufficiently smooth, Lv and Xu [20] (2016) proved by using energy method that the corresponding numerical method for solving time fractional partial differential equation has the convergence rate O (k 3 - α), 0 equation has low regularity and in this case the numerical method fails to have the convergence rate O (k 3 - α), 0 quadratic interpolation polynomials. Based on this scheme, we introduce a time discretization scheme to approximate the time fractional partial differential equation and show by using Laplace transform methods that the time discretization scheme has the convergence rate O (k 3 - α), 0 0 for smooth and nonsmooth data in both homogeneous and inhomogeneous cases. Numerical examples are given to show that the theoretical results are consistent with the numerical results.

A Broyden numerical Kutta condition for an unsteady panel method

International Nuclear Information System (INIS)

Liu, P.; Bose, N.; Colbourne, B.

2003-01-01

In panel methods, numerical Kutta conditions are applied in order to ensure that pressure differences between the surfaces at the trailing edges of lifting surface elements are close to zero. Previous numerical Kutta conditions for 3-D panel methods have focused on use of the Newton-Raphson iterative procedure. For extreme unsteady motions, such as for oscillating hydrofoils or for a propeller behind a blockage, the Newton-Raphson procedure can have severe convergence difficulties. The Broyden iteration, a modified Newton-Raphson iteration procedure, is applied here to obtain improved convergence behavior. Using the Broyden iteration increases the reliability, robustness and in many cases computing efficiency for unsteady, multi-body interactive flows. This method was tested in a time domain code for an ice class propeller in both open water flow and during interaction with a nearby ice blockage. Predictions showed that the method was effective in these extreme flows. (author)

Comparative study on triangular and quadrilateral meshes by a finite-volume method with a central difference scheme

KAUST Repository

Yu, Guojun

2012-10-01

In this article, comparative studies on computational accuracies and convergence rates of triangular and quadrilateral meshes are carried out in the frame work of the finite-volume method. By theoretical analysis, we conclude that the number of triangular cells needs to be 4/3 times that of quadrilateral cells to obtain similar accuracy. The conclusion is verified by a number of numerical examples. In addition, the convergence rates of the triangular meshes are found to be slower than those of the quadrilateral meshes when the same accuracy is obtained with these two mesh types. © 2012 Taylor and Francis Group, LLC.

Comparative study on triangular and quadrilateral meshes by a finite-volume method with a central difference scheme

KAUST Repository

Yu, Guojun; Yu, Bo; Sun, Shuyu; Tao, Wenquan

2012-01-01

In this article, comparative studies on computational accuracies and convergence rates of triangular and quadrilateral meshes are carried out in the frame work of the finite-volume method. By theoretical analysis, we conclude that the number of triangular cells needs to be 4/3 times that of quadrilateral cells to obtain similar accuracy. The conclusion is verified by a number of numerical examples. In addition, the convergence rates of the triangular meshes are found to be slower than those of the quadrilateral meshes when the same accuracy is obtained with these two mesh types. © 2012 Taylor and Francis Group, LLC.

Numerical Verification Methods for Spherical $t$-Designs

OpenAIRE

Chen, Xiaojun

2009-01-01

The construction of spherical $t$-designs with $(t+1)^2$ points on the unit sphere $S^2$ in $\\mathbb{R}^3$ can be reformulated as an underdetermined system of nonlinear equations. This system is highly nonlinear and involves the evaluation of a degree $t$ polynomial in $(t+1)^4$ arguments. This paper reviews numerical verification methods using the Brouwer fixed point theorem and Krawczyk interval operator for solutions of the underdetermined system of nonlinear equations...

Numerical investigations on contactless methods for measuring critical current density in HTS: application of modified constitutive-relation method

International Nuclear Information System (INIS)

Kamitani, A.; Takayama, T.; Itoh, T.; Ikuno, S.

2011-01-01

A fast method is proposed for calculating the shielding current density in an HTS. The J-E constitutive relation is modified so as not to change the solution. A numerical code is developed on the basis of the proposed method. The permanent magnet method is successfully simulated by means of the code. A fast method has been proposed for calculating the shielding current density in a high-temperature superconducting thin film. An initial-boundary-value problem of the shielding current density cannot be always solved by means of the Runge-Kutta method even when an adaptive step-size control algorithm is incorporated to the method. In order to suppress an overflow in the algorithm, the J-E constitutive relation is modified so that its solution may satisfy the original constitutive relation. A numerical code for analyzing the shielding current density has been developed on the basis of this method and, as an application of the code, the permanent magnet method for measuring the critical current density has been investigated numerically.

On a method of numerical calculation of nonlinear radial pulsations of stars

International Nuclear Information System (INIS)

Kosovichev, A.G.

1984-01-01

Some features of using the finite difference method for numerical investigation of nonradial pulsations of stars were considered. The mathematical model of these pulsations is described by time-dependent gasdynaMic equations with gravity. A one-dimentional (spherically-symmetric) case is considered. It was obtained a two-parametric family of ultimate conservative difference schemes where the diffepence analogy of the main conservative laws as well as the additional relations for the balance to individual kinds of energy are performed. Such difference schemes provide more exact calculation of nonlinear flows with shocks as compared with the other difference schemes with the same order of approximation. The methods of numerical solution of implicit (absolute stable) difference schemes for a given family were considered. The coupled equations are solved through iterative Newton method Using martrix and separate successive eliminations. Numerical method can be used for calculation of large amplitude radial pulsations of stars

Development of numerical methods for reactive transport; Developpement de methodes numeriques pour le transport reactif

Energy Technology Data Exchange (ETDEWEB)

Bouillard, N

2006-12-15

When a radioactive waste is stored in deep geological disposals, it is expected that the waste package will be damaged under water action (concrete leaching, iron corrosion). Then, to understand these damaging processes, chemical reactions and solutes transport are modelled. Numerical simulations of reactive transport can be done sequentially by the coupling of several codes. This is the case of the software platform ALLIANCES which is developed jointly with CEA, ANDRA and EDF. Stiff reactions like precipitation-dissolution are crucial for the radioactive waste storage applications, but standard sequential iterative approaches like Picard's fail in solving rapidly reactive transport simulations with such stiff reactions. In the first part of this work, we focus on a simplified precipitation and dissolution process: a system made up with one solid species and two aqueous species moving by diffusion is studied mathematically. It is assumed that a precipitation dissolution reaction occurs in between them, and it is modelled by a discontinuous kinetics law of unknown sign. By using monotonicity properties, the convergence of a finite volume scheme on admissible mesh is proved. Existence of a weak solution is obtained as a by-product of the convergence of the scheme. The second part is dedicated to coupling algorithms which improve Picard's method and can be easily used in an existing coupling code. By extending previous works, we propose a general and adaptable framework to solve nonlinear systems. Indeed by selecting special options, we can either recover well known methods, like nonlinear conjugate gradient methods, or design specific method. This algorithm has two main steps, a preconditioning one and an acceleration one. This algorithm is tested on several examples, some of them being rather academical and others being more realistic. We test it on the 'three species model'' example. Other reactive transport simulations use an external

Development of numerical methods for reactive transport; Developpement de methodes numeriques pour le transport reactif

Energy Technology Data Exchange (ETDEWEB)

Bouillard, N

2006-12-15

When a radioactive waste is stored in deep geological disposals, it is expected that the waste package will be damaged under water action (concrete leaching, iron corrosion). Then, to understand these damaging processes, chemical reactions and solutes transport are modelled. Numerical simulations of reactive transport can be done sequentially by the coupling of several codes. This is the case of the software platform ALLIANCES which is developed jointly with CEA, ANDRA and EDF. Stiff reactions like precipitation-dissolution are crucial for the radioactive waste storage applications, but standard sequential iterative approaches like Picard's fail in solving rapidly reactive transport simulations with such stiff reactions. In the first part of this work, we focus on a simplified precipitation and dissolution process: a system made up with one solid species and two aqueous species moving by diffusion is studied mathematically. It is assumed that a precipitation dissolution reaction occurs in between them, and it is modelled by a discontinuous kinetics law of unknown sign. By using monotonicity properties, the convergence of a finite volume scheme on admissible mesh is proved. Existence of a weak solution is obtained as a by-product of the convergence of the scheme. The second part is dedicated to coupling algorithms which improve Picard's method and can be easily used in an existing coupling code. By extending previous works, we propose a general and adaptable framework to solve nonlinear systems. Indeed by selecting special options, we can either recover well known methods, like nonlinear conjugate gradient methods, or design specific method. This algorithm has two main steps, a preconditioning one and an acceleration one. This algorithm is tested on several examples, some of them being rather academical and others being more realistic. We test it on the 'three species model'' example. Other reactive transport simulations use an external chemical code CHESS. For a

Development of the numerical method for liquid metal magnetohydrodynamics (I). Investigation of the method and development of the 2D method

International Nuclear Information System (INIS)

Ohira, H.; Ara, K.

2002-11-01

Advanced electromagnetic components are investigated in Feasibility Studies on Commercialized FR Cycle System to apply to the main cooling systems of Liquid Metal Fast Reactor. Although a lot of experiments and numerical analysis were carried out on both high Reynolds numbers and high magnetic Reynolds numbers, the complex phenomena could not be evaluated in detail. As the first step of the development of the numerical methods for the liquid metal magnetohydrodynamics, we investigated numerical methods that could be applied to the electromagnetic components with both complex structures and high magnetic turbulent field. As a result, we selected GSMAC (Generalized-Simplified MArker and Cell) method for calculating the liquid metal fluid dynamics because it could be easily applied to the complex flow field. We also selected the vector-FEM for calculating the magnetic field of the large components because the method had no interaction procedure. In the high magnetic turbulent field, the dynamic-SGS models would be also a promising model for the good estimation, because it could calculate the field directly without any experimental constant. In order to verify the GSMAC and the vector-FEM, we developed the 2D numerical models and calculated the magnetohydrodynamics in the large electromagnetic pump. It was estimated from these results that the methods were basically reasonable, because the calculated pressure differences had the similar tendencies to the experimental ones. (author)

Post-Processing in the Material-Point Method

DEFF Research Database (Denmark)

Andersen, Søren; Andersen, Lars Vabbersgaard

The material-point method (MPM) is a numerical method for dynamic or static analysis of solids using a discretization in time and space. The method has shown to be successful in modelling physical problems involving large deformations, which are difficult to model with traditional numerical tools...... such as the finite element method. In the material-point method, a set of material points is utilized to track the problem in time and space, while a computational background grid is utilized to obtain spatial derivatives relevant to the physical problem. Currently, the research within the material-point method......-point method. The first idea involves associating a volume with each material point and displaying the deformation of this volume. In the discretization process, the physical domain is divided into a number of smaller volumes each represented by a simple shape; here quadrilaterals are chosen for the presented...

Calculating regional tissue volume for hyperthermic isolated limb perfusion: Four methods compared.

Science.gov (United States)

Cecchin, D; Negri, A; Frigo, A C; Bui, F; Zucchetta, P; Bodanza, V; Gregianin, M; Campana, L G; Rossi, C R; Rastrelli, M

2016-12-01

Hyperthermic isolated limb perfusion (HILP) can be performed as an alternative to amputation for soft tissue sarcomas and melanomas of the extremities. Melphalan and tumor necrosis factor-alpha are used at a dosage that depends on the volume of the limb. Regional tissue volume is traditionally measured for the purposes of HILP using water displacement volumetry (WDV). Although this technique is considered the gold standard, it is time-consuming and complicated to implement, especially in obese and elderly patients. The aim of the present study was to compare the different methods described in the literature for calculating regional tissue volume in the HILP setting, and to validate an open source software. We reviewed the charts of 22 patients (11 males and 11 females) who had non-disseminated melanoma with in-transit metastases or sarcoma of the lower limb. We calculated the volume of the limb using four different methods: WDV, tape measurements and segmentation of computed tomography images using Osirix and Oncentra Masterplan softwares. The overall comparison provided a concordance correlation coefficient (CCC) of 0.92 for the calculations of whole limb volume. In particular, when Osirix was compared with Oncentra (validated for volume measures and used in radiotherapy), the concordance was near-perfect for the calculation of the whole limb volume (CCC = 0.99). With methods based on CT the user can choose a reliable plane for segmentation purposes. CT-based methods also provides the opportunity to separate the whole limb volume into defined tissue volumes (cortical bone, fat and water). Copyright © 2016 Elsevier Ltd. All rights reserved.

Advances in the discrete ordinates and finite volume methods for the solution of radiative heat transfer problems in participating media

International Nuclear Information System (INIS)

Coelho, Pedro J.

2014-01-01

Many methods are available for the solution of radiative heat transfer problems in participating media. Among these, the discrete ordinates method (DOM) and the finite volume method (FVM) are among the most widely used ones. They provide a good compromise between accuracy and computational requirements, and they are relatively easy to integrate in CFD codes. This paper surveys recent advances on these numerical methods. Developments concerning the grid structure (e.g., new formulations for axisymmetrical geometries, body-fitted structured and unstructured meshes, embedded boundaries, multi-block grids, local grid refinement), the spatial discretization scheme, and the angular discretization scheme are described. Progress related to the solution accuracy, solution algorithm, alternative formulations, such as the modified DOM and FVM, even-parity formulation, discrete-ordinates interpolation method and method of lines, and parallelization strategies is addressed. The application to non-gray media, variable refractive index media, and transient problems is also reviewed. - Highlights: • We survey recent advances in the discrete ordinates and finite volume methods. • Developments in spatial and angular discretization schemes are described. • Progress in solution algorithms and parallelization methods is reviewed. • Advances in the transient solution of the radiative transfer equation are appraised. • Non-gray media and variable refractive index media are briefly addressed

New numerical method for iterative or perturbative solution of quantum field theory

International Nuclear Information System (INIS)

Hahn, S.C.; Guralnik, G.S.

1999-01-01

A new computational idea for continuum quantum Field theories is outlined. This approach is based on the lattice source Galerkin methods developed by Garcia, Guralnik and Lawson. The method has many promising features including treating fermions on a relatively symmetric footing with bosons. As a spin-off of the technology developed for 'exact' solutions, the numerical methods used have a special case application to perturbation theory. We are in the process of developing an entirely numerical approach to evaluating graphs to high perturbative order. (authors)

Development of Numerical Technique to Analyze the Flow Characteristics of Porous Media Using Lattice Boltzmann Method

Energy Technology Data Exchange (ETDEWEB)

Kim, Hyung Min [Kyonggi Univ., Suwon (Korea, Republic of)

2016-11-15

The performance of proton exchange membrane fuel cells (PEMFC) is strongly related to the water flow and accumulation in the gas diffusion layer (GDL) and catalyst layer. Understanding the behavior of fluid from the characteristics of the media is crucial for the improvement of the performance and design of the GDL. In this paper, a numerical method is proposed to calculate the design parameters of the GDL, i.e., permeability, tortuosity, and effective diffusivity. The fluid flow in a channel filled with randomly packed hard spheres is simulated to validate the method. The flow simulation was performed by lattice Boltzmann method with bounce back condition for the solid volume fraction in the porous media, with different values of porosities. Permeability, which affects the flow, was calculated from the average pressure drop and the velocity in the porous media. Tortuosity, calculated by the ratio the average path length of the randomly injected massless particles to the thickness of the porous media, and the resultant effective diffusivity were in good agreement with the theoretical model. The suggested method can be used to calculate the parameters of real GDL accurately without any modification.

A two-dimensional adaptive numerical grids generation method and its realization

International Nuclear Information System (INIS)

Xu Tao; Shui Hongshou

1998-12-01

A two-dimensional adaptive numerical grids generation method and its particular realization is discussed. This method is effective and easy to realize if the control functions are given continuously, and the grids for some regions is showed in this case. For Computational Fluid Dynamics, because the control values of adaptive grids-numerical solution is given in dispersed form, it is needed to interpolate these values to get the continuous control functions. These interpolation techniques are discussed, and some efficient adaptive grids are given. A two-dimensional fluid dynamics example was also given

Discrete variational derivative method a structure-preserving numerical method for partial differential equations

CERN Document Server

Furihata, Daisuke

2010-01-01

Nonlinear Partial Differential Equations (PDEs) have become increasingly important in the description of physical phenomena. Unlike Ordinary Differential Equations, PDEs can be used to effectively model multidimensional systems. The methods put forward in Discrete Variational Derivative Method concentrate on a new class of ""structure-preserving numerical equations"" which improves the qualitative behaviour of the PDE solutions and allows for stable computing. The authors have also taken care to present their methods in an accessible manner, which means that the book will be useful to engineer

Finite Volume Method for Unstructured Grid

International Nuclear Information System (INIS)

Casmara; Kardana, N.D.

1997-01-01

The success of a computational method depends on the solution algorithm and mesh generation techniques. cell distributions are needed, which allow the solution to be calculated over the entire body surface with sufficient accuracy. to handle the mesh generation for multi-connected region such as multi-element bodies, the unstructured finite volume method will be applied. the advantages of the unstructured meshes are it provides a great deal more flexibility for generating meshes about complex geometries and provides a natural setting for the use of adaptive meshing. the governing equations to be discretized are inviscid and rotational euler equations. Applications of the method will be evaluated on flow around single and multi-component bodies

A numeric-analytic method for approximating the chaotic Chen system

International Nuclear Information System (INIS)

Mossa Al-sawalha, M.; Noorani, M.S.M.

2009-01-01

The epitome of this paper centers on the application of the differential transformation method (DTM) the renowned Chen system which is described as a three-dimensional system of ODEs with quadratic nonlinearities. Numerical comparisons are made between the DTM and the classical fourth-order Runge-Kutta method (RK4). Our work showcases the precision of the DTM as the Chen system transforms from a non-chaotic system to a chaotic one. Since the Lyapunov exponent for this system is much higher compared to other chaotic systems, we shall highlight the difficulties of the simulations with respect to its accuracy. We wrap up our investigations to reveal that this direct symbolic-numeric scheme is effective and accurate.

Numerical analysis for multi-group neutron-diffusion equation using Radial Point Interpolation Method (RPIM)

International Nuclear Information System (INIS)

Kim, Kyung-O; Jeong, Hae Sun; Jo, Daeseong

2017-01-01

Highlights: • Employing the Radial Point Interpolation Method (RPIM) in numerical analysis of multi-group neutron-diffusion equation. • Establishing mathematical formation of modified multi-group neutron-diffusion equation by RPIM. • Performing the numerical analysis for 2D critical problem. - Abstract: A mesh-free method is introduced to overcome the drawbacks (e.g., mesh generation and connectivity definition between the meshes) of mesh-based (nodal) methods such as the finite-element method and finite-difference method. In particular, the Point Interpolation Method (PIM) using a radial basis function is employed in the numerical analysis for the multi-group neutron-diffusion equation. The benchmark calculations are performed for the 2D homogeneous and heterogeneous problems, and the Multiquadrics (MQ) and Gaussian (EXP) functions are employed to analyze the effect of the radial basis function on the numerical solution. Additionally, the effect of the dimensionless shape parameter in those functions on the calculation accuracy is evaluated. According to the results, the radial PIM (RPIM) can provide a highly accurate solution for the multiplication eigenvalue and the neutron flux distribution, and the numerical solution with the MQ radial basis function exhibits the stable accuracy with respect to the reference solutions compared with the other solution. The dimensionless shape parameter directly affects the calculation accuracy and computing time. Values between 1.87 and 3.0 for the benchmark problems considered in this study lead to the most accurate solution. The difference between the analytical and numerical results for the neutron flux is significantly increased in the edge of the problem geometry, even though the maximum difference is lower than 4%. This phenomenon seems to arise from the derivative boundary condition at (x,0) and (0,y) positions, and it may be necessary to introduce additional strategy (e.g., the method using fictitious points and

Errors of the backextrapolation method in determination of the blood volume

Science.gov (United States)

Schröder, T.; Rösler, U.; Frerichs, I.; Hahn, G.; Ennker, J.; Hellige, G.

1999-01-01

Backextrapolation is an empirical method to calculate the central volume of distribution (for example the blood volume). It is based on the compartment model, which says that after an injection the substance is distributed instantaneously in the central volume with no time delay. The occurrence of recirculation is not taken into account. The change of concentration with time of indocyanine green (ICG) was observed in an in vitro model, in which the volume was recirculating in 60 s and the clearance of the ICG could be varied. It was found that the higher the elimination of ICG, the higher was the error of the backextrapolation method. The theoretical consideration of Schröder et al ( Biomed. Tech. 42 (1997) 7-11) was proved. If the injected substance is eliminated somewhere in the body (i.e. not by radioactive decay), the backextrapolation method produces large errors.

Appraisal of numerical methods in predicting the aerodynamics of forward-swept wings

CSIR Research Space (South Africa)

Lombardi, G

1998-07-01

Full Text Available The capabilities of different numerical methods in evaluating the aerodynamic characteristics of a forward-swept wing in subsonic and transonic now are analyzed. The numerical results, obtained by means of potential, Euler, and Navier-Stokes solvers...

Nonlinear bioheat transfer models and multi-objective numerical optimization of the cryosurgery operations

Energy Technology Data Exchange (ETDEWEB)

Kudryashov, Nikolay A.; Shilnikov, Kirill E. [National Research Nuclear University MEPhI, Department of Applied Mathematics, Moscow (Russian Federation)

2016-06-08

Numerical computation of the three dimensional problem of the freezing interface propagation during the cryosurgery coupled with the multi-objective optimization methods is used in order to improve the efficiency and safety of the cryosurgery operations performing. Prostate cancer treatment and cutaneous cryosurgery are considered. The heat transfer in soft tissue during the thermal exposure to low temperature is described by the Pennes bioheat model and is coupled with an enthalpy method for blurred phase change computations. The finite volume method combined with the control volume approximation of the heat fluxes is applied for the cryosurgery numerical modeling on the tumor tissue of a quite arbitrary shape. The flux relaxation approach is used for the stability improvement of the explicit finite difference schemes. The method of the additional heating elements mounting is studied as an approach to control the cellular necrosis front propagation. Whereas the undestucted tumor tissue and destucted healthy tissue volumes are considered as objective functions, the locations of additional heating elements in cutaneous cryosurgery and cryotips in prostate cancer cryotreatment are considered as objective variables in multi-objective problem. The quasi-gradient method is proposed for the searching of the Pareto front segments as the multi-objective optimization problem solutions.

Method of phase-Doppler anemometry free from the measurement-volume effect.

Science.gov (United States)

Qiu, H; Hsu, C T

1999-05-01

A novel method is developed to improve the accuracy of particle sizing in laser phase-Doppler anemometry (PDA). In this method the vector sum of refractive and reflective rays is taken into consideration in describing a dual-mechanism-scattering model caused by a nonuniformly illuminated PDA measurement volume. The constraint of the single-mechanism-scattering model in the conventional PDA is removed. As a result the error caused by the measurement-volume effect, which consists of a Gaussian-beam defect and a slit effect, can be eliminated. This new method can be easily implemented with minimal modification of the conventional PDA system. The results of simulation based on the generalized Lorenz-Mie theory show that the new method can provide a PDA system free from the measurement-volume effect.

Numerical simulation of the regularized long wave equation by He's homotopy perturbation method

International Nuclear Information System (INIS)

Inc, Mustafa; Ugurlu, Yavuz

2007-01-01

In this Letter, we present the homotopy perturbation method (shortly HPM) for obtaining the numerical solution of the RLW equation. We obtain the exact and numerical solutions of the Regularized Long Wave (RLW) equation for certain initial condition. The initial approximation can be freely chosen with possible unknown constants which can be determined by imposing the boundary and initial conditions. Comparison of the results with those of other methods have led us to significant consequences. The numerical solutions are compared with the known analytical solutions

Numerical methods in dynamic fracture mechanics

International Nuclear Information System (INIS)

Beskos, D.E.

1987-01-01

A review of numerical methods for the solution of dynamic problems of fracture mechanics is presented. Finite difference, finite element and boundary element methods as applied to linear elastic or viscoelastic and non-linear elastoplastic or elastoviscoplastic dynamic fracture mechanics problems are described and critically evaluated. Both cases of stationary cracks and rapidly propagating cracks of simple I, II, III or mixed modes are considered. Harmonically varying with time or general transient dynamic disturbances in the form of external loading or incident waves are taken into account. Determination of the dynamic stress intensity factor for stationary cracks or moving cracks with known velocity history as well as determination of the crack-tip propagation history for given dynamic fracture toughness versus crack velocity relation are described and illustrated by means of certain representative examples. Finally, a brief assessment of the present state of knowledge is made and research needs are identified

A Three-Dimensional, Immersed Boundary, Finite Volume Method for the Simulation of Incompressible Heat Transfer Flows around Complex Geometries

Directory of Open Access Journals (Sweden)

Hassan Badreddine

2017-01-01

Full Text Available The current work focuses on the development and application of a new finite volume immersed boundary method (IBM to simulate three-dimensional fluid flows and heat transfer around complex geometries. First, the discretization of the governing equations based on the second-order finite volume method on Cartesian, structured, staggered grid is outlined, followed by the description of modifications which have to be applied to the discretized system once a body is immersed into the grid. To validate the new approach, the heat conduction equation with a source term is solved inside a cavity with an immersed body. The approach is then tested for a natural convection flow in a square cavity with and without circular cylinder for different Rayleigh numbers. The results computed with the present approach compare very well with the benchmark solutions. As a next step in the validation procedure, the method is tested for Direct Numerical Simulation (DNS of a turbulent flow around a surface-mounted matrix of cubes. The results computed with the present method compare very well with Laser Doppler Anemometry (LDA measurements of the same case, showing that the method can be used for scale-resolving simulations of turbulence as well.

Calculating the tumor volume of acoustic neuromas: comparison of ABC/2 formula with planimetry method.

Science.gov (United States)

Yu, Yi-Lin; Lee, Meei-Shyuan; Juan, Chun-Jung; Hueng, Dueng-Yuan

2013-08-01

The ABC/2 equation is commonly applied to measure the volume of intracranial hematoma. However, the precision of ABC/2 equation in estimating the tumor volume of acoustic neuromas is less addressed. The study is to evaluate the accuracy of the ABC/2 formula by comparing with planimetry method for estimating the tumor volumes. Thirty-two patients diagnosed with acoustic neuroma received contrast-enhanced magnetic resonance imaging of brain were recruited. The volume was calculated by the ABC/2 equation and planimetry method (defined as exact volume) at the same time. The 32 patients were divided into three groups by tumor volume to avoid volume-dependent overestimation (6 ml). The tumor volume by ABC/2 method was highly correlated to that calculated by planimetry method using linear regression analysis (R2=0.985). Pearson correlation coefficient (r=0.993, pABC/2 formula is an easy method in estimating the tumor volume of acoustic neuromas that is not inferior to planimetry method. Copyright © 2013 Elsevier B.V. All rights reserved.

Finite elements volumes methods: applications to the Navier-Stokes equations and convergence results

International Nuclear Information System (INIS)

Emonot, P.

1992-01-01

In the first chapter are described the equations modeling incompressible fluid flow and a quick presentation of finite volumes method. The second chapter is an introduction to the finite elements volumes method. The box model is described and a method adapted to Navier-Stokes problems is proposed. The third chapter shows a fault analysis of the finite elements volumes method for the Laplacian problem and some examples in one, two, three dimensional calculations. The fourth chapter is an extension of the error analysis of the method for the Navier-Stokes problem

Efficient numerical method for district heating system hydraulics

International Nuclear Information System (INIS)

Stevanovic, Vladimir D.; Prica, Sanja; Maslovaric, Blazenka; Zivkovic, Branislav; Nikodijevic, Srdjan

2007-01-01

An efficient method for numerical simulation and analyses of the steady state hydraulics of complex pipeline networks is presented. It is based on the loop model of the network and the method of square roots for solving the system of linear equations. The procedure is presented in the comprehensive mathematical form that could be straightforwardly programmed into a computer code. An application of the method to energy efficiency analyses of a real complex district heating system is demonstrated. The obtained results show a potential for electricity savings in pumps operation. It is shown that the method is considerably more effective than the standard Hardy Cross method still widely used in engineering practice. Because of the ease of implementation and high efficiency, the method presented in this paper is recommended for hydraulic steady state calculations of complex networks

NUMERICAL METHODS FOR SOLVING THE MULTI-TERM TIME-FRACTIONAL WAVE-DIFFUSION EQUATION

OpenAIRE

Liu, F.; Meerschaert, M.M.; McGough, R.J.; Zhuang, P.; Liu, Q.

2013-01-01

In this paper, the multi-term time-fractional wave-diffusion equations are considered. The multi-term time fractional derivatives are defined in the Caputo sense, whose orders belong to the intervals [0,1], [1,2), [0,2), [0,3), [2,3) and [2,4), respectively. Some computationally effective numerical methods are proposed for simulating the multi-term time-fractional wave-diffusion equations. The numerical results demonstrate the effectiveness of theoretical analysis. These methods and technique...

A different approach to estimate nonlinear regression model using numerical methods

Science.gov (United States)

Mahaboob, B.; Venkateswarlu, B.; Mokeshrayalu, G.; Balasiddamuni, P.

2017-11-01

This research paper concerns with the computational methods namely the Gauss-Newton method, Gradient algorithm methods (Newton-Raphson method, Steepest Descent or Steepest Ascent algorithm method, the Method of Scoring, the Method of Quadratic Hill-Climbing) based on numerical analysis to estimate parameters of nonlinear regression model in a very different way. Principles of matrix calculus have been used to discuss the Gradient-Algorithm methods. Yonathan Bard [1] discussed a comparison of gradient methods for the solution of nonlinear parameter estimation problems. However this article discusses an analytical approach to the gradient algorithm methods in a different way. This paper describes a new iterative technique namely Gauss-Newton method which differs from the iterative technique proposed by Gorden K. Smyth [2]. Hans Georg Bock et.al [10] proposed numerical methods for parameter estimation in DAE’s (Differential algebraic equation). Isabel Reis Dos Santos et al [11], Introduced weighted least squares procedure for estimating the unknown parameters of a nonlinear regression metamodel. For large-scale non smooth convex minimization the Hager and Zhang (HZ) conjugate gradient Method and the modified HZ (MHZ) method were presented by Gonglin Yuan et al [12].

A convergent numerical method for the Navier-Stokes-Fourier system

Czech Academy of Sciences Publication Activity Database

Feireisl, Eduard; Karper, T.; Novotný, A.

2016-01-01

Roč. 36, č. 4 (2016), s. 1477-1535 ISSN 0272-4979 EU Projects: European Commission(XE) 320078 - MATHEF Institutional support: RVO:67985840 Keywords : Navier-Stokes-Fourier system * Crouzeix-Raviart finite element method * finite volume method Subject RIV: BA - General Mathematics Impact factor: 1.703, year: 2016 http://imajna.oxfordjournals.org/content/36/4/1477

New numerical method for solving the solute transport equation

International Nuclear Information System (INIS)

Ross, B.; Koplik, C.M.

1978-01-01

The solute transport equation can be solved numerically by approximating the water flow field by a network of stream tubes and using a Green's function solution within each stream tube. Compared to previous methods, this approach permits greater computational efficiency and easier representation of small discontinuities, and the results are easier to interpret physically. The method has been used to study hypothetical sites for disposal of high-level radioactive waste

International Winter Workshop on Differential Equations and Numerical Analysis

CERN Document Server

Miller, John; Narasimhan, Ramanujam; Mathiazhagan, Paramasivam; Victor, Franklin

2016-01-01

This book offers an ideal introduction to singular perturbation problems, and a valuable guide for researchers in the field of differential equations. It also includes chapters on new contributions to both fields: differential equations and singular perturbation problems. Written by experts who are active researchers in the related fields, the book serves as a comprehensive source of information on the underlying ideas in the construction of numerical methods to address different classes of problems with solutions of different behaviors, which will ultimately help researchers to design and assess numerical methods for solving new problems. All the chapters presented in the volume are complemented by illustrations in the form of tables and graphs.

Approximate Analytic and Numerical Solutions to Lane-Emden Equation via Fuzzy Modeling Method

Directory of Open Access Journals (Sweden)

De-Gang Wang

2012-01-01

Full Text Available A novel algorithm, called variable weight fuzzy marginal linearization (VWFML method, is proposed. This method can supply approximate analytic and numerical solutions to Lane-Emden equations. And it is easy to be implemented and extended for solving other nonlinear differential equations. Numerical examples are included to demonstrate the validity and applicability of the developed technique.

6th international symposium on finite volumes for complex applications

CERN Document Server

Halama, Jan; Herbin, Raphaèle; Hubert, Florence; Fort, Jaroslav; FVCA 6; Finite Volumes for Complex Applications VI : Problems and perspectives

2011-01-01

Finite volume methods are used for various applications in fluid dynamics, magnetohydrodynamics, structural analysis or nuclear physics. A closer look reveals many interesting phenomena and mathematical or numerical difficulties, such as true error analysis and adaptivity, modelling of multi-phase phenomena or fitting problems, stiff terms in convection/diffusion equations and sources. To overcome existing problems and to find solution methods for future applications requires many efforts and always new developments. The goal of The International Symposium on Finite Volumes for Complex Applica

Development of CAD implementing the algorithm of boundary elements’ numerical analytical method

Directory of Open Access Journals (Sweden)

Yulia V. Korniyenko

2015-03-01

Full Text Available Up to recent days the algorithms for numerical-analytical boundary elements method had been implemented with programs written in MATLAB environment language. Each program had a local character, i.e. used to solve a particular problem: calculation of beam, frame, arch, etc. Constructing matrices in these programs was carried out “manually” therefore being time-consuming. The research was purposed onto a reasoned choice of programming language for new CAD development, allows to implement algorithm of numerical analytical boundary elements method and to create visualization tools for initial objects and calculation results. Research conducted shows that among wide variety of programming languages the most efficient one for CAD development, employing the numerical analytical boundary elements method algorithm, is the Java language. This language provides tools not only for development of calculating CAD part, but also to build the graphic interface for geometrical models construction and calculated results interpretation.

Fictitious domain methods for elliptic problems with general boundary conditions with an application to the numerical simulation of two phase flows

International Nuclear Information System (INIS)

Ramiere, I.

2006-09-01

This work is dedicated to the introduction of two original fictitious domain methods for the resolution of elliptic problems (mainly convection-diffusion problems) with general and eventually mixed boundary conditions: Dirichlet, Robin or Neumann. The originality lies in the approximation of the immersed boundary by an approximate interface derived from the fictitious domain Cartesian mesh, which is generally not boundary-fitted to the physical domain. The same generic numerical scheme is used to impose the embedded boundary conditions. Hence, these methods require neither a surface mesh of the immersed boundary nor the local modification of the numerical scheme. We study two modelling of the immersed boundary. In the first one, called spread interface, the approximate immersed boundary is the union of the cells crossed by the physical immersed boundary. In the second one, called thin interface, the approximate immersed boundary lies on sides of mesh cells. Additional algebraic transmission conditions linking both flux and solution jumps through the thin approximate interface are introduced. The fictitious problem to solve as well as the treatment of the embedded boundary conditions are detailed for the two methods. A Q1 finite element scheme is implemented for the numerical validation of the spread interface approach while a new cell-centered finite volume scheme is derived for the thin interface approach with immersed jumps. Each method is then combined to multilevel local mesh refinement algorithms (with solution or flux residual) to increase the precision of the solution in the vicinity of the immersed interface. A convergence analysis of a Q1 finite element method with non-boundary fitted meshes is also presented. This study proves the convergence rates of the present methods. Among the various industrial applications, the simulation on a model of heat exchanger in french nuclear power plants enables us to appreciate the performances of the fictitious domain

A method for bubble volume calculating in vertical two-phase flow

International Nuclear Information System (INIS)

Wang, H Y; Dong, F

2009-01-01

The movement of bubble is a basic subject in gas-liquid two-phase flow research. A method for calculating bubble volume which is one of the most important characters in bubble motion research was proposed. A suit of visualized experimental device was designed and set up. Single bubble rising in stagnant liquid in a rectangular tank was studied using the high-speed video system. Bubbles generated by four orifice with different diameter (1mm, 2mm, 3mm, 4mm) were recorded respectively. Sequences of recorded high-speed images were processed by digital image processing method, such as image noise remove, binary image transform, bubble filling, and so on. then, Several parameters could be obtained from the processed image. Bubble area, equivalent diameter, bubble velocity, bubble acceleration are all indispensable in bubble volume calculating. In order to get the force balance equation, forces that work on bubble along vertical direction, including drag force, virtual mass force, buoyancy, gravity and liquid thrust, were analyzed. Finally, the bubble volume formula could be derived from the force balance equation and bubble parameters. Examples were given to shown the computing process and results. Comparison of the bubble volume calculated by geomettic method and the present method have shown the superiority of the proposed method in this paper.

Hybrid Multiscale Finite Volume method for multiresolution simulations of flow and reactive transport in porous media

Science.gov (United States)

Barajas-Solano, D. A.; Tartakovsky, A. M.

2017-12-01

We present a multiresolution method for the numerical simulation of flow and reactive transport in porous, heterogeneous media, based on the hybrid Multiscale Finite Volume (h-MsFV) algorithm. The h-MsFV algorithm allows us to couple high-resolution (fine scale) flow and transport models with lower resolution (coarse) models to locally refine both spatial resolution and transport models. The fine scale problem is decomposed into various "local'' problems solved independently in parallel and coordinated via a "global'' problem. This global problem is then coupled with the coarse model to strictly ensure domain-wide coarse-scale mass conservation. The proposed method provides an alternative to adaptive mesh refinement (AMR), due to its capacity to rapidly refine spatial resolution beyond what's possible with state-of-the-art AMR techniques, and the capability to locally swap transport models. We illustrate our method by applying it to groundwater flow and reactive transport of multiple species.

A volume of fluid method based on multidimensional advection and spline interface reconstruction

International Nuclear Information System (INIS)

Lopez, J.; Hernandez, J.; Gomez, P.; Faura, F.

2004-01-01

A new volume of fluid method for tracking two-dimensional interfaces is presented. The method involves a multidimensional advection algorithm based on the use of edge-matched flux polygons to integrate the volume fraction evolution equation, and a spline-based reconstruction algorithm. The accuracy and efficiency of the proposed method are analyzed using different tests, and the results are compared with those obtained recently by other authors. Despite its simplicity, the proposed method represents a significant improvement, and compares favorably with other volume of fluid methods as regards the accuracy and efficiency of both the advection and reconstruction steps

Numerical method for wave forces acting on partially perforated caisson

Science.gov (United States)

Jiang, Feng; Tang, Xiao-cheng; Jin, Zhao; Zhang, Li; Chen, Hong-zhou

2015-04-01

The perforated caisson is widely applied to practical engineering because of its great advantages in effectively wave energy consumption and cost reduction. The attentions of many scientists were paid to the fluid-structure interaction between wave and perforated caisson studies, but until now, most concerns have been put on theoretical analysis and experimental model set up. In this paper, interaction between the wave and the partial perforated caisson in a 2D numerical wave flume is investigated by means of the renewed SPH algorithm, and the mathematical equations are in the form of SPH numerical approximation based on Navier-Stokes equations. The validity of the SPH mathematical method is examined and the simulated results are compared with the results of theoretical models, meanwhile the complex hydrodynamic characteristics when the water particles flow in or out of a wave absorbing chamber are analyzed and the wave pressure distribution of the perforated caisson is also addressed here. The relationship between the ratio of total horizontal force acting on caisson under regular waves and its influence factors is examined. The data show that the numerical calculation of the ratio of total horizontal force meets the empirical regression equation very well. The simulations of SPH about the wave nonlinearity and breaking are briefly depicted in the paper, suggesting that the advantages and great potentiality of the SPH method is significant compared with traditional methods.

An energy stable evolution method for simulating two-phase equilibria of multi-component fluids at constant moles, volume and temperature

KAUST Repository

Kou, Jisheng

2016-02-25

In this paper, we propose an energy-stable evolution method for the calculation of the phase equilibria under given volume, temperature, and moles (VT-flash). An evolution model for describing the dynamics of two-phase fluid system is based on Fick’s law of diffusion for multi-component fluids and the Peng-Robinson equation of state. The mobility is obtained from diffusion coefficients by relating the gradient of chemical potential to the gradient of molar density. The evolution equation for moles of each component is derived using the discretization of diffusion equations, while the volume evolution equation is constructed based on the mechanical mechanism and the Peng-Robinson equation of state. It is proven that the proposed evolution system can well model the VT-flash problem, and moreover, it possesses the property of total energy decay. By using the Euler time scheme to discretize this evolution system, we develop an energy stable algorithm with an adaptive choice strategy of time steps, which allows us to calculate the suitable time step size to guarantee the physical properties of moles and volumes, including positivity, maximum limits, and correct definition of the Helmhotz free energy function. The proposed evolution method is also proven to be energy-stable under the proposed time step choice. Numerical examples are tested to demonstrate efficiency and robustness of the proposed method.

Numerical Solution of the Blasius Viscous Flow Problem by Quartic B-Spline Method

Directory of Open Access Journals (Sweden)

Hossein Aminikhah

2016-01-01

Full Text Available A numerical method is proposed to study the laminar boundary layer about a flat plate in a uniform stream of fluid. The presented method is based on the quartic B-spline approximations with minimizing the error L2-norm. Theoretical considerations are discussed. The computed results are compared with some numerical results to show the efficiency of the proposed approach.

Steady-state transport equation resolution by particle methods, and numerical results

International Nuclear Information System (INIS)

Mercier, B.

1985-10-01

A method to solve steady-state transport equation has been given. Principles of the method are given. The method is studied in two different cases; estimations given by the theory are compared to numerical results. Results got in 1-D (spherical geometry) and in 2-D (axisymmetric geometry) are given [fr

Adaptive moving grid methods for two-phase flow in porous media

KAUST Repository

Dong, Hao; Qiao, Zhonghua; Sun, Shuyu; Tang, Tao

2014-01-01

In this paper, we present an application of the moving mesh method for approximating numerical solutions of the two-phase flow model in porous media. The numerical schemes combine a mixed finite element method and a finite volume method, which can

Numerical method for solving the three-dimensional time-dependent neutron diffusion equation

International Nuclear Information System (INIS)

Khaled, S.M.; Szatmary, Z.

2005-01-01

A numerical time-implicit method has been developed for solving the coupled three-dimensional time-dependent multi-group neutron diffusion and delayed neutron precursor equations. The numerical stability of the implicit computation scheme and the convergence of the iterative associated processes have been evaluated. The computational scheme requires the solution of large linear systems at each time step. For this purpose, the point over-relaxation Gauss-Seidel method was chosen. A new scheme was introduced instead of the usual source iteration scheme. (author)

Two split cell numerical methods for solving 2-D non-equilibrium radiation transport equations

International Nuclear Information System (INIS)

Feng Tinggui

2004-11-01

Two numerically positive methods, the step characteristic integral method and subcell balance method, for solving radiative transfer equations on quadrilateral grids are presented. Numerical examples shows that the schemes presented are feasible on non-rectangle grid computation, and that the computing results by the schemes presented are comparative to that by the discrete ordinate diamond scheme on rectangle grid. (author)

Lung lesion doubling times: values and variability based on method of volume determination

International Nuclear Information System (INIS)

Eisenbud Quint, Leslie; Cheng, Joan; Schipper, Matthew; Chang, Andrew C.; Kalemkerian, Gregory

2008-01-01

Purpose: To determine doubling times (DTs) of lung lesions based on volumetric measurements from thin-section CT imaging. Methods: Previously untreated patients with ≥ two thin-section CT scans showing a focal lung lesion were identified. Lesion volumes were derived using direct volume measurements and volume calculations based on lesion area and diameter. Growth rates (GRs) were compared by tissue diagnosis and measurement technique. Results: 54 lesions were evaluated including 8 benign lesions, 10 metastases, 3 lymphomas, 15 adenocarcinomas, 11 squamous carcinomas, and 7 miscellaneous lung cancers. Using direct volume measurements, median DTs were 453, 111, 15, 181, 139 and 137 days, respectively. Lung cancer DTs ranged from 23-2239 days. There were no significant differences in GRs among the different lesion types. There was considerable variability among GRs using different volume determination methods. Conclusions: Lung cancer doubling times showed a substantial range, and different volume determination methods gave considerably different DTs

Application of numerical inverse method in calculation of composition-dependent interdiffusion coefficients in finite diffusion couples

DEFF Research Database (Denmark)

Liu, Yuanrong; Chen, Weimin; Zhong, Jing

2017-01-01

The previously developed numerical inverse method was applied to determine the composition-dependent interdiffusion coefficients in single-phase finite diffusion couples. The numerical inverse method was first validated in a fictitious binary finite diffusion couple by pre-assuming four standard...... sets of interdiffusion coefficients. After that, the numerical inverse method was then adopted in a ternary Al-Cu-Ni finite diffusion couple. Based on the measured composition profiles, the ternary interdiffusion coefficients along the entire diffusion path of the target ternary diffusion couple were...... obtained by using the numerical inverse approach. The comprehensive comparisons between the computations and the experiments indicate that the numerical inverse method is also applicable to high-throughput determination of the composition-dependent interdiffusion coefficients in finite diffusion couples....

Confining dyon gas with finite-volume effects under control

Energy Technology Data Exchange (ETDEWEB)

Bruckmann, Falk [Regensburg Univ. (Germany). Institut fuer Theoretische Physik; Dinter, Simon [Deutsches Elektronen-Synchrotron (DESY), Zeuthen (Germany). John von Neumann-Inst. fuer Computing NIC; Humboldt-Universitaet, Berlin (Germany). Inst. fuer Physik; Ilgenfritz, Ernst-Michael [Humboldt-Universitaet, Berlin (Germany). Inst. fuer Physik; Joint Institute for Nuclear Research, VBLHEP, Dubna (Russian Federation); Maier, Benjamin; Mueller-Preussker, Michael [Humboldt-Universitaet, Berlin (Germany). Inst. fuer Physik; Wagner, Marc [Humboldt-Universitaet, Berlin (Germany). Inst. fuer Physik; Frankfurt Univ. (Germany). Inst. fuer Theoretische Physik

2011-11-15

As an approach to describe the long-range properties of non-Abelian gauge theories at non-zero temperature Tnumerical treatments, the long-range tails of the dyon fields cause severe finite-volume effects. Therefore, we demonstrate the application of Ewald's summation method to this system. Finite-volume effects are shown to be under control, which is a crucial requirement for numerical studies of interacting dyon ensembles. (orig.)

Confining dyon gas with finite-volume effects under control

International Nuclear Information System (INIS)

Bruckmann, Falk; Maier, Benjamin; Mueller-Preussker, Michael; Wagner, Marc; Frankfurt Univ.

2011-11-01

As an approach to describe the long-range properties of non-Abelian gauge theories at non-zero temperature T c , we consider a non-interacting ensemble of dyons (magnetic monopoles) with non-trivial holonomy. We show analytically, that the quark-antiquark free energy from the Polyakov loop correlator grows linearly with the distance, and how the string tension scales with the dyon density. In numerical treatments, the long-range tails of the dyon fields cause severe finite-volume effects. Therefore, we demonstrate the application of Ewald's summation method to this system. Finite-volume effects are shown to be under control, which is a crucial requirement for numerical studies of interacting dyon ensembles. (orig.)

Numerical Method for Darcy Flow Derived Using Discrete Exterior Calculus

Science.gov (United States)

Hirani, A. N.; Nakshatrala, K. B.; Chaudhry, J. H.

2015-05-01

We derive a numerical method for Darcy flow, and also for Poisson's equation in mixed (first order) form, based on discrete exterior calculus (DEC). Exterior calculus is a generalization of vector calculus to smooth manifolds and DEC is one of its discretizations on simplicial complexes such as triangle and tetrahedral meshes. DEC is a coordinate invariant discretization, in that it does not depend on the embedding of the simplices or the whole mesh. We start by rewriting the governing equations of Darcy flow using the language of exterior calculus. This yields a formulation in terms of flux differential form and pressure. The numerical method is then derived by using the framework provided by DEC for discretizing differential forms and operators that act on forms. We also develop a discretization for a spatially dependent Hodge star that varies with the permeability of the medium. This also allows us to address discontinuous permeability. The matrix representation for our discrete non-homogeneous Hodge star is diagonal, with positive diagonal entries. The resulting linear system of equations for flux and pressure are saddle type, with a diagonal matrix as the top left block. The performance of the proposed numerical method is illustrated on many standard test problems. These include patch tests in two and three dimensions, comparison with analytically known solutions in two dimensions, layered medium with alternating permeability values, and a test with a change in permeability along the flow direction. We also show numerical evidence of convergence of the flux and the pressure. A convergence experiment is included for Darcy flow on a surface. A short introduction to the relevant parts of smooth and discrete exterior calculus is included in this article. We also include a discussion of the boundary condition in terms of exterior calculus.

Development of a set of benchmark problems to verify numerical methods for solving burnup equations

International Nuclear Information System (INIS)

Lago, Daniel; Rahnema, Farzad

2017-01-01

Highlights: • Description transmutation chain benchmark problems. • Problems for validating numerical methods for solving burnup equations. • Analytical solutions for the burnup equations. • Numerical solutions for the burnup equations. - Abstract: A comprehensive set of transmutation chain benchmark problems for numerically validating methods for solving burnup equations was created. These benchmark problems were designed to challenge both traditional and modern numerical methods used to solve the complex set of ordinary differential equations used for tracking the change in nuclide concentrations over time due to nuclear phenomena. Given the development of most burnup solvers is done for the purpose of coupling with an established transport solution method, these problems provide a useful resource in testing and validating the burnup equation solver before coupling for use in a lattice or core depletion code. All the relevant parameters for each benchmark problem are described. Results are also provided in the form of reference solutions generated by the Mathematica tool, as well as additional numerical results from MATLAB.

Gastrointestinal tract volume measurement method using a compound eye type endoscope

Science.gov (United States)

Yoshimoto, Kayo; Yamada, Kenji; Watabe, Kenji; Kido, Michiko; Nagakura, Toshiaki; Takahashi, Hideya; Nishida, Tsutomu; Iijima, Hideki; Tsujii, Masahiko; Takehara, Tetsuo; Ohno, Yuko

2015-03-01

We propose an intestine volume measurement method using a compound eye type endoscope. This method aims at assessment of the gastrointestinal function. Gastrointestinal diseases are mainly based on morphological abnormalities. However, gastrointestinal symptoms are sometimes apparent without visible abnormalities. Such diseases are called functional gastrointestinal disorder, for example, functional dyspepsia, and irritable bowel syndrome. One of the major factors for these diseases is abnormal gastrointestinal motility. For the diagnosis of the gastrointestinal tract, both aspects of organic and functional assessment is important. While endoscopic diagnosis is essential for assessment of organic abnormalities, three-dimensional information is required for assessment of the functional abnormalities. Thus, we proposed the three dimensional endoscope system using compound eye. In this study, we forces on the volume of gastrointestinal tract. The volume of the gastrointestinal tract is thought to related its function. In our system, we use a compound eye type endoscope system to obtain three-dimensional information of the tract. The volume can be calculated by integrating the slice data of the intestine tract shape using the obtained three-dimensional information. First, we evaluate the proposed method by known-shape tube. Then, we confirm that the proposed method can measure the tract volume using the tract simulated model. Our system can assess the wall of gastrointestinal tract directly in a three-dimensional manner. Our system can be used for examination of gastric morphological and functional abnormalities.

Annual review of numerical fluid mechanics and heat transfer. Volume 1

International Nuclear Information System (INIS)

Chawla, T.C.

1987-01-01

Numerical techniqes for the analysis of problems in fluid mechanics and heat transfer are discussed, reviewing the results of recent investigations. Topics addressed include thermal radiation in particulate media with dependent and independent scattering, pressure-velocity coupling in incompressiblefluid flow, new explicit methods for diffusion problems, and one-dimensional reaction-diffusion equations in combustion theory. Consideration is given to buckling flows, multidimensional radiative-transfer analysis in participating media, freezing and melting problems, and complex heat-transfer processes in heat-generating horizontal fluid layers

Numerical simulation of microdroplet dynamics in microfluidics using finite element and level set methods

CSIR Research Space (South Africa)

Mbanjwa, MB

2012-10-01

Full Text Available the dispersed phase (droplets). Miniaturised flow systems are predominately dominated by surface forces due to the Scaling Law (Figure 1). Figure 1: Surface forces dominate over volume forces in microsystems (image credit: DreamWorks[3]) Table 1: Important...?1065. 3. Antz Movie. 1998. DreamWorks Animation SKG Inc. 4. Bruus, H. 2008. Theoretical microfluidics. Oxford University Press, Oxford. K-10032 [www.kashan.co.za] Alongside experimental work, numerical tools, such as computational fl uid dynamics...

Flow studies in canine artery bifurcations using a numerical simulation method.

Science.gov (United States)

Xu, X Y; Collins, M W; Jones, C J

1992-11-01

Three-dimensional flows through canine femoral bifurcation models were predicted under physiological flow conditions by solving numerically the time-dependent three-dimensional Navier-stokes equations. In the calculations, two models were assumed for the blood, those of (a) a Newtonian fluid, and (b) a non-Newtonian fluid obeying the power law. The blood vessel wall was assumed to be rigid this being the only approximation to the prediction model. The numerical procedure utilized a finite volume approach on a finite element mesh to discretize the equations, and the code used (ASTEC) incorporated the SIMPLE velocity-pressure algorithm in performing the calculations. The predicted velocity profiles were in good qualitative agreement with the in vivo measurements recently obtained by Jones et al. The non-Newtonian effects on the bifurcation flow field were also investigated, and no great differences in velocity profiles were observed. This indicated that the non-Newtonian characteristics of the blood might not be an important factor in determining the general flow patterns for these bifurcations, but could have local significance. Current work involves modeling wall distensibility in an empirically valid manner. Predictions accommodating these will permit a true quantitative comparison with experiment.

The discrete ordinate method in association with the finite-volume method in non-structured mesh; Methode des ordonnees discretes associee a la methode des volumes finis en maillage non structure

Energy Technology Data Exchange (ETDEWEB)

Le Dez, V; Lallemand, M [Ecole Nationale Superieure de Mecanique et d` Aerotechnique (ENSMA), 86 - Poitiers (France); Sakami, M; Charette, A [Quebec Univ., Chicoutimi, PQ (Canada). Dept. des Sciences Appliquees

1997-12-31

The description of an efficient method of radiant heat transfer field determination in a grey semi-transparent environment included in a 2-D polygonal cavity with surface boundaries that reflect the radiation in a purely diffusive manner is proposed, at the equilibrium and in radiation-conduction coupling situation. The technique uses simultaneously the finite-volume method in non-structured triangular mesh, the discrete ordinate method and the ray shooting method. The main mathematical developments and comparative results with the discrete ordinate method in orthogonal curvilinear coordinates are included. (J.S.) 10 refs.

The discrete ordinate method in association with the finite-volume method in non-structured mesh; Methode des ordonnees discretes associee a la methode des volumes finis en maillage non structure

Energy Technology Data Exchange (ETDEWEB)

Le Dez, V.; Lallemand, M. [Ecole Nationale Superieure de Mecanique et d`Aerotechnique (ENSMA), 86 - Poitiers (France); Sakami, M.; Charette, A. [Quebec Univ., Chicoutimi, PQ (Canada). Dept. des Sciences Appliquees

1996-12-31

The description of an efficient method of radiant heat transfer field determination in a grey semi-transparent environment included in a 2-D polygonal cavity with surface boundaries that reflect the radiation in a purely diffusive manner is proposed, at the equilibrium and in radiation-conduction coupling situation. The technique uses simultaneously the finite-volume method in non-structured triangular mesh, the discrete ordinate method and the ray shooting method. The main mathematical developments and comparative results with the discrete ordinate method in orthogonal curvilinear coordinates are included. (J.S.) 10 refs.

ATHENA code manual. Volume 1. Code structure, system models, and solution methods

International Nuclear Information System (INIS)

Carlson, K.E.; Roth, P.A.; Ransom, V.H.

1986-09-01

The ATHENA (Advanced Thermal Hydraulic Energy Network Analyzer) code has been developed to perform transient simulation of the thermal hydraulic systems which may be found in fusion reactors, space reactors, and other advanced systems. A generic modeling approach is utilized which permits as much of a particular system to be modeled as necessary. Control system and secondary system components are included to permit modeling of a complete facility. Several working fluids are available to be used in one or more interacting loops. Different loops may have different fluids with thermal connections between loops. The modeling theory and associated numerical schemes are documented in Volume I in order to acquaint the user with the modeling base and thus aid effective use of the code. The second volume contains detailed instructions for input data preparation

Research on volume metrology method of large vertical energy storage tank based on internal electro-optical distance-ranging method

Science.gov (United States)

Hao, Huadong; Shi, Haolei; Yi, Pengju; Liu, Ying; Li, Cunjun; Li, Shuguang

2018-01-01

A Volume Metrology method based on Internal Electro-optical Distance-ranging method is established for large vertical energy storage tank. After analyzing the vertical tank volume calculation mathematical model, the key processing algorithms, such as gross error elimination, filtering, streamline, and radius calculation are studied for the point cloud data. The corresponding volume values are automatically calculated in the different liquids by calculating the cross-sectional area along the horizontal direction and integrating from vertical direction. To design the comparison system, a vertical tank which the nominal capacity is 20,000 m3 is selected as the research object, and there are shown that the method has good repeatability and reproducibility. Through using the conventional capacity measurement method as reference, the relative deviation of calculated volume is less than 0.1%, meeting the measurement requirements. And the feasibility and effectiveness are demonstrated.

Numerical Simulation of Partially-Coherent Broadband Optical Imaging Using the FDTD Method

Science.gov (United States)

Çapoğlu, İlker R.; White, Craig A.; Rogers, Jeremy D.; Subramanian, Hariharan; Taflove, Allen; Backman, Vadim

2012-01-01

Rigorous numerical modeling of optical systems has attracted interest in diverse research areas ranging from biophotonics to photolithography. We report the full-vector electromagnetic numerical simulation of a broadband optical imaging system with partially-coherent and unpolarized illumination. The scattering of light from the sample is calculated using the finite-difference time-domain (FDTD) numerical method. Geometrical optics principles are applied to the scattered light to obtain the intensity distribution at the image plane. Multilayered object spaces are also supported by our algorithm. For the first time, numerical FDTD calculations are directly compared to and shown to agree well with broadband experimental microscopy results. PMID:21540939

Efficient Numerical Methods for Stochastic Differential Equations in Computational Finance

KAUST Repository

Happola, Juho

2017-09-19

Stochastic Differential Equations (SDE) offer a rich framework to model the probabilistic evolution of the state of a system. Numerical approximation methods are typically needed in evaluating relevant Quantities of Interest arising from such models. In this dissertation, we present novel effective methods for evaluating Quantities of Interest relevant to computational finance when the state of the system is described by an SDE.

Efficient Numerical Methods for Stochastic Differential Equations in Computational Finance

KAUST Repository

Happola, Juho

2017-01-01

Stochastic Differential Equations (SDE) offer a rich framework to model the probabilistic evolution of the state of a system. Numerical approximation methods are typically needed in evaluating relevant Quantities of Interest arising from such models. In this dissertation, we present novel effective methods for evaluating Quantities of Interest relevant to computational finance when the state of the system is described by an SDE.

Advances in Reactor physics, mathematics and computation. Volume 3

Energy Technology Data Exchange (ETDEWEB)

1987-01-01

These proceedings of the international topical meeting on advances in reactor physics, mathematics and computation, volume 3, are divided into sessions bearing on: - poster sessions on benchmark and codes: 35 conferences - review of status of assembly spectrum codes: 9 conferences - Numerical methods in fluid mechanics and thermal hydraulics: 16 conferences - stochastic transport and methods: 7 conferences.

The influence of numerical models on determining the drag coefficient

Directory of Open Access Journals (Sweden)

Dobeš Josef

2014-03-01

Full Text Available The paper deals with numerical modelling of body aerodynamic drag coefficient in the transition from laminar to turbulent flow regimes, where the selection of a suitable numerical model is problematic. On the basic problem of flow around a simple body – sphere selected computational models are tested. The values obtained by numerical simulations of drag coefficients of each model are compared with the graph of dependency of the drag coefficient vs. Reynolds number for a sphere. Next the dependency of Strouhal number vs. Reynolds number is evaluated, where the vortex shedding frequency values for given speed are obtained numerically and experimentally and then the values are compared for each numerical model and experiment. The aim is to specify trends for the selection of appropriate numerical model for flow around bodies problem in which the precise description of the flow field around the obstacle is used to define the acoustic noise source. Numerical modelling is performed by finite volume method using CFD code.

Domain of composition and finite volume schemes on non-matching grids; Decomposition de domaine et schemas volumes finis sur maillages non-conformes

Energy Technology Data Exchange (ETDEWEB)

Saas, L.

2004-05-01

This Thesis deals with sedimentary basin modeling whose goal is the prediction through geological times of the localizations and appraisal of hydrocarbons quantities present in the ground. Due to the natural and evolutionary decomposition of the sedimentary basin in blocks and stratigraphic layers, domain decomposition methods are requested to simulate flows of waters and of hydrocarbons in the ground. Conservations laws are used to model the flows in the ground and form coupled partial differential equations which must be discretized by finite volume method. In this report we carry out a study on finite volume methods on non-matching grids solved by domain decomposition methods. We describe a family of finite volume schemes on non-matching grids and we prove that the associated global discretized problem is well posed. Then we give an error estimate. We give two examples of finite volume schemes on non matching grids and the corresponding theoretical results (Constant scheme and Linear scheme). Then we present the resolution of the global discretized problem by a domain decomposition method using arbitrary interface conditions (for example Robin conditions). Finally we give numerical results which validate the theoretical results and study the use of finite volume methods on non-matching grids for basin modeling. (author)

Three-body unitarity in the finite volume

Energy Technology Data Exchange (ETDEWEB)

Mai, M. [The George Washington University, Washington, DC (United States); Doering, M. [The George Washington University, Washington, DC (United States); Thomas Jefferson National Accelerator Facility, Newport News, VA (United States)

2017-12-15

The physical interpretation of lattice QCD simulations, performed in a small volume, requires an extrapolation to the infinite volume. A method is proposed to perform such an extrapolation for three interacting particles at energies above threshold. For this, a recently formulated relativistic 3 → 3 amplitude based on the isobar formulation is adapted to the finite volume. The guiding principle is two- and three-body unitarity that imposes the imaginary parts of the amplitude in the infinite volume. In turn, these imaginary parts dictate the leading power-law finite-volume effects. It is demonstrated that finite-volume poles arising from the singular interaction, from the external two-body sub-amplitudes, and from the disconnected topology cancel exactly leaving only the genuine three-body eigenvalues. The corresponding quantization condition is derived for the case of three identical scalar-isoscalar particles and its numerical implementation is demonstrated. (orig.)

The pressure equation arising in reservoir simulation. Mathematical properties, numerical methods and upscaling

Energy Technology Data Exchange (ETDEWEB)

Nielsen, Bjoern Fredrik

1997-12-31

The main purpose of this thesis has been to analyse self-adjoint second order elliptic partial differential equations arising in reservoir simulation. It studies several mathematical and numerical problems for the pressure equation arising in models of fluid flow in porous media. The theoretical results obtained have been illustrated by a series of numerical experiments. The influence of large variations in the mobility tensor upon the solution of the pressure equation is analysed. The performance of numerical methods applied to such problems have been studied. A new upscaling technique for one-phase flow in heterogeneous reservoirs is developed. The stability of the solution of the pressure equation with respect to small perturbations of the mobility tensor is studied. The results are used to develop a new numerical method for a model of fully nonlinear water waves. 158 refs, 39 figs., 12 tabs.

The pressure equation arising in reservoir simulation. Mathematical properties, numerical methods and upscaling

Energy Technology Data Exchange (ETDEWEB)

Nielsen, Bjoern Fredrik

1998-12-31

The main purpose of this thesis has been to analyse self-adjoint second order elliptic partial differential equations arising in reservoir simulation. It studies several mathematical and numerical problems for the pressure equation arising in models of fluid flow in porous media. The theoretical results obtained have been illustrated by a series of numerical experiments. The influence of large variations in the mobility tensor upon the solution of the pressure equation is analysed. The performance of numerical methods applied to such problems have been studied. A new upscaling technique for one-phase flow in heterogeneous reservoirs is developed. The stability of the solution of the pressure equation with respect to small perturbations of the mobility tensor is studied. The results are used to develop a new numerical method for a model of fully nonlinear water waves. 158 refs, 39 figs., 12 tabs.

Application of Numerical Integration and Data Fusion in Unit Vector Method

Science.gov (United States)

Zhang, J.

2012-01-01

The Unit Vector Method (UVM) is a series of orbit determination methods which are designed by Purple Mountain Observatory (PMO) and have been applied extensively. It gets the conditional equations for different kinds of data by projecting the basic equation to different unit vectors, and it suits for weighted process for different kinds of data. The high-precision data can play a major role in orbit determination, and accuracy of orbit determination is improved obviously. The improved UVM (PUVM2) promoted the UVM from initial orbit determination to orbit improvement, and unified the initial orbit determination and orbit improvement dynamically. The precision and efficiency are improved further. In this thesis, further research work has been done based on the UVM: Firstly, for the improvement of methods and techniques for observation, the types and decision of the observational data are improved substantially, it is also asked to improve the decision of orbit determination. The analytical perturbation can not meet the requirement. So, the numerical integration for calculating the perturbation has been introduced into the UVM. The accuracy of dynamical model suits for the accuracy of the real data, and the condition equations of UVM are modified accordingly. The accuracy of orbit determination is improved further. Secondly, data fusion method has been introduced into the UVM. The convergence mechanism and the defect of weighted strategy have been made clear in original UVM. The problem has been solved in this method, the calculation of approximate state transition matrix is simplified and the weighted strategy has been improved for the data with different dimension and different precision. Results of orbit determination of simulation and real data show that the work of this thesis is effective: (1) After the numerical integration has been introduced into the UVM, the accuracy of orbit determination is improved obviously, and it suits for the high-accuracy data of

Finite-State Mean-Field Games, Crowd Motion Problems, and its Numerical Methods

KAUST Repository

Machado Velho, Roberto

2017-09-10

In this dissertation, we present two research projects, namely finite-state mean-field games and the Hughes model for the motion of crowds. In the first part, we describe finite-state mean-field games and some applications to socio-economic sciences. Examples include paradigm shifts in the scientific community and the consumer choice behavior in a free market. The corresponding finite-state mean-field game models are hyperbolic systems of partial differential equations, for which we propose and validate a new numerical method. Next, we consider the dual formulation to two-state mean-field games, and we discuss numerical methods for these problems. We then depict different computational experiments, exhibiting a variety of behaviors, including shock formation, lack of invertibility, and monotonicity loss. We conclude the first part of this dissertation with an investigation of the shock structure for two-state problems. In the second part, we consider a model for the movement of crowds proposed by R. Hughes in [56] and describe a numerical approach to solve it. This model comprises a Fokker-Planck equation coupled with an Eikonal equation with Dirichlet or Neumann data. We first establish a priori estimates for the solutions. Next, we consider radial solutions, and we identify a shock formation mechanism. Subsequently, we illustrate the existence of congestion, the breakdown of the model, and the trend to the equilibrium. We also propose a new numerical method for the solution of Fokker-Planck equations and then to systems of PDEs composed by a Fokker-Planck equation and a potential type equation. Finally, we illustrate the use of the numerical method both to the Hughes model and mean-field games. We also depict cases such as the evacuation of a room and the movement of persons around Kaaba (Saudi Arabia).

Experimental Results and Numerical Simulation of the Target RCS using Gaussian Beam Summation Method

Directory of Open Access Journals (Sweden)

Ghanmi Helmi

2018-05-01

Full Text Available This paper presents a numerical and experimental study of Radar Cross Section (RCS of radar targets using Gaussian Beam Summation (GBS method. The purpose GBS method has several advantages over ray method, mainly on the caustic problem. To evaluate the performance of the chosen method, we started the analysis of the RCS using Gaussian Beam Summation (GBS and Gaussian Beam Launching (GBL, the asymptotic models Physical Optic (PO, Geometrical Theory of Diffraction (GTD and the rigorous Method of Moment (MoM. Then, we showed the experimental validation of the numerical results using experimental measurements which have been executed in the anechoic chamber of Lab-STICC at ENSTA Bretagne. The numerical and experimental results of the RCS are studied and given as a function of various parameters: polarization type, target size, Gaussian beams number and Gaussian beams width.

NUMERICAL METHODS FOR SOLVING THE MULTI-TERM TIME-FRACTIONAL WAVE-DIFFUSION EQUATION.

Science.gov (United States)

Liu, F; Meerschaert, M M; McGough, R J; Zhuang, P; Liu, Q

2013-03-01

In this paper, the multi-term time-fractional wave-diffusion equations are considered. The multi-term time fractional derivatives are defined in the Caputo sense, whose orders belong to the intervals [0,1], [1,2), [0,2), [0,3), [2,3) and [2,4), respectively. Some computationally effective numerical methods are proposed for simulating the multi-term time-fractional wave-diffusion equations. The numerical results demonstrate the effectiveness of theoretical analysis. These methods and techniques can also be extended to other kinds of the multi-term fractional time-space models with fractional Laplacian.

Droplet spreading and capillary imbibition in a porous medium: A coupled IB-VOF method based numerical study

Science.gov (United States)

Das, Saurish; Patel, H. V.; Milacic, E.; Deen, N. G.; Kuipers, J. A. M.

2018-01-01

We investigate the dynamics of a liquid droplet in contact with a surface of a porous structure by means of the pore-scale level, fully resolved numerical simulations. The geometrical details of the solid porous matrix are resolved by a sharp interface immersed boundary method on a Cartesian computational grid, whereas the motion of the gas-liquid interface is tracked by a mass conservative volume of fluid method. The numerical simulations are performed considering a model porous structure that is approximated by a 3D cubical scaffold with cylindrical struts. The effect of the porosity and the equilibrium contact angle (between the gas-liquid interface and the solid struts) on the spreading behavior, liquid imbibition, and apparent contact angle (between the gas-liquid interface and the porous base) are studied. We also perform several simulations for droplet spreading on a flat surface as a reference case. Gas-liquid systems of the Laplace number, La = 45 and La = 144 × 103 are considered neglecting the effect of gravity. We report the time exponent (n) and pre-factor (C) of the power law describing the evolution of the spreading diameter (S = Ctn) for different equilibrium contact angles and porosity. Our simulations reveal that the apparent or macroscopic contact angle varies linearly with the equilibrium contact angle and increases with porosity. Not necessarily for all the wetting porous structures, a continuous capillary drainage occurs, and we find that the rate of the capillary drainage very much depends on the fluid inertia. At La = 144 × 103, numerically we capture the capillary wave induced pinch-off and daughter droplet ejection. We observe that on the porous structure the pinch-off is weak compared to that on a flat plate.

Radionuclide method for blood volume determination in kidneys

International Nuclear Information System (INIS)

Trindev, P.; Nikolov, D.; Shejretova, E.; Garcheva-Tsacheva, M.

1989-01-01

The method is applied in nephrology for diagnosing changes in blood circulation of the kidneys. The blood volume of each kidney is determined separately by perfusion angioscintigraphy (PAS) with improved accuracy. The method consists in intravenous injection of 300-450 MBq 99m Tc for in-vivo labelling of the erythrocytes. About 30 images are registered every 2 sec, and through zones of interest perfusion histograms of kidneys are derived. Ten minutes later kidneys images (one full-face and two profiles) are registered. Correction coefficients for kidneys depth are derived and the activities registered according to full-face images and amplitudes of perfusion histograms are corrected. The activity of 1 ml blood is determined from blood sample of the patient. The blood volume of each kidney is expressed as a ratio of the activity corrected for background and depth and the activity of 1 ml blood of the sample. 1 claim

Numerical simulation of long-term radiation effects for MOSFETs

International Nuclear Information System (INIS)

Wei Yuan; Xie Honggang; Gong Ding; Zhu Jinhui; Niu Shengli; Huang Liuxing

2013-01-01

A coupled algorithm is introduced to simulate the long-term radiation effects of MOSFETs, which combines particle transport with semiconductor governing equations. The former is dealt with Monte-Carlo method, and the latter is solved by finite-volume method. The trapped charge in SiO 2 and the free charge in Si are both described by the drift-diffusion model, and the deposited energy by incident particles can be coupled with the continuous equations of charge, acting as a source item. The discrete form of governing equations is obtained using the finite-volume method, and the numerical solutions of these equations are the long-term radiation response result of MOSFETs. The threshold voltage shift and off-state leakage current of an irradiated MOSFET are simulated with the coupled algorithm respectively, showing a good accordance with results by other calculations. (authors)

Simulation of Intra-Aneurysmal Blood Flow by Different Numerical Methods

Directory of Open Access Journals (Sweden)

Frank Weichert

2013-01-01

Full Text Available The occlusional performance of sole endoluminal stenting of intracranial aneurysms is controversially discussed in the literature. Simulation of blood flow has been studied to shed light on possible causal attributions. The outcome, however, largely depends on the numerical method and various free parameters. The present study is therefore conducted to find ways to define parameters and efficiently explore the huge parameter space with finite element methods (FEMs and lattice Boltzmann methods (LBMs. The goal is to identify both the impact of different parameters on the results of computational fluid dynamics (CFD and their advantages and disadvantages. CFD is applied to assess flow and aneurysmal vorticity in 2D and 3D models. To assess and compare initial simulation results, simplified 2D and 3D models based on key features of real geometries and medical expert knowledge were used. A result obtained from this analysis indicates that a combined use of the different numerical methods, LBM for fast exploration and FEM for a more in-depth look, may result in a better understanding of blood flow and may also lead to more accurate information about factors that influence conditions for stenting of intracranial aneurysms.

On a New Method for Computing the Numerical Solution of Systems of Nonlinear Equations

Directory of Open Access Journals (Sweden)

H. Montazeri

2012-01-01

Full Text Available We consider a system of nonlinear equations F(x=0. A new iterative method for solving this problem numerically is suggested. The analytical discussions of the method are provided to reveal its sixth order of convergence. A discussion on the efficiency index of the contribution with comparison to the other iterative methods is also given. Finally, numerical tests illustrate the theoretical aspects using the programming package Mathematica.

Novel Parallel Numerical Methods for Radiation and Neutron Transport

International Nuclear Information System (INIS)

Brown, P N

2001-01-01