#### Sample records for advanced numerical schemes

1. Using Steffe's Advanced Fraction Schemes

ERIC Educational Resources Information Center

McCloskey, Andrea V.; Norton, Anderson H.

2009-01-01

Recognizing schemes, which are different from strategies, can help teachers understand their students' thinking about fractions. Using Steffe's advanced fraction schemes, the authors describe a progression of development that upper elementary and middle school students might follow in understanding fractions. Each scheme can be viewed as a…

2. Development of a numerical scheme to predict geomagnetic storms after intense solar events and geomagnetic activity 27 days in advance. Final report, 6 Aug 86-16 Nov 90

SciTech Connect

Akasofu, S.I.; Lee, L.H.

1991-02-01

The modern geomagnetic storm prediction scheme should be based on a numerical simulation method, rather than on a statistical result. Furthermore, the scheme should be able to predict the geomagnetic storm indices, such as the Dst and AE indices, as a function of time. By recognizing that geomagnetic storms are powered by the solar wind-magnetosphere generator and that its power is given in terms of the solar wind speed, the interplanetary magnetic field (IMF) magnitude and polar angle, the authors have made a major advance in predicting both flare-induced storms and recurrent storms. Furthermore, it is demonstrated that the prediction scheme can be calibrated using the interplanetary scintillation (IPS) observation, when the solar disturbance advances about half-way to the earth. It is shown, however, that we are still far from a reliable prediction scheme. The prediction of the IMF polar angle requires future advance in understanding characteristics of magnetic clouds.

3. A numerical scheme for ionizing shock waves

SciTech Connect

Aslan, Necdet . E-mail: naslan@yeditepe.edu.tr; Mond, Michael

2005-12-10

