Advanced numerical methods in mesh generation and mesh adaptation
Lipnikov, Konstantine; Danilov, A; Vassilevski, Y; Agonzal, A
2010-01-01
Numerical solution of partial differential equations requires appropriate meshes, efficient solvers and robust and reliable error estimates. Generation of high-quality meshes for complex engineering models is a non-trivial task. This task is made more difficult when the mesh has to be adapted to a problem solution. This article is focused on a synergistic approach to the mesh generation and mesh adaptation, where best properties of various mesh generation methods are combined to build efficiently simplicial meshes. First, the advancing front technique (AFT) is combined with the incremental Delaunay triangulation (DT) to build an initial mesh. Second, the metric-based mesh adaptation (MBA) method is employed to improve quality of the generated mesh and/or to adapt it to a problem solution. We demonstrate with numerical experiments that combination of all three methods is required for robust meshing of complex engineering models. The key to successful mesh generation is the high-quality of the triangles in the initial front. We use a black-box technique to improve surface meshes exported from an unattainable CAD system. The initial surface mesh is refined into a shape-regular triangulation which approximates the boundary with the same accuracy as the CAD mesh. The DT method adds robustness to the AFT. The resulting mesh is topologically correct but may contain a few slivers. The MBA uses seven local operations to modify the mesh topology. It improves significantly the mesh quality. The MBA method is also used to adapt the mesh to a problem solution to minimize computational resources required for solving the problem. The MBA has a solid theoretical background. In the first two experiments, we consider the convection-diffusion and elasticity problems. We demonstrate the optimal reduction rate of the discretization error on a sequence of adaptive strongly anisotropic meshes. The key element of the MBA method is construction of a tensor metric from hierarchical edge
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.
Advanced numerical methods and software approaches for semiconductor device simulation
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.
Advanced Numerical Methods and Software Approaches for Semiconductor Device Simulation
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
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)
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
NASA Astrophysics Data System (ADS)
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
NASA Astrophysics Data System (ADS)
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.
Advanced numerical methods for three dimensional two-phase flow calculations
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.
Recent advances in methods for numerical solution of O.D.E. initial value problems
NASA Technical Reports Server (NTRS)
Bui, T. D.; Oppenheim, A. K.; Pratt, D. T.
1984-01-01
In the mathematical modeling of physical systems, it is often necessary to solve an initial value problem (IVP), consisting of a system of ordinary differential equations (ODE). A typical program produces approximate solutions at certain mesh points. Almost all existing codes try to control the local truncation error, while the user is really interested in controlling the true or global error. The present investigation provides a review of recent advances regarding the solution of the IVP, giving particular attention to stiff systems. Stiff phenomena are customarily defined in terms of the eigenvalues of the Jacobian. There are, however, some difficulties connected with this approach. It is pointed out that an estimate of the Lipschitz constant proves to be a very practical way to determine the stiffness of a problem.
NASA Technical Reports Server (NTRS)
Atwater, James; Wheeler, Richard, Jr.; Akse, James; Jovanovic, Goran; Reed, Brian
2013-01-01
To support long-duration manned missions in space such as a permanent lunar base, Mars transit, or Mars Surface Mission, improved methods for the treatment of solid wastes, particularly methods that recover valuable resources, are needed. The ability to operate under microgravity and hypogravity conditions is essential to meet this objective. The utilization of magnetic forces to manipulate granular magnetic media has provided the means to treat solid wastes under variable gravity conditions by filtration using a consolidated magnetic media bed followed by thermal processing of the solid wastes in a fluidized bed reactor. Non-uniform magnetic fields will produce a magnetic field gradient in a bed of magnetically susceptible media toward the distributor plate of a fluidized bed reactor. A correctly oriented magnetic field gradient will generate a downward direct force on magnetic media that can substitute for gravitational force in microgravity, or which may augment low levels of gravity, such as on the Moon or Mars. This approach is termed Gradient Magnetically Assisted Fluidization (G-MAFB), in which the magnitude of the force on the fluidized media depends upon the intensity of the magnetic field (H), the intensity of the field gradient (dH/dz), and the magnetic susceptibility of the media. Fluidized beds based on the G-MAFB process can operate in any gravitational environment by tuning the magnetic field appropriately. Magnetic materials and methods have been developed that enable G-MAFB operation under variable gravity conditions.
NASA Astrophysics Data System (ADS)
Blum, Volker
This talk describes recent advances of a general, efficient, accurate all-electron electronic theory approach based on numeric atom-centered orbitals; emphasis is placed on developments related to materials for energy conversion and their discovery. For total energies and electron band structures, we show that the overall accuracy is on par with the best benchmark quality codes for materials, but scalable to large system sizes (1,000s of atoms) and amenable to both periodic and non-periodic simulations. A recent localized resolution-of-identity approach for the Coulomb operator enables O (N) hybrid functional based descriptions of the electronic structure of non-periodic and periodic systems, shown for supercell sizes up to 1,000 atoms; the same approach yields accurate results for many-body perturbation theory as well. For molecular systems, we also show how many-body perturbation theory for charged and neutral quasiparticle excitation energies can be efficiently yet accurately applied using basis sets of computationally manageable size. Finally, the talk highlights applications to the electronic structure of hybrid organic-inorganic perovskite materials, as well as to graphene-based substrates for possible future transition metal compound based electrocatalyst materials. All methods described here are part of the FHI-aims code. VB gratefully acknowledges contributions by numerous collaborators at Duke University, Fritz Haber Institute Berlin, TU Munich, USTC Hefei, Aalto University, and many others around the globe.
NASA Astrophysics Data System (ADS)
Zhang, Weizhong; Yoshida, Hiroyuki; Ose, Yasuo; Ohnuki, Akira; Akimoto, Hajime; Hotta, Akitoshi; Fujimura, Ken
In relation to the design of an innovative FLexible-fuel-cycle Water Reactor (FLWR), investigation of thermal-hydraulic performance in tight-lattice rod bundles of the FLWR is being carried out at Japan Atomic Energy Agency (JAEA). The FLWR core adopts a tight triangular lattice arrangement with about 1 mm gap clearance between adjacent fuel rods. In view of importance of accurate prediction of cross flow between subchannels in the evaluation of the boiling transition (BT) in the FLWR core, this study presents a statistical evaluation of numerical simulation results obtained by a detailed two-phase flow simulation code, TPFIT, which employs an advanced interface tracking method. In order to clarify mechanisms of cross flow in such tight lattice rod bundles, the TPFIT is applied to simulate water-steam two-phase flow in two modeled subchannels. Attention is focused on instantaneous fluctuation characteristics of cross flow. With the calculation of correlation coefficients between differential pressure and gas/liquid mixing coefficients, time scales of cross flow are evaluated, and effects of mixing section length, flow pattern and gap spacing on correlation coefficients are investigated. Differences in mechanism between gas and liquid cross flows are pointed out.
Numerical methods for turbulent flow
NASA Astrophysics Data System (ADS)
Turner, James C., Jr.
1988-09-01
It has generally become accepted that the Navier-Strokes equations predict the dynamic behavior of turbulent as well as laminar flows of a fluid at a point in space away form a discontinuity such as a shock wave. Turbulence is also closely related to the phenomena of non-uniqueness of solutions of the Navier-Strokes equations. These second order, nonlinear partial differential equations can be solved analytically for only a few simple flows. Turbulent flow fields are much to complex to lend themselves to these few analytical methods. Numerical methods, therefore, offer the only possibility of achieving a solution of turbulent flow equations. In spite of recent advances in computer technology, the direct solution, by discrete methods, of the Navier-Strokes equations for turbulent flow fields is today, and in the foreseeable future, impossible. Thus the only economically feasible way to solve practical turbulent flow problems numerically is to use statistically averaged equations governing mean-flow quantities. The objective is to study some recent developments relating to the use of numerical methods to study turbulent flow.
NASA Astrophysics Data System (ADS)
Wang, Ying; Krafczyk, Manfred; Geier, Martin; Schönherr, Martin
2014-05-01
The quantification of soil evaporation and of soil water content dynamics near the soil surface are critical in the physics of land-surface processes on many scales and are dominated by multi-component and multi-phase mass and energy fluxes between the ground and the atmosphere. Although it is widely recognized that both liquid and gaseous water movement are fundamental factors in the quantification of soil heat flux and surface evaporation, their computation has only started to be taken into account using simplified macroscopic models. As the flow field over the soil can be safely considered as turbulent, it would be natural to study the detailed transient flow dynamics by means of Large Eddy Simulation (LES [1]) where the three-dimensional flow field is resolved down to the laminar sub-layer. Yet this requires very fine resolved meshes allowing a grid resolution of at least one order of magnitude below the typical grain diameter of the soil under consideration. In order to gain reliable turbulence statistics, up to several hundred eddy turnover times have to be simulated which adds up to several seconds of real time. Yet, the time scale of the receding saturated water front dynamics in the soil is on the order of hours. Thus we are faced with the task of solving a transient turbulent flow problem including the advection-diffusion of water vapour over the soil-atmospheric interface represented by a realistic tomographic reconstruction of a real porous medium taken from laboratory probes. Our flow solver is based on the Lattice Boltzmann method (LBM) [2] which has been extended by a Cumulant approach similar to the one described in [3,4] to minimize the spurious coupling between the degrees of freedom in previous LBM approaches and can be used as an implicit LES turbulence model due to its low numerical dissipation and increased stability at high Reynolds numbers. The kernel has been integrated into the research code Virtualfluids [5] and delivers up to 30% of the
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.
Numerical methods in structural mechanics
NASA Astrophysics Data System (ADS)
Obraztsov, I. F.
The papers contained in this volume focus on numerical, numerical-analytical, and theoretical methods for dealing with strength, stability, and dynamics problems in the design of the structural elements of flight vehicles. Topics discussed include the solution of homogeneous boundary value problems for systems of ordinary differential equations modified by a difference factorization method, a study of the rupture strength of a welded joint between plates, singular solutions in mixed problems for a wedge and a half-strip, and a thermoelasticity problem for an open-profile cylindrical shell with a localized temperature field.
Brush seal numerical simulation: Concepts and advances
NASA Astrophysics Data System (ADS)
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.
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.
Numerical methods for molecular dynamics
Skeel, R.D.
1991-01-01
This report summarizes our research progress to date on the use of multigrid methods for three-dimensional elliptic partial differential equations, with particular emphasis on application to the Poisson-Boltzmann equation of molecular biophysics. This research is motivated by the need for fast and accurate numerical solution techniques for three-dimensional problems arising in physics and engineering. In many applications these problems must be solved repeatedly, and the extremely large number of discrete unknowns required to accurately approximate solutions to partial differential equations in three-dimensional regions necessitates the use of efficient solution methods. This situation makes clear the importance of developing methods which are of optimal order (or nearly so), meaning that the number of operations required to solve the discrete problem is on the order of the number of discrete unknowns. Multigrid methods are generally regarded as being in this class of methods, and are in fact provably optimal order for an increasingly large class of problems. The fundamental goal of this research is to develop a fast and accurate numerical technique, based on multi-level principles, for the solutions of the Poisson-Boltzmann equation of molecular biophysics and similar equations occurring in other applications. An outline of the report is as follows. We first present some background material, followed by a survey of the literature on the use of multigrid methods for solving problems similar to the Poisson-Boltzmann equation. A short description of the software we have developed so far is then given, and numerical results are discussed. Finally, our research plans for the coming year are presented.
Numerical methods for multibody systems
NASA Technical Reports Server (NTRS)
Glowinski, Roland; Nasser, Mahmoud G.
1994-01-01
This article gives a brief summary of some results obtained by Nasser on modeling and simulation of inequality problems in multibody dynamics. In particular, the augmented Lagrangian method discussed here is applied to a constrained motion problem with impulsive inequality constraints. A fundamental characteristic of the multibody dynamics problem is the lack of global convexity of its Lagrangian. The problem is transformed into a convex analysis problem by localization (piecewise linearization), where the augmented Lagrangian has been successfully used. A model test problem is considered and a set of numerical experiments is presented.
Some recent advances in the numerical solution of differential equations
NASA Astrophysics Data System (ADS)
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.
NASA Astrophysics Data System (ADS)
Yoshida, Hiroyuki; Nagayoshi, Takuji; Takase, Kazuyuki; Akimoto, Hajime
Thermal-hydraulic design of the current boiling water reactor (BWR) is performed by correlations with empirical results of actual-size tests. However, 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. 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, detailed Two-Phase Flow simulation code using advanced Interface Tracking method: TPFIT is developed to calculate the detailed information of the two-phase flow. We tried to verify the TPFIT code by comparing it with the 2-channel air-water and steam-water mixing experimental results. The predicted result agrees well the observed results and bubble dynamics through the gap and cross flow behavior could be effectively predicted by the TPFIT code, and pressure difference between fluid channels is responsible for the fluid mixing.
Numerical methods used in fusion science numerical modeling
NASA Astrophysics Data System (ADS)
Yagi, M.
2015-04-01
The dynamics of burning plasma is very complicated physics, which is dominated by multi-scale and multi-physics phenomena. To understand such phenomena, numerical simulations are indispensable. Fundamentals of numerical methods used in fusion science numerical modeling are briefly discussed in this paper. In addition, the parallelization technique such as open multi processing (OpenMP) and message passing interface (MPI) parallel programing are introduced and the loop-level parallelization is shown as an example.
RELAP-7 Numerical Stabilization: Entropy Viscosity Method
R. A. Berry; M. O. Delchini; J. Ragusa
2014-06-01
The RELAP-7 code is the next generation nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL). The code is based on the INL's modern scientific software development framework, MOOSE (Multi-Physics Object Oriented Simulation Environment). The overall design goal of RELAP-7 is to take advantage of the previous thirty years of advancements in computer architecture, software design, numerical integration methods, and physical models. The end result will be a reactor systems analysis capability that retains and improves upon RELAP5's capability and extends the analysis capability for all reactor system simulation scenarios. RELAP-7 utilizes a single phase and a novel seven-equation two-phase flow models as described in the RELAP-7 Theory Manual (INL/EXT-14-31366). The basic equation systems are hyperbolic, which generally require some type of stabilization (or artificial viscosity) to capture nonlinear discontinuities and to suppress advection-caused oscillations. This report documents one of the available options for this stabilization in RELAP-7 -- a new and novel approach known as the entropy viscosity method. Because the code is an ongoing development effort in which the physical sub models, numerics, and coding are evolving, so too must the specific details of the entropy viscosity stabilization method. Here the fundamentals of the method in their current state are presented.
Numerical methods for engine-airframe integration
Murthy, S.N.B.; Paynter, G.C.
1986-01-01
Various papers on numerical methods for engine-airframe integration are presented. The individual topics considered include: scientific computing environment for the 1980s, overview of prediction of complex turbulent flows, numerical solutions of the compressible Navier-Stokes equations, elements of computational engine/airframe integrations, computational requirements for efficient engine installation, application of CAE and CFD techniques to complete tactical missile design, CFD applications to engine/airframe integration, and application of a second-generation low-order panel methods to powerplant installation studies. Also addressed are: three-dimensional flow analysis of turboprop inlet and nacelle configurations, application of computational methods to the design of large turbofan engine nacelles, comparison of full potential and Euler solution algorithms for aeropropulsive flow field computations, subsonic/transonic, supersonic nozzle flows and nozzle integration, subsonic/transonic prediction capabilities for nozzle/afterbody configurations, three-dimensional viscous design methodology of supersonic inlet systems for advanced technology aircraft, and a user's technology assessment.
Advanced Numerical Model for Irradiated Concrete
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
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.
Advanced Numerical Modeling of Turbulent Atmospheric Flows
NASA Astrophysics Data System (ADS)
Kühnlein, Christian; Dörnbrack, Andreas; Gerz, Thomas
The present chapter introduces the method of computational simulation to predict and study turbulent atmospheric flows. This includes a description of the fundamental approach to computational simulation and the practical implementation using the technique of large-eddy simulation. In addition, selected contributions from IPA scientists to computational model development and various examples for applications are given. These examples include homogeneous turbulence, convective boundary layers, heated forest canopy, buoyant thermals, and large-scale flows with baroclinic wave instability.
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.
Numerical Methods for Radiation Magnetohydrodynamics in Astrophysics
Klein, R I; Stone, J M
2007-11-20
We describe numerical methods for solving the equations of radiation magnetohydrodynamics (MHD) for astrophysical fluid flow. Such methods are essential for the investigation of the time-dependent and multidimensional dynamics of a variety of astrophysical systems, although our particular interest is motivated by problems in star formation. Over the past few years, the authors have been members of two parallel code development efforts, and this review reflects that organization. In particular, we discuss numerical methods for MHD as implemented in the Athena code, and numerical methods for radiation hydrodynamics as implemented in the Orion code. We discuss the challenges introduced by the use of adaptive mesh refinement in both codes, as well as the most promising directions for future developments.
25 Years of Self-organized Criticality: Numerical Detection Methods
NASA Astrophysics Data System (ADS)
McAteer, R. T. James; Aschwanden, Markus J.; Dimitropoulou, Michaila; Georgoulis, Manolis K.; Pruessner, Gunnar; Morales, Laura; Ireland, Jack; Abramenko, Valentyna
2016-01-01
The detection and characterization of self-organized criticality (SOC), in both real and simulated data, has undergone many significant revisions over the past 25 years. The explosive advances in the many numerical methods available for detecting, discriminating, and ultimately testing, SOC have played a critical role in developing our understanding of how systems experience and exhibit SOC. In this article, methods of detecting SOC are reviewed; from correlations to complexity to critical quantities. A description of the basic autocorrelation method leads into a detailed analysis of application-oriented methods developed in the last 25 years. In the second half of this manuscript space-based, time-based and spatial-temporal methods are reviewed and the prevalence of power laws in nature is described, with an emphasis on event detection and characterization. The search for numerical methods to clearly and unambiguously detect SOC in data often leads us outside the comfort zone of our own disciplines—the answers to these questions are often obtained by studying the advances made in other fields of study. In addition, numerical detection methods often provide the optimum link between simulations and experiments in scientific research. We seek to explore this boundary where the rubber meets the road, to review this expanding field of research of numerical detection of SOC systems over the past 25 years, and to iterate forwards so as to provide some foresight and guidance into developing breakthroughs in this subject over the next quarter of a century.
On Numerical Methods For Hypersonic Turbulent Flows
NASA Astrophysics Data System (ADS)
Yee, H. C.; Sjogreen, B.; Shu, C. W.; Wang, W.; Magin, T.; Hadjadj, A.
2011-05-01
Proper control of numerical dissipation in numerical methods beyond the standard shock-capturing dissipation at discontinuities is an essential element for accurate and stable simulation of hypersonic turbulent flows, including combustion, and thermal and chemical nonequilibrium flows. Unlike rapidly developing shock interaction flows, turbulence computations involve long time integrations. Improper control of numerical dissipation from one time step to another would be compounded over time, resulting in the smearing of turbulent fluctuations to an unrecognizable form. Hypersonic turbulent flows around re- entry space vehicles involve mixed steady strong shocks and turbulence with unsteady shocklets that pose added computational challenges. Stiffness of the source terms and material mixing in combustion pose yet other types of numerical challenges. A low dissipative high order well- balanced scheme, which can preserve certain non-trivial steady solutions of the governing equations exactly, may help minimize some of these difficulties. For stiff reactions it is well known that the wrong propagation speed of discontinuities occurs due to the under-resolved numerical solutions in both space and time. Schemes to improve the wrong propagation speed of discontinuities for systems of stiff reacting flows remain a challenge for algorithm development. Some of the recent algorithm developments for direct numerical simulations (DNS) and large eddy simulations (LES) for the subject physics, including the aforementioned numerical challenges, will be discussed.
A numerical method of detecting singularity
NASA Technical Reports Server (NTRS)
Laporte, M.; Vignes, J.
1978-01-01
A numerical method is reported which determines a value C for the degree of conditioning of a matrix. This value is C = 0 for a singular matrix and has progressively larger values for matrices which are increasingly well-conditioned. This value is C sub = C max sub max (C defined by the precision of the computer) when the matrix is perfectly well conditioned.
A numerical method for predicting hypersonic flowfields
NASA Technical Reports Server (NTRS)
Maccormack, Robert W.; Candler, Graham V.
1989-01-01
The flow about a body traveling at hypersonic speed is energetic enough to cause the atmospheric gases to chemically react and reach states in thermal nonequilibrium. The prediction of hypersonic flowfields requires a numerical method capable of solving the conservation equations of fluid flow, the chemical rate equations for specie formation and dissociation, and the transfer of energy relations between translational and vibrational temperature states. Because the number of equations to be solved is large, the numerical method should also be as efficient as possible. The proposed paper presents a fully implicit method that fully couples the solution of the fluid flow equations with the gas physics and chemistry relations. The method flux splits the inviscid flow terms, central differences of the viscous terms, preserves element conservation in the strong chemistry source terms, and solves the resulting block matrix equation by Gauss Seidel line relaxation.
The TAB method for numerical calculation of spray droplet breakup
NASA Astrophysics Data System (ADS)
Orourke, P. J.; Amsden, A. A.
A short history is given of the major milestones in the development of the stochastic particle method for calculating liquid fuel sprays. The most recent advance has been the discovery of the importance of drop breakup in engine sprays. A new method, called TAB, for calculating drop breakup is presented. Some theoretical properties of the method are derived; its numerical implementation in the computer program KIVA is described; and comparisons are presented between TAB-method calculations and experiments and calculations using another breakup model.
Numerical Analysis of the Symmetric Methods
NASA Astrophysics Data System (ADS)
Xu, Ji-Hong; Zhang, A.-Li
1995-03-01
Aimed at the initial value problem of the particular second-order ordinary differential equations,y ″=f(x, y), the symmetric methods (Quinlan and Tremaine, 1990) and our methods (Xu and Zhang, 1994) have been compared in detail by integrating the artificial earth satellite orbits in this paper. In the end, we point out clearly that the integral accuracy of numerical integration of the satellite orbits by applying our methods is obviously higher than that by applying the same order formula of the symmetric methods when the integration time-interval is not greater than 12000 periods.
Recent advances in two-phase flow numerics
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.
Numerical Methods for Stochastic Partial Differential Equations
Sharp, D.H.; Habib, S.; Mineev, M.B.
1999-07-08
This is the final report of a Laboratory Directed Research and Development (LDRD) project at the Los Alamos National laboratory (LANL). The objectives of this proposal were (1) the development of methods for understanding and control of spacetime discretization errors in nonlinear stochastic partial differential equations, and (2) the development of new and improved practical numerical methods for the solutions of these equations. The authors have succeeded in establishing two methods for error control: the functional Fokker-Planck equation for calculating the time discretization error and the transfer integral method for calculating the spatial discretization error. In addition they have developed a new second-order stochastic algorithm for multiplicative noise applicable to the case of colored noises, and which requires only a single random sequence generation per time step. All of these results have been verified via high-resolution numerical simulations and have been successfully applied to physical test cases. They have also made substantial progress on a longstanding problem in the dynamics of unstable fluid interfaces in porous media. This work has lead to highly accurate quasi-analytic solutions of idealized versions of this problem. These may be of use in benchmarking numerical solutions of the full stochastic PDEs that govern real-world problems.
Numerical methods for supersonic astrophysical jets
NASA Astrophysics Data System (ADS)
Ha, Youngsoo
2003-09-01
The Euler equations of gas dynamics are used for the simulation of general astrophysical fluid flows including high Mach number astrophysical jets with radiative cooling. To accurately compute supersonic jet solutions with sharp resolution of shock waves, three modern numerical methods for gas dynamics were used: (1)a second-order Godunov method in LeVeque's software package CLAWPACK, (2)the Nessyahu-Tadmor-Kurganov (NTK) central hyperbolic scheme, and (3)the WENO-LF (Weighted Essentially Non-Oscillatory Lax-Friedrichs) scheme. Then simulations of supersonic astrophysical jets were compared, first without and then with radiative cooling. CLAWPACK consists of routines for solving time-dependent nonlinear hyperbolic conservation laws based on higher order Godunov methods and approximate Riemann problem solutions; the NTK scheme solves conservation laws using a modified Lax-Friedrichs central difference method without appealing to Riemann problem solutions; and the WENO-LF finite difference scheme is based on the Essentially Non-Oscillatory (ENO) idea by using Lax- Friedrichs flux splitting. The ENO method constructs a solution using the smoothness of the interpolating polynomial on given stencils; on the other hand, the WENO scheme uses a convex combination of the interpolate functions on all candidate stencils. The third-order and fifth-order WENO-LF methods were used to simulate the high Mach number jets. Appropriate numerical methods for incorporating radiative cooling in these numerical methods are also discussed. Interactions of supersonic jets with their environments (jet-“blob” interactions) are shown after modifying the codes to handle high Mach numbers and radiative cooling.
Hyperbolic conservation laws and numerical methods
NASA Technical Reports Server (NTRS)
Leveque, Randall J.
1990-01-01
The mathematical structure of hyperbolic systems and the scalar equation case of conservation laws are discussed. Linear, nonlinear systems and the Riemann problem for the Euler equations are also studied. The numerical methods for conservation laws are presented in a nonstandard manner which leads to large time steps generalizations and computations on irregular grids. The solution of conservation laws with stiff source terms is examined.
Numerical methods for molecular dynamics. Progress report
Skeel, R.D.
1991-12-31
This report summarizes our research progress to date on the use of multigrid methods for three-dimensional elliptic partial differential equations, with particular emphasis on application to the Poisson-Boltzmann equation of molecular biophysics. This research is motivated by the need for fast and accurate numerical solution techniques for three-dimensional problems arising in physics and engineering. In many applications these problems must be solved repeatedly, and the extremely large number of discrete unknowns required to accurately approximate solutions to partial differential equations in three-dimensional regions necessitates the use of efficient solution methods. This situation makes clear the importance of developing methods which are of optimal order (or nearly so), meaning that the number of operations required to solve the discrete problem is on the order of the number of discrete unknowns. Multigrid methods are generally regarded as being in this class of methods, and are in fact provably optimal order for an increasingly large class of problems. The fundamental goal of this research is to develop a fast and accurate numerical technique, based on multi-level principles, for the solutions of the Poisson-Boltzmann equation of molecular biophysics and similar equations occurring in other applications. An outline of the report is as follows. We first present some background material, followed by a survey of the literature on the use of multigrid methods for solving problems similar to the Poisson-Boltzmann equation. A short description of the software we have developed so far is then given, and numerical results are discussed. Finally, our research plans for the coming year are presented.
Advanced probabilistic method of development
NASA Technical Reports Server (NTRS)
Wirsching, P. H.
1987-01-01
Advanced structural reliability methods are utilized on the Probabilistic Structural Analysis Methods (PSAM) project to provide a tool for analysis and design of space propulsion system hardware. The role of the effort at the University of Arizona is to provide reliability technology support to this project. PSAM computer programs will provide a design tool for analyzing uncertainty associated with thermal and mechanical loading, material behavior, geometry, and the analysis methods used. Specifically, reliability methods are employed to perform sensitivity analyses, to establish the distribution of a critical response variable (e.g., stress, deflection), to perform reliability assessment, and ultimately to produce a design which will minimize cost and/or weight. Uncertainties in the design factors of space propulsion hardware are described by probability models constructed using statistical analysis of data. Statistical methods are employed to produce a probability model, i.e., a statistical synthesis or summary of each design variable in a format suitable for reliability analysis and ultimately, design decisions.
Numerical methods for finding stationary gravitational solutions
NASA Astrophysics Data System (ADS)
Dias, Óscar J. C.; Santos, Jorge E.; Way, Benson
2016-07-01
The wide applications of higher dimensional gravity and gauge/gravity duality have fuelled the search for new stationary solutions of the Einstein equation (possibly coupled to matter). In this topical review, we explain the mathematical foundations and give a practical guide for the numerical solution of gravitational boundary value problems. We present these methods by way of example: resolving asymptotically flat black rings, singly spinning lumpy black holes in anti-de Sitter (AdS), and the Gregory–Laflamme zero modes of small rotating black holes in AdS{}5× {S}5. We also include several tools and tricks that have been useful throughout the literature.
Numerical analysis method for linear induction machines.
NASA Technical Reports Server (NTRS)
Elliott, D. G.
1972-01-01
A numerical analysis method has been developed for linear induction machines such as liquid metal MHD pumps and generators and linear motors. Arbitrary phase currents or voltages can be specified and the moving conductor can have arbitrary velocity and conductivity variations from point to point. The moving conductor is divided into a mesh and coefficients are calculated for the voltage induced at each mesh point by unit current at every other mesh point. Combining the coefficients with the mesh resistances yields a set of simultaneous equations which are solved for the unknown currents.
Automatic numerical integration methods for Feynman integrals through 3-loop
NASA Astrophysics Data System (ADS)
de Doncker, E.; Yuasa, F.; Kato, K.; Ishikawa, T.; Olagbemi, O.
2015-05-01
We give numerical integration results for Feynman loop diagrams through 3-loop such as those covered by Laporta [1]. The methods are based on automatic adaptive integration, using iterated integration and extrapolation with programs from the QUADPACK package, or multivariate techniques from the ParInt package. The Dqags algorithm from QuadPack accommodates boundary singularities of fairly general types. PARINT is a package for multivariate integration layered over MPI (Message Passing Interface), which runs on clusters and incorporates advanced parallel/distributed techniques such as load balancing among processes that may be distributed over a network of nodes. Results are included for 3-loop self-energy diagrams without IR (infra-red) or UV (ultra-violet) singularities. A procedure based on iterated integration and extrapolation yields a novel method of numerical regularization for integrals with UV terms, and is applied to a set of 2-loop self-energy diagrams with UV singularities.
Application of numerical methods to elasticity imaging.
Castaneda, Benjamin; Ormachea, Juvenal; Rodríguez, Paul; Parker, Kevin J
2013-03-01
Elasticity imaging can be understood as the intersection of the study of biomechanical properties, imaging sciences, and physics. It was mainly motivated by the fact that pathological tissue presents an increased stiffness when compared to surrounding normal tissue. In the last two decades, research on elasticity imaging has been an international and interdisciplinary pursuit aiming to map the viscoelastic properties of tissue in order to provide clinically useful information. As a result, several modalities of elasticity imaging, mostly based on ultrasound but also on magnetic resonance imaging and optical coherence tomography, have been proposed and applied to a number of clinical applications: cancer diagnosis (prostate, breast, liver), hepatic cirrhosis, renal disease, thyroiditis, arterial plaque evaluation, wall stiffness in arteries, evaluation of thrombosis in veins, and many others. In this context, numerical methods are applied to solve forward and inverse problems implicit in the algorithms in order to estimate viscoelastic linear and nonlinear parameters, especially for quantitative elasticity imaging modalities. In this work, an introduction to elasticity imaging modalities is presented. The working principle of qualitative modalities (sonoelasticity, strain elastography, acoustic radiation force impulse) and quantitative modalities (Crawling Waves Sonoelastography, Spatially Modulated Ultrasound Radiation Force (SMURF), Supersonic Imaging) will be explained. Subsequently, the areas in which numerical methods can be applied to elasticity imaging are highlighted and discussed. Finally, we present a detailed example of applying total variation and AM-FM techniques to the estimation of elasticity. PMID:24010245
Mathematica with a Numerical Methods Course
NASA Astrophysics Data System (ADS)
Varley, Rodney
2003-04-01
An interdisciplinary "Numerical Methods" course has been shared between physics, mathematics and computer science since 1992 at Hunter C. Recently, the lectures and workshops for this course have become formalized and placed on the internet at http://www.ph.hunter.cuny.edu (follow the links "Course Listings and Websites" >> "PHYS385 (Numerical Methods)". Mathematica notebooks for the lectures are available for automatic download (by "double clicking" the lecture icon) for student use in the classroom or at home. AOL (or Netscape/Explorer) can be used provided Mathematica (or the "free" MathReader) has been made a "helper application". Using Mathematica has the virtue that mathematical equations (no LaTex required) can easily be included with the text and Mathematica's graphing is easy to use. Computational cells can be included within the notebook and students may easily modify the calculation to see the result of "what if..." questions. Homework is sent as Mathematica notebooks to the instructor via the internet and the corrected workshops are returned in the same manner. Most exam questions require computational solutions.
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.
Numerical methods for problems in computational aeroacoustics
NASA Astrophysics Data System (ADS)
Mead, Jodi Lorraine
1998-12-01
A goal of computational aeroacoustics is the accurate calculation of noise from a jet in the far field. This work concerns the numerical aspects of accurately calculating acoustic waves over large distances and long time. More specifically, the stability, efficiency, accuracy, dispersion and dissipation in spatial discretizations, time stepping schemes, and absorbing boundaries for the direct solution of wave propagation problems are determined. Efficient finite difference methods developed by Tam and Webb, which minimize dispersion and dissipation, are commonly used for the spatial and temporal discretization. Alternatively, high order pseudospectral methods can be made more efficient by using the grid transformation introduced by Kosloff and Tal-Ezer. Work in this dissertation confirms that the grid transformation introduced by Kosloff and Tal-Ezer is not spectrally accurate because, in the limit, the grid transformation forces zero derivatives at the boundaries. If a small number of grid points are used, it is shown that approximations with the Chebyshev pseudospectral method with the Kosloff and Tal-Ezer grid transformation are as accurate as with the Chebyshev pseudospectral method. This result is based on the analysis of the phase and amplitude errors of these methods, and their use for the solution of a benchmark problem in computational aeroacoustics. For the grid transformed Chebyshev method with a small number of grid points it is, however, more appropriate to compare its accuracy with that of high- order finite difference methods. This comparison, for an order of accuracy 10-3 for a benchmark problem in computational aeroacoustics, is performed for the grid transformed Chebyshev method and the fourth order finite difference method of Tam. Solutions with the finite difference method are as accurate. and the finite difference method is more efficient than, the Chebyshev pseudospectral method with the grid transformation. The efficiency of the Chebyshev
Numerical Forming Simulations and Optimisation in Advanced Materials
NASA Astrophysics Data System (ADS)
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.
An Advanced Manipulator For Poisson Series With Numerical Coefficients
NASA Astrophysics Data System (ADS)
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).
Numerical methods for analyzing electromagnetic scattering
NASA Technical Reports Server (NTRS)
Lee, S. W.; Lo, Y. T.; Chuang, S. L.; Lee, C. S.
1985-01-01
Numerical methods to analyze electromagnetic scattering are presented. The dispersions and attenuations of the normal modes in a circular waveguide coated with lossy material were completely analyzed. The radar cross section (RCS) from a circular waveguide coated with lossy material was calculated. The following is observed: (1) the interior irradiation contributes to the RCS much more than does the rim diffraction; (2) at low frequency, the RCS from the circular waveguide terminated by a perfect electric conductor (PEC) can be reduced more than 13 dB down with a coating thickness less than 1% of the radius using the best lossy material available in a 6 radius-long cylinder; (3) at high frequency, a modal separation between the highly attenuated and the lowly attenuated modes is evident if the coating material is too lossy, however, a large RCS reduction can be achieved for a small incident angle with a thin layer of coating. It is found that the waveguide coated with a lossy magnetic material can be used as a substitute for a corrugated waveguide to produce a circularly polarized radiation yield.
A numerical method for cardiac mechanoelectric simulations.
Pathmanathan, Pras; Whiteley, Jonathan P
2009-05-01
Much effort has been devoted to developing numerical techniques for solving the equations that describe cardiac electrophysiology, namely the monodomain equations and bidomain equations. Only a limited selection of publications, however, address the development of numerical techniques for mechanoelectric simulations where cardiac electrophysiology is coupled with deformation of cardiac tissue. One problem commonly encountered in mechanoelectric simulations is instability of the coupled numerical scheme. In this study, we develop a stable numerical scheme for mechanoelectric simulations. A number of convergence tests are carried out using this stable technique for simulations where deformations are of the magnitude typically observed in a beating heart. These convergence tests demonstrate that accurate computation of tissue deformation requires a nodal spacing of around 1 mm in the mesh used to calculate tissue deformation. This is a much finer computational grid than has previously been acknowledged, and has implications for the computational efficiency of the resulting numerical scheme. PMID:19263223
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.
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
Advanced surface paneling method for subsonic and supersonic flow
NASA Technical Reports Server (NTRS)
Erickson, L. L.; Johnson, F. T.; Ehlers, F. E.
1976-01-01
Numerical results illustrating the capabilities of an advanced aerodynamic surface paneling method are presented. The method is applicable to both subsonic and supersonic flow, as represented by linearized potential flow theory. The method is based on linearly varying sources and quadratically varying doublets which are distributed over flat or curved panels. These panels are applied to the true surface geometry of arbitrarily shaped three dimensional aerodynamic configurations.
Novel Parallel Numerical Methods for Radiation& Neutron Transport
Brown, P N
2001-03-06
In many of the multiphysics simulations performed at LLNL, transport calculations can take up 30 to 50% of the total run time. If Monte Carlo methods are used, the percentage can be as high as 80%. Thus, a significant core competence in the formulation, software implementation, and solution of the numerical problems arising in transport modeling is essential to Laboratory and DOE research. In this project, we worked on developing scalable solution methods for the equations that model the transport of photons and neutrons through materials. Our goal was to reduce the transport solve time in these simulations by means of more advanced numerical methods and their parallel implementations. These methods must be scalable, that is, the time to solution must remain constant as the problem size grows and additional computer resources are used. For iterative methods, scalability requires that (1) the number of iterations to reach convergence is independent of problem size, and (2) that the computational cost grows linearly with problem size. We focused on deterministic approaches to transport, building on our earlier work in which we performed a new, detailed analysis of some existing transport methods and developed new approaches. The Boltzmann equation (the underlying equation to be solved) and various solution methods have been developed over many years. Consequently, many laboratory codes are based on these methods, which are in some cases decades old. For the transport of x-rays through partially ionized plasmas in local thermodynamic equilibrium, the transport equation is coupled to nonlinear diffusion equations for the electron and ion temperatures via the highly nonlinear Planck function. We investigated the suitability of traditional-solution approaches to transport on terascale architectures and also designed new scalable algorithms; in some cases, we investigated hybrid approaches that combined both.
Numerical Methods for Forward and Inverse Problems in Discontinuous Media
Chartier, Timothy P.
2011-03-08
The research emphasis under this grant's funding is in the area of algebraic multigrid methods. The research has two main branches: 1) exploring interdisciplinary applications in which algebraic multigrid can make an impact and 2) extending the scope of algebraic multigrid methods with algorithmic improvements that are based in strong analysis.The work in interdisciplinary applications falls primarily in the field of biomedical imaging. Work under this grant demonstrated the effectiveness and robustness of multigrid for solving linear systems that result from highly heterogeneous finite element method models of the human head. The results in this work also give promise to medical advances possible with software that may be developed. Research to extend the scope of algebraic multigrid has been focused in several areas. In collaboration with researchers at the University of Colorado, Lawrence Livermore National Laboratory, and Los Alamos National Laboratory, the PI developed an adaptive multigrid with subcycling via complementary grids. This method has very cheap computing costs per iterate and is showing promise as a preconditioner for conjugate gradient. Recent work with Los Alamos National Laboratory concentrates on developing algorithms that take advantage of the recent advances in adaptive multigrid research. The results of the various efforts in this research could ultimately have direct use and impact to researchers for a wide variety of applications, including, astrophysics, neuroscience, contaminant transport in porous media, bi-domain heart modeling, modeling of tumor growth, and flow in heterogeneous porous media. This work has already led to basic advances in computational mathematics and numerical linear algebra and will continue to do so into the future.
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.
Advanced accelerator methods: The cyclotrino
Welch, J.J.; Bertsche, K.J.; Friedman, P.G.; Morris, D.E.; Muller, R.A.
1987-04-01
Several new and unusual, advanced techniques in the small cyclotron are described. The cyclotron is run at low energy, using negative ions and at high harmonics. Electrostatic focusing is used exclusively. The ion source and injection system is in the center, which unfortunately does not provide enough current, but the new system design should solve this problem. An electrostatic extractor that runs at low voltage, under 5 kV, and a microchannel plate detector which is able to discriminate low energy ions from the /sup 14/C are used. The resolution is sufficient for /sup 14/C dating and a higher intensity source should allow dating of a milligram size sample of 30,000 year old material with less than 10% uncertainty.
Advanced reliability methods - A review
NASA Astrophysics Data System (ADS)
Forsyth, David S.
2016-02-01
There are a number of challenges to the current practices for Probability of Detection (POD) assessment. Some Nondestructive Testing (NDT) methods, especially those that are image-based, may not provide a simple relationship between a scalar NDT response and a damage size. Some damage types are not easily characterized by a single scalar metric. Other sensing paradigms, such as structural health monitoring, could theoretically replace NDT but require a POD estimate. And the cost of performing large empirical studies to estimate POD can be prohibitive. The response of the research community has been to develop new methods that can be used to generate the same information, POD, in a form that can be used by engineering designers. This paper will highlight approaches to image-based data and complex defects, Model Assisted POD estimation, and Bayesian methods for combining information. This paper will also review the relationship of the POD estimate, confidence bounds, tolerance bounds, and risk assessment.
Finite element methods in numerical relativity.
NASA Astrophysics Data System (ADS)
Mann, P. J.
The finite element method is very successful in Newtonian fluid simulations, and can be extended to relativitstic fluid flows. This paper describes the general method, and then outlines some preliminary results for spherically symmetric geometries. The mixed finite element - finite difference scheme is introduced, and used for the description of spherically symmetric collapse. Baker's (Newtonian) shock modelling method and Miller's moving finite element method are also mentioned. Collapse in double-null coordinates requires non-constant time slicing, so the full finite element method in space and time is described.
Advanced Fine Particulate Characterization Methods
Steven Benson; Lingbu Kong; Alexander Azenkeng; Jason Laumb; Robert Jensen; Edwin Olson; Jill MacKenzie; A.M. Rokanuzzaman
2007-01-31
The characterization and control of emissions from combustion sources are of significant importance in improving local and regional air quality. Such emissions include fine particulate matter, organic carbon compounds, and NO{sub x} and SO{sub 2} gases, along with mercury and other toxic metals. This project involved four activities including Further Development of Analytical Techniques for PM{sub 10} and PM{sub 2.5} Characterization and Source Apportionment and Management, Organic Carbonaceous Particulate and Metal Speciation for Source Apportionment Studies, Quantum Modeling, and High-Potassium Carbon Production with Biomass-Coal Blending. The key accomplishments included the development of improved automated methods to characterize the inorganic and organic components particulate matter. The methods involved the use of scanning electron microscopy and x-ray microanalysis for the inorganic fraction and a combination of extractive methods combined with near-edge x-ray absorption fine structure to characterize the organic fraction. These methods have direction application for source apportionment studies of PM because they provide detailed inorganic analysis along with total organic and elemental carbon (OC/EC) quantification. Quantum modeling using density functional theory (DFT) calculations was used to further elucidate a recently developed mechanistic model for mercury speciation in coal combustion systems and interactions on activated carbon. Reaction energies, enthalpies, free energies and binding energies of Hg species to the prototype molecules were derived from the data obtained in these calculations. Bimolecular rate constants for the various elementary steps in the mechanism have been estimated using the hard-sphere collision theory approximation, and the results seem to indicate that extremely fast kinetics could be involved in these surface reactions. Activated carbon was produced from a blend of lignite coal from the Center Mine in North Dakota and
A Hybrid Numerical Analysis Method for Structural Health Monitoring
NASA Technical Reports Server (NTRS)
Forth, Scott C.; Staroselsky, Alexander
2001-01-01
A new hybrid surface-integral-finite-element numerical scheme has been developed to model a three-dimensional crack propagating through a thin, multi-layered coating. The finite element method was used to model the physical state of the coating (far field), and the surface integral method was used to model the fatigue crack growth. The two formulations are coupled through the need to satisfy boundary conditions on the crack surface and the external boundary. The coupling is sufficiently weak that the surface integral mesh of the crack surface and the finite element mesh of the uncracked volume can be set up independently. Thus when modeling crack growth, the finite element mesh can remain fixed for the duration of the simulation as the crack mesh is advanced. This method was implemented to evaluate the feasibility of fabricating a structural health monitoring system for real-time detection of surface cracks propagating in engine components. In this work, the authors formulate the hybrid surface-integral-finite-element method and discuss the mechanical issues of implementing a structural health monitoring system in an aircraft engine environment.
Numerical matrix method for quantum periodic potentials
NASA Astrophysics Data System (ADS)
Le Vot, Felipe; Meléndez, Juan J.; Yuste, Santos B.
2016-06-01
A numerical matrix methodology is applied to quantum problems with periodic potentials. The procedure consists essentially in replacing the true potential by an alternative one, restricted by an infinite square well, and in expressing the wave functions as finite superpositions of eigenfunctions of the infinite well. A matrix eigenvalue equation then yields the energy levels of the periodic potential within an acceptable accuracy. The methodology has been successfully used to deal with problems based on the well-known Kronig-Penney (KP) model. Besides the original model, these problems are a dimerized KP solid, a KP solid containing a surface, and a KP solid under an external field. A short list of additional problems that can be solved with this procedure is presented.
Method for numerical simulations of metastable states
Heller, U.M.; Seiberg, N.
1983-06-15
We present a numerical simulation of metastable states near a first-order phase transition in the example of a U(1) lattice gauge theory with a generalized action. In order to make measurements in these states possible their decay has to be prevented. We achieve this by using a microcanonical simulation for a finite system. We then obtain the coupling constant (inverse temperature) as a function of the action density. It turns out to be nonmonotonic and hence not uniquely invertible. From it we derive the effective potential for the action density. This effective potential is not always convex, a property that seems to be in contradiction with the standard lore about its convexity. This apparent ''paradox'' is resolved in a discussion about different definitions of the effective potential.
Recent advances in lattice Boltzmann methods
Chen, S.; Doolen, G.D.; He, X.; Nie, X.; Zhang, R.
1998-12-31
In this paper, the authors briefly present the basic principles of lattice Boltzmann method and summarize recent advances of the method, including the application of the lattice Boltzmann method for fluid flows in MEMS and simulation of the multiphase mixing and turbulence.
Numerical methods in Markov chain modeling
NASA Technical Reports Server (NTRS)
Philippe, Bernard; Saad, Youcef; Stewart, William J.
1989-01-01
Several methods for computing stationary probability distributions of Markov chains are described and compared. The main linear algebra problem consists of computing an eigenvector of a sparse, usually nonsymmetric, matrix associated with a known eigenvalue. It can also be cast as a problem of solving a homogeneous singular linear system. Several methods based on combinations of Krylov subspace techniques are presented. The performance of these methods on some realistic problems are compared.
Interpolation Method Needed for Numerical Uncertainty
NASA Technical Reports Server (NTRS)
Groves, Curtis E.; Ilie, Marcel; Schallhorn, Paul A.
2014-01-01
Using Computational Fluid Dynamics (CFD) to predict a flow field is an approximation to the exact problem and uncertainties exist. There is a method to approximate the errors in CFD via Richardson's Extrapolation. This method is based off of progressive grid refinement. To estimate the errors, the analyst must interpolate between at least three grids. This paper describes a study to find an appropriate interpolation scheme that can be used in Richardson's extrapolation or other uncertainty method to approximate errors.
An advanced numerical model for phase change problems in complicated geometries
NASA Astrophysics Data System (ADS)
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.
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…
Simultaneous source-mask optimization: a numerical combining method
NASA Astrophysics Data System (ADS)
Mülders, Thomas; Domnenko, Vitaliy; Küchler, Bernd; Klimpel, Thomas; Stock, Hans-Jürgen; Poonawala, Amyn A.; Taravade, Kunal N.; Stanton, William A.
2010-09-01
A new method for simultaneous Source-Mask Optimization (SMO) is presented. In order to produce optimum imaging fidelity with respect to exposure lattitude, depth of focus (DoF) and mask error enhancement factor (MEEF) the presented method aims to leverage both, the available degrees of freedom of a pixelated source and those available for the mask layout. The approach described in this paper is designed as to work with dissected mask polygons. The dissection of the mask patterns is to be performed in advance (before SMO) with the Synopsys Proteus OPC engine, providing the available degrees of freedom for mask pattern optimization. This is similar to mask optimization done for optical proximity correction (OPC). Additionally, however, the illumination source will be simultaneously optimized. The SMO approach borrows many of the performance enhancement methods of OPC software for mask correction, but is especially designed as to simultaneously optimize a pixelated source shape as nowadays available in production environments. Designed as a numerical optimization approach the method is able to assess in acceptable times several hundreds of thousands source-mask combinations for small, critical layout snippets. This allows a global optimization scheme to be applied to the SMO problem which is expected to better explore the optimization space and thus to yield an improved solution quality compared to local optimizations methods. The method is applied to an example system for investigating the impact of source constraints on the SMO results. Also, it is investigated how well possibly conflicting goals of low MEEF and large DoF can be balanced.
Modelling asteroid brightness variations. I - Numerical methods
NASA Technical Reports Server (NTRS)
Karttunen, H.
1989-01-01
A method for generating lightcurves of asteroid models is presented. The effects of the shape of the asteroid and the scattering law of a surface element are distinctly separable, being described by chosen functions that can easily be changed. The shape is specified by means of two functions that yield the length of the radius vector and the normal vector of the surface at a given point. The general shape must be convex, but spherical concavities producing macroscopic shadowing can also be modeled.
The Role of Numerical Simulation in Advancing Plasma Propulsion
NASA Astrophysics Data System (ADS)
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.
Numerical Simulations and Optimisation in Forming of Advanced Materials
NASA Astrophysics Data System (ADS)
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.
A numerical method for power plant simulations
Carcasci, C.; Facchini, B.
1996-03-01
This paper describes a highly flexible computerized method of calculating operating data in a power cycle. The computerized method presented here permits the study of steam, gas and combined plants. Its flexibility is not restricted by any defined cycle scheme. A power plant consists of simple elements (turbine, compressor, combustor chamber, pump, etc.). Each power plant component is represented by its typical equations relating to fundamental mechanical and thermodynamic laws, so a power plant system is represented by algebraic equations, which are the typical equations of components, continuity equations, and data concerning plant conditions. This equation system is not linear, but can be reduced to a linear equation system with variable coefficients. The solution is simultaneous for each component and it is determined by an iterative process. An example of a simple gas turbine cycle demonstrates the applied technique. This paper also presents the user interface based on MS-Windows. The input data, the results, and any characteristic parameters of a complex cycle scheme are also shown.
Numerical methods for analyzing electromagnetic scattering
NASA Technical Reports Server (NTRS)
Lee, S. W.; Lo, Y. T.; Chuang, S. L.; Lee, C. S.
1985-01-01
Attenuation properties of the normal modes in an overmoded waveguide coated with a lossy material were analyzed. It is found that the low-order modes, can be significantly attenuated even with a thin layer of coating if the coating material is not too lossy. A thinner layer of coating is required for large attenuation of the low-order modes if the coating material is magnetic rather than dielectric. The Radar Cross Section (RCS) from an uncoated circular guide terminated by a perfect electric conductor was calculated and compared with available experimental data. It is confirmed that the interior irradiation contributes to the RCS. The equivalent-current method based on the geometrical theory of diffraction (GTD) was chosen for the calculation of the contribution from the rim diffraction. The RCS reduction from a coated circular guide terminated by a PEC are planned schemes for the experiments are included. The waveguide coated with a lossy magnetic material is suggested as a substitute for the corrugated waveguide.
Advanced analysis methods in particle physics
Bhat, Pushpalatha C.; /Fermilab
2010-10-01
Each generation of high energy physics experiments is grander in scale than the previous - more powerful, more complex and more demanding in terms of data handling and analysis. The spectacular performance of the Tevatron and the beginning of operations of the Large Hadron Collider, have placed us at the threshold of a new era in particle physics. The discovery of the Higgs boson or another agent of electroweak symmetry breaking and evidence of new physics may be just around the corner. The greatest challenge in these pursuits is to extract the extremely rare signals, if any, from huge backgrounds arising from known physics processes. The use of advanced analysis techniques is crucial in achieving this goal. In this review, I discuss the concepts of optimal analysis, some important advanced analysis methods and a few examples. The judicious use of these advanced methods should enable new discoveries and produce results with better precision, robustness and clarity.
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
Cost estimating methods for advanced space systems
NASA Technical Reports Server (NTRS)
Cyr, Kelley
1988-01-01
The development of parametric cost estimating methods for advanced space systems in the conceptual design phase is discussed. The process of identifying variables which drive cost and the relationship between weight and cost are discussed. A theoretical model of cost is developed and tested using a historical data base of research and development projects.
An advanced probabilistic structural analysis method for implicit performance functions
NASA Technical Reports Server (NTRS)
Wu, Y.-T.; Millwater, H. R.; Cruse, T. A.
1989-01-01
In probabilistic structural analysis, the performance or response functions usually are implicitly defined and must be solved by numerical analysis methods such as finite element methods. In such cases, the most commonly used probabilistic analysis tool is the mean-based, second-moment method which provides only the first two statistical moments. This paper presents a generalized advanced mean value (AMV) method which is capable of establishing the distributions to provide additional information for reliability design. The method requires slightly more computations than the second-moment method but is highly efficient relative to the other alternative methods. In particular, the examples show that the AMV method can be used to solve problems involving non-monotonic functions that result in truncated distributions.
Hybrid Numerical Methods for Multiscale Simulations of Subsurface Biogeochemical Processes
Scheibe, Timothy D.; Tartakovsky, Alexandre M.; Tartakovsky, Daniel M.; Redden, George D.; Meakin, Paul
2007-08-01
Many subsurface flow and transport problems of importance today involve coupled non-linear flow, transport, and reaction in media exhibiting complex heterogeneity. In particular, problems involving biological mediation of reactions fall into this class of problems. Recent experimental research has revealed important details about the physical, chemical, and biological mechanisms involved in these processes at a variety of scales ranging from molecular to laboratory scales. However, it has not been practical or possible to translate detailed knowledge at small scales into reliable predictions of field-scale phenomena important for environmental management applications. A large assortment of numerical simulation tools have been developed, each with its own characteristic scale including molecular (e.g., molecular dynamics), microbial (e.g., cellular automata or particle individual-based models), pore (e.g., lattice-Boltzmann, pore network models, and discrete particle methods such as smoothed particle hydrodynamics) and continuum scales (e.g., traditional partial differential equations solved by finite difference or finite element methods). While many problems can be effectively addressed by one of these models at a single scale, some problems may require explicit integration of models across multiple scales. We are developing a hybrid multi-scale subsurface reactive transport modeling framework that integrates models with diverse representations of physics, chemistry and biology at different scales (sub-pore, pore and continuum). The modeling framework is being designed to take advantage of advanced computational technologies including parallel code components using the Common Component Architecture, parallel solvers, gridding, data and workflow management, and visualization. This paper describes the specific methods/codes being used at each scale, techniques used to directly and adaptively couple across model scales, and preliminary results of application to a
A Novel Numerical Method for Fuzzy Boundary Value Problems
NASA Astrophysics Data System (ADS)
Can, E.; Bayrak, M. A.; Hicdurmaz
2016-05-01
In the present paper, a new numerical method is proposed for solving fuzzy differential equations which are utilized for the modeling problems in science and engineering. Fuzzy approach is selected due to its important applications on processing uncertainty or subjective information for mathematical models of physical problems. A second-order fuzzy linear boundary value problem is considered in particular due to its important applications in physics. Moreover, numerical experiments are presented to show the effectiveness of the proposed numerical method on specific physical problems such as heat conduction in an infinite plate and a fin.
Asymptotic-induced numerical methods for conservation laws
NASA Technical Reports Server (NTRS)
Garbey, Marc; Scroggs, Jeffrey S.
1990-01-01
Asymptotic-induced methods are presented for the numerical solution of hyperbolic conservation laws with or without viscosity. The methods consist of multiple stages. The first stage is to obtain a first approximation by using a first-order method, such as the Godunov scheme. Subsequent stages of the method involve solving internal-layer problems identified by using techniques derived via asymptotics. Finally, a residual correction increases the accuracy of the scheme. The method is derived and justified with singular perturbation techniques.
Advanced electromagnetic methods for aerospace vehicles
NASA Technical Reports Server (NTRS)
Balanis, Constantine A.; El-Sharawy, El-Budawy; Hashemi-Yeganeh, Shahrokh; Aberle, James T.; Birtcher, Craig R.
1991-01-01
The Advanced Helicopter Electromagnetics is centered on issues that advance technology related to helicopter electromagnetics. Progress was made on three major topics: composite materials; precipitation static corona discharge; and antenna technology. In composite materials, the research has focused on the measurements of their electrical properties, and the modeling of material discontinuities and their effect on the radiation pattern of antennas mounted on or near material surfaces. The electrical properties were used to model antenna performance when mounted on composite materials. Since helicopter platforms include several antenna systems at VHF and UHF bands, measuring techniques are being explored that can be used to measure the properties at these bands. The effort on corona discharge and precipitation static was directed toward the development of a new two dimensional Voltage Finite Difference Time Domain computer program. Results indicate the feasibility of using potentials for simulating electromagnetic problems in the cases where potentials become primary sources. In antenna technology the focus was on Polarization Diverse Conformal Microstrip Antennas, Cavity Backed Slot Antennas, and Varactor Tuned Circular Patch Antennas. Numerical codes were developed for the analysis of two probe fed rectangular and circular microstrip patch antennas fed by resistive and reactive power divider networks.
Parallel processing numerical method for confined vortex dynamics and applications
NASA Astrophysics Data System (ADS)
Bistrian, Diana Alina
2013-10-01
This paper explores a combined analytical and numerical technique to investigate the hydrodynamic instability of confined swirling flows, with application to vortex rope dynamics in a Francis turbine diffuser, in condition of sophisticated boundary constraints. We present a new approach based on the method of orthogonal decomposition in the Hilbert space, implemented with a spectral descriptor scheme in discrete space. A parallel implementation of the numerical scheme is conducted reducing the computational time compared to other techniques.
Numerical design of advanced multi-element airfoils
NASA Technical Reports Server (NTRS)
Mathias, Donovan L.; Cummings, Russell M.
1994-01-01
The current study extends the application of computational fluid dynamics to three-dimensional high-lift systems. Structured, overset grids are used in conjunction with an incompressible Navier-Stokes flow solver to investigate flow over a two-element high-lift configuration. The computations were run in a fully turbulent mode using the one-equation Baldwin-Barth turbulence model. The geometry consisted of an unswept wing which spanned a wind tunnel test section. Flows over full and half-span Fowler flap configurations were computed. Grid resolution issues were investigated in two dimensional studies of the flapped airfoil. Results of the full-span flap wing agreed well with experimental data and verified the method. Flow over the wing with the half-span was computed to investigate the details of the flow at the free edge of the flap. The results illustrated changes in flow streamlines, separation locations, and surface pressures due to the vortex shed from the flap edge.
Investigating Convergence Patterns for Numerical Methods Using Data Analysis
ERIC Educational Resources Information Center
Gordon, Sheldon P.
2013-01-01
The article investigates the patterns that arise in the convergence of numerical methods, particularly those in the errors involved in successive iterations, using data analysis and curve fitting methods. In particular, the results obtained are used to convey a deeper level of understanding of the concepts of linear, quadratic, and cubic…
Collocation Method for Numerical Solution of Coupled Nonlinear Schroedinger Equation
Ismail, M. S.
2010-09-30
The coupled nonlinear Schroedinger equation models several interesting physical phenomena presents a model equation for optical fiber with linear birefringence. In this paper we use collocation method to solve this equation, we test this method for stability and accuracy. Numerical tests using single soliton and interaction of three solitons are used to test the resulting scheme.
A numerical method for solving singular De`s
Mahaver, W.T.
1996-12-31
A numerical method is developed for solving singular differential equations using steepest descent based on weighted Sobolev gradients. The method is demonstrated on a variety of first and second order problems, including linear constrained, unconstrained, and partially constrained first order problems, a nonlinear first order problem with irregular singularity, and two second order variational problems.
A new numerical method of total solar eclipse photography processing
NASA Astrophysics Data System (ADS)
Druckmüller, M.; Rušin, V.; Minarovjech, M.
2006-10-01
A new numerical method of image processing suitable for visualization of corona images taken during total solar eclipses is presented. This method allows us to study both small- and large-scale coronal structures that remain invisible on original images because of their very high dynamic range of the coronal brightness. The method is based on the use of adaptive filters inspired by human vision and the sensitivity of resulting images is thus very close to that of the human eye during an eclipse. A high precision alignment method for white-light corona images is also discussed. The proposed method highly improves a widely used unsharp masking method employing a radially blurred mask. The results of these numerical image processing techniques are illustrated by a series of images taken during eclipses of the last decade. The method minimizes the risk of processing artifacts.
Comparison of methods for numerical calculation of continuum damping
Bowden, G. W.; Hole, M. J.; Dennis, G. R.; Könies, A.; Gorelenkov, N. N.
2014-05-15
Continuum resonance damping is an important factor in determining the stability of certain global modes in fusion plasmas. A number of analytic and numerical approaches have been developed to compute this damping, particularly, in the case of the toroidicity-induced shear Alfvén eigenmode. This paper compares results obtained using an analytical perturbative approach with those found using resistive and complex contour numerical approaches. It is found that the perturbative method does not provide accurate agreement with reliable numerical methods for the range of parameters examined. This discrepancy exists even in the limit where damping approaches zero. When the perturbative technique is implemented using a standard finite element method, the damping estimate fails to converge with radial grid resolution. The finite elements used cannot accurately represent the eigenmode in the region of the continuum resonance, regardless of the number of radial grid points used.
Advanced Bayesian Method for Planetary Surface Navigation
NASA Technical Reports Server (NTRS)
Center, Julian
2015-01-01
Autonomous Exploration, Inc., has developed an advanced Bayesian statistical inference method that leverages current computing technology to produce a highly accurate surface navigation system. The method combines dense stereo vision and high-speed optical flow to implement visual odometry (VO) to track faster rover movements. The Bayesian VO technique improves performance by using all image information rather than corner features only. The method determines what can be learned from each image pixel and weighs the information accordingly. This capability improves performance in shadowed areas that yield only low-contrast images. The error characteristics of the visual processing are complementary to those of a low-cost inertial measurement unit (IMU), so the combination of the two capabilities provides highly accurate navigation. The method increases NASA mission productivity by enabling faster rover speed and accuracy. On Earth, the technology will permit operation of robots and autonomous vehicles in areas where the Global Positioning System (GPS) is degraded or unavailable.
NASA Astrophysics Data System (ADS)
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.
Advances in numerical solutions to integral equations in liquid state theory
NASA Astrophysics Data System (ADS)
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
Transonic wing analysis using advanced computational methods
NASA Technical Reports Server (NTRS)
Henne, P. A.; Hicks, R. M.
1978-01-01
This paper discusses the application of three-dimensional computational transonic flow methods to several different types of transport wing designs. The purpose of these applications is to evaluate the basic accuracy and limitations associated with such numerical methods. The use of such computational methods for practical engineering problems can only be justified after favorable evaluations are completed. The paper summarizes a study of both the small-disturbance and the full potential technique for computing three-dimensional transonic flows. Computed three-dimensional results are compared to both experimental measurements and theoretical results. Comparisons are made not only of pressure distributions but also of lift and drag forces. Transonic drag rise characteristics are compared. Three-dimensional pressure distributions and aerodynamic forces, computed from the full potential solution, compare reasonably well with experimental results for a wide range of configurations and flow conditions.
Advanced Analysis Methods in High Energy Physics
Pushpalatha C. Bhat
2001-10-03
During the coming decade, high energy physics experiments at the Fermilab Tevatron and around the globe will use very sophisticated equipment to record unprecedented amounts of data in the hope of making major discoveries that may unravel some of Nature's deepest mysteries. The discovery of the Higgs boson and signals of new physics may be around the corner. The use of advanced analysis techniques will be crucial in achieving these goals. The author discusses some of the novel methods of analysis that could prove to be particularly valuable for finding evidence of any new physics, for improving precision measurements and for exploring parameter spaces of theoretical models.
Numerical modeling of magnetic induction tomography using the impedance method.
Ramos, Airton; Wolff, Julia G B
2011-02-01
This article discusses the impedance method in the forward calculation in magnetic induction tomography (MIT). Magnetic field and eddy current distributions were obtained numerically for a sphere in the field of a coil and were compared with an analytical model. Additionally, numerical and experimental results for phase sensitivity in MIT were obtained and compared for a cylindrical object in a planar array of sensors. The results showed that the impedance method provides results that agree very well with reality in the frequency range from 100 kHz to 20 MHz and for low conductivity objects (10 S/m or less). This opens the possibility of using this numerical approach in image reconstruction in MIT. PMID:21229327
Numerical Methods of Computational Electromagnetics for Complex Inhomogeneous Systems
Cai, Wei
2014-05-15
Understanding electromagnetic phenomena is the key in many scientific investigation and engineering designs such as solar cell designs, studying biological ion channels for diseases, and creating clean fusion energies, among other things. The objectives of the project are to develop high order numerical methods to simulate evanescent electromagnetic waves occurring in plasmon solar cells and biological ion-channels, where local field enhancement within random media in the former and long range electrostatic interactions in the latter are of major challenges for accurate and efficient numerical computations. We have accomplished these objectives by developing high order numerical methods for solving Maxwell equations such as high order finite element basis for discontinuous Galerkin methods, well-conditioned Nedelec edge element method, divergence free finite element basis for MHD, and fast integral equation methods for layered media. These methods can be used to model the complex local field enhancement in plasmon solar cells. On the other hand, to treat long range electrostatic interaction in ion channels, we have developed image charge based method for a hybrid model in combining atomistic electrostatics and continuum Poisson-Boltzmann electrostatics. Such a hybrid model will speed up the molecular dynamics simulation of transport in biological ion-channels.
COMPARING NUMERICAL METHODS FOR ISOTHERMAL MAGNETIZED SUPERSONIC TURBULENCE
Kritsuk, Alexei G.; Collins, David; Norman, Michael L.; Xu Hao E-mail: dccollins@lanl.gov
2011-08-10
Many astrophysical applications involve magnetized turbulent flows with shock waves. Ab initio star formation simulations require a robust representation of supersonic turbulence in molecular clouds on a wide range of scales imposing stringent demands on the quality of numerical algorithms. We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine popular astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. These applications employ a variety of numerical approaches, including both split and unsplit, finite difference and finite volume, divergence preserving and divergence cleaning, a variety of Riemann solvers, and a range of spatial reconstruction and time integration techniques. We present a comprehensive set of statistical measures designed to quantify the effects of numerical dissipation in these MHD solvers. We compare power spectra for basic fields to determine the effective spectral bandwidth of the methods and rank them based on their relative effective Reynolds numbers. We also compare numerical dissipation for solenoidal and dilatational velocity components to check for possible impacts of the numerics on small-scale density statistics. Finally, we discuss the convergence of various characteristics for the turbulence decay test and the impact of various components of numerical schemes on the accuracy of solutions. The nine codes gave qualitatively the same results, implying that they are all performing reasonably well and are useful for scientific applications. We show that the best performing codes employ a consistently high order of accuracy for spatial reconstruction of the evolved fields, transverse gradient interpolation, conservation law update step, and Lorentz force computation. The best results are achieved with divergence-free evolution of the
A gyrokinetic continuum code based on the numerical Lie transform (NLT) method
NASA Astrophysics Data System (ADS)
Ye, Lei; Xu, Yingfeng; Xiao, Xiaotao; Dai, Zongliang; Wang, Shaojie
2016-07-01
In this work, we report a novel gyrokinetic simulation method named numerical Lie transform (NLT), which depends on a new physical model derived from the I-transform theory. In this model, the perturbed motion of a particle is decoupled from the unperturbed motion. Due to this property, the unperturbed orbit can be computed in advance and saved as numerical tables for real-time computation. A 4D tensor B-spline interpolation module is developed and applied with the semi-Lagrangian scheme to avoid operator splitting. The NLT code is verified by the Rosenbluth-Hinton test and the linear ITG Cyclone test.
A novel gas-droplet numerical method for spray combustion
NASA Technical Reports Server (NTRS)
Chen, C. P.; Shang, H. M.; Jiang, Y.
1991-01-01
This paper presents a non-iterative numerical technique for computing time-dependent gas-droplet flows. The method is a fully-interacting combination of Eulerian fluid and Lagrangian particle calculation. The interaction calculations between the two phases are formulated on a pressure-velocity coupling procedure based on the operator-splitting technique. This procedure eliminates the global iterations required in the conventional particle-source-in-cell (PSIC) procedure. Turbulent dispersion calculations are treated by a stochastic procedure. Numerical calculations and comparisons with available experimental data, as well as efficiency assessments are given for some sprays typical of spray combustion applications.
Simple numerical method for predicting steady compressible flows
NASA Technical Reports Server (NTRS)
Vonlavante, Ernst; Nelson, N. Duane
1986-01-01
A numerical method for solving the isenthalpic form of the governing equations for compressible viscous and inviscid flows was developed. The method was based on the concept of flux vector splitting in its implicit form. The method was tested on several demanding inviscid and viscous configurations. Two different forms of the implicit operator were investigated. The time marching to steady state was accelerated by the implementation of the multigrid procedure. Its various forms very effectively increased the rate of convergence of the present scheme. High quality steady state results were obtained in most of the test cases; these required only short computational times due to the relative efficiency of the basic method.
Singularity Preserving Numerical Methods for Boundary Integral Equations
NASA Technical Reports Server (NTRS)
Kaneko, Hideaki (Principal Investigator)
1996-01-01
In the past twelve months (May 8, 1995 - May 8, 1996), under the cooperative agreement with Division of Multidisciplinary Optimization at NASA Langley, we have accomplished the following five projects: a note on the finite element method with singular basis functions; numerical quadrature for weakly singular integrals; superconvergence of degenerate kernel method; superconvergence of the iterated collocation method for Hammersteion equations; and singularity preserving Galerkin method for Hammerstein equations with logarithmic kernel. This final report consists of five papers describing these projects. Each project is preceeded by a brief abstract.
The development of accurate and efficient methods of numerical quadrature
NASA Technical Reports Server (NTRS)
Feagin, T.
1973-01-01
Some new methods for performing numerical quadrature of an integrable function over a finite interval are described. Each method provides a sequence of approximations of increasing order to the value of the integral. Each approximation makes use of all previously computed values of the integrand. The points at which new values of the integrand are computed are selected in such a way that the order of the approximation is maximized. The methods are compared with the quadrature methods of Clenshaw and Curtis, Gauss, Patterson, and Romberg using several examples.
A numerical method for interface problems in elastodynamics
NASA Technical Reports Server (NTRS)
Mcghee, D. S.
1984-01-01
The numerical implementation of a formulation for a class of interface problems in elastodynamics is discussed. This formulation combines the use of the finite element and boundary integral methods to represent the interior and the exteriro regions, respectively. In particular, the response of a semicylindrical alluvial valley in a homogeneous halfspace to incident antiplane SH waves is considered to determine the accuracy and convergence of the numerical procedure. Numerical results are obtained from several combinations of the incidence angle, frequency of excitation, and relative stiffness between the inclusion and the surrounding halfspace. The results tend to confirm the theoretical estimates that the convergence is of the order H(2) for the piecewise linear elements used. It was also observed that the accuracy descreases as the frequency of excitation increases or as the relative stiffness of the inclusion decreases.
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.
Numerical simulation methods for the Rouse model in flow
NASA Astrophysics Data System (ADS)
Howard, Michael P.; Milner, Scott T.
2011-11-01
Simulation of the Rouse model in flow underlies a great variety of numerical investigations of polymer dynamics, in both entangled melts and solutions and in dilute solution. Typically a simple explicit stochastic Euler method is used to evolve the Rouse model. Here we compare this approach to an operator splitting method, which splits the evolution operator into stochastic linear and deterministic nonlinear parts and takes advantage of an analytical solution for the linear Rouse model in terms of the noise history. We show that this splitting method has second-order weak convergence, whereas the Euler method has only first-order weak convergence. Furthermore, the splitting method is unconditionally stable, in contrast to the limited stability range of the Euler method. Similar splitting methods are applicable to a broad class of problems in stochastic dynamics in which noise competes with ordering and flow to determine steady-state order parameter structures.
Advancing MODFLOW Applying the Derived Vector Space Method
NASA Astrophysics Data System (ADS)
Herrera, G. S.; Herrera, I.; Lemus-García, M.; Hernandez-Garcia, G. D.
2015-12-01
The most effective domain decomposition methods (DDM) are non-overlapping DDMs. Recently a new approach, the DVS-framework, based on an innovative discretization method that uses a non-overlapping system of nodes (the derived-nodes), was introduced and developed by I. Herrera et al. [1, 2]. Using the DVS-approach a group of four algorithms, referred to as the 'DVS-algorithms', which fulfill the DDM-paradigm (i.e. the solution of global problems is obtained by resolution of local problems exclusively) has been derived. Such procedures are applicable to any boundary-value problem, or system of such equations, for which a standard discretization method is available and then software with a high degree of parallelization can be constructed. In a parallel talk, in this AGU Fall Meeting, Ismael Herrera will introduce the general DVS methodology. The application of the DVS-algorithms has been demonstrated in the solution of several boundary values problems of interest in Geophysics. Numerical examples for a single-equation, for the cases of symmetric, non-symmetric and indefinite problems were demonstrated before [1,2]. For these problems DVS-algorithms exhibited significantly improved numerical performance with respect to standard versions of DDM algorithms. In view of these results our research group is in the process of applying the DVS method to a widely used simulator for the first time, here we present the advances of the application of this method for the parallelization of MODFLOW. Efficiency results for a group of tests will be presented. References [1] I. Herrera, L.M. de la Cruz and A. Rosas-Medina. Non overlapping discretization methods for partial differential equations, Numer Meth Part D E, (2013). [2] Herrera, I., & Contreras Iván "An Innovative Tool for Effectively Applying Highly Parallelized Software To Problems of Elasticity". Geofísica Internacional, 2015 (In press)
Numerical Polynomial Homotopy Continuation Method and String Vacua
Mehta, Dhagash
2011-01-01
Finding vmore » acua for the four-dimensional effective theories for supergravity which descend from flux compactifications and analyzing them according to their stability is one of the central problems in string phenomenology. Except for some simple toy models, it is, however, difficult to find all the vacua analytically. Recently developed algorithmic methods based on symbolic computer algebra can be of great help in the more realistic models. However, they suffer from serious algorithmic complexities and are limited to small system sizes. In this paper, we review a numerical method called the numerical polynomial homotopy continuation (NPHC) method, first used in the areas of lattice field theories, which by construction finds all of the vacua of a given potential that is known to have only isolated solutions. The NPHC method is known to suffer from no major algorithmic complexities and is embarrassingly parallelizable , and hence its applicability goes way beyond the existing symbolic methods. We first solve a simple toy model as a warm-up example to demonstrate the NPHC method at work. We then show that all the vacua of a more complicated model of a compactified M theory model, which has an S U ( 3 ) structure, can be obtained by using a desktop machine in just about an hour, a feat which was reported to be prohibitively difficult by the existing symbolic methods. Finally, we compare the various technicalities between the two methods.« less
Numerical methods for characterization of synchrotron radiation based on the Wigner function method
NASA Astrophysics Data System (ADS)
Tanaka, Takashi
2014-06-01
Numerical characterization of synchrotron radiation based on the Wigner function method is explored in order to accurately evaluate the light source performance. A number of numerical methods to compute the Wigner functions for typical synchrotron radiation sources such as bending magnets, undulators and wigglers, are presented, which significantly improve the computation efficiency and reduce the total computation time. As a practical example of the numerical characterization, optimization of betatron functions to maximize the brilliance of undulator radiation is discussed.
Projected discrete ordinates methods for numerical transport problems
Larsen, E.W.
1985-01-01
A class of Projected Discrete-Ordinates (PDO) methods is described for obtaining iterative solutions of discrete-ordinates problems with convergence rates comparable to those observed using Diffusion Synthetic Acceleration (DSA). The spatially discretized PDO solutions are generally not equal to the DSA solutions, but unlike DSA, which requires great care in the use of spatial discretizations to preserve stability, the PDO solutions remain stable and rapidly convergent with essentially arbitrary spatial discretizations. Numerical results are presented which illustrate the rapid convergence and the accuracy of solutions obtained using PDO methods with commonplace differencing methods.
Computational methods for aerodynamic design using numerical optimization
NASA Technical Reports Server (NTRS)
Peeters, M. F.
1983-01-01
Five methods to increase the computational efficiency of aerodynamic design using numerical optimization, by reducing the computer time required to perform gradient calculations, are examined. The most promising method consists of drastically reducing the size of the computational domain on which aerodynamic calculations are made during gradient calculations. Since a gradient calculation requires the solution of the flow about an airfoil whose geometry was slightly perturbed from a base airfoil, the flow about the base airfoil is used to determine boundary conditions on the reduced computational domain. This method worked well in subcritical flow.
Fast and stable numerical method for neuronal modelling
NASA Astrophysics Data System (ADS)
Hashemi, Soheil; Abdolali, Ali
2016-11-01
Excitable cell modelling is of a prime interest in predicting and targeting neural activity. Two main limits in solving related equations are speed and stability of numerical method. Since there is a tradeoff between accuracy and speed, most previously presented methods for solving partial differential equations (PDE) are focused on one side. More speed means more accurate simulations and therefore better device designing. By considering the variables in finite differenced equation in proper time and calculating the unknowns in the specific sequence, a fast, stable and accurate method is introduced in this paper for solving neural partial differential equations. Propagation of action potential in giant axon is studied by proposed method and traditional methods. Speed, consistency and stability of the methods are compared and discussed. The proposed method is as fast as forward methods and as stable as backward methods. Forward methods are known as fastest methods and backward methods are stable in any circumstances. Complex structures can be simulated by proposed method due to speed and stability of the method.
A numerical method to study the dynamics of capillary fluid systems
NASA Astrophysics Data System (ADS)
Herrada, M. A.; Montanero, J. M.
2016-02-01
We propose a numerical approach to study both the nonlinear dynamics and linear stability of capillary fluid systems. In the nonlinear analysis, the time-dependent fluid region is mapped onto a fixed numerical domain through a coordinate transformation. The hydrodynamic equations are spatially discretized with the Chebyshev spectral collocation technique, while an implicit time advancement is performed using second-order backward finite differences. The resulting algebraic equations are solved with the iterative Newton-Raphson technique. The most novel aspect of the method is the fact that the elements of the Jacobian of the discretized system of equations are symbolic functions calculated before running the simulation. These functions are evaluated numerically in the Newton-Raphson iterations to find the solution at each time step, which reduces considerably the computing time. Besides, this numerical procedure can be easily adapted to solve the eigenvalue problem which determines the linear global modes of the capillary system. Therefore, both the nonlinear dynamics and the linear stability analysis can be conducted with essentially the same algorithm. We validate this numerical approach by studying the dynamics of a liquid bridge close to its minimum volume stability limit. The results are virtually the same as those obtained with other methods. The proposed approach proves to be much more computationally efficient than those other methods. Finally, we show the versatility of the method by calculating the linear global modes of a gravitational jet.
Numerical method for shear bands in ductile metal with inclusions
Plohr, Jee Yeon N; Plohr, Bradley J
2010-01-01
A numerical method for mesoscale simulation of high strain-rate loading of ductile metal containing inclusions is described. Because of small-scale inhomogeneities, such a composite material is prone to localized shear deformation (adiabatic shear bands). The modeling framework is the Generalized Method of Cells of Paley and Aboudi [Mech. Materials, vol. 14, pp. /27-139, 1992], which ensures that the micromechanical response of the material is reflected in the behavior of the composite at the mesoscale. To calculate the effective plastic strain rate when shear bands are present, the analytic and numerical analysis of shear bands by Glimm, Plohr, and Sharp [Mech. Materials, vol. 24, pp. 31-41, 1996] is adapted and extended.
Numerical Method for the Astronomical Almanac and Orbit Calculations
NASA Astrophysics Data System (ADS)
Kim, Kap-Sung
1993-12-01
We have calculated the astronomical almanac 1994 and simulated the trajectory of a satellite orbit considering all perturbative forces with various initial conditions. In this work, Gauss Jackson multistep integration method has been used to calculate our basic equation of motion with high numerical accuracy. It has been found that our results agree well with the Astronomical Almanac Data distributed by JPL of NASA and the orbit simulations have been carried out with fast speed, stability and excellent round-off error accumulation, comparing with other numerical methods. In order to be carried out our works on almanac and orbit calculations easily by anyone who uses a personal computer, we have made a computer program on graphical user interface to provide various menus for detail works selected by a mouse.
Asymptotic and Numerical Methods for Rapidly Rotating Buoyant Flow
NASA Astrophysics Data System (ADS)
Grooms, Ian G.
This thesis documents three investigations carried out in pursuance of a doctoral degree in applied mathematics at the University of Colorado (Boulder). The first investigation concerns the properties of rotating Rayleigh-Benard convection -- thermal convection in a rotating infinite plane layer between two constant-temperature boundaries. It is noted that in certain parameter regimes convective Taylor columns appear which dominate the dynamics, and a semi-analytical model of these is presented. Investigation of the columns and of various other properties of the flow is ongoing. The second investigation concerns the interactions between planetary-scale and mesoscale dynamics in the oceans. Using multiple-scale asymptotics the possible connections between planetary geostrophic and quasigeostrophic dynamics are investigated, and three different systems of coupled equations are derived. Possible use of these equations in conjunction with the method of superparameterization, and extension of the asymptotic methods to the interactions between mesoscale and submesoscale dynamics is ongoing. The third investigation concerns the linear stability properties of semi-implicit methods for the numerical integration of ordinary differential equations, focusing in particular on the linear stability of IMEX (Implicit-Explicit) methods and exponential integrators applied to systems of ordinary differential equations arising in the numerical solution of spatially discretized nonlinear partial differential equations containing both dispersive and dissipative linear terms. While these investigations may seem unrelated at first glance, some reflection shows that they are in fact closely linked. The investigation of rotating convection makes use of single-space, multiple-time-scale asymptotics to deal with dynamics strongly constrained by rotation. Although the context of thermal convection in an infinite layer seems somewhat removed from large-scale ocean dynamics, the asymptotic
Numerical optimization methods for critical currents in superconductors
NASA Astrophysics Data System (ADS)
Kimmel, Gregory; Sadovskyy, Ivan; Koshelev, Alex; Glatz, Andreas
In this work, I present optimization methods for maximizing the critical current in high-temperature superconductors for energy applications. The critical current in the presence of an external magnetic field is mostly defined by the pinning landscape (pinscape) within the superconductor, which prevents magnetic vortices from moving and, therefore, increases its critical current. Our approach is to generate different pinscapes and obtain the resulting critical current by large-scale time-dependent Ginzburg-Landau equations. Pinning centers could be any combination of defects, including spherical and columnar defects. The parameters controlling the pinscape are adaptively adjusted in order to find the optimal parameter set, which maximizes the critical current. Here, we compare different optimization methods and discuss their performance. Work was supported by the Scientific Discovery through Advanced Computing (SciDAC) program funded by U.S. Department of Energy, Office of Science, Advanced Scientific Computing Research and Basic Energy Sciences.
An alternative numerical method for the stationary pulsar magnetosphere
NASA Astrophysics Data System (ADS)
Takamori, Yohsuke; Okawa, Hirotada; Takamoto, Makoto; Suwa, Yudai
2014-02-01
Stationary pulsar magnetospheres in the force-free system are governed by the pulsar equation. In 1999, Contopoulos, Kazanas, and Fendt (hereafter CKF) numerically solved the pulsar equation and obtained a pulsar magnetosphere model called the CKF solution that has both closed and open magnetic field lines. The CKF solution is a successful solution, but it contains a poloidal current sheet that flows along the last open field line. This current sheet is artificially added to make the current system closed. In this paper, we suggest an alternative method to solve the pulsar equation and construct pulsar magnetosphere models without a current sheet. In our method, the pulsar equation is decomposed into Ampère's law and the force-free condition. We numerically solve these equations simultaneously with a fixed poloidal current. As a result, we obtain a pulsar magnetosphere model without a current sheet, which is similar to the CKF solution near the neutron star and has a jet-like structure at a distance along the pole. In addition, we discuss physical properties of the model and find that the force-free condition breaks down in a vicinity of the light cylinder due to dissipation that is included implicitly in the numerical method.
[Numerical methods for multi-fluid flows]. Final progress report
Pozrikidis, C.
1998-07-21
The central objective of this research has been to develop efficient numerical methods for computing multi-fluid flows with large interfacial deformations, and apply these methods to study the rheology of suspensions of deformable particles with viscous and non-Newtonian interfacial behavior. The mathematical formulation employs boundary-integral, immersed-boundary, and related numerical methods. Particles of interest include liquid drops with constant surface tension and capsules whose interfaces exhibit viscoelastic and incompressible characteristics. In one family of problems, the author has considered the shear-driven and pressure-driven flow of a suspension of two-dimensional liquid drops with ordered and random structure. In a second series of investigations, the author carried out dynamic simulations of two-dimensional, unbounded, doubly-periodic shear flows with random structure. Another family of problems addresses the deformation of three-dimensional capsules whose interfaces exhibit isotropic surface tension, viscous, elastic, or incompressible behavior, in simple shear flow. The numerical results extend previous asymptotic theories for small deformations and illuminate the mechanism of membrane rupture.
Numerical method for wave forces acting on partially perforated caisson
NASA Astrophysics Data System (ADS)
Jiang, Feng; Tang, Xiao-cheng; Jin, Zhao; Zhang, Li; Chen, Hong-zhou
2015-04-01
The perforated caisson is widely applied to practical engineering because of its great advantages in effectively wave energy consumption and cost reduction. The attentions of many scientists were paid to the fluid-structure interaction between wave and perforated caisson studies, but until now, most concerns have been put on theoretical analysis and experimental model set up. In this paper, interaction between the wave and the partial perforated caisson in a 2D numerical wave flume is investigated by means of the renewed SPH algorithm, and the mathematical equations are in the form of SPH numerical approximation based on Navier-Stokes equations. The validity of the SPH mathematical method is examined and the simulated results are compared with the results of theoretical models, meanwhile the complex hydrodynamic characteristics when the water particles flow in or out of a wave absorbing chamber are analyzed and the wave pressure distribution of the perforated caisson is also addressed here. The relationship between the ratio of total horizontal force acting on caisson under regular waves and its influence factors is examined. The data show that the numerical calculation of the ratio of total horizontal force meets the empirical regression equation very well. The simulations of SPH about the wave nonlinearity and breaking are briefly depicted in the paper, suggesting that the advantages and great potentiality of the SPH method is significant compared with traditional methods.
The instanton method and its numerical implementation in fluid mechanics
NASA Astrophysics Data System (ADS)
Grafke, Tobias; Grauer, Rainer; Schäfer, Tobias
2015-08-01
A precise characterization of structures occurring in turbulent fluid flows at high Reynolds numbers is one of the last open problems of classical physics. In this review we discuss recent developments related to the application of instanton methods to turbulence. Instantons are saddle point configurations of the underlying path integrals. They are equivalent to minimizers of the related Freidlin-Wentzell action and known to be able to characterize rare events in such systems. While there is an impressive body of work concerning their analytical description, this review focuses on the question on how to compute these minimizers numerically. In a short introduction we present the relevant mathematical and physical background before we discuss the stochastic Burgers equation in detail. We present algorithms to compute instantons numerically by an efficient solution of the corresponding Euler-Lagrange equations. A second focus is the discussion of a recently developed numerical filtering technique that allows to extract instantons from direct numerical simulations. In the following we present modifications of the algorithms to make them efficient when applied to two- or three-dimensional (2D or 3D) fluid dynamical problems. We illustrate these ideas using the 2D Burgers equation and the 3D Navier-Stokes equations.
Dielectric Boundary Forces in Numerical Poisson-Boltzmann Methods: Theory and Numerical Strategies.
Cai, Qin; Ye, Xiang; Wang, Jun; Luo, Ray
2011-10-01
Continuum modeling of electrostatic interactions based upon the numerical solutions of the Poisson-Boltzmann equation has been widely adopted in biomolecular applications. To extend their applications to molecular dynamics and energy minimization, robust and efficient methodologies to compute solvation forces must be developed. In this study, we have first reviewed the theory for the computation of dielectric boundary forces based on the definition of the Maxwell stress tensor. This is followed by a new formulation of the dielectric boundary force suitable for the finite-difference Poisson-Boltzmann methods. We have validated the new formulation with idealized analytical systems and realistic molecular systems. PMID:22125339
Dielectric boundary force in numerical Poisson-Boltzmann methods: Theory and numerical strategies
NASA Astrophysics Data System (ADS)
Cai, Qin; Ye, Xiang; Wang, Jun; Luo, Ray
2011-10-01
Continuum modeling of electrostatic interactions based upon the numerical solutions of the Poisson-Boltzmann equation has been widely adopted in biomolecular applications. To extend their applications to molecular dynamics and energy minimization, robust and efficient methodologies to compute solvation forces must be developed. In this study, we have first reviewed the theory for the computation of dielectric boundary force based on the definition of the Maxwell stress tensor. This is followed by a new formulation of the dielectric boundary force suitable for the finite-difference Poisson-Boltzmann methods. We have validated the new formulation with idealized analytical systems and realistic molecular systems.
Numerical Method for Darcy Flow Derived Using Discrete Exterior Calculus
NASA Astrophysics Data System (ADS)
Hirani, A. N.; Nakshatrala, K. B.; Chaudhry, J. H.
2015-05-01
We derive a numerical method for Darcy flow, and also for Poisson's equation in mixed (first order) form, based on discrete exterior calculus (DEC). Exterior calculus is a generalization of vector calculus to smooth manifolds and DEC is one of its discretizations on simplicial complexes such as triangle and tetrahedral meshes. DEC is a coordinate invariant discretization, in that it does not depend on the embedding of the simplices or the whole mesh. We start by rewriting the governing equations of Darcy flow using the language of exterior calculus. This yields a formulation in terms of flux differential form and pressure. The numerical method is then derived by using the framework provided by DEC for discretizing differential forms and operators that act on forms. We also develop a discretization for a spatially dependent Hodge star that varies with the permeability of the medium. This also allows us to address discontinuous permeability. The matrix representation for our discrete non-homogeneous Hodge star is diagonal, with positive diagonal entries. The resulting linear system of equations for flux and pressure are saddle type, with a diagonal matrix as the top left block. The performance of the proposed numerical method is illustrated on many standard test problems. These include patch tests in two and three dimensions, comparison with analytically known solutions in two dimensions, layered medium with alternating permeability values, and a test with a change in permeability along the flow direction. We also show numerical evidence of convergence of the flux and the pressure. A convergence experiment is included for Darcy flow on a surface. A short introduction to the relevant parts of smooth and discrete exterior calculus is included in this article. We also include a discussion of the boundary condition in terms of exterior calculus.
Numerical methods for scattering from electrically large objects
NASA Astrophysics Data System (ADS)
Enguist, Bjorn; Murphy, W. D.; Rokhlin, Vladimir; Vassiliou, Marius S.
1991-05-01
A new and computationally very efficient integral equation numerical method for computing electromagnetic scattering and radar cross section (RCS) was developed. A theory of higher order impedance boundary conditions was derived to handle single and multiple dielectric coatings around conductors. The method was tested in two dimensions using a 14,000-line FORTRAN program and was found to be very promising for electrically large objects. Initial ideas for extensions to three dimensions were explored. Treatments of trailing edge and corner singularities were developed.
Teaching Thermal Hydraulics & Numerical Methods: An Introductory Control Volume Primer
Lucas, D.S.
2004-10-03
This paper covers the basics of the implementation of the control volume method in the context of the Homogeneous Equilibrium Model (HEM)(T/H) code using the conservation equations of mass, momentum, and energy. This primer uses the advection equation as a template. The discussion will cover the basic equations of the control volume portion of the course in the primer, which includes the advection equation, numerical methods, along with the implementation of the various equations via FORTRAN into computer programs and the final result for a three equation HEM code and its validation.
Integrated numerical methods for hypersonic aircraft cooling systems analysis
NASA Technical Reports Server (NTRS)
Petley, Dennis H.; Jones, Stuart C.; Dziedzic, William M.
1992-01-01
Numerical methods have been developed for the analysis of hypersonic aircraft cooling systems. A general purpose finite difference thermal analysis code is used to determine areas which must be cooled. Complex cooling networks of series and parallel flow can be analyzed using a finite difference computer program. Both internal fluid flow and heat transfer are analyzed, because increased heat flow causes a decrease in the flow of the coolant. The steady state solution is a successive point iterative method. The transient analysis uses implicit forward-backward differencing. Several examples of the use of the program in studies of hypersonic aircraft and rockets are provided.
Optimization methods and silicon solar cell numerical models
NASA Technical Reports Server (NTRS)
Girardini, K.; Jacobsen, S. E.
1986-01-01
An optimization algorithm for use with numerical silicon solar cell models was developed. By coupling an optimization algorithm with a solar cell model, it is possible to simultaneously vary design variables such as impurity concentrations, front junction depth, back junction depth, and cell thickness to maximize the predicted cell efficiency. An optimization algorithm was developed and interfaced with the Solar Cell Analysis Program in 1 Dimension (SCAP1D). SCAP1D uses finite difference methods to solve the differential equations which, along with several relations from the physics of semiconductors, describe mathematically the performance of a solar cell. A major obstacle is that the numerical methods used in SCAP1D require a significant amount of computer time, and during an optimization the model is called iteratively until the design variables converge to the values associated with the maximum efficiency. This problem was alleviated by designing an optimization code specifically for use with numerically intensive simulations, to reduce the number of times the efficiency has to be calculated to achieve convergence to the optimal solution.
Numerical methods for the Poisson-Fermi equation in electrolytes
NASA Astrophysics Data System (ADS)
Liu, Jinn-Liang
2013-08-01
The Poisson-Fermi equation proposed by Bazant, Storey, and Kornyshev [Phys. Rev. Lett. 106 (2011) 046102] for ionic liquids is applied to and numerically studied for electrolytes and biological ion channels in three-dimensional space. This is a fourth-order nonlinear PDE that deals with both steric and correlation effects of all ions and solvent molecules involved in a model system. The Fermi distribution follows from classical lattice models of configurational entropy of finite size ions and solvent molecules and hence prevents the long and outstanding problem of unphysical divergence predicted by the Gouy-Chapman model at large potentials due to the Boltzmann distribution of point charges. The equation reduces to Poisson-Boltzmann if the correlation length vanishes. A simplified matched interface and boundary method exhibiting optimal convergence is first developed for this equation by using a gramicidin A channel model that illustrates challenging issues associated with the geometric singularities of molecular surfaces of channel proteins in realistic 3D simulations. Various numerical methods then follow to tackle a range of numerical problems concerning the fourth-order term, nonlinearity, stability, efficiency, and effectiveness. The most significant feature of the Poisson-Fermi equation, namely, its inclusion of steric and correlation effects, is demonstrated by showing good agreement with Monte Carlo simulation data for a charged wall model and an L type calcium channel model.
Numerical integration of population models satisfying conservation laws: NSFD methods.
Mickens, Ronald E
2007-10-01
Population models arising in ecology, epidemiology and mathematical biology may involve a conservation law, i.e. the total population is constant. In addition to these cases, other situations may occur for which the total population, asymptotically in time, approach a constant value. Since it is rarely the situation that the equations of motion can be analytically solved to obtain exact solutions, it follows that numerical techniques are needed to provide solutions. However, numerical procedures are only valid if they can reproduce fundamental properties of the differential equations modeling the phenomena of interest. We show that for population models, involving a dynamical conservation law the use of nonstandard finite difference (NSFD) methods allows the construction of discretization schemes such that they are dynamically consistent (DC) with the original differential equations. The paper will briefly discuss the NSFD methodology, the concept of DC, and illustrate their application to specific problems for population models. PMID:22876826
7 CFR 27.92 - Method of payment; advance deposit.
Code of Federal Regulations, 2010 CFR
2010-01-01
... 7 Agriculture 2 2010-01-01 2010-01-01 false Method of payment; advance deposit. 27.92 Section 27... Micronaire § 27.92 Method of payment; advance deposit. Any payment or advance deposit under this subpart...,” and may not be made in cash except in cases where the total payment or deposit does not exceed...
A Numerical Method for Solving Elasticity Equations with Interfaces
Li, Zhilin; Wang, Liqun; Wang, Wei
2012-01-01
Solving elasticity equations with interfaces is a challenging problem for most existing methods. Nonetheless, it has wide applications in engineering and science. An accurate and efficient method is desired. In this paper, an efficient non-traditional finite element method with non-body-fitting grids is proposed to solve elasticity equations with interfaces. The main idea is to choose the test function basis to be the standard finite element basis independent of the interface and to choose the solution basis to be piecewise linear satisfying the jump conditions across the interface. The resulting linear system of equations is shown to be positive definite under certain assumptions. Numerical experiments show that this method is second order accurate in the L∞ norm for piecewise smooth solutions. More than 1.5th order accuracy is observed for solution with singularity (second derivative blows up) on the sharp-edged interface corner. PMID:22707984
Time-dependent corona models - A numerical method
NASA Astrophysics Data System (ADS)
Korevaar, P.; van Leer, B.
1988-07-01
A time-dependent numerical method for calculating gas flows is described. The method is implicit and especially suitable for finding stationary flow solutions. Although the method is general in its application to ideal compressible fluids, this paper applies it to a stellar atmosphere, heated to coronal temperatures by dissipation of mechanical energy. The integration scheme is based on conservative upwind spatial differencing. The upwind switching is provided by Van Leer's method of differentiable flux-splitting. It is shown that the code can handle large differences in density: up to 14 orders of magnitude. Special attention is paid to the boundary conditions, which are made completely transparent to disturbances. Besides some test-results, converged solutions for various values of the initial mechanical flux are presented which are in good agreement with previous time-independent calculations.
Advanced Electromagnetic Methods for Aerospace Vehicles
NASA Technical Reports Server (NTRS)
Balanis, Constantine A.; Polycarpou, Anastasis; Birtcher, Craig R.; Georgakopoulos, Stavros; Han, Dong-Ho; Ballas, Gerasimos
1999-01-01
The imminent destructive threats of Lightning on helicopters and other airborne systems has always been a topic of great interest to this research grant. Previously, the lightning induced currents on the surface of the fuselage and its interior were predicted using the finite-difference time-domain (FDTD) method as well as the NEC code. The limitations of both methods, as applied to lightning, were identified and extensively discussed in the last meeting. After a thorough investigation of the capabilities of the FDTD, it was decided to incorporate into the numerical method a subcell model to accurately represent current diffusion through conducting materials of high conductivity and finite thickness. Because of the complexity of the model, its validity will be first tested for a one-dimensional FDTD problem. Although results are not available yet, the theory and formulation of the subcell model are presented and discussed here to a certain degree. Besides lightning induced currents in the interior of an aircraft, penetration of electromagnetic fields through apertures (e.g., windows and cracks) could also be devastating for the navigation equipment, electronics, and communications systems in general. The main focus of this study is understanding and quantifying field penetration through apertures. The simulation is done using the FDTD method and the predictions are compared with measurements and moment method solutions obtained from the NASA Langley Research Center. Cavity-backed slot (CBS) antennas or slot antennas in general have many applications in aircraft-satellite type of communications. These can be flushmounted on the surface of the fuselage and, therefore, they retain the aerodynamic shape of the aircraft. In the past, input impedance and radiation patterns of CBS antennas were computed using a hybrid FEM/MoM code. The analysis is now extended to coupling between two identical slot antennas mounted on the same structure. The predictions are performed
Advanced electromagnetic methods for aerospace vehicles
NASA Technical Reports Server (NTRS)
Balanis, Constantine A.; Choi, Jachoon; El-Sharawy, El-Budawy; Hashemi-Yeganeh, Shahrokh; Birtcher, Craig R.
1990-01-01
High- and low-frequency methods to analyze various radiation elements located on aerospace vehicles with combinations of conducting, nonconducting, and energy absorbing surfaces and interfaces. The focus was on developing fundamental concepts, techniques, and algorithms which would remove some of the present limitations in predicting radiation characteristics of antennas on complex aerospace vehicles. In order to accomplish this, the following subjects were examined: (1) the development of techniques for rigorous analysis of surface discontinuities of metallic and nonmetallic surfaces using the equivalent surface impedance concept and Green's function; (2) the effects of anisotropic material on antenna radiation patterns through the use of an equivalent surface impedance concept which is incorporated into the existing numerical electromagnetics computer codes; and (3) the fundamental concepts of precipitation static (P-Static), such as formulations and analytical models. A computer code was used to model the P-Static process on a simple structure. Measurement techniques were also developed to characterized the electrical properties at microwave frequencies. Samples of typical materials used in airframes were tested and the results are included.
Advanced continuous cultivation methods for systems microbiology.
Adamberg, Kaarel; Valgepea, Kaspar; Vilu, Raivo
2015-09-01
Increasing the throughput of systems biology-based experimental characterization of in silico-designed strains has great potential for accelerating the development of cell factories. For this, analysis of metabolism in the steady state is essential as only this enables the unequivocal definition of the physiological state of cells, which is needed for the complete description and in silico reconstruction of their phenotypes. In this review, we show that for a systems microbiology approach, high-resolution characterization of metabolism in the steady state--growth space analysis (GSA)--can be achieved by using advanced continuous cultivation methods termed changestats. In changestats, an environmental parameter is continuously changed at a constant rate within one experiment whilst maintaining cells in the physiological steady state similar to chemostats. This increases the resolution and throughput of GSA compared with chemostats, and, moreover, enables following of the dynamics of metabolism and detection of metabolic switch-points and optimal growth conditions. We also describe the concept, challenge and necessary criteria of the systematic analysis of steady-state metabolism. Finally, we propose that such systematic characterization of the steady-state growth space of cells using changestats has value not only for fundamental studies of metabolism, but also for systems biology-based metabolic engineering of cell factories. PMID:26220303
Advanced electromagnetic methods for aerospace vehicles
NASA Technical Reports Server (NTRS)
Balanis, Constantine A.; Sun, Weimin; El-Sharawy, El-Budawy; Aberle, James T.; Birtcher, Craig R.; Peng, Jian; Tirkas, Panayiotis A.; Andrew, William V.; Kokotoff, David; Zavosh, Frank
1993-01-01
The Advanced Helicopter Electromagnetics (AHE) Industrial Associates Program has fruitfully completed its fourth year. Under the support of the AHE members and the joint effort of the research team, new and significant progress has been achieved in the year. Following the recommendations by the Advisory Task Force, the research effort is placed on more practical helicopter electromagnetic problems, such as HF antennas, composite materials, and antenna efficiencies. In this annual report, the main topics to be addressed include composite materials and antenna technology. The research work on each topic has been driven by the AHE consortium members' interests and needs. The remarkable achievements and progresses in each subject is reported respectively in individual sections of the report. The work in the area of composite materials includes: modeling of low conductivity composite materials by using Green's function approach; guidelines for composite material modeling by using the Green's function approach in the NEC code; development of 3-D volume mesh generator for modeling thick and volumetric dielectrics by using FD-TD method; modeling antenna elements mounted on a composite Comanche tail stabilizer; and antenna pattern control and efficiency estimate for a horn antenna loaded with composite dielectric materials.
Advanced electromagnetic methods for aerospace vehicles
NASA Technical Reports Server (NTRS)
Balanis, Constantine A.; Sun, Weimin; El-Sharawy, El-Budawy; Aberle, James T.; Birtcher, Craig R.; Peng, Jian; Tirkas, Panayiotis A.; Kokotoff, David; Zavosh, Frank
1993-01-01
The Advanced Helicopter Electromagnetics (AHE) Industrial Associates Program has continuously progressed with its research effort focused on subjects identified and recommended by the Advisory Task Force of the program. The research activities in this reporting period have been steered toward practical helicopter electromagnetic problems, such as HF antenna problems and antenna efficiencies, recommended by the AHE members at the annual conference held at Arizona State University on 28-29 Oct. 1992 and the last biannual meeting held at the Boeing Helicopter on 19-20 May 1993. The main topics addressed include the following: Composite Materials and Antenna Technology. The research work on each topic is closely tied with the AHE Consortium members' interests. Significant progress in each subject is reported. Special attention in the area of Composite Materials has been given to the following: modeling of material discontinuity and their effects on towel-bar antenna patterns; guidelines for composite material modeling by using the Green's function approach in the NEC code; measurements of towel-bar antennas grounded with a partially material-coated plate; development of 3-D volume mesh generator for modeling thick and volumetric dielectrics by using FD-TD method; FDTD modeling of horn antennas with composite E-plane walls; and antenna efficiency analysis for a horn antenna loaded with composite dielectric materials.
Advanced electromagnetic methods for aerospace vehicles
NASA Astrophysics Data System (ADS)
Balanis, Constantine A.; Sun, Weimin; El-Sharawy, El-Budawy; Aberle, James T.; Birtcher, Craig R.; Peng, Jian; Tirkas, Panayiotis A.; Kokotoff, David; Zavosh, Frank
1993-06-01
The Advanced Helicopter Electromagnetics (AHE) Industrial Associates Program has continuously progressed with its research effort focused on subjects identified and recommended by the Advisory Task Force of the program. The research activities in this reporting period have been steered toward practical helicopter electromagnetic problems, such as HF antenna problems and antenna efficiencies, recommended by the AHE members at the annual conference held at Arizona State University on 28-29 Oct. 1992 and the last biannual meeting held at the Boeing Helicopter on 19-20 May 1993. The main topics addressed include the following: Composite Materials and Antenna Technology. The research work on each topic is closely tied with the AHE Consortium members' interests. Significant progress in each subject is reported. Special attention in the area of Composite Materials has been given to the following: modeling of material discontinuity and their effects on towel-bar antenna patterns; guidelines for composite material modeling by using the Green's function approach in the NEC code; measurements of towel-bar antennas grounded with a partially material-coated plate; development of 3-D volume mesh generator for modeling thick and volumetric dielectrics by using FD-TD method; FDTD modeling of horn antennas with composite E-plane walls; and antenna efficiency analysis for a horn antenna loaded with composite dielectric materials.
Comparison of four stable numerical methods for Abel's integral equation
NASA Technical Reports Server (NTRS)
Murio, Diego A.; Mejia, Carlos E.
1991-01-01
The 3-D image reconstruction from cone-beam projections in computerized tomography leads naturally, in the case of radial symmetry, to the study of Abel-type integral equations. If the experimental information is obtained from measured data, on a discrete set of points, special methods are needed in order to restore continuity with respect to the data. A new combined Regularized-Adjoint-Conjugate Gradient algorithm, together with two different implementations of the Mollification Method (one based on a data filtering technique and the other on the mollification of the kernal function) and a regularization by truncation method (initially proposed for 2-D ray sample schemes and more recently extended to 3-D cone-beam image reconstruction) are extensively tested and compared for accuracy and numerical stability as functions of the level of noise in the data.
Numerical Analysis of a Finite Element/Volume Penalty Method
NASA Astrophysics Data System (ADS)
Maury, Bertrand
The penalty method makes it possible to incorporate a large class of constraints in general purpose Finite Element solvers like freeFEM++. We present here some contributions to the numerical analysis of this method. We propose an abstract framework for this approach, together with some general error estimates based on the discretization parameter ɛ and the space discretization parameter h. As this work is motivated by the possibility to handle constraints like rigid motion for fluid-particle flows, we shall pay a special attention to a model problem of this kind, where the constraint is prescribed over a subdomain. We show how the abstract estimate can be applied to this situation, in the case where a non-body-fitted mesh is used. In addition, we describe how this method provides an approximation of the Lagrange multiplier associated to the constraint.
Numerical methods for high-dimensional probability density function equations
NASA Astrophysics Data System (ADS)
Cho, H.; Venturi, D.; Karniadakis, G. E.
2016-01-01
In this paper we address the problem of computing the numerical solution to kinetic partial differential equations involving many phase variables. These types of equations arise naturally in many different areas of mathematical physics, e.g., in particle systems (Liouville and Boltzmann equations), stochastic dynamical systems (Fokker-Planck and Dostupov-Pugachev equations), random wave theory (Malakhov-Saichev equations) and coarse-grained stochastic systems (Mori-Zwanzig equations). We propose three different classes of new algorithms addressing high-dimensionality: The first one is based on separated series expansions resulting in a sequence of low-dimensional problems that can be solved recursively and in parallel by using alternating direction methods. The second class of algorithms relies on truncation of interaction in low-orders that resembles the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) framework of kinetic gas theory and it yields a hierarchy of coupled probability density function equations. The third class of algorithms is based on high-dimensional model representations, e.g., the ANOVA method and probabilistic collocation methods. A common feature of all these approaches is that they are reducible to the problem of computing the solution to high-dimensional equations via a sequence of low-dimensional problems. The effectiveness of the new algorithms is demonstrated in numerical examples involving nonlinear stochastic dynamical systems and partial differential equations, with up to 120 variables.
Improved numerical method for subchannel cross-flow calculations
Kaya, S.; Anghaie, S.
1986-01-01
COBRA-OSU is a fast running computer code for coupled kinetic and thermal-hydraulic analysis of nuclear reactor core subchannels, currently under development at Oregon State University. This code is a modified version of COBRA-IV with two major improved features. First, COBRA-OSU uses the Gaussian elimination method instead of Gauss-Seidel iteration for subchannel cross-flow calculation. Second, COBRA-OSU has an additional model for regionwise point reactor kinetics which includes all major feedback reactivity effects on calculation of the axial power profile during the course of a transient. This paper summarizes the improved numerical features of the COBRA-OSU code.
Optimization methods and silicon solar cell numerical models
NASA Technical Reports Server (NTRS)
Girardini, K.
1986-01-01
The goal of this project is the development of an optimization algorithm for use with a solar cell model. It is possible to simultaneously vary design variables such as impurity concentrations, front junction depth, back junctions depth, and cell thickness to maximize the predicted cell efficiency. An optimization algorithm has been developed and interfaced with the Solar Cell Analysis Program in 1 Dimension (SCAPID). SCAPID uses finite difference methods to solve the differential equations which, along with several relations from the physics of semiconductors, describe mathematically the operation of a solar cell. A major obstacle is that the numerical methods used in SCAPID require a significant amount of computer time, and during an optimization the model is called iteratively until the design variables converge to the value associated with the maximum efficiency. This problem has been alleviated by designing an optimization code specifically for use with numerically intensive simulations, to reduce the number of times the efficiency has to be calculated to achieve convergence to the optimal solution. Adapting SCAPID so that it could be called iteratively by the optimization code provided another means of reducing the cpu time required to complete an optimization. Instead of calculating the entire I-V curve, as is usually done in SCAPID, only the efficiency is calculated (maximum power voltage and current) and the solution from previous calculations is used to initiate the next solution.
Teaching Thermal Hydraulics & Numerical Methods: An Introductory Control Volume Primer
D. S. Lucas
2004-10-01
A graduate level course for Thermal Hydraulics (T/H) was taught through Idaho State University in the spring of 2004. A numerical approach was taken for the content of this course since the students were employed at the Idaho National Laboratory and had been users of T/H codes. The majority of the students had expressed an interest in learning about the Courant Limit, mass error, semi-implicit and implicit numerical integration schemes in the context of a computer code. Since no introductory text was found the author developed notes taught from his own research and courses taught for Westinghouse on the subject. The course started with a primer on control volume methods and the construction of a Homogeneous Equilibrium Model (HEM) (T/H) code. The primer was valuable for giving the students the basics behind such codes and their evolution to more complex codes for Thermal Hydraulics and Computational Fluid Dynamics (CFD). The course covered additional material including the Finite Element Method and non-equilibrium (T/H). The control volume primer and the construction of a three-equation (mass, momentum and energy) HEM code are the subject of this paper . The Fortran version of the code covered in this paper is elementary compared to its descendants. The steam tables used are less accurate than the available commercial version written in C Coupled to a Graphical User Interface (GUI). The Fortran version and input files can be downloaded at www.microfusionlab.com.
Validation of a numerical method for unsteady flow calculations
Giles, M.; Haimes, R. . Dept. of Aeronautics and Astronautics)
1993-01-01
This paper describes and validates a numerical method for the calculation of unsteady inviscid and viscous flows. A companion paper compares experimental measurements of unsteady heat transfer on a transonic rotor with the corresponding computational results. The mathematical model is the Reynolds-averaged unsteady Navier-Stokes equations for a compressible ideal gas. Quasi-three-dimensionality is included through the use of a variable streamtube thickness. The numerical algorithm is unusual in two respects: (a) For reasons of efficiency and flexibility, it uses a hybrid Navier-Stokes/Euler method, and (b) to allow for the computation of stator/rotor combinations with arbitrary pitch ratio, a novel space-time coordinate transformation is used. Several test cases are presented to validate the performance of the computer program, UNSFLO. These include: (a) unsteady, inviscid flat plate cascade flows (b) steady and unsteady, viscous flat plate cascade flows, (c) steady turbine heat transfer and loss prediction. In the first two sets of cases comparisons are made with theory, and in the third the comparison is with experimental data.
Asymmetric MRI magnet design using a hybrid numerical method.
Zhao, H; Crozier, S; Doddrell, D M
1999-12-01
This paper describes a hybrid numerical method for the design of asymmetric magnetic resonance imaging magnet systems. The problem is formulated as a field synthesis and the desired current density on the surface of a cylinder is first calculated by solving a Fredholm equation of the first kind. Nonlinear optimization methods are then invoked to fit practical magnet coils to the desired current density. The field calculations are performed using a semi-analytical method. A new type of asymmetric magnet is proposed in this work. The asymmetric MRI magnet allows the diameter spherical imaging volume to be positioned close to one end of the magnet. The main advantages of making the magnet asymmetric include the potential to reduce the perception of claustrophobia for the patient, better access to the patient by attending physicians, and the potential for reduced peripheral nerve stimulation due to the gradient coil configuration. The results highlight that the method can be used to obtain an asymmetric MRI magnet structure and a very homogeneous magnetic field over the central imaging volume in clinical systems of approximately 1.2 m in length. Unshielded designs are the focus of this work. This method is flexible and may be applied to magnets of other geometries. PMID:10579958
A method for improving time-stepping numerics
NASA Astrophysics Data System (ADS)
Williams, P. D.
2012-04-01
In contemporary numerical simulations of the atmosphere, evidence suggests that time-stepping errors may be a significant component of total model error, on both weather and climate time-scales. This presentation will review the available evidence, and will then suggest a simple but effective method for substantially improving the time-stepping numerics at no extra computational expense. The most common time-stepping method is the leapfrog scheme combined with the Robert-Asselin (RA) filter. This method is used in the following atmospheric models (and many more): ECHAM, MAECHAM, MM5, CAM, MESO-NH, HIRLAM, KMCM, LIMA, SPEEDY, IGCM, PUMA, COSMO, FSU-GSM, FSU-NRSM, NCEP-GFS, NCEP-RSM, NSEAM, NOGAPS, RAMS, and CCSR/NIES-AGCM. Although the RA filter controls the time-splitting instability in these models, it also introduces non-physical damping and reduces the accuracy. This presentation proposes a simple modification to the RA filter. The modification has become known as the RAW filter (Williams 2011). When used in conjunction with the leapfrog scheme, the RAW filter eliminates the non-physical damping and increases the amplitude accuracy by two orders, yielding third-order accuracy. (The phase accuracy remains second-order.) The RAW filter can easily be incorporated into existing models, typically via the insertion of just a single line of code. Better simulations are obtained at no extra computational expense. Results will be shown from recent implementations of the RAW filter in various atmospheric models, including SPEEDY and COSMO. For example, in SPEEDY, the skill of weather forecasts is found to be significantly improved. In particular, in tropical surface pressure predictions, five-day forecasts made using the RAW filter have approximately the same skill as four-day forecasts made using the RA filter (Amezcua, Kalnay & Williams 2011). These improvements are encouraging for the use of the RAW filter in other models.
Arkundato, Artoto; Su'ud, Zaki; Abdullah, Mikrajuddin; Widayani,; Celino, Massimo
2012-06-06
In this present work, we report numerical results of iron (cladding) corrosion study in interaction with lead-bismuth eutectic coolant of advanced nuclear reactors. The goal of this work is to study how the oxygen can be used to reduce the corrosion rate of cladding. The molecular dynamics method was applied to simulate corrosion process. By evaluating the diffusion coefficients, RDF functions, MSD curves of the iron and also observed the crystal structure of iron before and after oxygen injection to the coolant then we concluded that a significant and effective reduction can be achieved by issuing about 2% number of oxygen atoms to lead-bismuth eutectic coolant.
Unsaturated Shear Strength and Numerical Analysis Methods for Unsaturated Soils
NASA Astrophysics Data System (ADS)
Kim, D.; Kim, G.; Kim, D.; Baek, H.; Kang, S.
2011-12-01
The angles of shearing resistance(φb) and internal friction(φ') appear to be identical in low suction range, but the angle of shearing resistance shows non-linearity as suction increases. In most numerical analysis however, a fixed value for the angle of shearing resistance is applied even in low suction range for practical reasons, often leading to a false conclusion. In this study, a numerical analysis has been undertaken employing the estimated shear strength curve of unsaturated soils from the residual water content of SWCC proposed by Vanapalli et al.(1996). The result was also compared with that from a fixed value of φb. It is suggested that, in case it is difficult to measure the unsaturated shear strength curve through the triaxial soil tests, the estimated shear strength curve using the residual water content can be a useful alternative. This result was applied for analyzing the slope stablity of unsaturated soils. The effects of a continuous rainfall on slope stability were analyzed using a commercial program "SLOPE/W", with the coupled infiltration analysis program "SEEP/W" from the GEO-SLOPE International Ltd. The results show that, prior to the infiltration by the intensive rainfall, the safety factors using the estimated shear strength curve were substantially higher than that from the fixed value of φb at all time points. After the intensive infiltration, both methods showed a similar behavior.
Libration Orbit Mission Design: Applications of Numerical & Dynamical Methods
NASA Technical Reports Server (NTRS)
Bauer, Frank (Technical Monitor); Folta, David; Beckman, Mark
2002-01-01
Sun-Earth libration point orbits serve as excellent locations for scientific investigations. These orbits are often selected to minimize environmental disturbances and maximize observing efficiency. Trajectory design in support of libration orbits is ever more challenging as more complex missions are envisioned in the next decade. Trajectory design software must be further enabled to incorporate better understanding of the libration orbit solution space and thus improve the efficiency and expand the capabilities of current approaches. The Goddard Space Flight Center (GSFC) is currently supporting multiple libration missions. This end-to-end support consists of mission operations, trajectory design, and control. It also includes algorithm and software development. The recently launched Microwave Anisotropy Probe (MAP) and upcoming James Webb Space Telescope (JWST) and Constellation-X missions are examples of the use of improved numerical methods for attaining constrained orbital parameters and controlling their dynamical evolution at the collinear libration points. This paper presents a history of libration point missions, a brief description of the numerical and dynamical design techniques including software used, and a sample of future GSFC mission designs.
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.
Numerical study of Alfvén eigenmodes in the Experimental Advanced Superconducting Tokamak
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.
Space-time adaptive numerical methods for geophysical applications.
Castro, C E; Käser, M; Toro, E F
2009-11-28
In this paper we present high-order formulations of the finite volume and discontinuous Galerkin finite-element methods for wave propagation problems with a space-time adaptation technique using unstructured meshes in order to reduce computational cost without reducing accuracy. Both methods can be derived in a similar mathematical framework and are identical in their first-order version. In their extension to higher order accuracy in space and time, both methods use spatial polynomials of higher degree inside each element, a high-order solution of the generalized Riemann problem and a high-order time integration method based on the Taylor series expansion. The static adaptation strategy uses locally refined high-resolution meshes in areas with low wave speeds to improve the approximation quality. Furthermore, the time step length is chosen locally adaptive such that the solution is evolved explicitly in time by an optimal time step determined by a local stability criterion. After validating the numerical approach, both schemes are applied to geophysical wave propagation problems such as tsunami waves and seismic waves comparing the new approach with the classical global time-stepping technique. The problem of mesh partitioning for large-scale applications on multi-processor architectures is discussed and a new mesh partition approach is proposed and tested to further reduce computational cost. PMID:19840984
Numerical method of characteristics for one-dimensional blood flow
NASA Astrophysics Data System (ADS)
Acosta, Sebastian; Puelz, Charles; Rivière, Béatrice; Penny, Daniel J.; Rusin, Craig G.
2015-08-01
Mathematical modeling at the level of the full cardiovascular system requires the numerical approximation of solutions to a one-dimensional nonlinear hyperbolic system describing flow in a single vessel. This model is often simulated by computationally intensive methods like finite elements and discontinuous Galerkin, while some recent applications require more efficient approaches (e.g. for real-time clinical decision support, phenomena occurring over multiple cardiac cycles, iterative solutions to optimization/inverse problems, and uncertainty quantification). Further, the high speed of pressure waves in blood vessels greatly restricts the time step needed for stability in explicit schemes. We address both cost and stability by presenting an efficient and unconditionally stable method for approximating solutions to diagonal nonlinear hyperbolic systems. Theoretical analysis of the algorithm is given along with a comparison of our method to a discontinuous Galerkin implementation. Lastly, we demonstrate the utility of the proposed method by implementing it on small and large arterial networks of vessels whose elastic and geometrical parameters are physiologically relevant.
Development of numerical methods to problems of micromechanics
NASA Astrophysics Data System (ADS)
Garcia-Martinez, Jose Ramon
In this dissertation we utilize the finite element method to investigate three micromechanical problems. In Chapter 2, we study the compliance contribution tensor H of multiple branched cracks. The cracks grow from a deltoid pore at their center into a triple crack. For plain strain conditions, two-dimensional models of the branched crack are modeled and solved in ABAQUS. The displacement field over the surface of the branched crack and the deltoid is curve fitted to carry out the integral surface of the compliance contribution tensor H. The predicted values are in good agreement with analytical solution. In Chapter 3 a three-dimensional finite element program using an unaligned mesh with an eight-node isoparametric element is developed to study the compliance contribution tensor H of cavities with superellipsoid shapes. A mesh clustering algorithm to increase the number of elements inside and near the superellipsoid surface to obtain a mesh independent solution is used. The numerical results are compared with the analytical solution of a sphere; the error of the numerical approximation varied from 8 to 11%. It is found that the number of elements inside the superellipsoid are insufficient. An algorithm to mesh independently the volumes inside and outside the cube is proposed to increase the accuracy in the calculation of H. As n1 and n2 increase, the numerical solutions show that, H1111 → 0 and H2211 → 0. Although, for these concave shapes no analytical solution exists a bound of 0 for the terms H 1111 and H2211 is suggested. Finally, in Chapter 4 a numerical verification of the cross-property connection between the effective fluid permeability and the effective electrical conductivity is study. A molecular dynamics algorithm is used to generate a set of different microstructural patterns. The volumetric average over a cubic volume is used to obtain the effective electrical conductivity and the effective fluid permeability. The tortuosity of the porous phase
Why Video? How Technology Advances Method
ERIC Educational Resources Information Center
Downing, Martin J., Jr.
2008-01-01
This paper reports on the use of video to enhance qualitative research. Advances in technology have improved our ability to capture lived experiences through visual means. I reflect on my previous work with individuals living with HIV/AIDS, the results of which are described in another paper, to evaluate the effectiveness of video as a medium that…
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.
a Numerical Method for Stability Analysis of Pinned Flexible Mechanisms
NASA Astrophysics Data System (ADS)
Beale, D. G.; Lee, S. W.
1996-05-01
A technique is presented to investigate the stability of mechanisms with pin-jointed flexible members. The method relies on a special floating frame from which elastic link co-ordinates are defined. Energies are easily developed for use in a Lagrange equation formulation, leading to a set of non-linear and mixed ordinary differential-algebraic equations of motion with constraints. Stability and bifurcation analysis is handled using a numerical procedure (generalized co-ordinate partitioning) that avoids the tedious and difficult task of analytically reducing the system of equations to a number equalling the system degrees of freedom. The proposed method was then applied to (1) a slider-crank mechanism with a flexible connecting rod and crank of constant rotational speed, and (2) a four-bar linkage with a flexible coupler with a constant speed crank. In both cases, a single pinned-pinned beam bending mode is employed to develop resonance curves and stability boundaries in the crank length-crank speed parameter plane. Flip and fold bifurcations are common occurrences in both mechanisms. The accuracy of the proposed method was also verified by comparison with previous experimental results [1].
Introduction to finite-difference methods for numerical fluid dynamics
Scannapieco, E.; Harlow, F.H.
1995-09-01
This work is intended to be a beginner`s exercise book for the study of basic finite-difference techniques in computational fluid dynamics. It is written for a student level ranging from high-school senior to university senior. Equations are derived from basic principles using algebra. Some discussion of partial-differential equations is included, but knowledge of calculus is not essential. The student is expected, however, to have some familiarity with the FORTRAN computer language, as the syntax of the computer codes themselves is not discussed. Topics examined in this work include: one-dimensional heat flow, one-dimensional compressible fluid flow, two-dimensional compressible fluid flow, and two-dimensional incompressible fluid flow with additions of the equations of heat flow and the {Kappa}-{epsilon} model for turbulence transport. Emphasis is placed on numerical instabilities and methods by which they can be avoided, techniques that can be used to evaluate the accuracy of finite-difference approximations, and the writing of the finite-difference codes themselves. Concepts introduced in this work include: flux and conservation, implicit and explicit methods, Lagrangian and Eulerian methods, shocks and rarefactions, donor-cell and cell-centered advective fluxes, compressible and incompressible fluids, the Boussinesq approximation for heat flow, Cartesian tensor notation, the Boussinesq approximation for the Reynolds stress tensor, and the modeling of transport equations. A glossary is provided which defines these and other terms.
Numerical performance of projection methods in finite element consolidation models
NASA Astrophysics Data System (ADS)
Gambolati, Giuseppe; Pini, Giorgio; Ferronato, Massimiliano
2001-12-01
Projection, or conjugate gradient like, methods are becoming increasingly popular for the efficient solution of large sparse sets of unsymmetric indefinite equations arising from the numerical integration of (initial) boundary value problems. One such problem is soil consolidation coupling a flow and a structural model, typically solved by finite elements (FE) in space and a marching scheme in time (e.g. the Crank-Nicolson scheme). The attraction of a projection method stems from a number of factors, including the ease of implementation, the requirement of limited core memory and the low computational cost if a cheap and effective matrix preconditioner is available. In the present paper, biconjugate gradient stabilized (Bi- CGSTAB) is used to solve FE consolidation equations in 2-D and 3-D settings with variable time integration steps. Three different nodal orderings are selected along with the preconditioner ILUT based on incomplete triangular factorization and variable fill-in. The overall cost of the solver is made up of the preconditioning cost plus the cost to converge which is in turn related to the number of iterations and the elementary operations required by each iteration. The results show that nodal ordering affects the perfor mance of Bi-CGSTAB. For normally conditioned consolidation problems Bi-CGSTAB with the best ILUT preconditioner may converge in a number of iterations up to two order of magnitude smaller than the size of the FE model and proves an accurate, cost-effective and robust alternative to direct methods.
Numerical Simulations of Granular Dynamics: Method and Tests
NASA Astrophysics Data System (ADS)
Richardson, Derek C.; Walsh, K. J.; Murdoch, N.; Michel, P.; Schwartz, S. R.
2010-10-01
We present a new particle-based numerical method for the simulation of granular dynamics, with application to motions of particles (regolith) on small solar system bodies and planetary surfaces [1]. The method employs the parallel N-body tree code pkdgrav [2] to search for collisions and compute particle trajectories. Particle confinement is achieved by combining arbitrary combinations of four provided wall primitives, namely infinite plane, finite disk, infinite cylinder, and finite cylinder, and degenerate cases of these. Various wall movements, including translation, oscillation, and rotation, are supported. Several tests of the method are presented, including a model granular "atmosphere” that achieves correct energy equipartition, and a series of tumbler simulations that compare favorably with actual laboratory experiments [3]. DCR and SRS acknowledge NASA Grant No. NNX08AM39G and NSF Grant No. AST0524875; KJW, the Poincaré Fellowship at OCA; NM, Thales Alenia Space and The Open University; and PM and NM, the French Programme National de Planétologie. References: [1] Richardson et al. (2010), Icarus, submitted; [2] Cf. Richardson et al. (2009), P&SS 57, 183 and references therein; [3] Brucks et al. (2007), PRE 75, 032301-1-032301-4.
Controlling template erosion with advanced cleaning methods
NASA Astrophysics Data System (ADS)
Singh, SherJang; Yu, Zhaoning; Wähler, Tobias; Kurataka, Nobuo; Gauzner, Gene; Wang, Hongying; Yang, Henry; Hsu, Yautzong; Lee, Kim; Kuo, David; Dress, Peter
2012-03-01
We studied the erosion and feature stability of fused silica patterns under different template cleaning conditions. The conventional SPM cleaning is compared with an advanced non-acid process. Spectroscopic ellipsometry optical critical dimension (SE-OCD) measurements were used to characterize the changes in pattern profile with good sensitivity. This study confirmed the erosion of the silica patterns in the traditional acid-based SPM cleaning mixture (H2SO4+H2O2) at a rate of ~0.1nm per cleaning cycle. The advanced non-acid clean process however only showed CD shift of ~0.01nm per clean. Contamination removal & pattern integrity of sensitive 20nm features under MegaSonic assisted cleaning is also demonstrated.
Advanced electromagnetic methods for aerospace vehicles
NASA Technical Reports Server (NTRS)
Balanis, Constantine A.; Sun, Weimin; El-Sharawy, El-Budawy; Aberle, James T.; Birtcher, Craig R.; Peng, Jian; Tirkas, Panayiotis A.
1992-01-01
The Advanced Helicopter Electromagnetics (AHE) Industrial Associates Program continues its research on variety of main topics identified and recommended by the Advisory Task Force of the program. The research activities center on issues that advance technology related to helicopter electromagnetics. While most of the topics are a continuation of previous works, special effort has been focused on some of the areas due to recommendations from the last annual conference. The main topics addressed in this report are: composite materials, and antenna technology. The area of composite materials continues getting special attention in this period. The research has focused on: (1) measurements of the electrical properties of low-conductivity materials; (2) modeling of material discontinuity and their effects on the scattering patterns; (3) preliminary analysis on interaction of electromagnetic fields with multi-layered graphite fiberglass plates; and (4) finite difference time domain (FDTD) modeling of fields penetration through composite panels of a helicopter.
Method and apparatus for advancing tethers
Zollinger, W. Thor
1998-01-01
A tether puller for advancing a tether through a channel may include a bellows assembly having a leading end fixedly attached to the tether at a first position and a trailing end fixedly attached to the tether at a second position so that the leading and trailing ends of the bellows assembly are located a substantially fixed distance apart. The bellows assembly includes a plurality of independently inflatable elements each of which may be separately inflated to an extended position and deflated to a retracted position. Each of the independently inflatable elements expands radially and axially upon inflation. An inflation system connected to the independently inflatable elements inflates and deflates selected ones of the independently inflatable elements to cause the bellows assembly to apply a tractive force to the tether and advance it in the channel.
Method and apparatus for advancing tethers
Zollinger, W.T.
1998-06-02
A tether puller for advancing a tether through a channel may include a bellows assembly having a leading end fixedly attached to the tether at a first position and a trailing end fixedly attached to the tether at a second position so that the leading and trailing ends of the bellows assembly are located a substantially fixed distance apart. The bellows assembly includes a plurality of independently inflatable elements each of which may be separately inflated to an extended position and deflated to a retracted position. Each of the independently inflatable elements expands radially and axially upon inflation. An inflation system connected to the independently inflatable elements inflates and deflates selected ones of the independently inflatable elements to cause the bellows assembly to apply a tractive force to the tether and advance it in the channel. 9 figs.
Numerical Comparison of Periodic MoM (Method of Moments) and BMIA (Banded Matrix Iteration Method)
NASA Technical Reports Server (NTRS)
Kim, Y.; Rodriguez, E.; Michel, T.
1995-01-01
The most popular numerical technique in rough surface scattering is the Method of Moments (MoM). Since the scattering patch size is finite, the edge current must be suppressed to obtain accurate scattering cross sections. Two standard ways to minimize the edge current are periodic boundary conditions and incident wave tapering. We compare the accuracy & computational requirements of these methods.
A survey of numerical methods for shock physics applications
Hertel, E.S. Jr.
1997-10-01
Hydrocodes or more accurately, shock physics analysis packages, have been widely used in the US Department of Energy (DOE) laboratories and elsewhere around the world for over 30 years. Initial applications included weapons effects studies where the pressure levels were high enough to disregard the material strength, hence the term hydrocode. Over the last 30 years, Sandia has worked extensively to develop and apply advanced hydrocodes to armor/anti-armor interactions, warhead design, high explosive initiation, and nuclear weapon safety issues. The needs of the DOE have changed over the last 30 years, especially over the last decade. A much stronger emphasis is currently placed on the details of material deformation and high explosive initiation phenomena. The hydrocodes of 30 years ago have now evolved into sophisticated analysis tools that can replace testing in some situations and complement it in all situations. A brief history of the development of hydrocodes in the US will be given. The author also discusses and compares the four principal methods in use today for the solution of the conservation equations of mass, momentum, and energy for shock physics applications. The techniques discussed are the Eulerian methods currently employed by the Sandia multi-dimensional shock physics analysis package known as CTH; the element based Lagrangian method currently used by codes like DYNA; the element free Lagrangian method (also known as smooth particle hydrodynamics) used by codes like the Los Alamos code SPHINX; and the Arbitrary Lagrangian Eulerian methods used by codes like the Lawrence Livermore code CALE or the Sandia code ALEGRA.
A review on recent advances in the numerical simulation for coalbed-methane-recovery process
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.
NATO PILOT STUDY ON ADVANCED CANCER RISK ASSESSMENT METHODS
NCEA scientists are participating in a study of advanced cancer risk assessment methods, conducted under the auspices of NATO's Committee on the Challenges of Modern Society. The product will be a book of case studies that illustrate advanced cancer risk assessment methods, avail...
Advanced methods of structural and trajectory analysis for transport aircraft
NASA Technical Reports Server (NTRS)
Ardema, Mark D.
1995-01-01
This report summarizes the efforts in two areas: (1) development of advanced methods of structural weight estimation, and (2) development of advanced methods of trajectory optimization. The majority of the effort was spent in the structural weight area. A draft of 'Analytical Fuselage and Wing Weight Estimation of Transport Aircraft', resulting from this research, is included as an appendix.
Advances in Analytical and Numerical Dispersion Modeling of Pollutants Releasing from an Area-source
NASA Astrophysics Data System (ADS)
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
Numerical methods for determining interstitial oxygen in silicon
Stevenson, J.O.; Medernach, J.W.
1995-01-01
The interstitial oxygen (O{sub i}) concentration in Czochralski silicon and the subsequent SiO{sub x} precipitation are important parameters for integrated circuit fabrication. Uncontrolled SiO{sub x} precipitation during processing can create detrimental mechanical and electrical effects that contribute to poor performance. An inability to consistently and accurately measure the initial O{sub i} concentration in heavily doped silicon has led to contradictory results regarding the effects of dopant type and concentration on SiO{sub x} precipitation. The authors have developed a software package for reliably determining and comparing O{sub i} in heavily doped silicon. The SiFTIR{copyright} code implements three independent oxygen analysis methods in a single integrated package. Routine oxygen measurements are desirable over a wide range of silicon resistivities, but there has been confusion concerning which of the three numerical methods is most suitable for the low resistivity portion of the continuum. A major strength of the software is an ability to rapidly produce results for all three methods using only a single Fourier Transform Infrared Spectroscopy (FTIR) spectrum as input. This ability to perform three analyses on a single data set allows a detailed comparison of the three methods across the entire range of resistivities in question. Integrated circuit manufacturers could use the enabling technology provided by SiFTIR{copyright} to monitor O{sub i} content. Early detection of O{sub i} using this diagnostic could be beneficial in controlling SiO{sub x} precipitation during integrated circuit processing.
Numerical Methods for Two-Dimensional Stem Cell Tissue Growth.
Ovadia, Jeremy; Nie, Qing
2014-01-01
Growth of developing and regenerative biological tissues of different cell types is usually driven by stem cells and their local environment. Here, we present a computational framework for continuum tissue growth models consisting of stem cells, cell lineages, and diffusive molecules that regulate proliferation and differentiation through feedback. To deal with the moving boundaries of the models in both open geometries and closed geometries (through polar coordinates) in two dimensions, we transform the dynamic domains and governing equations to fixed domains, followed by solving for the transformation functions to track the interface explicitly. Clustering grid points in local regions for better efficiency and accuracy can be achieved by appropriate choices of the transformation. The equations resulting from the incompressibility of the tissue is approximated by high-order finite difference schemes and is solved using the multigrid algorithms. The numerical tests demonstrate an overall spatiotemporal second-order accuracy of the methods and their capability in capturing large deformations of the tissue boundaries. The methods are applied to two biological systems: stratified epithelia for studying the effects of two different types of stem cell niches and the scaling of a morphogen gradient with the size of the Drosophila imaginal wing disc during growth. Direct simulations of both systems suggest that that the computational framework is robust and accurate, and it can incorporate various biological processes critical to stem cell dynamics and tissue growth. PMID:24415847
Numerical Weather Predictions Evaluation Using Spatial Verification Methods
NASA Astrophysics Data System (ADS)
Tegoulias, I.; Pytharoulis, I.; Kotsopoulos, S.; Kartsios, S.; Bampzelis, D.; Karacostas, T.
2014-12-01
During the last years high-resolution numerical weather prediction simulations have been used to examine meteorological events with increased convective activity. Traditional verification methods do not provide the desired level of information to evaluate those high-resolution simulations. To assess those limitations new spatial verification methods have been proposed. In the present study an attempt is made to estimate the ability of the WRF model (WRF -ARW ver3.5.1) to reproduce selected days with high convective activity during the year 2010 using those feature-based verification methods. Three model domains, covering Europe, the Mediterranean Sea and northern Africa (d01), the wider area of Greece (d02) and central Greece - Thessaly region (d03) are used at horizontal grid-spacings of 15km, 5km and 1km respectively. By alternating microphysics (Ferrier, WSM6, Goddard), boundary layer (YSU, MYJ) and cumulus convection (Kain--Fritsch, BMJ) schemes, a set of twelve model setups is obtained. The results of those simulations are evaluated against data obtained using a C-Band (5cm) radar located at the centre of the innermost domain. Spatial characteristics are well captured but with a variable time lag between simulation results and radar data. Acknowledgements: This research is cofinanced by the European Union (European Regional Development Fund) and Greek national funds, through the action "COOPERATION 2011: Partnerships of Production and Research Institutions in Focused Research and Technology Sectors" (contract number 11SYN_8_1088 - DAPHNE) in the framework of the operational programme "Competitiveness and Entrepreneurship" and Regions in Transition (OPC II, NSRF 2007--2013).
Trigonometrically fitted two step hybrid method for the numerical integration of second order IVPs
NASA Astrophysics Data System (ADS)
Monovasilis, Th.; Kalogiratou, Z.; Simos, T. E.
2016-06-01
In this work we consider the numerical integration of second order ODEs where the first derivative is missing. We construct trigonometrically fitted two step hybrid methods. We apply the new methods on the numerical integration of several test problems.
A Numerical Method for Determining Diffusivity from Annealing Experiments
NASA Astrophysics Data System (ADS)
Harris-Kuhlman, K. R.; Kulcinski, G. L.
1998-12-01
Terrestrial analogs of lunar ilmenite (FeTiO3) have been implanted with solar-wind energy 4He at 4 keV and 3He at 3 keV using Plasma Source Ion Implantation (PSII). Isochronal annealing of the samples revealed thermally induced 4He evolution similar to the helium release of the Apollo 11 regoliths reported by Pepin, et. al., [1970]. These annealing experiments are analyzed with a three dimensional numerical method based on Fick's law for diffusion. An iterative method is used to calculate the diffusivity. The code uses an assumed diffusivity to calculate the amount of gas released during a temperature step. The initial depth profile of the implanted species is generated using the TRIM electronic stopping code [Ziegler, 1996]. The calculated value is compared to the measured value and a linear regression is used to calculate a new diffusivity until there is convergence within a specified tolerance level. The diffusivity as a function of temperature is then fitted to an Arrhenius equation. Analysis of results for 4 keV 4He on ilmenite shows two distinct regions of Arrehnius behavior with activation energies of 0.5 +/- 0.1 eV at emperatures below 800 deg C and 1.5 +/- 0.2 eV at temperatures from 800 deg C to 1100 deg C. Pepin, R. O., L. E. Nyquist, D. Phinney, and D. C. Black (1970) "Rare Gases in Apollo 11 Lunar Material," Proceedings of the Apollo 11 Lunar Science Conference, 2, pp. 1435-1454. Ziegler, J. P. (1996) SRIM Instruction Manual: The Stopping and Range of Ions in Matter, (Yorktown, New York: IBM - Research); based on Ziegler, J. P., J. P. Biersack and U. Littmark, The Stopping and Range of Ions in Solids, (New York: Pergamon Press, 1985).
Advanced particulate matter control apparatus and methods
Miller, Stanley J.; Zhuang, Ye; Almlie, Jay C.
2012-01-10
Apparatus and methods for collection and removal of particulate matter, including fine particulate matter, from a gas stream, comprising a unique combination of high collection efficiency and ultralow pressure drop across the filter. The apparatus and method utilize simultaneous electrostatic precipitation and membrane filtration of a particular pore size, wherein electrostatic collection and filtration occur on the same surface.
Asuero, A G; Navas, M J; Jiminez-Trillo, J L
1986-02-01
The spectrophotometric methods applicable to the numerical evaluation of acidity constants of monobasic acids are briefly reviewed. The equations are presented in a form suitable for easy calculation with a programmable pocket calculator. The aim of this paper is to cover a gap in the education analytical literature. PMID:18964064
Active Problem Solving and Applied Research Methods in a Graduate Course on Numerical Methods
ERIC Educational Resources Information Center
Maase, Eric L.; High, Karen A.
2008-01-01
"Chemical Engineering Modeling" is a first-semester graduate course traditionally taught in a lecture format at Oklahoma State University. The course as taught by the author for the past seven years focuses on numerical and mathematical methods as necessary skills for incoming graduate students. Recent changes to the course have included Visual…
Numerical approach for the voloxidation process of an advanced spent fuel conditioning process (ACP)
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)
Advanced spectral methods for climatic time series
Ghil, M.; Allen, M.R.; Dettinger, M.D.; Ide, K.; Kondrashov, D.; Mann, M.E.; Robertson, A.W.; Saunders, A.; Tian, Y.; Varadi, F.; Yiou, P.
2002-01-01
The analysis of univariate or multivariate time series provides crucial information to describe, understand, and predict climatic variability. The discovery and implementation of a number of novel methods for extracting useful information from time series has recently revitalized this classical field of study. Considerable progress has also been made in interpreting the information so obtained in terms of dynamical systems theory. In this review we describe the connections between time series analysis and nonlinear dynamics, discuss signal- to-noise enhancement, and present some of the novel methods for spectral analysis. The various steps, as well as the advantages and disadvantages of these methods, are illustrated by their application to an important climatic time series, the Southern Oscillation Index. This index captures major features of interannual climate variability and is used extensively in its prediction. Regional and global sea surface temperature data sets are used to illustrate multivariate spectral methods. Open questions and further prospects conclude the review.
PRATHAM: Parallel Thermal Hydraulics Simulations using Advanced Mesoscopic Methods
Joshi, Abhijit S; Jain, Prashant K; Mudrich, Jaime A; Popov, Emilian L
2012-01-01
At the Oak Ridge National Laboratory, efforts are under way to develop a 3D, parallel LBM code called PRATHAM (PaRAllel Thermal Hydraulic simulations using Advanced Mesoscopic Methods) to demonstrate the accuracy and scalability of LBM for turbulent flow simulations in nuclear applications. The code has been developed using FORTRAN-90, and parallelized using the message passing interface MPI library. Silo library is used to compact and write the data files, and VisIt visualization software is used to post-process the simulation data in parallel. Both the single relaxation time (SRT) and multi relaxation time (MRT) LBM schemes have been implemented in PRATHAM. To capture turbulence without prohibitively increasing the grid resolution requirements, an LES approach [5] is adopted allowing large scale eddies to be numerically resolved while modeling the smaller (subgrid) eddies. In this work, a Smagorinsky model has been used, which modifies the fluid viscosity by an additional eddy viscosity depending on the magnitude of the rate-of-strain tensor. In LBM, this is achieved by locally varying the relaxation time of the fluid.
Advanced verification methods for OVI security ink
NASA Astrophysics Data System (ADS)
Coombs, Paul G.; McCaffery, Shaun F.; Markantes, Tom
2006-02-01
OVI security ink +, incorporating OVP security pigment* microflakes, enjoys a history of effective document protection. This security feature provides not only first-line recognition by the person on the street, but also facilitates machine-readability. This paper explores the evolution of OVI reader technology from proof-of-concept to miniaturization. Three different instruments have been built to advance the technology of OVI machine verification. A bench-top unit has been constructed which allows users to automatically verify a multitude of different banknotes and OVI images. In addition, high speed modules were fabricated and tested in a state of the art banknote sorting machine. Both units demonstrate the ability of modern optical components to illuminate and collect light reflected from the interference platelets within OVI ink. Electronic hardware and software convert and process the optical information in milliseconds to accurately determine the authenticity of the security feature. Most recently, OVI ink verification hardware has been miniaturized and simplified providing yet another platform for counterfeit protection. These latest devices provide a tool for store clerks and bank tellers to unambiguously determine the validity of banknotes in the time period it takes the cash drawer to be opened.
Indentation Methods in Advanced Materials Research Introduction
Pharr, George Mathews; Cheng, Yang-Tse; Hutchings, Ian; Sakai, Mototsugu; Moody, Neville; Sundararajan, G.; Swain, Michael V.
2009-01-01
Since its commercialization early in the 20th century, indentation testing has played a key role in the development of new materials and understanding their mechanical behavior. Progr3ess in the field has relied on a close marriage between research in the mechanical behavior of materials and contact mechanics. The seminal work of Hertz laid the foundations for bringing these two together, with his contributions still widely utilized today in examining elastic behavior and the physics of fracture. Later, the pioneering work of Tabor, as published in his classic text 'The Hardness of Metals', exapdned this understanding to address the complexities of plasticity. Enormous progress in the field has been achieved in the last decade, made possible both by advances in instrumentation, for example, load and depth-sensing indentation and scanning electron microscopy (SEM) and transmission electron microscopy (TEM) based in situ testing, as well as improved modeling capabilities that use computationally intensive techniques such as finite element analysis and molecular dynamics simulation. The purpose of this special focus issue is to present recent state of the art developments in the field.
Advances in methods for deepwater TLP installations
Wybro, P.G.
1995-10-01
This paper describes a method suitable for installing deepwater TLP structures in water depths beyond 3,000 ft. An overview is presented of previous TLP installation, wherein an evaluation is made of the various methods and their suitability to deepwater applications. A novel method for installation of deepwater TLP`s is described. This method of installation is most suitable for deepwater and/or large TLP structures, but can also be used in moderate water depth as well. The tendon installation method utilizes the so-called Platform Arrestor Concept (PAC), wherein tendon sections are transported by barges to site, and assembled vertically using a dynamically position crane vessel. The tendons are transferred to the platform where they are hung off until there are a full complement of tendons. The hull lock off operation is performed on all tendons simultaneously, avoiding dangerous platform resonant behavior. The installation calls for relatively simple installation equipment, and also enables the use of simple tendon tie-off equipment, such as a single piece nut.
NUMERICAL METHODS FOR THE SIMULATION OF HIGH INTENSITY HADRON SYNCHROTRONS.
LUCCIO, A.; D'IMPERIO, N.; MALITSKY, N.
2005-09-12
Numerical algorithms for PIC simulation of beam dynamics in a high intensity synchrotron on a parallel computer are presented. We introduce numerical solvers of the Laplace-Poisson equation in the presence of walls, and algorithms to compute tunes and twiss functions in the presence of space charge forces. The working code for the simulation here presented is SIMBAD, that can be run as stand alone or as part of the UAL (Unified Accelerator Libraries) package.
Advanced reliability method for fatigue analysis
NASA Technical Reports Server (NTRS)
Wu, Y.-T.; Wirsching, P. H.
1984-01-01
When design factors are considered as random variables and the failure condition cannot be expressed by a closed form algebraic inequality, computations of risk (or probability of failure) may become extremely difficult or very inefficient. This study suggests using a simple and easily constructed second degree polynomial to approximate the complicated limit state in the neighborhood of the design point; a computer analysis relates the design variables at selected points. Then a fast probability integration technique (i.e., the Rackwitz-Fiessler algorithm) can be used to estimate risk. The capability of the proposed method is demonstrated in an example of a low cycle fatigue problem for which a computer analysis is required to perform local strain analysis to relate the design variables. A comparison of the performance of this method is made with a far more costly Monte Carlo solution. Agreement of the proposed method with Monte Carlo is considered to be good.
The role of numerical simulation for the development of an advanced HIFU system
NASA Astrophysics Data System (ADS)
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.
Numerical divergence effects of equivalence theory in the nodal expansion method
Zika, M.R.; Downar, T.J. )
1993-11-01
Accurate solutions of the advanced nodal equations require the use of discontinuity factors (DFs) to account for the homogenization errors that are inherent in all coarse-mesh nodal methods. During the last several years, nodal equivalence theory (NET) has successfully been implemented for the Cartesian geometry and has received widespread acceptance in the light water reactor industry. The extension of NET to other reactor types has had limited success. Recent efforts to implement NET within the framework of the nodal expansion method have successfully been applied to the fast breeder reactor. However, attempts to apply the same methods to thermal reactors such as the Modular High-Temperature Gas Reactor (MHTGR) have led to numerical divergence problems that can be attributed directly to the magnitude of the DFs. In the work performed here, it was found that the numerical problems occur in the inner and upscatter iterations of the solution algorithm. These iterations use a Gauss-Seidel iterative technique that is always convergent for problems with unity DFs. However, for an MHTGR model that requires large DFs, both the inner and upscatter iterations were divergent. Initial investigations into methods for bounding the DFs have proven unsatisfactory as a means of remedying the convergence problems. Although the DFs could be bounded to yield a convergent solution, several cases were encountered where the resulting flux solution was less accurate than the solution without DFs. For the specific case of problems without upscattering, an alternate numerical method for the inner iteration, an LU decomposition, was identified and shown to be feasible.
Advanced method for making vitreous waste forms
Pope, J.M.; Harrison, D.E.
1980-01-01
A process is described for making waste glass that circumvents the problems of dissolving nuclear waste in molten glass at high temperatures. Because the reactive mixing process is independent of the inherent viscosity of the melt, any glass composition can be prepared with equal facility. Separation of the mixing and melting operations permits novel glass fabrication methods to be employed.
A Fast Numerical Method for a Nonlinear Black-Scholes Equation
NASA Astrophysics Data System (ADS)
Koleva, Miglena N.; Vulkov, Lubin G.
2009-11-01
In this paper we will present an effective numerical method for the Black-Scholes equation with transaction costs for the limiting price u(s, t;a). The technique combines the Rothe method with a two-grid (coarse-fine) algorithm for computation of numerical solutions to initial boundary-values problems to this equation. Numerical experiments for comparison the accuracy ant the computational cost of the method with other known numerical schemes are discussed.
Numerical Analysis on the Vortex Pattern and Flux Particle Dispersion in KR Method Using MPS Method
NASA Astrophysics Data System (ADS)
Hirata, N.; Xu, Y.; Anzai, K.
2015-06-01
The mechanically-stirring vessel is widely used in many fields, such as chemical reactor, bioreactor, and metallurgy, etc. The type of vortex mode that formed during impeller stirring has great effect on stirring efficiency, chemical reacting rate and air entrapment. Many efforts have been made to numerically simulate the fluid flow in the stirring vessel with classical Eulerian method. However, it is difficult to directly investigate the vortex mode and flux particle dispersion. Therefore, moving particle semi-implicit (MPS) method, which is based on Lagrangian method, is applied to simulate the fluid flow in a KR method in this practice. Top height and bottom heights of vortex surface in a steady state under several rotation speed was taken as key parameters to compare the results of numerical and published results. Flux particle dispersion behaviour under a rotation speed range from 80 to 480 rpm was also compared with the past study. The result shows that the numerical calculation has high consistency with experimental results. It is confirmed that the calculation using MPS method well reflected the vortex mode and flux particle dispersion in a mechanically-stirring vessel.
Advancing-layers method for generation of unstructured viscous grids
NASA Technical Reports Server (NTRS)
Pirzadeh, Shahyar
1993-01-01
A novel approach for generating highly stretched grids which is based on a modified advancing-front technique and benefits from the generality, flexibility, and grid quality of the conventional advancing-front-based Euler grid generators is presented. The method is self-sufficient for the insertion of grid points in the boundary layer and beyond. Since it is based on a totally unstructured grid strategy, the method alleviates the difficulties stemming from the structural limitations of the prismatic techniques.
SAMSAN- MODERN NUMERICAL METHODS FOR CLASSICAL SAMPLED SYSTEM ANALYSIS
NASA Technical Reports Server (NTRS)
Frisch, H. P.
1994-01-01
SAMSAN was developed to aid the control system analyst by providing a self consistent set of computer algorithms that support large order control system design and evaluation studies, with an emphasis placed on sampled system analysis. Control system analysts have access to a vast array of published algorithms to solve an equally large spectrum of controls related computational problems. The analyst usually spends considerable time and effort bringing these published algorithms to an integrated operational status and often finds them less general than desired. SAMSAN reduces the burden on the analyst by providing a set of algorithms that have been well tested and documented, and that can be readily integrated for solving control system problems. Algorithm selection for SAMSAN has been biased toward numerical accuracy for large order systems with computational speed and portability being considered important but not paramount. In addition to containing relevant subroutines from EISPAK for eigen-analysis and from LINPAK for the solution of linear systems and related problems, SAMSAN contains the following not so generally available capabilities: 1) Reduction of a real non-symmetric matrix to block diagonal form via a real similarity transformation matrix which is well conditioned with respect to inversion, 2) Solution of the generalized eigenvalue problem with balancing and grading, 3) Computation of all zeros of the determinant of a matrix of polynomials, 4) Matrix exponentiation and the evaluation of integrals involving the matrix exponential, with option to first block diagonalize, 5) Root locus and frequency response for single variable transfer functions in the S, Z, and W domains, 6) Several methods of computing zeros for linear systems, and 7) The ability to generate documentation "on demand". All matrix operations in the SAMSAN algorithms assume non-symmetric matrices with real double precision elements. There is no fixed size limit on any matrix in any
NASA Technical Reports Server (NTRS)
Przekwas, A. J.; Yang, H. Q.
1989-01-01
The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.
Cost estimating methods for advanced space systems
NASA Technical Reports Server (NTRS)
Cyr, Kelley
1988-01-01
Parametric cost estimating methods for space systems in the conceptual design phase are developed. The approach is to identify variables that drive cost such as weight, quantity, development culture, design inheritance, and time. The relationship between weight and cost is examined in detail. A theoretical model of cost is developed and tested statistically against a historical data base of major research and development programs. It is concluded that the technique presented is sound, but that it must be refined in order to produce acceptable cost estimates.
Advances in organometallic synthesis with mechanochemical methods.
Rightmire, Nicholas R; Hanusa, Timothy P
2016-02-14
Solvent-based syntheses have long been normative in all areas of chemistry, although mechanochemical methods (specifically grinding and milling) have been used to good effect for decades in organic, and to a lesser but growing extent, inorganic coordination chemistry. Organometallic synthesis, in contrast, represents a relatively underdeveloped area for mechanochemical research, and the potential benefits are considerable. From access to new classes of unsolvated complexes, to control over stoichiometries that have not been observed in solution routes, mechanochemical (or 'M-chem') approaches have much to offer the synthetic chemist. It has already become clear that removing the solvent from an organometallic reaction can change reaction pathways considerably, so that prediction of the outcome is not always straightforward. This Perspective reviews recent developments in the field, and describes equipment that can be used in organometallic synthesis. Synthetic chemists are encouraged to add mechanochemical methods to their repertoire in the search for new and highly reactive metal complexes and novel types of organometallic transformations. PMID:26763151
Advancements in Research Synthesis Methods: From a Methodologically Inclusive Perspective
ERIC Educational Resources Information Center
Suri, Harsh; Clarke, David
2009-01-01
The dominant literature on research synthesis methods has positivist and neo-positivist origins. In recent years, the landscape of research synthesis methods has changed rapidly to become inclusive. This article highlights methodologically inclusive advancements in research synthesis methods. Attention is drawn to insights from interpretive,…
Current methods and advances in bone densitometry
NASA Technical Reports Server (NTRS)
Guglielmi, G.; Gluer, C. C.; Majumdar, S.; Blunt, B. A.; Genant, H. K.
1995-01-01
Bone mass is the primary, although not the only, determinant of fracture. Over the past few years a number of noninvasive techniques have been developed to more sensitively quantitate bone mass. These include single and dual photon absorptiometry (SPA and DPA), single and dual X-ray absorptiometry (SXA and DXA) and quantitative computed tomography (QCT). While differing in anatomic sites measured and in their estimates of precision, accuracy, and fracture discrimination, all of these methods provide clinically useful measurements of skeletal status. It is the intent of this review to discuss the pros and cons of these techniques and to present the new applications of ultrasound (US) and magnetic resonance (MRI) in the detection and management of osteoporosis.
Advances in LC: bioanalytical method transfer.
Wright, Patricia; Wright, Adrian
2016-09-01
There are three main reasons for transferring from an existing bioanalytical assay to an alternative chromatographic method: speed, cost and sensitivity. These represent a challenge to the analyst in that there is an interplay between these three considerations and one factor is often improved at the expense of another. These three factors act as drivers to encourage technology development and support its uptake. The more recently introduced chromatographic technologies may show significant improvements against one of more of these factors relative to conventional 4.6-mm id reversed-phase HPLC. In this article, some of these new chromatographic approaches will be considered in terms of what they can offer the bioanalysts. PMID:27491842
Current methods and advances in bone densitometry.
Guglielmi, G; Gluer, C C; Majumdar, S; Blunt, B A; Genant, H K
1995-01-01
Bone mass is the primary, although not the only, determinant of fracture. Over the past few years a number of noninvasive techniques have been developed to more sensitively quantitate bone mass. These include single and dual photon absorptiometry (SPA and DPA), single and dual X-ray absorptiometry (SXA and DXA) and quantitative computed tomography (QCT). While differing in anatomic sites measured and in their estimates of precision, accuracy, and fracture discrimination, all of these methods provide clinically useful measurements of skeletal status. It is the intent of this review to discuss the pros and cons of these techniques and to present the new applications of ultrasound (US) and magnetic resonance (MRI) in the detection and management of osteoporosis. PMID:11539928
Advanced Fuzzy Potential Field Method for Mobile Robot Obstacle Avoidance.
Park, Jong-Wook; Kwak, Hwan-Joo; Kang, Young-Chang; Kim, Dong W
2016-01-01
An advanced fuzzy potential field method for mobile robot obstacle avoidance is proposed. The potential field method primarily deals with the repulsive forces surrounding obstacles, while fuzzy control logic focuses on fuzzy rules that handle linguistic variables and describe the knowledge of experts. The design of a fuzzy controller--advanced fuzzy potential field method (AFPFM)--that models and enhances the conventional potential field method is proposed and discussed. This study also examines the rule-explosion problem of conventional fuzzy logic and assesses the performance of our proposed AFPFM through simulations carried out using a mobile robot. PMID:27123001
Advanced Fuzzy Potential Field Method for Mobile Robot Obstacle Avoidance
Park, Jong-Wook; Kwak, Hwan-Joo; Kang, Young-Chang; Kim, Dong W.
2016-01-01
An advanced fuzzy potential field method for mobile robot obstacle avoidance is proposed. The potential field method primarily deals with the repulsive forces surrounding obstacles, while fuzzy control logic focuses on fuzzy rules that handle linguistic variables and describe the knowledge of experts. The design of a fuzzy controller—advanced fuzzy potential field method (AFPFM)—that models and enhances the conventional potential field method is proposed and discussed. This study also examines the rule-explosion problem of conventional fuzzy logic and assesses the performance of our proposed AFPFM through simulations carried out using a mobile robot. PMID:27123001
Cost estimating methods for advanced space systems
NASA Technical Reports Server (NTRS)
Cyr, Kelley
1994-01-01
NASA is responsible for developing much of the nation's future space technology. Cost estimates for new programs are required early in the planning process so that decisions can be made accurately. Because of the long lead times required to develop space hardware, the cost estimates are frequently required 10 to 15 years before the program delivers hardware. The system design in conceptual phases of a program is usually only vaguely defined and the technology used is so often state-of-the-art or beyond. These factors combine to make cost estimating for conceptual programs very challenging. This paper describes an effort to develop parametric cost estimating methods for space systems in the conceptual design phase. The approach is to identify variables that drive cost such as weight, quantity, development culture, design inheritance and time. The nature of the relationships between the driver variables and cost will be discussed. In particular, the relationship between weight and cost will be examined in detail. A theoretical model of cost will be developed and tested statistically against a historical database of major research and development projects.
An Improved Numerical Integration Method for Springback Predictions
NASA Astrophysics Data System (ADS)
Ibrahim, R.; Smith, L. M.; Golovashchenko, Sergey F.
2011-08-01
In this investigation, the focus is on the springback of steel sheets in V-die air bending. A full replication to a numerical integration algorithm presented rigorously in [1] to predict the springback in air bending was performed and confirmed successfully. Algorithm alteration and extensions were proposed here. The altered approach used in solving the moment equation numerically resulted in springback values much closer to the trend presented by the experimental data, Although investigation here extended to use a more realistic work-hardening model, the differences in the springback values obtained by both hardening models were almost negligible. The algorithm was extended to be applied on thin sheets down to 0.8 mm. Results show that this extension is possible as verified by FEA and other published experiments on TRIP steel sheets.
Some variance reduction methods for numerical stochastic homogenization.
Blanc, X; Le Bris, C; Legoll, F
2016-04-28
We give an overview of a series of recent studies devoted to variance reduction techniques for numerical stochastic homogenization. Numerical homogenization requires that a set of problems is solved at the microscale, the so-called corrector problems. In a random environment, these problems are stochastic and therefore need to be repeatedly solved, for several configurations of the medium considered. An empirical average over all configurations is then performed using the Monte Carlo approach, so as to approximate the effective coefficients necessary to determine the macroscopic behaviour. Variance severely affects the accuracy and the cost of such computations. Variance reduction approaches, borrowed from other contexts in the engineering sciences, can be useful. Some of these variance reduction techniques are presented, studied and tested here. PMID:27002065
Algorithms for the Fractional Calculus: A Selection of Numerical Methods
NASA Technical Reports Server (NTRS)
Diethelm, K.; Ford, N. J.; Freed, A. D.; Luchko, Yu.
2003-01-01
Many recently developed models in areas like viscoelasticity, electrochemistry, diffusion processes, etc. are formulated in terms of derivatives (and integrals) of fractional (non-integer) order. In this paper we present a collection of numerical algorithms for the solution of the various problems arising in this context. We believe that this will give the engineer the necessary tools required to work with fractional models in an efficient way.
On numerical methods in non-Newtonian flows
NASA Astrophysics Data System (ADS)
Fileas, G.
1982-12-01
The constitutive equations for non-Newtonian flows are presented and the various flow models derived from continuum mechanics and molecular theories are considered and evaluated. Detailed account is given of numerical simulation employing differential and integral models of different kinds of non-Newtonian flows using finite difference and finite element techniques. Procedures for computer set ups are described and references are given for finite difference, finite element and molecular theory based programs for several kinds of flow. Achievements and unreached goals in the field of numerical simulation of non-Newtonian flows are discussed and the lack of numerical work in the fields of suspension flows and heat transfer is pointed out. Finally, FFOCUS is presented as a newly built computer program which can simulate freezing flows of Newtonian fluids through various geometries and is aimed to be further developed to handle non-Newtonian freezing flows and certain types of suspension phenomena involved in corium flow after a hypothetical core melt down accident in a pressurized water reactor.
Unstructured viscous grid generation by advancing-front method
NASA Technical Reports Server (NTRS)
Pirzadeh, Shahyar
1993-01-01
A new method of generating unstructured triangular/tetrahedral grids with high-aspect-ratio cells is proposed. The method is based on new grid-marching strategy referred to as 'advancing-layers' for construction of highly stretched cells in the boundary layer and the conventional advancing-front technique for generation of regular, equilateral cells in the inviscid-flow region. Unlike the existing semi-structured viscous grid generation techniques, the new procedure relies on a totally unstructured advancing-front grid strategy resulting in a substantially enhanced grid flexibility and efficiency. The method is conceptually simple but powerful, capable of producing high quality viscous grids for complex configurations with ease. A number of two-dimensional, triangular grids are presented to demonstrate the methodology. The basic elements of the method, however, have been primarily designed with three-dimensional problems in mind, making it extendible for tetrahedral, viscous grid generation.
Advanced Ablative Insulators and Methods of Making Them
NASA Technical Reports Server (NTRS)
Congdon, William M.
2005-01-01
Advanced ablative (more specifically, charring) materials that provide temporary protection against high temperatures, and advanced methods of designing and manufacturing insulators based on these materials, are undergoing development. These materials and methods were conceived in an effort to replace the traditional thermal-protection systems (TPSs) of re-entry spacecraft with robust, lightweight, better-performing TPSs that can be designed and manufactured more rapidly and at lower cost. These materials and methods could also be used to make improved TPSs for general aerospace, military, and industrial applications.
Strategy to Promote Active Learning of an Advanced Research Method
ERIC Educational Resources Information Center
McDermott, Hilary J.; Dovey, Terence M.
2013-01-01
Research methods courses aim to equip students with the knowledge and skills required for research yet seldom include practical aspects of assessment. This reflective practitioner report describes and evaluates an innovative approach to teaching and assessing advanced qualitative research methods to final-year psychology undergraduate students. An…
Magnetohydrodynamic (MHD) modelling of solar active phenomena via numerical methods
NASA Technical Reports Server (NTRS)
Wu, S. T.
1988-01-01
Numerical ideal MHD models for the study of solar active phenomena are summarized. Particular attention is given to the following physical phenomena: (1) local heating of a coronal loop in an isothermal and stratified atmosphere, and (2) the coronal dynamic responses due to magnetic field movement. The results suggest that local heating of a magnetic loop will lead to the enhancement of the density of the neighboring loops through MHD wave compression. It is noted that field lines can be pinched off and may form a self-contained magnetized plasma blob that may move outward into interplanetary space.
Achieving better cooling of turbine blades using numerical simulation methods
NASA Astrophysics Data System (ADS)
Inozemtsev, A. A.; Tikhonov, A. S.; Sendyurev, C. I.; Samokhvalov, N. Yu.
2013-02-01
A new design of the first-stage nozzle vane for the turbine of a prospective gas-turbine engine is considered. The blade's thermal state is numerically simulated in conjugate statement using the ANSYS CFX 13.0 software package. Critical locations in the blade design are determined from the distribution of heat fluxes, and measures aimed at achieving more efficient cooling are analyzed. Essentially lower (by 50-100°C) maximal temperature of metal has been achieved owing to the results of the performed work.
Path Integrals and Exotic Options:. Methods and Numerical Results
NASA Astrophysics Data System (ADS)
Bormetti, G.; Montagna, G.; Moreni, N.; Nicrosini, O.
2005-09-01
In the framework of Black-Scholes-Merton model of financial derivatives, a path integral approach to option pricing is presented. A general formula to price path dependent options on multidimensional and correlated underlying assets is obtained and implemented by means of various flexible and efficient algorithms. As an example, we detail the case of Asian call options. The numerical results are compared with those obtained with other procedures used in quantitative finance and found to be in good agreement. In particular, when pricing at the money (ATM) and out of the money (OTM) options, path integral exhibits competitive performances.
Numerical methods for a general class of porous medium equations
Rose, M. E.
1980-03-01
The partial differential equation par. deltau/par. deltat + par. delta(f(u))/par. deltax = par. delta(g(u)par. deltau/par. deltax)/par. deltax, where g(u) is a non-negative diffusion coefficient that may vanish for one or more values of u, was used to model fluid flow through a porous medium. Error estimates for a numerical procedure to approximate the solution are derived. A revised version of this report will appear in Computers and Mathematics with Applications.
A Primer In Advanced Fatigue Life Prediction Methods
NASA Technical Reports Server (NTRS)
Halford, Gary R.
2000-01-01
Metal fatigue has plagued structural components for centuries, and it remains a critical durability issue in today's aerospace hardware. This is true despite vastly improved and advanced materials, increased mechanistic understanding, and development of accurate structural analysis and advanced fatigue life prediction tools. Each advance is quickly taken advantage of to produce safer, more reliable more cost effective, and better performing products. In other words, as the envelop is expanded, components are then designed to operate just as close to the newly expanded envelop as they were to the initial one. The problem is perennial. The economic importance of addressing structural durability issues early in the design process is emphasized. Tradeoffs with performance, cost, and legislated restrictions are pointed out. Several aspects of structural durability of advanced systems, advanced materials and advanced fatigue life prediction methods are presented. Specific items include the basic elements of durability analysis, conventional designs, barriers to be overcome for advanced systems, high-temperature life prediction for both creep-fatigue and thermomechanical fatigue, mean stress effects, multiaxial stress-strain states, and cumulative fatigue damage accumulation assessment.
Projection methods for the numerical solution of Markov chain models
NASA Technical Reports Server (NTRS)
Saad, Youcef
1989-01-01
Projection methods for computing stationary probability distributions for Markov chain models are presented. A general projection method is a method which seeks an approximation from a subspace of small dimension to the original problem. Thus, the original matrix problem of size N is approximated by one of dimension m, typically much smaller than N. A particularly successful class of methods based on this principle is that of Krylov subspace methods which utilize subspaces of the form span(v,av,...,A(exp m-1)v). These methods are effective in solving linear systems and eigenvalue problems (Lanczos, Arnoldi,...) as well as nonlinear equations. They can be combined with more traditional iterative methods such as successive overrelaxation, symmetric successive overrelaxation, or with incomplete factorization methods to enhance convergence.
Advanced Numerical Imaging Procedure Accounting for Non-Ideal Effects in GPR Scenarios
NASA Astrophysics Data System (ADS)
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.
Method for numerical simulation of two-term exponentially correlated colored noise
Yilmaz, B.; Ayik, S.; Abe, Y.; Gokalp, A.; Yilmaz, O.
2006-04-15
A method for numerical simulation of two-term exponentially correlated colored noise is proposed. The method is an extension of traditional method for one-term exponentially correlated colored noise. The validity of the algorithm is tested by comparing numerical simulations with analytical results in two physical applications.
Numerical conformal mapping methods for exterior and doubly connected regions
DeLillo, T.K.; Pfaltzgraff, J.A.
1996-12-31
Methods are presented and analyzed for approximating the conformal map from the exterior of the disk to the exterior a smooth, simple closed curve and from an annulus to a bounded, doubly connected region with smooth boundaries. The methods are Newton-like methods for computing the boundary correspondences and conformal moduli similar to Fornberg`s method for the interior of the disk. We show that the linear systems are discretizations of the identity plus a compact operator and, hence, that the conjugate gradient method converges superlinearly.
Advanced stress analysis methods applicable to turbine engine structures
NASA Technical Reports Server (NTRS)
Pian, T. H. H.
1985-01-01
Advanced stress analysis methods applicable to turbine engine structures are investigated. Constructions of special elements which containing traction-free circular boundaries are investigated. New versions of mixed variational principle and version of hybrid stress elements are formulated. A method is established for suppression of kinematic deformation modes. semiLoof plate and shell elements are constructed by assumed stress hybrid method. An elastic-plastic analysis is conducted by viscoplasticity theory using the mechanical subelement model.
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-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
Recent advances in theoretical and numerical studies of wire array Z-pinch in the IAPCM
NASA Astrophysics Data System (ADS)
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
Analysis of free turbulent shear flows by numerical methods
NASA Technical Reports Server (NTRS)
Korst, H. H.; Chow, W. L.; Hurt, R. F.; White, R. A.; Addy, A. L.
1973-01-01
Studies are described in which the effort was essentially directed to classes of problems where the phenomenologically interpreted effective transport coefficients could be absorbed by, and subsequently extracted from (by comparison with experimental data), appropriate coordinate transformations. The transformed system of differential equations could then be solved without further specifications or assumptions by numerical integration procedures. An attempt was made to delineate different regimes for which specific eddy viscosity models could be formulated. In particular, this would account for the carryover of turbulence from attached boundary layers, the transitory adjustment, and the asymptotic behavior of initially disturbed mixing regions. Such models were subsequently used in seeking solutions for the prescribed two-dimensional test cases, yielding a better insight into overall aspects of the exchange mechanisms.
Numerical method to digital photoelasticity using plane polariscope.
Júnior, P A A Magalhães; Vieira, F G; Magalhães, C A; Ribeiro, J S; Rios, I G
2016-06-13
This research aims to find a new way to get the intensity equations for the phase-shifting model in digital photoelasticity. The procedure is based on the rotation of the analyzer itself. From the intensity equations, the isoclinic and isochromatic equations parameters are deduced by applying a new numerical technique. This approach can be done to calculate how many images allow the resolution of the polariscope. Each image indicates the stress forces in the object. In this study the plane polariscope was used. The amount of images will determinate the number of errors and uncertainties of the study, due to the observation that the veracity of the equations increases considerably with a large amout of images. Several analyses are performed with different amounts of photographic images. The results showed the possibility to measure stress forces with high precision using plane polariscopes. PMID:27410283
MODELING COLLISIONAL CASCADES IN DEBRIS DISKS: THE NUMERICAL METHOD
Gaspar, Andras; Psaltis, Dimitrios; Oezel, Feryal; Rieke, George H.; Cooney, Alan E-mail: dpsaltis@as.arizona.edu E-mail: grieke@as.arizona.edu
2012-04-10
We develop a new numerical algorithm to model collisional cascades in debris disks. Because of the large dynamical range in particle masses, we solve the integro-differential equations describing erosive and catastrophic collisions in a particle-in-a-box approach, while treating the orbital dynamics of the particles in an approximate fashion. We employ a new scheme for describing erosive (cratering) collisions that yields a continuous set of outcomes as a function of colliding masses. We demonstrate the stability and convergence characteristics of our algorithm and compare it with other treatments. We show that incorporating the effects of erosive collisions results in a decay of the particle distribution that is significantly faster than with purely catastrophic collisions.
On a New Numerical Method for Solving General Variational Inequalities
NASA Astrophysics Data System (ADS)
Bnouhachem, Abdellah; Noor, Muhammad Aslam; Khalfaoui, Mohamed; Sheng, Zhaohan
In this paper, we suggest and analyze a new extragradient method for solving the general variational inequalities involving two operators. We also prove the global convergence of the proposed modified method under certain mild conditions. We used a self-adaptive technique to adjust parameter ρ at each iteration. It is proved theoretically that the lower-bound of the progress obtained by the proposed method is greater than that by the extragradient method. An example is given to illustrate the efficiency and its comparison with the extragradient method. Since the general variational inequalities include the classical variational inequalities and complementarity problems as special cases, our results obtained in this paper continue to hold for these problems. Results obtained in this paper may be viewed as an improvement and refinement of the previously known results in this field.
Exploring the Use of Discontinuous Galerkin Methods for Numerical Relativity
NASA Astrophysics Data System (ADS)
Hebert, Francois; Kidder, Lawrence; Teukolsky, Saul; SXS Collaboration
2015-04-01
The limited accuracy of relativistic hydrodynamic simulations constrains our insight into several important research problems, including among others our ability to generate accurate template waveforms for black hole-neutron star mergers, or our understanding of supernova explosion mechanisms. In many codes the algorithms used to evolve the matter, based on the finite volume method, struggle to reach the desired accuracy. We aim to show improved accuracy by using a discontinuous Galerkin method. This method's attractiveness comes from its combination of spectral convergence properties for smooth solutions and robust stability properties for shocks. We present the status of our work implementing a testbed GR-hydro code using discontinuous Galerkin.
A numerical method for eigenvalue problems in modeling liquid crystals
Baglama, J.; Farrell, P.A.; Reichel, L.; Ruttan, A.; Calvetti, D.
1996-12-31
Equilibrium configurations of liquid crystals in finite containments are minimizers of the thermodynamic free energy of the system. It is important to be able to track the equilibrium configurations as the temperature of the liquid crystals decreases. The path of the minimal energy configuration at bifurcation points can be computed from the null space of a large sparse symmetric matrix. We describe a new variant of the implicitly restarted Lanczos method that is well suited for the computation of extreme eigenvalues of a large sparse symmetric matrix, and we use this method to determine the desired null space. Our implicitly restarted Lanczos method determines adoptively a polynomial filter by using Leja shifts, and does not require factorization of the matrix. The storage requirement of the method is small, and this makes it attractive to use for the present application.
Numerical Stability and Convergence of Approximate Methods for Conservation Laws
NASA Astrophysics Data System (ADS)
Galkin, V. A.
We present the new approach to background of approximate methods convergence based on functional solutions theory for conservation laws. The applications to physical kinetics, gas and fluid dynamics are considered.
NASA Astrophysics Data System (ADS)
Min, Xiaoyi
This thesis first presents the study of the interaction of electromagnetic waves with three-dimensional heterogeneous, dielectric, magnetic, and lossy bodies by surface integral equation modeling. Based on the equivalence principle, a set of coupled surface integral equations is formulated and then solved numerically by the method of moments. Triangular elements are used to model the interfaces of the heterogeneous body, and vector basis functions are defined to expand the unknown current in the formulation. The validity of this formulation is verified by applying it to concentric spheres for which an exact solution exists. The potential applications of this formulation to a partially coated sphere and a homogeneous human body are discussed. Next, this thesis also introduces an efficient new set of integral equations for treating the scattering problem of a perfectly conducting body coated with a thin magnetically lossy layer. These electric field integral equations and magnetic field integral equations are numerically solved by the method of moments (MoM). To validate the derived integral equations, an alternative method to solve the scattering problem of an infinite circular cylinder coated with a thin magnetic lossy layer has also been developed, based on the eigenmode expansion. Results for the radar cross section and current densities via the MoM and the eigenmode expansion method are compared. The agreement is excellent. The finite difference time domain method is subsequently implemented to solve a metallic object coated with a magnetic thin layer and numerical results are compared with that by the MoM. Finally, this thesis presents an application of the finite-difference time-domain approach to the problem of electromagnetic receiving and scattering by a cavity -backed antenna situated on an infinite conducting plane. This application involves modifications of Yee's model, which applies the difference approximations of field derivatives to differential
METHODS ADVANCEMENT FOR MILK ANALYSIS: THE MAMA STUDY
The Methods Advancement for Milk Analysis (MAMA) study was designed by US EPA and CDC investigators to provide data to support the technological and study design needs of the proposed National Children=s Study (NCS). The NCS is a multi-Agency-sponsored study, authorized under the...
Advanced boundary layer transition measurement methods for flight applications
NASA Technical Reports Server (NTRS)
Holmes, B. J.; Croom, C. C.; Gail, P. D.; Manuel, G. S.; Carraway, D. L.
1986-01-01
In modern laminar flow flight research, it is important to understand the specific cause(s) of laminar to turbulent boundary-layer transition. Such information is crucial to the exploration of the limits of practical application of laminar flow for drag reduction on aircraft. The transition modes of interest in current flight investigations include the viscous Tollmien-Schlichting instability, the inflectional instability at laminar separation, and the crossflow inflectional instability, as well as others. This paper presents the results to date of research on advanced devices and methods used for the study of laminar boundary-layer transition phenomena in the flight environment. Recent advancements in the development of arrayed hot-film devices and of a new flow visualization method are discussed. Arrayed hot-film devices have been designed to detect the presence of laminar separation, and of crossflow vorticity. The advanced flow visualization method utilizes color changes in liquid-crystal coatings to detect boundary-layer transition at high altitude flight conditions. Flight and wind tunnel data are presented to illustrate the design and operation of these advanced methods. These new research tools provide information on disturbance growth and transition mode which is essential to furthering our understanding of practical design limits for applications of laminar flow technology.
Domain Decomposition By the Advancing-Partition Method
NASA Technical Reports Server (NTRS)
Pirzadeh, Shahyar Z.
2008-01-01
A new method of domain decomposition has been developed for generating unstructured grids in subdomains either sequentially or using multiple computers in parallel. Domain decomposition is a crucial and challenging step for parallel grid generation. Prior methods are generally based on auxiliary, complex, and computationally intensive operations for defining partition interfaces and usually produce grids of lower quality than those generated in single domains. The new technique, referred to as "Advancing Partition," is based on the Advancing-Front method, which partitions a domain as part of the volume mesh generation in a consistent and "natural" way. The benefits of this approach are: 1) the process of domain decomposition is highly automated, 2) partitioning of domain does not compromise the quality of the generated grids, and 3) the computational overhead for domain decomposition is minimal. The new method has been implemented in NASA's unstructured grid generation code VGRID.
NASA Technical Reports Server (NTRS)
Stern, Martin O.
1992-01-01
This report describes a study to evaluate the benefits of advanced propulsion technologies for transporting materials between low Earth orbit and the Moon. A relatively conventional reference transportation system, and several other systems, each of which includes one advanced technology component, are compared in terms of how well they perform a chosen mission objective. The evaluation method is based on a pairwise life-cycle cost comparison of each of the advanced systems with the reference system. Somewhat novel and economically important features of the procedure are the inclusion not only of mass payback ratios based on Earth launch costs, but also of repair and capital acquisition costs, and of adjustments in the latter to reflect the technological maturity of the advanced technologies. The required input information is developed by panels of experts. The overall scope and approach of the study are presented in the introduction. The bulk of the paper describes the evaluation method; the reference system and an advanced transportation system, including a spinning tether in an eccentric Earth orbit, are used to illustrate it.
Validation of a Numerical Method for Determining Liner Impedance
NASA Technical Reports Server (NTRS)
Watson, Willie R.; Jones, Michael G.; Tanner, Sharon E.; Parrott, Tony L.
1996-01-01
This paper reports the initial results of a test series to evaluate a method for determining the normal incidence impedance of a locally reacting acoustically absorbing liner, located on the lower wall of a duct in a grazing incidence, multi-modal, non-progressive acoustic wave environment without flow. This initial evaluation is accomplished by testing the methods' ability to converge to the known normal incidence impedance of a solid steel plate, and to the normal incidence impedance of an absorbing test specimen whose impedance was measured in a conventional normal incidence tube. The method is shown to converge to the normal incident impedance values and thus to be an adequate tool for determining the impedance of specimens in a grazing incidence, multi-modal, nonprogressive acoustic wave environment for a broad range of source frequencies.
NASA Astrophysics Data System (ADS)
Schuster, Jonathan
Infrared (IR) detectors are well established as a vital sensor technology for military, defense and commercial applications. Due to the expense and effort required to fabricate pixel arrays, it is imperative to develop numerical simulation models to perform predictive device simulations which assess device characteristics and design considerations. Towards this end, we have developed a robust three-dimensional (3D) numerical simulation model for IR detector pixel arrays. We used the finite-difference time-domain technique to compute the optical characteristics including the reflectance and the carrier generation rate in the device. Subsequently, we employ the finite element method to solve the drift-diffusion equations to compute the electrical characteristics including the I(V) characteristics, quantum efficiency, crosstalk and modulation transfer function. We use our 3D numerical model to study a new class of detector based on the nBn-architecture. This detector is a unipolar unity-gain barrier device consisting of a narrow-gap absorber layer, a wide-gap barrier layer, and a narrow-gap collector layer. We use our model to study the underlying physics of these devices and to explain the anomalously long lateral collection lengths for photocarriers measured experimentally. Next, we investigate the crosstalk in HgCdTe photovoltaic pixel arrays employing a photon-trapping (PT) structure realized with a periodic array of pillars intended to provide broadband operation. The PT region drastically reduces the crosstalk; making the use of the PT structures not only useful to obtain broadband operation, but also desirable for reducing crosstalk, especially in small pitch detector arrays. Then, the power and flexibility of the nBn architecture is coupled with a PT structure to engineer spectrally filtering detectors. Last, we developed a technique to reduce the cost of large-format, high performance HgCdTe detectors by nondestructively screen-testing detector arrays prior
Advances and future directions of research on spectral methods
NASA Technical Reports Server (NTRS)
Patera, A. T.
1986-01-01
Recent advances in spectral methods are briefly reviewed and characterized with respect to their convergence and computational complexity. Classical finite element and spectral approaches are then compared, and spectral element (or p-type finite element) approximations are introduced. The method is applied to the full Navier-Stokes equations, and examples are given of the application of the technique to several transitional flows. Future directions of research in the field are outlined.
Improved numerical methods for turbulent viscous recirculating flows
NASA Technical Reports Server (NTRS)
Turan, A.; Vandoormaal, J. P.
1988-01-01
The performance of discrete methods for the prediction of fluid flows can be enhanced by improving the convergence rate of solvers and by increasing the accuracy of the discrete representation of the equations of motion. This report evaluates the gains in solver performance that are available when various acceleration methods are applied. Various discretizations are also examined and two are recommended because of their accuracy and robustness. Insertion of the improved discretization and solver accelerator into a TEACH mode, that has been widely applied to combustor flows, illustrates the substantial gains to be achieved.
Numerical solution of 2D-vector tomography problem using the method of approximate inverse
NASA Astrophysics Data System (ADS)
Svetov, Ivan; Maltseva, Svetlana; Polyakova, Anna
2016-08-01
We propose a numerical solution of reconstruction problem of a two-dimensional vector field in a unit disk from the known values of the longitudinal and transverse ray transforms. The algorithm is based on the method of approximate inverse. Numerical simulations confirm that the proposed method yields good results of reconstruction of vector fields.
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.
A finite volume method for numerical grid generation
NASA Astrophysics Data System (ADS)
Beale, S. B.
1999-07-01
A novel method to generate body-fitted grids based on the direct solution for three scalar functions is derived. The solution for scalar variables , and is obtained with a conventional finite volume method based on a physical space formulation. The grid is adapted or re-zoned to eliminate the residual error between the current solution and the desired solution, by means of an implicit grid-correction procedure. The scalar variables are re-mapped and the process is reiterated until convergence is obtained. Calculations are performed for a variety of problems by assuming combined Dirichlet-Neumann and pure Dirichlet boundary conditions involving the use of transcendental control functions, as well as functions designed to effect grid control automatically on the basis of boundary values. The use of dimensional analysis to build stable exponential functions and other control functions is demonstrated. Automatic procedures are implemented: one based on a finite difference approximation to the Cristoffel terms assuming local-boundary orthogonality, and another designed to procure boundary orthogonality. The performance of the new scheme is shown to be comparable with that of conventional inverse methods when calculations are performed on benchmark problems through the application of point-by-point and whole-field solution schemes. Advantages and disadvantages of the present method are critically appraised. Copyright
Evaluating numerical ODE/DAE methods, algorithms and software
NASA Astrophysics Data System (ADS)
Soderlind, Gustaf; Wang, Lina
2006-01-01
Until recently, the testing of ODE/DAE software has been limited to simple comparisons and benchmarking. The process of developing software from a mathematically specified method is complex: it entails constructing control structures and objectives, selecting iterative methods and termination criteria, choosing norms and many more decisions. Most software constructors have taken a heuristic approach to these design choices, and as a consequence two different implementations of the same method may show significant differences in performance. Yet it is common to try to deduce from software comparisons that one method is better than another. Such conclusions are not warranted, however, unless the testing is carried out under true ceteris paribus conditions. Moreover, testing is an empirical science and as such requires a formal test protocol; without it conclusions are questionable, invalid or even false.We argue that ODE/DAE software can be constructed and analyzed by proven, "standard" scientific techniques instead of heuristics. The goals are computational stability, reproducibility, and improved software quality. We also focus on different error criteria and norms, and discuss modifications to DASPK and RADAU5. Finally, some basic principles of a test protocol are outlined and applied to testing these codes on a variety of problems.
Advances in subtyping methods of foodborne disease pathogens.
Boxrud, Dave
2010-04-01
Current subtyping methods for the detection of foodborne disease outbreaks have limitations that reduce their use by public health laboratories. Recent advances in subtyping of foodborne disease pathogens utilize techniques that identify nucleic acid polymorphisms. Recent methods of nucleic acid characterization such as microarrays and mass spectrometry (MS) may provide improvements such as increasing speed and data portability while decreasing labor compared to current methods. This article discusses multiple-locus variable-number tandem-repeat analysis, single-nucleotide polymorphisms, nucleic acid sequencing, whole genome sequencing, variable absent or present loci, microarrays and MS as potential subtyping methods to enhance our ability to detect foodborne disease outbreaks. PMID:20299203
Applications of numerical methods to simulate the movement of contaminants in groundwater.
Sun, N Z
1989-01-01
This paper reviews mathematical models and numerical methods that have been extensively used to simulate the movement of contaminants through the subsurface. The major emphasis is placed on the numerical methods of advection-dominated transport problems and inverse problems. Several mathematical models that are commonly used in field problems are listed. A variety of numerical solutions for three-dimensional models are introduced, including the multiple cell balance method that can be considered a variation of the finite element method. The multiple cell balance method is easy to understand and convenient for solving field problems. When the advection transport dominates the dispersion transport, two kinds of numerical difficulties, overshoot and numerical dispersion, are always involved in solving standard, finite difference methods and finite element methods. To overcome these numerical difficulties, various numerical techniques are developed, such as upstream weighting methods and moving point methods. A complete review of these methods is given and we also mention the problems of parameter identification, reliability analysis, and optimal-experiment design that are absolutely necessary for constructing a practical model. PMID:2695327
NASA Technical Reports Server (NTRS)
Srinivasan, G. R.; Baeder, J. D.
1991-01-01
This paper outlines some recent advances in the application of the Euler and Navier-Stokes computational fluid dynamics methods to analyze nonlinear problems of helicopter aerodynamics and acoustics. A complete flowfield simulation of helicopters is currently not feasible with these methods. However, the use of the state-of-the-art numerical algorithms in conjunction with powerful supercomputers, like the Cray-2, have enabled notable progress to be made in modeling several individual components of this complex flow in hover and forward flight.