A two-dimensional (2D) visual computer code to solve the steady state (SS) or transient shock problems including partially ionizing plasma is presented. Since the flows considered are hypersonic and the resulting temperatures are high, the plasma is partially ionized. Hence the plasma constituents are electrons, ions and neutral atoms. It is assumed that all the above species are in thermal equilibrium, namely, that they all have the same temperature. The ionization degree is calculated from Saha equation as a function of electron density and pressure by means of a nonlinear Newton type root finding algorithms. The code utilizes a wave model and numerical fluctuation distribution (FD) scheme that runs on structured or unstructured triangular meshes. This scheme is based on evaluating the mesh averaged fluctuations arising from a number of waves and distributing them to the nodes of these meshes in an upwind manner. The physical properties (directions, strengths, etc.) of these wave patterns are obtained by a new wave model: ION-A developed from the eigen-system of the flux Jacobian matrices. Since the equation of state (EOS) which is used to close up the conservation laws includes electronic effects, it is a nonlinear function and it must be inverted by iterations to determine the ionization degree as a function of density and temperature. For the time advancement, the scheme utilizes a multi-stage Runge-Kutta (RK) algorithm with time steps carefully evaluated from the maximum possible propagation speed in the solution domain. The code runs interactively with the user and allows to create different meshes to use different initial and boundary conditions and to see changes of desired physical quantities in the form of color and vector graphics. The details of the visual properties of the code has been published before (see [N. Aslan, A visual fluctuation splitting scheme for magneto-hydrodynamics with a new sonic fix and Euler limit, J. Comput. Phys. 197 (2004) 1

4. Multidimensional numerical scheme for resistive relativistic magnetohydrodynamics

Komissarov, Serguei S.

2007-12-01

The paper describes a new upwind conservative numerical scheme for special relativistic resistive magnetohydrodynamics with scalar resistivity. The magnetic field is kept approximately divergence free and the divergence of the electric field is kept consistent with the electric charge distribution via the method of Generalized Lagrange Multiplier. The hyperbolic fluxes are computed using the Harten-Lax-van Leer (HLL) prescription and the source terms are accounted via the time-splitting technique. The results of test simulations show that the scheme can handle equally well both resistive current sheets and shock waves, and thus can be a useful tool for studying phenomena of relativistic astrophysics that involve both colliding supersonic flows and magnetic reconnection.

5. Numerical Schemes for Rough Parabolic Equations

SciTech Connect

Deya, Aurelien

2012-04-15

This paper is devoted to the study of numerical approximation schemes for a class of parabolic equations on (0,1) perturbed by a non-linear rough signal. It is the continuation of Deya (Electron. J. Probab. 16:1489-1518, 2011) and Deya et al. (Probab. Theory Relat. Fields, to appear), where the existence and uniqueness of a solution has been established. The approach combines rough paths methods with standard considerations on discretizing stochastic PDEs. The results apply to a geometric 2-rough path, which covers the case of the multidimensional fractional Brownian motion with Hurst index H>1/3.

6. Numerical viscosity and the entropy condition for conservative difference schemes

NASA Technical Reports Server (NTRS)

1983-01-01

Consider a scalar, nonlinear conservative difference scheme satisfying the entropy condition. It is shown that difference schemes containing more numerical viscosity will necessarily converge to the unique, physically relevant weak solution of the approximated conservation equation. In particular, entropy satisfying convergence follows for E schemes - those containing more numerical viscosity than Godunov's scheme.

7. Simple Numerical Schemes for the Korteweg-deVries Equation

SciTech Connect

C. J. McKinstrie; M. V. Kozlov

2000-12-01

Two numerical schemes, which simulate the propagation of dispersive non-linear waves, are described. The first is a split-step Fourier scheme for the Korteweg-de Vries (KdV) equation. The second is a finite-difference scheme for the modified KdV equation. The stability and accuracy of both schemes are discussed. These simple schemes can be used to study a wide variety of physical processes that involve dispersive nonlinear waves.

8. Study on the numerical schemes for hypersonic flow simulation

Nagdewe, S. P.; Shevare, G. R.; Kim, Heuy-Dong

2009-10-01

Hypersonic flow is full of complex physical and chemical processes, hence its investigation needs careful analysis of existing schemes and choosing a suitable scheme or designing a brand new scheme. The present study deals with two numerical schemes Harten, Lax, and van Leer with Contact (HLLC) and advection upstream splitting method (AUSM) to effectively simulate hypersonic flow fields, and accurately predict shock waves with minimal diffusion. In present computations, hypersonic flows have been modeled as a system of hyperbolic equations with one additional equation for non-equilibrium energy and relaxing source terms. Real gas effects, which appear typically in hypersonic flows, have been simulated through energy relaxation method. HLLC and AUSM methods are modified to incorporate the conservation laws for non-equilibrium energy. Numerical implementation have shown that non-equilibrium energy convect with mass, and hence has no bearing on the basic numerical scheme. The numerical simulation carried out shows good comparison with experimental data available in literature. Both numerical schemes have shown identical results at equilibrium. Present study has demonstrated that real gas effects in hypersonic flows can be modeled through energy relaxation method along with either AUSM or HLLC numerical scheme.

9. High order parallel numerical schemes for solving incompressible flows

NASA Technical Reports Server (NTRS)

Lin, Avi; Milner, Edward J.; Liou, May-Fun; Belch, Richard A.

1992-01-01

The use of parallel computers for numerically solving flow fields has gained much importance in recent years. This paper introduces a new high order numerical scheme for computational fluid dynamics (CFD) specifically designed for parallel computational environments. A distributed MIMD system gives the flexibility of treating different elements of the governing equations with totally different numerical schemes in different regions of the flow field. The parallel decomposition of the governing operator to be solved is the primary parallel split. The primary parallel split was studied using a hypercube like architecture having clusters of shared memory processors at each node. The approach is demonstrated using examples of simple steady state incompressible flows. Future studies should investigate the secondary split because, depending on the numerical scheme that each of the processors applies and the nature of the flow in the specific subdomain, it may be possible for a processor to seek better, or higher order, schemes for its particular subcase.

10. Numerical simulation of shock wave diffraction by TVD schemes

NASA Technical Reports Server (NTRS)

Young, Victor Y. C.; Yee, H. C.

1987-01-01

An upwind total variation diminishing (TVD) scheme and a predictor-corrector symmetric TVD scheme were used to numerically simulate the blast wave diffraction on a stationary object. The objective is to help design an optimum configuration so that lateral motion is minimized and at the same time vortex shedding and flow separation are reduced during a blast wave encounter. Results are presented for a generic configuration for both a coarse grid and a fine grid to illustrate the global and local diffraction flow fields. Numerical experiments for the shock wave reflection on a wedge are also included to validate the current approach. Numerical study indicated that these TVD schemes are more stable and produced higher shock resolution than classical shock capturing methods such as the explicit MacCormack scheme.

11. A hybrid numerical scheme for the numerical solution of the Burgers' equation

Jiwari, Ram

2015-03-01

In this article, a hybrid numerical scheme based on Euler implicit method, quasilinearization and uniform Haar wavelets has been developed for the numerical solutions of Burgers' equation. Most of the numerical methods available in the literature fail to capture the physical behavior of the equations when viscosity ν → 0. In Jiwari (2012), the author presented the numerical results up to ν = 0.003 and the scheme failed for values smaller than ν = 0.003. The main aim in the development of the present scheme is to overcome the drawback of the scheme developed in Jiwari (2012). Lastly, three test problems are chosen to check the accuracy of the proposed scheme. The approximated results are compared with existing numerical and exact solutions found in literature. The use of uniform Haar wavelet is found to be accurate, simple, fast, flexible, convenient and at small computation costs.

12. Suppressing the numerical Cherenkov radiation in the Yee numerical scheme

2016-01-01

The next generation of laser facilities will routinely produce relativistic particle beams from the interaction of intense laser pulses with solids and/or gases. Their modeling with Particle-In-Cell (PIC) codes needs dispersion-free Maxwell solvers in order to properly describe the interaction of electromagnetic waves with relativistic particles. A particular attention is devoted to the suppression of the numerical Cherenkov instability, responsible for the noise generation. It occurs when the electromagnetic wave is artificially slowed down because of the finite mesh size, thus allowing for the high energy particles to propagate with super-luminous velocities. In the present paper, we show how a slight increase of the light velocity in the Maxwell's equations enables to suppress this instability while keeping a good overall precision of calculations.

13. Fundamental Limitations in Advanced LC Schemes

SciTech Connect

Mikhailichenko, A. A.

2010-11-04

Fundamental limitations in acceleration gradient, emittance, alignment and polarization in acceleration schemes are considered in application for novel schemes of acceleration, including laser-plasma and structure-based schemes. Problems for each method are underlined whenever it is possible. Main attention is paid to the scheme with a tilted laser bunch.

14. A factored implicit scheme for numerical weather prediction

NASA Technical Reports Server (NTRS)

Augenbaum, J. M.; Cohn, S. E.; Isaacson, E.; Dee, D. P.; Marchesin, D.

1985-01-01

An implicit method is proposed to factor the nonlinear partial differential equations governing fast and slow modes of dynamic motion in numerical weather prediction schemes. The method permits separate factorization of the slow and fast modes of the implicit operator. A simple two-dimensional version of the system of three-dimensional equations governing atmospheric dynamics over shallow water was analyzed to assess the accuracy of the proposed method. It is shown that the method has a small error which is comparable to other discretization errors in the overall scheme.

15. A multidimensional numerical scheme for two-fluid relativistic magnetohydrodynamics

Barkov, Maxim; Komissarov, Serguei S.; Korolev, Vitaly; Zankovich, Andrey

2014-02-01

This paper describes an explicit multidimensional numerical scheme for special relativistic two-fluid magnetohydrodynamics of electron-positron plasma and a suit of test problems. The scheme utilizes Cartesian grid and the third-order weighted essentially non-oscillatory interpolation. Time integration is carried out using the third-order total variation diminishing method of Runge-Kutta type, thus ensuring overall third-order accuracy on smooth solutions. The magnetic field is kept near divergence-free by means of the method of generalized Lagrange multiplier. The test simulations, which include linear and non-linear continuous plasma waves, shock waves, strong explosions and the tearing instability, show that the scheme is sufficiently robust and confirm its accuracy.

16. Numerical analysis of boosting scheme for scalable NMR quantum computation

SciTech Connect

SaiToh, Akira; Kitagawa, Masahiro

2005-02-01

Among initialization schemes for ensemble quantum computation beginning at thermal equilibrium, the scheme proposed by Schulman and Vazirani [in Proceedings of the 31st ACM Symposium on Theory of Computing (STOC'99) (ACM Press, New York, 1999), pp. 322-329] is known for the simple quantum circuit to redistribute the biases (polarizations) of qubits and small time complexity. However, our numerical simulation shows that the number of qubits initialized by the scheme is rather smaller than expected from the von Neumann entropy because of an increase in the sum of the binary entropies of individual qubits, which indicates a growth in the total classical correlation. This result--namely, that there is such a significant growth in the total binary entropy--disagrees with that of their analysis.

17. Unsteady boundary layers with an intelligent numerical scheme

Cebeci, T.

1986-02-01

A numerical method has been developed to represent unsteady boundary layers with large flow reversal. It makes use of the characteristic box scheme which examines the finite-difference grid in relation to the magnitude and direction of local velocity and reaches and implements a decision to ensure that the Courant, Friedricks and Lewey stability criterion is not violated. The method has been applied to the problem of an impulsively started circular cylinder and the results, though generally consistent with those of van Dommelen and Shen obtained with a Lagrangian method, show some differences. The time step is identified as very important and, with the present intelligent numerical scheme, the results were readily extended to times far beyond those previously achieved with Eulerian methods. Extrapolation of the results suggests that the much-discussed singularity for this unsteady flow is the same as that of the corresponding steady flow.

18. Advances in numerical and applied mathematics

NASA Technical Reports Server (NTRS)

South, J. C., Jr. (Editor); Hussaini, M. Y. (Editor)

1986-01-01

This collection of papers covers some recent developments in numerical analysis and computational fluid dynamics. Some of these studies are of a fundamental nature. They address basic issues such as intermediate boundary conditions for approximate factorization schemes, existence and uniqueness of steady states for time dependent problems, and pitfalls of implicit time stepping. The other studies deal with modern numerical methods such as total variation diminishing schemes, higher order variants of vortex and particle methods, spectral multidomain techniques, and front tracking techniques. There is also a paper on adaptive grids. The fluid dynamics papers treat the classical problems of imcompressible flows in helically coiled pipes, vortex breakdown, and transonic flows.

19. Numerical Methods in Quantum Mechanics: Analysis of Numerical Schemes on One-Dimensional Schrodinger Wave Problems

Jones, Marvin Quenten, Jr.

The motion and behavior of quantum processes can be described by the Schrodinger equation using the wave function, Psi(x,t). The use of the Schrodinger equation to study quantum phenomena is known as Quantum Mechanics, akin to classical mechanics being the tool to study classical physics. This research focuses on the emphasis of numerical techniques: Finite-Difference, Fast Fourier Transform (spectral method), finite difference schemes such as the Leapfrog method and the Crank-Nicolson scheme and second quantization to solve and analyze the Schrodinger equation for the infinite square well problem, the free particle with periodic boundary conditions, the barrier problem, tight-binding hamiltonians and a potential wall problem. We discuss these techniques and the problems created to test how these different techniques draw both physical and numerical conclusions in a tabular summary. We observed both numerical stability and quantum stability (conservation of energy, probability, momentum, etc.). We found in our results that the Crank-Nicolson scheme is an unconditionally stable scheme and conserves probability (unitary), and momentum, though dissipative with energy. The time-independent problems conserved energy, momentum and were unitary, which is of interest, but we found when time-dependence was introduced, quantum stability (i.e. conservation of mass, momentum, etc.) was not implied by numerical stability. Hence, we observed schemes that were numerically stable, but not quantum stable as well as schemes that were quantum stable, but not numerically stable for all of time, t. We also observed that second quantization removed the issues with stability as the problem was transformed into a discrete problem. Moreover, all quantum information is conserved in second quantization. This method, however, does not work universally for all problems.

20. Numerical solution of a semilinear elliptic equation via difference scheme

2016-08-01

We consider the Bitsadze-Samarskii type nonlocal boundary value problem { -d/2v (t ) d t2 +B v (t ) =h (t ,v (t ) ) ,0 scheme. The numerical results are computed by MATLAB.

1. COMPARISON OF NUMERICAL SCHEMES FOR SOLVING A SPHERICAL PARTICLE DIFFUSION EQUATION

EPA Science Inventory

A new robust iterative numerical scheme was developed for a nonlinear diffusive model that described sorption dynamics in spherical particle suspensions. he numerical scheme had been applied to finite difference and finite element models that showed rapid convergence and stabilit...

2. Energy balanced numerical schemes with very high order. The Augmented Roe Flux ADER scheme. Application to the shallow water equations

Navas-Montilla, A.; Murillo, J.

2015-06-01

In this work, an ADER type finite volume numerical scheme is proposed as an extension of a first order solver based on weak solutions of RPs with source terms. The type of source terms considered here are a special but relevant type of source terms: their spatial integral is discontinuous. The relevant difference with other previously defined ADER schemes is that it considers the presence of the source term in the solutions of the DRP. Unlike the original ADER schemes, the proposed numerical scheme computes the RPs of the high order terms of the DRP departing from time derivatives of the fluxes as initial conditions for these RPs. Weak solutions of the RPs defined for the DRP are computed using an augmented version of the Roe solver that includes an extra wave that accounts for the contribution of the source term. The discretization done over the source term leads to an energy balanced numerical scheme that allows to obtain the exact solution for steady cases with independence of the grid refinement. In unsteady problems, the numerical scheme ensures the convergence to the exact solution. The numerical scheme is constructed with an arbitrary order of accuracy, and has no theoretical barrier. Numerical results for the Burger's equation and the shallow water equations are presented in this work and indicate that the proposed numerical scheme is able to converge with the expected order of accuracy.

3. A numerical scheme for coastal morphodynamic modelling on unstructured grids

Guerin, Thomas; Bertin, Xavier; Dodet, Guillaume

2016-08-01

Over the last decade, modelling systems based on unstructured grids have been appearing increasingly attractive to investigate the dynamics of coastal zones. However, the resolution of the sediment continuity equation to simulate bed evolution is a complex problem which often leads to the development of numerical oscillations. To overcome this problem, addition of artificial diffusion or bathymetric filters are commonly employed methods, although these techniques can potentially over-smooth the bathymetry. This study aims to present a numerical scheme based on the Weighted Essentially Non-Oscillatory (WENO) formalism to solve the bed continuity equation on unstructured grids in a finite volume formulation. The new solution is compared against a classical method, which combines a basic node-centered finite volume method with artificial diffusion, for three idealized test cases. This comparison reveals that a higher accuracy is obtained with our new method while the addition of diffusion appears inappropriate mainly due to the arbitrary choice of the diffusion coefficient. Moreover, the increased computation time associated with the WENO-based method to solve the bed continuity equation is negligible when considering a fully-coupled simulation with tides and waves. Finally, the application of the new method to the pluri-monthly evolution of an idealized inlet subjected to tides and waves shows the development of realistic bed features (e.g. secondary flood channels, ebb-delta sandbars, or oblique sandbars at the adjacent beaches), that are smoothed or nonexistent when using additional diffusion.

4. Variationally consistent discretization schemes and numerical algorithms for contact problems

Wohlmuth, Barbara

We consider variationally consistent discretization schemes for mechanical contact problems. Most of the results can also be applied to other variational inequalities, such as those for phase transition problems in porous media, for plasticity or for option pricing applications from finance. The starting point is to weakly incorporate the constraint into the setting and to reformulate the inequality in the displacement in terms of a saddle-point problem. Here, the Lagrange multiplier represents the surface forces, and the constraints are restricted to the boundary of the simulation domain. Having a uniform inf-sup bound, one can then establish optimal low-order a priori convergence rates for the discretization error in the primal and dual variables. In addition to the abstract framework of linear saddle-point theory, complementarity terms have to be taken into account. The resulting inequality system is solved by rewriting it equivalently by means of the non-linear complementarity function as a system of equations. Although it is not differentiable in the classical sense, semi-smooth Newton methods, yielding super-linear convergence rates, can be applied and easily implemented in terms of a primal-dual active set strategy. Quite often the solution of contact problems has a low regularity, and the efficiency of the approach can be improved by using adaptive refinement techniques. Different standard types, such as residual- and equilibrated-based a posteriori error estimators, can be designed based on the interpretation of the dual variable as Neumann boundary condition. For the fully dynamic setting it is of interest to apply energy-preserving time-integration schemes. However, the differential algebraic character of the system can result in high oscillations if standard methods are applied. A possible remedy is to modify the fully discretized system by a local redistribution of the mass. Numerical results in two and three dimensions illustrate the wide range of

5. Comparative study of numerical schemes of TVD3, UNO3-ACM and optimized compact scheme

NASA Technical Reports Server (NTRS)

Lee, Duck-Joo; Hwang, Chang-Jeon; Ko, Duck-Kon; Kim, Jae-Wook

1995-01-01

Three different schemes are employed to solve the benchmark problem. The first one is a conventional TVD-MUSCL (Monotone Upwind Schemes for Conservation Laws) scheme. The second scheme is a UNO3-ACM (Uniformly Non-Oscillatory Artificial Compression Method) scheme. The third scheme is an optimized compact finite difference scheme modified by us: the 4th order Runge Kutta time stepping, the 4th order pentadiagonal compact spatial discretization with the maximum resolution characteristics. The problems of category 1 are solved by using the second (UNO3-ACM) and third (Optimized Compact) schemes. The problems of category 2 are solved by using the first (TVD3) and second (UNO3-ACM) schemes. The problem of category 5 is solved by using the first (TVD3) scheme. It can be concluded from the present calculations that the Optimized Compact scheme and the UN03-ACM show good resolutions for category 1 and category 2 respectively.

6. The hybrid Eulerian Lagrangian numerical scheme tested with Chemistry

Hansen, A. B.; Sørensen, B.; Tarning-Andersen, P.; Christensen, J. H.; Brandt, J.; Kaas, E.

2012-11-01

A newly developed advection scheme, the Hybrid Eulerian Lagrangian (HEL) scheme, has been tested, including a module for atmospheric chemistry, including 58 chemical species, and compared to two other traditional advection schemes; a classical pseudospectral Eulerian method the Accurate Space Derivative (ASD) scheme and the bi-cubic semi-Lagrangian (SL) scheme using classical rotation tests. The rotation tests have been designed to test and compare the advection schemes for different spatial and temporal resolutions in different chemical conditions (rural and urban) and for different shapes (cone and slotted cylinder) giving the advection schemes different challenges with respect to relatively slow or fast chemistry and smooth or sharp gradients, respectively. In every test, error measures have been calculated and used for ranking of the advection schemes with respect to performance, i.e. lowest overall errors for all chemical species. Furthermore, the HEL and SL schemes have been compared in a shallow water model, demonstrating the performance in a more realistic non-linear deformation flow. The results in this paper show that the new advection scheme, HEL, by far outperforms both the Eulerian and semi-Lagrangian schemes with very low error estimates compared to the two other schemes. Although no analytic solution can be obtained for the performance in the non-linear shallow water model flow, the tracer distribution appears realistic as compared to LMCSL when a mixing between local parcel concentrations is introduced in HEL.

7. Numerical study of chemically reacting flows using an LU scheme

NASA Technical Reports Server (NTRS)

Shuen, Jian Shun; Yoon, Seokkwan

1988-01-01

A new computational fluid dynamic code has been developed for the study of mixing and chemical reactions in the flow fields of ramjets and scramjets. The code employs an implicit finite volume, lower-upper symmetric successive overrelaxation scheme for solving the complete two-dimensional Navier-Stokes equations and species transport equations in a fully-coupled and very efficient manner. The combustion processes are modeled by an 8-species, 14-step finite rate chemistry model whereas turbulence is simulated by a Baldwin-Lomax algebraic model. The validity of the code is demonstrated by comparing the numerical calculations with both experimental data and previous calculations of a cold flow helium injection into a straight channel and premixed hydrogen-air reacting flows in a ramped duct. The code is then used to calculate the mixing and chemical reactions of a hydrogen jet transversely injected into a supersonic airstream. Results are presented describing the flow field, the recirculation regions in front and behind the injector, and the chemical reactions.

8. The hybrid Eulerian Lagrangian numerical scheme tested with Chemistry

Hansen, A. B.; Sørensen, B.; Tarning-Andersen, P.; Christensen, J. H.; Brandt, J.; Kaas, E.

2012-12-01

A newly developed transport scheme, the Hybrid Eulerian Lagrangian (HEL) scheme, has been tested using a module for atmospheric chemistry, including 58 chemical species, and compared to two other traditional advection schemes; a classical pseudospectral Eulerian method the Accurate Space Derivative (ASD) scheme and the bi-cubic semi-Lagrangian (SL) scheme using classical rotation tests. The rotation tests have been designed to test and compare the advection schemes for different spatial and temporal resolutions in different chemical conditions (rural and urban) and for different shapes (cone and slotted cylinder). This gives the advection schemes different challenges with respect to relatively slow or fast chemistry and smooth or sharp gradients. In every test, error measures have been calculated and used for ranking of the advection schemes with respect to performance, i.e. lowest overall errors for all chemical species. The results presented show that the new transport scheme, HEL, by far outperforms both the Eulerian and semi-Lagrangian schemes with very low error estimates compared to the two other schemes.

9. Advanced numerics for multi-dimensional fluid flow calculations

NASA Technical Reports Server (NTRS)

Vanka, S. P.

1984-01-01

In recent years, there has been a growing interest in the development and use of mathematical models for the simulation of fluid flow, heat transfer and combustion processes in engineering equipment. The equations representing the multi-dimensional transport of mass, momenta and species are numerically solved by finite-difference or finite-element techniques. However despite the multiude of differencing schemes and solution algorithms, and the advancement of computing power, the calculation of multi-dimensional flows, especially three-dimensional flows, remains a mammoth task. The following discussion is concerned with the author's recent work on the construction of accurate discretization schemes for the partial derivatives, and the efficient solution of the set of nonlinear algebraic equations resulting after discretization. The present work has been jointly supported by the Ramjet Engine Division of the Wright Patterson Air Force Base, Ohio, and the NASA Lewis Research Center.

10. A faster numerical scheme for a coupled system modeling soil erosion and sediment transport

Le, M.-H.; Cordier, S.; Lucas, C.; Cerdan, O.

2015-02-01

Overland flow and soil erosion play an essential role in water quality and soil degradation. Such processes, involving the interactions between water flow and the bed sediment, are classically described by a well-established system coupling the shallow water equations and the Hairsine-Rose model. Numerical approximation of this coupled system requires advanced methods to preserve some important physical and mathematical properties; in particular, the steady states and the positivity of both water depth and sediment concentration. Recently, finite volume schemes based on Roe's solver have been proposed by Heng et al. (2009) and Kim et al. (2013) for one and two-dimensional problems. In their approach, an additional and artificial restriction on the time step is required to guarantee the positivity of sediment concentration. This artificial condition can lead the computation to be costly when dealing with very shallow flow and wet/dry fronts. The main result of this paper is to propose a new and faster scheme for which only the CFL condition of the shallow water equations is sufficient to preserve the positivity of sediment concentration. In addition, the numerical procedure of the erosion part can be used with any well-balanced and positivity preserving scheme of the shallow water equations. The proposed method is tested on classical benchmarks and also on a realistic configuration.

11. Numerical Schemes for the Hamilton-Jacobi and Level Set Equations on Triangulated Domains

NASA Technical Reports Server (NTRS)

Barth, Timothy J.; Sethian, James A.

1997-01-01

Borrowing from techniques developed for conservation law equations, numerical schemes which discretize the Hamilton-Jacobi (H-J), level set, and Eikonal equations on triangulated domains are presented. The first scheme is a provably monotone discretization for certain forms of the H-J equations. Unfortunately, the basic scheme lacks proper Lipschitz continuity of the numerical Hamiltonian. By employing a virtual edge flipping technique, Lipschitz continuity of the numerical flux is restored on acute triangulations. Next, schemes are introduced and developed based on the weaker concept of positive coefficient approximations for homogeneous Hamiltonians. These schemes possess a discrete maximum principle on arbitrary triangulations and naturally exhibit proper Lipschitz continuity of the numerical Hamiltonian. Finally, a class of Petrov-Galerkin approximations are considered. These schemes are stabilized via a least-squares bilinear form. The Petrov-Galerkin schemes do not possess a discrete maximum principle but generalize to high order accuracy.

12. Verification and comparison of four numerical schemes for a 1D viscoelastic blood flow model.

PubMed

Wang, Xiaofei; Fullana, Jose-Maria; Lagrée, Pierre-Yves

2015-01-01

A reliable and fast numerical scheme is crucial for the 1D simulation of blood flow in compliant vessels. In this paper, a 1D blood flow model is incorporated with a Kelvin-Voigt viscoelastic arterial wall. This leads to a nonlinear hyperbolic-parabolic system, which is then solved with four numerical schemes, namely: MacCormack, Taylor-Galerkin, monotonic upwind scheme for conservation law and local discontinuous Galerkin. The numerical schemes are tested on a single vessel, a simple bifurcation and a network with 55 arteries. The numerical solutions are checked favorably against analytical, semi-analytical solutions or clinical observations. Among the numerical schemes, comparisons are made in four important aspects: accuracy, ability to capture shock-like phenomena, computational speed and implementation complexity. The suitable conditions for the application of each scheme are discussed. PMID:25145651

13. An Advanced Leakage Scheme for Neutrino Treatment in Astrophysical Simulations

Perego, A.; Cabezón, R. M.; Käppeli, R.

2016-04-01

We present an Advanced Spectral Leakage (ASL) scheme to model neutrinos in the context of core-collapse supernovae (CCSNe) and compact binary mergers. Based on previous gray leakage schemes, the ASL scheme computes the neutrino cooling rates by interpolating local production and diffusion rates (relevant in optically thin and thick regimes, respectively) separately for discretized values of the neutrino energy. Neutrino trapped components are also modeled, based on equilibrium and timescale arguments. The better accuracy achieved by the spectral treatment allows a more reliable computation of neutrino heating rates in optically thin conditions. The scheme has been calibrated and tested against Boltzmann transport in the context of Newtonian spherically symmetric models of CCSNe. ASL shows a very good qualitative and a partial quantitative agreement for key quantities from collapse to a few hundreds of milliseconds after core bounce. We have proved the adaptability and flexibility of our ASL scheme, coupling it to an axisymmetric Eulerian and to a three-dimensional smoothed particle hydrodynamics code to simulate core collapse. Therefore, the neutrino treatment presented here is ideal for large parameter-space explorations, parametric studies, high-resolution tests, code developments, and long-term modeling of asymmetric configurations, where more detailed neutrino treatments are not available or are currently computationally too expensive.

14. LES of turbulent heat transfer: proper convection numerical schemes for temperature transport

Châtelain, A.; Ducros, F.; Métais, O.

2004-03-01

Large eddy simulations of two basic configurations (decay of isotropic turbulence, and the academic plane channel flow) with heat transfer have been performed comparing several convection numerical schemes, in order to discuss their ability to evaluate temperature fluctuations properly. Results are compared with the available incompressible heat transfer direct numerical simulation data. It is shown that the use of regularizing schemes (such as high order upwind type schemes) for the temperature transport equation in combination with centered schemes for momentum transport equation gives better results than the use of centred schemes for both equations.

15. Analyzing numerics of bulk microphysics schemes in community models: warm rain processes

Sednev, I.; Menon, S.

2012-08-01

Implementation of bulk cloud microphysics (BLK) parameterizations in atmospheric models of different scales has gained momentum in the last two decades. Utilization of these parameterizations in cloud-resolving models when timesteps used for the host model integration are a few seconds or less is justified from the point of view of cloud physics. However, mechanistic extrapolation of the applicability of BLK schemes to the regional or global scales and the utilization of timesteps of hundreds up to thousands of seconds affect both physics and numerics. We focus on the mathematical aspects of BLK schemes, such as stability and positive-definiteness. We provide a strict mathematical definition for the problem of warm rain formation. We also derive a general analytical condition (SM-criterion) that remains valid regardless of parameterizations for warm rain processes in an explicit Eulerian time integration framework used to advanced finite-difference equations, which govern warm rain formation processes in microphysics packages in the Community Atmosphere Model and the Weather Research and Forecasting model. The SM-criterion allows for the existence of a unique positive-definite stable mass-conserving numerical solution, imposes an additional constraint on the timestep permitted due to the microphysics (like the Courant-Friedrichs-Lewy condition for the advection equation), and prohibits use of any additional assumptions not included in the strict mathematical definition of the problem under consideration. By analyzing the numerics of warm rain processes in source codes of BLK schemes implemented in community models we provide general guidelines regarding the appropriate choice of time steps in these models.

16. Improvements to the RELAP5-3D Nearly-Implicit Numerical Scheme

SciTech Connect

Richard A. Riemke; Walter L. Weaver; RIchard R. Schultz

2005-05-01

The RELAP5-3D computer program has been improved with regard to its nearly-implicit numerical scheme for twophase flow and single-phase flow. Changes were made to the nearly-implicit numerical scheme finite difference momentum equations as follows: (1) added the velocity flip-flop mass/energy error mitigation logic, (2) added the modified Henry-Fauske choking model, (3) used the new time void fraction in the horizontal stratification force terms and gravity head, and (4) used an implicit form of the artificial viscosity. The code modifications allow the nearly-implicit numerical scheme to be more implicit and lead to enhanced numerical stability.

17. Evaluating two numerical advection schemes in HYCOM for eddy-resolving modelling of the Agulhas Current

Backeberg, B. C.; Bertino, L.; Johannessen, J. A.

2009-06-01

A 4th order advection scheme is applied in a nested eddy-resolving Hybrid Coordinate Ocean Model (HYCOM) of the greater Agulhas Current system for the purpose of testing advanced numerics as a means for improving the model simulation for eventual operational implementation. Model validation techniques comparing sea surface height variations, sea level skewness and variogram analyses to satellite altimetry measurements quantify that generally the 4th order advection scheme improves the realism of the model simulation. The most striking improvement over the standard 2nd order momentum advection scheme, is that the southern Agulhas Current is simulated as a well-defined meandering current, rather than a train of successive eddies. A better vertical structure and stronger poleward transports in the Agulhas Current core contribute toward a better southwestward penetration of the current, and its temperature field, implying a stronger Indo-Atlantic inter-ocean exchange. It is found that the transport, and hence this exchange, is sensitive to the occurrences of mesoscale features originating upstream in the Mozambique Channel and southern East Madagascar Current, and that the improved HYCOM simulation is well suited for further studies of these inter-actions.

18. Evaluating two numerical advection schemes in HYCOM for eddy-resolving modelling of the Agulhas Current

Backeberg, B. C.; Bertino, L.; Johannessen, J. A.

2009-02-01

A 4th order advection scheme is applied in a nested eddy-resolving Hybrid Coordinate Ocean Model (HYCOM) of the greater Agulhas Current system for the purpose of testing advanced numerics as a means for improving the model simulation for eventual operational implementation. Model validation techniques comparing sea surface height variations, sea level skewness and variogram analyses to satellite altimetry measurements quantify that generally the 4th order advection scheme improves the realism of the model simulation. The most striking improvement over the standard 2nd order momentum advection scheme, is that the Southern Agulhas Current is simulated as a well-defined meandering current, rather than a train of successive eddies. A better vertical structure and stronger poleward transports in the Agulhas Current core contribute toward a better southwestward penetration of the current, and its temperature field, implying a stronger Indo-Atlantic inter-ocean exchange. It is found that the transport, and hence this exchange, is sensitive to the occurrences of mesoscale features originating upstream in the Mozambique Channel and Southern East Madagascar Current, and that the improved HYCOM simulation is well suited for further studies of these inter-actions.

19. A benchmark study of numerical schemes for one-dimensional arterial blood flow modelling.

PubMed

Boileau, Etienne; Nithiarasu, Perumal; Blanco, Pablo J; Müller, Lucas O; Fossan, Fredrik Eikeland; Hellevik, Leif Rune; Donders, Wouter P; Huberts, Wouter; Willemet, Marie; Alastruey, Jordi

2015-10-01

Haemodynamical simulations using one-dimensional (1D) computational models exhibit many of the features of the systemic circulation under normal and diseased conditions. Recent interest in verifying 1D numerical schemes has led to the development of alternative experimental setups and the use of three-dimensional numerical models to acquire data not easily measured in vivo. In most studies to date, only one particular 1D scheme is tested. In this paper, we present a systematic comparison of six commonly used numerical schemes for 1D blood flow modelling: discontinuous Galerkin, locally conservative Galerkin, Galerkin least-squares finite element method, finite volume method, finite difference MacCormack method and a simplified trapezium rule method. Comparisons are made in a series of six benchmark test cases with an increasing degree of complexity. The accuracy of the numerical schemes is assessed by comparison with theoretical results, three-dimensional numerical data in compatible domains with distensible walls or experimental data in a network of silicone tubes. Results show a good agreement among all numerical schemes and their ability to capture the main features of pressure, flow and area waveforms in large arteries. All the information used in this study, including the input data for all benchmark cases, experimental data where available and numerical solutions for each scheme, is made publicly available online, providing a comprehensive reference data set to support the development of 1D models and numerical schemes. PMID:26100764

20. Multi-dimensional high-order numerical schemes for Lagrangian hydrodynamics

SciTech Connect

Dai, William W; Woodward, Paul R

2009-01-01

An approximate solver for multi-dimensional Riemann problems at grid points of unstructured meshes, and a numerical scheme for multi-dimensional hydrodynamics have been developed in this paper. The solver is simple, and is developed only for the use in numerical schemes for hydrodynamics. The scheme is truely multi-dimensional, is second order accurate in both space and time, and satisfies conservation laws exactly for mass, momentum, and total energy. The scheme has been tested through numerical examples involving strong shocks. It has been shown that the scheme offers the principle advantages of high-order Codunov schemes; robust operation in the presence of very strong shocks and thin shock fronts.

1. Eulerian-Lagrangian numerical scheme for simulating advection, dispersion, and transient storage in streams and a comparison of numerical methods

USGS Publications Warehouse

Cox, T.J.; Runkel, R.L.

2008-01-01

Past applications of one-dimensional advection, dispersion, and transient storage zone models have almost exclusively relied on a central differencing, Eulerian numerical approximation to the nonconservative form of the fundamental equation. However, there are scenarios where this approach generates unacceptable error. A new numerical scheme for this type of modeling is presented here that is based on tracking Lagrangian control volumes across a fixed (Eulerian) grid. Numerical tests are used to provide a direct comparison of the new scheme versus nonconservative Eulerian numerical methods, in terms of both accuracy and mass conservation. Key characteristics of systems for which the Lagrangian scheme performs better than the Eulerian scheme include: nonuniform flow fields, steep gradient plume fronts, and pulse and steady point source loadings in advection-dominated systems. A new analytical derivation is presented that provides insight into the loss of mass conservation in the nonconservative Eulerian scheme. This derivation shows that loss of mass conservation in the vicinity of spatial flow changes is directly proportional to the lateral inflow rate and the change in stream concentration due to the inflow. While the nonconservative Eulerian scheme has clearly worked well for past published applications, it is important for users to be aware of the scheme's limitations. ?? 2008 ASCE.

2. Brush seal numerical simulation: Concepts and advances

NASA Technical Reports Server (NTRS)

Braun, M. J.; Kudriavtsev, V. V.

1994-01-01

The development of the brush seal is considered to be most promising among the advanced type seals that are presently in use in the high speed turbomachinery. The brush is usually mounted on the stationary portions of the engine and has direct contact with the rotating element, in the process of limiting the 'unwanted' leakage flows between stages, or various engine cavities. This type of sealing technology is providing high (in comparison with conventional seals) pressure drops due mainly to the high packing density (around 100 bristles/sq mm), and brush compliance with the rotor motions. In the design of modern aerospace turbomachinery leakage flows between the stages must be minimal, thus contributing to the higher efficiency of the engine. Use of the brush seal instead of the labyrinth seal reduces the leakage flow by one order of magnitude. Brush seals also have been found to enhance dynamic performance, cost less, and are lighter than labyrinth seals. Even though industrial brush seals have been successfully developed through extensive experimentation, there is no comprehensive numerical methodology for the design or prediction of their performance. The existing analytical/numerical approaches are based on bulk flow models and do not allow the investigation of the effects of brush morphology (bristle arrangement), or brushes arrangement (number of brushes, spacing between them), on the pressure drops and flow leakage. An increase in the brush seal efficiency is clearly a complex problem that is closely related to the brush geometry and arrangement, and can be solved most likely only by means of a numerically distributed model.

3. Brush seal numerical simulation: Concepts and advances

Braun, M. J.; Kudriavtsev, V. V.

1994-07-01

The development of the brush seal is considered to be most promising among the advanced type seals that are presently in use in the high speed turbomachinery. The brush is usually mounted on the stationary portions of the engine and has direct contact with the rotating element, in the process of limiting the 'unwanted' leakage flows between stages, or various engine cavities. This type of sealing technology is providing high (in comparison with conventional seals) pressure drops due mainly to the high packing density (around 100 bristles/sq mm), and brush compliance with the rotor motions. In the design of modern aerospace turbomachinery leakage flows between the stages must be minimal, thus contributing to the higher efficiency of the engine. Use of the brush seal instead of the labyrinth seal reduces the leakage flow by one order of magnitude. Brush seals also have been found to enhance dynamic performance, cost less, and are lighter than labyrinth seals. Even though industrial brush seals have been successfully developed through extensive experimentation, there is no comprehensive numerical methodology for the design or prediction of their performance. The existing analytical/numerical approaches are based on bulk flow models and do not allow the investigation of the effects of brush morphology (bristle arrangement), or brushes arrangement (number of brushes, spacing between them), on the pressure drops and flow leakage. An increase in the brush seal efficiency is clearly a complex problem that is closely related to the brush geometry and arrangement, and can be solved most likely only by means of a numerically distributed model.

4. Positivity-preserving numerical schemes for multidimensional advection

NASA Technical Reports Server (NTRS)

Leonard, B. P.; Macvean, M. K.; Lock, A. P.

1993-01-01

This report describes the construction of an explicit, single time-step, conservative, finite-volume method for multidimensional advective flow, based on a uniformly third-order polynomial interpolation algorithm (UTOPIA). Particular attention is paid to the problem of flow-to-grid angle-dependent, anisotropic distortion typical of one-dimensional schemes used component-wise. The third-order multidimensional scheme automatically includes certain cross-difference terms that guarantee good isotropy (and stability). However, above first-order, polynomial-based advection schemes do not preserve positivity (the multidimensional analogue of monotonicity). For this reason, a multidimensional generalization of the first author's universal flux-limiter is sought. This is a very challenging problem. A simple flux-limiter can be found; but this introduces strong anisotropic distortion. A more sophisticated technique, limiting part of the flux and then restoring the isotropy-maintaining cross-terms afterwards, gives more satisfactory results. Test cases are confined to two dimensions; three-dimensional extensions are briefly discussed.

Energy Science and Technology Software Center (ESTSC)

2006-05-12

The model was developed as a part of the study, "Advanced Fuel Cycles and Waste Management", which was performed during 20032005 by an ad-hoc expert group under the Nuclear Development Committee in the OECD/NEA. The model was designed for an efficient conduct of nuclear fuel cycle scheme cost analyses. It is simple, transparent and offers users the capability to track down the cost analysis results. All the fuel cycle schemes considered in the model aremore » represented in a graphic format and all values related to a fuel cycle step are shown in the graphic interface, i.e., there are no hidden values embedded in the calculations. All data on the fuel cycle schemes considered in the study including mass flows, waste generation, cost data, and other data such as activities, decay heat and neutron sources of spent fuel and highlevel waste along time are included in the model and can be displayed. The user can modify easily the values of mass flows and/or cost parameters and see the corresponding changes in the results. The model calculates: frontend fuel cycle mass flows such as requirements of enrichment and conversion services and natural uranium; mass of waste based on the waste generation parameters and the mass flow; and all costs. It performs Monte Carlo simulations with changing the values of all unit costs within their respective ranges (from lower to upper bounds).« less

Energy Science and Technology Software Center (ESTSC)

2008-03-17

Version 00 SMAFS was developed as a part of the study, "Advanced Fuel Cycles and Waste Management", which was performed during 2003-2005 by an ad-hoc expert group under the Nuclear Development Committee in the OECD/NEA. The model was designed for an efficient conduct of nuclear fuel cycle scheme cost analyses. It is simple, transparent and offers users the capability to track down cost analysis results. All the fuel cycle schemes considered in the model aremore » represented in a graphic format and all values related to a fuel cycle step are shown in the graphic interface, i.e., there are no hidden values embedded in the calculations. All data on the fuel cycle schemes considered in the study including mass flows, waste generation, cost data, and other data such as activities, decay heat and neutron sources of spent fuel and high-level waste along time are included in the model and can be displayed. The user can easily modify values of mass flows and/or cost parameters and see corresponding changes in the results. The model calculates: front-end fuel cycle mass flows such as requirements of enrichment and conversion services and natural uranium; mass of waste based on the waste generation parameters and the mass flow; and all costs.« less

7. Analysis and design of numerical schemes for gas dynamics. 2: Artificial diffusion and discrete shock structure

NASA Technical Reports Server (NTRS)

Jameson, Antony

1994-01-01

The effect of artificial diffusion on discrete shock structures is examined for a family of schemes which includes scalar diffusion, convective upwind and split pressure (CUSP) schemes, and upwind schemes with characteristics splitting. The analysis leads to conditions on the diffusive flux such that stationary discrete shocks can contain a single interior point. The simplest formulation which meets these conditions is a CUSP scheme in which the coefficients of the pressure differences is fully determined by the coefficient of convective diffusion. It is also shown how both the characteristic and CUSP schemes can be modified to preserve constant stagnation enthalpy in steady flow, leading to four variants, the E and H-characteristic schemes, and the E and H-CUSP schemes. Numerical results are presented which confirm the properties of these schemes.

8. On some numerical scheme of solving diffraction problem on open and closed screens

SciTech Connect

Ryzhakov, Gleb V.

2015-03-10

In the paper, the problem of diffraction on thin ideally conductive screens is reduced to vector hypersingular integral equation with integral treated in the sense of finite Hadamard value. An numerical scheme to solve the equation is introduced. The scheme is based on piecewise approximation of unknown function. The advantage of the scheme is that integral of singular part is reduced to contour integral which can be analytically calculated so numerical calculation are significantly accelerated. Several examples of resulting numerical experiments are given in comparison with known theoretical and experimental data.

9. The numerical viscosity of entropy stable schemes for systems of conservation laws. I

NASA Technical Reports Server (NTRS)

1987-01-01

Discrete approximations to hyperbolic systems of conservation laws are studied. The amount of numerical viscosity present in such schemes is quantified and related to their entropy stability by means of comparison.To this end, conservative schemes which are also entropy-conservative are constructed. These entropy-conservative schemes enjoy second-order accuracy; moreover, they can be interpreted as piecewise-linear finite-element methods, and hence can be formulated on various mesh configurations. It is then shown that conservative schemes are entropy stable, if and (for three-point schemes) only they contain more viscosity than that present in the above-mentioned entropy-conservative ones.

10. SEAWAT 2000: modelling unstable flow and sensitivity to discretization levels and numerical schemes

Al-Maktoumi, A.; Lockington, D. A.; Volker, R. E.

2007-09-01

A systematic analysis shows how results from the finite difference code SEAWAT are sensitive to choice of grid dimension, time step, and numerical scheme for unstable flow problems. Guidelines to assist in selecting appropriate combinations of these factors are suggested. While the SEAWAT code has been tested for a wide range of problems, the sensitivity of results to spatial and temporal discretization levels and numerical schemes has not been studied in detail for unstable flow problems. Here, the Elder-Voss-Souza benchmark problem has been used to systematically explore the sensitivity of SEAWAT output to spatio-temporal resolution and numerical solver choice. A grid size of 0.38 and 0.60% of the total domain length and depth respectively is found to be fine enough to deliver results with acceptable accuracy for most of the numerical schemes when Courant number (Cr) is 0.1. All numerical solvers produced similar results for extremely fine meshes; however, some schemes converged faster than others. For instance, the 3rd-order total variation-diminishing method (TVD3) scheme converged at a much coarser mesh than the standard finite difference methods (SFDM) upstream weighting (UW) scheme. The sensitivity of the results to Cr number depends on the numerical scheme as expected.

11. Fourth-order compact schemes for the numerical simulation of coupled Burgers' equation

Bhatt, H. P.; Khaliq, A. Q. M.

2016-03-01

This paper introduces two new modified fourth-order exponential time differencing Runge-Kutta (ETDRK) schemes in combination with a global fourth-order compact finite difference scheme (in space) for direct integration of nonlinear coupled viscous Burgers' equations in their original form without using any transformations or linearization techniques. One scheme is a modification of the Cox and Matthews ETDRK4 scheme based on (1 , 3) -Padé approximation and other is a modification of Krogstad's ETDRK4-B scheme based on (2 , 2) -Padé approximation. Efficient versions of the proposed schemes are obtained by using a partial fraction splitting technique of rational functions. The stability properties of the proposed schemes are studied by plotting the stability regions, which provide an explanation of their behavior for dispersive and dissipative problems. The order of convergence of the schemes is examined empirically and found that the modification of ETDRK4 converges with the expected rate even if the initial data are nonsmooth. On the other hand, modification of ETDRK4-B suffers with order reduction if the initial data are nonsmooth. Several numerical experiments are carried out in order to demonstrate the performance and adaptability of the proposed schemes. The numerical results indicate that the proposed schemes provide better accuracy than other schemes available in the literature. Moreover, the results show that the modification of ETDRK4 is reliable and yields more accurate results than modification of ETDRK4-B, while solving problems with nonsmooth data or with high Reynolds number.

12. Adaptive Numerical Dissipative Control in High Order Schemes for Multi-D Non-Ideal MHD

NASA Technical Reports Server (NTRS)

Yee, H. C.; Sjoegreen, B.

2004-01-01

The goal is to extend our adaptive numerical dissipation control in high order filter schemes and our new divergence-free methods for ideal MHD to non-ideal MHD that include viscosity and resistivity. The key idea consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and leave the rest of the region free of numerical dissipation contamination. These scheme-independent detectors are capable of distinguishing shocks/shears, flame sheets, turbulent fluctuations and spurious high-frequency oscillations. The detection algorithm is based on an artificial compression method (ACM) (for shocks/shears), and redundant multi-resolution wavelets (WAV) (for the above types of flow feature). These filter approaches also provide a natural and efficient way for the minimization of Div(B) numerical error. The filter scheme consists of spatially sixth order or higher non-dissipative spatial difference operators as the base scheme for the inviscid flux derivatives. If necessary, a small amount of high order linear dissipation is used to remove spurious high frequency oscillations. For example, an eighth-order centered linear dissipation (AD8) might be included in conjunction with a spatially sixth-order base scheme. The inviscid difference operator is applied twice for the viscous flux derivatives. After the completion of a full time step of the base scheme step, the solution is adaptively filtered by the product of a 'flow detector' and the 'nonlinear dissipative portion' of a high-resolution shock-capturing scheme. In addition, the scheme independent wavelet flow detector can be used in conjunction with spatially compact, spectral or spectral element type of base schemes. The ACM and wavelet filter schemes using the dissipative portion of a second-order shock-capturing scheme with sixth-order spatial central base scheme for both the inviscid and viscous MHD flux

13. Finite-difference scheme for the numerical solution of the Schroedinger equation

NASA Technical Reports Server (NTRS)

1992-01-01

A finite-difference scheme for numerical integration of the Schroedinger equation is constructed. Asymptotically (r goes to infinity), the method gives the exact solution correct to terms of order r exp -2.

14. A Continuing Search for a Near-Perfect Numerical Flux Scheme. Part 1; [AUSM+

NASA Technical Reports Server (NTRS)

Liou, Meng-Sing

1994-01-01

While enjoying demonstrated improvement in accuracy, efficiency, and robustness over existing schemes, the Advection Upstream Splitting Scheme (AUSM) was found to have some deficiencies in extreme cases. This recent progress towards improving the AUSM while retaining its advantageous features is described. The new scheme, termed AUSM+, features: unification of velocity and Mach number splitting; exact capture of a single stationary shock; and improvement in accuracy. A general construction of the AUSM+ scheme is layed out and then focus is on the analysis of the a scheme and its mathematical properties, heretofore unreported. Monotonicity and positivity are proved, and a CFL-like condition is given for first and second order schemes and for generalized curvilinear co-ordinates. Finally, results of numerical tests on many problems are given to confirm the capability and improvements on a variety of problems including those failed by prominent schemes.

15. The numerical viscosity of entropy stable schemes for systems of conservation laws

NASA Technical Reports Server (NTRS)

1985-01-01

Discrete approximations to hyperbolic systems of conservation laws are studied. The amount of numerical viscosity present in such schemes, is quantified and related to their entropy stability by means of comparison. To this end, conservative schemes which are also entropy conservative are constructed. These entropy conservative schemes enjoy second-order accuracy; moreover, they admit a particular interpretation within the finite-element frameworks, and hence can be formulated on various mesh configurations. It is then shown that conservative schemes are entropy stable if and only if they contain more viscosity than the mentioned above entropy conservative ones.

16. Determination of Solution Accuracy of Numerical Schemes as Part of Code and Calculation Verification

SciTech Connect

Blottner, F.G.; Lopez, A.R.

1998-10-01

This investigation is concerned with the accuracy of numerical schemes for solving partial differential equations used in science and engineering simulation codes. Richardson extrapolation methods for steady and unsteady problems with structured meshes are presented as part of the verification procedure to determine code and calculation accuracy. The local truncation error de- termination of a numerical difference scheme is shown to be a significant component of the veri- fication procedure as it determines the consistency of the numerical scheme, the order of the numerical scheme, and the restrictions on the mesh variation with a non-uniform mesh. Genera- tion of a series of co-located, refined meshes with the appropriate variation of mesh cell size is in- vestigated and is another important component of the verification procedure. The importance of mesh refinement studies is shown to be more significant than just a procedure to determine solu- tion accuracy. It is suggested that mesh refinement techniques can be developed to determine con- sistency of numerical schemes and to determine if governing equations are well posed. The present investigation provides further insight into the conditions and procedures required to effec- tively use Richardson extrapolation with mesh refinement studies to achieve confidence that sim- ulation codes are producing accurate numerical solutions.

17. Numerical Compression Schemes for Proteomics Mass Spectrometry Data*

PubMed Central

Teleman, Johan; Dowsey, Andrew W.; Gonzalez-Galarza, Faviel F.; Perkins, Simon; Pratt, Brian; Röst, Hannes L.; Malmström, Lars; Malmström, Johan; Jones, Andrew R.; Deutsch, Eric W.; Levander, Fredrik

2014-01-01

The open XML format mzML, used for representation of MS data, is pivotal for the development of platform-independent MS analysis software. Although conversion from vendor formats to mzML must take place on a platform on which the vendor libraries are available (i.e. Windows), once mzML files have been generated, they can be used on any platform. However, the mzML format has turned out to be less efficient than vendor formats. In many cases, the naïve mzML representation is fourfold or even up to 18-fold larger compared with the original vendor file. In disk I/O limited setups, a larger data file also leads to longer processing times, which is a problem given the data production rates of modern mass spectrometers. In an attempt to reduce this problem, we here present a family of numerical compression algorithms called MS-Numpress, intended for efficient compression of MS data. To facilitate ease of adoption, the algorithms target the binary data in the mzML standard, and support in main proteomics tools is already available. Using a test set of 10 representative MS data files we demonstrate typical file size decreases of 90% when combined with traditional compression, as well as read time decreases of up to 50%. It is envisaged that these improvements will be beneficial for data handling within the MS community. PMID:24677029

18. Numerical study of read scheme in one-selector one-resistor crossbar array

Kim, Sungho; Kim, Hee-Dong; Choi, Sung-Jin

2015-12-01

A comprehensive numerical circuit analysis of read schemes of a one selector-one resistance change memory (1S1R) crossbar array is carried out. Three schemes-the ground, V/2, and V/3 schemes-are compared with each other in terms of sensing margin and power consumption. Without the aid of a complex analytical approach or SPICE-based simulation, a simple numerical iteration method is developed to simulate entire current flows and node voltages within a crossbar array. Understanding such phenomena is essential in successfully evaluating the electrical specifications of selectors for suppressing intrinsic drawbacks of crossbar arrays, such as sneaky current paths and series line resistance problems. This method provides a quantitative tool for the accurate analysis of crossbar arrays and provides guidelines for developing an optimal read scheme, array configuration, and selector device specifications.

19. Advanced unrepeatered systems using novel Raman amplification schemes

Chang, Do-il; Pelouch, Wayne; Burtsev, Sergey; Perrier, Philippe; Fevrier, Herve

2015-01-01

Unrepeatered transmission systems provide a cost-effective solution to transmit high capacity channels in submarine networks to communicate between coastal population centers or in terrestrial networks to connect remote areas where service access is difficult. The main goal of unrepeatered systems has traditionally been to achieve the longest reach, however, increasing traffic demands now require unrepeatered systems to support both longer reach and higher transport capacity. As a result, transmission rate of unrepeatered systems has quickly moved from 10 Gb/s to 40 Gb/s or 100 Gb/s. This paper reviews the key basic technologies, with a specific focus on Raman amplification, required for long-reach, high-capacity unrepeatered optical transmission systems. We will discuss novel Raman amplification schemes, enhanced remote optically pumped amplifiers (ROPA), ultra-low loss / large effective area fibers, and coherent transmission with advanced modulation format and high FEC coding gain. We will also report recent experimental demonstrations that show how these technologies have been combined to achieve industry's leading capacity and reach transmission.

20. A practical numerical scheme for the ternary Cahn-Hilliard system with a logarithmic free energy

Jeong, Darae; Kim, Junseok

2016-01-01

We consider a practically stable finite difference method for the ternary Cahn-Hilliard system with a logarithmic free energy modeling the phase separation of a three-component mixture. The numerical scheme is based on a linear unconditionally gradient stable scheme by Eyre and is solved by an efficient and accurate multigrid method. The logarithmic function has a singularity at zero. To remove the singularity, we regularize the function near zero by using a quadratic polynomial approximation. We perform a convergence test, a linear stability analysis, and a robustness test of the ternary Cahn-Hilliard equation. We observe that our numerical solutions are convergent, consistent with the exact solutions of linear stability analysis, and stable with practically large enough time steps. Using the proposed numerical scheme, we also study the temporal evolution of morphology patterns during phase separation in one-, two-, and three-dimensional spaces.

1. Ancient numerical daemons of conceptual hydrological modeling: 1. Fidelity and efficiency of time stepping schemes

Clark, Martyn P.; Kavetski, Dmitri

2010-10-01

A major neglected weakness of many current hydrological models is the numerical method used to solve the governing model equations. This paper thoroughly evaluates several classes of time stepping schemes in terms of numerical reliability and computational efficiency in the context of conceptual hydrological modeling. Numerical experiments are carried out using 8 distinct time stepping algorithms and 6 different conceptual rainfall-runoff models, applied in a densely gauged experimental catchment, as well as in 12 basins with diverse physical and hydroclimatic characteristics. Results show that, over vast regions of the parameter space, the numerical errors of fixed-step explicit schemes commonly used in hydrology routinely dwarf the structural errors of the model conceptualization. This substantially degrades model predictions, but also, disturbingly, generates fortuitously adequate performance for parameter sets where numerical errors compensate for model structural errors. Simply running fixed-step explicit schemes with shorter time steps provides a poor balance between accuracy and efficiency: in some cases daily-step adaptive explicit schemes with moderate error tolerances achieved comparable or higher accuracy than 15 min fixed-step explicit approximations but were nearly 10 times more efficient. From the range of simple time stepping schemes investigated in this work, the fixed-step implicit Euler method and the adaptive explicit Heun method emerge as good practical choices for the majority of simulation scenarios. In combination with the companion paper, where impacts on model analysis, interpretation, and prediction are assessed, this two-part study vividly highlights the impact of numerical errors on critical performance aspects of conceptual hydrological models and provides practical guidelines for robust numerical implementation.

2. Numerical Forming Simulations and Optimisation in Advanced Materials

Huétink, J.; van den Boogaard, A. H.; Geijselears, H. J. M.; Meinders, T.

2007-05-01

With the introduction of new materials as high strength steels, metastable steels and fibre reinforced composites, the need for advanced physically valid constitutive models arises. In finite deformation problems constitutive relations are commonly formulated in terms the Cauchy stress as a function of the elastic Finger tensor and an objective rate of the Cauchy stress as a function of the rate of deformation tensor. For isotropic materials models this is rather straightforward, but for anisotropic material models, including elastic anisotropy as well as plastic anisotropy, this may lead to confusing formulations. It will be shown that it is more convenient to define the constitutive relations in terms of invariant tensors referred to the deformed metric. Experimental results are presented that show new combinations of strain rate and strain path sensitivity. An adaptive through- thickness integration scheme for plate elements is developed, which improves the accuracy of spring back prediction at minimal costs. A procedure is described to automatically compensate the CAD tool shape numerically to obtain the desired product shape. Forming processes need to be optimized for cost saving and product improvement. Until recently, a trial-and-error process in the factory primarily did this optimization. An optimisation strategy is proposed that assists an engineer to model an optimization problem that suits his needs, including an efficient algorithm for solving the problem.

3. A numerical study of ENO and TVD schemes for shock capturing

NASA Technical Reports Server (NTRS)

Chang, Shih-Hung; Liou, Meng-Sing

1988-01-01

The numerical performance of a second-order upwind-based total variation diminishing (TVD) scheme and that of a uniform second-order essentially non-oscillatory (ENO) scheme for shock capturing are compared. The TVD scheme used is a modified version of Liou, using the flux-difference splitting (FDS) of Roe and his superbee function as the limiter. The construction of the basic ENO scheme is based on Harten, Engquist, Osher, and Chakravarthy, and the 2-D extensions are obtained by using a Strang-type of fractional-step time-splitting method. Numerical results presented include both steady and unsteady, 1-D and 2-D calculations. All the chosen test problems have exact solutions so that numerical performance can be measured by comparing the computer results to them. For 1-D calculations, the standard shock-tube problems of Sod and Lax are chosen. A very strong shock-tube problem, with the initial density ratio of 400 to 1 and pressure ratio of 500 to 1, is also used to study the behavior of the two schemes. For 2-D calculations, the shock wave reflection problems are adopted for testing. The cases presented in this report include flows with Mach numbers of 2.9, 5.0, and 10.0.

4. Quantitative evaluation of numerical integration schemes for Lagrangian particle dispersion models

Ramli, Huda Mohd.; Esler, J. Gavin

2016-07-01

A rigorous methodology for the evaluation of integration schemes for Lagrangian particle dispersion models (LPDMs) is presented. A series of one-dimensional test problems are introduced, for which the Fokker-Planck equation is solved numerically using a finite-difference discretisation in physical space and a Hermite function expansion in velocity space. Numerical convergence errors in the Fokker-Planck equation solutions are shown to be much less than the statistical error associated with a practical-sized ensemble (N = 106) of LPDM solutions; hence, the former can be used to validate the latter. The test problems are then used to evaluate commonly used LPDM integration schemes. The results allow for optimal time-step selection for each scheme, given a required level of accuracy. The following recommendations are made for use in operational models. First, if computational constraints require the use of moderate to long time steps, it is more accurate to solve the random displacement model approximation to the LPDM rather than use existing schemes designed for long time steps. Second, useful gains in numerical accuracy can be obtained, at moderate additional computational cost, by using the relatively simple "small-noise" scheme of Honeycutt.

5. A block interface flux reconstruction method for numerical simulation with high order finite difference scheme

Gao, Junhui

2013-05-01

Overlap grid is usually used in numerical simulation of flow with complex geometry by high order finite difference scheme. It is difficult to generate overlap grid and the connectivity information between adjacent blocks, especially when interpolation is required for non-coincident overlap grids. In this study, an interface flux reconstruction (IFR) method is proposed for numerical simulation using high order finite difference scheme with multi-block structured grids. In this method the neighboring blocks share a common face, and the fluxes on each block are matched to set the boundary conditions for each interior block. Therefore this method has the promise of allowing discontinuous grids on either side of an interior block interface. The proposed method is proven to be stable for 7-point central DRP scheme coupled with 4-point and 5-point boundary closure schemes, as well as the 4th order compact scheme coupled with 3rd order boundary closure scheme. Four problems are numerically solved with the developed code to validate the interface flux reconstruction method in this study. The IFR method coupled with the 4th order DRP scheme or compact scheme is validated to be 4th order accuracy with one and two dimensional waves propagation problems. Two dimensional pulse propagation in mean flow is computed with wavy mesh to demonstrate the ability of the proposed method for non-uniform grid. To demonstrate the ability of the proposed method for complex geometry, sound scattering by two cylinders is simulated and the numerical results are compared with the analytical data. It is shown that the numerical results agree well with the analytical data. Finally the IFR method is applied to simulate viscous flow pass a cylinder at Reynolds number 150 to show its capability for viscous problem. The computed pressure coefficient on the cylinder surface, the frequency of vortex shedding, the lift and drag coefficients are presented. The numerical results are compared with the data

6. RELAP5 two-phase fluid model and numerical scheme for economic LWR system simulation

SciTech Connect

Ransom, V.H.; Wagner, R.J.; Trapp, J.A.

1981-01-01

The RELAP5 two-phase fluid model and the associated numerical scheme are summarized. The experience accrued in development of a fast running light water reactor system transient analysis code is reviewed and example of the code application are given.

SciTech Connect

Giorla, Alain B.

2015-03-01

In this report, we establish a numerical model for concrete exposed to irradiation to address these three critical points. The model accounts for creep in the cement paste and its coupling with damage, temperature and relative humidity. The shift in failure mode with the loading rate is also properly represented. The numerical model for creep has been validated and calibrated against different experiments in the literature [Wittmann, 1970, Le Roy, 1995]. Results from a simplified model are shown to showcase the ability of numerical homogenization to simulate irradiation effects in concrete. In future works, the complete model will be applied to the analysis of the irradiation experiments of Elleuch et al. [1972] and Kelly et al. [1969]. This requires a careful examination of the experimental environmental conditions as in both cases certain critical information are missing, including the relative humidity history. A sensitivity analysis will be conducted to provide lower and upper bounds of the concrete expansion under irradiation, and check if the scatter in the simulated results matches the one found in experiments. The numerical and experimental results will be compared in terms of expansion and loss of mechanical stiffness and strength. Both effects should be captured accordingly by the model to validate it. Once the model has been validated on these two experiments, it can be applied to simulate concrete from nuclear power plants. To do so, the materials used in these concrete must be as well characterized as possible. The main parameters required are the mechanical properties of each constituent in the concrete (aggregates, cement paste), namely the elastic modulus, the creep properties, the tensile and compressive strength, the thermal expansion coefficient, and the drying shrinkage. These can be either measured experimentally, estimated from the initial composition in the case of cement paste, or back-calculated from mechanical tests on concrete. If some

8. Advanced in turbulence physics and modeling by direct numerical simulations

NASA Technical Reports Server (NTRS)

Reynolds, W. C.

1987-01-01

The advent of direct numerical simulations of turbulence has opened avenues for research on turbulence physics and turbulence modeling. Direct numerical simulation provides values for anything that the scientist or modeler would like to know about the flow. An overview of some recent advances in the physical understanding of turbulence and in turbulence modeling obtained through such simulations is presented.

9. Numerical calculation of tidal current with UTOPIA scheme for advection and application to Osaka Bay

Komoda, Jun; Matsuyama, Masaji

UTOPIA scheme was applied to advection term for the numerical calculation of tide and tidal current to reproduce the strong tidal current realistically. Numerical model is constructed by boundary-fitted coordinate method vertically using Arakawa A grid in space. The new method is designed to suppress a numerical oscillation usually induced by Arakawa A grid. UTOPIA scheme was confirmed to be suitable to express a strong current around complicated topography. This model was applied to the tidal calculation for M2 constituent in Osaka Bay with two narrow straits, i.e., Akashi and Tomogashima straits. The tidal currents obtained in this model agree with them observed at monitoring stations, and the four eddies in the bay were also reproduced as the residual currents, i.e., tide induced transient eddy (TITE). The generation, growth and lifetime of the eddies also were investigated.

10. Numerical pricing of options using high-order compact finite difference schemes

Tangman, D. Y.; Gopaul, A.; Bhuruth, M.

2008-09-01

We consider high-order compact (HOC) schemes for quasilinear parabolic partial differential equations to discretise the Black-Scholes PDE for the numerical pricing of European and American options. We show that for the heat equation with smooth initial conditions, the HOC schemes attain clear fourth-order convergence but fail if non-smooth payoff conditions are used. To restore the fourth-order convergence, we use a grid stretching that concentrates grid nodes at the strike price for European options. For an American option, an efficient procedure is also described to compute the option price, Greeks and the optimal exercise curve. Comparisons with a fourth-order non-compact scheme are also done. However, fourth-order convergence is not experienced with this strategy. To improve the convergence rate for American options, we discuss the use of a front-fixing transformation with the HOC scheme. We also show that the HOC scheme with grid stretching along the asset price dimension gives accurate numerical solutions for European options under stochastic volatility.

11. 3 Lectures: "Lagrangian Models", "Numerical Transport Schemes", and "Chemical and Transport Models"

NASA Technical Reports Server (NTRS)

Douglass, A.

2005-01-01

The topics for the three lectures for the Canadian Summer School are Lagrangian Models, numerical transport schemes, and chemical and transport models. In the first lecture I will explain the basic components of the Lagrangian model (a trajectory code and a photochemical code), the difficulties in using such a model (initialization) and show some applications in interpretation of aircraft and satellite data. If time permits I will show some results concerning inverse modeling which is being used to evaluate sources of tropospheric pollutants. In the second lecture I will discuss one of the core components of any grid point model, the numerical transport scheme. I will explain the basics of shock capturing schemes, and performance criteria. I will include an example of the importance of horizontal resolution to polar processes. We have learned from NASA's global modeling initiative that horizontal resolution matters for predictions of the future evolution of the ozone hole. The numerical scheme will be evaluated using performance metrics based on satellite observations of long-lived tracers. The final lecture will discuss the evolution of chemical transport models over the last decade. Some of the problems with assimilated winds will be demonstrated, using satellite data to evaluate the simulations.

12. Numerical schemes for dynamically orthogonal equations of stochastic fluid and ocean flows

SciTech Connect

Ueckermann, M.P.; Lermusiaux, P.F.J.; Sapsis, T.P.

2013-01-15

The quantification of uncertainties is critical when systems are nonlinear and have uncertain terms in their governing equations or are constrained by limited knowledge of initial and boundary conditions. Such situations are common in multiscale, intermittent and non-homogeneous fluid and ocean flows. The dynamically orthogonal (DO) field equations provide an adaptive methodology to predict the probability density functions of such flows. The present work derives efficient computational schemes for the DO methodology applied to unsteady stochastic Navier-Stokes and Boussinesq equations, and illustrates and studies the numerical aspects of these schemes. Semi-implicit projection methods are developed for the mean and for the DO modes, and time-marching schemes of first to fourth order are used for the stochastic coefficients. Conservative second-order finite-volumes are employed in physical space with new advection schemes based on total variation diminishing methods. Other results include: (i) the definition of pseudo-stochastic pressures to obtain a number of pressure equations that is linear in the subspace size instead of quadratic; (ii) symmetric advection schemes for the stochastic velocities; (iii) the use of generalized inversion to deal with singular subspace covariances or deterministic modes; and (iv) schemes to maintain orthonormal modes at the numerical level. To verify our implementation and study the properties of our schemes and their variations, a set of stochastic flow benchmarks are defined including asymmetric Dirac and symmetric lock-exchange flows, lid-driven cavity flows, and flows past objects in a confined channel. Different Reynolds number and Grashof number regimes are employed to illustrate robustness. Optimal convergence under both time and space refinements is shown as well as the convergence of the probability density functions with the number of stochastic realizations.

13. A New Framework to Compare Mass-Flux Schemes Within the AROME Numerical Weather Prediction Model

Riette, Sébastien; Lac, Christine

2016-08-01

In the Application of Research to Operations at Mesoscale (AROME) numerical weather forecast model used in operations at Météo-France, five mass-flux schemes are available to parametrize shallow convection at kilometre resolution. All but one are based on the eddy-diffusivity-mass-flux approach, and differ in entrainment/detrainment, the updraft vertical velocity equation and the closure assumption. The fifth is based on a more classical mass-flux approach. Screen-level scores obtained with these schemes show few discrepancies and are not sufficient to highlight behaviour differences. Here, we describe and use a new experimental framework, able to compare and discriminate among different schemes. For a year, daily forecast experiments were conducted over small domains centred on the five French metropolitan radio-sounding locations. Cloud base, planetary boundary-layer height and normalized vertical profiles of specific humidity, potential temperature, wind speed and cloud condensate were compared with observations, and with each other. The framework allowed the behaviour of the different schemes in and above the boundary layer to be characterized. In particular, the impact of the entrainment/detrainment formulation, closure assumption and cloud scheme were clearly visible. Differences mainly concerned the transport intensity thus allowing schemes to be separated into two groups, with stronger or weaker updrafts. In the AROME model (with all interactions and the possible existence of compensating errors), evaluation diagnostics gave the advantage to the first group.

14. A New Framework to Compare Mass-Flux Schemes Within the AROME Numerical Weather Prediction Model

Riette, Sébastien; Lac, Christine

2016-03-01

In the Application of Research to Operations at Mesoscale (AROME) numerical weather forecast model used in operations at Météo-France, five mass-flux schemes are available to parametrize shallow convection at kilometre resolution. All but one are based on the eddy-diffusivity-mass-flux approach, and differ in entrainment/detrainment, the updraft vertical velocity equation and the closure assumption. The fifth is based on a more classical mass-flux approach. Screen-level scores obtained with these schemes show few discrepancies and are not sufficient to highlight behaviour differences. Here, we describe and use a new experimental framework, able to compare and discriminate among different schemes. For a year, daily forecast experiments were conducted over small domains centred on the five French metropolitan radio-sounding locations. Cloud base, planetary boundary-layer height and normalized vertical profiles of specific humidity, potential temperature, wind speed and cloud condensate were compared with observations, and with each other. The framework allowed the behaviour of the different schemes in and above the boundary layer to be characterized. In particular, the impact of the entrainment/detrainment formulation, closure assumption and cloud scheme were clearly visible. Differences mainly concerned the transport intensity thus allowing schemes to be separated into two groups, with stronger or weaker updrafts. In the AROME model (with all interactions and the possible existence of compensating errors), evaluation diagnostics gave the advantage to the first group.

15. Numerical study of fourth-order linearized compact schemes for generalized NLS equations

Liao, Hong-lin; Shi, Han-sheng; Zhao, Ying

2014-08-01

The fourth-order compact approximation for the spatial second-derivative and several linearized approaches, including the time-lagging method of Zhang et al. (1995), the local-extrapolation technique of Chang et al. (1999) and the recent scheme of Dahlby et al. (2009), are considered in constructing fourth-order linearized compact difference (FLCD) schemes for generalized NLS equations. By applying a new time-lagging linearized approach, we propose a symmetric fourth-order linearized compact difference (SFLCD) scheme, which is shown to be more robust in long-time simulations of plane wave, breather, periodic traveling-wave and solitary wave solutions. Numerical experiments suggest that the SFLCD scheme is a little more accurate than some other FLCD schemes and the split-step compact difference scheme of Dehghan and Taleei (2010). Compared with the time-splitting pseudospectral method of Bao et al. (2003), our SFLCD method is more suitable for oscillating solutions or the problems with a rapidly varying potential.

16. Efficient Schemes for Reducing Numerical Dispersion in ModelingMultiphase Transport through Porous and Fractured Media

SciTech Connect

Wu, Yu-Shu; Forsyth, Peter A.

2006-04-13

Numerical issues with modeling transport of chemicals or solute in realistic large-scale subsurface systems have been a serious concern, even with the continual progress made in both simulation algorithms and computer hardware in the past few decades. The problem remains and becomes even more difficult when dealing with chemical transport in a multiphase flow system using coarse, multidimensional regular or irregular grids, because of the known effects of numerical dispersion associated with moving plume fronts. We have investigated several total-variation-diminishing (TVD) or flux-limiter schemes by implementing and testing them in the T2R3D code, one of the TOUGH2 family of codes. The objectives of this paper are (1) to investigate the possibility of applying these TVD schemes, using multi-dimensional irregular unstructured grids, and (2) to help select more accurate spatial averaging methods for simulating chemical transport given a numerical grid or spatial discretization. We present an application example to show that such TVD schemes are able to effectively reduce numerical dispersion.

17. Numerical scheme for riser motion calculation during 3-D VIV simulation

Huang, Kevin; Chen, Hamn-Ching; Chen, Chia-Rong

2011-10-01

This paper presents a numerical scheme for riser motion calculation and its application to riser VIV simulations. The discretisation of the governing differential equation is studied first. The top tensioned risers are simplified as tensioned beams. A centered space and forward time finite difference scheme is derived from the governing equations of motion. Then an implicit method is adopted for better numerical stability. The method meets von Neumann criteria and is shown to be unconditionally stable. The discretized linear algebraic equations are solved using a LU decomposition method. This approach is then applied to a series of benchmark cases with known solutions. The comparisons show good agreement. Finally the method is applied to practical riser VIV simulations. The studied cases cover a wide range of riser VIV problems, i.e. different riser outer diameter, length, tensioning conditions, and current profiles. Reasonable agreement is obtained between the numerical simulations and experimental data on riser motions and cross-flow VIV a/D . These validations and comparisons confirm that the present numerical scheme for riser motion calculation is valid and effective for long riser VIV simulation.

18. Numerical dissipation control in high order shock-capturing schemes for LES of low speed flows

Kotov, D. V.; Yee, H. C.; Wray, A. A.; Sjögreen, B.; Kritsuk, A. G.

2016-02-01

The Yee & Sjögreen adaptive numerical dissipation control in high order scheme (High Order Filter Methods for Wide Range of Compressible Flow Speeds, ICOSAHOM 09, 2009) is further improved for DNS and LES of shock-free turbulence and low speed turbulence with shocklets. There are vastly different requirements in the minimization of numerical dissipation for accurate turbulence simulations of different compressible flow types and flow speeds. Traditionally, the method of choice for shock-free turbulence and low speed turbulence are by spectral, high order central or high order compact schemes with high order linear filters. With a proper control of a local flow sensor, appropriate amount of numerical dissipation in high order shock-capturing schemes can have spectral-like accuracy for compressible low speed turbulent flows. The development of the method includes an adaptive flow sensor with automatic selection on the amount of numerical dissipation needed at each flow location for more accurate DNS and LES simulations with less tuning of parameters for flows with a wide range of flow speed regime during the time-accurate evolution, e.g., time varying random forcing. An automatic selection of the different flow sensors catered to the different flow types is constructed. A Mach curve and high-frequency oscillation indicators are used to reduce the tuning of parameters in controlling the amount of shock-capturing numerical dissipation to be employed for shock-free turbulence, low speed turbulence and turbulence with strong shocks. In Kotov et al. (High Order Numerical Methods for LES of Turbulent Flows with Shocks, ICCFD8, Chengdu, Sichuan, China, July 14-18, 2014) the LES of a turbulent flow with a strong shock by the Yee & Sjögreen scheme indicated a good agreement with the filtered DNS data. A work in progress for the application of the adaptive flow sensor for compressible turbulence with time-varying random forcing is forthcoming. The present study examines the

19. Adaptive Numerical Dissipation Control in High Order Schemes for Multi-D Non-Ideal MHD

NASA Technical Reports Server (NTRS)

Yee, H. C.; Sjoegreen, B.

2005-01-01

The required type and amount of numerical dissipation/filter to accurately resolve all relevant multiscales of complex MHD unsteady high-speed shock/shear/turbulence/combustion problems are not only physical problem dependent, but also vary from one flow region to another. In addition, proper and efficient control of the divergence of the magnetic field (Div(B)) numerical error for high order shock-capturing methods poses extra requirements for the considered type of CPU intensive computations. The goal is to extend our adaptive numerical dissipation control in high order filter schemes and our new divergence-free methods for ideal MHD to non-ideal MHD that include viscosity and resistivity. The key idea consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and leave the rest of the region free from numerical dissipation contamination. These scheme-independent detectors are capable of distinguishing shocks/shears, flame sheets, turbulent fluctuations and spurious high-frequency oscillations. The detection algorithm is based on an artificial compression method (ACM) (for shocks/shears), and redundant multiresolution wavelets (WAV) (for the above types of flow feature). These filters also provide a natural and efficient way for the minimization of Div(B) numerical error.

20. Higher-Order Compact Schemes for Numerical Simulation of Incompressible Flows

NASA Technical Reports Server (NTRS)

Wilson, Robert V.; Demuren, Ayodeji O.; Carpenter, Mark

1998-01-01

A higher order accurate numerical procedure has been developed for solving incompressible Navier-Stokes equations for 2D or 3D fluid flow problems. It is based on low-storage Runge-Kutta schemes for temporal discretization and fourth and sixth order compact finite-difference schemes for spatial discretization. The particular difficulty of satisfying the divergence-free velocity field required in incompressible fluid flow is resolved by solving a Poisson equation for pressure. It is demonstrated that for consistent global accuracy, it is necessary to employ the same order of accuracy in the discretization of the Poisson equation. Special care is also required to achieve the formal temporal accuracy of the Runge-Kutta schemes. The accuracy of the present procedure is demonstrated by application to several pertinent benchmark problems.

ERIC Educational Resources Information Center

Sparks, Alan; Ingram, Hadyn; Phillips, Sunny

2009-01-01

Purpose: The purpose of this paper is to evaluate an innovative way to train adult apprentices for the construction industry. Design/methodology/approach: The paper emphasizes that, in order to address skills shortages for international construction, training methods must be improved. It looks at the example of an adult apprenticeship scheme in…

2. A numerical scheme for optimal transition paths of stochastic chemical kinetic systems

Liu, Di

2008-10-01

We present a new framework for finding the optimal transition paths of metastable stochastic chemical kinetic systems with large system size. The optimal transition paths are identified to be the most probable paths according to the Large Deviation Theory of stochastic processes. Dynamical equations for the optimal transition paths are derived using the variational principle. A modified Minimum Action Method (MAM) is proposed as a numerical scheme to solve the optimal transition paths. Applications to Gene Regulatory Networks such as the toggle switch model and the Lactose Operon Model in Escherichia coli are presented as numerical examples.

3. Numerical integration of the stochastic Landau-Lifshitz-Gilbert equation in generic time-discretization schemes.

PubMed

Romá, Federico; Cugliandolo, Leticia F; Lozano, Gustavo S

2014-08-01

We introduce a numerical method to integrate the stochastic Landau-Lifshitz-Gilbert equation in spherical coordinates for generic discretization schemes. This method conserves the magnetization modulus and ensures the approach to equilibrium under the expected conditions. We test the algorithm on a benchmark problem: the dynamics of a uniformly magnetized ellipsoid. We investigate the influence of various parameters, and in particular, we analyze the efficiency of the numerical integration, in terms of the number of steps needed to reach a chosen long time with a given accuracy. PMID:25215839

4. Numerical Simulations and Optimisation in Forming of Advanced Materials

Huétink, J.

2007-04-01

With the introduction of new materials as high strength steels, metastable steels and fiber reinforce composites, the need for advanced physically valid constitutive models arises. A biaxial test equipment is developed and applied for the determination of material data as well as for validation of material models. An adaptive through- thickness integration scheme for plate elements is developed, which improves the accuracy of spring back prediction at minimal costs. An optimization strategy is proposed that assists an engineer to model an optimization problem.

5. Numerical investigation of a modified family of centered schemes applied to multiphase equations with nonconservative sources

Crochet, M. W.; Gonthier, K. A.

2013-12-01

Systems of hyperbolic partial differential equations are frequently used to model the flow of multiphase mixtures. These equations often contain sources, referred to as nozzling terms, that cannot be posed in divergence form, and have proven to be particularly challenging in the development of finite-volume methods. Upwind schemes have recently shown promise in properly resolving the steady wave solution of the associated multiphase Riemann problem. However, these methods require a full characteristic decomposition of the system eigenstructure, which may be either unavailable or computationally expensive. Central schemes, such as the Kurganov-Tadmor (KT) family of methods, require minimal characteristic information, which makes them easily applicable to systems with an arbitrary number of phases. However, the proper implementation of nozzling terms in these schemes has been mathematically ambiguous. The primary objectives of this work are twofold: first, an extension of the KT family of schemes is proposed that formally accounts for the nonconservative nozzling sources. This modification results in a semidiscrete form that retains the simplicity of its predecessor and introduces little additional computational expense. Second, this modified method is applied to multiple, but equivalent, forms of the multiphase equations to perform a numerical study by solving several one-dimensional test problems. Both ideal and Mie-Grüneisen equations of state are used, with the results compared to an analytical solution. This study demonstrates that the magnitudes of the resulting numerical errors are sensitive to the form of the equations considered, and suggests an optimal form to minimize these errors. Finally, a separate modification of the wave propagation speeds used in the KT family is also suggested that can reduce the extent of numerical diffusion in multiphase flows.

6. Numerical Modeling of Deep Mantle Convection: Advection and Diffusion Schemes for Marker Methods

Mulyukova, Elvira; Dabrowski, Marcin; Steinberger, Bernhard

2013-04-01

7. Advances in Numerical Boundary Conditions for Computational Aeroacoustics

NASA Technical Reports Server (NTRS)

Tam, Christopher K. W.

1997-01-01

Advances in Computational Aeroacoustics (CAA) depend critically on the availability of accurate, nondispersive, least dissipative computation algorithm as well as high quality numerical boundary treatments. This paper focuses on the recent developments of numerical boundary conditions. In a typical CAA problem, one often encounters two types of boundaries. Because a finite computation domain is used, there are external boundaries. On the external boundaries, boundary conditions simulating the solution outside the computation domain are to be imposed. Inside the computation domain, there may be internal boundaries. On these internal boundaries, boundary conditions simulating the presence of an object or surface with specific acoustic characteristics are to be applied. Numerical boundary conditions, both external or internal, developed for simple model problems are reviewed and examined. Numerical boundary conditions for real aeroacoustic problems are also discussed through specific examples. The paper concludes with a description of some much needed research in numerical boundary conditions for CAA.

8. Automatic optimization of metrology sampling scheme for advanced process control

Chue, Chuei-Fu; Huang, Chun-Yen; Shih, Chiang-Lin

2011-03-01

In order to ensure long-term profitability, driving the operational costs down and improving the yield of a DRAM manufacturing process are continuous efforts. This includes optimal utilization of the capital equipment. The costs of metrology needed to ensure yield are contributing to the overall costs. As the shrinking of device dimensions continues, the costs of metrology are increasing because of the associated tightening of the on-product specifications requiring more metrology effort. The cost-of-ownership reduction is tackled by increasing the throughput and availability of metrology systems. However, this is not the only way to reduce metrology effort. In this paper, we discuss how the costs of metrology can be improved by optimizing the recipes in terms of the sampling layout, thereby eliminating metrology that does not contribute to yield. We discuss results of sampling scheme optimization for on-product overlay control of two DRAM manufacturing processes at Nanya Technology Corporation. For a 6x DRAM production process, we show that the reduction of metrology waste can be as high as 27% and overlay can be improved by 36%, comparing with a baseline sampling scheme. For a 4x DRAM process, having tighter overlay specs, a gain of ca. 0.5nm on-product overlay could be achieved, without increasing the metrology effort relative to the original sampling plan.

9. Building fast well-balanced two-stage numerical schemes for a model of two-phase flows

Thanh, Mai Duc

2014-06-01

We present a set of well-balanced two-stage schemes for an isentropic model of two-phase flows arisen from the modeling of deflagration-to-detonation transition in granular materials. The first stage is to absorb the source term in nonconservative form into equilibria. Then in the second stage, these equilibria will be composed into a numerical flux formed by using a convex combination of the numerical flux of a stable Lax-Friedrichs-type scheme and the one of a higher-order Richtmyer-type scheme. Numerical schemes constructed in such a way are expected to get the interesting property: they are fast and stable. Tests show that the method works out until the parameter takes on the value CFL, and so any value of the parameter between zero and this value is expected to work as well. All the schemes in this family are shown to capture stationary waves and preserves the positivity of the volume fractions. The special values of the parameter 0,1/2,1/(1+CFL), and CFL in this family define the Lax-Friedrichs-type, FAST1, FAST2, and FAST3 schemes, respectively. These schemes are shown to give a desirable accuracy. The errors and the CPU time of these schemes and the Roe-type scheme are calculated and compared. The constructed schemes are shown to be well-balanced and faster than the Roe-type scheme.

10. Efficient numerical schemes for viscoplastic avalanches. Part 1: The 1D case

SciTech Connect

Fernández-Nieto, Enrique D.

2014-05-01

This paper deals with the numerical resolution of a shallow water viscoplastic flow model. Viscoplastic materials are characterized by the existence of a yield stress: below a certain critical threshold in the imposed stress, there is no deformation and the material behaves like a rigid solid, but when that yield value is exceeded, the material flows like a fluid. In the context of avalanches, it means that after going down a slope, the material can stop and its free surface has a non-trivial shape, as opposed to the case of water (Newtonian fluid). The model involves variational inequalities associated with the yield threshold: finite-volume schemes are used together with duality methods (namely Augmented Lagrangian and Bermúdez–Moreno) to discretize the problem. To be able to accurately simulate the stopping behavior of the avalanche, new schemes need to be designed, involving the classical notion of well-balancing. In the present context, it needs to be extended to take into account the viscoplastic nature of the material as well as general bottoms with wet/dry fronts which are encountered in geophysical geometries. We derived such schemes and numerical experiments are presented to show their performances.

11. Toward a consistent framework for high order mesh refinement schemes in numerical relativity

Mongwane, Bishop

2015-05-01

It has now become customary in the field of numerical relativity to couple high order finite difference schemes to mesh refinement algorithms. To this end, different modifications to the standard Berger-Oliger adaptive mesh refinement algorithm have been proposed. In this work we present a fourth order stable mesh refinement scheme with sub-cycling in time for numerical relativity. We do not use buffer zones to deal with refinement boundaries but explicitly specify boundary data for refined grids. We argue that the incompatibility of the standard mesh refinement algorithm with higher order Runge Kutta methods is a manifestation of order reduction phenomena, caused by inconsistent application of boundary data in the refined grids. Our scheme also addresses the problem of spurious reflections that are generated when propagating waves cross mesh refinement boundaries. We introduce a transition zone on refined levels within which the phase velocity of propagating modes is allowed to decelerate in order to smoothly match the phase velocity of coarser grids. We apply the method to test problems involving propagating waves and show a significant reduction in spurious reflections.

12. von Neumann Stability Analysis of Numerical Solution Schemes for 1D and 2D Euler Equations

Konangi, Santosh; Palakurthi, Nikhil Kumar; Ghia, Urmila

2014-11-01

A von Neumann stability analysis is conducted for numerical schemes for the full system of coupled, density-based 1D and 2D Euler equations, closed by an isentropic equation of state. The governing equations are discretized on a staggered grid, which permits equivalence to finite-volume discretization. Presently, first-order accurate spatial and temporal finite-difference techniques are analyzed. The momentum convection term is treated as explicit, semi-implicit or implicit. Density upwind bias is included in the spatial operator of the continuity equation. By combining the discretization techniques, ten solution schemes are formulated. For each scheme, unstable and stable regimes are identified through the stability analysis, and the maximum allowable CFL number is predicted. The predictions are verified for selected schemes, using the Riemann problem at incompressible and compressible Mach numbers. Very good agreement is obtained between the analytically predicted and ``experimentally'' observed CFL values for all cases, thereby validating the analysis. The demonstrated analysis provides an accurate indication of stability conditions for the Euler equations, in contrast to the simplistic conditions arising from model equations, such as the wave equation.

13. On the validation of cloud parametrization schemes in numerical atmospheric models with satellite data from ISCCP

Meinke, I.

2003-04-01

A new method is presented to validate cloud parametrization schemes in numerical atmospheric models with satellite data of scanning radiometers. This method is applied to the regional atmospheric model HRM (High Resolution Regional Model) using satellite data from ISCCP (International Satellite Cloud Climatology Project). Due to the limited reliability of former validations there has been a need for developing a new validation method: Up to now differences between simulated and measured cloud properties are mostly declared as deficiencies of the cloud parametrization scheme without further investigation. Other uncertainties connected with the model or with the measurements have not been taken into account. Therefore changes in the cloud parametrization scheme based on such kind of validations might not be realistic. The new method estimates uncertainties of the model and the measurements. Criteria for comparisons of simulated and measured data are derived to localize deficiencies in the model. For a better specification of these deficiencies simulated clouds are classified regarding their parametrization. With this classification the localized model deficiencies are allocated to a certain parametrization scheme. Applying this method to the regional model HRM the quality of forecasting cloud properties is estimated in detail. The overestimation of simulated clouds in low emissivity heights especially during the night is localized as model deficiency. This is caused by subscale cloudiness. As the simulation of subscale clouds in the regional model HRM is described by a relative humidity parametrization these deficiencies are connected with this parameterization.

14. Advanced numerical methods and software approaches for semiconductor device simulation

SciTech Connect

CAREY,GRAHAM F.; PARDHANANI,A.L.; BOVA,STEVEN W.

2000-03-23

In this article the authors concisely present several modern strategies that are applicable to drift-dominated carrier transport in higher-order deterministic models such as the drift-diffusion, hydrodynamic, and quantum hydrodynamic systems. The approaches include extensions of upwind and artificial dissipation schemes, generalization of the traditional Scharfetter-Gummel approach, Petrov-Galerkin and streamline-upwind Petrov Galerkin (SUPG), entropy variables, transformations, least-squares mixed methods and other stabilized Galerkin schemes such as Galerkin least squares and discontinuous Galerkin schemes. The treatment is representative rather than an exhaustive review and several schemes are mentioned only briefly with appropriate reference to the literature. Some of the methods have been applied to the semiconductor device problem while others are still in the early stages of development for this class of applications. They have included numerical examples from the recent research tests with some of the methods. A second aspect of the work deals with algorithms that employ unstructured grids in conjunction with adaptive refinement strategies. The full benefits of such approaches have not yet been developed in this application area and they emphasize the need for further work on analysis, data structures and software to support adaptivity. Finally, they briefly consider some aspects of software frameworks. These include dial-an-operator approaches such as that used in the industrial simulator PROPHET, and object-oriented software support such as those in the SANDIA National Laboratory framework SIERRA.

15. Advanced Numerical Methods and Software Approaches for Semiconductor Device Simulation

DOE PAGESBeta

Carey, Graham F.; Pardhanani, A. L.; Bova, S. W.

2000-01-01

In this article we concisely present several modern strategies that are applicable to driftdominated carrier transport in higher-order deterministic models such as the driftdiffusion, hydrodynamic, and quantum hydrodynamic systems. The approaches include extensions of “upwind” and artificial dissipation schemes, generalization of the traditional Scharfetter – Gummel approach, Petrov – Galerkin and streamline-upwind Petrov Galerkin (SUPG), “entropy” variables, transformations, least-squares mixed methods and other stabilized Galerkin schemes such as Galerkin least squares and discontinuous Galerkin schemes. The treatment is representative rather than an exhaustive review and several schemes are mentioned only briefly with appropriate reference to the literature. Some of themore » methods have been applied to the semiconductor device problem while others are still in the early stages of development for this class of applications. We have included numerical examples from our recent research tests with some of the methods. A second aspect of the work deals with algorithms that employ unstructured grids in conjunction with adaptive refinement strategies. The full benefits of such approaches have not yet been developed in this application area and we emphasize the need for further work on analysis, data structures and software to support adaptivity. Finally, we briefly consider some aspects of software frameworks. These include dial-an-operator approaches such as that used in the industrial simulator PROPHET, and object-oriented software support such as those in the SANDIA National Laboratory framework SIERRA.« less

16. Multiecho scheme advances surface NMR for aquifer characterization

Grunewald, Elliot; Walsh, David

2013-12-01

nuclear magnetic resonance (NMR) is increasingly used as a method to noninvasively characterize aquifers. This technology follows a successful history of NMR logging, applied over decades to estimate hydrocarbon reservoir properties. In contrast to logging, however, surface methods have utilized relatively simple acquisition sequences, from which pore-scale properties may not be reliably and efficiently estimated. We demonstrate for the first time the capability of sophisticated multiecho measurements to rapidly record a surface NMR response that more directly reflects aquifer characteristics. Specifically, we develop an adaptation of the multipulse Carr-Purcell-Meiboom-Gill (CPMG) sequence, widely used in logging, to measure the T2 relaxation response in a single scan. We validate this approach in a field surface NMR data set and by direct comparison with an NMR log. Adoption of the CPMG marked a landmark advancement in the history of logging NMR; we have now realized this same advancement in the surface NMR method.

17. Application of an efficient hybrid scheme for aeroelastic analysis of advanced propellers

NASA Technical Reports Server (NTRS)

Srivastava, R.; Sankar, N. L.; Reddy, T. S. R.; Huff, D. L.

1989-01-01

An efficient 3-D hybrid scheme is applied for solving Euler equations to analyze advanced propellers. The scheme treats the spanwise direction semi-explicitly and the other two directions implicitly, without affecting the accuracy, as compared to a fully implicit scheme. This leads to a reduction in computer time and memory requirement. The calculated power coefficients for two advanced propellers, SR3 and SR7L, and various advanced ratios showed good correlation with experiment. Spanwise distribution of elemental power coefficient and steady pressure coefficient differences also showed good agreement with experiment. A study of the effect of structural flexibility on the performance of the advanced propellers showed that structural deformation due to centrifugal and aero loading should be included for better correlation.

18. A numerical study of a class of TVD schemes for compressible mixing layers

NASA Technical Reports Server (NTRS)

Sandham, N. D.; Yee, H. C.

1989-01-01

At high Mach numbers the two-dimensional time-developing mixing layer develops shock waves, positioned around large-scale vortical structures. A suitable numerical method has to be able to capture the inherent instability of the flow, leading to the roll-up of vortices, and also must be able to capture shock waves when they develop. Standard schemes for low speed turbulent flows, for example spectral methods, rely on resolution of all flow-features and cannot handle shock waves, which become too thin at any realistic Reynolds number. The performance of a class of second-order explicit total variation diminishing (TVD) schemes on a compressible mixing layer problem was studied. The basic idea is to capture the physics of the flow correctly, by resolving down to the smallest turbulent length scales, without resorting to turbulence or sub-grid scale modeling, and at the same time capture shock waves without spurious oscillations. The present study indicates that TVD schemes can capture the shocks accurately when they form, but (without resorting to a finer grid) have poor accuracy in computing the vortex growth. The solution accuracy depends on the choice of limiter. However a larger number of grid points are in general required to resolve the correct vortex growth. The low accuracy in computing time-dependent problems containing shock waves as well as vortical structures is partly due to the inherent shock-capturing property of all TVD schemes. In order to capture shock waves without spurious oscillations these schemes reduce to first-order near extrema and indirectly produce clipping phenomena, leading to inaccuracy in the computation of vortex growth. Accurate simulation of unsteady turbulent fluid flows with shock waves will require further development of efficient, uniformly higher than second-order accurate, shock-capturing methods.

19. Numerical Schemes for the Hamilton-Jacobi and Level Set Equations on Triangulated Domains

NASA Technical Reports Server (NTRS)

Barth, Timothy J.; Sethian, James A.

2006-01-01

Borrowing from techniques developed for conservation law equations, we have developed both monotone and higher order accurate numerical schemes which discretize the Hamilton-Jacobi and level set equations on triangulated domains. The use of unstructured meshes containing triangles (2D) and tetrahedra (3D) easily accommodates mesh adaptation to resolve disparate level set feature scales with a minimal number of solution unknowns. The minisymposium talk will discuss these algorithmic developments and present sample calculations using our adaptive triangulation algorithm applied to various moving interface problems such as etching, deposition, and curvature flow.

20. Thermoplastic polymer patterning without residual layer by advanced nanoimprinting schemes

Park, Hyunsoo; Cheng, Xing

2009-06-01

Nanoimprinting is a fast-growing technique for nanoscale patterning. One of the remaining issues in nanoimprinting is the removal of the residual layer after nanoimprinting. Traditionally the residual layer is removed by an oxygen reactive-ion etching (RIE) step. The need for a vacuum environment and dedicated equipment in this step lowers the throughput and increases the cost of the nanoimprinting process. It also prevents the possibility of patterning isolated functional polymers because oxygen RIE destroys the functional materials. In this work, novel nanoimprinting schemes are developed to nondestructively remove the residual layer in thermal nanoimprinting by solvent developing and dewetting. Combined with a transfer-bonding technique, three-dimensional polymer scaffolds are achieved. The techniques developed here eliminate the RIE step in thermal nanoimprinting and are compatible with roller nanoimprinting for large-scale patterning of polymer micro- or nanostructures. The technique also opens up new applications for nanoimprinting in patterning isolated conjugated polymers for organic electronic devices and circuits.

1. Thermoplastic polymer patterning without residual layer by advanced nanoimprinting schemes.

PubMed

Park, Hyunsoo; Cheng, Xing

2009-06-17

Nanoimprinting is a fast-growing technique for nanoscale patterning. One of the remaining issues in nanoimprinting is the removal of the residual layer after nanoimprinting. Traditionally the residual layer is removed by an oxygen reactive-ion etching (RIE) step. The need for a vacuum environment and dedicated equipment in this step lowers the throughput and increases the cost of the nanoimprinting process. It also prevents the possibility of patterning isolated functional polymers because oxygen RIE destroys the functional materials. In this work, novel nanoimprinting schemes are developed to nondestructively remove the residual layer in thermal nanoimprinting by solvent developing and dewetting. Combined with a transfer-bonding technique, three-dimensional polymer scaffolds are achieved. The techniques developed here eliminate the RIE step in thermal nanoimprinting and are compatible with roller nanoimprinting for large-scale patterning of polymer micro- or nanostructures. The technique also opens up new applications for nanoimprinting in patterning isolated conjugated polymers for organic electronic devices and circuits. PMID:19468173

2. A hybrid numerical prediction scheme for solar radiation estimation in un-gauged catchments.

Shamim, M. A.; Bray, M.; Ishak, A. M.; Remesan, R.; Han, D.

2009-09-01

The importance of solar radiation on earth's surface is depicted in its wide range of applications in the fields of meteorology, agricultural sciences, engineering, hydrology, crop water requirements, climatic changes and energy assessment. It is quite random in nature as it has to go through different processes of assimilation and dispersion while on its way to earth. Compared to other meteorological parameters, solar radiation is quite infrequently measured, for example, the worldwide ratio of stations collecting solar radiation to those collecting temperature is 1:500 (Badescu, 2008). Researchers, therefore, have to rely on indirect techniques of estimation that include nonlinear models, artificial intelligence (e.g. neural networks), remote sensing and numerical weather predictions (NWP). This study proposes a hybrid numerical prediction scheme for solar radiation estimation in un-gauged catchments. It uses the PSU/NCAR's Mesoscale Modelling system (MM5) (Grell et al., 1995) to parameterise the cloud effect on extraterrestrial radiation by dividing the atmosphere into four layers of very high (6-12 km), high (3-6 km), medium (1.5-3) and low (0-1.5) altitudes from earth. It is believed that various cloud forms exist within each of these layers. An hourly time series of upper air pressure and relative humidity data sets corresponding to all of these layers is determined for the Brue catchment, southwest UK, using MM5. Cloud Index (CI) was then determined using (Yang and Koike, 2002): 1 p?bi [ (Rh - Rh )] ci =------- max 0.0,---------cri dp pbi - ptipti (1- Rhcri) where, pbi and pti represent the air pressure at the top and bottom of each layer and Rhcri is the critical value of relative humidity at which a certain cloud type is formed. Output from a global clear sky solar radiation model (MRM v-5) (Kambezidis and Psiloglu, 2008) is used along with meteorological datasets of temperature and precipitation and astronomical information. The analysis is aided by the

3. Finite volume numerical scheme for high-resolution gravity field modelling and its parallel implementation

Fašková, Z.; Macák, M.; Čunderlík, R.; Mikula, K.

2012-04-01

The paper discusses a numerical solution of the geodetic boundary value problem (GBVP) by the finite volume method (FVM). The FVM is a numerical method where numerical flux is conserved from one discretization cell to its neighbour, so it's very appropriate for solving GBVP with the Neumann and the Dirichlet BCs. Our numerical scheme is developed for 3D computational domain above an ellipsoid. It is shown that a refinement of the discretization in height's direction leads to more precise numerical results. In order to achieve high-resolution numerical results, parallel implementations of algorithms using the MPI procedures were developed and computations on parallel computers were successfully performed. This basis includes the splitting of all arrays in meridian's direction, usage of an implementation of the Bi-CGSTAB non-stationary iterative solver instead of the standard SOR and an optimization of communications on parallel computers with the NUMA architecture. This gives us higher speed up in comparison to standard approaches and enables us to develop an efficient tool for high-resolution global or regional gravity field modelling in huge areas. Numerical experiments present global modelling with the resolution comparable with EGM2008 and detailed regional modelling in the Pacific Ocean with the resolution 2x2 arc min. Input gravity disturbances are generated from the DTU10-GRAV gravity field model and the disturbing potential is computed from the GOCE_DIR2 satellite geopotential model up to degree 240. Finally, the obtained disturbing potential is used to evaluate the geopotential on the DTU10 mean sea surface and the achieved mean dynamic topography is compared with the ECCO oceanographic model.

4. Advanced scheme for high-yield laser driven nuclear reactions

Margarone, D.; Picciotto, A.; Velyhan, A.; Krasa, J.; Kucharik, M.; Mangione, A.; Szydlowsky, A.; Malinowska, A.; Bertuccio, G.; Shi, Y.; Crivellari, M.; Ullschmied, J.; Bellutti, P.; Korn, G.

2015-01-01

The use of a low contrast nanosecond laser pulse with a relatively low intensity (3  ×  1016 W cm-2) allowed the enhancing of the yield of induced nuclear reactions in advanced solid targets. In particular the ‘ultraclean’ proton-boron fusion reaction, producing energetic alpha particles without neutron generation, was chosen. A spatially well-defined layer of boron dopants in a hydrogen-enriched silicon substrate was used as a target. A combination of the specific target composition and the laser pulse temporal shape allowed the enhancing of the yield of alpha particles up to 109 per steradian. This result can be ascribed to the interaction of the long-laser pre-pulse with the target and to the optimal target geometry and composition.

5. Advanced receiver autonomous integrity monitoring (ARAIM) schemes with GNSS time offsets

Wu, Yun; Wang, Jinling; Jiang, Yiping

2013-07-01

Within the current Advanced Receiver Integrity Monitoring (ARAIM) scheme, the time offsets between different Global Navigation Satellite System (GNSS) constellations are estimated along with a position solution and the GNSS receiver clock error. This scheme is called the Time-offsets Estimated ARAIM, or the TOE ARAIM. In order to enhance the interoperability and compatibility between different constellations, the time offsets are expected to be broadcast to users in future multi-GNSS positioning and navigation applications. This paper describes two new ARAIM schemes to make use of the Broadcast Time Offsets (BTOs): Time Offset Observed ARAIM (TOO ARAIM) and Time Offset Synchronized ARAIM (TOS ARAIM). It has been shown that the VPL performances of these two new ARAIM schemes rely strongly on the accuracy of BTOs. By varying the error model of the BTOs, the simulation results also demonstrate that the proposed new TOO ARAIM scheme can outperform the existing TOE ARAIM scheme-even if the accuracy of BTOs for integrity is degraded to 4.5 m and the probability of a BTO fault is relaxed to 10-2 h-1. In addition, the new Time Offset Synchronized ARAIM scheme (TOS ARAIM) can also perform better than the existing TOE ARAIM scheme if the accuracy of BTOs for integrity can reach 0.75 m. As the TOO ARAIM also has a very relaxed requirement on BTOs and better VPL performance, the TOO ARAIM is regarded as a superior ARAIM scheme for multi-GNSS with BTOs available.

6. Numerical simulation of transonic limit cycle oscillations using high-order low-diffusion schemes

Wang, Baoyuan; Zha, Ge-Cheng

2010-05-01

This paper simulates the NLR7301 airfoil limit cycle oscillation (LCO) caused by fluid-structure interaction (FSI) using Reynolds averaged Navier-Stokes equations (RANS) coupled with Spalart-Allmaras (S-A) one-equation turbulence model. A low diffusion E-CUSP (LDE) scheme with 5th order weighted essentially nonoscillatory scheme (WENO) is employed to calculate the inviscid fluxes. A fully conservative 4th order central differencing is used for the viscous terms. A fully coupled fluid-structural interaction model is employed. For the case computed in this paper, the predicted LCO frequency, amplitudes, averaged lift and moment, all agree excellently with the experiment performed by Schewe et al. The solutions appear to have bifurcation and are dependent on the initial fields or initial perturbation. The developed computational fluid dynamics (CFD)/computational structure dynamics (CSD) simulation is able to capture the LCO with very small amplitudes measured in the experiment. This is attributed to the high order low diffusion schemes, fully coupled FSI model, and the turbulence model used. This research appears to be the first time that a numerical simulation of LCO matches the experiment. The simulation confirms several observations of the experiment.

7. Orbital Advection by Interpolation: A Fast and Accurate Numerical Scheme for Super-Fast MHD Flows

SciTech Connect

Johnson, B M; Guan, X; Gammie, F

2008-04-11

In numerical models of thin astrophysical disks that use an Eulerian scheme, gas orbits supersonically through a fixed grid. As a result the timestep is sharply limited by the Courant condition. Also, because the mean flow speed with respect to the grid varies with position, the truncation error varies systematically with position. For hydrodynamic (unmagnetized) disks an algorithm called FARGO has been developed that advects the gas along its mean orbit using a separate interpolation substep. This relaxes the constraint imposed by the Courant condition, which now depends only on the peculiar velocity of the gas, and results in a truncation error that is more nearly independent of position. This paper describes a FARGO-like algorithm suitable for evolving magnetized disks. Our method is second order accurate on a smooth flow and preserves {del} {center_dot} B = 0 to machine precision. The main restriction is that B must be discretized on a staggered mesh. We give a detailed description of an implementation of the code and demonstrate that it produces the expected results on linear and nonlinear problems. We also point out how the scheme might be generalized to make the integration of other supersonic/super-fast flows more efficient. Although our scheme reduces the variation of truncation error with position, it does not eliminate it. We show that the residual position dependence leads to characteristic radial variations in the density over long integrations.

8. Multiscale/fractional step schemes for the numerical simulation of the rotating shallow water flows with complex periodic topography

Jauberteau, F.; Temam, R. M.; Tribbia, J.

2014-08-01

In this paper, we study several multiscale/fractional step schemes for the numerical solution of the rotating shallow water equations with complex topography. We consider the case of periodic boundary conditions (f-plane model). Spatial discretization is obtained using a Fourier spectral Galerkin method. For the schemes presented in this paper we consider two approaches. The first approach (multiscale schemes) is based on topography scale separation and the numerical time integration is function of the scales. The second approach is based on a splitting of the operators, and the time integration method is function of the operator considered (fractional step schemes). The numerical results obtained are compared with the explicit reference scheme (Leap-Frog scheme). With these multiscale/fractional step schemes the objective is to propose new schemes giving numerical results similar to those obtained using only one uniform fine grid N×N and a time step Δt, but with a CPU time near the CPU time needed when using only one coarse grid N1×N1, N1Δt.

9. Resource Allocation and Power Management Schemes in an LTE-Advanced Femtocell Network

Lee, Byung-Bog; Yu, Jae-Hak; Lee, In-Hwan; Pyo, Cheol-Sig; Kim, Se-Jin

In this letter, we introduce two different resource allocation and Tx power management schemes, called resource control and fixed power (RCFP) and fixed resource and power control (FRPC), in an LTE-Advanced femtocell network. We analyze and compare the two schemes in terms of the system throughput for downlink and energy consumption of home evolved NodeB (HeNB) Tx power according to the number of HeNBs and home user equipment (HUE)'s user traffic density (C). The simulation results show that the FRPC scheme has better performance in terms of system throughput for macro user equipments (MUEs) and energy consumption in low C.

10. A new hybrid-Lagrangian numerical scheme for gyrokinetic simulation of tokamak edge plasma

Ku, S.; Hager, R.; Chang, C. S.; Kwon, J. M.; Parker, S. E.

2016-06-01

In order to enable kinetic simulation of non-thermal edge plasmas at a reduced computational cost, a new hybrid-Lagrangian δf scheme has been developed that utilizes the phase space grid in addition to the usual marker particles, taking advantage of the computational strengths from both sides. The new scheme splits the particle distribution function of a kinetic equation into two parts. Marker particles contain the fast space-time varying, δf, part of the distribution function and the coarse-grained phase-space grid contains the slow space-time varying part. The coarse-grained phase-space grid reduces the memory-requirement and the computing cost, while the marker particles provide scalable computing ability for the fine-grained physics. Weights of the marker particles are determined by a direct weight evolution equation instead of the differential form weight evolution equations that the conventional delta-f schemes use. The particle weight can be slowly transferred to the phase space grid, thereby reducing the growth of the particle weights. The non-Lagrangian part of the kinetic equation - e.g., collision operation, ionization, charge exchange, heat-source, radiative cooling, and others - can be operated directly on the phase space grid. Deviation of the particle distribution function on the velocity grid from a Maxwellian distribution function - driven by ionization, charge exchange and wall loss - is allowed to be arbitrarily large. The numerical scheme is implemented in the gyrokinetic particle code XGC1, which specializes in simulating the tokamak edge plasma that crosses the magnetic separatrix and is in contact with the material wall.

11. Some recent advances in the numerical solution of differential equations

D'Ambrosio, Raffaele

2016-06-01

The purpose of the talk is the presentation of some recent advances in the numerical solution of differential equations, with special emphasis to reaction-diffusion problems, Hamiltonian problems and ordinary differential equations with discontinuous right-hand side. As a special case, in this short paper we focus on the solution of reaction-diffusion problems by means of special purpose numerical methods particularly adapted to the problem: indeed, following a problem oriented approach, we propose a modified method of lines based on the employ of finite differences shaped on the qualitative behavior of the solutions. Constructive issues and a brief analysis are presented, together with some numerical experiments showing the effectiveness of the approach and a comparison with existing solvers.

12. Modeling of Convective-Stratiform Precipitation Processes: Sensitivity to Partitioning Methods and Numerical Advection Schemes

NASA Technical Reports Server (NTRS)

Lang, Steve; Tao, W.-K.; Simpson, J.; Ferrier, B.; Einaudi, Franco (Technical Monitor)

2001-01-01

Six different convective-stratiform separation techniques, including a new technique that utilizes the ratio of vertical and terminal velocities, are compared and evaluated using two-dimensional numerical simulations of a tropical [Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE)] and midlatitude continental [Preliminary Regional Experiment for STORM-Central (PRESTORM)] squall line. The simulations are made using two different numerical advection schemes: 4th order and positive definite advection. Comparisons are made in terms of rainfall, cloud coverage, mass fluxes, apparent heating and moistening, mean hydrometeor profiles, CFADs (Contoured Frequency with Altitude Diagrams), microphysics, and latent heating retrieval. Overall, it was found that the different separation techniques produced results that qualitatively agreed. However, the quantitative differences were significant. Observational comparisons were unable to conclusively evaluate the performance of the techniques. Latent heating retrieval was shown to be sensitive to the use of separation technique mainly due to the stratiform region for methods that found very little stratiform rain. The midlatitude PRESTORM simulation was found to be nearly invariant with respect to advection type for most quantities while for TOGA COARE fourth order advection produced numerous shallow convective cores and positive definite advection fewer cells that were both broader and deeper penetrating above the freezing level.

13. A benchmark for numerical scheme validation of airborne particle exposure in street canyons.

PubMed

Marini, S; Buonanno, G; Stabile, L; Avino, P

2015-02-01

Measurements of particle concentrations and distributions in terms of number, surface area, and mass were performed simultaneously at eight sampling points within a symmetric street canyon of an Italian city. The aim was to obtain a useful benchmark for validation of wind tunnel experiments and numerical schemes: to this purpose, the influence of wind directions and speeds was considered. Particle number concentrations (PNCs) were higher on the leeward side than the windward side of the street canyon due to the wind vortex effect. Different vertical PNC profiles were observed between the two canyon sides depending on the wind direction and speed at roof level. A decrease in particle concentrations was observed with increasing rooftop wind speed, except for the coarse fraction indicating a possible particle resuspension due to the traffic and wind motion. This study confirms that particle concentration fields in urban street canyons are strongly influenced by traffic emissions and meteorological parameters, especially wind direction and speed. PMID:25167823

14. Numerical characterization of DNA sequences in a 2-D graphical representation scheme of low degeneracy

Guo, Xiaofeng; Nandy, Ashesh

2003-02-01

Some 2-D and 3-D graphical representations of DNA sequences have been given by Gate, Nandy, Leong, Randic, and Guo et al. Based on 2-D graphical representation of DNA sequences, Raychaudhury and Nandy introduced the first-order moments of the x and y coordinates and the radius of the plot of a DNA sequence for indexing scheme and similarity measures of DNA sequences. In this Letter, based on Guo's novel 2-D graphical representation of DNA sequences of low degeneracy, we introduce the improved first-order moments of the x and y coordinates and the radius of DNA sequences, and the distance of two DNA sequences. The new descriptors of DNA sequences give a good numerical characterization of DNA sequences, which have lower degeneracy.

15. Parallel solution of high-order numerical schemes for solving incompressible flows

NASA Technical Reports Server (NTRS)

Milner, Edward J.; Lin, Avi; Liou, May-Fun; Blech, Richard A.

1993-01-01

A new parallel numerical scheme for solving incompressible steady-state flows is presented. The algorithm uses a finite-difference approach to solving the Navier-Stokes equations. The algorithms are scalable and expandable. They may be used with only two processors or with as many processors as are available. The code is general and expandable. Any size grid may be used. Four processors of the NASA LeRC Hypercluster were used to solve for steady-state flow in a driven square cavity. The Hypercluster was configured in a distributed-memory, hypercube-like architecture. By using a 50-by-50 finite-difference solution grid, an efficiency of 74 percent (a speedup of 2.96) was obtained.

16. A numerical scheme and some theoretical aspects for the cylindrically and spherically symmetric sine-Gordon equations

Nguyen, Lu Trong Khiem

2016-07-01

A finite difference formula based on the predictor-corrector technique is presented to integrate the cylindrically and spherically symmetric sine-Gordon equations numerically. Based on various numerical observations, one property of the waves of kink type is conjectured and used to explain their returning effect. Several numerical experiments are carried out and they are in excellent agreement with the existing results. In addition, the corresponding modulation solution for the two-dimensional ring-shaped kink is extended to that in three-dimension. Both numerical and theoretical aspects are utilized to verify the reliability of the proposed numerical scheme and thus the analytical modulation solutions.

17. An implicit numerical scheme for the simulation of internal viscous flows on unstructured grids

NASA Technical Reports Server (NTRS)

Jorgenson, Philip C. E.; Pletcher, Richard H.

1994-01-01

The Navier-Stokes equations are solved numerically for two-dimensional steady viscous laminar flows. The grids are generated based on the method of Delaunay triangulation. A finite-volume approach is used to discretize the conservation law form of the compressible flow equations written in terms of primitive variables. A preconditioning matrix is added to the equations so that low Mach number flows can be solved economically. The equations are time marched using either an implicit Gauss-Seidel iterative procedure or a solver based on a conjugate gradient like method. A four color scheme is employed to vectorize the block Gauss-Seidel relaxation procedure. This increases the memory requirements minimally and decreases the computer time spent solving the resulting system of equations substantially. A factor of 7.6 speed up in the matrix solver is typical for the viscous equations. Numerical results are obtained for inviscid flow over a bump in a channel at subsonic and transonic conditions for validation with structured solvers. Viscous results are computed for developing flow in a channel, a symmetric sudden expansion, periodic tandem cylinders in a cross-flow, and a four-port valve. Comparisons are made with available results obtained by other investigators.

18. Numerical investigation of complex flooding schemes for surfactant polymer based enhanced oil recovery

Dutta, Sourav; Daripa, Prabir

2015-11-01

Surfactant-polymer flooding is a widely used method of chemical enhanced oil recovery (EOR) in which an array of complex fluids containing suitable and varying amounts of surfactant or polymer or both mixed with water is injected into the reservoir. This is an example of multiphase, multicomponent and multiphysics porous media flow which is characterized by the spontaneous formation of complex viscous fingering patterns and is modeled by a system of strongly coupled nonlinear partial differential equations with appropriate initial and boundary conditions. Here we propose and discuss a modern, hybrid method based on a combination of a discontinuous, multiscale finite element formulation and the method of characteristics to accurately solve the system. Several types of flooding schemes and rheological properties of the injected fluids are used to numerically study the effectiveness of various injection policies in minimizing the viscous fingering and maximizing oil recovery. Numerical simulations are also performed to investigate the effect of various other physical and model parameters such as heterogeneity, relative permeability and residual saturation on the quantities of interest like cumulative oil recovery, sweep efficiency, fingering intensity to name a few. Supported by the grant NPRP 08-777-1-141 from the Qatar National Research Fund (a member of The Qatar Foundation).

19. First order comparison of numerical calculation and two different turtle input schemes to represent a SLC defocusing magnet

SciTech Connect

Jaeger, J.

1983-07-14

Correcting the dispersion function in the SLC north arc it turned out that backleg-windings (BLW) acting horizontally as well as BLW acting vertically have to be used. In the latter case the question arose what is the best representation of a defocusing magnet with excited BLW acting in the vertical plane for the computer code TURTLE. Two different schemes, the 14.-scheme and the 20.-scheme were studied and the TURTLE output for one ray through such a magnet compared with the numerical solution of the equation of motion; only terms of first order have been taken into account.

20. Sensitivity analysis of numerical weather prediction radiative schemes to forecast direct solar radiation over Australia

Mukkavilli, S. K.; Kay, M. J.; Taylor, R.; Prasad, A. A.; Troccoli, A.

2014-12-01

The Australian Solar Energy Forecasting System (ASEFS) project requires forecasting timeframes which range from nowcasting to long-term forecasts (minutes to two years). As concentrating solar power (CSP) plant operators are one of the key stakeholders in the national energy market, research and development enhancements for direct normal irradiance (DNI) forecasts is a major subtask. This project involves comparing different radiative scheme codes to improve day ahead DNI forecasts on the national supercomputing infrastructure running mesoscale simulations on NOAA's Weather Research & Forecast (WRF) model. ASEFS also requires aerosol data fusion for improving accurate representation of spatio-temporally variable atmospheric aerosols to reduce DNI bias error in clear sky conditions over southern Queensland & New South Wales where solar power is vulnerable to uncertainities from frequent aerosol radiative events such as bush fires and desert dust. Initial results from thirteen years of Bureau of Meteorology's (BOM) deseasonalised DNI and MODIS NASA-Terra aerosol optical depth (AOD) anomalies demonstrated strong negative correlations in north and southeast Australia along with strong variability in AOD (~0.03-0.05). Radiative transfer schemes, DNI and AOD anomaly correlations will be discussed for the population and transmission grid centric regions where current and planned CSP plants dispatch electricity to capture peak prices in the market. Aerosol and solar irradiance datasets include satellite and ground based assimilations from the national BOM, regional aerosol researchers and agencies. The presentation will provide an overview of this ASEFS project task on WRF and results to date. The overall goal of this ASEFS subtask is to develop a hybrid numerical weather prediction (NWP) and statistical/machine learning multi-model ensemble strategy that meets future operational requirements of CSP plant operators.

SciTech Connect

Lipnikov, Konstantine; Danilov, A; Vassilevski, Y; Agonzal, A

2010-01-01

2. Comparison of ice-phase microphysical parameterization schemes using numerical simulations of tropical convection

NASA Technical Reports Server (NTRS)

Mccumber, Michael; Tao, Wei-Kuo; Simpson, Joanne; Penc, Richard; Soong, Su-Tzai

1991-01-01

The performance of several ice parameterizations has been evaluated through a numerical cloud model. Ice effects using different schemes are contrasted with each other and with an ice-free control by incorporating them into the cloud model and by applying them to simulations of tropical squall systems. The latter are simulated in 2D so that a large domain can be used to incorporate a complete anvil. Nonsquall-type convective lines are simulated in 3D owing to their smaller horizontal scale. It is concluded that inclusion of ice microphysics in the cloud model enhanced the agreement of the simulated convection with some features of observed convection, including the proportion of surface rainfall in the anvil region and the intensity and structure of the radar brightband near the melting level in the anvil. In the experiments with bulk microphysics, three ice categories produced much better results than two ice categories, which in turn was better than no ice. For the tropical squall-type and nonsquall-type systems the optimal mix was ice, snow, and graupel.

3. Analysis and Dynamically Consistent Numerical Schemes for the SIS Model and Related Reaction Diffusion Equation

Lubuma, J. M.-S.; Mureithi, E.; Terefe, Y. A.

2011-11-01

The classical SIS epidemiological model is extended in two directions: (a) The number of adequate contacts per infective in unit time is assumed to be a function of the total population in such a way that this number grows less rapidly as the total population increases; (b) A diffusion term is added to the SIS model and this leads to a reaction diffusion equation, which governs the spatial spread of the disease. With the parameter R0 representing the basic reproduction number, it is shown that R0 = 1 is a forward bifurcation for the model (a), with the disease-free equilibrium being globally asymptotic stable when R0 is less than 1. In the case when R0 is greater than 1, traveling wave solutions are found for the model (b). Nonstandard finite difference (NSFD) schemes that replicate the dynamics of the continuous models are presented. In particular, for the model (a), a nonstandard version of the Runge-Kutta method having high order of convergence is investigated. Numerical experiments that support the theory are provided.

4. Gas Evolution Dynamics in Godunov-Type Schemes and Analysis of Numerical Shock Instability

NASA Technical Reports Server (NTRS)

Xu, Kun

1999-01-01

In this paper we are going to study the gas evolution dynamics of the exact and approximate Riemann solvers, e.g., the Flux Vector Splitting (FVS) and the Flux Difference Splitting (FDS) schemes. Since the FVS scheme and the Kinetic Flux Vector Splitting (KFVS) scheme have the same physical mechanism and similar flux function, based on the analysis of the discretized KFVS scheme the weakness and advantage of the FVS scheme are closely observed. The subtle dissipative mechanism of the Godunov method in the 2D case is also analyzed, and the physical reason for shock instability, i.e., carbuncle phenomena and odd-even decoupling, is presented.

5. Numerical simulation of the debris flow dynamics with an upwind scheme and specific friction treatment

Sánchez Burillo, Guillermo; Beguería, Santiago; Latorre, Borja; Burguete, Javier

2014-05-01

Debris flows, snow and rock avalanches, mud and earth flows are often modeled by means of a particular realization of the so called shallow water equations (SWE). Indeed, a number of simulation models have been already developed [1], [2], [3], [4], [5], [6], [7]. Debris flow equations differ from shallow water equations in two main aspects. These are (a) strong bed gradient and (b) rheology friction terms that differ from the traditional SWE. A systematic analysis of the numerical solution of the hyperbolic system of equations rising from the shallow water equations with different rheological laws has not been done. Despite great efforts have been done to deal with friction expressions common in hydraulics (such as Manning friction), landslide rheologies are characterized by more complicated expressions that may deal to unphysical solutions if not treated carefully. In this work, a software that solves the time evolution of sliding masses over complex bed configurations is presented. The set of non- linear equations is treated by means of a first order upwind explicit scheme, and the friction contribution to the dynamics is treated with a suited numerical scheme [8]. In addition, the software incorporates various rheological models to accommodate for different flow types, such as the Voellmy frictional model [9] for rock and debris avalanches, or the Herschley-Bulkley model for debris and mud flows. The aim of this contribution is to release this code as a free, open source tool for the simulation of mass movements, and to encourage the scientific community to make use of it. The code uses as input data the friction coefficients and two input files: the topography of the bed and the initial (pre-failure) position of the sliding mass. In addition, another file with the final (post-event) position of the sliding mass, if desired, can be introduced to be compared with the simulation obtained result. If the deposited mass is given, an error estimation is computed by

6. Numerical simulation of the debris flow dynamics with an upwind scheme and specific friction treatment

Sánchez Burillo, Guillermo; Beguería, Santiago; Latorre, Borja; Burguete, Javier

2014-05-01

Debris flows, snow and rock avalanches, mud and earth flows are often modeled by means of a particular realization of the so called shallow water equations (SWE). Indeed, a number of simulation models have been already developed [1], [2], [3], [4], [5], [6], [7]. Debris flow equations differ from shallow water equations in two main aspects. These are (a) strong bed gradient and (b) rheology friction terms that differ from the traditional SWE. A systematic analysis of the numerical solution of the hyperbolic system of equations rising from the shallow water equations with different rheological laws has not been done. Despite great efforts have been done to deal with friction expressions common in hydraulics (such as Manning friction), landslide rheologies are characterized by more complicated expressions that may deal to unphysical solutions if not treated carefully. In this work, a software that solves the time evolution of sliding masses over complex bed configurations is presented. The set of non- linear equations is treated by means of a first order upwind explicit scheme, and the friction contribution to the dynamics is treated with a suited numerical scheme [8]. In addition, the software incorporates various rheological models to accommodate for different flow types, such as the Voellmy frictional model [9] for rock and debris avalanches, or the Herschley-Bulkley model for debris and mud flows. The aim of this contribution is to release this code as a free, open source tool for the simulation of mass movements, and to encourage the scientific community to make use of it. The code uses as input data the friction coefficients and two input files: the topography of the bed and the initial (pre-failure) position of the sliding mass. In addition, another file with the final (post-event) position of the sliding mass, if desired, can be introduced to be compared with the simulation obtained result. If the deposited mass is given, an error estimation is computed by

7. An Improved Transformation and Optimized Sampling Scheme for the Numerical Evaluation of Singular and Near-Singular Potentials

NASA Technical Reports Server (NTRS)

Khayat, Michael A.; Wilton, Donald R.; Fink, Patrick W.

2007-01-01

Simple and efficient numerical procedures using singularity cancellation methods are presented for evaluating singular and near-singular potential integrals. Four different transformations are compared and the advantages of the Radial-angular transform are demonstrated. A method is then described for optimizing this integration scheme.

8. On the numerical simulation of particle dynamics in the radiation belt: 1. Implicit and semi-implicit schemes

Camporeale, E.; Delzanno, G. L.; Zaharia, S.; Koller, J.

2013-06-01

The particle dynamics in the Earth's radiation belt is generally modeled by means of a two-dimensional diffusion equation for the particle distribution function in energy and pitch angle. The goal of this paper is to survey and compare different numerical schemes for the solution of the diffusion equation, and to outline the optimal strategy from a numerical point of view. We focus on the general (and more computationally challenging) case where the mixed terms in the diffusion tensor are retained. In Part 1, we compare fully implicit and semi-implicit schemes. For the former, we have analyzed a direct solver based on a LU decomposition routine for sparse matrices, and an iterative incomplete LU preconditioned Generalized Minimal REsidual solver. For the semi-implicit scheme, we have studied an alternating direction implicit scheme. We present a convergence study for a realistic case that shows that the time step and grid size are strongly constrained by the desired accuracy of the solution. We show that the fully implicit scheme is to be preferred in most cases as the more computationally efficient.

9. A numerical scheme for modelling reacting flow with detailed chemistry and transport.

SciTech Connect

Knio, Omar M.; Najm, Habib N.; Paul, Phillip H.

2003-09-01

An efficient projection scheme is developed for the simulation of reacting flow with detailed kinetics and transport. The scheme is based on a zero-Mach-number formulation of the compressible conservation equations for an ideal gas mixture. It is a modified version of the stiff operator-split scheme developed by Knio, Najm & Wyckoff (1999, J. Comput. Phys. 154, 428). Similar to its predecessor, the new scheme relies on Strang splitting of the discrete evolution equations, where diffusion is integrated in two half steps that are symmetrically distributed around a single stiff step for the reaction source terms. The diffusive half-step is integrated using an explicit single-step, multistage, Runge-Kutta-Chebyshev (RKC) method, which replaces the explicit, multi-step, fractional sub-step approach used in the previous formulation. This modification maintains the overall second-order convergence properties of the scheme and enhances the efficiency of the computations by taking advantage of the extended real-stability region of the RKC scheme. Two additional efficiency-enhancements are also explored, based on an extrapolation procedure for the transport coefficients and on the use of approximate Jacobian data evaluated on a coarse mesh. By including these enhancement schemes, performance tests using 2D computations with a detailed C{sub 1}C{sub 2} methane-air mechanism and a detailed mixture-averaged transport model indicate that speedup factors of about 15 are achieved over the previous split-stiff scheme.

10. Nonstandard Finite Difference Schemes: Relations Between Time and Space Step-Sizes in Numerical Schemes for PDE's That Follow from Positivity Condition

NASA Technical Reports Server (NTRS)

Mickens, Ronald E.

1996-01-01

A large class of physical phenomena can be modeled by evolution and wave type Partial Differential Equations (PDE). Few of these equations have known explicit exact solutions. Finite-difference techniques are a popular method for constructing discrete representations of these equations for the purpose of numerical integration. However, the solutions to the difference equations often contain so called numerical instabilities; these are solutions to the difference equations that do not correspond to any solution of the PDE's. For explicit schemes, the elimination of this behavior requires functional relations to exist between the time and space steps-sizes. We show that such functional relations can be obtained for certain PDE's by use of a positivity condition. The PDE's studied are the Burgers, Fisher, and linearized Euler equations.

11. An Advanced Manipulator For Poisson Series With Numerical Coefficients

Biscani, Francesco; Casotto, S.

2006-06-01

The availability of an efficient and featureful manipulator for Poisson deries with numerical coefficients is a standard need for celestial mechanicians and has arisen during our work on the analytical development of the Tide-Generating-Potential (TGP). In the harmonic expansion of the TGP the Poisson series appearing in the theories of motion of the celestial bodies are subjected to a wide set of mathematical operations, ranging from simple additions and multiplications to more sophisticated operations on Legendre polynomials and spherical harmonics with Poisson series as arguments. To perform these operations we have developed an algebraic manipulator, called Piranha, structured as an object-oriented multi-platform C++ library. Piranha handles series with real and complex coefficients, and operates with an arbitrary degree of precision. It supports advanced features such as trigonometric operations and the generation of special functions from Poisson series. Piranha is provided with a proof-of-concept, multi-platform GUI, which serves as a testbed and benchmark for the library. We describe Piranha's architecture and characteristics, what it accomplishes currently and how it will be extended in the future (e.g., to handle series with symbolic coefficients in a consistent fashion with its current design).

12. Dynamic design, numerical solution and effective verification of acceleration-level obstacle-avoidance scheme for robot manipulators

Xiao, Lin; Zhang, Yunong

2016-03-01

For avoiding obstacles and joint physical constraints of robot manipulators, this paper proposes and investigates a novel obstacle avoidance scheme (termed the acceleration-level obstacle-avoidance scheme). The scheme is based on a new obstacle-avoidance criterion that is designed by using the gradient neural network approach for the first time. In addition, joint physical constraints such as joint-angle limits, joint-velocity limits and joint-acceleration limits are incorporated into such a scheme, which is further reformulated as a quadratic programming (QP). Two important 'bridge' theorems are established so that such a QP can be converted equivalently to a linear variational inequality and then equivalently to a piecewise-linear projection equation (PLPE). A numerical algorithm based on a PLPE is thus developed and applied for an online solution of the resultant QP. Four path-tracking tasks based on the PA10 robot in the presence of point and window-shaped obstacles demonstrate and verify the effectiveness and accuracy of the acceleration-level obstacle-avoidance scheme. Besides, the comparisons between the non-obstacle-avoidance and obstacle-avoidance results further validate the superiority of the proposed scheme.

13. A hybrid numerical technique for predicting the aerodynamic and acoustic fields of advanced turboprops

NASA Technical Reports Server (NTRS)

Homicz, G. F.; Moselle, J. R.

1985-01-01

A hybrid numerical procedure is presented for the prediction of the aerodynamic and acoustic performance of advanced turboprops. A hybrid scheme is proposed which in principle leads to a consistent simultaneous prediction of both fields. In the inner flow a finite difference method, the Approximate-Factorization Alternating-Direction-Implicit (ADI) scheme, is used to solve the nonlinear Euler equations. In the outer flow the linearized acoustic equations are solved via a Boundary-Integral Equation (BIE) method. The two solutions are iteratively matched across a fictitious interface in the flow so as to maintain continuity. At convergence the resulting aerodynamic load prediction will automatically satisfy the appropriate free-field boundary conditions at the edge of the finite difference grid, while the acoustic predictions will reflect the back-reaction of the radiated field on the magnitude of the loading source terms, as well as refractive effects in the inner flow. The equations and logic needed to match the two solutions are developed and the computer program implementing the procedure is described. Unfortunately, no converged solutions were obtained, due to unexpectedly large running times. The reasons for this are discussed and several means to alleviate the situation are suggested.

14. Numerical solution of 3D Navier-Stokes equations with upwind implicit schemes

NASA Technical Reports Server (NTRS)

Marx, Yves P.

1990-01-01

An upwind MUSCL type implicit scheme for the three-dimensional Navier-Stokes equations is presented. Comparison between different approximate Riemann solvers (Roe and Osher) are performed and the influence of the reconstructions schemes on the accuracy of the solution as well as on the convergence of the method is studied. A new limiter is introduced in order to remove the problems usually associated with non-linear upwind schemes. The implementation of a diagonal upwind implicit operator for the three-dimensional Navier-Stokes equations is also discussed. Finally the turbulence modeling is assessed. Good prediction of separated flows are demonstrated if a non-equilibrium turbulence model is used.

15. An implicit factored scheme for the compressible Navier-Stokes equations. II - The numerical ODE connection

NASA Technical Reports Server (NTRS)

Beam, R. M.; Warming, R. F.

1979-01-01

An attempt is made to establish a connection between linear multistep methods for applications to ordinary differential equations and their extension (by approximate factorization) to alternating direction implicit methods for partial differential equations. An earlier implicit factored scheme for the compressible Navier-Stokes equations is generalized by innovations that (1) increase the class of temporal difference schemes to include all linear multistep methods, (2) optimize the class of unconditionally stable factored schemes by a new choice of unknown variable, and (3) improve the computational efficiency by the introduction of quasi-one-leg methods.

16. Recent advances in the development of implicit schemes for the equations of fluid dynamics

NASA Technical Reports Server (NTRS)

Warming, R. F.; Beam, R. M.

1981-01-01

Innovations and extensions of implicit schemes for equations of fluid dynamics are presented. The notation and theory for linear multistep methods are reviewed, and extensions of work by Beam and Warming (1979) include the implementation of one-leg methods, ADI methods for equations with mixed derivatives, flux vector splitting, the P-dimensional wave equation, and boundary conditions. Numerical experiments indicate that implicit treatment of the boundary conditions is necessary for unconditional stability, and the improvement and implementation of the boundary condition theory should improve the implicit algorithms for gas dynamic equations.

17. Enhanced numerical inviscid and viscous fluxes for cell centered finite volume schemes

Eberle, Albrecht

1991-07-01

The most attractive features of cell centered finite volume schemes seem to be the easy introduction of the solid body boundary condition and the implementation of characteristic based methods for evaluating the convective fluxes at the cell faces of the finite volumes. For the viscous parts of the fluxes, however, cell centered finite volume schemes are not as well suited as cell vertex based discretizations since in a general grid, cell centered schemes usually are not linear flow preserving concerning the viscous terms. That means that the viscous stress tensor and the heat flux vector may spuriously vary in a flow field with linear velocity and/or temperature distribution. Several enhancements of the flux formulations for cell centered finite volume schemes are described.

18. Numerical Simulation of Chennai Heavy Rainfall Using MM5 Mesoscale Model with Different Cumulus Parameterization Schemes

Litta, A. J.; Chakrapani, B.; Mohankumar, K.

2007-07-01

Heavy rainfall events become significant in human affairs when they are combined with hydrological elements. The problem of forecasting heavy precipitation is especially difficult since it involves making a quantitative precipitation forecast, a problem well recognized as challenging. Chennai (13.04°N and 80.17°E) faced incessant and heavy rain about 27 cm in 24 hours up to 8.30 a.m on 27th October 2005 completely threw life out of gear. This torrential rain caused by deep depression which lay 150km east of Chennai city in Bay of Bengal intensified and moved west north-west direction and crossed north Tamil Nadu and south Andhra Pradesh coast on 28th morning. In the present study, we investigate the predictability of the MM5 mesoscale model using different cumulus parameterization schemes for the heavy rainfall event over Chennai. MM5 Version 3.7 (PSU/NCAR) is run with two-way triply nested grids using Lambert Conformal Coordinates (LCC) with a nest ratio of 3:1 and 23 vertical layers. Grid sizes of 45, 15 and 5 km are used for domains 1, 2 and 3 respectively. The cumulus parameterization schemes used in this study are Anthes-Kuo scheme (AK), the Betts-Miller scheme (BM), the Grell scheme (GR) and the Kain-Fritsch scheme (KF). The present study shows that the prediction of heavy rainfall is sensitive to cumulus parameterization schemes. In the time series of rainfall, Grell scheme is in good agreement with observation. The ideal combination of the nesting domains, horizontal resolution and cloud parameterization is able to simulate the heavy rainfall event both qualitatively and quantitatively.

19. The numerical prediction of planar viscoelastic contraction flows using the pom pom model and higher-order finite volume schemes

Aguayo, J. P.; Phillips, P. M.; Phillips, T. N.; Tamaddon-Jahromi, H. R.; Snigerev, B. A.; Webster, M. F.

2007-01-01

This study investigates the numerical solution of viscoelastic flows using two contrasting high-order finite volume schemes. We extend our earlier work for Poiseuille flow in a planar channel and the single equation form of the extended pom-pom (SXPP) model [M. Aboubacar, J.P. Aguayo, P.M. Phillips, T.N. Phillips, H.R. Tamaddon-Jahromi, B.A. Snigerev, M.F. Webster, Modelling pom-pom type models with high-order finite volume schemes, J. Non-Newtonian Fluid Mech. 126 (2005) 207-220], to determine steady-state solutions for planar 4:1 sharp contraction flows. The numerical techniques employed are time-stepping algorithms: one of hybrid finite element/volume type, the other of pure finite volume form. The pure finite volume scheme is a staggered-grid cell-centred scheme based on area-weighting and a semi-Lagrangian formulation. This may be implemented on structured or unstructured rectangular grids, utilising backtracking along the solution characteristics in time. For the hybrid scheme, we solve the momentum-continuity equations by a fractional-staged Taylor-Galerkin pressure-correction procedure and invoke a cell-vertex finite volume scheme for the constitutive law. A comparison of the two finite volume approaches is presented, concentrating upon the new features posed by the pom-pom class of models in this context of non-smooth flows. Here, the dominant feature of larger shear and extension in the entry zone influences both stress and stretch, so that larger stretch develops around the re-entrant corner zone as Weissenberg number increases, whilst correspondingly stress levels decline.

20. Accelerating development of advanced inverters : evaluation of anti-islanding schemes with grid support functions and preliminary laboratory demonstration.

SciTech Connect

Neely, Jason C.; Gonzalez, Sigifredo; Ropp, Michael; Schutz, Dustin

2013-11-01

The high penetration of utility interconnected photovoltaic (PV) systems is causing heightened concern over the effect that variable renewable generation will have on the electrical power system (EPS). These concerns have initiated the need to amend the utility interconnection standard to allow advanced inverter control functionalities that provide: (1) reactive power control for voltage support, (2) real power control for frequency support and (3) better tolerance of grid disturbances. These capabilities are aimed at minimizing the negative impact distributed PV systems may have on EPS voltage and frequency. Unfortunately, these advanced control functions may interfere with island detection schemes, and further development of advanced inverter functions requires a study of the effect of advanced functions on the efficacy of antiislanding schemes employed in industry. This report summarizes the analytical, simulation and experimental work to study interactions between advanced inverter functions and anti-islanding schemes being employed in distributed PV systems.

1. Numerical simulation by TVD schemes of complex shock reflections from airfoils at high angle of attack. [Total Variation Diminishing

NASA Technical Reports Server (NTRS)

Moon, Young J.; Yee, H. C.

1987-01-01

The shock-capturing capability of total variation diminishing (TVD) schemes is demonstrated for a more realistic complex shock-diffraction problem for which the experimental data are available. Second-order explicit upwind and symmetric TVD schemes are used to solve the time-dependent Euler equations of gas dynamics for the interaction of a blast wave with an airfoil at high angle-of-attack. The test cases considered are a time-dependent moving curved-shock wave and a contant moving planar-shock wave impinging at an angle-of-attack 30 deg on a NACA 0018 airfoil. Good agreement is obtained between isopycnic contours computed by the TVD schemes and those from experimental interferograms. No drastic difference in flow-field structure is found between the curved- and planar-shock wave cases, except for a difference in density level near the lower surface of the airfoil. Computation for cases with higher shock Mach numbers is also possible. Numerical experiments show that the symmetric TVD scheme is less sensitive to the boundary conditions treatment than the upwind scheme.

2. Conservative numerical schemes for unsteady one-dimensional two phase flow

Garcia Cascales, Jose Ramon

The thesis is devoted to the modelization of non steady two phase mixtures of liquid and vapour. It has been motivated by the great amount of industrial applications in which we find these phenomena. Transient two phase flow is a very important issue in nuclear, chemical and industrial applications. In the case of the nuclear industry due to the importance of preventing loss of coolant accidents (LOCA) and guaranteeing a good performance of the coolant system in power plants. We justify the present development by means of the introduction of the most important codes developed during the last two decades and their associated mesh techniques. It is basically focused on the extension of some conservative and explicit schemes to obtain approximate solutions of the system of equations in one dimensional one pressure two phase flow. They have been centred and upwind schemes to solve multiphase flow problems, most of them based on the exact or approximate solution of Riemann problems using Godunov's like methods such as Approximate Riemann solvers or Flux Splitting methods. We have studied mainly TVD schemes, Adapted TVD schemes (ATVD) and the AUSM family of schemes. Firstly we introduce the 1D two phase flow system of equations with which we will work. We consider the systems of equations more used depending on the model. Thus we introduce the homogeneous model, the isentropic model and the separated model will be treated in some detail. The evaluation of the eigenstructure of the homogeneous and the separated two phase flow is studied. Different methods to determine the eigenvalues are presented. A general method to determine the eigenvectors is studied as well. We extend different conservative schemes to two phase flow whose good behaviour in single phase has been well proved. They are basically TVD schemes, the Adapted TVD schemes developed by Gascon and Corberan and the AUSM family of schemes, firstly introduced by Steffen and Liou. Most of the extensions developed

3. A numerical resolution study of high order essentially non-oscillatory schemes applied to incompressible flow

NASA Technical Reports Server (NTRS)

Weinan, E.; Shu, Chi-Wang

1992-01-01

High order essentially non-oscillatory (ENO) schemes, originally designed for compressible flow and in general for hyperbolic conservation laws, are applied to incompressible Euler and Navier-Stokes equations with periodic boundary conditions. The projection to divergence-free velocity fields is achieved by fourth order central differences through Fast Fourier Transforms (FFT) and a mild high-order filtering. The objective of this work is to assess the resolution of ENO schemes for large scale features of the flow when a coarse grid is used and small scale features of the flow, such as shears and roll-ups, are not fully resolved. It is found that high-order ENO schemes remain stable under such situations and quantities related to large-scale features, such as the total circulation around the roll-up region, are adequately resolved.

4. Numerical Speed of Sound and its Application to Schemes for all Speeds

NASA Technical Reports Server (NTRS)

Liou, Meng-Sing; Edwards, Jack R.

1999-01-01

The concept of "numerical speed of sound" is proposed in the construction of numerical flux. It is shown that this variable is responsible for the accurate resolution of' discontinuities, such as contacts and shocks. Moreover, this concept can he readily extended to deal with low speed and multiphase flows. As a results, the numerical dissipation for low speed flows is scaled with the local fluid speed, rather than the sound speed. Hence, the accuracy is enhanced the correct solution recovered, and the convergence rate improved. We also emphasize the role of mass flux and analyze the behavior of this flux. Study of mass flux is important because the numerical diffusivity introduced in it can be identified. In addition, it is the term common to all conservation equations. We show calculated results for a wide variety of flows to validate the effectiveness of using the numerical speed of sound concept in constructing the numerical flux. We especially aim at achieving these two goals: (1) improving accuracy and (2) gaining convergence rates for all speed ranges. We find that while the performance at high speed range is maintained, the flux now has the capability of performing well even with the low: speed flows. Thanks to the new numerical speed of sound, the convergence is even enhanced for the flows outside of the low speed range. To realize the usefulness of the proposed method in engineering problems, we have also performed calculations for complex 3D turbulent flows and the results are in excellent agreement with data.

5. A novel type II relay-assisted retransmission scheme for uplink of LTE-advanced system

Li, Anxin; Nagata, Satoshi; Harada, Atsushi; Suda, Hirohito

2013-12-01

Relay, which enables coverage extension and throughput enhancement, is a very promising technique for future wireless communication systems. Among different types of relay, type II relay is one kind of inband relays and is hotly discussed in LTE-Advanced system for throughput enhancement. In order to support type II relay, many challenges must be overcome. In this article, we focus on relay-assisted uplink data retransmission and propose a novel joint design of reference signal and data precoding for type II relay. The proposed method not only solves the problem of channel estimation mismatch for control information, but also achieves cooperative diversity gain for data transmission. Simulation results show the effectiveness of the proposed method over existing schemes.

6. Recent advances in two-phase flow numerics

SciTech Connect

Mahaffy, J.H.; Macian, R.

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.

NASA Technical Reports Server (NTRS)

Usab, William J., Jr.; Jiang, Yi-Tsann

1991-01-01

The objective of the present research is to develop a general solution adaptive scheme for the accurate prediction of inviscid quasi-three-dimensional flow in advanced compressor and turbine designs. The adaptive solution scheme combines an explicit finite-volume time-marching scheme for unstructured triangular meshes and an advancing front triangular mesh scheme with a remeshing procedure for adapting the mesh as the solution evolves. The unstructured flow solver has been tested on a series of two-dimensional airfoil configurations including a three-element analytic test case presented here. Mesh adapted quasi-three-dimensional Euler solutions are presented for three spanwise stations of the NASA rotor 67 transonic fan. Computed solutions are compared with available experimental data.

8. On one numerical scheme of the solution of a three-dimensional problem of diffraction of an electromagnetic wave on thin ideally conductive screens

SciTech Connect

Ryzhakov, G. V.

2014-11-12

In the paper, the problem of diffraction on thin ideally conductive screens is reduced to vector hypersingular integral equation with integral treated in the sense of finite Hadamard value. An numerical scheme to solve the equation is introduced. The scheme is based on piecewise approximation of unknown function. The advantage of the scheme is that integral of singular part is reduced to contour integral which can be analytically calculated so numerical calculation are significantly accelerated. Several examples of resulting numerical experiments are given in comparison with known theoretical and experimental data.

9. Comparison of Numerical Schemes for a Realistic Computational Aeroacoustics Benchmark Problem

NASA Technical Reports Server (NTRS)

Hixon, R.; Wu, J.; Nallasamy, M.; Sawyer, S.; Dyson, R.

2004-01-01

In this work, a nonlinear structured-multiblock CAA solver, the NASA GRC BASS code, will be tested on a realistic CAA benchmark problem. The purpose of this test is to ascertain what effect the high-accuracy solution methods used in CAA have on a realistic test problem, where both the mean flow and the unsteady waves are simultaneously computed on a fully curvilinear grid from a commercial grid generator. The proposed test will compare the solutions obtained using several finite-difference methods on identical grids to determine whether high-accuracy schemes have advantages for this benchmark problem.

10. Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes

NASA Technical Reports Server (NTRS)

Jameson, A.; Schmidt, Wolfgang; Turkel, Eli

1981-01-01

A new combination of a finite volume discretization in conjunction with carefully designed dissipative terms of third order, and a Runge Kutta time stepping scheme, is shown to yield an effective method for solving the Euler equations in arbitrary geometric domains. The method has been used to determine the steady transonic flow past an airfoil using an O mesh. Convergence to a steady state is accelerated by the use of a variable time step determined by the local Courant member, and the introduction of a forcing term proportional to the difference between the local total enthalpy and its free stream value.

11. Application of vortex identification schemes to direct numerical simulation data of a transitional boundary layer

Pierce, Brian; Moin, Parviz; Sayadi, Taraneh

2013-01-01

We have demonstrated how various vortex identification and visualization criteria perform using direct numerical simulation data from a transitional and turbulent boundary layer by Sayadi, Hamman, and Moin ["Direct numerical simulation of complete transition to turbulence via h-type and k-type secondary instabilities," Technical Report, Stanford University, CTR Annual Research Briefs, 2011]. The presence of well-known Λ vortices in the transitional region provides a well defined and yet realistic benchmark for evaluation of various criteria. We investigate the impact of changing the threshold used for iso-surface plotting.

12. System-level performance of LTE-Advanced with joint transmission and dynamic point selection schemes

Määttänen, Helka-Liina; Hämäläinen, Kari; Venäläinen, Juha; Schober, Karol; Enescu, Mihai; Valkama, Mikko

2012-12-01

In this article, we present a practical coordinated multipoint (CoMP) system for LTE-Advanced. In this CoMP system, cooperation is enabled for cell-edge users via dynamic switching between the normal single-cell operation and CoMP. We first formulate a general CoMP system model of several CoMP schemes. We then investigate a practical finite-rate feedback design that simultaneously supports interference coordination, joint transmission (JT), and dynamic point selection (DPS) with a varying number of cooperating transmission points while operating a single-cell transmission as a fallback mode. We provide both link-level and system-level results for the evaluation of different feedback options for general CoMP operation. The results show that there are substantial performance gains in cell-edge throughputs for both JT and DPS CoMP over the baseline Release 10 LTE-Advanced with practical feedback options. We also show that CoMP can enable improved mobility management in real networks.

13. Generation of a composite grid for turbine flows and consideration of a numerical scheme

NASA Technical Reports Server (NTRS)

Choo, Y.; Yoon, S.; Reno, C.

1986-01-01

A composite grid was generated for flows in turbines. It consisted of the C-grid (or O-grid) in the immediate vicinity of the blade and the H-grid in the middle of the blade passage between the C-grids and in the upstream region. This new composite grid provides better smoothness, resolution, and orthogonality than any single grid for a typical turbine blade with a large camber and rounded leading and trailing edges. The C-H (or O-H) composite grid has an unusual grid point that is connected to more than four neighboring nodes in two dimensions (more than six neighboring nodes in three dimensions). A finite-volume lower-upper (LU) implicit scheme to be used on this grid poses no problem and requires no special treatment because each interior cell of this composite grid has only four neighboring cells in two dimensions (six cells in three dimensions). The LU implicit scheme was demonstrated to be efficient and robust for external flows in a broad flow regime and can be easily applied to internal flows and extended from two to three dimensions.

14. EVALUATION OF NUMERICAL SCHEMES FOR SOLVING A CONSERVATION OF SPECIES EQUATION WITH CHEMICAL TERMS

EPA Science Inventory

Numerical methods are investigated for solving a system of continuity equations that contain linear and nonlinear chemistry as source and sink terms. It is shown that implicit, finite-difference approximations, when applied to the chemical kinetic terms, yield accurate results wh...

15. Numerical study of unsteady shockwave reflections using an upwind TVD scheme

NASA Technical Reports Server (NTRS)

Hsu, Andrew T.; Liou, Meng-Sing

1990-01-01

An unsteady TVD Navier-Stokes solver was developed and applied to the problem of shock reflection on a circular cylinder. The obtained numerical results were compared with the Schlieren photos from an experimental study. These results show that the present computer code has the ability of capturing moving shocks.

16. Numerical simulation of Stokes flow around particles via a hybrid Finite Difference-Boundary Integral scheme

Bhattacharya, Amitabh

2013-11-01

An efficient algorithm for simulating Stokes flow around particles is presented here, in which a second order Finite Difference method (FDM) is coupled to a Boundary Integral method (BIM). This method utilizes the strong points of FDM (i.e. localized stencil) and BIM (i.e. accurate representation of particle surface). Specifically, in each iteration, the flow field away from the particles is solved on a Cartesian FDM grid, while the traction on the particle surface (given the the velocity of the particle) is solved using BIM. The two schemes are coupled by matching the solution in an intermediate region between the particle and surrounding fluid. We validate this method by solving for flow around an array of cylinders, and find good agreement with Hasimoto's (J. Fluid Mech. 1959) analytical results.

17. The Role of Numerical Simulation in Advancing Plasma Propulsion

Turchi, P. J.; Mikellides, P. G.; Mikellides, I. G.

1999-11-01

Plasma thrusters often involve a complex set of interactions among several distinct physical processes. While each process can yield to separate mathematical representation, their combination generally requires numerical simulation. We have extended and used the MACH2 code successfully to simulate both self-field and applied-field magnetoplasmadynamic thrusters and, more recently, ablation-fed pulsed plasma microthrusters. MACH2 provides a framework in which to compute 2-1/2 dimensional, unsteady, MHD flows in two-temperature LTE. It couples to several options for electrical circuitry and allows access to both analytic formulas and tabular values for material properties and transport coefficients, including phenomenological models for anomalous transport. Even with all these capabilities, however, successful modeling demands comparison with experiment and with analytic solutions in idealized limits, and careful combination of MACH2 results with separate physical reasoning. Although well understood elsewhere in plasma physics, the strengths and limitations of numerical simulation for plasma propulsion needs further discussion.

18. Numerical and experimental investigation into passive hydrogen recovery scheme using vacuum ejector

Hwang, Jenn-Jiang; Cho, Ching-Chang; Wu, Wei; Chiu, Ching-Huang; Chiu, Kuo-Ching; Lin, Chih-Hong

2015-02-01

The current work presents a numerical and experimental investigation into a passive ejector for recovering the anode off-gas in a proton exchange membrane fuel cell (PEMFC) system. The proposed ejector is consisted of a convergent-divergent channel and a suction channel, and it is connected with the anode outlet of PEMFC system for recovery the anode off-gas into the main gas supply. Numerical simulations based on a three-dimensional compressible steady-state k-ɛ turbulent model are performed to examine the effects of the inlet mass flow rate and nozzle throat diameter on the pressure, Mach number, temperature, suction channel mass flow rate, outlet channel mass flow rate, and suction channel entrainment ratio, respectively. The numerical results are confirmed by means of an experimental investigation. It is shown that supersonic flow conditions are induced in the ejector; resulting in the induction of a vacuum pressure in the suction channel and the subsequent recovery of the anode off-gas at the outlet of the main channel. In addition, it is shown that the mass flow rate in the suction channel increases with an increasing mass flow rate at the primary channel inlet. Finally, the results show that a higher entrainment ratio is obtained as the throat diameter of the nozzle in the ejector is reduced. Overall, the results presented in this study provide a useful source of reference for developing the ejector devices applied to fuel cell systems while simultaneously avoiding extra energy consumption.

19. A nonlocal finite difference scheme for simulation of wave propagation in 2D models with reduced numerical dispersion

Martowicz, A.; Ruzzene, M.; Staszewski, W. J.; Rimoli, J. J.; Uhl, T.

2014-03-01

The work deals with the reduction of numerical dispersion in simulations of wave propagation in solids. The phenomenon of numerical dispersion naturally results from time and spatial discretization present in a numerical model of mechanical continuum. Although discretization itself makes possible to model wave propagation in structures with complicated geometries and made of different materials, it inevitably causes simulation errors when improper time and length scales are chosen for the simulations domains. Therefore, by definition, any characteristic parameter for spatial and time resolution must create limitations on maximal wavenumber and frequency for a numerical model. It should be however noted that expected increase of the model quality and its functionality in terms of affordable wavenumbers, frequencies and speeds should not be achieved merely by denser mesh and reduced time integration step. The computational cost would be simply unacceptable. The authors present a nonlocal finite difference scheme with the coefficients calculated applying a Fourier series, which allows for considerable reduction of numerical dispersion. There are presented the results of analyses for 2D models, with isotropic and anisotropic materials, fulfilling the planar stress state. Reduced numerical dispersion is shown in the dispersion surfaces for longitudinal and shear waves propagating for different directions with respect to the mesh orientation and without dramatic increase of required number of nonlocal interactions. A case with the propagation of longitudinal wave in composite material is studied with given referential solution of the initial value problem for verification of the time-domain outcomes. The work gives a perspective of modeling of any type of real material dispersion according to measurements and with assumed accuracy.

20. Numerical scheme for a spatially inhomogeneous matrix-valued quantum Boltzmann equation

Lu, Jianfeng; Mendl, Christian B.

2015-06-01

We develop an efficient algorithm for a spatially inhomogeneous matrix-valued quantum Boltzmann equation derived from the Hubbard model. The distribution functions are 2 × 2 matrix-valued to accommodate the spin degree of freedom, and the scalar quantum Boltzmann equation is recovered as a special case when all matrices are proportional to the identity. We use Fourier discretization and fast Fourier transform to efficiently evaluate the collision kernel with spectral accuracy, and numerically investigate periodic, Dirichlet and Maxwell boundary conditions. Model simulations quantify the convergence to local and global thermal equilibrium.

1. Static structure of general relativistic partially degenerate, self-gravitational, lepton configurations I. Numerical scheme

SciTech Connect

Chau, W.Y.; Lake, K.; Stone, J.

1984-06-15

The Oppenheimer-Volkoff equations for isothermal, partially degenerate neutral lepton configurations have been solved numerically for arbitrary temperature and degree of degeneracy. For the specific case where the leptons are massive neutrions (approx.30 eV) with properties quite tightly constrained by cosmological considerations, we find that it possible to have a neutrino halo surrounding a normal galaxy with the right values of mass and radius required of the ''invisible halo'' in the missing mass problem. The density distribution, however, yields a rotational curve which, with the present chosen values of parameters, does not fit the observations for spiral galaxies.

2. Improved computational schemes for the numerical modeling of hydrothermal resources in Wyoming

SciTech Connect

Heasler, H.P.; George, J.H.; Allen, M.B.

1990-05-01

A new method, the Conjugate Gradient Squared (CGS) solution technique, is shown to be extremely effective when applied to the finite-difference solution of conductive and convective heat transfer in geologic systems. The CGS method is compared to the Successive Over/Under Relaxation schemes, a version of the Gaussian elimination method, and the Generalized Minimum Residual (GMRES) approach. The CGS procedure converges at least ten times faster than the nearest competitor. The model is applied to the Thermopolis hydrothermal system, located in northwestern Wyoming. Modeled results are compared with measured temperature-depth profiles and results from other studies. The temperature decrease from 72{degree}C to 54{degrees}C along the crest of the Thermopolis anticline is shown to result from cooling of the geothermal fluid as it moves to the southeast. Modeled results show correct general trends, however, a time-varying three-dimensional model will be needed to fully explain the effects of mixing within the aquifers along the crest of the anticline and thermal affects of surface surface topography. 29 refs., 18 figs., 2 tabs.

3. A discontinuous wave-in-cell numerical scheme for hyperbolic conservation laws

Thompson, Richard J.; Moeller, Trevor

2015-10-01

A new Riemann solver scheme for hyperbolic systems is introduced. The method consists of a discretization of the initial data into an approximate representation by discrete, discontinuous waves. Instead of calculating an intercell flux based on these waves, the discontinuous waves are propagated directly. Since the sum total of all discontinuous waves represents an extension of the linear Riemann problem, the solution is determined straightforwardly. For nonlinear systems, each timestep is considered a separate linear Riemann problem, and the projected waves are weighted to a background grid. This method is strikingly similar to the particle-in-cell approach, except that discontinuous waves are pushed around instead of macroparticles. The method is applied to Maxwell's equations and the equations of inviscid gasdynamics. For linear systems, exact solutions can be achieved without dissipation, and exact transmissive boundary condition treatments are trivial to implement. For inviscid gasdynamics, the nonlinear method tended to resolve discontinuities more sharply than the Roe, HLL or HLLC methods while requiring only between 30% and 50% of the execution time under identical conditions. The approach is also extremely robust, as it works for any Courant number. In the limit that the Courant number becomes infinite, the nonlinear solution approaches the linearized solution.

4. Basic and Advanced Numerical Performances Relate to Mathematical Expertise but Are Fully Mediated by Visuospatial Skills

PubMed Central

2016-01-01

Recent studies have highlighted the potential role of basic numerical processing in the acquisition of numerical and mathematical competences. However, it is debated whether high-level numerical skills and mathematics depends specifically on basic numerical representations. In this study mathematicians and nonmathematicians performed a basic number line task, which required mapping positive and negative numbers on a physical horizontal line, and has been shown to correlate with more advanced numerical abilities and mathematical achievement. We found that mathematicians were more accurate compared with nonmathematicians when mapping positive, but not negative numbers, which are considered numerical primitives and cultural artifacts, respectively. Moreover, performance on positive number mapping could predict whether one is a mathematician or not, and was mediated by more advanced mathematical skills. This finding might suggest a link between basic and advanced mathematical skills. However, when we included visuospatial skills, as measured by block design subtest, the mediation analysis revealed that the relation between the performance in the number line task and the group membership was explained by non-numerical visuospatial skills. These results demonstrate that relation between basic, even specific, numerical skills and advanced mathematical achievement can be artifactual and explained by visuospatial processing. PMID:26913930

5. High security chaotic multiple access scheme for visible light communication systems with advanced encryption standard interleaving

Qiu, Junchao; Zhang, Lin; Li, Diyang; Liu, Xingcheng

2016-06-01

Chaotic sequences can be applied to realize multiple user access and improve the system security for a visible light communication (VLC) system. However, since the map patterns of chaotic sequences are usually well known, eavesdroppers can possibly derive the key parameters of chaotic sequences and subsequently retrieve the information. We design an advanced encryption standard (AES) interleaving aided multiple user access scheme to enhance the security of a chaotic code division multiple access-based visible light communication (C-CDMA-VLC) system. We propose to spread the information with chaotic sequences, and then the spread information is interleaved by an AES algorithm and transmitted over VLC channels. Since the computation complexity of performing inverse operations to deinterleave the information is high, the eavesdroppers in a high speed VLC system cannot retrieve the information in real time; thus, the system security will be enhanced. Moreover, we build a mathematical model for the AES-aided VLC system and derive the theoretical information leakage to analyze the system security. The simulations are performed over VLC channels, and the results demonstrate the effectiveness and high security of our presented AES interleaving aided chaotic CDMA-VLC system.

6. Langevin spin dynamics based on ab initio calculations: numerical schemes and applications.

PubMed

Rózsa, L; Udvardi, L; Szunyogh, L

2014-05-28

A method is proposed to study the finite-temperature behaviour of small magnetic clusters based on solving the stochastic Landau-Lifshitz-Gilbert equations, where the effective magnetic field is calculated directly during the solution of the dynamical equations from first principles instead of relying on an effective spin Hamiltonian. Different numerical solvers are discussed in the case of a one-dimensional Heisenberg chain with nearest-neighbour interactions. We performed detailed investigations for a monatomic chain of ten Co atoms on top of a Au(0 0 1) surface. We found a spiral-like ground state of the spins due to Dzyaloshinsky-Moriya interactions, while the finite-temperature magnetic behaviour of the system was well described by a nearest-neighbour Heisenberg model including easy-axis anisotropy. PMID:24806308

7. Numerical simulation of the effect of dissipation and phase fluctuation in a direct communication scheme

Li, Fu; Zhang, Jun-Xiang; Zhu, Shi-Yao

2015-06-01

Recently, the direct counterfactual communication protocol, proposed by Salih et al (2013 Phys. Rev. Lett. 110 170502) using a single photon source under ideal conditions (no dissipation, no phase fluctuation and an infinite number of beam splitters), has attracted much interest from a broad range of scientists. In order to put the direct communication protocol into a realistic framework, we numerically simulate the effect of the dissipation and the phase fluctuation with a finite number of beam splitters. Our calculation shows that the dissipation and phase fluctuation will dramatically decrease the reliability and the efficiency of communication, and even corrupt the communication. To counteract the negative effect of dissipation, we propose the balanced dissipation method, which substantially improves the reliability of the protocol at the expense of decreasing communication efficiency. Meanwhile, our theoretical derivation shows that the reliability and efficiency of communication are independent of the input state: a single photon state or a coherent state.

8. Numerical modeling of spray combustion with an advanced VOF method

NASA Technical Reports Server (NTRS)

Chen, Yen-Sen; Shang, Huan-Min; Shih, Ming-Hsin; Liaw, Paul

1995-01-01

This paper summarizes the technical development and validation of a multiphase computational fluid dynamics (CFD) numerical method using the volume-of-fluid (VOF) model and a Lagrangian tracking model which can be employed to analyze general multiphase flow problems with free surface mechanism. The gas-liquid interface mass, momentum and energy conservation relationships are modeled by continuum surface mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed flow regimes. The objectives of the present study are to develop and verify the fractional volume-of-fluid cell partitioning approach into a predictor-corrector algorithm and to demonstrate the effectiveness of the present approach by simulating benchmark problems including laminar impinging jets, shear coaxial jet atomization and shear coaxial spray combustion flows.

9. Effect of spatial configuration of an extended nonlinear Kierstead-Slobodkin reaction-transport model with adaptive numerical scheme.

PubMed

Owolabi, Kolade M; Patidar, Kailash C

2016-01-01

In this paper, we consider the numerical simulations of an extended nonlinear form of Kierstead-Slobodkin reaction-transport system in one and two dimensions. We employ the popular fourth-order exponential time differencing Runge-Kutta (ETDRK4) schemes proposed by Cox and Matthew (J Comput Phys 176:430-455, 2002), that was modified by Kassam and Trefethen (SIAM J Sci Comput 26:1214-1233, 2005), for the time integration of spatially discretized partial differential equations. We demonstrate the supremacy of ETDRK4 over the existing exponential time differencing integrators that are of standard approaches and provide timings and error comparison. Numerical results obtained in this paper have granted further insight to the question 'What is the minimal size of the spatial domain so that the population persists?' posed by Kierstead and Slobodkin (J Mar Res 12:141-147, 1953), with a conclusive remark that the population size increases with the size of the domain. In attempt to examine the biological wave phenomena of the solutions, we present the numerical results in both one- and two-dimensional space, which have interesting ecological implications. Initial data and parameter values were chosen to mimic some existing patterns. PMID:27064984

10. A Numerical Scheme for Special Relativistic Radiation Magnetohydrodynamics Based on Solving the Time-dependent Radiative Transfer Equation

Ohsuga, Ken; Takahashi, Hiroyuki R.

2016-02-01

11. Advances in numerical solutions to integral equations in liquid state theory

Howard, Jesse J.

Solvent effects play a vital role in the accurate description of the free energy profile for solution phase chemical and structural processes. The inclusion of solvent effects in any meaningful theoretical model however, has proven to be a formidable task. Generally, methods involving Poisson-Boltzmann (PB) theory and molecular dynamic (MD) simulations are used, but they either fail to accurately describe the solvent effects or require an exhaustive computation effort to overcome sampling problems. An alternative to these methods are the integral equations (IEs) of liquid state theory which have become more widely applicable due to recent advancements in the theory of interaction site fluids and the numerical methods to solve the equations. In this work a new numerical method is developed based on a Newton-type scheme coupled with Picard/MDIIS routines. To extend the range of these numerical methods to large-scale data systems, the size of the Jacobian is reduced using basis functions, and the Newton steps are calculated using a GMRes solver. The method is then applied to calculate solutions to the 3D reference interaction site model (RISM) IEs of statistical mechanics, which are derived from first principles, for a solute model of a pair of parallel graphene plates at various separations in pure water. The 3D IEs are then extended to electrostatic models using an exact treatment of the long-range Coulomb interactions for negatively charged walls and DNA duplexes in aqueous electrolyte solutions to calculate the density profiles and solution thermodynamics. It is found that the 3D-IEs provide a qualitative description of the density distributions of the solvent species when compared to MD results, but at a much reduced computational effort in comparison to MD simulations. The thermodynamics of the solvated systems are also qualitatively reproduced by the IE results. The findings of this work show the IEs to be a valuable tool for the study and prediction of

12. Advanced numerical methods for three dimensional two-phase flow calculations

SciTech Connect

Toumi, I.; Caruge, D.

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.

13. Propagation of gravity waves through an SPH scheme with numerical diffusive terms

Antuono, M.; Colagrossi, A.; Marrone, S.; Lugni, C.

2011-04-01

Basing on the work by Antuono et al. (2010) [1], an SPH model with numerical diffusive terms (here denoted δ-SPH) is combined with an enhanced treatment of solid boundaries to simulate 2D gravity waves generated by a wave maker and propagating into a basin. Both regular and transient wave systems are considered. In the former, a large number of simulations is performed for different wave steepness and height-to-depth ratio and the results are compared with a BEM Mixed-Eulerian-Lagrangian solver (here denoted BEM-MEL solver). In the latter, the δ-SPH model has been compared with both the experimental measurements available in the literature and with the BEM-MEL solver, at least until the breaking event occurs. The results show a satisfactory agreement between the δ-SPH model, the BEM-MEL solver and the experiments. Finally, the influence of the weakly-compressibility assumption on the SPH results is inspected and a convergence analysis is provided in order to identify the minimal spatial resolution needed to get an accurate representation of gravity waves.

14. PROBABILISTIC SIMULATION OF SUBSURFACE FLUID FLOW: A STUDY USING A NUMERICAL SCHEME

SciTech Connect

Buscheck, Timothy Eric

1980-03-01

There has been an increasing interest in probabilistic modeling of hydrogeologic systems. The classical approach to groundwater modeling has been deterministic in nature, where individual layers and formations are assumed to be uniformly homogeneous. Even in the case of complex heterogeneous systems, the heterogeneities describe the differences in parameter values between various layers, but not within any individual layer. In a deterministic model a single-number is assigned to each hydrogeologic parameter, given a particular scale of interest. However, physically there is no such entity as a truly uniform and homogeneous unit. Single-number representations or deterministic predictions are subject to uncertainties. The approach used in this work models such uncertainties with probabilistic parameters. The resulting statistical distributions of output variables are analyzed. A numerical algorithm, based on axiomatic principles of probability theory, performs arithmetic operations between probability distributions. Two subroutines are developed from the algorithm and incorporated into the computer program TERZAGI, which solves groundwater flow problems in saturated, multi-dimensional systems. The probabilistic computer program is given the name, PROGRES. The algorithm has been applied to study the following problems: one-dimensional flow through homogeneous media, steady-state and transient flow conditions, one-dimensional flow through heterogeneous media, steady-state and transient flow conditions, and two-dimensional steady-stte flow through heterogeneous media. The results are compared with those available in the literature.

15. Eulerian methods for the description of soot: mathematical modeling and numerical scheme

Nguyen, T. T.; Wick, A.; Laurent, F.; Fox, R.; Pitsch, H.

2014-11-01

A development and comparison between numerical methods for soot modeling derived from the population balance equations (PBE) is presented. The soot mechanism includes nucleation, surface growth, oxidation, aggregation and breakage (Mueller et al., Proceed. Combust. Inst., 2009, 2011). For comparison, data from the ethylene premixed flame of Xu et al. (Combust. Flame 108, 1997) over a range of equivalence ratios are used. Two types of methods are introduced. The first is a moment method in which the closure is obtained through a reconstruction of the number density function (NDF). In particular, the NDF can be approximated by a sum of Gamma distribution functions (Yuan et al., J. Aero. Sci. 51, 2012). The second is Eulerian multi-fluid (MF), which is a size discretization method (Laurent et al., Combust. Theory Modelling 5, 2001) considering one or two moments per section. The case of one moment per section is also known as a sectional method. The accuracy of MF methods depends on the number of sections. Eventually, an extension of these two methods considering the surface area as a function of volume is taken into account to describe more precisely the geometry of soot particles. The solutions from these methods are compared with solutions from Monte Carlo method.

16. Efficient Conservative Numerical Schemes for 1D Nonlinear Spherical Diffusion Equations with Applications in Battery Modeling

SciTech Connect

Zeng, Y; Albertus, P; Klein, R; Chaturvedi, N; Kojic, A; Bazant, MZ; Christensen, J

2013-06-07

Mathematical models of batteries which make use of the intercalation of a species into a solid phase need to solve the corresponding mass transfer problem. Because solving this equation can significantly add to the computational cost of a model, various methods have been devised to reduce the computational time. In this paper we focus on a comparison of the formulation, accuracy, and order of the accuracy for two numerical methods of solving the spherical diffusion problem with a constant or non-constant diffusion coefficient: the finite volume method and the control volume method. Both methods provide perfect mass conservation and second order accuracy in mesh spacing, but the control volume method provides the surface concentration directly, has a higher accuracy for a given numbers of mesh points and can also be easily extended to variable mesh spacing. Variable mesh spacing can significantly reduce the number of points that are required to achieve a given degree of accuracy in the surface concentration (which is typically coupled to the other battery equations) by locating more points where the concentration gradients are highest. (C) 2013 The Electrochemical Society. All rights reserved.

17. Numerical optimization design of advanced transonic wing configurations

NASA Technical Reports Server (NTRS)

Cosentino, G. B.; Holst, T. L.

1984-01-01

A computationally efficient and versatile technique for use in the design of advanced transonic wing configurations has been developed. A reliable and fast transonic wing flow-field analysis program, TWING, has been coupled with a modified quasi-Newton method, unconstrained optimization algorithm, QNMDIF, to create a new design tool. Fully three-dimensional wing designs utilizing both specified wing pressure distributions and drag-to-lift ration minimization as design objectives are demonstrated. Because of the high computational efficiency of each of the components of the design code, in particular the vectorization of TWING and the high speed of the Cray X-MP vector computer, the computer time required for a typical wing design is reduced by approximately an order of magnitude over previous methods. In the results presented here, this computed wave drag has been used as the quantity to be optimized (minimized) with great success, yielding wing designs with nearly shock-free (zero wave drag) pressure distributions and very reasonable wing section shapes.

18. Immersed boundary Eulerian-Lagrangian 3D simulation of pyroclastic density currents: numerical scheme and experimental validation

Doronzo, Domenico Maria; de Tullio, Marco; Pascazio, Giuseppe; Dellino, Pierfrancesco

2010-05-01

Pyroclastic density currents are ground hugging, hot, gas-particle flows representing the most hazardous events of explosive volcanism. Their impact on structures is a function of dynamic pressure, which expresses the lateral load that such currents exert over buildings. In this paper we show how analog experiments can be matched with numerical simulations for capturing the essential physics of the multiphase flow. We used an immersed boundary scheme for the mesh generation, which helped in reconstructing the steep velocity and particle concentration gradients near the ground surface. Results show that the calculated values of dynamic pressure agree reasonably with the experimental measurements. These outcomes encourage future application of our method for the assessment of the impact of pyroclastic density currents at the natural scale.

19. Development of comprehensive numerical schemes for predicting evaporating gas-droplets flow processes of a liquid-fueled combustor

NASA Technical Reports Server (NTRS)

Chen, C. P.

1990-01-01

An existing Computational Fluid Dynamics code for simulating complex turbulent flows inside a liquid rocket combustion chamber was validated and further developed. The Advanced Rocket Injector/Combustor Code (ARICC) is simplified and validated against benchmark flow situations for laminar and turbulent flows. The numerical method used in ARICC Code is re-examined for incompressible flow calculations. For turbulent flows, both the subgrid and the two equation k-epsilon turbulence models are studied. Cases tested include idealized Burger's equation in complex geometries and boundaries, a laminar pipe flow, a high Reynolds number turbulent flow, and a confined coaxial jet with recirculations. The accuracy of the algorithm is examined by comparing the numerical results with the analytical solutions as well as experimented data with different grid sizes.

20. Technical note: Improving the AWAT filter with interpolation schemes for advanced processing of high resolution data

Peters, Andre; Nehls, Thomas; Wessolek, Gerd

2016-06-01

Weighing lysimeters with appropriate data filtering yield the most precise and unbiased information for precipitation (P) and evapotranspiration (ET). A recently introduced filter scheme for such data is the AWAT (Adaptive Window and Adaptive Threshold) filter (Peters et al., 2014). The filter applies an adaptive threshold to separate significant from insignificant mass changes, guaranteeing that P and ET are not overestimated, and uses a step interpolation between the significant mass changes. In this contribution we show that the step interpolation scheme, which reflects the resolution of the measuring system, can lead to unrealistic prediction of P and ET, especially if they are required in high temporal resolution. We introduce linear and spline interpolation schemes to overcome these problems. To guarantee that medium to strong precipitation events abruptly following low or zero fluxes are not smoothed in an unfavourable way, a simple heuristic selection criterion is used, which attributes such precipitations to the step interpolation. The three interpolation schemes (step, linear and spline) are tested and compared using a data set from a grass-reference lysimeter with 1 min resolution, ranging from 1 January to 5 August 2014. The selected output resolutions for P and ET prediction are 1 day, 1 h and 10 min. As expected, the step scheme yielded reasonable flux rates only for a resolution of 1 day, whereas the other two schemes are well able to yield reasonable results for any resolution. The spline scheme returned slightly better results than the linear scheme concerning the differences between filtered values and raw data. Moreover, this scheme allows continuous differentiability of filtered data so that any output resolution for the fluxes is sound. Since computational burden is not problematic for any of the interpolation schemes, we suggest always using the spline scheme.

1. Analytical and numerical schemes for a derivative with filtering property and no singular kernel with applications to diffusion

Doungmo Goufo, Emile Franc; Atangana, Abdon

2016-08-01

There have been numbers of conflicting and confusing situations, but also uniformity, in the application of the two most popular fractional derivatives, namely the classic Riemann-Liouville and Caputo fractional derivatives. The range of these issues is wide, including the initialization with the Caputo derivative and its observed difficulties compared to the Riemann-Liouville initialization conditions. In this paper, being aware of these issues and reacting to the newly introduced Caputo-Fabrizio fractional derivative (CFFD) without singular kernel, we introduce a new definition of fractional derivative called the new Riemann-Liouville fractional derivative (NRLFD) without singular kernel. The filtering property of the NRLFD is pointed out by showing it as the derivative of a convolution and contrary to the CFFD, it matches with the function when the order is zero. We also explore various scientific situations that may be conflicting and confusing in the applicability of both new derivatives. In particular, we show that both definitions still have some basic similarities, like not obeying the traditional chain rule. Furthermore, we provide the explicit formula for the Laplace transform of the NRLFD and we prove that, contrary to the CFFD, the NRLFD requires non-constant initial conditions and does not require the function f to be continuous or differentiable. Some simulations for the NRLFD are presented for different values of the derivative order. In the second part of this work, numerical approximations for the first- and second-order NRLFD are developped followed by a concrete application to diffusion. The stability of the numerical scheme is proved and numerical simulations are performed for different values of the derivative order α. They exhibit similar behavior for closed values of α.

2. Self-consistent field theory and numerical scheme for calculating the phase diagram of wormlike diblock copolymers.

PubMed

Jiang, Ying; Chen, Jeff Z Y

2013-10-01

This paper concerns establishing a theoretical basis and numerical scheme for studying the phase behavior of AB diblock copolymers made of wormlike chains. The general idea of a self-consistent field theory is the combination of the mean-field approach together with a statistical weight that describes the configurational properties of a polymer chain. In recent years, this approach has been extensively used for structural prediction of block copolymers, based on the Gaussian-model description of a polymer chain. The wormlike-chain model has played an important role in the description of polymer systems, covering the semiflexible-to-rod crossover of the polymer properties and the highly stretching regime, which the Gaussian-chain model has difficulties to describe. Although the idea of developing a self-consistent field theory for wormlike chains could be traced back to early development in polymer physics, the solution of such a theory has been limited due to technical difficulties. In particular, a challenge has been to develop a numerical algorithm enabling the calculation of the phase diagram containing three-dimensional structures for wormlike AB diblock copolymers. This paper describes a computational algorithm that combines a number of numerical tricks, which can be used for such a calculation. A phase diagram covering major parameter areas was constructed for the wormlike-chain system and reported by us, where the ratio between the total length and the persistence length of a constituent polymer is suggested as another tuning parameter for the microphase-separated structures; all detailed technical issues are carefully addressed in the current paper. PMID:24229202

3. BOOK REVIEW: Advanced Topics in Computational Partial Differential Equations: Numerical Methods and Diffpack Programming

Katsaounis, T. D.

2005-02-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. The first chapter is an introduction to parallel processing. It covers fundamentals of parallel processing in a simple and concrete way and no prior knowledge of the subject is required. Examples of parallel implementation of basic linear algebra operations are presented using the Message Passing Interface (MPI) programming environment. Here, some knowledge of MPI routines is required by the reader. Examples solving in parallel simple PDEs using

4. A New Cell-Centered Implicit Numerical Scheme for Ions in the 2-D Axisymmetric Code Hall2de

NASA Technical Reports Server (NTRS)

Lopez Ortega, Alejandro; Mikellides, Ioannis G.

2014-01-01

We present a new algorithm in the Hall2De code to simulate the ion hydrodynamics in the acceleration channel and near plume regions of Hall-effect thrusters. This implementation constitutes an upgrade of the capabilities built in the Hall2De code. The equations of mass conservation and momentum for unmagnetized ions are solved using a conservative, finite-volume, cell-centered scheme on a magnetic-field-aligned grid. Major computational savings are achieved by making use of an implicit predictor/multi-corrector algorithm for time evolution. Inaccuracies in the prediction of the motion of low-energy ions in the near plume in hydrodynamics approaches are addressed by implementing a multi-fluid algorithm that tracks ions of different energies separately. A wide range of comparisons with measurements are performed to validate the new ion algorithms. Several numerical experiments with the location and value of the anomalous collision frequency are also presented. Differences in the plasma properties in the near-plume between the single fluid and multi-fluid approaches are discussed. We complete our validation by comparing predicted erosion rates at the channel walls of the thruster with measurements. Erosion rates predicted by the plasma properties obtained from simulations replicate accurately measured rates of erosion within the uncertainty range of the sputtering models employed.

5. Recent advances in numerical simulation and control of asymmetric flows around slender bodies

NASA Technical Reports Server (NTRS)

Kandil, Osama A.; Wong, Tin-Chee; Sharaf, Hazem H.; Liu, C. H.

1992-01-01

The problems of asymmetric flow around slender bodies and its control are formulated using the unsteady, compressible, thin-layer or full Navier-Stokes equations which are solved using an implicit, flux-difference splitting, finite-volume scheme. The problem is numerically simulated for both locally-conical and three-dimensional flows. The numerical applications include studies of the effects of relative incidence, Mach number and Reynolds number on the flow asymmetry. For the control of flow asymmetry, the numerical simulation cover passive and active control methods. For the passive control, the effectiveness of vertical fins placed in the leeward plane of geometric symmetry and side strakes with different orientations is studied. For the active control, the effectiveness of normal and tangential flow injection and surface heating and a combination of these methods is studied.

6. A Straightforward Method for Advance Estimation of User Charges for Information in Numeric Databases.

ERIC Educational Resources Information Center

Jarvelin, Kalervo

1986-01-01

Describes a method for advance estimation of user charges for queries in relational data model-based numeric databases when charges are based on data retrieved. Use of this approach is demonstrated by sample queries to an imaginary marketing database. The principles and methods of this approach and its relevance are discussed. (MBR)

7. New energy-preserving schemes using Hamiltonian Boundary Value and Fourier pseudospectral methods for the numerical solution of the "good" Boussinesq equation

Yan, Jinliang; Zhang, Zhiyue

2016-04-01

Two energy-preserving schemes are proposed for the "good" Boussinesq (GBq) equation using the Hamiltonian Boundary Value and Fourier pseudospectral methods. The equation is discretized in space by Fourier pseudospectral method and in time by Hamiltonian Boundary Value methods (HBVMs). The outstanding advantages of the proposed schemes are that they can precisely conserve the global mass and energy, and provide highly accurate results. The single solitary wave, the interaction of two solitary waves and the birth of solitary waves are presented to validate the accuracy and conservation properties of the proposed schemes. In addition, we also compare our numerical results with other known studied methods in terms of numerical accuracy and conservation properties.

8. Numerical framework and performance of the new multiple-phase cloud microphysics scheme in RegCM4.5: precipitation, cloud microphysics, and cloud radiative effects

Nogherotto, Rita; Tompkins, Adrian Mark; Giuliani, Graziano; Coppola, Erika; Giorgi, Filippo

2016-07-01

We implement and evaluate a new parameterization scheme for stratiform cloud microphysics and precipitation within regional climate model RegCM4. This new parameterization is based on a multiple-phase one-moment cloud microphysics scheme built upon the implicit numerical framework recently developed and implemented in the ECMWF operational forecasting model. The parameterization solves five prognostic equations for water vapour, cloud liquid water, rain, cloud ice, and snow mixing ratios. Compared to the pre-existing scheme, it allows a proper treatment of mixed-phase clouds and a more physically realistic representation of cloud microphysics and precipitation. Various fields from a 10-year long integration of RegCM4 run in tropical band mode with the new scheme are compared with their counterparts using the previous cloud scheme and are evaluated against satellite observations. In addition, an assessment using the Cloud Feedback Model Intercomparison Project (CFMIP) Observational Simulator Package (COSP) for a 1-year sub-period provides additional information for evaluating the cloud optical properties against satellite data. The new microphysics parameterization yields an improved simulation of cloud fields, and in particular it removes the overestimation of upper level cloud characteristics of the previous scheme, increasing the agreement with observations and leading to an amelioration of a long-standing problem in the RegCM system. The vertical cloud profile produced by the new scheme leads to a considerably improvement of the representation of the longwave and shortwave components of the cloud radiative forcing.

9. Developing a TPCK-SRL Assessment Scheme for Conceptually Advancing Technology in Education

ERIC Educational Resources Information Center

Kohen, Zehavit; Kramarski, Bracha

2012-01-01

The present study aimed to: (a) develop a conceptual TPCK-SRL scheme for assessing teachers' integration of self-regulated learning (SRL) considerations while infusing technology into a TPCK classroom context (blending K = knowledge about T = technology, P = pedagogy, and C = content), which reflects all three knowledge components' dynamic…

10. Physical control oriented model of large scale refrigerators to synthesize advanced control schemes. Design, validation, and first control results

Bonne, François; Alamir, Mazen; Bonnay, Patrick

2014-01-01

In this paper, a physical method to obtain control-oriented dynamical models of large scale cryogenic refrigerators is proposed, in order to synthesize model-based advanced control schemes. These schemes aim to replace classical user experience designed approaches usually based on many independent PI controllers. This is particularly useful in the case where cryoplants are submitted to large pulsed thermal loads, expected to take place in the cryogenic cooling systems of future fusion reactors such as the International Thermonuclear Experimental Reactor (ITER) or the Japan Torus-60 Super Advanced Fusion Experiment (JT-60SA). Advanced control schemes lead to a better perturbation immunity and rejection, to offer a safer utilization of cryoplants. The paper gives details on how basic components used in the field of large scale helium refrigeration (especially those present on the 400W @1.8K helium test facility at CEA-Grenoble) are modeled and assembled to obtain the complete dynamic description of controllable subsystems of the refrigerator (controllable subsystems are namely the Joule-Thompson Cycle, the Brayton Cycle, the Liquid Nitrogen Precooling Unit and the Warm Compression Station). The complete 400W @1.8K (in the 400W @4.4K configuration) helium test facility model is then validated against experimental data and the optimal control of both the Joule-Thompson valve and the turbine valve is proposed, to stabilize the plant under highly variable thermals loads. This work is partially supported through the European Fusion Development Agreement (EFDA) Goal Oriented Training Program, task agreement WP10-GOT-GIRO.

11. Physical control oriented model of large scale refrigerators to synthesize advanced control schemes. Design, validation, and first control results

SciTech Connect

Bonne, François; Bonnay, Patrick

2014-01-29

In this paper, a physical method to obtain control-oriented dynamical models of large scale cryogenic refrigerators is proposed, in order to synthesize model-based advanced control schemes. These schemes aim to replace classical user experience designed approaches usually based on many independent PI controllers. This is particularly useful in the case where cryoplants are submitted to large pulsed thermal loads, expected to take place in the cryogenic cooling systems of future fusion reactors such as the International Thermonuclear Experimental Reactor (ITER) or the Japan Torus-60 Super Advanced Fusion Experiment (JT-60SA). Advanced control schemes lead to a better perturbation immunity and rejection, to offer a safer utilization of cryoplants. The paper gives details on how basic components used in the field of large scale helium refrigeration (especially those present on the 400W @1.8K helium test facility at CEA-Grenoble) are modeled and assembled to obtain the complete dynamic description of controllable subsystems of the refrigerator (controllable subsystems are namely the Joule-Thompson Cycle, the Brayton Cycle, the Liquid Nitrogen Precooling Unit and the Warm Compression Station). The complete 400W @1.8K (in the 400W @4.4K configuration) helium test facility model is then validated against experimental data and the optimal control of both the Joule-Thompson valve and the turbine valve is proposed, to stabilize the plant under highly variable thermals loads. This work is partially supported through the European Fusion Development Agreement (EFDA) Goal Oriented Training Program, task agreement WP10-GOT-GIRO.

12. The effect of a sub-grid statistical cloud-cover scheme applied to the COSMO local numerical weather prediction model over the wider geographical domain of Greece

Avgoustoglou, Euripides; Tzeferi, Theodora

2015-01-01

The COSMO model uses operationally two sub-grid schemes for the evaluation of stratus clouds. A semi-empirical scheme based on relative humidity is used in the radiation module while a statistical scheme is used in the turbulence module. The objective is to investigate the possibility of the implementation of the statistical scheme also in the radiation module. The relative impact is presented in reference to a spring test case with synoptic conditions that favor stratiform clouds. The domain considered is the wider Balkan region around the Hellenic geographical area and is characterized by comparable sea and land partitions. This particular domain choice gives rise to a strong coexistence of continental as well as marine clouds which is one of the most challenging features regarding the operational use of numerical weather prediction models by the Hellenic Meteorological Service. The results are evaluated through direct comparisons with satellite data as well as the observed 2-m temperatures for an approximate total of fifty Greek synoptic meteorological stations. The implementation of the statistical scheme led to an underestimation of low cloud-cover by the model in contrast to the implementation of the default relative-humidity scheme, while regarding medium cloud-cover the situation was reversed. Also, the daily 2-meter minimum and maximum temperatures were slightly better simulated, but not conclusively, when the statistical scheme was implemented in the radiation module. Although the statistical scheme cannot in its present form replace operationally the relative-humidity scheme in the radiation module, it is an important asset to COSMO model invoking valuable insight to the physics of the model and can be used as a basis to support the ongoing research in this crucial area of atmospheric sciences.

13. Particle exhaust schemes in the DIII-D advanced divertor configuration

SciTech Connect

Menon, M.M.; Mioduszewski, P.K.

1989-01-01

For density control in long-pulse operation, the open divertor on the DIII-D tokamak will be equipped with a baffled chamber and a pumping system. The throat of the baffle chamber is sized to provide optimal pumping for the typical plasma equilibrium configuration. Severe limitations on the toroidal conductance of this baffle chamber require the use of in-vessel pumping to achieve the desired particle exhaust of about 25 Torr{center dot}l/s. Two separate pumping schemes are considered: an array of titanium getter modules based on the design developed by the Tore Supra team and a cryocondensation pump. The merits and demerits of each scheme are analyzed, and the design considerations introduced by the tokamak environment are brought out. 3 refs., 5 figs.

14. Left Ventricular Flow Analysis: Recent Advances in Numerical Methods and Applications in Cardiac Ultrasound

PubMed Central

Borazjani, Iman; Westerdale, John; McMahon, Eileen M.; Rajaraman, Prathish K.; Heys, Jeffrey J.

2013-01-01

The left ventricle (LV) pumps oxygenated blood from the lungs to the rest of the body through systemic circulation. The efficiency of such a pumping function is dependent on blood flow within the LV chamber. It is therefore crucial to accurately characterize LV hemodynamics. Improved understanding of LV hemodynamics is expected to provide important clinical diagnostic and prognostic information. We review the recent advances in numerical and experimental methods for characterizing LV flows and focus on analysis of intraventricular flow fields by echocardiographic particle image velocimetry (echo-PIV), due to its potential for broad and practical utility. Future research directions to advance patient-specific LV simulations include development of methods capable of resolving heart valves, higher temporal resolution, automated generation of three-dimensional (3D) geometry, and incorporating actual flow measurements into the numerical solution of the 3D cardiovascular fluid dynamics. PMID:23690874

15. Recent advances in Runge-Kutta schemes for solving 3-D Navier-Stokes equations

NASA Technical Reports Server (NTRS)

Vatsa, Veer N.; Wedan, Bruce W.; Abid, Ridha

1989-01-01

A thin-layer Navier-Stokes has been developed for solving high Reynolds number, turbulent flows past aircraft components under transonic flow conditions. The computer code has been validated through data comparisons for flow past isolated wings, wing-body configurations, prolate spheroids and wings mounted inside wind-tunnels. The basic code employs an explicit Runge-Kutta time-stepping scheme to obtain steady state solution to the unsteady governing equations. Significant gain in the efficiency of the code has been obtained by implementing a multigrid acceleration technique to achieve steady-state solutions. The improved efficiency of the code has made it feasible to conduct grid-refinement and turbulence model studies in a reasonable amount of computer time. The non-equilibrium turbulence model of Johnson and King has been extended to three-dimensional flows and excellent agreement with pressure data has been obtained for transonic separated flow over a transport type of wing.

16. Advanced control schemes and kinematic analysis for a kinematically redundant 7 DOF manipulator

NASA Technical Reports Server (NTRS)

Nguyen, Charles C.; Zhou, Zhen-Lei

1990-01-01

The kinematic analysis and control of a kinematically redundant manipulator is addressed. The manipulator is the slave arm of a telerobot system recently built at Goddard Space Flight Center (GSFC) to serve as a testbed for investigating research issues in telerobotics. A forward kinematic transformation is developed in its most simplified form, suitable for real-time control applications, and the manipulator Jacobian is derived using the vector cross product method. Using the developed forward kinematic transformation and quaternion representation of orientation matrices, we perform computer simulation to evaluate the efficiency of the Jacobian in converting joint velocities into Cartesian velocities and to investigate the accuracy of Jacobian pseudo-inverse for various sampling times. The equivalence between Cartesian velocities and quaternion is also verified using computer simulation. Three control schemes are proposed and discussed for controlling the motion of the slave arm end-effector.

SciTech Connect

1995-12-31

The possibility to operate with a short pulse FEL (electron pulse length a {sigma} {approx_equal} 1mm) in far infrared region (with the radiation wavelength 10{mu}m {le}{lambda}{le}300{mu}m) gives rise to many new applications. The parameter {mu}{sub c} = N{lambda}/{sigma}, which characterises the interaction time between the electron and radiation pulses in the undulator with N periods, becomes rather large ({mu}{sub c}{approximately}10) even for {lambda}{approximately}100{mu}m. An increase in {lambda} leads to the reduction in output power owing to an increase in both the diffraction losses and {mu}{sub c}. We consider two options capable to improve the situation in the FEL with planar undulator geometry: I. Reduction of the diffraction losses by use of circular or rectangular hollow waveguides in the cavity. II. Effective reduction of {mu}{sub c} factor by increasing the FEL-interaction time by means of the optical klystron undulator scheme. We demonstrate that; I. Introduction of a waveguide into the cavity causes additional difficulties. In both circular and rectangular cases the waveguide dispersion leads to the undesirable variation of radiation pulse shape. For the circular case the effect of depolarization becomes important and leads to the additional losses (up to 30%) due to the violation of the axial symmetry and linear polarization. II. The optical klystron scheme consisting of several separated undulators, can indeed increase the FEL-interaction time. The right choice of electron bunch delay between these undulators gives the possibility to couple several successive electron pulses with one radiation pulse during one passage through the cavity. This results in considerable increase in FEL-interaction time and corresponding increase in gain.

18. Numerical simulation of the reactive flow in advanced (HSR) combustors using KIVA-2

NASA Technical Reports Server (NTRS)

Winowich, Nicholas S.

1991-01-01

Recent work has been done with the goal of establishing ultralow emission aircraft gas turbine combustors. A significant portion of the effort is the development of three dimensional computational combustor models. The KIVA-II computer code which is based on the Implicit Continuous Eulerian Difference mesh Arbitrary Lagrangian Eulerian (ICED-ALE) numerical scheme is one of the codes selected by NASA to achieve these goals. This report involves a simulation of jet injection through slanted slots within the Rich burn/Quick quench/Lean burn (RQL) baseline experimental rig. The RQL combustor distinguishes three regions of combustion. This work specifically focuses on modeling the quick quench mixer region in which secondary injection air is introduced radially through 12 equally spaced slots around the mixer circumference. Steady state solutions are achieved with modifications to the KIVA-II program. Work currently underway will evaluate thermal mixing as a function of injection air velocity and angle of inclination of the slots.

19. An efficient liner cooling scheme for advanced small gas turbine combustors

NASA Technical Reports Server (NTRS)

Paskin, Marc D.; Mongia, Hukam C.; Acosta, Waldo A.

1993-01-01

A joint Army/NASA program was conducted to design, fabricate, and test an advanced, small gas turbine, reverse-flow combustor utilizing a compliant metal/ceramic (CMC) wall cooling concept. The objectives of this effort were to develop a design method (basic design data base and analysis) for the CMC cooling technique and then demonstrate its application to an advanced cycle, small, reverse-flow combustor with 3000 F burner outlet temperature. The CMC concept offers significant improvements in wall cooling effectiveness resulting in a large reduction in cooling air requirements. Therefore, more air is available for control of burner outlet temperature pattern in addition to the benefits of improved efficiency, reduced emissions, and lower smoke levels. The program was divided into four tasks. Task 1 defined component materials and localized design of the composite wall structure in conjunction with development of basic design models for the analysis of flow and heat transfer through the wall. Task 2 included implementation of the selected materials and validated design models during combustor preliminary design. Detail design of the selected combustor concept and its refinement with 3D aerothermal analysis were completed in Task 3. Task 4 covered detail drawings, process development and fabrication, and a series of burner rig tests. The purpose of this paper is to provide details of the investigation into the fundamental flow and heat transfer characteristics of the CMC wall structure as well as implementation of the fundamental analysis method for full-scale combustor design.

20. Numerical simulation of flood inundation using a well-balanced kinetic scheme for the shallow water equations with bulk recharge and discharge

Ersoy, Mehmet; Lakkis, Omar; Townsend, Philip

2016-04-01

The flow of water in rivers and oceans can, under general assumptions, be efficiently modelled using Saint-Venant's shallow water system of equations (SWE). SWE is a hyperbolic system of conservation laws (HSCL) which can be derived from a starting point of incompressible Navier-Stokes. A common difficulty in the numerical simulation of HSCLs is the conservation of physical entropy. Work by Audusse, Bristeau, Perthame (2000) and Perthame, Simeoni (2001), proposed numerical SWE solvers known as kinetic schemes (KSs), which can be shown to have desirable entropy-consistent properties, and are thus called well-balanced schemes. A KS is derived from kinetic equations that can be integrated into the SWE. In flood risk assessment models the SWE must be coupled with other equations describing interacting meteorological and hydrogeological phenomena such as rain and groundwater flows. The SWE must therefore be appropriately modified to accommodate source and sink terms, so kinetic schemes are no longer valid. While modifications of SWE in this direction have been recently proposed, e.g., Delestre (2010), we depart from the extant literature by proposing a novel model that is "entropy-consistent" and naturally extends the SWE by respecting its kinetic formulation connections. This allows us to derive a system of partial differential equations modelling flow of a one-dimensional river with both a precipitation term and a groundwater flow model to account for potential infiltration and recharge. We exhibit numerical simulations of the corresponding kinetic schemes. These simulations can be applied to both real world flood prediction and the tackling of wider issues on how climate and societal change are affecting flood risk.

1. The HiPER project for inertial confinement fusion and some experimental results on advanced ignition schemes

Batani, D.; Koenig, M.; Baton, S.; Perez, F.; Gizzi, L. A.; Koester, P.; Labate, L.; Honrubia, J.; Antonelli, L.; Morace, A.; Volpe, L.; Santos, J.; Schurtz, G.; Hulin, S.; Ribeyre, X.; Fourment, C.; Nicolai, P.; Vauzour, B.; Gremillet, L.; Nazarov, W.; Pasley, J.; Richetta, M.; Lancaster, K.; Spindloe, Ch; Tolley, M.; Neely, D.; Kozlová, M.; Nejdl, J.; Rus, B.; Wolowski, J.; Badziak, J.; Dorchies, F.

2011-12-01

This paper presents the goals and some of the results of experiments conducted within the Working Package 10 (Fusion Experimental Programme) of the HiPER Project. These experiments concern the study of the physics connected to 'advanced ignition schemes', i.e. the fast ignition and the shock ignition approaches to inertial fusion. Such schemes are aimed at achieving a higher gain, as compared with the classical approach which is used in NIF, as required for future reactors, and make fusion possible with smaller facilities. In particular, a series of experiments related to fast ignition were performed at the RAL (UK) and LULI (France) Laboratories and studied the propagation of fast electrons (created by a short-pulse ultra-high-intensity beam) in compressed matter, created either by cylindrical implosions or by compression of planar targets by (planar) laser-driven shock waves. A more recent experiment was performed at PALS and investigated the laser-plasma coupling in the 1016 W cm-2 intensity regime of interest for shock ignition.

2. Numerical simulation of a dust event in northeastern Germany with a new dust emission scheme in COSMO-ART

Technology Transfer Automated Retrieval System (TEKTRAN)

The dust emission scheme of Shao (2004) has been implemented into the regional atmospheric model COSMO-ART and has been applied to a severe dust event in northeastern Germany on 8th April 2011. The model sensitivity to soil moisture and vegetation cover has been studied. Soil moisture has been found...

3. Recent numerical and algorithmic advances within the volume tracking framework for modeling interfacial flows

SciTech Connect

François, Marianne M.

2015-05-28

A review of recent advances made in numerical methods and algorithms within the volume tracking framework is presented. The volume tracking method, also known as the volume-of-fluid method has become an established numerical approach to model and simulate interfacial flows. Its advantage is its strict mass conservation. However, because the interface is not explicitly tracked but captured via the material volume fraction on a fixed mesh, accurate estimation of the interface position, its geometric properties and modeling of interfacial physics in the volume tracking framework remain difficult. Several improvements have been made over the last decade to address these challenges. In this study, the multimaterial interface reconstruction method via power diagram, curvature estimation via heights and mean values and the balanced-force algorithm for surface tension are highlighted.

4. Recent numerical and algorithmic advances within the volume tracking framework for modeling interfacial flows

DOE PAGESBeta

François, Marianne M.

2015-05-28

A review of recent advances made in numerical methods and algorithms within the volume tracking framework is presented. The volume tracking method, also known as the volume-of-fluid method has become an established numerical approach to model and simulate interfacial flows. Its advantage is its strict mass conservation. However, because the interface is not explicitly tracked but captured via the material volume fraction on a fixed mesh, accurate estimation of the interface position, its geometric properties and modeling of interfacial physics in the volume tracking framework remain difficult. Several improvements have been made over the last decade to address these challenges.more » In this study, the multimaterial interface reconstruction method via power diagram, curvature estimation via heights and mean values and the balanced-force algorithm for surface tension are highlighted.« less

5. An energy and potential enstrophy conserving numerical scheme for the multi-layer shallow water equations with complete Coriolis force

Stewart, Andrew L.; Dellar, Paul J.

2016-05-01

We present an energy- and potential enstrophy-conserving scheme for the non-traditional shallow water equations that include the complete Coriolis force and topography. These integral conservation properties follow from material conservation of potential vorticity in the continuous shallow water equations. The latter property cannot be preserved by a discretisation on a fixed Eulerian grid, but exact conservation of a discrete energy and a discrete potential enstrophy seems to be an effective substitute that prevents any distortion of the forward and inverse cascades in quasi-two dimensional turbulence through spurious sources and sinks of energy and potential enstrophy, and also increases the robustness of the scheme against nonlinear instabilities. We exploit the existing Arakawa-Lamb scheme for the traditional shallow water equations, reformulated by Salmon as a discretisation of the Hamiltonian and Poisson bracket for this system. The non-rotating, traditional, and our non-traditional shallow water equations all share the same continuous Hamiltonian structure and Poisson bracket, provided one distinguishes between the particle velocity and the canonical momentum per unit mass. We have determined a suitable discretisation of the non-traditional canonical momentum, which includes additional coupling between the layer thickness and velocity fields, and modified the discrete kinetic energy to suppress an internal symmetric computational instability that otherwise arises for multiple layers. The resulting scheme exhibits the expected second-order convergence under spatial grid refinement. We also show that the drifts in the discrete total energy and potential enstrophy due to temporal truncation error may be reduced to machine precision under suitable refinement of the timestep using the third-order Adams-Bashforth or fourth-order Runge-Kutta integration schemes.

6. Numerical study of Alfvén eigenmodes in the Experimental Advanced Superconducting Tokamak

SciTech Connect

Hu, Youjun; Li, Guoqiang; Yang, Wenjun; Zhou, Deng; Ren, Qilong; Gorelenkov, N. N.; Cai, Huishan

2014-05-15

Alfvén eigenmodes in up-down asymmetric tokamak equilibria are studied by a new magnetohydrodynamic eigenvalue code. The code is verified with the NOVA code for the Solovév equilibrium and then is used to study Alfvén eigenmodes in a up-down asymmetric equilibrium of the Experimental Advanced Superconducting Tokamak. The frequency and mode structure of toroidicity-induced Alfvén eigenmodes are calculated. It is demonstrated numerically that up-down asymmetry induces phase variation in the eigenfunction across the major radius on the midplane.

7. Simulation studies of the impact of advanced observing systems on numerical weather prediction

NASA Technical Reports Server (NTRS)

Atlas, R.; Kalnay, E.; Susskind, J.; Reuter, D.; Baker, W. E.; Halem, M.

1984-01-01

To study the potential impact of advanced passive sounders and lidar temperature, pressure, humidity, and wind observing systems on large-scale numerical weather prediction, a series of realistic simulation studies between the European Center for medium-range weather forecasts, the National Meteorological Center, and the Goddard Laboratory for Atmospheric Sciences is conducted. The project attempts to avoid the unrealistic character of earlier simulation studies. The previous simulation studies and real-data impact tests are reviewed and the design of the current simulation system is described. Consideration is given to the simulation of observations of space-based sounding systems.

8. Advanced integration schemes for high-functionality/high-performance photonic integrated circuits

Raring, James W.; Sysak, Matthew N.; Tauke-Pedretti, Anna; Dummer, Matthew; Skogen, Erik J.; Barton, Jonathon S.; DenBaars, S. P.; Coldren, Larry A.

2006-02-01

The evolution of optical communication systems has facilitated the required bandwidth to meet the increasing data rate demands. However, as the peripheral technologies have progressed to meet the requirements of advanced systems, an abundance of viable solutions and products have emerged. The finite market for these products will inevitably force a paradigm shift upon the communications industry. Monolithic integration is a key technology that will facilitate this shift as it will provide solutions at low cost with reduced power dissipation and foot-print in the form of highly functional optical components based on photonic integrated circuits (PICs). In this manuscript, we discuss the advantages, potential applications, and challenges of photonic integration. After a brief overview of various integration techniques, we present our novel approaches to increase the performance of the individual components comprising highly functional PICs.

9. Advanced scheme for high-yield laser driven proton-boron fusion reaction

Margarone, D.; Picciotto, A.; Velyhan, A.; Krasa, J.; Kucharik, M.; Morrissey, M.; Mangione, A.; Szydlowsky, A.; Malinowska, A.; Bertuccio, G.; Shi, Y.; Crivellari, M.; Ullschmied, J.; Bellutti, P.; Korn, G.

2015-02-01

A low contrast nanosecond laser pulse with relatively low intensity (3 × 1016 W cm-2) was used to enhance the yield of induced nuclear reactions in advanced solid targets. In particular the "ultraclean" proton-boron fusion reaction, producing energetic alpha-particles without neutron generation, was chosen. A spatially well-defined layer of boron dopants in a hydrogen-enriched silicon substrate was used as target. The combination of the specific target geometry and the laser pulse temporal shape allowed enhancing the yield of alpha-particles up to 109 per steradian, i.e 100 times higher than previous experimental achievements. Moreover the alpha particle stream presented a clearly peaked angular and energy distribution, which make this secondary source attractive for potential applications. This result can be ascribed to the interaction of the long laser pre-pulse with the target and to the optimal target geometry and composition.

10. A membrane-integrated advanced scheme for treatment of industrial wastewater: dynamic modeling towards scale up.

PubMed

Kumar, Ramesh; Pal, Parimal

2013-08-01

Modeling and simulation was carried out for an advanced membrane-integrated hybrid treatment process that ensures reuse of water with conversion and recovery of ammoniacal nitrogen as value-added struvite fertilizer from coke wastewater. While toxic cyanide was largely removed in a pre-chemical treatment unit using Fenton's reagents under optimized conditions, more than 95% of NH4(+)-N could be recovered as a valuable by-product called struvite through addition of appropriate doses of magnesium and phosphate salts. Water could be turned reusable through a polishing treatment by nanofiltration membranes in a largely fouling free membrane module following a biodegradation step. Mathematical modeling of such an integrated process was done with Haldane-Andrew approach for the associated microbial degradation of phenol by Pseudomonas putida. Residual NH4(+) was degraded by nitrification and denitrification following the modified Monod kinetics. The model could successfully predict the plant performance as reflected in reasonably low relative error (0.03-0.18) and high Willmott d-index (>0.98). PMID:23735488

11. A review on recent advances in the numerical simulation for coalbed-methane-recovery process

SciTech Connect

Wei, X.R.; Wang, G.X.; Massarotto, P.; Golding, S.D.; Rudolph, V.

2007-12-15

The recent advances in numerical simulation for primary coalbed methane (CBM) recovery and enhanced coalbed-methane recovery (ECBMR) processes are reviewed, primarily focusing on the progress that has occurred since the late 1980s. Two major issues regarding the numerical modeling will be discussed in this review: first, multicomponent gas transport in in-situ bulk coal and, second, changes of coal properties during methane (CH{sub 4}) production. For the former issues, a detailed review of more recent advances in modeling gas and water transport within a coal matrix is presented. Further, various factors influencing gas diffusion through the coal matrix will be highlighted as well, such as pore structure, concentration and pressure, and water effects. An ongoing bottleneck for evaluating total mass transport rate is developing a reasonable representation of multiscale pore space that considers coal type and rank. Moreover, few efforts have been concerned with modeling water-flow behavior in the coal matrix and its effects on CH{sub 4} production and on the exchange of carbon dioxide (CO{sub 2}) and CH{sub 4}. As for the second issue, theoretical coupled fluid-flow and geomechanical models have been proposed to describe the evolution of pore structure during CH{sub 4} production, instead of traditional empirical equations. However, there is currently no effective coupled model for engineering applications. Finally, perspectives on developing suitable simulation models for CBM production and for predicting CO{sub 2}-sequestration ECBMR are suggested.

12. Using the PPML approach for constructing a low-dissipation, operator-splitting scheme for numerical simulations of hydrodynamic flows

Kulikov, Igor; Vorobyov, Eduard

2016-07-01

An approach for constructing a low-dissipation numerical method is described. The method is based on a combination of the operator-splitting method, Godunov method, and piecewise-parabolic method on the local stencil. Numerical method was tested on a standard suite of hydrodynamic test problems. In addition, the performance of the method is demonstrated on a global test problem showing the development of a spiral structure in a gravitationally unstable gaseous galactic disk.

13. An advanced numerical model for phase change problems in complicated geometries

Khashan, Saud Abdel-Aziz

1998-11-01

An advanced fixed-grid enthalpy formulation based finite volume numerical method is developed to solve the phase change problems in complicated geometries. The numerical method is based on a general non-orthogonal grid structure and a colocated arrangement of variables. Second order discretizations and interpolations are used. The convergence rate is considerably accelerated by switching-off the velocity in the solidified region in an implicit way. This switching-off technique has a strong compatibility with SIMPLE-like methods. For all test cases conducted in this study, the rate of convergence using the new treatment exceeds that of the other enthalpy formulation-based methods and with less numerical stability constraints, when used in convection-diffusion phase change problems. For better run in vector computers, The Incomplete LU decomposition (ILU) matrix solver is partially vectorized. The Mflops (million floating point operation per second) number is raised from 60 to over 300. Water freezing in orthogonal and non-orthogonal geometry are studied under the effect of density inversion. All thermo-physical properties of the water are dealt with as temperature-dependent (no Boussinsq approximation). The results show a profound effect of density inversion on the flow/energy field and on the local as well as on the universal freezing rate.

14. Improved scheme for parametrization of convection in the Met Office's Numerical Atmospheric-dispersion Modelling Environment (NAME)

Meneguz, Elena; Thomson, David; Witham, Claire; Kusmierczyk-Michulec, Jolanta

2015-04-01

NAME is a Lagrangian atmospheric dispersion model used by the Met Office to predict the dispersion of both natural and man-made contaminants in the atmosphere, e.g. volcanic ash, radioactive particles and chemical species. Atmospheric convection is responsible for transport and mixing of air resulting in a large exchange of heat and energy above the boundary layer. Although convection can transport material through the whole troposphere, convective clouds have a small horizontal length scale (of the order of few kilometres). Therefore, for large-scale transport the horizontal scale on which the convection exists is below the global NWP resolution used as input to NAME and convection must be parametrized. Prior to the work presented here, the enhanced vertical mixing generated by non-resolved convection was reproduced by randomly redistributing Lagrangian particles between the cloud base and cloud top with probability equal to 1/25th of the NWP predicted convective cloud fraction. Such a scheme is essentially diffusive and it does not make optimal use of all the information provided by the driving meteorological model. To make up for these shortcomings and make the parametrization more physically based, the convection scheme has been recently revised. The resulting version, presented in this paper, is now based on the balance equation between upward, entrainment and detrainment fluxes. In particular, upward mass fluxes are calculated with empirical formulas derived from Cloud Resolving Models and using the NWP convective precipitation diagnostic as closure. The fluxes are used to estimate how many particles entrain, move upward and detrain. Lastly, the scheme is completed by applying a compensating subsidence flux. The performance of the updated convection scheme is benchmarked against available observational data of passive tracers. In particular, radioxenon is a noble gas that can undergo significant long range transport: this study makes use of observations of

15. Explicit robust schemes for implementation of a class of principal value-based constitutive models: Symbolic and numeric implementation

NASA Technical Reports Server (NTRS)

Arnold, S. M.; Saleeb, A. F.; Tan, H. Q.; Zhang, Y.

1993-01-01

The issue of developing effective and robust schemes to implement a class of the Ogden-type hyperelastic constitutive models is addressed. To this end, special purpose functions (running under MACSYMA) are developed for the symbolic derivation, evaluation, and automatic FORTRAN code generation of explicit expressions for the corresponding stress function and material tangent stiffness tensors. These explicit forms are valid over the entire deformation range, since the singularities resulting from repeated principal-stretch values have been theoretically removed. The required computational algorithms are outlined, and the resulting FORTRAN computer code is presented.

16. Numerical approach for the voloxidation process of an advanced spent fuel conditioning process (ACP)

SciTech Connect

Park, Byung Heung; Jeong, Sang Mun; Seo, Chung-Seok

2007-07-01

A voloxidation process is adopted as the first step of an advanced spent fuel conditioning process in order to prepare the SF oxide to be reduced in the following electrolytic reduction process. A semi-batch type voloxidizer was devised to transform a SF pellet into powder. In this work, a simple reactor model was developed for the purpose of correlating a gas phase flow rate with an operation time as a numerical approach. With an assumption that a solid phase and a gas phase are homogeneous in a reactor, a reaction rate for an oxidation was introduced into a mass balance equation. The developed equation can describe a change of an outlet's oxygen concentration including such a case that a gas flow is not sufficient enough to continue a reaction at its maximum reaction rate. (authors)

17. Direct numerical simulation of transitional and turbulent flow over a heated flat plate using finite-difference schemes

NASA Technical Reports Server (NTRS)

1995-01-01

This report deals with the direct numerical simulation of transitional and turbulent flow at low Mach numbers using high-order-accurate finite-difference techniques. A computation of transition to turbulence of the spatially-evolving boundary layer on a heated flat plate in the presence of relatively high freestream turbulence was performed. The geometry and flow conditions were chosen to match earlier experiments. The development of the momentum and thermal boundary layers was documented. Velocity and temperature profiles, as well as distributions of skin friction, surface heat transfer rate, Reynolds shear stress, and turbulent heat flux, were shown to compare well with experiment. The results indicate that the essential features of the transition process have been captured. The numerical method used here can be applied to complex geometries in a straightforward manner.

18. Direct Numerical Simulation of Transitional and Turbulent Flow Over a Heated Flat Plate Using Finite-Difference Schemes

NASA Technical Reports Server (NTRS)

1995-01-01

The work in this report was conducted at NASA Ames Research Center during the period from August 1993 to January 1995 deals with the direct numerical simulation of transitional and turbulent flow at low Mach numbers using high-order-accurate finite-difference techniques. A computation of transition to turbulence of the spatially-evolving boundary layer on a heated flat plate in the presence of relatively high freestream turbulence was performed. The geometry and flow conditions were chosen to match earlier experiments. The development of the momentum and thermal boundary layers was documented. Velocity and temperature profiles, as well as distributions of skin friction, surface heat transfer rate, Reynolds shear stress, and turbulent heat flux were shown to compare well with experiment. The numerical method used here can be applied to complex geometries in a straightforward manner.

19. A Low-order Coupled Chemistry Meteorology Model for Testing Online and Offline Advanced Data Assimilation Schemes

Bocquet, M.; Haussaire, J. M.

2015-12-01

Bocquet and Sakov have recently introduced a low-order model based on the coupling of thechaotic Lorenz-95 model which simulates winds along a mid-latitude circle, with thetransport of a tracer species advected by this wind field. It has been used to testadvanced data assimilation methods with an online model that couples meteorology andtracer transport. In the present study, the tracer subsystem of the model is replacedwith a reduced photochemistry module meant to emulate reactive air pollution. Thiscoupled chemistry meteorology model, the L95-GRS model, mimics continental andtranscontinental transport and photochemistry of ozone, volatile organic compounds andnitrogen dioxides.The L95-GRS is specially useful in testing advanced data assimilation schemes, such as theiterative ensemble Kalman smoother (IEnKS) that combines the best of ensemble andvariational methods. The model provides useful insights prior to any implementation ofthe data assimilation method on larger models. For instance, online and offline dataassimilation strategies based on the ensemble Kalman filter or the IEnKS can easily beevaluated with it. It allows to document the impact of species concentration observationson the wind estimation. The model also illustrates a long standing issue in atmosphericchemistry forecasting: the impact of the wind chaotic dynamics and of the chemical speciesnon-chaotic but highly nonlinear dynamics on the selected data assimilation approach.

20. The role of numerical simulation for the development of an advanced HIFU system

Okita, Kohei; Narumi, Ryuta; Azuma, Takashi; Takagi, Shu; Matumoto, Yoichiro

2014-10-01

High-intensity focused ultrasound (HIFU) has been used clinically and is under clinical trials to treat various diseases. An advanced HIFU system employs ultrasound techniques for guidance during HIFU treatment instead of magnetic resonance imaging in current HIFU systems. A HIFU beam imaging for monitoring the HIFU beam and a localized motion imaging for treatment validation of tissue are introduced briefly as the real-time ultrasound monitoring techniques. Numerical simulations have a great impact on the development of real-time ultrasound monitoring as well as the improvement of the safety and efficacy of treatment in advanced HIFU systems. A HIFU simulator was developed to reproduce ultrasound propagation through the body in consideration of the elasticity of tissue, and was validated by comparison with in vitro experiments in which the ultrasound emitted from the phased-array transducer propagates through the acrylic plate acting as a bone phantom. As the result, the defocus and distortion of the ultrasound propagating through the acrylic plate in the simulation quantitatively agree with that in the experimental results. Therefore, the HIFU simulator accurately reproduces the ultrasound propagation through the medium whose shape and physical properties are well known. In addition, it is experimentally confirmed that simulation-assisted focus control of the phased-array transducer enables efficient assignment of the focus to the target. Simulation-assisted focus control can contribute to design of transducers and treatment planning.

1. A new design of the LAPS land surface scheme for use over and through heterogeneous and non-heterogeneous surfaces: Numerical simulations and tests

Mihailovic, Dragutin T.; Lazic, Jelena; Leśny, Jacek; Olejnik, Janusz; Lalic, Branislava; Kapor, Darko; Cirisan, Ana

2010-05-01

Numerical simulations and tests with the recently redesigned land-air parameterization scheme (LAPS) are presented. In all experiments, supported either by one-point micrometeorological, 1D or 3D simulations, the attention has been directed to: (1) comparison of simulation outputs, expressing the energy transfer over and through heterogeneous and non-heterogeneous surfaces, versus observations and (2) analysis of uncertainties occurring in the solution of the energy balance equation at the land-air interface. To check the proposed method for aggregation of albedo, "propagating hole" sensitivity tests with LAPS over a sandstone rock grid cell have been performed with the forcing meteorological data for July 17, 1999 in Baxter site, Philadelphia (USA). Micrometeorological and biophysical measurements from the surface experiments conducted over crops and apple orchard in Serbia, Poland, Austria and France were used to test the operation of LAPS in calculating surface fluxes and canopy environment temperatures within and above plant covers of different densities. In addition, sensitivity tests with single canopy covers over the Central Europe region and comparison against the observations taken from SYNOP data using 3D simulations were made. Validation of LAPS performances over a solid surface has been done by comparison of 2 m air temperature observations against 5-day simulations over the Sahara Desert rocky ground using 3D model. To examine how realistically the LAPS simulates surface processes over a heterogeneous surface, we compared the air temperature measured at 2 m and that predicted by the 1D model with the LAPS as the surface scheme. Finally, the scheme behaviour over urban surface was tested by runs over different parts of a hypothetical urban area. The corresponding 1D simulations were carried out with an imposed meteorological dataset collected during HAPEX-MOBILHY experiment at Caumont (France). The quantities predicted by the LAPS compare well with the

2. Numerical simulation of a dust event in northeastern Germany with a new dust emission scheme in COSMO-ART

Deetz, K.; Klose, M.; Kirchner, I.; Cubasch, U.

2016-02-01

The dust emission scheme of Shao (2004) has been implemented into the regional atmospheric model COSMO-ART and has been applied to a severe dust event in northeastern Germany on 8th April 2011. The model sensitivity to soil moisture and vegetation cover has been studied. Soil moisture has been found to be relatively high in the model during the investigation period and has been reduced by different degree to investigate the resulting changes in dust emissions. Two different vegetation datasets have been tested as model input: the climatological vegetation cover data of COSMO-ART (ECOCLIMAP) and the SPOT5 remote sensing vegetation cover data for the time of the event. By varying soil moisture, vegetation cover and by restricting the potential emission area, a set of eleven simulations was generated. Vegetation cover during the event was about 24% lower on average compared to the climatological mean. Thus, dust emissions modeled with SPOT5 vegetation exceeded that modeled with climatological data by a factor of about 5. The modeled dust concentrations were compared with in-situ measurements of aerosol concentration. The temporal evolutions of simulations and observations have significant correlations (0.42-0.85) especially in rural backgrounds. The lower correlations at urban sites are attributed to anthropogenic PM10 sources, which are not included in the model. However, a verification of the magnitude of modeled dust concentrations is not possible due to the uncertainty in soil moisture and emission area.

3. Numerical Analysis on a Flow Field of Liquid Metals Under a Magnetic Field, Using a Spectral Finite Difference Scheme

Im, Kichang; Mochimaru, Yoshihiro

A steady-state axisymmetric flow field of a liquid metal in a coreless induction furnace under an axisymmetric magnetic field is analyzed numerically, using a spectral finite difference method. Vorticity-stream function formulation is used in conjunction with Maxwell's equations, in a boundary-fitted coordinate system. For boundary conditions, both no-slip on the wall and no shear stress tensor on the free surface are used as dynamic conditions, and a field equivalent to the magnetic field induced by external coils is adopted as an electromagnetic field condition. Presented are streamlines, magnetic streamlines, and radial profiles of the axial velocity component at two Reynolds numbers for various parameters. It is found that the flow field varies remarkably according to the Reynolds number, the dimensionless height of the liquid metal, and the dimensionless height of external coils.

4. A numerical scheme for the identification of hybrid systems describing the vibration of flexible beams with tip bodies

NASA Technical Reports Server (NTRS)

Rosen, I. G.

1984-01-01

A cubic spline based Galerkin-like method is developed for the identification of a class of hybrid systems which describe the transverse vibration to flexible beams with attached tip bodies. The identification problem is formulated as a least squares fit to data subject to the system dynamics given by a coupled system of ordnary and partial differential equations recast as an abstract evolution equation (AEE) in an appropriate infinite dimensional Hilbert space. Projecting the AEE into spline-based subspaces leads naturally to a sequence of approximating finite dimensional identification problems. The solutions to these problems are shown to exist, are relatively easily computed, and are shown to, in some sense, converge to solutions to the original identification problem. Numerical results for a variety of examples are discussed.

5. Impact of a variational objective analysis scheme on a regional area numerical model: The Italian Air Force Weather Service experience

Bonavita, M.; Torrisi, L.

2005-03-01

A new data assimilation system has been designed and implemented at the National Center for Aeronautic Meteorology and Climatology of the Italian Air Force (CNMCA) in order to improve its operational numerical weather prediction capabilities and provide more accurate guidance to operational forecasters. The system, which is undergoing testing before operational use, is based on an “observation space” version of the 3D-VAR method for the objective analysis component, and on the High Resolution Regional Model (HRM) of the Deutscher Wetterdienst (DWD) for the prognostic component. Notable features of the system include a completely parallel (MPI+OMP) implementation of the solution of analysis equations by a preconditioned conjugate gradient descent method; correlation functions in spherical geometry with thermal wind constraint between mass and wind field; derivation of the objective analysis parameters from a statistical analysis of the innovation increments.

6. Numerical Simulations of Optical Turbulence Using an Advanced Atmospheric Prediction Model: Implications for Adaptive Optics Design

Alliss, R.

2014-09-01

Optical turbulence (OT) acts to distort light in the atmosphere, degrading imagery from astronomical telescopes and reducing the data quality of optical imaging and communication links. Some of the degradation due to turbulence can be corrected by adaptive optics. However, the severity of optical turbulence, and thus the amount of correction required, is largely dependent upon the turbulence at the location of interest. Therefore, it is vital to understand the climatology of optical turbulence at such locations. In many cases, it is impractical and expensive to setup instrumentation to characterize the climatology of OT, so numerical simulations become a less expensive and convenient alternative. The strength of OT is characterized by the refractive index structure function Cn2, which in turn is used to calculate atmospheric seeing parameters. While attempts have been made to characterize Cn2 using empirical models, Cn2 can be calculated more directly from Numerical Weather Prediction (NWP) simulations using pressure, temperature, thermal stability, vertical wind shear, turbulent Prandtl number, and turbulence kinetic energy (TKE). In this work we use the Weather Research and Forecast (WRF) NWP model to generate Cn2 climatologies in the planetary boundary layer and free atmosphere, allowing for both point-to-point and ground-to-space seeing estimates of the Fried Coherence length (ro) and other seeing parameters. Simulations are performed using a multi-node linux cluster using the Intel chip architecture. The WRF model is configured to run at 1km horizontal resolution and centered on the Mauna Loa Observatory (MLO) of the Big Island. The vertical resolution varies from 25 meters in the boundary layer to 500 meters in the stratosphere. The model top is 20 km. The Mellor-Yamada-Janjic (MYJ) TKE scheme has been modified to diagnose the turbulent Prandtl number as a function of the Richardson number, following observations by Kondo and others. This modification

7. Development and validation of burnup dependent computational schemes for the analysis of assemblies with advanced lattice codes

Ramamoorthy, Karthikeyan

The main aim of this research is the development and validation of computational schemes for advanced lattice codes. The advanced lattice code which forms the primary part of this research is "DRAGON Version4". The code has unique features like self shielding calculation with capabilities to represent distributed and mutual resonance shielding effects, leakage models with space-dependent isotropic or anisotropic streaming effect, availability of the method of characteristics (MOC), burnup calculation with reaction-detailed energy production etc. Qualified reactor physics codes are essential for the study of all existing and envisaged designs of nuclear reactors. Any new design would require a thorough analysis of all the safety parameters and burnup dependent behaviour. Any reactor physics calculation requires the estimation of neutron fluxes in various regions of the problem domain. The calculation goes through several levels before the desired solution is obtained. Each level of the lattice calculation has its own significance and any compromise at any step will lead to poor final result. The various levels include choice of nuclear data library and energy group boundaries into which the multigroup library is cast; self shielding of nuclear data depending on the heterogeneous geometry and composition; tracking of geometry, keeping error in volume and surface to an acceptable minimum; generation of regionwise and groupwise collision probabilities or MOC-related information and their subsequent normalization thereof, solution of transport equation using the previously generated groupwise information and obtaining the fluxes and reaction rates in various regions of the lattice; depletion of fuel and of other materials based on normalization with constant power or constant flux. Of the above mentioned levels, the present research will mainly focus on two aspects, namely self shielding and depletion. The behaviour of the system is determined by composition of resonant

8. Implementation of Black Sea numerical model based on NEMO and 3DVAR data assimilation scheme for operational forecasting

Ciliberti, Stefania Angela; Peneva, Elisaveta; Storto, Andrea; Rostislav, Kandilarov; Lecci, Rita; Yang, Chunxue; Coppini, Giovanni; Masina, Simona; Pinardi, Nadia

2016-04-01

This study describes a new model implementation for the Black Sea, which uses data assimilation, towards operational forecasting, based on NEMO (Nucleus for European Modelling of the Ocean, Madec et al., 2012). The Black Sea domain is resolved with 1/27°×1/36° horizontal resolution (~3 km) and 31 z-levels with partial steps based on the GEBCO bathymetry data (Grayek et al., 2010). The model is forced by momentum, water and heat fluxes interactively computed by bulk formulae using high resolution atmospheric forcing provided by the European Centre for Medium-Range Forecast (ECMWF). The initial condition is calculated from long-term climatological temperature and salinity 3D fields. Precipitation field over the basin has been computed from the climatological GPCP rainfall monthly data (Adler et al., 2003; Huffman et al., 2009), while the evaporation is derived from the latent heat flux. The climatological monthly mean runoff of the major rivers in the Black Sea is computed using the hydrological dataset provided by SESAME project (Ludvig et al., 2009). The exchange with Mediterranean Sea through the Bosporus Straits is represented by a surface boundary condition taking into account the barotropic transport calculated to balance the fresh water fluxes on monthly bases (Stanev and Beckers, 1999, Peneva et al., 2001). A multi-annual run 2011-2015 has been completed in order to describe the main characteristics of the Black Sea circulation dynamics and thermohaline structure and the numerical results have been validated using in-situ (ARGO) and satellite (SST, SLA) data. The Black Sea model represents also the core of the new Black Sea Forecasting System, implemented at CMCC operationally since January 2016, which produces at daily frequency 10-day forecasts, 3-days analyses and 1-day simulation. Once a week, the system is run 15-day in the past in analysis mode to compute the new optimal initial condition for the forecast cycle. The assimilation is performed by a

9. Numerical Evaluation of Fluid Mixing Phenomena in Boiling Water Reactor Using Advanced Interface Tracking Method

Yoshida, Hiroyuki; Takase, Kazuyuki

Thermal-hydraulic design of the current boiling water reactor (BWR) is performed with the subchannel analysis codes which incorporated the correlations based on empirical results including actual-size tests. Then, for the Innovative Water Reactor for Flexible Fuel Cycle (FLWR) core, an actual size test of an embodiment of its design is required to confirm or modify such correlations. In this situation, development of a method that enables the thermal-hydraulic design of nuclear reactors without these actual size tests is desired, because these tests take a long time and entail great cost. For this reason, we developed an advanced thermal-hydraulic design method for FLWRs using innovative two-phase flow simulation technology. In this study, a detailed Two-Phase Flow simulation code using advanced Interface Tracking method: TPFIT is developed to calculate the detailed information of the two-phase flow. In this paper, firstly, we tried to verify the TPFIT code by comparing it with the existing 2-channel air-water mixing experimental results. Secondary, the TPFIT code was applied to simulation of steam-water two-phase flow in a model of two subchannels of a current BWRs and FLWRs rod bundle. The fluid mixing was observed at a gap between the subchannels. The existing two-phase flow correlation for fluid mixing is evaluated using detailed numerical simulation data. This data indicates that pressure difference between fluid channels is responsible for the fluid mixing, and thus the effects of the time average pressure difference and fluctuations must be incorporated in the two-phase flow correlation for fluid mixing. When inlet quality ratio of subchannels is relatively large, it is understood that evaluation precision of the existing two-phase flow correlations for fluid mixing are relatively low.

10. Investigation of upwind, multigrid, multiblock numerical schemes for three dimensional flows. Volume 1: Runge-Kutta methods for a thin layer Navier-Stokes solver

NASA Technical Reports Server (NTRS)

Cannizzaro, Frank E.; Ash, Robert L.

1992-01-01

A state-of-the-art computer code has been developed that incorporates a modified Runge-Kutta time integration scheme, upwind numerical techniques, multigrid acceleration, and multi-block capabilities (RUMM). A three-dimensional thin-layer formulation of the Navier-Stokes equations is employed. For turbulent flow cases, the Baldwin-Lomax algebraic turbulence model is used. Two different upwind techniques are available: van Leer's flux-vector splitting and Roe's flux-difference splitting. Full approximation multi-grid plus implicit residual and corrector smoothing were implemented to enhance the rate of convergence. Multi-block capabilities were developed to provide geometric flexibility. This feature allows the developed computer code to accommodate any grid topology or grid configuration with multiple topologies. The results shown in this dissertation were chosen to validate the computer code and display its geometric flexibility, which is provided by the multi-block structure.

11. Pollutant transport by shallow water equations on unstructured meshes: Hyperbolization of the model and numerical solution via a novel flux splitting scheme

Vanzo, Davide; Siviglia, Annunziato; Toro, Eleuterio F.

2016-09-01

The purpose of this paper is twofold. First, using the Cattaneo's relaxation approach, we reformulate the system of governing equations for the pollutant transport by shallow water flows over non-flat topography and anisotropic diffusion as hyperbolic balance laws with stiff source terms. The proposed relaxation system circumvents the infinite wave speed paradox which is inherent in standard advection-diffusion models. This turns out to give a larger stability range for the choice of the time step. Second, following a flux splitting approach, we derive a novel numerical method to discretise the resulting problem. In particular, we propose a new flux splitting and study the associated two systems of differential equations, called the "hydrodynamic" and the "relaxed diffusive" system, respectively. For the presented splitting we analyse the resulting two systems of differential equations and propose two discretisation schemes of the Godunov-type. These schemes are simple to implement, robust, accurate and fast when compared with existing methods. The resulting method is implemented on unstructured meshes and is systematically assessed for accuracy, robustness and efficiency on a carefully selected suite of test problems including non-flat topography and wetting and drying problems. Formal second order accuracy is assessed through convergence rates studies.

12. Advanced Numerical Imaging Procedure Accounting for Non-Ideal Effects in GPR Scenarios

Comite, Davide; Galli, Alessandro; Catapano, Ilaria; Soldovieri, Francesco

2015-04-01

advanced implementation have also been tested by introducing 'errors' on the knowledge of the background medium permittivity, by simulating the presence of one or more layers, and by choosing different models of the surface roughness. The impact of these issues on the performance of both the conventional procedure and the advanced one will be extensively highlighted and discussed at the conference. [1] G. Valerio et al., "GPR detectability of rocks in a Martian-like shallow subsoil: A numerical approach," Plan. Sp. Sci., vol. 62, pp. 31-40, 2012. [2] A. Galli et al., "3D imaging of buried dielectric targets with a tomographic microwave approach applied to GPR synthetic data," Int. J. Antennas Propag., art. ID 610389, 10 pp., 2013 [3] F. Soldovieri et al., "A linear inverse scattering algorithm for realistic GPR applications," Near Surface Geophysics, 5 (1), pp. 29-42, 2007.

13. Recent advances in theoretical and numerical studies of wire array Z-pinch in the IAPCM

SciTech Connect

Ding, Ning Zhang, Yang Xiao, Delong Wu, Jiming Huang, Jun Yin, Li Sun, Shunkai Xue, Chuang Dai, Zihuan Ning, Cheng Shu, Xiaojian Wang, Jianguo Li, Hua

2014-12-15

Fast Z-pinch has produced the most powerful X-ray radiation source in laboratory and also shows the possibility to drive inertial confinement fusion (ICF). Recent advances in wire-array Z-pinch researches at the Institute of Applied Physics and Computational Mathematics are presented in this paper. A typical wire array Z-pinch process has three phases: wire plasma formation and ablation, implosion and the MRT instability development, stagnation and radiation. A mass injection model with azimuthal modulation coefficient is used to describe the wire initiation, and the dynamics of ablated plasmas of wire-array Z-pinches in (r, θ) geometry is numerically studied. In the implosion phase, a two-dimensional(r, z) three temperature radiation MHD code MARED has been developed to investigate the development of the Magneto-Rayleigh-Taylor(MRT) instability. We also analyze the implosion modes of nested wire-array and find that the inner wire-array is hardly affected before the impaction of the outer wire-array. While the plasma accelerated to high speed in the implosion stage stagnates on the axis, abundant x-ray radiation is produced. The energy spectrum of the radiation and the production mechanism are investigated. The computational x-ray pulse shows a reasonable agreement with the experimental result. We also suggest that using alloyed wire-arrays can increase multi-keV K-shell yield by decreasing the opacity of K-shell lines. In addition, we use a detailed circuit model to study the energy coupling between the generator and the Z-pinch implosion. Recently, we are concentrating on the problems of Z-pinch driven ICF, such as dynamic hohlraum and capsule implosions. Our numerical investigations on the interaction of wire-array Z-pinches on foam convertors show qualitative agreements with experimental results on the “Qiangguang I” facility. An integrated two-dimensional simulation of dynamic hohlraum driven capsule implosion provides us the physical insights of wire

14. Recent advances in theoretical and numerical studies of wire array Z-pinch in the IAPCM

Ding, Ning; Zhang, Yang; Xiao, Delong; Wu, Jiming; Huang, Jun; Yin, Li; Sun, Shunkai; Xue, Chuang; Dai, Zihuan; Ning, Cheng; Shu, Xiaojian; Wang, Jianguo; Li, Hua

2014-12-01

Fast Z-pinch has produced the most powerful X-ray radiation source in laboratory and also shows the possibility to drive inertial confinement fusion (ICF). Recent advances in wire-array Z-pinch researches at the Institute of Applied Physics and Computational Mathematics are presented in this paper. A typical wire array Z-pinch process has three phases: wire plasma formation and ablation, implosion and the MRT instability development, stagnation and radiation. A mass injection model with azimuthal modulation coefficient is used to describe the wire initiation, and the dynamics of ablated plasmas of wire-array Z-pinches in (r, θ) geometry is numerically studied. In the implosion phase, a two-dimensional(r, z) three temperature radiation MHD code MARED has been developed to investigate the development of the Magneto-Rayleigh-Taylor(MRT) instability. We also analyze the implosion modes of nested wire-array and find that the inner wire-array is hardly affected before the impaction of the outer wire-array. While the plasma accelerated to high speed in the implosion stage stagnates on the axis, abundant x-ray radiation is produced. The energy spectrum of the radiation and the production mechanism are investigated. The computational x-ray pulse shows a reasonable agreement with the experimental result. We also suggest that using alloyed wire-arrays can increase multi-keV K-shell yield by decreasing the opacity of K-shell lines. In addition, we use a detailed circuit model to study the energy coupling between the generator and the Z-pinch implosion. Recently, we are concentrating on the problems of Z-pinch driven ICF, such as dynamic hohlraum and capsule implosions. Our numerical investigations on the interaction of wire-array Z-pinches on foam convertors show qualitative agreements with experimental results on the "Qiangguang I" facility. An integrated two-dimensional simulation of dynamic hohlraum driven capsule implosion provides us the physical insights of wire

15. An Advanced Temporal Credential-Based Security Scheme with Mutual Authentication and Key Agreement for Wireless Sensor Networks

PubMed Central

Li, Chun-Ta; Weng, Chi-Yao; Lee, Cheng-Chi

2013-01-01

Wireless sensor networks (WSNs) can be quickly and randomly deployed in any harsh and unattended environment and only authorized users are allowed to access reliable sensor nodes in WSNs with the aid of gateways (GWNs). Secure authentication models among the users, the sensor nodes and GWN are important research issues for ensuring communication security and data privacy in WSNs. In 2013, Xue et al. proposed a temporal-credential-based mutual authentication and key agreement scheme for WSNs. However, in this paper, we point out that Xue et al.'s scheme cannot resist stolen-verifier, insider, off-line password guessing, smart card lost problem and many logged-in users' attacks and these security weaknesses make the scheme inapplicable to practical WSN applications. To tackle these problems, we suggest a simple countermeasure to prevent proposed attacks while the other merits of Xue et al.'s authentication scheme are left unchanged. PMID:23887085

16. An advanced temporal credential-based security scheme with mutual authentication and key agreement for wireless sensor networks.

PubMed

Li, Chun-Ta; Weng, Chi-Yao; Lee, Cheng-Chi

2013-01-01

Wireless sensor networks (WSNs) can be quickly and randomly deployed in any harsh and unattended environment and only authorized users are allowed to access reliable sensor nodes in WSNs with the aid of gateways (GWNs). Secure authentication models among the users, the sensor nodes and GWN are important research issues for ensuring communication security and data privacy in WSNs. In 2013, Xue et al. proposed a temporal-credential-based mutual authentication and key agreement scheme for WSNs. However, in this paper, we point out that Xue et al.'s scheme cannot resist stolen-verifier, insider, off-line password guessing, smart card lost problem and many logged-in users' attacks and these security weaknesses make the scheme inapplicable to practical WSN applications. To tackle these problems, we suggest a simple countermeasure to prevent proposed attacks while the other merits of Xue et al.'s authentication scheme are left unchanged. PMID:23887085

17. Development and validation of burnup dependent computational schemes for the analysis of assemblies with advanced lattice codes

Ramamoorthy, Karthikeyan

The main aim of this research is the development and validation of computational schemes for advanced lattice codes. The advanced lattice code which forms the primary part of this research is "DRAGON Version4". The code has unique features like self shielding calculation with capabilities to represent distributed and mutual resonance shielding effects, leakage models with space-dependent isotropic or anisotropic streaming effect, availability of the method of characteristics (MOC), burnup calculation with reaction-detailed energy production etc. Qualified reactor physics codes are essential for the study of all existing and envisaged designs of nuclear reactors. Any new design would require a thorough analysis of all the safety parameters and burnup dependent behaviour. Any reactor physics calculation requires the estimation of neutron fluxes in various regions of the problem domain. The calculation goes through several levels before the desired solution is obtained. Each level of the lattice calculation has its own significance and any compromise at any step will lead to poor final result. The various levels include choice of nuclear data library and energy group boundaries into which the multigroup library is cast; self shielding of nuclear data depending on the heterogeneous geometry and composition; tracking of geometry, keeping error in volume and surface to an acceptable minimum; generation of regionwise and groupwise collision probabilities or MOC-related information and their subsequent normalization thereof, solution of transport equation using the previously generated groupwise information and obtaining the fluxes and reaction rates in various regions of the lattice; depletion of fuel and of other materials based on normalization with constant power or constant flux. Of the above mentioned levels, the present research will mainly focus on two aspects, namely self shielding and depletion. The behaviour of the system is determined by composition of resonant

18. Advances in Analytical and Numerical Dispersion Modeling of Pollutants Releasing from an Area-source

Nimmatoori, Praneeth

The air quality near agricultural activities such as tilling, plowing, harvesting, and manure application is of main concern because they release fine particulate matter into the atmosphere. These releases are modeled as area-sources in the air quality modeling research. None of the currently available dispersion models relate and incorporate physical characteristics and meteorological conditions for modeling the dispersion and deposition of particulates emitting from such area-sources. This knowledge gap was addressed by developing the advanced analytical and numerical methods for modeling the dispersion of particulate matter. The development, application, and evaluation of new dispersion modeling methods are discussed in detail in this dissertation. In the analytical modeling, a ground-level area source analytical dispersion model known as particulate matter deposition -- PMD was developed for predicting the concentrations of different particle sizes. Both the particle dynamics (particle physical characteristics) and meteorological conditions which have significant effect on the dispersion of particulates were related and incorporated in the PMD model using the formulations of particle gravitational settling and dry deposition velocities. The modeled particle size concentrations of the PMD model were evaluated statistically after applying it to particulates released from a biosolid applied agricultural field. The evaluation of the PMD model using the statistical criteria concluded effective and successful inclusion of dry deposition theory for modeling particulate matter concentrations. A comprehensive review of analytical area-source dispersion models, which do not account for dry deposition and treat pollutants as gases, was conducted and determined three models -- the Shear, the Parker, and the Smith. A statistical evaluation of these dispersion models was conducted after applying them to two different field data sets and the statistical results concluded that

19. Advanced material modelling in numerical simulation of primary acetabular press-fit cup stability.

PubMed

Souffrant, R; Zietz, C; Fritsche, A; Kluess, D; Mittelmeier, W; Bader, R

2012-01-01

Primary stability of artificial acetabular cups, used for total hip arthroplasty, is required for the subsequent osteointegration and good long-term clinical results of the implant. Although closed-cell polymer foams represent an adequate bone substitute in experimental studies investigating primary stability, correct numerical modelling of this material depends on the parameter selection. Material parameters necessary for crushable foam plasticity behaviour were originated from numerical simulations matched with experimental tests of the polymethacrylimide raw material. Experimental primary stability tests of acetabular press-fit cups consisting of static shell assembly with consecutively pull-out and lever-out testing were subsequently simulated using finite element analysis. Identified and optimised parameters allowed the accurate numerical reproduction of the raw material tests. Correlation between experimental tests and the numerical simulation of primary implant stability depended on the value of interference fit. However, the validated material model provides the opportunity for subsequent parametric numerical studies. PMID:22817471

ERIC Educational Resources Information Center

Iler, H. Darrell; Brown, Amber; Landis, Amanda; Schimke, Greg; Peters, George

2014-01-01

A numerical analysis of the free radical addition polymerization system is described that provides those teaching polymer, physical, or advanced organic chemistry courses the opportunity to introduce students to numerical methods in the context of a simple but mathematically stiff chemical kinetic system. Numerical analysis can lead students to an…

1. Springback Simulation: Impact of Some Advanced Constitutive Models and Numerical Parameters

2005-08-01

The impact of material models on the numerical simulation of springback is investigated. The study is focused on the strain-path sensitivity of two hardening models. While both models predict the Bauschinger effect, their response in the transient zone after a strain-path change is fairly different. Their respective predictions are compared in terms of sequential test response and of strip-drawing springback. For this purpose, an accurate and general time integration algorithm has been developed and implemented in the Abaqus code. The impact of several numerical parameters is also studied in order to assess the overall accuracy of the finite element prediction. For some test geometries, both material and numerical parameters are shown to clearly influence the springback behavior at a large extent. Moreover, a general trend cannot always be extracted, thus justifying the need for the finite element simulation of the stamping process.

2. A review of recent advances in numerical simulations of microscale fuel processor for hydrogen production

2015-05-01

Microscale (<5 W) reformers for hydrogen production have been investigated for over a decade. These devices are intended to provide hydrogen for small fuel cells. Due to the reformer's small size, numerical simulations are critical to understand heat and mass transfer phenomena occurring in the systems and help guide the further improvements. This paper reviews the development of the numerical codes and details the reaction equations used. The majority of the devices utilized methanol as the fuel due to methanol's low reforming temperature and high conversion, although, there are several methane fueled systems. The increased computational power and more complex codes have led to improved accuracy of numerical simulations. Initial models focused on the reformer, while more recently, the simulations began including other unit operations such as vaporizers, inlet manifolds, and combustors. These codes are critical for developing the next generation systems. The systems reviewed included plate reactors, microchannel reactors, and annulus reactors for both wash-coated and packed bed systems.

3. CNC Turning Center Advanced Operations. Computer Numerical Control Operator/Programmer. 444-332.

ERIC Educational Resources Information Center

Skowronski, Steven D.; Tatum, Kenneth

This student guide provides materials for a course designed to introduce the student to the operations and functions of a two-axis computer numerical control (CNC) turning center. The course consists of seven units. Unit 1 presents course expectations and syllabus, covers safety precautions, and describes the CNC turning center components, CNC…

4. A numerical technique for calculation of the noise of high-speed propellers with advanced blade geometry

NASA Technical Reports Server (NTRS)

Nystrom, P. A.; Farassat, F.

1980-01-01

A numerical technique and computer program were developed for the prediction of the noise of propellers with advanced geometry. The blade upper and lower surfaces are described by a curvilinear coordinate system, which was also used to divide the blade surfaces into panels. Two different acoustic formulations in the time domain were used to improve the speed and efficiency of the noise calculations: an acoustic formualtion with the Doppler factor singularity for panels moving at subsonic speeds and the collapsing sphere formulation for panels moving at transonic or supersonic speeds. This second formulation involves a sphere which is centered at the observer position and whose radius decreases at the speed of sound. The acoustic equation consisted of integrals over the curve of intersection for both the sphere and the panels on the blade. Algorithms used in some parts of the computer program are discussed. Comparisons with measured acoustic data for two model high speed propellers with advanced geometry are also presented.

5. Numerical assessment of radiation binary targeted therapy for HER-2 positive breast cancers: advanced calculations and radiation dosimetry

Sztejnberg Gonçalves-Carralves, Manuel L.; Jevremovic, Tatjana

2007-07-01

In our previous publication (Mundy et al 2006 Phys. Med. Biol. 51 1377) we have described the theoretical assessment of our novel approach in radiation binary targeted therapy for HER-2 positive breast cancers and summarized the future directions in this area of research. In this paper we advanced the numerical analysis to show the detailed radiation dose distribution for various neutron sources in combination with the required boron concentration and allowed radiation skin doses. We once again proved the feasibility of the concept and will use these data and conclusions to start with the experimental verifications.

6. A review of recent advances of numerical simulations of microscale fuel processors for hydrogen production

SciTech Connect

2015-05-01

Microscale (<5W) reformers for hydrogen production have been investigated for over a decade. These devices are intended to provide hydrogen for small fuel cells. Due to the reformer’s small size, numerical simulations are critical to understand heat and mass transfer phenomena occurring in the systems. This paper reviews the development of the numerical codes and details the reaction equations used. The majority of the devices utilized methanol as the fuel due to methanol’s low reforming temperature and high conversion, although, there are several methane fueled systems. As computational power has decreased in cost and increased in availability, the codes increased in complexity and accuracy. Initial models focused on the reformer, while more recently, the simulations began including other unit operations such as vaporizers, inlet manifolds, and combustors. These codes are critical for developing the next generation systems. The systems reviewed included, plate reactors, microchannel reactors, annulus reactors, wash-coated, packed bed systems.

7. Advanced Techniques for Seismic Protection of Historical Buildings: Experimental and Numerical Approach

SciTech Connect

Mazzolani, Federico M.

2008-07-08

The seismic protection of historical and monumental buildings, namely dating back from the ancient age up to the 20th Century, is being looked at with greater and greater interest, above all in the Euro-Mediterranean area, its cultural heritage being strongly susceptible to undergo severe damage or even collapse due to earthquake. The cultural importance of historical and monumental constructions limits, in many cases, the possibility to upgrade them from the seismic point of view, due to the fear of using intervention techniques which could have detrimental effects on their cultural value. Consequently, a great interest is growing in the development of sustainable methodologies for the use of Reversible Mixed Technologies (RMTs) in the seismic protection of the existing constructions. RMTs, in fact, are conceived for exploiting the peculiarities of innovative materials and special devices, and they allow ease of removal when necessary. This paper deals with the experimental and numerical studies, framed within the EC PROHITECH research project, on the application of RMTs to the historical and monumental constructions mainly belonging to the cultural heritage of the Euro-Mediterranean area. The experimental tests and the numerical analyses are carried out at five different levels, namely full scale models, large scale models, sub-systems, devices, materials and elements.

8. Numerical Study on Crossflow Printed Circuit Heat Exchanger for Advanced Small Modular Reactors

SciTech Connect

Yoon, Su-Jong; Sabharwall, Piyush; Kim, Eung-Soo

2014-03-01

Various fluids such as water, gases (helium), molten salts (FLiNaK, FLiBe) and liquid metal (sodium) are used as a coolant of advanced small modular reactors (SMRs). The printed circuit heat exchanger (PCHE) has been adopted as the intermediate and/or secondary heat exchanger of SMR systems because this heat exchanger is compact and effective. The size and cost of PCHE can be changed by the coolant type of each SMR. In this study, the crossflow PCHE analysis code for advanced small modular reactor has been developed for the thermal design and cost estimation of the heat exchanger. The analytical solution of single pass, both unmixed fluids crossflow heat exchanger model was employed to calculate a two dimensional temperature profile of a crossflow PCHE. The analytical solution of crossflow heat exchanger was simply implemented by using built in function of the MATLAB program. The effect of fluid property uncertainty on the calculation results was evaluated. In addition, the effect of heat transfer correlations on the calculated temperature profile was analyzed by taking into account possible combinations of primary and secondary coolants in the SMR systems. Size and cost of heat exchanger were evaluated for the given temperature requirement of each SMR.

9. Numerical simulation and intercomparison of boundary layer structure with different PBL schemes in WRF using experimental observations at a tropical site

Hariprasad, K. B. R. R.; Srinivas, C. V.; Singh, A. Bagavath; Vijaya Bhaskara Rao, S.; Baskaran, R.; Venkatraman, B.

2014-08-01

In this study the performance of seven PBL parameterizations in the Weather Research and Forecast (WRF-ARW) mesoscale model was tested at the tropical site Kalpakkam. Meteorological observations collected during an intense observation campaign for wind field modeling called Round Robin Exercise (RRE) were used for comparison. High resolution simulations were conducted for a warm summer condition on 22-24 September 2010. The observations included GPS Sonde vertical profiles, surface level data from meteorological towers and turbulent fluxes from sonic anemometers. Sensitivity experiments with seven PBL schemes [Mellor-Yamada-Janjic (MYJ), Mellor-Yamada-Nakanishi-Niino (MYNN), Quasi Normal Scale Elimination (QNSE), Yonsei University (YSU), Asymmetric Convective Model (ACM2), Bougeault-Lacarrére (BL), Bretherton-Park (UW)] indicated that while all the schemes similarly produced the stable boundary layer characteristics there were large differences in the convective daytime PBL. It has been found that while ACM2 and QNSE produced highly unstable and deep convective layers, the UW produced relatively shallow mixed layer and all other schemes (YSU, MYNN, MYJ, BL) produced intermediately deep convective layers. All the schemes well produced the vertical wind directional shear within the PBL. A wide variation in the eddy diffusivities was simulated by different PBL schemes in convective daytime condition. ACM2 and UW produced excessive diffusivities which led to relatively weaker winds, warmer and dryer mixed layers with these schemes. Overall the schemes MYNN and YSU simulated the various PBL quantities in better agreement with observations. The differences in the simulated PBL structures could be partly due to various surface layer formulations that produced variation in friction velocity and heat fluxes in each case.

10. Numerical modelling of the groundwater inflow to an advancing open pit mine: Kolahdarvazeh pit, Central Iran.

PubMed

Bahrami, Saeed; Doulati Ardejani, Faramarz; Aslani, Soheyla; Baafi, Ernest

2014-12-01

The groundwater inflow into a mine during its life and after ceasing operations is one of the most important concerns of the mining industry. This paper presents a hydrogeological assessment of the Irankuh Zn-Pb mine at 20 km south of Esfahan and 1 km northeast of Abnil in west-Central Iran. During mine excavation, the upper impervious bed of a confined aquifer was broken and water at high-pressure flowed into an open pit mine associated with the Kolahdarvazeh deposit. The inflow rates were 6.7 and 1.4 m(3)/s at the maximum and minimum quantities, respectively. Permeability, storage coefficient, thickness and initial head of the fully saturated confined aquifer were 3.5 × 10(-4) m/s, 0.2, 30 m and 60 m, respectively. The hydraulic heads as a function of time were monitored at four observation wells in the vicinity of the pit over 19 weeks and at an observation well near a test well over 21 h. In addition, by measuring the rate of pumping out from the pit sump, at a constant head (usually equal to height of the pit floor), the real inflow rates to the pit were monitored. The main innovations of this work were to make comparison between numerical modelling using a finite element software called SEEP/W and actual data related to inflow and extend the applicability of the numerical model. This model was further used to estimate the hydraulic heads at the observation wells around the pit over 19 weeks during mining operations. Data from a pump-out test and observation wells were used for model calibration and verification. In order to evaluate the model efficiency, the modelling results of inflow quantity and hydraulic heads were compared to those from analytical solutions, as well as the field data. The mean percent error in relation to field data for the inflow quantity was 0.108. It varied between 1.16 and 1.46 for hydraulic head predictions, which are much lower values than the mean percent errors resulted from the analytical solutions (from 1.8 to 5