ERIC Educational Resources Information Center
Sozio, Gerry
2009-01-01
Senior secondary students cover numerical integration techniques in their mathematics courses. In particular, students would be familiar with the "midpoint rule," the elementary "trapezoidal rule" and "Simpson's rule." This article derives these techniques by methods which secondary students may not be familiar with and an approach that…
NASA Astrophysics Data System (ADS)
Khabarov, Nikolay; Huggel, Christian; Obersteiner, Michael; Ramírez, Juan Manuel
2010-05-01
Mountain regions are typically characterized by rugged terrain which is susceptible to different types of landslides during high-intensity precipitation. Landslides account for billions of dollars of damage and many casualties, and are expected to increase in frequency in the future due to a projected increase of precipitation intensity. Early warning systems (EWS) are thought to be a primary tool for related disaster risk reduction and climate change adaptation to extreme climatic events and hydro-meteorological hazards, including landslides. An EWS for hazards such as landslides consist of different components, including environmental monitoring instruments (e.g. rainfall or flow sensors), physical or empirical process models to support decision-making (warnings, evacuation), data and voice communication, organization and logistics-related procedures, and population response. Considering this broad range, EWS are highly complex systems, and it is therefore difficult to understand the effect of the different components and changing conditions on the overall performance, ultimately being expressed as human lives saved or structural damage reduced. In this contribution we present a further development of our approach to assess a landslide EWS in an integral way, both at the system and component level. We utilize a numerical model using 6 hour rainfall data as basic input. A threshold function based on a rainfall-intensity/duration relation was applied as a decision criterion for evacuation. Damage to infrastructure and human lives was defined as a linear function of landslide magnitude, with the magnitude modelled using a power function of landslide frequency. Correct evacuation was assessed with a ‘true' reference rainfall dataset versus a dataset of artificially reduced quality imitating the observation system component. Performance of the EWS using these rainfall datasets was expressed in monetary terms (i.e. damage related to false and correct evacuation). We
Rythmos Numerical Integration Package
Coffey, Todd S.; Bartlett, Roscoe A.
2006-09-01
Rythmos numerically integrates transient differential equations. The differential equations can be explicit or implicit ordinary differential equations ofr formulated as fully implicit differential-algebraic equations. Methods include backward Euler, forward Euler, explicit Runge-Kutta, and implicit BDF at this time. Native support for operator split methods and strict modularity are strong design goals. Forward sensitivity computations will be included in the first release with adjoint sensitivities coming in the near future. Rythmos heavily relies on Thyra for linear algebra and nonlinear solver interfaces to AztecOO, Amesos, IFPack, and NOX in Tilinos. Rythmos is specially suited for stiff differential equations and thos applictions where operator split methods have a big advantage, e.g. Computational fluid dynamics, convection-diffusion equations, etc.
Rythmos Numerical Integration Package
Energy Science and Technology Software Center (ESTSC)
2006-09-01
Rythmos numerically integrates transient differential equations. The differential equations can be explicit or implicit ordinary differential equations ofr formulated as fully implicit differential-algebraic equations. Methods include backward Euler, forward Euler, explicit Runge-Kutta, and implicit BDF at this time. Native support for operator split methods and strict modularity are strong design goals. Forward sensitivity computations will be included in the first release with adjoint sensitivities coming in the near future. Rythmos heavily relies on Thyra formore » linear algebra and nonlinear solver interfaces to AztecOO, Amesos, IFPack, and NOX in Tilinos. Rythmos is specially suited for stiff differential equations and thos applictions where operator split methods have a big advantage, e.g. Computational fluid dynamics, convection-diffusion equations, etc.« less
Accelerated adaptive integration method.
Kaus, Joseph W; Arrar, Mehrnoosh; McCammon, J Andrew
2014-05-15
Conformational changes that occur upon ligand binding may be too slow to observe on the time scales routinely accessible using molecular dynamics simulations. The adaptive integration method (AIM) leverages the notion that when a ligand is either fully coupled or decoupled, according to λ, barrier heights may change, making some conformational transitions more accessible at certain λ values. AIM adaptively changes the value of λ in a single simulation so that conformations sampled at one value of λ seed the conformational space sampled at another λ value. Adapting the value of λ throughout a simulation, however, does not resolve issues in sampling when barriers remain high regardless of the λ value. In this work, we introduce a new method, called Accelerated AIM (AcclAIM), in which the potential energy function is flattened at intermediate values of λ, promoting the exploration of conformational space as the ligand is decoupled from its receptor. We show, with both a simple model system (Bromocyclohexane) and the more complex biomolecule Thrombin, that AcclAIM is a promising approach to overcome high barriers in the calculation of free energies, without the need for any statistical reweighting or additional processors. PMID:24780083
Accelerated Adaptive Integration Method
2015-01-01
Conformational changes that occur upon ligand binding may be too slow to observe on the time scales routinely accessible using molecular dynamics simulations. The adaptive integration method (AIM) leverages the notion that when a ligand is either fully coupled or decoupled, according to λ, barrier heights may change, making some conformational transitions more accessible at certain λ values. AIM adaptively changes the value of λ in a single simulation so that conformations sampled at one value of λ seed the conformational space sampled at another λ value. Adapting the value of λ throughout a simulation, however, does not resolve issues in sampling when barriers remain high regardless of the λ value. In this work, we introduce a new method, called Accelerated AIM (AcclAIM), in which the potential energy function is flattened at intermediate values of λ, promoting the exploration of conformational space as the ligand is decoupled from its receptor. We show, with both a simple model system (Bromocyclohexane) and the more complex biomolecule Thrombin, that AcclAIM is a promising approach to overcome high barriers in the calculation of free energies, without the need for any statistical reweighting or additional processors. PMID:24780083
Adaptive Numerical Algorithms in Space Weather Modeling
NASA Technical Reports Server (NTRS)
Toth, Gabor; vanderHolst, Bart; Sokolov, Igor V.; DeZeeuw, Darren; Gombosi, Tamas I.; Fang, Fang; Manchester, Ward B.; Meng, Xing; Nakib, Dalal; Powell, Kenneth G.; Stout, Quentin F.; Glocer, Alex; Ma, Ying-Juan; Opher, Merav
2010-01-01
Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different physics in different domains. A multi-physics system can be modeled by a software framework comprising of several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solar wind Roe Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamics (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical
Adaptive numerical algorithms in space weather modeling
NASA Astrophysics Data System (ADS)
Tóth, Gábor; van der Holst, Bart; Sokolov, Igor V.; De Zeeuw, Darren L.; Gombosi, Tamas I.; Fang, Fang; Manchester, Ward B.; Meng, Xing; Najib, Dalal; Powell, Kenneth G.; Stout, Quentin F.; Glocer, Alex; Ma, Ying-Juan; Opher, Merav
2012-02-01
Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different relevant physics in different domains. A multi-physics system can be modeled by a software framework comprising several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solarwind Roe-type Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamic (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit
Numerical integration of subtraction terms
NASA Astrophysics Data System (ADS)
Seth, Satyajit; Weinzierl, Stefan
2016-06-01
Numerical approaches to higher-order calculations often employ subtraction terms, both for the real emission and the virtual corrections. These subtraction terms have to be added back. In this paper we show that at NLO the real subtraction terms, the virtual subtraction terms, the integral representations of the field renormalization constants and—in the case of initial-state partons—the integral representation for the collinear counterterm can be grouped together to give finite integrals, which can be evaluated numerically. This is useful for an extension towards next-to-next-to-leading order.
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.
Fresnel Integral Equations: Numerical Properties
Adams, R J; Champagne, N J II; Davis, B A
2003-07-22
A spatial-domain solution to the problem of electromagnetic scattering from a dielectric half-space is outlined. The resulting half-space operators are referred to as Fresnel surface integral operators. When used as preconditioners for nonplanar geometries, the Fresnel operators yield surface Fresnel integral equations (FIEs) which are stable with respect to dielectric constant, discretization, and frequency. Numerical properties of the formulations are discussed.
Cuba: Multidimensional numerical integration library
NASA Astrophysics Data System (ADS)
Hahn, Thomas
2016-08-01
The Cuba library offers four independent routines for multidimensional numerical integration: Vegas, Suave, Divonne, and Cuhre. The four algorithms work by very different methods, and can integrate vector integrands and have very similar Fortran, C/C++, and Mathematica interfaces. Their invocation is very similar, making it easy to cross-check by substituting one method by another. For further safeguarding, the output is supplemented by a chi-square probability which quantifies the reliability of the error estimate.
Recursive adaptive frame integration limited
NASA Astrophysics Data System (ADS)
Rafailov, Michael K.
2006-05-01
Recursive Frame Integration Limited was proposed as a way to improve frame integration performance and mitigate issues related to high data rate needed for conventional frame integration. The technique applies two thresholds - one tuned for optimum probability of detection, the other to manage required false alarm rate - and allows a non-linear integration process that, along with Signal-to-Noise Ratio (SNR) gain, provides system designers more capability where cost, weight, or power considerations limit system data rate, processing, or memory capability. However, Recursive Frame Integration Limited may have performance issues when single frame SNR is really low. Recursive Adaptive Frame Integration Limited is proposed as a means to improve limited integration performance with really low single frame SNR. It combines the benefits of nonlinear recursive limited frame integration and adaptive thresholds with a kind of conventional frame integration.
Adaptive Urban Dispersion Integrated Model
Wissink, A; Chand, K; Kosovic, B; Chan, S; Berger, M; Chow, F K
2005-11-03
Numerical simulations represent a unique predictive tool for understanding the three-dimensional flow fields and associated concentration distributions from contaminant releases in complex urban settings (Britter and Hanna 2003). Utilization of the most accurate urban models, based on fully three-dimensional computational fluid dynamics (CFD) that solve the Navier-Stokes equations with incorporated turbulence models, presents many challenges. We address two in this work; first, a fast but accurate way to incorporate the complex urban terrain, buildings, and other structures to enforce proper boundary conditions in the flow solution; second, ways to achieve a level of computational efficiency that allows the models to be run in an automated fashion such that they may be used for emergency response and event reconstruction applications. We have developed a new integrated urban dispersion modeling capability based on FEM3MP (Gresho and Chan 1998, Chan and Stevens 2000), a CFD model from Lawrence Livermore National Lab. The integrated capability incorporates fast embedded boundary mesh generation for geometrically complex problems and full three-dimensional Cartesian adaptive mesh refinement (AMR). Parallel AMR and embedded boundary gridding support are provided through the SAMRAI library (Wissink et al. 2001, Hornung and Kohn 2002). Embedded boundary mesh generation has been demonstrated to be an automatic, fast, and efficient approach for problem setup. It has been used for a variety of geometrically complex applications, including urban applications (Pullen et al. 2005). The key technology we introduce in this work is the application of AMR, which allows the application of high-resolution modeling to certain important features, such as individual buildings and high-resolution terrain (including important vegetative and land-use features). It also allows the urban scale model to be readily interfaced with coarser resolution meso or regional scale models. This talk
Numerical design of an adaptive aileron
NASA Astrophysics Data System (ADS)
Amendola, Gianluca; Dimino, Ignazio; Concilio, Antonio; Magnifico, Marco; Pecora, Rosario
2016-04-01
The study herein described is aimed at investigating the feasibility of an innovative full-scale camber morphing aileron device. In the framework of the "Adaptive Aileron" project, an international cooperation between Italy and Canada, this goal was carried out with the integration of different morphing concepts in a wing-tip prototype. As widely demonstrated in recent European projects such as Clean Sky JTI and SARISTU, wing trailing edge morphing may lead to significant drag reduction (up to 6%) in off-design flight points by adapting chord-wise camber variations in cruise to compensate A/C weight reduction following fuel consumption. Those researches focused on the flap region as the most immediate solution to implement structural adaptations. However, there is also a growing interest in extending morphing functionalities to the aileron region preserving its main functionality in controlling aircraft directional stability. In fact, the external region of the wing seems to be the most effective in producing "lift over drag" improvements by morphing. Thus, the objective of the presented research is to achieve a certain drag reduction in off-design flight points by adapting wing shape and lift distribution following static deflections. In perspective, the developed device could also be used as a load alleviation system to reduce gust effects, augmenting its frequency bandwidth. In this paper, the preliminary design of the adaptive aileron is first presented, assessed on the base of the external aerodynamic loads. The primary structure is made of 5 segmented ribs, distributed along 4 bays, each splitted into three consecutive parts, connected with spanwise stringers. The aileron shape modification is then implemented by means of an actuation system, based on a classical quick-return mechanism, opportunely suited for the presented application. Finite element analyses were assessed for properly sizing the load-bearing structure and actuation systems and for
GRChombo: Numerical relativity with adaptive mesh refinement
NASA Astrophysics Data System (ADS)
Clough, Katy; Figueras, Pau; Finkel, Hal; Kunesch, Markus; Lim, Eugene A.; Tunyasuvunakool, Saran
2015-12-01
In this work, we introduce {\\mathtt{GRChombo}}: a new numerical relativity code which incorporates full adaptive mesh refinement (AMR) using block structured Berger-Rigoutsos grid generation. The code supports non-trivial 'many-boxes-in-many-boxes' mesh hierarchies and massive parallelism through the message passing interface. {\\mathtt{GRChombo}} evolves the Einstein equation using the standard BSSN formalism, with an option to turn on CCZ4 constraint damping if required. The AMR capability permits the study of a range of new physics which has previously been computationally infeasible in a full 3 + 1 setting, while also significantly simplifying the process of setting up the mesh for these problems. We show that {\\mathtt{GRChombo}} can stably and accurately evolve standard spacetimes such as binary black hole mergers and scalar collapses into black holes, demonstrate the performance characteristics of our code, and discuss various physics problems which stand to benefit from the AMR technique.
Numerical Integration: One Step at a Time
ERIC Educational Resources Information Center
Yang, Yajun; Gordon, Sheldon P.
2016-01-01
This article looks at the effects that adding a single extra subdivision has on the level of accuracy of some common numerical integration routines. Instead of automatically doubling the number of subdivisions for a numerical integration rule, we investigate what happens with a systematic method of judiciously selecting one extra subdivision for…
DANA: distributed numerical and adaptive modelling framework.
Rougier, Nicolas P; Fix, Jérémy
2012-01-01
DANA is a python framework ( http://dana.loria.fr ) whose computational paradigm is grounded on the notion of a unit that is essentially a set of time dependent values varying under the influence of other units via adaptive weighted connections. The evolution of a unit's value are defined by a set of differential equations expressed in standard mathematical notation which greatly ease their definition. The units are organized into groups that form a model. Each unit can be connected to any other unit (including itself) using a weighted connection. The DANA framework offers a set of core objects needed to design and run such models. The modeler only has to define the equations of a unit as well as the equations governing the training of the connections. The simulation is completely transparent to the modeler and is handled by DANA. This allows DANA to be used for a wide range of numerical and distributed models as long as they fit the proposed framework (e.g. cellular automata, reaction-diffusion system, decentralized neural networks, recurrent neural networks, kernel-based image processing, etc.). PMID:22994650
An Integrative Theory of Numerical Development
ERIC Educational Resources Information Center
Siegler, Robert; Lortie-Forgues, Hugues
2014-01-01
Understanding of numerical development is growing rapidly, but the volume and diversity of findings can make it difficult to perceive any coherence in the process. The integrative theory of numerical development posits that a coherent theme is present, however--progressive broadening of the set of numbers whose magnitudes can be accurately…
Adaptive numerical methods for partial differential equations
Cololla, P.
1995-07-01
This review describes a structured approach to adaptivity. The Automated Mesh Refinement (ARM) algorithms developed by M Berger are described, touching on hyperbolic and parabolic applications. Adaptivity is achieved by overlaying finer grids only in areas flagged by a generalized error criterion. The author discusses some of the issues involved in abutting disparate-resolution grids, and demonstrates that suitable algorithms exist for dissipative as well as hyperbolic systems.
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
Adapting Inspection Data for Computer Numerical Control
NASA Technical Reports Server (NTRS)
Hutchison, E. E.
1986-01-01
Machining time for repetitive tasks reduced. Program converts measurements of stub post locations by coordinate-measuring machine into form used by numerical-control computer. Work time thus reduced by 10 to 15 minutes for each post. Since there are 600 such posts on each injector, time saved per injector is 100 to 150 hours. With modifications this approach applicable to machining of many precise holes on large machine frames and similar objects.
Effects of aliasing on numerical integration.
Edwards, Timothy S.
2005-02-01
During the course of processing acceleration data from mechanical systems it is often desirable to integrate the data to obtain velocity or displacement waveforms. However, those who have attempted these operations may be painfully aware that the integrated records often yield unrealistic residual values. This is true whether the data has been obtained experimentally or through numerical simulation such as Runge-Kutta integration or the explicit finite element method. In the case of experimentally obtained data, the integration errors are usually blamed on accelerometer zero shift or amplifier saturation. In the case of simulation data, incorrect integrations are often incorrectly blamed on the integration algorithm itself. This work demonstrates that seemingly small aliased content can cause appreciable errors in the integrated waveforms and explores the unavoidable source of aliasing in both experiment and simulation-the sampling operation. Numerical analysts are often puzzled as to why the integrated acceleration from their simulation does not match the displacement output from the same simulation. This work shows that these strange results can be caused by aliasing induced by interpolation of the model output during sampling regularization.
Highly Parallel, High-Precision Numerical Integration
Bailey, David H.; Borwein, Jonathan M.
2005-04-22
This paper describes a scheme for rapidly computing numerical values of definite integrals to very high accuracy, ranging from ordinary machine precision to hundreds or thousands of digits, even for functions with singularities or infinite derivatives at endpoints. Such a scheme is of interest not only in computational physics and computational chemistry, but also in experimental mathematics, where high-precision numerical values of definite integrals can be used to numerically discover new identities. This paper discusses techniques for a parallel implementation of this scheme, then presents performance results for 1-D and 2-D test suites. Results are also given for a certain problem from mathematical physics, which features a difficult singularity, confirming a conjecture to 20,000 digit accuracy. The performance rate for this latter calculation on 1024 CPUs is 690 Gflop/s. We believe that this and one other 20,000-digit integral evaluation that we report are the highest-precision non-trivial numerical integrations performed to date.
Numerical multi-loop integrals and applications
NASA Astrophysics Data System (ADS)
Freitas, A.
2016-09-01
Higher-order radiative corrections play an important role in precision studies of the electroweak and Higgs sector, as well as for the detailed understanding of large backgrounds to new physics searches. For corrections beyond the one-loop level and involving many independent mass and momentum scales, it is in general not possible to find analytic results, so that one needs to resort to numerical methods instead. This article presents an overview of a variety of numerical loop integration techniques, highlighting their range of applicability, suitability for automatization, and numerical precision and stability. In a second part of this article, the application of numerical loop integration methods in the area of electroweak precision tests is illustrated. Numerical methods were essential for obtaining full two-loop predictions for the most important precision observables within the Standard Model. The theoretical foundations for these corrections will be described in some detail, including aspects of the renormalization, resummation of leading log contributions, and the evaluation of the theory uncertainty from missing higher orders.
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.
Efficient numerical evaluation of Feynman integrals
NASA Astrophysics Data System (ADS)
Li, Zhao; Wang, Jian; Yan, Qi-Shu; Zhao, Xiaoran
2016-03-01
Feynman loop integrals are a key ingredient for the calculation of higher order radiation effects, and are responsible for reliable and accurate theoretical prediction. We improve the efficiency of numerical integration in sector decomposition by implementing a quasi-Monte Carlo method associated with the CUDA/GPU technique. For demonstration we present the results of several Feynman integrals up to two loops in both Euclidean and physical kinematic regions in comparison with those obtained from FIESTA3. It is shown that both planar and non-planar two-loop master integrals in the physical kinematic region can be evaluated in less than half a minute with accuracy, which makes the direct numerical approach viable for precise investigation of higher order effects in multi-loop processes, e.g. the next-to-leading order QCD effect in Higgs pair production via gluon fusion with a finite top quark mass. Supported by the Natural Science Foundation of China (11305179 11475180), Youth Innovation Promotion Association, CAS, IHEP Innovation (Y4545170Y2), State Key Lab for Electronics and Particle Detectors, Open Project Program of State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, China (Y4KF061CJ1), Cluster of Excellence Precision Physics, Fundamental Interactions and Structure of Matter (PRISMA-EXC 1098)
Integrated optical circuits for numerical computation
NASA Technical Reports Server (NTRS)
Verber, C. M.; Kenan, R. P.
1983-01-01
The development of integrated optical circuits (IOC) for numerical-computation applications is reviewed, with a focus on the use of systolic architectures. The basic architecture criteria for optical processors are shown to be the same as those proposed by Kung (1982) for VLSI design, and the advantages of IOCs over bulk techniques are indicated. The operation and fabrication of electrooptic grating structures are outlined, and the application of IOCs of this type to an existing 32-bit, 32-Mbit/sec digital correlator, a proposed matrix multiplier, and a proposed pipeline processor for polynomial evaluation is discussed. The problems arising from the inherent nonlinearity of electrooptic gratings are considered. Diagrams and drawings of the application concepts are provided.
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 integration of asymptotic solutions of ordinary differential equations
NASA Technical Reports Server (NTRS)
Thurston, Gaylen A.
1989-01-01
Classical asymptotic analysis of ordinary differential equations derives approximate solutions that are numerically stable. However, the analysis also leads to tedious expansions in powers of the relevant parameter for a particular problem. The expansions are replaced with integrals that can be evaluated by numerical integration. The resulting numerical solutions retain the linear independence that is the main advantage of asymptotic solutions. Examples, including the Falkner-Skan equation from laminar boundary layer theory, illustrate the method of asymptotic analysis with numerical integration.
Ying, Wenjun; Henriquez, Craig S.
2015-01-01
A both space and time adaptive algorithm is presented for simulating electrical wave propagation in the Purkinje system of the heart. The equations governing the distribution of electric potential over the system are solved in time with the method of lines. At each timestep, by an operator splitting technique, the space-dependent but linear diffusion part and the nonlinear but space-independent reactions part in the partial differential equations are integrated separately with implicit schemes, which have better stability and allow larger timesteps than explicit ones. The linear diffusion equation on each edge of the system is spatially discretized with the continuous piecewise linear finite element method. The adaptive algorithm can automatically recognize when and where the electrical wave starts to leave or enter the computational domain due to external current/voltage stimulation, self-excitation, or local change of membrane properties. Numerical examples demonstrating efficiency and accuracy of the adaptive algorithm are presented. PMID:26581455
Numerical simulation of immiscible viscous fingering using adaptive unstructured meshes
NASA Astrophysics Data System (ADS)
Adam, A.; Salinas, P.; Percival, J. R.; Pavlidis, D.; Pain, C.; Muggeridge, A. H.; Jackson, M.
2015-12-01
Displacement of one fluid by another in porous media occurs in various settings including hydrocarbon recovery, CO2 storage and water purification. When the invading fluid is of lower viscosity than the resident fluid, the displacement front is subject to a Saffman-Taylor instability and is unstable to transverse perturbations. These instabilities can grow, leading to fingering of the invading fluid. Numerical simulation of viscous fingering is challenging. The physics is controlled by a complex interplay of viscous and diffusive forces and it is necessary to ensure physical diffusion dominates numerical diffusion to obtain converged solutions. This typically requires the use of high mesh resolution and high order numerical methods. This is computationally expensive. We demonstrate here the use of a novel control volume - finite element (CVFE) method along with dynamic unstructured mesh adaptivity to simulate viscous fingering with higher accuracy and lower computational cost than conventional methods. Our CVFE method employs a discontinuous representation for both pressure and velocity, allowing the use of smaller control volumes (CVs). This yields higher resolution of the saturation field which is represented CV-wise. Moreover, dynamic mesh adaptivity allows high mesh resolution to be employed where it is required to resolve the fingers and lower resolution elsewhere. We use our results to re-examine the existing criteria that have been proposed to govern the onset of instability.Mesh adaptivity requires the mapping of data from one mesh to another. Conventional methods such as consistent interpolation do not readily generalise to discontinuous fields and are non-conservative. We further contribute a general framework for interpolation of CV fields by Galerkin projection. The method is conservative, higher order and yields improved results, particularly with higher order or discontinuous elements where existing approaches are often excessively diffusive.
Numerical simulations of the LBT adaptive secondary mirror
NASA Astrophysics Data System (ADS)
Del Vecchio, Ciro; Gallieni, Daniele
2000-07-01
In this paper we describe the design of the deformable mirror of the Large Binocular Telescope adaptive secondary unit. Starting from the optical design, a numerical model of the ultra-thin, aspherical glass shell, accommodating the 918 magnets on the selected actuator geometry, has been run. As a result, we can evaluate the response of this crucial component of the telescope optics with great accuracy. The DM is analyzed from the mechanical standpoint -- gravity deformations, wavefront residue, residue of low-order Zernike aberrations, corrections of magnetic interactions -- in order to compute the optical performances in the most demanding operational circumstances.
Fuzzy Adaptive Cubature Kalman Filter for Integrated Navigation Systems.
Tseng, Chien-Hao; Lin, Sheng-Fuu; Jwo, Dah-Jing
2016-01-01
This paper presents a sensor fusion method based on the combination of cubature Kalman filter (CKF) and fuzzy logic adaptive system (FLAS) for the integrated navigation systems, such as the GPS/INS (Global Positioning System/inertial navigation system) integration. The third-degree spherical-radial cubature rule applied in the CKF has been employed to avoid the numerically instability in the system model. In processing navigation integration, the performance of nonlinear filter based estimation of the position and velocity states may severely degrade caused by modeling errors due to dynamics uncertainties of the vehicle. In order to resolve the shortcoming for selecting the process noise covariance through personal experience or numerical simulation, a scheme called the fuzzy adaptive cubature Kalman filter (FACKF) is presented by introducing the FLAS to adjust the weighting factor of the process noise covariance matrix. The FLAS is incorporated into the CKF framework as a mechanism for timely implementing the tuning of process noise covariance matrix based on the information of degree of divergence (DOD) parameter. The proposed FACKF algorithm shows promising accuracy improvement as compared to the extended Kalman filter (EKF), unscented Kalman filter (UKF), and CKF approaches. PMID:27472336
Numerical modeling of seismic waves using frequency-adaptive meshes
NASA Astrophysics Data System (ADS)
Hu, Jinyin; Jia, Xiaofeng
2016-08-01
An improved modeling algorithm using frequency-adaptive meshes is applied to meet the computational requirements of all seismic frequency components. It automatically adopts coarse meshes for low-frequency computations and fine meshes for high-frequency computations. The grid intervals are adaptively calculated based on a smooth inversely proportional function of grid size with respect to the frequency. In regular grid-based methods, the uniform mesh or non-uniform mesh is used for frequency-domain wave propagators and it is fixed for all frequencies. A too coarse mesh results in inaccurate high-frequency wavefields and unacceptable numerical dispersion; on the other hand, an overly fine mesh may cause storage and computational overburdens as well as invalid propagation angles of low-frequency wavefields. Experiments on the Padé generalized screen propagator indicate that the Adaptive mesh effectively solves these drawbacks of regular fixed-mesh methods, thus accurately computing the wavefield and its propagation angle in a wide frequency band. Several synthetic examples also demonstrate its feasibility for seismic modeling and migration.
Systems integration of innate and adaptive immunity.
Zak, Daniel E; Aderem, Alan
2015-09-29
The pathogens causing AIDS, malaria, and tuberculosis have proven too complex to be overcome by classical approaches to vaccination. The complexities of human immunology and pathogen-induced modulation of the immune system mandate new approaches to vaccine discovery and design. A new field, systems vaccinology, weds holistic analysis of innate and adaptive immunity within a quantitative framework to enable rational design of new vaccines that elicit tailored protective immune responses. A key step in the approach is to discover relationships between the earliest innate inflammatory responses to vaccination and the subsequent vaccine-induced adaptive immune responses and efficacy. Analysis of these responses in clinical studies is complicated by the inaccessibility of relevant tissue compartments (such as the lymph node), necessitating reliance upon peripheral blood responses as surrogates. Blood transcriptomes, although indirect to vaccine mechanisms, have proven very informative in systems vaccinology studies. The approach is most powerful when innate and adaptive immune responses are integrated with vaccine efficacy, which is possible for malaria with the advent of a robust human challenge model. This is more difficult for AIDS and tuberculosis, given that human challenge models are lacking and efficacy observed in clinical trials has been low or highly variable. This challenge can be met by appropriate clinical trial design for partially efficacious vaccines and by analysis of natural infection cohorts. Ultimately, systems vaccinology is an iterative approach in which mechanistic hypotheses-derived from analysis of clinical studies-are evaluated in model systems, and then used to guide the development of new vaccine strategies. In this review, we will illustrate the above facets of the systems vaccinology approach with case studies. PMID:26102534
Deb, M.K.; Kennon, S.R.
1998-04-01
A cooperative R&D effort between industry and the US government, this project, under the HPPP (High Performance Parallel Processing) initiative of the Dept. of Energy, started the investigations into parallel object-oriented (OO) numerics. The basic goal was to research and utilize the emerging technologies to create a physics-independent computational kernel for applications using adaptive finite element method. The industrial team included Computational Mechanics Co., Inc. (COMCO) of Austin, TX (as the primary contractor), Scientific Computing Associates, Inc. (SCA) of New Haven, CT, Texaco and CONVEX. Sandia National Laboratory (Albq., NM) was the technology partner from the government side. COMCO had the responsibility of the main kernel design and development, SCA had the lead in parallel solver technology and guidance on OO technologies was Sandia`s main expertise in this venture. CONVEX and Texaco supported the partnership by hardware resource and application knowledge, respectively. As such, a minimum of fifty-percent cost-sharing was provided by the industry partnership during this project. This report describes the R&D activities and provides some details about the prototype kernel and example applications.
Numerical study of Taylor bubbles with adaptive unstructured meshes
NASA Astrophysics Data System (ADS)
Xie, Zhihua; Pavlidis, Dimitrios; Percival, James; Pain, Chris; Matar, Omar; Hasan, Abbas; Azzopardi, Barry
2014-11-01
The Taylor bubble is a single long bubble which nearly fills the entire cross section of a liquid-filled circular tube. This type of bubble flow regime often occurs in gas-liquid slug flows in many industrial applications, including oil-and-gas production, chemical and nuclear reactors, and heat exchangers. The objective of this study is to investigate the fluid dynamics of Taylor bubbles rising in a vertical pipe filled with oils of extremely high viscosity (mimicking the ``heavy oils'' found in the oil-and-gas industry). A modelling and simulation framework is presented here which can modify and adapt anisotropic unstructured meshes to better represent the underlying physics of bubble rise and reduce the computational effort without sacrificing accuracy. The numerical framework consists of a mixed control-volume and finite-element formulation, a ``volume of fluid''-type method for the interface capturing based on a compressive control volume advection method, and a force-balanced algorithm for the surface tension implementation. Numerical examples of some benchmark tests and the dynamics of Taylor bubbles are presented to show the capability of this method. EPSRC Programme Grant, MEMPHIS, EP/K0039761/1.
Error Estimates for Numerical Integration Rules
ERIC Educational Resources Information Center
Mercer, Peter R.
2005-01-01
The starting point for this discussion of error estimates is the fact that integrals that arise in Fourier series have properties that can be used to get improved bounds. This idea is extended to more general situations.
Numerical evaluation of Feynman path integrals
NASA Astrophysics Data System (ADS)
Baird, William Hugh
1999-11-01
The notion of path integration developed by Feynman, while an incredibly successful method of solving quantum mechanical problems, leads to frequently intractable integrations over an infinite number of paths. Two methods now exist which sidestep this difficulty by defining "densities" of actions which give the relative number of paths found at different values of the action. These densities are sampled by computer generation of paths and the propagators are found to a high degree of accuracy for the case of a particle on the infinite half line and in a finite square well in one dimension. The problem of propagation within a two dimensional radial well is also addressed as the precursor to the problem of a particle in a stadium (quantum billiard).
Adaptive integral dynamic surface control of a hypersonic flight vehicle
NASA Astrophysics Data System (ADS)
Aslam Butt, Waseem; Yan, Lin; Amezquita S., Kendrick
2015-07-01
In this article, non-linear adaptive dynamic surface air speed and flight path angle control designs are presented for the longitudinal dynamics of a flexible hypersonic flight vehicle. The tracking performance of the control design is enhanced by introducing a novel integral term that caters to avoiding a large initial control signal. To ensure feasibility, the design scheme incorporates magnitude and rate constraints on the actuator commands. The uncertain non-linear functions are approximated by an efficient use of the neural networks to reduce the computational load. A detailed stability analysis shows that all closed-loop signals are uniformly ultimately bounded and the ? tracking performance is guaranteed. The robustness of the design scheme is verified through numerical simulations of the flexible flight vehicle model.
Numerical integration of ordinary differential equations of various orders
NASA Technical Reports Server (NTRS)
Gear, C. W.
1969-01-01
Report describes techniques for the numerical integration of differential equations of various orders. Modified multistep predictor-corrector methods for general initial-value problems are discussed and new methods are introduced.
A Numerical Study of Hypersonic Forebody/Inlet Integration Problem
NASA Technical Reports Server (NTRS)
Kumar, Ajay
1991-01-01
A numerical study of hypersonic forebody/inlet integration problem is presented in the form of the view-graphs. The following topics are covered: physical/chemical modeling; solution procedure; flow conditions; mass flow rate at inlet face; heating and skin friction loads; 3-D forebogy/inlet integration model; and sensitivity studies.
Translation and integration of numerical atomic orbitals in linear molecules.
Heinäsmäki, Sami
2014-02-14
We present algorithms for translation and integration of atomic orbitals for LCAO calculations in linear molecules. The method applies to arbitrary radial functions given on a numerical mesh. The algorithms are based on pseudospectral differentiation matrices in two dimensions and the corresponding two-dimensional Gaussian quadratures. As a result, multicenter overlap and Coulomb integrals can be evaluated effectively. PMID:24527905
Translation and integration of numerical atomic orbitals in linear molecules
NASA Astrophysics Data System (ADS)
Heinäsmäki, Sami
2014-02-01
We present algorithms for translation and integration of atomic orbitals for LCAO calculations in linear molecules. The method applies to arbitrary radial functions given on a numerical mesh. The algorithms are based on pseudospectral differentiation matrices in two dimensions and the corresponding two-dimensional Gaussian quadratures. As a result, multicenter overlap and Coulomb integrals can be evaluated effectively.
Numerical solution of boundary-integral equations for molecular electrostatics.
Bardhan, Jaydeep P
2009-03-01
Numerous molecular processes, such as ion permeation through channel proteins, are governed by relatively small changes in energetics. As a result, theoretical investigations of these processes require accurate numerical methods. In the present paper, we evaluate the accuracy of two approaches to simulating boundary-integral equations for continuum models of the electrostatics of solvation. The analysis emphasizes boundary-element method simulations of the integral-equation formulation known as the apparent-surface-charge (ASC) method or polarizable-continuum model (PCM). In many numerical implementations of the ASC/PCM model, one forces the integral equation to be satisfied exactly at a set of discrete points on the boundary. We demonstrate in this paper that this approach to discretization, known as point collocation, is significantly less accurate than an alternative approach known as qualocation. Furthermore, the qualocation method offers this improvement in accuracy without increasing simulation time. Numerical examples demonstrate that electrostatic part of the solvation free energy, when calculated using the collocation and qualocation methods, can differ significantly; for a polypeptide, the answers can differ by as much as 10 kcal/mol (approximately 4% of the total electrostatic contribution to solvation). The applicability of the qualocation discretization to other integral-equation formulations is also discussed, and two equivalences between integral-equation methods are derived. PMID:19275391
Fast Quantum Algorithms for Numerical Integrals and Stochastic Processes
NASA Technical Reports Server (NTRS)
Abrams, D.; Williams, C.
1999-01-01
We discuss quantum algorithms that calculate numerical integrals and descriptive statistics of stochastic processes. With either of two distinct approaches, one obtains an exponential speed increase in comparison to the fastest known classical deterministic algotithms and a quadratic speed increase incomparison to classical Monte Carlo methods.
Integrated product definition representation for agile numerical control applications
Simons, W.R. Jr.; Brooks, S.L.; Kirk, W.J. III; Brown, C.W.
1994-11-01
Realization of agile manufacturing capabilities for a virtual enterprise requires the integration of technology, management, and work force into a coordinated, interdependent system. This paper is focused on technology enabling tools for agile manufacturing within a virtual enterprise specifically relating to Numerical Control (N/C) manufacturing activities and product definition requirements for these activities.
Monograph - The Numerical Integration of Ordinary Differential Equations.
ERIC Educational Resources Information Center
Hull, T. E.
The materials presented in this monograph are intended to be included in a course on ordinary differential equations at the upper division level in a college mathematics program. These materials provide an introduction to the numerical integration of ordinary differential equations, and they can be used to supplement a regular text on this…
Designing Adaptive Low-Dissipative High Order Schemes for Long-Time Integrations. Chapter 1
NASA Technical Reports Server (NTRS)
Yee, Helen C.; Sjoegreen, B.; Mansour, Nagi N. (Technical Monitor)
2001-01-01
A general framework for the design of adaptive low-dissipative high order schemes is presented. It encompasses a rather complete treatment of the numerical approach based on four integrated design criteria: (1) For stability considerations, condition the governing equations before the application of the appropriate numerical scheme whenever it is possible; (2) For consistency, compatible schemes that possess stability properties, including physical and numerical boundary condition treatments, similar to those of the discrete analogue of the continuum are preferred; (3) For the minimization of numerical dissipation contamination, efficient and adaptive numerical dissipation control to further improve nonlinear stability and accuracy should be used; and (4) For practical considerations, the numerical approach should be efficient and applicable to general geometries, and an efficient and reliable dynamic grid adaptation should be used if necessary. These design criteria are, in general, very useful to a wide spectrum of flow simulations. However, the demand on the overall numerical approach for nonlinear stability and accuracy is much more stringent for long-time integration of complex multiscale viscous shock/shear/turbulence/acoustics interactions and numerical combustion. Robust classical numerical methods for less complex flow physics are not suitable or practical for such applications. The present approach is designed expressly to address such flow problems, especially unsteady flows. The minimization of employing very fine grids to overcome the production of spurious numerical solutions and/or instability due to under-resolved grids is also sought. The incremental studies to illustrate the performance of the approach are summarized. Extensive testing and full implementation of the approach is forthcoming. The results shown so far are very encouraging.
Technical Report: Scalable Parallel Algorithms for High Dimensional Numerical Integration
Masalma, Yahya; Jiao, Yu
2010-10-01
We implemented a scalable parallel quasi-Monte Carlo numerical high-dimensional integration for tera-scale data points. The implemented algorithm uses the Sobol s quasi-sequences to generate random samples. Sobol s sequence was used to avoid clustering effects in the generated random samples and to produce low-discrepancy random samples which cover the entire integration domain. The performance of the algorithm was tested. Obtained results prove the scalability and accuracy of the implemented algorithms. The implemented algorithm could be used in different applications where a huge data volume is generated and numerical integration is required. We suggest using the hyprid MPI and OpenMP programming model to improve the performance of the algorithms. If the mixed model is used, attention should be paid to the scalability and accuracy.
Numerical integration of ordinary differential equations on manifolds
NASA Astrophysics Data System (ADS)
Crouch, P. E.; Grossman, R.
1993-12-01
This paper is concerned with the problem of developing numerical integration algorithms for differential equations that, when viewed as equations in some Euclidean space, naturally evolve on some embedded submanifold. It is desired to construct algorithms whose iterates also evolve on the same manifold. These algorithms can therefore be viewed as integrating ordinary differential equations on manifolds. The basic method “decouples” the computation of flows on the submanifold from the numerical integration process. It is shown that two classes of single-step and multistep algorithms can be posed and analyzed theoretically, using the concept of “freezing” the coefficients of differential operators obtained from the defining vector field. Explicit third-order algorithms are derived, with additional equations augmenting those of their classical counterparts, obtained from “obstructions” defined by nonvanishing Lie brackets.
Complexity and network dynamics in physiological adaptation: an integrated view.
Baffy, György; Loscalzo, Joseph
2014-05-28
Living organisms constantly interact with their surroundings and sustain internal stability against perturbations. This dynamic process follows three fundamental strategies (restore, explore, and abandon) articulated in historical concepts of physiological adaptation such as homeostasis, allostasis, and the general adaptation syndrome. These strategies correspond to elementary forms of behavior (ordered, chaotic, and static) in complex adaptive systems and invite a network-based analysis of the operational characteristics, allowing us to propose an integrated framework of physiological adaptation from a complex network perspective. Applicability of this concept is illustrated by analyzing molecular and cellular mechanisms of adaptation in response to the pervasive challenge of obesity, a chronic condition resulting from sustained nutrient excess that prompts chaotic exploration for system stability associated with tradeoffs and a risk of adverse outcomes such as diabetes, cardiovascular disease, and cancer. Deconstruction of this complexity holds the promise of gaining novel insights into physiological adaptation in health and disease. PMID:24751342
Ensemble-type numerical uncertainty information from single model integrations
Rauser, Florian Marotzke, Jochem; Korn, Peter
2015-07-01
We suggest an algorithm that quantifies the discretization error of time-dependent physical quantities of interest (goals) for numerical models of geophysical fluid dynamics. The goal discretization error is estimated using a sum of weighted local discretization errors. The key feature of our algorithm is that these local discretization errors are interpreted as realizations of a random process. The random process is determined by the model and the flow state. From a class of local error random processes we select a suitable specific random process by integrating the model over a short time interval at different resolutions. The weights of the influences of the local discretization errors on the goal are modeled as goal sensitivities, which are calculated via automatic differentiation. The integration of the weighted realizations of local error random processes yields a posterior ensemble of goal approximations from a single run of the numerical model. From the posterior ensemble we derive the uncertainty information of the goal discretization error. This algorithm bypasses the requirement of detailed knowledge about the models discretization to generate numerical error estimates. The algorithm is evaluated for the spherical shallow-water equations. For two standard test cases we successfully estimate the error of regional potential energy, track its evolution, and compare it to standard ensemble techniques. The posterior ensemble shares linear-error-growth properties with ensembles of multiple model integrations when comparably perturbed. The posterior ensemble numerical error estimates are of comparable size as those of a stochastic physics ensemble.
Stability of numerical integration techniques for transient rotor dynamics
NASA Technical Reports Server (NTRS)
Kascak, A. F.
1977-01-01
A finite element model of a rotor bearing system was analyzed to determine the stability limits of the forward, backward, and centered Euler; Runge-Kutta; Milne; and Adams numerical integration techniques. The analysis concludes that the highest frequency mode determines the maximum time step for a stable solution. Thus, the number of mass elements should be minimized. Increasing the damping can sometimes cause numerical instability. For a uniform shaft, with 10 mass elements, operating at approximately the first critical speed, the maximum time step for the Runge-Kutta, Milne, and Adams methods is that which corresponds to approximately 1 degree of shaft movement. This is independent of rotor dimensions.
Numerical computation of 2D Sommerfeld integrals - Decomposition of the angular integral
NASA Astrophysics Data System (ADS)
Dvorak, Steven L.; Kuester, Edward F.
1992-02-01
The computational efficiency of the 2D Sommerfeld integrals is shown to undergo improvement through the discovery of novel ways to compute the inner angular integral in polar representations. It is shown that the angular integral can be decomposed into a finite number of incomplete Lipschitz-Hankel integrals; these can in turn be calculated through a series of expansions, so that the angular integral can be computed by summing a series rather than applying a standard numerical integration algorithm. The technique is most efficient and accurate when piecewise-sinusoidal basis functions are employed to analyze a printed strip-dipole antenna in a layered medium.
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.
Adaptive Spontaneous Transitions between Two Mechanisms of Numerical Averaging.
Brezis, Noam; Bronfman, Zohar Z; Usher, Marius
2015-01-01
We investigated the mechanism with which humans estimate numerical averages. Participants were presented with 4, 8 or 16 (two-digit) numbers, serially and rapidly (2 numerals/second) and were instructed to convey the sequence average. As predicted by a dual, but not a single-component account, we found a non-monotonic influence of set-size on accuracy. Moreover, we observed a marked decrease in RT as set-size increases and RT-accuracy tradeoff in the 4-, but not in the 16-number condition. These results indicate that in accordance with the normative directive, participants spontaneously employ analytic/sequential thinking in the 4-number condition and intuitive/holistic thinking in the 16-number condition. When the presentation rate is extreme (10 items/sec) we find that, while performance still remains high, the estimations are now based on intuitive processing. The results are accounted for by a computational model postulating population-coding underlying intuitive-averaging and working-memory-mediated symbolic procedures underlying analytical-averaging, with flexible allocation between the two. PMID:26041580
Adaptive Spontaneous Transitions between Two Mechanisms of Numerical Averaging
Brezis, Noam; Bronfman, Zohar Z.; Usher, Marius
2015-01-01
We investigated the mechanism with which humans estimate numerical averages. Participants were presented with 4, 8 or 16 (two-digit) numbers, serially and rapidly (2 numerals/second) and were instructed to convey the sequence average. As predicted by a dual, but not a single-component account, we found a non-monotonic influence of set-size on accuracy. Moreover, we observed a marked decrease in RT as set-size increases and RT-accuracy tradeoff in the 4-, but not in the 16-number condition. These results indicate that in accordance with the normative directive, participants spontaneously employ analytic/sequential thinking in the 4-number condition and intuitive/holistic thinking in the 16-number condition. When the presentation rate is extreme (10 items/sec) we find that, while performance still remains high, the estimations are now based on intuitive processing. The results are accounted for by a computational model postulating population-coding underlying intuitive-averaging and working-memory-mediated symbolic procedures underlying analytical-averaging, with flexible allocation between the two. PMID:26041580
Simple and Efficient Numerical Evaluation of Near-Hypersingular Integrals
NASA Technical Reports Server (NTRS)
Fink, Patrick W.; Wilton, Donald R.; Khayat, Michael A.
2007-01-01
Recently, significant progress has been made in the handling of singular and nearly-singular potential integrals that commonly arise in the Boundary Element Method (BEM). To facilitate object-oriented programming and handling of higher order basis functions, cancellation techniques are favored over techniques involving singularity subtraction. However, gradients of the Newton-type potentials, which produce hypersingular kernels, are also frequently required in BEM formulations. As is the case with the potentials, treatment of the near-hypersingular integrals has proven more challenging than treating the limiting case in which the observation point approaches the surface. Historically, numerical evaluation of these near-hypersingularities has often involved a two-step procedure: a singularity subtraction to reduce the order of the singularity, followed by a boundary contour integral evaluation of the extracted part. Since this evaluation necessarily links basis function, Green s function, and the integration domain (element shape), the approach ill fits object-oriented programming concepts. Thus, there is a need for cancellation-type techniques for efficient numerical evaluation of the gradient of the potential. Progress in the development of efficient cancellation-type procedures for the gradient potentials was recently presented. To the extent possible, a change of variables is chosen such that the Jacobian of the transformation cancels the singularity. However, since the gradient kernel involves singularities of different orders, we also require that the transformation leaves remaining terms that are analytic. The terms "normal" and "tangential" are used herein with reference to the source element. Also, since computational formulations often involve the numerical evaluation of both potentials and their gradients, it is highly desirable that a single integration procedure efficiently handles both.
System integration of pattern recognition, adaptive aided, upper limb prostheses
NASA Technical Reports Server (NTRS)
Lyman, J.; Freedy, A.; Solomonow, M.
1975-01-01
The requirements for successful integration of a computer aided control system for multi degree of freedom artificial arms are discussed. Specifications are established for a system which shares control between a human amputee and an automatic control subsystem. The approach integrates the following subsystems: (1) myoelectric pattern recognition, (2) adaptive computer aiding; (3) local reflex control; (4) prosthetic sensory feedback; and (5) externally energized arm with the functions of prehension, wrist rotation, elbow extension and flexion and humeral rotation.
Integrating Learning Styles into Adaptive E-Learning System
ERIC Educational Resources Information Center
Truong, Huong May
2015-01-01
This paper provides an overview and update on my PhD research project which focuses on integrating learning styles into adaptive e-learning system. The project, firstly, aims to develop a system to classify students' learning styles through their online learning behaviour. This will be followed by a study on the complex relationship between…
Integrating Numerical Groundwater Modeling Results With Geographic Information Systems
NASA Astrophysics Data System (ADS)
Witkowski, M. S.; Robinson, B. A.; Linger, S. P.
2001-12-01
Many different types of data are used to create numerical models of flow and transport of groundwater in the vadose zone. Results from water balance studies, infiltration models, hydrologic properties, and digital elevation models (DEMs) are examples of such data. Because input data comes in a variety of formats, for consistency the data need to be assembled in a coherent fashion on a single platform. Through the use of a geographic information system (GIS), all data sources can effectively be integrated on one platform to store, retrieve, query, and display data. In our vadoze zone modeling studies in support of Los Alamos National Laboratory's Environmental Restoration Project, we employ a GIS comprised of a Raid storage device, an Oracle database, ESRI's spatial database engine (SDE), ArcView GIS, and custom GIS tools for three-dimensional (3D) analysis. We store traditional GIS data, such as, contours, historical building footprints, and study area locations, as points, lines, and polygons with attributes. Numerical flow and transport model results from the Finite Element Heat and Mass Transfer Code (FEHM) are stored as points with attributes, such as fluid saturation, or pressure, or contaminant concentration at a given location. We overlay traditional types of GIS data with numerical model results, thereby allowing us to better build conceptual models and perform spatial analyses. We have also developed specialized analysis tools to assist in the data and model analysis process. This approach provides an integrated framework for performing tasks such as comparing the model to data and understanding the relationship of model predictions to existing contaminant source locations and water supply wells. Our process of integrating GIS and numerical modeling results allows us to answer a wide variety of questions about our conceptual model design: - Which set of locations should be identified as contaminant sources based on known historical building operations
Comparison of integrated numerical experiments with accelerator and FEL experiments
Thode, L.E.; Carlsten, B.E.; Chan, K.C.D.; Cooper, R.K.; Elliott, J.C.; Gitomer, S.J.; Goldstein, J.C.; Jones, M.E.; McVey, B.D.; Schmitt, M.J.; Takeda, H.; Tokar, R.L.; Wang, T.S.; Young, L.M.
1991-01-01
Even at the conceptual level the strong coupling between the laser subsystem elements, such as the accelerator, wiggler, optics, and control, greatly complicates the understanding and design of an FEL. Given the requirements for a high-performance FEL, the coupling between the laser subsystems must be included in the design approach. To address the subsystem coupling the concept of an integrated numerical experiment (INEX) has been implemented. Unique features of the INEX approach are consistency and numerical equivalence of experimental diagnostic. The equivalent numerical diagnostics mitigates the major problem of misinterpretation that often occurs when theoretical and experimental data are compared. A complete INEX model has been applied to the 10{mu}m high-extraction-efficiency experiment at Los Alamos and the 0.6-{mu}m Burst Mode experiment at Boeing Aerospace. In addition, various subsets of the INEX model have been compared with a number of other experiments. Overall, the agreement between INEX and the experiments is very good. With the INEX approach, it now appears possible to design high-performance FELS for numerous applications. The first full-scale test of the INEX approach is the Los Alamos HIBAF experiment. The INEX concept, implementation, and validation with experiments are discussed. 28 refs., 13 figs., 1 tab.
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.
INEX (integrated numerical experiment) simulations of the Boeing FEL system
Tokar, R.L.; Young, L.M.; Lumpkin, A.H.; McVey, B.D.; Thode, L.E.; Bender, S.C.; Chan, K.C.D. ); Yeremian, A.D.; Dowell, D.H.; Lowrey, A.R. )
1989-01-01
The INEX (integrated numerical experiment) numerical model is applied to the 0.6 {mu}m FEL oscillator at Boeing Aerospace and Electronics Company in Seattle, WA. This system consists of a 110 MeV L-band rf linac, a beam transport line from the accelerator to the entrance of the wiggler, the 5.0 meter THUNDER variable taper wiggler, and a near concentric two mirror optical oscillator. Many aspects of the model for the electron beam accelerator and transport line agree with experimental measurements. Predictions for lasing performance are compared with data obtained in May and June 1989 using a mild tapered wiggler. We obtain good agreement with the achieved extraction efficiency, while 1D pulse simulations reproduce the observed sideband instability. 15 refs., 11 figs.
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.
A wavelet-optimized, very high order adaptive grid and order numerical method
NASA Technical Reports Server (NTRS)
Jameson, Leland
1996-01-01
Differencing operators of arbitrarily high order can be constructed by interpolating a polynomial through a set of data followed by differentiation of this polynomial and finally evaluation of the polynomial at the point where a derivative approximation is desired. Furthermore, the interpolating polynomial can be constructed from algebraic, trigonometric, or, perhaps exponential polynomials. This paper begins with a comparison of such differencing operator construction. Next, the issue of proper grids for high order polynomials is addressed. Finally, an adaptive numerical method is introduced which adapts the numerical grid and the order of the differencing operator depending on the data. The numerical grid adaptation is performed on a Chebyshev grid. That is, at each level of refinement the grid is a Chebvshev grid and this grid is refined locally based on wavelet analysis.
McDowell, Julia Z.; Luber, George
2011-01-01
Background: Climate change is expected to have a range of health impacts, some of which are already apparent. Public health adaptation is imperative, but there has been little discussion of how to increase adaptive capacity and resilience in public health systems. Objectives: We explored possible explanations for the lack of work on adaptive capacity, outline climate–health challenges that may lie outside public health’s coping range, and consider changes in practice that could increase public health’s adaptive capacity. Methods: We conducted a substantive, interdisciplinary literature review focused on climate change adaptation in public health, social learning, and management of socioeconomic systems exhibiting dynamic complexity. Discussion: There are two competing views of how public health should engage climate change adaptation. Perspectives differ on whether climate change will primarily amplify existing hazards, requiring enhancement of existing public health functions, or present categorically distinct threats requiring innovative management strategies. In some contexts, distinctly climate-sensitive health threats may overwhelm public health’s adaptive capacity. Addressing these threats will require increased emphasis on institutional learning, innovative management strategies, and new and improved tools. Adaptive management, an iterative framework that embraces uncertainty, uses modeling, and integrates learning, may be a useful approach. We illustrate its application to extreme heat in an urban setting. Conclusions: Increasing public health capacity will be necessary for certain climate–health threats. Focusing efforts to increase adaptive capacity in specific areas, promoting institutional learning, embracing adaptive management, and developing tools to facilitate these processes are important priorities and can improve the resilience of local public health systems to climate change. PMID:21997387
NASA Astrophysics Data System (ADS)
Rastigejev, Y.; Semakin, A. N.
2012-12-01
In this work we present a multilevel Wavelet-based Adaptive Mesh Refinement (WAMR) method for numerical modeling of global atmospheric chemical transport problems. An accurate numerical simulation of such problems presents an enormous challenge. Atmospheric Chemical Transport Models (CTMs) combine chemical reactions with meteorologically predicted atmospheric advection and turbulent mixing. The resulting system of multi-scale advection-reaction-diffusion equations is extremely stiff, nonlinear and involves a large number of chemically interacting species. As a consequence, the need for enormous computational resources for solving these equations imposes severe limitations on the spatial resolution of the CTMs implemented on uniform or quasi-uniform grids. In turn, this relatively crude spatial resolution results in significant numerical diffusion introduced into the system. This numerical diffusion is shown to noticeably distort the pollutant mixing and transport dynamics for typically used grid resolutions. The developed WAMR method for numerical modeling of atmospheric chemical evolution equations presented in this work provides a significant reduction in the computational cost, without upsetting numerical accuracy, therefore it addresses the numerical difficulties described above. WAMR method introduces a fine grid in the regions where sharp transitions occur and cruder grid in the regions of smooth solution behavior. Therefore WAMR results in much more accurate solutions than conventional numerical methods implemented on uniform or quasi-uniform grids. The algorithm allows one to provide error estimates of the solution that are used in conjunction with appropriate threshold criteria to adapt the non-uniform grid. The method has been tested for a variety of problems including numerical simulation of traveling pollution plumes. It was shown that pollution plumes in the remote troposphere can propagate as well-defined layered structures for two weeks or more as
Bergeron, Bryan; Cline, Andrew; Shipley, Jaime
2012-01-01
We have developed a distributed, standards-based architecture that enables simulation and simulator designers to leverage adaptive learning systems. Our approach, which incorporates an electronic competency record, open source LMS, and open source microcontroller hardware, is a low-cost, pragmatic option to integrating simulators with traditional courseware. PMID:22356955
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.
Quantum Calisthenics: Gaussians, The Path Integral and Guided Numerical Approximations
Weinstein, Marvin; /SLAC
2009-02-12
It is apparent to anyone who thinks about it that, to a large degree, the basic concepts of Newtonian physics are quite intuitive, but quantum mechanics is not. My purpose in this talk is to introduce you to a new, much more intuitive way to understand how quantum mechanics works. I begin with an incredibly easy way to derive the time evolution of a Gaussian wave-packet for the case free and harmonic motion without any need to know the eigenstates of the Hamiltonian. This discussion is completely analytic and I will later use it to relate the solution for the behavior of the Gaussian packet to the Feynman path-integral and stationary phase approximation. It will be clear that using the information about the evolution of the Gaussian in this way goes far beyond what the stationary phase approximation tells us. Next, I introduce the concept of the bucket brigade approach to dealing with problems that cannot be handled totally analytically. This approach combines the intuition obtained in the initial discussion, as well as the intuition obtained from the path-integral, with simple numerical tools. My goal is to show that, for any specific process, there is a simple Hilbert space interpretation of the stationary phase approximation. I will then argue that, from the point of view of numerical approximations, the trajectory obtained from my generalization of the stationary phase approximation specifies that subspace of the full Hilbert space that is needed to compute the time evolution of the particular state under the full Hamiltonian. The prescription I will give is totally non-perturbative and we will see, by the grace of Maple animations computed for the case of the anharmonic oscillator Hamiltonian, that this approach allows surprisingly accurate computations to be performed with very little work. I think of this approach to the path-integral as defining what I call a guided numerical approximation scheme. After the discussion of the anharmonic oscillator I will
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.
Analysis of adaptive algorithms for an integrated communication network
NASA Technical Reports Server (NTRS)
Reed, Daniel A.; Barr, Matthew; Chong-Kwon, Kim
1985-01-01
Techniques were examined that trade communication bandwidth for decreased transmission delays. When the network is lightly used, these schemes attempt to use additional network resources to decrease communication delays. As the network utilization rises, the schemes degrade gracefully, still providing service but with minimal use of the network. Because the schemes use a combination of circuit and packet switching, they should respond to variations in the types and amounts of network traffic. Also, a combination of circuit and packet switching to support the widely varying traffic demands imposed on an integrated network was investigated. The packet switched component is best suited to bursty traffic where some delays in delivery are acceptable. The circuit switched component is reserved for traffic that must meet real time constraints. Selected packet routing algorithms that might be used in an integrated network were simulated. An integrated traffic places widely varying workload demands on a network. Adaptive algorithms were identified, ones that respond to both the transient and evolutionary changes that arise in integrated networks. A new algorithm was developed, hybrid weighted routing, that adapts to workload changes.
Progress in integrated analysis with adaptive unstructured meshing
NASA Technical Reports Server (NTRS)
Dechaumphai, Pramote
1992-01-01
Design of lightweight structures and thermal protection systems for hypersonic vehicles depend on accurate prediction of aerothermal loads, structural temperatures and their gradients, and structural deformations and stresses. Concentration is on an alternative meshing technique which generates an entirely new adaptive unstructured mesh based on the solution obtained from the earlier mesh. The technique combined with the finite element method has been shown to significantly improve the efficiency and accuracy of the fluid, thermal, and structural analyses. Current capability of the adaptive unstructured meshing technique for the integrated fluid-thermal-structural analysis is described first. The technique was extended to transient thermal analysis of structures with time-dependent adaptive meshing to capture the detailed temperature response with a minimum number of unknowns and computational cost. Both linear and higher-order finite elements are implemented to demonstrate the generality of the technique and to investigate their solution accuracy. Currently, the adaptive meshing technique is being developed for plane structures that can be modeled with membrane elements and built-up structures modeled with membrane and bending elements. The capability of the technique to these different disciplinary problems is demonstrated by several examples.
Integrated modeling of the GMT laser tomography adaptive optics system
NASA Astrophysics Data System (ADS)
Piatrou, Piotr
2014-08-01
Laser Tomography Adaptive Optics (LTAO) is one of adaptive optics systems planned for the Giant Magellan Telescope (GMT). End-to-end simulation tools that are able to cope with the complexity and computational burden of the AO systems to be installed on the extremely large telescopes such as GMT prove to be an integral part of the GMT LTAO system development endeavors. SL95, the Fortran 95 Simulation Library, is one of the software tools successfully used for the LTAO system end-to-end simulations. The goal of SL95 project is to provide a complete set of generic, richly parameterized mathematical models for key elements of the segmented telescope wavefront control systems including both active and adaptive optics as well as the models for atmospheric turbulence, extended light sources like Laser Guide Stars (LGS), light propagation engines and closed-loop controllers. The library is implemented as a hierarchical collection of classes capable of mutual interaction, which allows one to assemble complex wavefront control system configurations with multiple interacting control channels. In this paper we demonstrate the SL95 capabilities by building an integrated end-to-end model of the GMT LTAO system with 7 control channels: LGS tomography with Adaptive Secondary and on-instrument deformable mirrors, tip-tilt and vibration control, LGS stabilization, LGS focus control, truth sensor-based dynamic noncommon path aberration rejection, pupil position control, SLODAR-like embedded turbulence profiler. The rich parameterization of the SL95 classes allows to build detailed error budgets propagating through the system multiple errors and perturbations such as turbulence-, telescope-, telescope misalignment-, segment phasing error-, non-common path-induced aberrations, sensor noises, deformable mirror-to-sensor mis-registration, vibration, temporal errors, etc. We will present a short description of the SL95 architecture, as well as the sample GMT LTAO system simulation
Shaping the Cities of Tomorrow: Integrating Local Urban Adaptation within an Environmental Framework
NASA Astrophysics Data System (ADS)
Georgescu, M.
2014-12-01
Contemporary methods focused on increasing urban sustainability are largely based on the reduction of greenhouse gas emissions. While these efforts are essential steps forward, continued characterization of urban sustainability solely within a biogeochemical framework, with neglect of the biophysical impact of the built environment, omits regional hydroclimatic forcing of the same order of magnitude as greenhouse gas emissions. Using a suite of continuous, multi-year and multi-member continental scale numerical simulations with the WRF model for the U.S., we examine hydroclimatic impacts for a variety of U.S. urban expansion scenarios (for the year 2100) and urban adaptation futures (cool roofs, green roofs, and a hypothetical hybrid approach integrating biophysical properties of both cool and green roofs), and compare those to experiments utilizing a contemporary urban extent. Widespread adoption of adaptation strategies exhibits regionally and seasonally dependent hydroclimatic impacts. For some regions and seasons, urban-induced warming in excess of 3°C can be completely offset by all adaptation approaches examined. For other regions, widespread adoption of some adaptation approaches leads to significant rainfall decline. Sustainable urban expansion therefore requires an integrated assessment that also incorporates biophysically induced urban impacts, and demands tradeoff assessment of various strategies aimed to ameliorate deleterious consequences of growth (e.g., urban heat island reduction).
Integrated Decision Support for Global Environmental Change Adaptation
NASA Astrophysics Data System (ADS)
Kumar, S.; Cantrell, S.; Higgins, G. J.; Marshall, J.; VanWijngaarden, F.
2011-12-01
Environmental changes are happening now that has caused concern in many parts of the world; particularly vulnerable are the countries and communities with limited resources and with natural environments that are more susceptible to climate change impacts. Global leaders are concerned about the observed phenomena and events such as Amazon deforestation, shifting monsoon patterns affecting agriculture in the mountain slopes of Peru, floods in Pakistan, water shortages in Middle East, droughts impacting water supplies and wildlife migration in Africa, and sea level rise impacts on low lying coastal communities in Bangladesh. These environmental changes are likely to get exacerbated as the temperatures rise, the weather and climate patterns change, and sea level rise continues. Large populations and billions of dollars of infrastructure could be affected. At Northrop Grumman, we have developed an integrated decision support framework for providing necessary information to stakeholders and planners to adapt to the impacts of climate variability and change at the regional and local levels. This integrated approach takes into account assimilation and exploitation of large and disparate weather and climate data sets, regional downscaling (dynamic and statistical), uncertainty quantification and reduction, and a synthesis of scientific data with demographic and economic data to generate actionable information for the stakeholders and decision makers. Utilizing a flexible service oriented architecture and state-of-the-art visualization techniques, this information can be delivered via tailored GIS portals to meet diverse set of user needs and expectations. This integrated approach can be applied to regional and local risk assessments, predictions and decadal projections, and proactive adaptation planning for vulnerable communities. In this paper we will describe this comprehensive decision support approach with selected applications and case studies to illustrate how this
On the Computation of Integral Curves in Adaptive Mesh Refinement Vector Fields
Deines, Eduard; Weber, Gunther H.; Garth, Christoph; Van Straalen, Brian; Borovikov, Sergey; Martin, Daniel F.; Joy, Kenneth I.
2011-06-27
Integral curves, such as streamlines, streaklines, pathlines, and timelines, are an essential tool in the analysis of vector field structures, offering straightforward and intuitive interpretation of visualization results. While such curves have a long-standing tradition in vector field visualization, their application to Adaptive Mesh Refinement (AMR) simulation results poses unique problems. AMR is a highly effective discretization method for a variety of physical simulation problems and has recently been applied to the study of vector fields in flow and magnetohydrodynamic applications. The cell-centered nature of AMR data and discontinuities in the vector field representation arising from AMR level boundaries complicate the application of numerical integration methods to compute integral curves. In this paper, we propose a novel approach to alleviate these problems and show its application to streamline visualization in an AMR model of the magnetic field of the solar system as well as to a simulation of two incompressible viscous vortex rings merging.
Optimizing aircraft performance with adaptive, integrated flight/propulsion control
NASA Technical Reports Server (NTRS)
Smith, R. H.; Chisholm, J. D.; Stewart, J. F.
1991-01-01
The Performance-Seeking Control (PSC) integrated flight/propulsion adaptive control algorithm presented was developed in order to optimize total aircraft performance during steady-state engine operation. The PSC multimode algorithm minimizes fuel consumption at cruise conditions, while maximizing excess thrust during aircraft accelerations, climbs, and dashes, and simultaneously extending engine service life through reduction of fan-driving turbine inlet temperature upon engagement of the extended-life mode. The engine models incorporated by the PSC are continually upgraded, using a Kalman filter to detect anomalous operations. The PSC algorithm will be flight-demonstrated by an F-15 at NASA-Dryden.
Wavelet-based adaptive numerical simulation of unsteady 3D flow around a bluff body
NASA Astrophysics Data System (ADS)
de Stefano, Giuliano; Vasilyev, Oleg
2012-11-01
The unsteady three-dimensional flow past a two-dimensional bluff body is numerically simulated using a wavelet-based method. The body is modeled by exploiting the Brinkman volume-penalization method, which results in modifying the governing equations with the addition of an appropriate forcing term inside the spatial region occupied by the obstacle. The volume-penalized incompressible Navier-Stokes equations are numerically solved by means of the adaptive wavelet collocation method, where the non-uniform spatial grid is dynamically adapted to the flow evolution. The combined approach is successfully applied to the simulation of vortex shedding flow behind a stationary prism with square cross-section. The computation is conducted at transitional Reynolds numbers, where fundamental unstable three-dimensional vortical structures exist, by well-predicting the unsteady forces arising from fluid-structure interaction.
Integrated flight/propulsion control - Adaptive engine control system mode
NASA Technical Reports Server (NTRS)
Yonke, W. A.; Terrell, L. A.; Meyers, L. P.
1985-01-01
The adaptive engine control system mode (ADECS) which is developed and tested on an F-15 aircraft with PW1128 engines, using the NASA sponsored highly integrated digital electronic control program, is examined. The operation of the ADECS mode, as well as the basic control logic, the avionic architecture, and the airframe/engine interface are described. By increasing engine pressure ratio (EPR) additional thrust is obtained at intermediate power and above. To modulate the amount of EPR uptrim and to prevent engine stall, information from the flight control system is used. The performance benefits, anticipated from control integration are shown for a range of flight conditions and power settings. It is found that at higher altitudes, the ADECS mode can increase thrust as much as 12 percent, which is used for improved acceleration, improved turn rate, or sustained turn angle.
Earthquake Rupture Dynamics using Adaptive Mesh Refinement and High-Order Accurate Numerical Methods
NASA Astrophysics Data System (ADS)
Kozdon, J. E.; Wilcox, L.
2013-12-01
Our goal is to develop scalable and adaptive (spatial and temporal) numerical methods for coupled, multiphysics problems using high-order accurate numerical methods. To do so, we are developing an opensource, parallel library known as bfam (available at http://bfam.in). The first application to be developed on top of bfam is an earthquake rupture dynamics solver using high-order discontinuous Galerkin methods and summation-by-parts finite difference methods. In earthquake rupture dynamics, wave propagation in the Earth's crust is coupled to frictional sliding on fault interfaces. This coupling is two-way, required the simultaneous simulation of both processes. The use of laboratory-measured friction parameters requires near-fault resolution that is 4-5 orders of magnitude higher than that needed to resolve the frequencies of interest in the volume. This, along with earlier simulations using a low-order, finite volume based adaptive mesh refinement framework, suggest that adaptive mesh refinement is ideally suited for this problem. The use of high-order methods is motivated by the high level of resolution required off the fault in earlier the low-order finite volume simulations; we believe this need for resolution is a result of the excessive numerical dissipation of low-order methods. In bfam spatial adaptivity is handled using the p4est library and temporal adaptivity will be accomplished through local time stepping. In this presentation we will present the guiding principles behind the library as well as verification of code against the Southern California Earthquake Center dynamic rupture code validation test problems.
Integrated Framework for an Urban Climate Adaptation Tool
NASA Astrophysics Data System (ADS)
Omitaomu, O.; Parish, E. S.; Nugent, P.; Mei, R.; Sylvester, L.; Ernst, K.; Absar, M.
2015-12-01
Cities have an opportunity to become more resilient to future climate change through investments made in urban infrastructure today. However, most cities lack access to credible high-resolution climate change projection information needed to assess and address potential vulnerabilities from future climate variability. Therefore, we present an integrated framework for developing an urban climate adaptation tool (Urban-CAT). Urban-CAT consists of four modules. Firstly, it provides climate projections at different spatial resolutions for quantifying urban landscape. Secondly, this projected data is combined with socio-economic data using leading and lagging indicators for assessing landscape vulnerability to climate extremes (e.g., urban flooding). Thirdly, a neighborhood scale modeling approach is presented for identifying candidate areas for adaptation strategies (e.g., green infrastructure as an adaptation strategy for urban flooding). Finally, all these capabilities are made available as a web-based tool to support decision-making and communication at the neighborhood and city levels. In this paper, we present some of the methods that drive each of the modules and demo some of the capabilities available to-date using the City of Knoxville in Tennessee as a case study.
The adaptive significance of adult neurogenesis: an integrative approach
Konefal, Sarah; Elliot, Mick; Crespi, Bernard
2013-01-01
Adult neurogenesis in mammals is predominantly restricted to two brain regions, the dentate gyrus (DG) of the hippocampus and the olfactory bulb (OB), suggesting that these two brain regions uniquely share functions that mediate its adaptive significance. Benefits of adult neurogenesis across these two regions appear to converge on increased neuronal and structural plasticity that subserves coding of novel, complex, and fine-grained information, usually with contextual components that include spatial positioning. By contrast, costs of adult neurogenesis appear to center on potential for dysregulation resulting in higher risk of brain cancer or psychological dysfunctions, but such costs have yet to be quantified directly. The three main hypotheses for the proximate functions and adaptive significance of adult neurogenesis, pattern separation, memory consolidation, and olfactory spatial, are not mutually exclusive and can be reconciled into a simple general model amenable to targeted experimental and comparative tests. Comparative analysis of brain region sizes across two major social-ecological groups of primates, gregarious (mainly diurnal haplorhines, visually-oriented, and in large social groups) and solitary (mainly noctural, territorial, and highly reliant on olfaction, as in most rodents) suggest that solitary species, but not gregarious species, show positive associations of population densities and home range sizes with sizes of both the hippocampus and OB, implicating their functions in social-territorial systems mediated by olfactory cues. Integrated analyses of the adaptive significance of adult neurogenesis will benefit from experimental studies motivated and structured by ecologically and socially relevant selective contexts. PMID:23882188
Are integral controllers adapted to the new era of ELT adaptive optics?
NASA Astrophysics Data System (ADS)
Conan, J.-M.; Raynaud, H.-F.; Kulcsár, C.; Meimon, S.
2011-09-01
With ELTs we are now entering a new era in adaptive optics developments. Meeting unprecedented level of performance with incredibly complex systems implies reconsidering AO concepts at all levels, including controller design. Concentrating mainly on temporal aspects, one may wonder if integral controllers remain an adequate solution. This question is all the more important that, with ever larger degrees of freedom, one may be tempted to discard more sophisticated approaches because they are deemed too complex to implement. The respective performance of integrator versus LQG control should therefore be carefully evaluated in the ELT context. We recall for instance the impressive correction improvement brought by such controllers for the rejection of windshake and vibration components. LQG controller significantly outperforms the integrator because its disturbance rejection transfer function closely matches the energy concentration, respectively at low temporal frequencies for windshake, and around localized resonant peaks for vibrations. The application to turbulent modes should also be investigated, especially for very low spatial frequencies now explored on the huge ELT pupil. The questions addressed here are: 1/ How do integral and LQG controllers compare in terms of performance for a given sampling frequency and noise level?; 2/ Could we relax sampling frequency with LQG control?; 3/ Does a mode to mode adaptation of temporal rejection bring significant performance improvement?; 4/ Which modes particularly benefit from this fine tuning of the rejection transfer function? Based on a simplified ELT AO configuration, and through a simple analytical formulation, performance is evaluated for several control approaches. Various assumptions concerning the perturbation parameters (seeing and outer-scale value, windshake amplitude) are considered. Bode's integral theorem allows intuitive understanding of the results. Practical implementation and computation complexity
NASA Technical Reports Server (NTRS)
Sjoegreen, B.; Yee, H. C.
2001-01-01
The recently developed essentially fourth-order or higher low dissipative shock-capturing scheme of Yee, Sandham and Djomehri (1999) aimed at minimizing nu- merical dissipations for high speed compressible viscous flows containing shocks, shears and turbulence. To detect non smooth behavior and control the amount of numerical dissipation to be added, Yee et al. employed an artificial compression method (ACM) of Harten (1978) but utilize it in an entirely different context than Harten originally intended. The ACM sensor consists of two tuning parameters and is highly physical problem dependent. To minimize the tuning of parameters and physical problem dependence, new sensors with improved detection properties are proposed. The new sensors are derived from utilizing appropriate non-orthogonal wavelet basis functions and they can be used to completely switch to the extra numerical dissipation outside shock layers. The non-dissipative spatial base scheme of arbitrarily high order of accuracy can be maintained without compromising its stability at all parts of the domain where the solution is smooth. Two types of redundant non-orthogonal wavelet basis functions are considered. One is the B-spline wavelet (Mallat & Zhong 1992) used by Gerritsen and Olsson (1996) in an adaptive mesh refinement method, to determine regions where re nement should be done. The other is the modification of the multiresolution method of Harten (1995) by converting it to a new, redundant, non-orthogonal wavelet. The wavelet sensor is then obtained by computing the estimated Lipschitz exponent of a chosen physical quantity (or vector) to be sensed on a chosen wavelet basis function. Both wavelet sensors can be viewed as dual purpose adaptive methods leading to dynamic numerical dissipation control and improved grid adaptation indicators. Consequently, they are useful not only for shock-turbulence computations but also for computational aeroacoustics and numerical combustion. In addition, these
Adaptive Numerical Dissipation Control in High Order Schemes for Multi-D Non-Ideal MHD
NASA Technical Reports Server (NTRS)
Yee, H. C.; Sjoegreen, B.
2005-01-01
The required type and amount of numerical dissipation/filter to accurately resolve all relevant multiscales of complex MHD unsteady high-speed shock/shear/turbulence/combustion problems are not only physical problem dependent, but also vary from one flow region to another. In addition, proper and efficient control of the divergence of the magnetic field (Div(B)) numerical error for high order shock-capturing methods poses extra requirements for the considered type of CPU intensive computations. The goal is to extend our adaptive numerical dissipation control in high order filter schemes and our new divergence-free methods for ideal MHD to non-ideal MHD that include viscosity and resistivity. The key idea consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and leave the rest of the region free from numerical dissipation contamination. These scheme-independent detectors are capable of distinguishing shocks/shears, flame sheets, turbulent fluctuations and spurious high-frequency oscillations. The detection algorithm is based on an artificial compression method (ACM) (for shocks/shears), and redundant multiresolution wavelets (WAV) (for the above types of flow feature). These filters also provide a natural and efficient way for the minimization of Div(B) numerical error.
Numerical Relations and Skill Level Constrain Co-Adaptive Behaviors of Agents in Sports Teams
Silva, Pedro; Travassos, Bruno; Vilar, Luís; Aguiar, Paulo; Davids, Keith; Araújo, Duarte; Garganta, Júlio
2014-01-01
Similar to other complex systems in nature (e.g., a hunting pack, flocks of birds), sports teams have been modeled as social neurobiological systems in which interpersonal coordination tendencies of agents underpin team swarming behaviors. Swarming is seen as the result of agent co-adaptation to ecological constraints of performance environments by collectively perceiving specific possibilities for action (affordances for self and shared affordances). A major principle of invasion team sports assumed to promote effective performance is to outnumber the opposition (creation of numerical overloads) during different performance phases (attack and defense) in spatial regions adjacent to the ball. Such performance principles are assimilated by system agents through manipulation of numerical relations between teams during training in order to create artificially asymmetrical performance contexts to simulate overloaded and underloaded situations. Here we evaluated effects of different numerical relations differentiated by agent skill level, examining emergent inter-individual, intra- and inter-team coordination. Groups of association football players (national – NLP and regional-level – RLP) participated in small-sided and conditioned games in which numerical relations between system agents were manipulated (5v5, 5v4 and 5v3). Typical grouping tendencies in sports teams (major ranges, stretch indices, distances of team centers to goals and distances between the teams' opposing line-forces in specific team sectors) were recorded by plotting positional coordinates of individual agents through continuous GPS tracking. Results showed that creation of numerical asymmetries during training constrained agents' individual dominant regions, the underloaded teams' compactness and each team's relative position on-field, as well as distances between specific team sectors. We also observed how skill level impacted individual and team coordination tendencies. Data revealed
Numerical relations and skill level constrain co-adaptive behaviors of agents in sports teams.
Silva, Pedro; Travassos, Bruno; Vilar, Luís; Aguiar, Paulo; Davids, Keith; Araújo, Duarte; Garganta, Júlio
2014-01-01
Similar to other complex systems in nature (e.g., a hunting pack, flocks of birds), sports teams have been modeled as social neurobiological systems in which interpersonal coordination tendencies of agents underpin team swarming behaviors. Swarming is seen as the result of agent co-adaptation to ecological constraints of performance environments by collectively perceiving specific possibilities for action (affordances for self and shared affordances). A major principle of invasion team sports assumed to promote effective performance is to outnumber the opposition (creation of numerical overloads) during different performance phases (attack and defense) in spatial regions adjacent to the ball. Such performance principles are assimilated by system agents through manipulation of numerical relations between teams during training in order to create artificially asymmetrical performance contexts to simulate overloaded and underloaded situations. Here we evaluated effects of different numerical relations differentiated by agent skill level, examining emergent inter-individual, intra- and inter-team coordination. Groups of association football players (national--NLP and regional-level--RLP) participated in small-sided and conditioned games in which numerical relations between system agents were manipulated (5v5, 5v4 and 5v3). Typical grouping tendencies in sports teams (major ranges, stretch indices, distances of team centers to goals and distances between the teams' opposing line-forces in specific team sectors) were recorded by plotting positional coordinates of individual agents through continuous GPS tracking. Results showed that creation of numerical asymmetries during training constrained agents' individual dominant regions, the underloaded teams' compactness and each team's relative position on-field, as well as distances between specific team sectors. We also observed how skill level impacted individual and team coordination tendencies. Data revealed emergence of
Numerical modelling of tsunami generation by deformable submarine slides using mesh adaptivity
NASA Astrophysics Data System (ADS)
Smith, Rebecca; Parkinson, Samuel; Hill, Jon; Collins, Gareth; Piggott, Matthew
2014-05-01
Tsunamis generated by submarine slides are often under considered in comparison to earthquake generated tsunami, despite several recent examples. Tsunamigenic slides have generated waves that have caused significant damage and loss of life, for example the 1998 Papua New Guinea submarine mass failure resulted in a tsunami that devastated coastal villages and killed over 2,100 people. Numerical simulations of submarine slide generated waves can help us understand the nature of the waves that are generated, and identify the important factors in determining wave characteristics. There have not been many studies of tsunami generation by deformable submarine slides, largely because of the complexities and computational expense involved in modelling these large scale events. At large, real world, scales modelling of tsunami waves by the generation of slides is computationally challenging. Fluidity is an open source finite element code that is ideally suited to tackle this type of problem as it uses unstructured, adaptive meshes, which help to reduce the computational expense without losing accuracy in the results. Adaptive meshes change topology and resolution based on the current simulation state and as such can focus or reduce resolution when and where it is required. The model also allows a number of different numerical approaches to be taken to simulate the same problem within the same numerical framework. In this example we use multi-material approach, with both two materials (slide and water) and three materials (slide, water and air), alongside a density-driven sediment model approach. We will present results of validating Fluidity against benchmarks from experimental and other numerical studies, at different scales, for deformable underwater slides, and consider the utility of mesh adaptivity. We show good agreement to both laboratory results and other numerical models, both with a fixed mesh and a dynamically adaptive mesh, tracking important features of the
Carbon Dioxide Dispersion in the Combustion Integrated Rack Simulated Numerically
NASA Technical Reports Server (NTRS)
Wu, Ming-Shin; Ruff, Gary A.
2004-01-01
When discharged into an International Space Station (ISS) payload rack, a carbon dioxide (CO2) portable fire extinguisher (PFE) must extinguish a fire by decreasing the oxygen in the rack by 50 percent within 60 sec. The length of time needed for this oxygen reduction throughout the rack and the length of time that the CO2 concentration remains high enough to prevent the fire from reigniting is important when determining the effectiveness of the response and postfire procedures. Furthermore, in the absence of gravity, the local flow velocity can make the difference between a fire that spreads rapidly and one that self-extinguishes after ignition. A numerical simulation of the discharge of CO2 from PFE into the Combustion Integrated Rack (CIR) in microgravity was performed to obtain the local velocity and CO2 concentration. The complicated flow field around the PFE nozzle exits was modeled by sources of equivalent mass and momentum flux at a location downstream of the nozzle. The time for the concentration of CO2 to reach a level that would extinguish a fire anywhere in the rack was determined using the Fire Dynamics Simulator (FDS), a computational fluid dynamics code developed by the National Institute of Standards and Technology specifically to evaluate the development of a fire and smoke transport. The simulation shows that CO2, as well as any smoke and combustion gases produced by a fire, would be discharged into the ISS cabin through the resource utility panel at the bottom of the rack. These simulations will be validated by comparing the results with velocity and CO2 concentration measurements obtained during the fire suppression system verification tests conducted on the CIR in March 2003. Once these numerical simulations are validated, portions of the ISS labs and living areas will be modeled to determine the local flow conditions before, during, and after a fire event. These simulations can yield specific information about how long it takes for smoke and
Black shale weathering: An integrated field and numerical modeling study
NASA Astrophysics Data System (ADS)
Bolton, E. W.; Wildman, R. A., Jr.; Berner, R. A.; Eckert, J. O., Jr.; Petsch, S. T.; Mok, U.; Evans, B.
2003-04-01
We present an integrated study of black shale weathering in a near surface environment. Implications of this study contribute to our understanding of organic matter oxidation in uplifted sediments, along with erosion and reburial of ancient unoxidized organic matter, as major controls on atmospheric oxygen levels over geologic time. The field study used to launch the modeling effort is based on core samples from central-eastern Kentucky near Clay City (Late Devonian New Albany/Ohio Shale), where the strata are essentially horizontal. Samples from various depth intervals (up to 12 m depth) were analyzed for texture (SEM images), porosity fraction (0.02 to 0.1), and horizontal and vertical permeability (water and air permeabilities differ due to the fine-grained nature of the sediments, but are on the order of 0.01 to 1. millidarcies, respectively). Chemical analyses were also performed for per cent C, N, S, and basic mineralogy was determined (clays, quartz, pyrite, in addition to organic matter). The samples contained from 2 to 15 per cent ancient (non-modern soil) organic matter. These results were used in the creation of a numerical model for kinetically controlled oxidation of the organic matter within the shale (based on kinetics from Chang and Berner, 1999). The one-dimensional model includes erosion, oxygen diffusion in the partially saturated vadose zone as well as water percolation and solute transport. This study extends the studies of Petsch (2000) and the weathering component of Lasaga and Ohmoto (2002) to include more reactions (e.g., pyrite oxidation to sulfuric acid and weathering of silicates due to low pH) and to resolve the near-surface boundary layer. The model provides a convenient means of exploring the influence of variable rates of erosion, oxygen level, rainfall, as well as physical and chemical characteristics of the shale on organic matter oxidation.
Liu, Fang; Lin, Lin; Vigil-Fowler, Derek; Lischner, Johannes; Kemper, Alexander F.; Sharifzadeh, Sahar; Jornada, Felipe H. da; Deslippe, Jack; Yang, Chao; and others
2015-04-01
We present a numerical integration scheme for evaluating the convolution of a Green's function with a screened Coulomb potential on the real axis in the GW approximation of the self energy. Our scheme takes the zero broadening limit in Green's function first, replaces the numerator of the integrand with a piecewise polynomial approximation, and performs principal value integration on subintervals analytically. We give the error bound of our numerical integration scheme and show by numerical examples that it is more reliable and accurate than the standard quadrature rules such as the composite trapezoidal rule. We also discuss the benefit of using different self energy expressions to perform the numerical convolution at different frequencies.
NASA Astrophysics Data System (ADS)
Zeff, H. B.; Characklis, G. W.; Reed, P. M.; Herman, J. D.
2015-12-01
Water supply policies that integrate portfolios of short-term management decisions with long-term infrastructure development enable utilities to adapt to a range of future scenarios. An effective mix of short-term management actions can augment existing infrastructure, potentially forestalling new development. Likewise, coordinated expansion of infrastructure such as regional interconnections and shared treatment capacity can increase the effectiveness of some management actions like water transfers. Highly adaptable decision pathways that mix long-term infrastructure options and short-term management actions require decision triggers capable of incorporating the impact of these time-evolving decisions on growing water supply needs. Here, we adapt risk-based triggers to sequence a set of potential infrastructure options in combination with utility-specific conservation actions and inter-utility water transfers. Individual infrastructure pathways can be augmented with conservation or water transfers to reduce the cost of meeting utility objectives, but they can also include cooperatively developed, shared infrastructure that expands regional capacity to transfer water. This analysis explores the role of cooperation among four water utilities in the 'Research Triangle' region of North Carolina by formulating three distinct categories of adaptive policy pathways: independent action (utility-specific conservation and supply infrastructure only), weak cooperation (utility-specific conservation and infrastructure development with regional transfers), and strong cooperation (utility specific conservation and jointly developed of regional infrastructure that supports transfers). Results suggest that strong cooperation aids the utilities in meeting their individual objections at substantially lower costs and with fewer irreversible infrastructure options.
Adaptive multi-sensor integration for mine detection
Baker, J.E.
1997-05-01
State-of-the-art in multi-sensor integration (MSI) application involves extensive research and development time to understand and characterize the application domain; to determine and define the appropriate sensor suite; to analyze, characterize, and calibrate the individual sensor systems; to recognize and accommodate the various sensor interactions; and to develop and optimize robust merging code. Much of this process can benefit from adaptive learning, i.e., an output-based system can take raw sensor data and desired merged results as input and adaptively develop/determine an effective method if interpretation and merger. This approach significantly reduces the time required to apply MSI to a given application, while increasing the quality of the final result and provides a quantitative measure for comparing competing MSI techniques and sensor suites. The ability to automatically develop and optimize MSI techniques for new sensor suites and operating environments makes this approach well suited to the detection of mines and mine-like targets. Perhaps more than any other, this application domain is characterized by diverse, innovative, and dynamic sensor suites, whose nature and interactions are not yet well established. This paper presents such an outcome-based multi-image analysis system. An empirical evaluation of its performance and its application, sensor and domain robustness is presented.
Integrating Numerical Computation into the Modeling Instruction Curriculum
ERIC Educational Resources Information Center
Caballero, Marcos D.; Burk, John B.; Aiken, John M.; Thoms, Brian D.; Douglas, Scott S.; Scanlon, Erin M.; Schatz, Michael F.
2014-01-01
Numerical computation (the use of a computer to solve, simulate, or visualize a physical problem) has fundamentally changed the way scientific research is done. Systems that are too difficult to solve in closed form are probed using computation. Experiments that are impossible to perform in the laboratory are studied numerically. Consequently, in…
Applying integrals of motion to the numerical solution of differential equations
NASA Technical Reports Server (NTRS)
Jezewski, D. J.
1979-01-01
A method is developed for using the integrals of systems of nonlinear, ordinary differential equations in a numerical integration process to control the local errors in these integrals and reduce the global errors of the solution. The method is general and can be applied to either scaler or vector integrals. A number of example problems, with accompanying numerical results, are used to verify the analysis and support the conjecture of global error reduction.
Numerical Integration with GeoGebra in High School
ERIC Educational Resources Information Center
Herceg, Dorde; Herceg, Dragoslav
2010-01-01
The concept of definite integral is almost always introduced as the Riemann integral, which is defined in terms of the Riemann sum, and its geometric interpretation. This definition is hard to understand for high school students. With the aid of mathematical software for visualisation and computation of approximate integrals, the notion of…
Integrated Power Adapter: Isolated Converter with Integrated Passives and Low Material Stress
2010-09-01
ADEPT Project: CPES at Virginia Tech is developing an extremely efficient power converter that could be used in power adapters for small, lightweight laptops and other types of mobile electronic devices. Power adapters convert electrical energy into useable power for an electronic device, and they currently waste a lot of energy when they are plugged into an outlet to power up. CPES at Virginia Tech is integrating high-density capacitors, new magnetic materials, high-frequency integrated circuits, and a constant-flux transformer to create its efficient power converter. The high-density capacitors enable the power adapter to store more energy. The new magnetic materials also increase energy storage, and they can be precisely dispensed using a low-cost ink-jet printer which keeps costs down. The high-frequency integrated circuits can handle more power, and they can handle it more efficiently. And, the constant-flux transformer processes a consistent flow of electrical current, which makes the converter more efficient.
NASA Astrophysics Data System (ADS)
Burago, N. G.; Nikitin, I. S.; Yakushev, V. L.
2016-06-01
Techniques that improve the accuracy of numerical solutions and reduce their computational costs are discussed as applied to continuum mechanics problems with complex time-varying geometry. The approach combines shock-capturing computations with the following methods: (1) overlapping meshes for specifying complex geometry; (2) elastic arbitrarily moving adaptive meshes for minimizing the approximation errors near shock waves, boundary layers, contact discontinuities, and moving boundaries; (3) matrix-free implementation of efficient iterative and explicit-implicit finite element schemes; (4) balancing viscosity (version of the stabilized Petrov-Galerkin method); (5) exponential adjustment of physical viscosity coefficients; and (6) stepwise correction of solutions for providing their monotonicity and conservativeness.
A Diffusion Approximation and Numerical Methods for Adaptive Neuron Models with Stochastic Inputs
Rosenbaum, Robert
2016-01-01
Characterizing the spiking statistics of neurons receiving noisy synaptic input is a central problem in computational neuroscience. Monte Carlo approaches to this problem are computationally expensive and often fail to provide mechanistic insight. Thus, the field has seen the development of mathematical and numerical approaches, often relying on a Fokker-Planck formalism. These approaches force a compromise between biological realism, accuracy and computational efficiency. In this article we develop an extension of existing diffusion approximations to more accurately approximate the response of neurons with adaptation currents and noisy synaptic currents. The implementation refines existing numerical schemes for solving the associated Fokker-Planck equations to improve computationally efficiency and accuracy. Computer code implementing the developed algorithms is made available to the public. PMID:27148036
A Diffusion Approximation and Numerical Methods for Adaptive Neuron Models with Stochastic Inputs.
Rosenbaum, Robert
2016-01-01
Characterizing the spiking statistics of neurons receiving noisy synaptic input is a central problem in computational neuroscience. Monte Carlo approaches to this problem are computationally expensive and often fail to provide mechanistic insight. Thus, the field has seen the development of mathematical and numerical approaches, often relying on a Fokker-Planck formalism. These approaches force a compromise between biological realism, accuracy and computational efficiency. In this article we develop an extension of existing diffusion approximations to more accurately approximate the response of neurons with adaptation currents and noisy synaptic currents. The implementation refines existing numerical schemes for solving the associated Fokker-Planck equations to improve computationally efficiency and accuracy. Computer code implementing the developed algorithms is made available to the public. PMID:27148036
Implicit numerical integration for periodic solutions of autonomous nonlinear systems
NASA Technical Reports Server (NTRS)
Thurston, G. A.
1982-01-01
A change of variables that stabilizes numerical computations for periodic solutions of autonomous systems is derived. Computation of the period is decoupled from the rest of the problem for conservative systems of any order and for any second-order system. Numerical results are included for a second-order conservative system under a suddenly applied constant load. Near the critical load for the system, a small increment in load amplitude results in a large increase in amplitude of the response.
Managing Climate Risk. Integrating Adaptation into World Bank Group Operations
Van Aalst, M.
2006-08-15
Climate change is already taking place, and further changes are inevitable. Developing countries, and particularly the poorest people in these countries, are most at risk. The impacts result not only from gradual changes in temperature and sea level but also, in particular, from increased climate variability and extremes, including more intense floods, droughts, and storms. These changes are already having major impacts on the economic performance of developing countries and on the lives and livelihoods of millions of poor people around the world. Climate change thus directly affects the World Bank Group's mission of eradicating poverty. It also puts at risk many projects in a wide range of sectors, including infrastructure, agriculture, human health, water resources, and environment. The risks include physical threats to the investments, potential underperformance, and the possibility that projects will indirectly contribute to rising vulnerability by, for example, triggering investment and settlement in high-risk areas. The way to address these concerns is not to separate climate change adaptation from other priorities but to integrate comprehensive climate risk management into development planning, programs, and projects. While there is a great need to heighten awareness of climate risk in Bank work, a large body of experience on climate risk management is already available, in analytical work, in country dialogues, and in a growing number of investment projects. This operational experience highlights the general ingredients for successful integration of climate risk management into the mainstream development agenda: getting the right sectoral departments and senior policy makers involved; incorporating risk management into economic planning; engaging a wide range of nongovernmental actors (businesses, nongovernmental organizations, communities, and so on); giving attention to regulatory issues; and choosing strategies that will pay off immediately under current
Adaptive Numerical Dissipative Control in High Order Schemes for Multi-D Non-Ideal MHD
NASA Technical Reports Server (NTRS)
Yee, H. C.; Sjoegreen, B.
2004-01-01
The goal is to extend our adaptive numerical dissipation control in high order filter schemes and our new divergence-free methods for ideal MHD to non-ideal MHD that include viscosity and resistivity. The key idea consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and leave the rest of the region free of numerical dissipation contamination. These scheme-independent detectors are capable of distinguishing shocks/shears, flame sheets, turbulent fluctuations and spurious high-frequency oscillations. The detection algorithm is based on an artificial compression method (ACM) (for shocks/shears), and redundant multi-resolution wavelets (WAV) (for the above types of flow feature). These filter approaches also provide a natural and efficient way for the minimization of Div(B) numerical error. The filter scheme consists of spatially sixth order or higher non-dissipative spatial difference operators as the base scheme for the inviscid flux derivatives. If necessary, a small amount of high order linear dissipation is used to remove spurious high frequency oscillations. For example, an eighth-order centered linear dissipation (AD8) might be included in conjunction with a spatially sixth-order base scheme. The inviscid difference operator is applied twice for the viscous flux derivatives. After the completion of a full time step of the base scheme step, the solution is adaptively filtered by the product of a 'flow detector' and the 'nonlinear dissipative portion' of a high-resolution shock-capturing scheme. In addition, the scheme independent wavelet flow detector can be used in conjunction with spatially compact, spectral or spectral element type of base schemes. The ACM and wavelet filter schemes using the dissipative portion of a second-order shock-capturing scheme with sixth-order spatial central base scheme for both the inviscid and viscous MHD flux
Research on the Evolutionary Strategy Based on AIS and Its Application on Numerical Integration
NASA Astrophysics Data System (ADS)
Bei, Li
Based on the features of artificial immune system, a new evolutionary strategy is proposed in order to calculate the numerical integration of functions. This evolutionary strategy includes the mechanisms of swarm searching and constructing the fitness function. Finally, numerical examples are given for verifying the effectiveness of evolutionary strategy. The results show that the performance of evolutionary strategy is satisfactory and more accurate than traditional methods of numerical integration, such as trapezoid formula and Simpson formula.
Three-dimensional numerical simulations of falling films using an adaptive unstructured mesh
NASA Astrophysics Data System (ADS)
Pain, Chris; Xie, Zhihua; Matar, Omar
2015-11-01
Falling liquid films have rich wave dynamics, often occurring in many industrial applications, such as condensers, evaporators and chemical reactors. A number of numerical studies featuring falling liquid films are available in the literature; the majority of them, however, have focused on two-dimensional falling films. Far fewer studies have considered three-dimensional falling films, and those that have only studied the flow in a periodic domain. The objective of this study is to investigate flow dynamics of developing three-dimensional falling films using the Navier-Stokes equations coupled with interface capturing approach over extended domains. An adaptive, unstructured mesh modelling framework is employed here to study this problem, which can modify and adapt three-dimensional meshes to better represent the underlying physics of multiphase problems and reduce computational effort without sacrificing accuracy. Numerical examples of three-dimensional falling films in a long domain are presented and discussed. EPSRC Programme Grant, MEMPHIS, EP/K0039761/1.
A Numerical Study of Mesh Adaptivity in Multiphase Flows with Non-Newtonian Fluids
NASA Astrophysics Data System (ADS)
Percival, James; Pavlidis, Dimitrios; Xie, Zhihua; Alberini, Federico; Simmons, Mark; Pain, Christopher; Matar, Omar
2014-11-01
We present an investigation into the computational efficiency benefits of dynamic mesh adaptivity in the numerical simulation of transient multiphase fluid flow problems involving Non-Newtonian fluids. Such fluids appear in a range of industrial applications, from printing inks to toothpastes and introduce new challenges for mesh adaptivity due to the additional ``memory'' of viscoelastic fluids. Nevertheless, the multiscale nature of these flows implies huge potential benefits for a successful implementation. The study is performed using the open source package Fluidity, which couples an unstructured mesh control volume finite element solver for the multiphase Navier-Stokes equations to a dynamic anisotropic mesh adaptivity algorithm, based on estimated solution interpolation error criteria, and conservative mesh-to-mesh interpolation routine. The code is applied to problems involving rheologies ranging from simple Newtonian to shear-thinning to viscoelastic materials and verified against experimental data for various industrial and microfluidic flows. This work was undertaken as part of the EPSRC MEMPHIS programme grant EP/K003976/1.
Construction of the two-electron contribution to the Fock matrix by numerical integration
NASA Astrophysics Data System (ADS)
Losilla, Sergio A.; Mehine, Mooses M.; Sundholm, Dage
2012-10-01
A novel method to numerically calculate the Fock matrix is presented. The Coulomb operator is re-expressed as an integral identity, which is discretized. The discretization of the auxiliary t dimension separates the x, y, and z dependencies transforming the two-electron Coulomb integrals of Gaussian-type orbitals (GTO) to a linear sum of products of two-dimensional integrals. The s-type integrals are calculated analytically and integrals of the higher angular-momentum functions are obtained using recursion formulae. The contributions to the two-body Coulomb integrals obtained for each discrete t value can be evaluated independently. The two-body Fock matrix elements can be integrated numerically, using common sets of quadrature points and weights. The aim is to calculate Fock matrices of enough accuracy for electronic structure calculations. Preliminary calculations indicate that it is possible to achieve an overall accuracy of at least 10-12 E h using the numerical approach.
NASA Technical Reports Server (NTRS)
Hickey, Michael Philip
1988-01-01
A proposed replacement scheme for the integration of the barometric and diffusion equations in the NASA Marshall Engineering Thermosphere (MET) model is presented. This proposed integration scheme is based on Gaussian Quadrature. Extensive numerical testing reveals it to be faster, more accurate and more reliable than the present integration scheme (a modified form of Simpson's Rule) used in the MET model. Numerous graphical examples are provided, along with a listing of a modified form of the MET model in which subroutine INTEGRATE (using Simpson's Rule) is replaced by subroutine GAUSS (which uses Gaussian Quadrature). It is recommended that the Gaussian Quadrature integration scheme, as used here, be used in the MET model.
NASA Astrophysics Data System (ADS)
Jiang, Shidong; Luo, Li-Shi
2016-07-01
The integral equation for the flow velocity u (x ; k) in the steady Couette flow derived from the linearized Bhatnagar-Gross-Krook-Welander kinetic equation is studied in detail both theoretically and numerically in a wide range of the Knudsen number k between 0.003 and 100.0. First, it is shown that the integral equation is a Fredholm equation of the second kind in which the norm of the compact integral operator is less than 1 on Lp for any 1 ≤ p ≤ ∞ and thus there exists a unique solution to the integral equation via the Neumann series. Second, it is shown that the solution is logarithmically singular at the endpoints. More precisely, if x = 0 is an endpoint, then the solution can be expanded as a double power series of the form ∑n=0∞∑m=0∞cn,mxn(xln x) m about x = 0 on a small interval x ∈ (0 , a) for some a > 0. And third, a high-order adaptive numerical algorithm is designed to compute the solution numerically to high precision. The solutions for the flow velocity u (x ; k), the stress Pxy (k), and the half-channel mass flow rate Q (k) are obtained in a wide range of the Knudsen number 0.003 ≤ k ≤ 100.0; and these solutions are accurate for at least twelve significant digits or better, thus they can be used as benchmark solutions.
NASA Astrophysics Data System (ADS)
Roig, Jaume; Stefanov, Evgueniy; Morancho, Frédéric
2007-07-01
The use of super-junction (SJ) techniques in PIN photodiodes is proposed in this letter for the first time with the objective to assist the optoelectronic integrated circuits (OEICs) implementation in complementary metal oxide semiconductor (CMOS), bipolar CMOS (BiCMOS) and bipolar-CMOS-double diffused MOS (BCD) technologies. Its technological viability is also discussed to make it credible as an alternative to other OEICs approaches. Numerical simulation of realistic SJ-PIN devices, widely used in high power electronics, demonstrates the possibility to integrate high-performance CMOS-based OEICs in epitaxial layers with doping concentrations above 1× 1015 cm-3. The induced lateral depletion at low reverse biased voltage, assisted by the alternated N and P-doped pillars, allows high-speed transient response in SJ-PIN detecting wavelengths between 400 and 800 nm. Moreover, other important parameters as the responsivity and the dark current are not degraded in respect to the conventional PIN (C-PIN) structures.
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
Impact of numerical integration on gas curtain simulations
Rider, W.; Kamm, J.
2000-11-01
In recent years, we have presented a less than glowing experimental comparison of hydrodynamic codes with the gas curtain experiment (e.g., Kamm et al. 1999a). Here, we discuss the manner in which the details of the hydrodynamic integration techniques may conspire to produce poor results. This also includes some progress in improving the results and agreement with experimental results. Because our comparison was conducted on the details of the experimental images (i.e., their detailed structural information), our results do not conflict with previously published results of good agreement with Richtmyer-Meshkov instabilities based on the integral scale of mixing. New experimental and analysis techniques are also discussed.
NASA Astrophysics Data System (ADS)
Blanes, Sergio; Budd, Chris J.
2004-05-01
We present a generalisation of the Levi-Civita and Kustaanheimo-Stiefel regularisation. This allows the use of more general time rescalings. In particular, it is possible to find a regularisation which removes the singularity of the equations and preserves scaling invariance. In addition, these equations can, in certain cases, be integrated with explicit symplectic Runge-Kutta-Nyström methods. The combination of both techniques gives an explicit adaptive symplectic (EASY) integrator. We apply those methods to some perturbations of the Kepler problem and illustrate, by means of some numerical examples, when scaling invariant regularisations are more efficient that the LC/KS regularisation.
Advantages of vertically adaptive coordinates in numerical models of stratified shelf seas
NASA Astrophysics Data System (ADS)
Gräwe, Ulf; Holtermann, Peter; Klingbeil, Knut; Burchard, Hans
2015-08-01
Shelf seas such as the North Sea and the Baltic Sea are characterised by spatially and temporally varying stratification that is highly relevant for their physical dynamics and the evolution of their ecosystems. Stratification may vary from unstably stratified (e.g., due to convective surface cooling) to strongly stratified with density jumps of up to 10 kg/m3 per m (e.g., in overflows into the Baltic Sea). Stratification has a direct impact on vertical turbulent transports (e.g., of nutrients) and influences the entrainment rate of ambient water into dense bottom currents which in turn determine the stratification of and oxygen supply to, e.g., the central Baltic Sea. Moreover, the suppression of the vertical diffusivity at the summer thermocline is one of the limiting factors for the vertical exchange of nutrients in the North Sea. Due to limitations of computational resources and since the locations of such density jumps (either by salinity or temperature) are predicted by the model simulation itself, predefined vertical coordinates cannot always reliably resolve these features. Thus, all shelf sea models with a predefined vertical coordinate distribution are inherently subject to under-resolution of the density structure. To solve this problem, Burchard and Beckers (2004) and Hofmeister et al. (2010) developed the concept of vertically adaptive coordinates for ocean models, where zooming of vertical coordinates at locations of strong stratification (and shear) is imposed. This is achieved by solving a diffusion equation for the position of the coordinates (with the diffusivity being proportional to the stratification or shear frequencies). We will show for a coupled model system of the North Sea and the Baltic Sea (resolution ˜ 1.8 km) how numerical mixing is substantially reduced and model results become significantly more realistic when vertically adaptive coordinates are applied. We additionally demonstrate that vertically adaptive coordinates perform well
A time-accurate adaptive grid method and the numerical simulation of a shock-vortex interaction
NASA Technical Reports Server (NTRS)
Bockelie, Michael J.; Eiseman, Peter R.
1990-01-01
A time accurate, general purpose, adaptive grid method is developed that is suitable for multidimensional steady and unsteady numerical simulations. The grid point movement is performed in a manner that generates smooth grids which resolve the severe solution gradients and the sharp transitions in the solution gradients. The temporal coupling of the adaptive grid and the PDE solver is performed with a grid prediction correction method that is simple to implement and ensures the time accuracy of the grid. Time accurate solutions of the 2-D Euler equations for an unsteady shock vortex interaction demonstrate the ability of the adaptive method to accurately adapt the grid to multiple solution features.
Numerical implications of stabilization by the use of integrals
NASA Technical Reports Server (NTRS)
Beaudet, P. R.
1975-01-01
Liapunov or energy restraint methods for dynamic stabilization in two body motion perturbation problems are considered. Results of computerized orbital stabilization estimates show that the application of energy restraint prevents the occurrence of consistent timing errors in the stepwise integration of equations of motion for a nearly circular orbit.
Daeva, S.G.; Setukha, A.V.
2015-03-10
A numerical method for solving a problem of diffraction of acoustic waves by system of solid and thin objects based on the reduction the problem to a boundary integral equation in which the integral is understood in the sense of finite Hadamard value is proposed. To solve this equation we applied piecewise constant approximations and collocation methods numerical scheme. The difference between the constructed scheme and earlier known is in obtaining approximate analytical expressions to appearing system of linear equations coefficients by separating the main part of the kernel integral operator. The proposed numerical scheme is tested on the solution of the model problem of diffraction of an acoustic wave by inelastic sphere.
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.
Numerical simulation of diffusion MRI signals using an adaptive time-stepping method.
Li, Jing-Rebecca; Calhoun, Donna; Poupon, Cyril; Le Bihan, Denis
2014-01-20
The effect on the MRI signal of water diffusion in biological tissues in the presence of applied magnetic field gradient pulses can be modelled by a multiple compartment Bloch-Torrey partial differential equation. We present a method for the numerical solution of this equation by coupling a standard Cartesian spatial discretization with an adaptive time discretization. The time discretization is done using the explicit Runge-Kutta-Chebyshev method, which is more efficient than the forward Euler time discretization for diffusive-type problems. We use this approach to simulate the diffusion MRI signal from the extra-cylindrical compartment in a tissue model of the brain gray matter consisting of cylindrical and spherical cells and illustrate the effect of cell membrane permeability. PMID:24351275
Numerical simulation of diffusion MRI signals using an adaptive time-stepping method
NASA Astrophysics Data System (ADS)
Li, Jing-Rebecca; Calhoun, Donna; Poupon, Cyril; Le Bihan, Denis
2014-01-01
The effect on the MRI signal of water diffusion in biological tissues in the presence of applied magnetic field gradient pulses can be modelled by a multiple compartment Bloch-Torrey partial differential equation. We present a method for the numerical solution of this equation by coupling a standard Cartesian spatial discretization with an adaptive time discretization. The time discretization is done using the explicit Runge-Kutta-Chebyshev method, which is more efficient than the forward Euler time discretization for diffusive-type problems. We use this approach to simulate the diffusion MRI signal from the extra-cylindrical compartment in a tissue model of the brain gray matter consisting of cylindrical and spherical cells and illustrate the effect of cell membrane permeability.
Integration of AdaptiSPECT, a small-animal adaptive SPECT imaging system
Chaix, Cécile; Kovalsky, Stephen; Kosmider, Matthew; Barrett, Harrison H.; Furenlid, Lars R.
2015-01-01
AdaptiSPECT is a pre-clinical adaptive SPECT imaging system under final development at the Center for Gamma-ray Imaging. The system incorporates multiple adaptive features: an adaptive aperture, 16 detectors mounted on translational stages, and the ability to switch between a non-multiplexed and a multiplexed imaging configuration. In this paper, we review the design of AdaptiSPECT and its adaptive features. We then describe the on-going integration of the imaging system. PMID:26347197
Integration, heterochrony, and adaptation in pedal digits of syndactylous marsupials
2008-01-01
Background Marsupial syndactyly is a curious morphology of the foot found in all species of diprotodontian and peramelemorph marsupials. It is traditionally defined as a condition in which digits II and III of the foot are bound by skin and are reduced. Past treatments of marsupial syndactyly have not considered the implications of this unique morphology for broader issues of digit development and evolution, and the ongoing debate regarding its phylogenetic meaning lacks a broad empirical basis. This study undertakes the first interdisciplinary characterisation of syndactyly, using variance/covariance matrix comparisons of morphometric measurements, locomotor indices, ossification sequences, and re-assessment of the largely anecdotal data on the phylogenetic distribution of tarsal/metatarsal articulations and "incipient syndactyly". Results Syndactylous digits have virtually identical variance/covariance matrices and display heterochronic ossification timing with respect to digits IV/V. However, this does not impact on overall locomotor adaptation patterns in the syndactylous foot as determined by analysis of locomotor predictor ratios. Reports of incipient syndactyly in some marsupial clades could not be confirmed; contrary to previous claims, syndactyly does not appear to impact on tarsal bone arrangement. Conclusion The results suggest that marsupial syndactyly originates from a constraint that is rooted in early digit ontogeny and results in evolution of the syndactylous digits as a highly integrated unit. Although convergent evolution appears likely, syndactyly in Diprotodontia and Peramelemorpha may occur through homologous developmental processes. We argue that the term "syndactyly" is a misnomer because the marsupial condition only superficially resembles its name-giving human soft-tissue syndactyly. PMID:18501017
NASA Astrophysics Data System (ADS)
von Sydow, Lina
2013-10-01
The discontinuous Galerkin method for time integration of the Black-Scholes partial differential equation for option pricing problems is studied and compared with more standard time-integrators. In space an adaptive finite difference discretization is employed. The results show that the dG method are in most cases at least comparable to standard time-integrators and in some cases superior to them. Together with adaptive spatial grids the suggested pricing method shows great qualities.
Parallel implementation of an adaptive and parameter-free N-body integrator
NASA Astrophysics Data System (ADS)
Pruett, C. David; Ingham, William H.; Herman, Ralph D.
2011-05-01
Previously, Pruett et al. (2003) [3] described an N-body integrator of arbitrarily high order M with an asymptotic operation count of O(MN). The algorithm's structure lends itself readily to data parallelization, which we document and demonstrate here in the integration of point-mass systems subject to Newtonian gravitation. High order is shown to benefit parallel efficiency. The resulting N-body integrator is robust, parameter-free, highly accurate, and adaptive in both time-step and order. Moreover, it exhibits linear speedup on distributed parallel processors, provided that each processor is assigned at least a handful of bodies. Program summaryProgram title: PNB.f90 Catalogue identifier: AEIK_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEIK_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC license, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 3052 No. of bytes in distributed program, including test data, etc.: 68 600 Distribution format: tar.gz Programming language: Fortran 90 and OpenMPI Computer: All shared or distributed memory parallel processors Operating system: Unix/Linux Has the code been vectorized or parallelized?: The code has been parallelized but has not been explicitly vectorized. RAM: Dependent upon N Classification: 4.3, 4.12, 6.5 Nature of problem: High accuracy numerical evaluation of trajectories of N point masses each subject to Newtonian gravitation. Solution method: Parallel and adaptive extrapolation in time via power series of arbitrary degree. Running time: 5.1 s for the demo program supplied with the package.
NASA Astrophysics Data System (ADS)
Bargatze, L. F.
2015-12-01
Active Data Archive Product Tracking (ADAPT) is a collection of software routines that permits one to generate XML metadata files to describe and register data products in support of the NASA Heliophysics Virtual Observatory VxO effort. ADAPT is also a philosophy. The ADAPT concept is to use any and all available metadata associated with scientific data to produce XML metadata descriptions in a consistent, uniform, and organized fashion to provide blanket access to the full complement of data stored on a targeted data server. In this poster, we present an application of ADAPT to describe all of the data products that are stored by using the Common Data File (CDF) format served out by the CDAWEB and SPDF data servers hosted at the NASA Goddard Space Flight Center. These data servers are the primary repositories for NASA Heliophysics data. For this purpose, the ADAPT routines have been used to generate data resource descriptions by using an XML schema named Space Physics Archive, Search, and Extract (SPASE). SPASE is the designated standard for documenting Heliophysics data products, as adopted by the Heliophysics Data and Model Consortium. The set of SPASE XML resource descriptions produced by ADAPT includes high-level descriptions of numerical data products, display data products, or catalogs and also includes low-level "Granule" descriptions. A SPASE Granule is effectively a universal access metadata resource; a Granule associates an individual data file (e.g. a CDF file) with a "parent" high-level data resource description, assigns a resource identifier to the file, and lists the corresponding assess URL(s). The CDAWEB and SPDF file systems were queried to provide the input required by the ADAPT software to create an initial set of SPASE metadata resource descriptions. Then, the CDAWEB and SPDF data repositories were queried subsequently on a nightly basis and the CDF file lists were checked for any changes such as the occurrence of new, modified, or deleted
Feischl, Michael; Gantner, Gregor; Praetorius, Dirk
2015-01-01
We consider the Galerkin boundary element method (BEM) for weakly-singular integral equations of the first-kind in 2D. We analyze some residual-type a posteriori error estimator which provides a lower as well as an upper bound for the unknown Galerkin BEM error. The required assumptions are weak and allow for piecewise smooth parametrizations of the boundary, local mesh-refinement, and related standard piecewise polynomials as well as NURBS. In particular, our analysis gives a first contribution to adaptive BEM in the frame of isogeometric analysis (IGABEM), for which we formulate an adaptive algorithm which steers the local mesh-refinement and the multiplicity of the knots. Numerical experiments underline the theoretical findings and show that the proposed adaptive strategy leads to optimal convergence. PMID:26085698
NASA Astrophysics Data System (ADS)
Calvisi, Michael; Manmi, Kawa; Wang, Qianxi
2014-11-01
Ultrasound contrast agents (UCAs) are microbubbles stabilized with a shell typically of lipid, polymer, or protein and are emerging as a unique tool for noninvasive therapies ranging from gene delivery to tumor ablation. The nonspherical dynamics of contrast agents are thought to play an important role in both diagnostic and therapeutic applications, for example, causing the emission of subharmonic frequency components and enhancing the uptake of therapeutic agents across cell membranes and tissue interfaces. A three-dimensional model for nonspherical contrast agent dynamics based on the boundary integral method is presented. The effects of the encapsulating shell are approximated by adapting Hoff's model for thin-shell, spherical contrast agents to the nonspherical case. A high-quality mesh of the bubble surface is maintained by implementing a hybrid approach of the Lagrangian method and elastic mesh technique. Numerical analyses for the dynamics of UCAs in an infinite liquid and near a rigid wall are performed in parameter regimes of clinical relevance. The results show that the presence of a coating significantly reduces the oscillation amplitude and period, increases the ultrasound pressure amplitude required to incite jetting, and reduces the jet width and velocity.
On the accuracy of numerical integration over the unit sphere applied to full network models
NASA Astrophysics Data System (ADS)
Itskov, Mikhail
2016-05-01
This paper is motivated by a recent study by Verron (Mecha Mater 89:216-228, 2015) which revealed huge errors of the numerical integration over the unit sphere in application to large strain problems. For the verification of numerical integration schemes we apply here other analytical integrals over the unit sphere which demonstrate much more accurate results. Relative errors of these integrals with respect to corresponding analytical solutions are evaluated also for a full network model of rubber elasticity based on a Padé approximation of the inverse Langevin function as the chain force. According to the results of our study, the numerical integration over the unit sphere can still be considered as a reliable and accurate tool for full network models.
An adaptive procedure for the numerical parameters of a particle simulation
NASA Astrophysics Data System (ADS)
Galitzine, Cyril; Boyd, Iain D.
2015-01-01
In this article, a computational procedure that automatically determines the optimum time step, cell weight and species weights for steady-state multi-species DSMC (direct simulation Monte Carlo) simulations is presented. The time step is required to satisfy the basic requirements of the DSMC method while the weight and relative weights fields are chosen so as to obtain a user-specified average number of particles in all cells of the domain. The procedure allows the conduct of efficient DSMC simulations with minimal user input and is integrable into existing DSMC codes. The adaptive method is used to simulate a test case consisting of two counterflowing jets at a Knudsen number of 0.015. Large accuracy gains for sampled number densities and velocities over a standard simulation approach for the same number of particles are observed.
NASA Astrophysics Data System (ADS)
Grothmann, T.; Grecksch, K.; Winges, M.; Siebenhüner, B.
2013-12-01
Several case studies show that social factors like institutions, perceptions and social capital strongly affect social capacities to adapt to climate change. Together with economic and technological development they are important for building social capacities. However, there are almost no methodologies for the systematic assessment of social factors. After reviewing existing methodologies we identify the Adaptive Capacity Wheel (ACW) by Gupta et al. (2010), developed for assessing the adaptive capacity of institutions, as the most comprehensive and operationalised framework to assess social factors. The ACW differentiates 22 criteria to assess 6 dimensions: variety, learning capacity, room for autonomous change, leadership, availability of resources, fair governance. To include important psychological factors we extended the ACW by two dimensions: "adaptation motivation" refers to actors' motivation to realise, support and/or promote adaptation to climate; "adaptation belief" refers to actors' perceptions of realisability and effectiveness of adaptation measures. We applied the extended ACW to assess adaptive capacities of four sectors - water management, flood/coastal protection, civil protection and regional planning - in northwestern Germany. The assessments of adaptation motivation and belief provided a clear added value. The results also revealed some methodological problems in applying the ACW (e.g. overlap of dimensions), for which we propose methodological solutions.
NASA Astrophysics Data System (ADS)
Grothmann, T.; Grecksch, K.; Winges, M.; Siebenhüner, B.
2013-03-01
Several case studies show that "soft social factors" (e.g. institutions, perceptions, social capital) strongly affect social capacities to adapt to climate change. Many soft social factors can probably be changed faster than "hard social factors" (e.g. economic and technological development) and are therefore particularly important for building social capacities. However, there are almost no methodologies for the systematic assessment of soft social factors. Gupta et al. (2010) have developed the Adaptive Capacity Wheel (ACW) for assessing the adaptive capacity of institutions. The ACW differentiates 22 criteria to assess six dimensions: variety, learning capacity, room for autonomous change, leadership, availability of resources, fair governance. To include important psychological factors we extended the ACW by two dimensions: "adaptation motivation" refers to actors' motivation to realise, support and/or promote adaptation to climate. "Adaptation belief" refers to actors' perceptions of realisability and effectiveness of adaptation measures. We applied the extended ACW to assess adaptive capacities of four sectors - water management, flood/coastal protection, civil protection and regional planning - in North Western Germany. The assessments of adaptation motivation and belief provided a clear added value. The results also revealed some methodological problems in applying the ACW (e.g. overlap of dimensions), for which we propose methodological solutions.
Zhao, Guoliang; Li, Hongxing
2013-01-01
This paper proposes new methodologies for the design of adaptive integral-sliding mode control. A tensor product model transformation based adaptive integral-sliding mode control law with respect to uncertainties and perturbations is studied, while upper bounds on the perturbations and uncertainties are assumed to be unknown. The advantage of proposed controllers consists in having a dynamical adaptive control gain to establish a sliding mode right at the beginning of the process. Gain dynamics ensure a reasonable adaptive gain with respect to the uncertainties. Finally, efficacy of the proposed controller is verified by simulations on an uncertain nonlinear system model. PMID:24453897
NASA Astrophysics Data System (ADS)
Moniem, T. A.
2016-05-01
This article presents a methodology for an integrated Bragg grating using an alloy of GaAs, AlGaAs, and InGaAs with a controllable refractive index to obtain an adaptive Bragg grating suitable for many applications on optical processing and adaptive control systems, such as limitation and filtering. The refractive index of a Bragg grating is controlled by using an external electric field for controlling periodic modulation of the refractive index of the active waveguide region. The designed Bragg grating has refractive indices programmed by using that external electric field. This article presents two approaches for designing the controllable refractive indices active region of a Bragg grating. The first approach is based on the modification of a planar micro-strip structure of the iGaAs traveling wave as the active region, and the second is based on the modification of self-assembled InAs/GaAs quantum dots of an alloy from GaAs and InGaAs with a GaP traveling wave. The overall design and results are discussed through numerical simulation by using the finite-difference time-domain, plane wave expansion, and opto-wave simulation methods to confirm its operation and feasibility.
NASA Astrophysics Data System (ADS)
Plakhov, Iu. V.; Mytsenko, A. V.; Shel'Pov, V. A.
A numerical integration method is developed that is more accurate than Everhart's (1974) implicit single-sequence approach for integrating orbits. This method can be used to solve problems of space geodesy based on the use of highly precise laser observations.
An efficient step-size control method in numerical integration for astrodynamical equations
NASA Astrophysics Data System (ADS)
Liu, C. Z.; Cui, D. X.
2002-11-01
Using the curvature of the integral curve, a step-size control method is introduced in this paper. This method will prove to be the efficient scheme in the sense that it saves computation time and improve accuracy of numerical integration.
Internal Models in Sensorimotor Integration: Perspectives from Adaptive Control Theory
Tin, Chung; Poon, Chi-Sang
2007-01-01
Internal model and adaptive control are empirical and mathematical paradigms that have evolved separately to describe learning control processes in brain systems and engineering systems, respectively. This paper presents a comprehensive appraisal of the correlation between these paradigms with a view to forging a unified theoretical framework that may benefit both disciplines. It is suggested that the classic equilibrium-point theory of impedance control of arm movement is analogous to continuous gain-scheduling or high-gain adaptive control within or across movement trials, respectively, and that the recently proposed inverse internal model is akin to adaptive sliding control originally for robotic manipulator applications. Modular internal models architecture for multiple motor tasks is a form of multi-model adaptive control. Stochastic methods such as generalized predictive control, reinforcement learning, Bayesian learning and Hebbian feedback covariance learning are reviewed and their possible relevance to motor control is discussed. Possible applicability of Luenberger observer and extended Kalman filter to state estimation problems such as sensorimotor prediction or the resolution of vestibular sensory ambiguity is also discussed. The important role played by vestibular system identification in postural control suggests an indirect adaptive control scheme whereby system states or parameters are explicitly estimated prior to the implementation of control. This interdisciplinary framework should facilitate the experimental elucidation of the mechanisms of internal model in sensorimotor systems and the reverse engineering of such neural mechanisms into novel brain-inspired adaptive control paradigms in future. PMID:16135881
Numerical modeling of landslide-generated tsunami using adaptive unstructured meshes
NASA Astrophysics Data System (ADS)
Wilson, Cian; Collins, Gareth; Desousa Costa, Patrick; Piggott, Matthew
2010-05-01
Landslides impacting into or occurring under water generate waves, which can have devastating environmental consequences. Depending on the characteristics of the landslide the waves can have significant amplitude and potentially propagate over large distances. Linear models of classical earthquake-generated tsunamis cannot reproduce the highly nonlinear generation mechanisms required to accurately predict the consequences of landslide-generated tsunamis. Also, laboratory-scale experimental investigation is limited to simple geometries and short time-scales before wave reflections contaminate the data. Computational fluid dynamics models based on the nonlinear Navier-Stokes equations can simulate landslide-tsunami generation at realistic scales. However, traditional chessboard-like structured meshes introduce superfluous resolution and hence the computing power required for such a simulation can be prohibitively high, especially in three dimensions. Unstructured meshes allow the grid spacing to vary rapidly from high resolution in the vicinity of small scale features to much coarser, lower resolution in other areas. Combining this variable resolution with dynamic mesh adaptivity allows such high resolution zones to follow features like the interface between the landslide and the water whilst minimising the computational costs. Unstructured meshes are also better suited to representing complex geometries and bathymetries allowing more realistic domains to be simulated. Modelling multiple materials, like water, air and a landslide, on an unstructured adaptive mesh poses significant numerical challenges. Novel methods of interface preservation must be considered and coupled to a flow model in such a way that ensures conservation of the different materials. Furthermore this conservation property must be maintained during successive stages of mesh optimisation and interpolation. In this paper we validate a new multi-material adaptive unstructured fluid dynamics model
Integrated numeric and symbolic signal processing using a heterogeneous design environment
NASA Astrophysics Data System (ADS)
Mani, Ramamurthy; Nawab, S. Hamid; Winograd, Joseph M.; Evans, Brian L.
1996-10-01
We present a solution to a complex multi-tone transient detection problem to illustrate the integrated use of symbolic and numeric processing techniques which are supported by well-established underlying models. Examples of such models include synchronous dataflow for numeric processing and the blackboard paradigm for symbolic heuristic search. Our transient detection solution serves to emphasize the importance of developing system design methods and tools which can support the integrated use of well- established symbolic and numerical models of computation. Recently, we incorporated a blackboard-based model of computation underlying the Integrated Processing and Understanding of Signals (IPUS) paradigm into a system-level design environment for numeric processing called Ptolemy. Using the IPUS/Ptolemy environment, we are implementing our solution to the multi-tone transient detection problem.
Numerical solution of a class of integral equations arising in two-dimensional aerodynamics
NASA Technical Reports Server (NTRS)
Fromme, J.; Golberg, M. A.
1978-01-01
We consider the numerical solution of a class of integral equations arising in the determination of the compressible flow about a thin airfoil in a ventilated wind tunnel. The integral equations are of the first kind with kernels having a Cauchy singularity. Using appropriately chosen Hilbert spaces, it is shown that the kernel gives rise to a mapping which is the sum of a unitary operator and a compact operator. This allows the problem to be studied in terms of an equivalent integral equation of the second kind. A convergent numerical algorithm for its solution is derived by using Galerkin's method. It is shown that this algorithm is numerically equivalent to Bland's collocation method, which is then used as the method of computation. Extensive numerical calculations are presented establishing the validity of the theory.
Adaptive Wavelet-Based Direct Numerical Simulations of Rayleigh-Taylor Instability
NASA Astrophysics Data System (ADS)
Reckinger, Scott J.
The compressible Rayleigh-Taylor instability (RTI) occurs when a fluid of low molar mass supports a fluid of higher molar mass against a gravity-like body force or in the presence of an accelerating front. Intrinsic to the problem are highly stratified background states, acoustic waves, and a wide range of physical scales. The objective of this thesis is to develop a specialized computational framework that addresses these challenges and to apply the advanced methodologies for direct numerical simulations of compressible RTI. Simulations are performed using the Parallel Adaptive Wavelet Collocation Method (PAWCM). Due to the physics-based adaptivity and direct error control of the method, PAWCM is ideal for resolving the wide range of scales present in RTI growth. Characteristics-based non-reflecting boundary conditions are developed for highly stratified systems to be used in conjunction with PAWCM. This combination allows for extremely long domains, which is necessary for observing the late time growth of compressible RTI. Initial conditions that minimize acoustic disturbances are also developed. The initialization is consistent with linear stability theory, where the background state consists of two diffusively mixed stratified fluids of differing molar masses. The compressibility effects on the departure from the linear growth, the onset of strong non-linear interactions, and the late-time behavior of the fluid structures are investigated. It is discovered that, for the thermal equilibrium case, the background stratification acts to suppress the instability growth when the molar mass difference is small. A reversal in this monotonic behavior is observed for large molar mass differences, where stratification enhances the bubble growth. Stratification also affects the vortex creation and the associated induced velocities. The enhancement and suppression of the RTI growth has important consequences for a detailed understanding of supernovae flame front
Vogel, Stephan E; Goffin, Celia; Ansari, Daniel
2015-04-01
The way the human brain constructs representations of numerical symbols is poorly understood. While increasing evidence from neuroimaging studies has indicated that the intraparietal sulcus (IPS) becomes increasingly specialized for symbolic numerical magnitude representation over developmental time, the extent to which these changes are associated with age-related differences in symbolic numerical magnitude representation or with developmental changes in non-numerical processes, such as response selection, remains to be uncovered. To address these outstanding questions we investigated developmental changes in the cortical representation of symbolic numerical magnitude in 6- to 14-year-old children using a passive functional magnetic resonance imaging adaptation design, thereby mitigating the influence of response selection. A single-digit Arabic numeral was repeatedly presented on a computer screen and interspersed with the presentation of novel digits deviating as a function of numerical ratio (smaller/larger number). Results demonstrated a correlation between age and numerical ratio in the left IPS, suggesting an age-related increase in the extent to which numerical symbols are represented in the left IPS. Brain activation of the right IPS was modulated by numerical ratio but did not correlate with age, indicating hemispheric differences in IPS engagement during the development of symbolic numerical representation. PMID:25555264
NASA Technical Reports Server (NTRS)
Sidi, A.; Israeli, M.
1986-01-01
High accuracy numerical quadrature methods for integrals of singular periodic functions are proposed. These methods are based on the appropriate Euler-Maclaurin expansions of trapezoidal rule approximations and their extrapolations. They are used to obtain accurate quadrature methods for the solution of singular and weakly singular Fredholm integral equations. Such periodic equations are used in the solution of planar elliptic boundary value problems, elasticity, potential theory, conformal mapping, boundary element methods, free surface flows, etc. The use of the quadrature methods is demonstrated with numerical examples.
NASA Astrophysics Data System (ADS)
Imada, Masatoshi; Kashima, Tsuyoshi
2000-09-01
A numerical algorithm for studying strongly correlated electron systems is proposed. The groundstate wavefunction is projected out after a numerical renormalization procedure in the path integral formalism. The wavefunction is expressed from the optimized linear combination of retained states in the truncated Hilbert space with a numerically chosen basis. This algorithm does not suffer from the negative sign problem and can be applied to any type of Hamiltonian in any dimension. The efficiency is tested in examples of the Hubbard model where the basis of Slater determinants is numerically optimized. We show results on fast convergence and accuracy achieved with a small number of retained states.
ERIC Educational Resources Information Center
Mitsuhara, Hiroyuki; Kurose, Yoshinobu; Ochi, Youji; Yano, Yoneo
The authors developed a Web-based Adaptive Educational System (Web-based AES) named ITMS (Individualized Teaching Material System). ITMS adaptively integrates knowledge on the distributed Web pages and generates individualized teaching material that has various contents. ITMS also presumes the learners' knowledge levels from the states of their…
Planning for an uncertain future - Monitoring, integration, and adaptation
Webb, Richard M. T., (Edited By); Semmens, Darius J.
2009-01-01
The 6.7 billion human inhabitants of the earth have the ability to drastically alter ecosystems and the populations of species that have taken eons to evolve. By better understanding how our actions affect the environment, we stand a better chance of designing successful strategies to manage ecosystems sustainably. Toward this end, the Third Interagency Conference on Research in the Watersheds (ICRW) was convened in Estes Park, CO, on September 8-11, 2008. The Conference provided a forum to present adaptive management as a practical tool for learning how to manage complex ecosystems more sustainably. Further complexity introduced by spatially variable and continuously changing environmental drivers favors this management approach because of its emphasis on adaptation in response to changing conditions or ineffective actions. For climate change in particular, an adaptive approach can more effectively accommodate the uncertainty in future climate scenarios. Scenarios compiled by the Intergovernmental Panel on Climate Change are built on distinct economic, energy, and societal models. The scenarios predict potential changes in greenhouse gases, temperature, precipitation, and atmospheric aerosols, which would have direct or indirect impacts on the timing, volume, and quality of runoff, vegetation, snowpack, stream temperature, groundwater, thawing permafrost, and icecaps. Through presentations and field trips, researchers and stakeholders described how their findings and issues fit into the adaptive management 'learning by doing' paradigm of Assess > Design > Implement > Monitor > Evaluate > Adjust > Assess.
Dissociating Conflict Adaptation from Feature Integration: A Multiple Regression Approach
ERIC Educational Resources Information Center
Notebaert, Wim; Verguts, Tom
2007-01-01
Congruency effects are typically smaller after incongruent than after congruent trials. One explanation is in terms of higher levels of cognitive control after detection of conflict (conflict adaptation; e.g., M. M. Botvinick, T. S. Braver, D. M. Barch, C. S. Carter, & J. D. Cohen, 2001). An alternative explanation for these results is based on…
A comparison of the efficiency of numerical methods for integrating chemical kinetic rate equations
NASA Technical Reports Server (NTRS)
Radhakrishnan, K.
1984-01-01
The efficiency of several algorithms used for numerical integration of stiff ordinary differential equations was compared. The methods examined included two general purpose codes EPISODE and LSODE and three codes (CHEMEQ, CREK1D and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes were applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code available for the integration of combustion kinetic rate equations. It is shown that an iterative solution of the algebraic energy conservation equation to compute the temperature can be more efficient then evaluating the temperature by integrating its time-derivative.
Controlled time integration for the numerical simulation of meteor radar reflections
NASA Astrophysics Data System (ADS)
Räbinä, Jukka; Mönkölä, Sanna; Rossi, Tuomo; Markkanen, Johannes; Gritsevich, Maria; Muinonen, Karri
2016-07-01
We model meteoroids entering the Earth's atmosphere as objects surrounded by non-magnetized plasma, and consider efficient numerical simulation of radar reflections from meteors in the time domain. Instead of the widely used finite difference time domain method (FDTD), we use more generalized finite differences by applying the discrete exterior calculus (DEC) and non-uniform leapfrog-style time discretization. The computational domain is presented by convex polyhedral elements. The convergence of the time integration is accelerated by the exact controllability method. The numerical experiments show that our code is efficiently parallelized. The DEC approach is compared to the volume integral equation (VIE) method by numerical experiments. The result is that both methods are competitive in modelling non-magnetized plasma scattering. For demonstrating the simulation capabilities of the DEC approach, we present numerical experiments of radar reflections and vary parameters in a wide range.
INCORPORATING CATASTROPHES INTO INTEGRATED ASSESSMENT: SCIENCE, IMPACTS, AND ADAPTATION
Incorporating potential catastrophic consequences into integrated assessment models of climate change has been a top priority of policymakers and modelers alike. We review the current state of scientific understanding regarding three frequently mentioned geophysical catastrophes,...
NASA Astrophysics Data System (ADS)
Alliss, R.
2014-09-01
Optical turbulence (OT) acts to distort light in the atmosphere, degrading imagery from astronomical telescopes and reducing the data quality of optical imaging and communication links. Some of the degradation due to turbulence can be corrected by adaptive optics. However, the severity of optical turbulence, and thus the amount of correction required, is largely dependent upon the turbulence at the location of interest. Therefore, it is vital to understand the climatology of optical turbulence at such locations. In many cases, it is impractical and expensive to setup instrumentation to characterize the climatology of OT, so numerical simulations become a less expensive and convenient alternative. The strength of OT is characterized by the refractive index structure function Cn2, which in turn is used to calculate atmospheric seeing parameters. While attempts have been made to characterize Cn2 using empirical models, Cn2 can be calculated more directly from Numerical Weather Prediction (NWP) simulations using pressure, temperature, thermal stability, vertical wind shear, turbulent Prandtl number, and turbulence kinetic energy (TKE). In this work we use the Weather Research and Forecast (WRF) NWP model to generate Cn2 climatologies in the planetary boundary layer and free atmosphere, allowing for both point-to-point and ground-to-space seeing estimates of the Fried Coherence length (ro) and other seeing parameters. Simulations are performed using a multi-node linux cluster using the Intel chip architecture. The WRF model is configured to run at 1km horizontal resolution and centered on the Mauna Loa Observatory (MLO) of the Big Island. The vertical resolution varies from 25 meters in the boundary layer to 500 meters in the stratosphere. The model top is 20 km. The Mellor-Yamada-Janjic (MYJ) TKE scheme has been modified to diagnose the turbulent Prandtl number as a function of the Richardson number, following observations by Kondo and others. This modification
Data rate management and real time operation: recursive adaptive frame integration of limited data
NASA Astrophysics Data System (ADS)
Rafailov, Michael K.
2006-08-01
Recursive Limited Frame Integration was proposed as a way to improve frame integration performance and mitigate issues related to high data rate needed to support conventional frame integration. The technique uses two thresholds -one tuned for optimum probability of detection, the other to manage required false alarm rate, and places integration process between those thresholds. This configuration allows a non-linear integration process that, along with Signal-to-Noise Ratio (SNR) gain, provides system designers more capability where cost, weight, or power considerations limit system data rate, processing, or memory capability. However, Recursive Frame Integration Limited may have performance issues when single-frame SNR is really low. Recursive Adaptive Limited Frame Integration was proposed as a means to improve limited integration performance with really low single-frame SNR. It combines the benefits of nonlinear recursive limited frame integration and adaptive thresholds with a kind of conventional frame integration. Adding the third threshold may help in managing real time operations. In the paper the Recursive Frame Integration is presented in form of multiple parallel recursive integration. Such an approach can help not only in data rate management but in mitigation of low single frame SNR issue for Recursive Integration as well as in real time operations with frame integration.
Romá, Federico; Cugliandolo, Leticia F; Lozano, Gustavo S
2014-08-01
We introduce a numerical method to integrate the stochastic Landau-Lifshitz-Gilbert equation in spherical coordinates for generic discretization schemes. This method conserves the magnetization modulus and ensures the approach to equilibrium under the expected conditions. We test the algorithm on a benchmark problem: the dynamics of a uniformly magnetized ellipsoid. We investigate the influence of various parameters, and in particular, we analyze the efficiency of the numerical integration, in terms of the number of steps needed to reach a chosen long time with a given accuracy. PMID:25215839
Toyoda, Masayuki; Ozaki, Taisuke
2009-03-28
A numerical method to calculate the four-center electron-repulsion integrals for strictly localized pseudoatomic orbital basis sets has been developed. Compared to the conventional Gaussian expansion method, this method has an advantage in the ease of combination with O(N) density functional calculations. Additional mathematical derivations are also presented including the analytic derivatives of the integrals with respect to atomic positions and spatial damping of the Coulomb interaction due to the screening effect. In the numerical test for a simple molecule, the convergence up to 10(-5) hartree in energy is successfully obtained with a feasible cost of computation. PMID:19334815
An integrated architecture of adaptive neural network control for dynamic systems
Ke, Liu; Tokar, R.; Mcvey, B.
1994-07-01
In this study, an integrated neural network control architecture for nonlinear dynamic systems is presented. Most of the recent emphasis in the neural network control field has no error feedback as the control input which rises the adaptation problem. The integrated architecture in this paper combines feed forward control and error feedback adaptive control using neural networks. The paper reveals the different internal functionality of these two kinds of neural network controllers for certain input styles, e.g., state feedback and error feedback. Feed forward neural network controllers with state feedback establish fixed control mappings which can not adapt when model uncertainties present. With error feedbacks, neural network controllers learn the slopes or the gains respecting to the error feedbacks, which are error driven adaptive control systems. The results demonstrate that the two kinds of control scheme can be combined to realize their individual advantages. Testing with disturbances added to the plant shows good tracking and adaptation.
NASA Astrophysics Data System (ADS)
Sagert, I.; Fann, G. I.; Fattoyev, F. J.; Postnikov, S.; Horowitz, C. J.
2016-05-01
Background: Neutron star and supernova matter at densities just below the nuclear matter saturation density is expected to form a lattice of exotic shapes. These so-called nuclear pasta phases are caused by Coulomb frustration. Their elastic and transport properties are believed to play an important role for thermal and magnetic field evolution, rotation, and oscillation of neutron stars. Furthermore, they can impact neutrino opacities in core-collapse supernovae. Purpose: In this work, we present proof-of-principle three-dimensional (3D) Skyrme Hartree-Fock (SHF) simulations of nuclear pasta with the Multi-resolution ADaptive Numerical Environment for Scientific Simulations (MADNESS). Methods: We perform benchmark studies of 16O, 208Pb, and 238U nuclear ground states and calculate binding energies via 3D SHF simulations. Results are compared with experimentally measured binding energies as well as with theoretically predicted values from an established SHF code. The nuclear pasta simulation is initialized in the so-called waffle geometry as obtained by the Indiana University Molecular Dynamics (IUMD) code. The size of the unit cell is 24 fm with an average density of about ρ =0.05 fm-3 , proton fraction of Yp=0.3 , and temperature of T =0 MeV. Results: Our calculations reproduce the binding energies and shapes of light and heavy nuclei with different geometries. For the pasta simulation, we find that the final geometry is very similar to the initial waffle state. We compare calculations with and without spin-orbit forces. We find that while subtle differences are present, the pasta phase remains in the waffle geometry. Conclusions: Within the MADNESS framework, we can successfully perform calculations of inhomogeneous nuclear matter. By using pasta configurations from IUMD it is possible to explore different geometries and test the impact of self-consistent calculations on the latter.
NASA Astrophysics Data System (ADS)
Parkinson, S. D.; Hill, J.; Piggott, M. D.; Allison, P. A.
2014-05-01
High resolution direct numerical simulations (DNS) are an important tool for the detailed analysis of turbidity current dynamics. Models that resolve the vertical structure and turbulence of the flow are typically based upon the Navier-Stokes equations. Two-dimensional simulations are known to produce unrealistic cohesive vortices that are not representative of the real three-dimensional physics. The effect of this phenomena is particularly apparent in the later stages of flow propagation. The ideal solution to this problem is to run the simulation in three dimensions but this is computationally expensive. This paper presents a novel finite-element (FE) DNS turbidity current model that has been built within Fluidity, an open source, general purpose, computational fluid dynamics code. The model is validated through re-creation of a lock release density current at a Grashof number of 5 × 106 in two, and three-dimensions. Validation of the model considers the flow energy budget, sedimentation rate, head speed, wall normal velocity profiles and the final deposit. Conservation of energy in particular is found to be a good metric for measuring mesh performance in capturing the range of dynamics. FE models scale well over many thousands of processors and do not impose restrictions on domain shape, but they are computationally expensive. Use of discontinuous discretisations and adaptive unstructured meshing technologies, which reduce the required element count by approximately two orders of magnitude, results in high resolution DNS models of turbidity currents at a fraction of the cost of traditional FE models. The benefits of this technique will enable simulation of turbidity currents in complex and large domains where DNS modelling was previously unachievable.
Numerical integration of discontinuities on arbitrary domains based on moment fitting
NASA Astrophysics Data System (ADS)
Joulaian, Meysam; Hubrich, Simeon; Düster, Alexander
2016-03-01
Discretization methods based on meshes that do not conform to the geometry of the problem under consideration require special treatment when it comes to the integration of finite elements that are broken by the boundary or internal interfaces. To this end, we propose a numerical approach suitable for integrating broken elements with a low number of integration points. In this method, which is based on the moment fitting approach, an individual quadrature rule is set up for each cut element. The approach requires a B-rep representation of the broken element, which can be either achieved by processing a triangulated surface obtained from a CAD software or by taking advantage of a voxel model resulting from computed tomography. The numerical examples presented in this paper reveal that the proposed method delivers for a wide variety of geometrical situations very accurate results and requires a rather low number of integration points.
Numerical solutions to ill-posed and well-posed impedance boundary condition integral equations
NASA Astrophysics Data System (ADS)
Rogers, J. R.
1983-11-01
Exterior scattering from a three-dimensional impedance body can be formulated in terms of various integral equations derived from the Leontovich impedance boundary condition (IBC). The electric and magnetic field integral equations are ill-posed because they theoretically admit spurious solutions at the frequencies of interior perfect conductor cavity resonances. A combined field formulation is well-posed because it does not allow the spurious solutions. This report outlines the derivation of IBC integral equations and describes a procedure for constructing moment-method solutions for bodies of revolution. Numerical results for scattering from impedance spheres are presented which contrast the stability and accuracy of solutions to the ill-posed equations with those of the well-posed equation. The results show that numerical solutions for exterior scattering to the electric and magnetic field integral equations can be severely contaminated by spurious resonant solutions regardless of whether the surface impedance of the body is lossy or lossless.
Numerical integration of discontinuities on arbitrary domains based on moment fitting
NASA Astrophysics Data System (ADS)
Joulaian, Meysam; Hubrich, Simeon; Düster, Alexander
2016-06-01
Discretization methods based on meshes that do not conform to the geometry of the problem under consideration require special treatment when it comes to the integration of finite elements that are broken by the boundary or internal interfaces. To this end, we propose a numerical approach suitable for integrating broken elements with a low number of integration points. In this method, which is based on the moment fitting approach, an individual quadrature rule is set up for each cut element. The approach requires a B-rep representation of the broken element, which can be either achieved by processing a triangulated surface obtained from a CAD software or by taking advantage of a voxel model resulting from computed tomography. The numerical examples presented in this paper reveal that the proposed method delivers for a wide variety of geometrical situations very accurate results and requires a rather low number of integration points.
ERIC Educational Resources Information Center
Rule, Audrey C.; Barrera, Manuel T., III
2008-01-01
Integration of subject areas with technology and thinking skills is a way to help teachers cope with today's overloaded curriculum and to help students see the connectedness of different curriculum areas. This study compares three authentic approaches to teaching a science unit on bird adaptations for habitat that integrate thinking skills and…
ERIC Educational Resources Information Center
Yu, Baohua; Downing, Kevin
2012-01-01
This study examined the influence of integrative motivation, instrumental motivation and second language (L2) proficiency on socio-cultural/academic adaptation in a sample of two groups of international students studying Chinese in China. Results revealed that the non-Asian student group reported higher levels of integrative motivation,…
NASA Astrophysics Data System (ADS)
Ma, Shaokang; Wu, Peijun; Ji, Jinhu; Li, Xuchun
2016-02-01
This article presents a sensorless control approach of salient PMSM with an online parameter identifier. Adaptive Integrator is proposed and utilised for the estimation of active flux and rotor position. As a result, integrator overflow caused by DC offset is avoided. Meanwhile, an online stator resistance identification algorithm using strong tracking filter is employed, and the identified stator resistance is fed back to the estimating algorithm. Thus, the estimating algorithm can calculate the rotor position correctly. Simulations and experimental results validate the feasibility of both adaptive integrator and the parameter identification method.
High-performance Integrated numerical methods for Two-phase Flow in Heterogeneous Porous Media
NASA Astrophysics Data System (ADS)
Chueh, Chih-Che; Djilali, Ned; Bangerth, Wolfgang
2010-11-01
Modelling of two-phase flow in heterogeneous porous media has been playing a decisive role in a variety of areas. However, how to efficiently and accurately solve the governing equation in the flow in porous media remains a challenge. In order to ensure the accurate representative flow field and simultaneously increase the computational efficiency, we incorporate a number of state-of-the-art techniques into a numerical framework on which more complicated models in the field of multi-phase flow in porous media will be based. Such a numerical framework consists of a h-adaptive refinement method, an entropy-based artificial diffusive term, a new adaptive operator splitting method and efficient preconditioners. In particular, it is emphasized that we propose a new efficient adaptive operator splitting to avoid solving a time-consuming pressure-velocity part every saturation time step and, most importantly, we also provide a theoretically numerical analysis as well as proof. A few benchmarks will be demonstrated in the presentation.
NASA Astrophysics Data System (ADS)
Rodriguez Bueno, R. A.; Byrne, J. M.
2015-12-01
The Environment Service Center of Matanzas (ESCM), Cuba and the University of Lethbridge are collaborating on the development of climate mitigation and adaptation programs in Matanzas province. Tourism is the largest industry in Matanzas. Protecting that industry means protecting coastal zones and conservation areas of value to tourism. These same areas are critical to protecting the landscape from global environmental change: enhanced tropical cyclones, flooding, drought and a range of other environmental change impacts. Byrne (2014) adapted a multidisciplinary methodology for climate adaptation capacity definition for the population of Nicaragua. A wide array of adaptive capacity skills and resources were integrated with agricultural crop modeling to define regions of the country where adaptive capacity development were weakest and should be improved. In Matanzas province, we are developing a series of multidisciplinary mitigation and adaptation programs that builds social science and science knowledge to expand capacity within the ESCM and the provincial population. We will be exploring increased risk due to combined watershed and tropical cyclone flooding, stresses on crops, and defining a range of possibilities in shifting from fossil fuels to renewable energy. The program will build ongoing interactions with thousands of Matanzas citizens through site visits carried out by numerous Cuban and visiting students participating in a four-month education semester with a number of Lethbridge and Matanzas faculty. These visits will also provide local citizens with better access to web-based interactions. We will evaluate mitigation and adaptive capacities in three municipalities and some rural areas across the province. Furthermore, we will explore better ways and means to communicate between the research and conservation staff and the larger population of the province.
Peskin, Michael E
2003-02-13
In upper-division undergraduate physics courses, it is desirable to give numerical problem-solving exercises integrated naturally into weekly problem sets. I explain a method for doing this that makes use of the built-in class structure of the Java programming language. I also supply a Java class library that can assist instructors in writing programs of this type.
NUMERICAL APPROXIMATION OF SEMI-INTEGRALS AND SEMIDERIVATIVES BY PRODUCT QUADRATURE RULES
This paper is concerned with the numerical calculation of the semi-integral and semiderivative of a function f, whose values f (xj) are known on a discrete set of abscissas 0 = x(1) < x(2) < ... < x(n). A family of product quadrature rules is developed to approximate the semi-int...
Some numerical methods for integrating systems of first-order ordinary differential equations
NASA Technical Reports Server (NTRS)
Clark, N. W.
1969-01-01
Report on numerical methods of integration includes the extrapolation methods of Bulirsch-Stoer and Neville. A comparison is made nith the Runge-Kutta and Adams-Moulton methods, and circumstances are discussed under which the extrapolation method may be preferred.
On the stability of numerical integration routines for ordinary differential equations.
NASA Technical Reports Server (NTRS)
Glover, K.; Willems, J. C.
1973-01-01
Numerical integration methods for the solution of initial value problems for ordinary vector differential equations may be modelled as discrete time feedback systems. The stability criteria discovered in modern control theory are applied to these systems and criteria involving the routine, the step size and the differential equation are derived. Linear multistep, Runge-Kutta, and predictor-corrector methods are all investigated.
NASA Astrophysics Data System (ADS)
Tang, Xiaojun
2016-04-01
The main purpose of this work is to provide multiple-interval integral Gegenbauer pseudospectral methods for solving optimal control problems. The latest developed single-interval integral Gauss/(flipped Radau) pseudospectral methods can be viewed as special cases of the proposed methods. We present an exact and efficient approach to compute the mesh pseudospectral integration matrices for the Gegenbauer-Gauss and flipped Gegenbauer-Gauss-Radau points. Numerical results on benchmark optimal control problems confirm the ability of the proposed methods to obtain highly accurate solutions.
Rapid adaptation of multisensory integration in vestibular pathways.
Carriot, Jerome; Jamali, Mohsen; Cullen, Kathleen E
2015-01-01
Sensing gravity is vital for our perception of spatial orientation, the control of upright posture, and generation of our everyday activities. When an astronaut transitions to microgravity or returns to earth, the vestibular input arising from self-motion will not match the brain's expectation. Our recent neurophysiological studies have provided insight into how the nervous system rapidly reorganizes when vestibular input becomes unreliable by both (1) updating its internal model of the sensory consequences of motion and (2) up-weighting more reliable extra-vestibular information. These neural strategies, in turn, are linked to improvements in sensorimotor performance (e.g., gaze and postural stability, locomotion, orienting) and perception characterized by similar time courses. We suggest that furthering our understanding of the neural mechanisms that underlie sensorimotor adaptation will have important implications for optimizing training programs for astronauts before and after space exploration missions and for the design of goal-oriented rehabilitation for patients. PMID:25932009
Rapid adaptation of multisensory integration in vestibular pathways
Carriot, Jerome; Jamali, Mohsen; Cullen, Kathleen E.
2015-01-01
Sensing gravity is vital for our perception of spatial orientation, the control of upright posture, and generation of our everyday activities. When an astronaut transitions to microgravity or returns to earth, the vestibular input arising from self-motion will not match the brain's expectation. Our recent neurophysiological studies have provided insight into how the nervous system rapidly reorganizes when vestibular input becomes unreliable by both (1) updating its internal model of the sensory consequences of motion and (2) up-weighting more reliable extra-vestibular information. These neural strategies, in turn, are linked to improvements in sensorimotor performance (e.g., gaze and postural stability, locomotion, orienting) and perception characterized by similar time courses. We suggest that furthering our understanding of the neural mechanisms that underlie sensorimotor adaptation will have important implications for optimizing training programs for astronauts before and after space exploration missions and for the design of goal-oriented rehabilitation for patients. PMID:25932009
Induction, adaptation and recovery of lysosomal integrity in green-lipped mussel Perna viridis.
Fang, J K H; Wu, R S S; Zheng, G J; Lam, P K S; Shin, P K S
2008-01-01
Biomarkers are generally applied to detect pollution in environmental monitoring. Such biological responses should accurately reflect the stress over time in a quantitative manner. As such, the initial and maximum responses induced by stress, as well as adaptation and recovery of these biomarkers, need to be fully understood or else erroneous false-negative or false-positive may be arrived. However, most of the biomarker studies only provided information on initially induced responses under different concentrations of toxicants, while biological adaptation and recovery were poorly known. In this study, the time required for induction, adaptation and recovery of lysosomal integrity in green-lipped mussel Perna viridis upon exposure to benzo[a]pyrene was investigated over a period of 62 days. Maximum induction occurred on day 6 when lysosomal integrity was significantly reduced by 51%, and no further change or adaptation was detected thereafter. When mussels were depurated in clean seawater after 18 days of exposure to benzo[a]pyrene, a gradual recovery was observed, with lysosomal integrity returning to its background level and showing a complete recovery after 20 days of depuration. Lysosomal integrity was significantly correlated with the body burden concentrations of benzo[a]pyrene and condition index of the mussels. The relatively fast induction (6 days) and recovery (20 days) without apparent adaptation suggested that lysosomal integrity in P. viridis can serve as a good biomarker in biomonitoring, as its response is not likely to generate both false-negative and false-positive results. PMID:18466928
NASA Astrophysics Data System (ADS)
Zhang, Feng; Duan, Guangren
2013-05-01
This paper handles the integrated translation and rotation tracking control problem of a rigid spacecraft with unknown mass property, actuator misalignment and control saturation. In view of the system natural coupling, the coupled translational and rotational dynamics of the spacecraft is developed, where a thruster configuration with installation misalignment is taken into account. By using anti-windup technique and backstepping philosophy, a robust adaptive integrated control scheme is proposed such that the spacecraft is able to track the command position and attitude signals in the presence of external disturbance, unknown mass property, thruster misalignment and control saturation. Within the Lyapunov framework, the uniformly ultimate boundedness of the system states is guaranteed. In particular, given the nominal case, the asymptotic convergence of the system states can be further ensured by the proposed control scheme. Finally, numerical simulation demonstrates the effect of the designed control strategy.
Magnitude Estimation with Noisy Integrators Linked by an Adaptive Reference
Thurley, Kay
2016-01-01
Judgments of physical stimuli show characteristic biases; relatively small stimuli are overestimated whereas relatively large stimuli are underestimated (regression effect). Such biases likely result from a strategy that seeks to minimize errors given noisy estimates about stimuli that itself are drawn from a distribution, i.e., the statistics of the environment. While being conceptually well described, it is unclear how such a strategy could be implemented neurally. The present paper aims toward answering this question. A theoretical approach is introduced that describes magnitude estimation as two successive stages of noisy (neural) integration. Both stages are linked by a reference memory that is updated with every new stimulus. The model reproduces the behavioral characteristics of magnitude estimation and makes several experimentally testable predictions. Moreover, the model identifies the regression effect as a means of minimizing estimation errors and explains how this optimality strategy depends on the subject's discrimination abilities and on the stimulus statistics. The latter influence predicts another property of magnitude estimation, the so-called range effect. Beyond being successful in describing decision-making, the present work suggests that noisy integration may also be important in processing magnitudes. PMID:26909028
Adapting line integral convolution for fabricating artistic virtual environment
NASA Astrophysics Data System (ADS)
Lee, Jiunn-Shyan; Wang, Chung-Ming
2003-04-01
Vector field occurs not only extensively in scientific applications but also in treasured art such as sculptures and paintings. Artist depicts our natural environment stressing valued directional feature besides color and shape information. Line integral convolution (LIC), developed for imaging vector field in scientific visualization, has potential of producing directional image. In this paper we present several techniques of exploring LIC techniques to generate impressionistic images forming artistic virtual environment. We take advantage of directional information given by a photograph, and incorporate many investigations to the work including non-photorealistic shading technique and statistical detail control. In particular, the non-photorealistic shading technique blends cool and warm colors into the photograph to imitate artists painting convention. Besides, we adopt statistical technique controlling integral length according to image variance to preserve details. Furthermore, we also propose method for generating a series of mip-maps, which revealing constant strokes under multi-resolution viewing and achieving frame coherence in an interactive walkthrough system. The experimental results show merits of emulating satisfyingly and computing efficiently, as a consequence, relying on the proposed technique successfully fabricates a wide category of non-photorealistic rendering (NPR) application such as interactive virtual environment with artistic perception.
Magnitude Estimation with Noisy Integrators Linked by an Adaptive Reference.
Thurley, Kay
2016-01-01
Judgments of physical stimuli show characteristic biases; relatively small stimuli are overestimated whereas relatively large stimuli are underestimated (regression effect). Such biases likely result from a strategy that seeks to minimize errors given noisy estimates about stimuli that itself are drawn from a distribution, i.e., the statistics of the environment. While being conceptually well described, it is unclear how such a strategy could be implemented neurally. The present paper aims toward answering this question. A theoretical approach is introduced that describes magnitude estimation as two successive stages of noisy (neural) integration. Both stages are linked by a reference memory that is updated with every new stimulus. The model reproduces the behavioral characteristics of magnitude estimation and makes several experimentally testable predictions. Moreover, the model identifies the regression effect as a means of minimizing estimation errors and explains how this optimality strategy depends on the subject's discrimination abilities and on the stimulus statistics. The latter influence predicts another property of magnitude estimation, the so-called range effect. Beyond being successful in describing decision-making, the present work suggests that noisy integration may also be important in processing magnitudes. PMID:26909028
NASA Astrophysics Data System (ADS)
Zhou, Yanlai; Guo, Shenglian; Xu, Chong-Yu; Liu, Dedi; Chen, Lu; Ye, Yushi
2015-12-01
Due to the adaption, dynamic and multi-objective characteristics of complex water resources system, it is a considerable challenge to manage water resources in an efficient, equitable and sustainable way. An integrated optimal allocation model is proposed for complex adaptive system of water resources management. The model consists of three modules: (1) an agent-based module for revealing evolution mechanism of complex adaptive system using agent-based, system dynamic and non-dominated sorting genetic algorithm II methods, (2) an optimal module for deriving decision set of water resources allocation using multi-objective genetic algorithm, and (3) a multi-objective evaluation module for evaluating the efficiency of the optimal module and selecting the optimal water resources allocation scheme using project pursuit method. This study has provided a theoretical framework for adaptive allocation, dynamic allocation and multi-objective optimization for a complex adaptive system of water resources management.
Numerous strategies but limited implementation guidance in US local adaptation plans
NASA Astrophysics Data System (ADS)
Woodruff, Sierra C.; Stults, Missy
2016-08-01
Adaptation planning offers a promising approach for identifying and devising solutions to address local climate change impacts. Yet there is little empirical understanding of the content and quality of these plans. We use content analysis to evaluate 44 local adaptation plans in the United States and multivariate regression to examine how plan quality varies across communities. We find that plans draw on multiple data sources to analyse future climate impacts and include a breadth of strategies. Most plans, however, fail to prioritize impacts and strategies or provide detailed implementation processes, raising concerns about whether adaptation plans will translate into on-the-ground reductions in vulnerability. Our analysis also finds that plans authored by the planning department and those that engaged elected officials in the planning process were of higher quality. The results provide important insights for practitioners, policymakers and scientists wanting to improve local climate adaptation planning and action.
Numerical Modeling of Long Bone Adaptation due to Mechanical Loading: Correlation with Experiments
Kumar, Natarajan Chennimalai; Dantzig, Jonathan A.; Jasiuk, Iwona M.; Robling, Alex G.; Turner, Charles H.
2011-01-01
The process of external bone adaptation in cortical bone is modeled mathematically using finite element (FE) stress analysis coupled with an evolution model, in which adaptation response is triggered by mechanical stimulus represented by strain energy density. The model is applied to experiments in which a rat ulna is subjected to cyclic loading, and the results demonstrate the ability of the model to predict the bone adaptation response. The FE mesh is generated from micro-computed tomography (μCT) images of the rat ulna, and the stress analysis is carried out using boundary and loading conditions on the rat ulna obtained from the experiments [Robling, A. G., F. M. Hinant, D. B. Burr, and C. H. Turner. J. Bone Miner. Res. 17:1545–1554, 2002]. The external adaptation process is implemented in the model by moving the surface nodes of the FE mesh based on an evolution law characterized by two parameters: one that captures the rate of the adaptation process (referred to as gain); and the other characterizing the threshold value of the mechanical stimulus required for adaptation (referred to as threshold-sensitivity). A parametric study is carried out to evaluate the effect of these two parameters on the adaptation response. We show, following comparison of results from the simulations to the experimental observations of Robling et al. (J. Bone Miner. Res. 17:1545–1554, 2002), that splitting the loading cycles into different number of bouts affects the threshold-sensitivity but not the rate of adaptation. We also show that the threshold-sensitivity parameter can quantify the mechanosensitivity of the osteocytes. PMID:20013156
INTEGRATING EVOLUTIONARY AND FUNCTIONAL APPROACHES TO INFER ADAPTATION AT SPECIFIC LOCI
Storz, Jay F.; Wheat, Christopher W.
2010-01-01
Inferences about adaptation at specific loci are often exclusively based on the static analysis of DNA sequence variation. Ideally, population-genetic evidence for positive selection serves as a stepping-off point for experimental studies to elucidate the functional significance of the putatively adaptive variation. We argue that inferences about adaptation at specific loci are best achieved by integrating the indirect, retrospective insights provided by population-genetic analyses with the more direct, mechanistic insights provided by functional experiments. Integrative studies of adaptive genetic variation may sometimes be motivated by experimental insights into molecular function, which then provide the impetus to perform population genetic tests to evaluate whether the functional variation is of adaptive significance. In other cases, studies may be initiated by genome scans of DNA variation to identify candidate loci for recent adaptation. Results of such analyses can then motivate experimental efforts to test whether the identified candidate loci do in fact contribute to functional variation in some fitness-related phenotype. Functional studies can provide corroborative evidence for positive selection at particular loci, and can potentially reveal specific molecular mechanisms of adaptation. PMID:20500215
NASA Astrophysics Data System (ADS)
Evans, W. A. B.; Torre, A.
2012-11-01
The paper focusses on the advantages of using high-order Gauss-Legendre quadratures for the precise evaluation of integrals with both smooth and rapidly changing integrands. Aspects of their precision are analysed in the light of Gauss' error formula. Some "test examples" are considered and evaluated in multiple precision to ≈ 200 significant decimal digits with David Bailey's multiprecision package to eliminate truncation/rounding errors. The increase of precision on doubling the number of subintervals is analysed, the relevant quadrature attribute being the precision increment. In order to exemplify the advantages that high-order quadrature afford, the technique is then used to evaluate several plots of the Rayleigh-Sommerfeld diffraction integral for axi-symmetric source fields defined on a planar aperture. A comparison of the high-order quadrature method against various FFT-based methods is finally given.
NASA Astrophysics Data System (ADS)
Pantano, Carlos
2005-11-01
We describe a hybrid finite difference method for large-eddy simulation (LES) of compressible flows with a low-numerical dissipation scheme and structured adaptive mesh refinement (SAMR). Numerical experiments and validation calculations are presented including a turbulent jet and the strongly shock-driven mixing of a Richtmyer-Meshkov instability. The approach is a conservative flux-based SAMR formulation and as such, it utilizes refinement to computational advantage. The numerical method for the resolved scale terms encompasses the cases of scheme alternation and internal mesh interfaces resulting from SAMR. An explicit centered scheme that is consistent with a skew-symmetric finite difference formulation is used in turbulent flow regions while a weighted essentially non-oscillatory (WENO) scheme is employed to capture shocks. The subgrid stresses and transports are calculated by means of the streched-vortex model, Misra & Pullin (1997)
NASA Astrophysics Data System (ADS)
Pantano, C.; Deiterding, R.; Hill, D. J.; Pullin, D. I.
2007-01-01
We present a methodology for the large-eddy simulation of compressible flows with a low-numerical dissipation scheme and structured adaptive mesh refinement (SAMR). A description of a conservative, flux-based hybrid numerical method that uses both centered finite-difference and a weighted essentially non-oscillatory (WENO) scheme is given, encompassing the cases of scheme alternation and internal mesh interfaces resulting from SAMR. In this method, the centered scheme is used in turbulent flow regions while WENO is employed to capture shocks. One-, two- and three-dimensional numerical experiments and example simulations are presented including homogeneous shock-free turbulence, a turbulent jet and the strongly shock-driven mixing of a Richtmyer-Meshkov instability.
NASA Technical Reports Server (NTRS)
Rosenbaum, J. S.
1976-01-01
If a system of ordinary differential equations represents a property conserving system that can be expressed linearly (e.g., conservation of mass), it is then desirable that the numerical integration method used conserve the same quantity. It is shown that both linear multistep methods and Runge-Kutta methods are 'conservative' and that Newton-type methods used to solve the implicit equations preserve the inherent conservation of the numerical method. It is further shown that a method used by several authors is not conservative.
A comparison of the efficiency of numerical methods for integrating chemical kinetic rate equations
NASA Technical Reports Server (NTRS)
Radhakrishnan, K.
1984-01-01
A comparison of the efficiency of several algorithms recently developed for the efficient numerical integration of stiff ordinary differential equations is presented. The methods examined include two general-purpose codes EPISODE and LSODE and three codes (CHEMEQ, CREK1D, and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes are applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code currently available for the integration of combustion kinetic rate equations. An important finding is that an iterative solution of the algebraic energy conservation equation to compute the temperature can be more efficient than evaluating the temperature by integrating its time-derivative.
Sea Extremes: Integrated impact assessment in coastal climate adaptation
NASA Astrophysics Data System (ADS)
Sorensen, Carlo; Knudsen, Per; Broge, Niels; Molgaard, Mads; Andersen, Ole
2016-04-01
We investigate effects of sea level rise and a change in precipitation pattern on coastal flooding hazards. Historic and present in situ and satellite data of water and groundwater levels, precipitation, vertical ground motion, geology, and geotechnical soil properties are combined with flood protection measures, topography, and infrastructure to provide a more complete picture of the water-related impact from climate change at an exposed coastal location. Results show that future sea extremes evaluated from extreme value statistics may, indeed, have a large impact. The integrated effects from future storm surges and other geo- and hydro-parameters need to be considered in order to provide for the best protection and mitigation efforts, however. Based on the results we present and discuss a simple conceptual model setup that can e.g. be used for 'translation' of regional sea level rise evidence and projections to concrete impact measures. This may be used by potentially affected stakeholders -often working in different sectors and across levels of governance, in a common appraisal of the challenges faced ahead. The model may also enter dynamic tools to evaluate local impact as sea level research advances and projections for the future are updated.
Comparison of symbolic and numerical integration methods for an assumed-stress hybrid shell element
NASA Technical Reports Server (NTRS)
Rengarajan, Govind; Knight, Norman F., Jr.; Aminpour, Mohammad A.
1993-01-01
Hybrid shell elements have long been regarded with reserve by the commercial finite element developers despite the high degree of reliability and accuracy associated with such formulations. The fundamental reason is the inherent higher computational cost of the hybrid approach as compared to the displacement-based formulations. However, a noteworthy factor in favor of hybrid elements is that numerical integration to generate element matrices can be entirely avoided by the use of symbolic integration. In this paper, the use of the symbolic computational approach is presented for an assumed-stress hybrid shell element with drilling degrees of freedom and the significant time savings achieved is demonstrated through an example.
NLOS UV channel modeling using numerical integration and an approximate closed-form path loss model
NASA Astrophysics Data System (ADS)
Gupta, Ankit; Noshad, Mohammad; Brandt-Pearce, Maïté
2012-10-01
In this paper we propose a simulation method using numerical integration, and develop a closed-form link loss model for physical layer channel characterization for non-line of sight (NLOS) ultraviolet (UV) communication systems. The impulse response of the channel is calculated by assuming both uniform and Gaussian profiles for transmitted beams and different geometries. The results are compared with previously published results. The accuracy of the integration approach is compared to the Monte Carlo simulation. Then the path loss using the simulation method and the suggested closed-form expression are presented for different link geometries. The accuracies are evaluated and compared to the results obtained using other methods.
Numerical evaluation of two-center integrals over Slater type orbitals
NASA Astrophysics Data System (ADS)
Kurt, S. A.; Yükçü, N.
2016-03-01
Slater Type Orbitals (STOs) which one of the types of exponential type orbitals (ETOs) are used usually as basis functions in the multicenter molecular integrals to better understand physical and chemical properties of matter. In this work, we develop algorithms for two-center overlap and two-center two-electron hybrid and Coulomb integrals which are calculated with help of translation method for STOs and some auxiliary functions by V. Magnasco's group. We use Mathematica programming language to produce algorithms for these calculations. Numerical results for some quantum numbers are presented in the tables. Consequently, we compare our obtained numerical results with the other known literature results and other details of evaluation method are discussed.
Extended RKN-type methods for numerical integration of perturbed oscillators
NASA Astrophysics Data System (ADS)
Yang, Hongli; Wu, Xinyuan; You, Xiong; Fang, Yonglei
2009-10-01
In this paper, extended Runge-Kutta-Nyström-type methods for the numerical integration of perturbed oscillators with low frequencies are presented, which inherit the framework of RKN methods and make full use of the special feature of the true flows for both the internal stages and the updates. Following the approach of J. Butcher, E. Hairer and G. Wanner, we develop a new kind of tree set to derive order conditions for the extended Runge-Kutta-Nyström-type methods. The numerical stability and phase properties of the new methods are analyzed. Numerical experiments are accompanied to show the applicability and efficiency of our new methods in comparison with some well-known high quality methods proposed in the scientific literature.
A two-dimensional depth-integrated non-hydrostatic numerical model for nearshore wave propagation
NASA Astrophysics Data System (ADS)
Lu, Xinhua; Dong, Bingjiang; Mao, Bing; Zhang, Xiaofeng
2015-12-01
In this study, we develop a shallow-water depth-integrated non-hydrostatic numerical model (SNH model) using a hybrid finite-volume and finite-difference method. Numerical discretization is performed using the non-incremental pressure-correction method on a collocated grid. We demonstrate that an extension can easily be made from an existing finite-volume method and collocated-grid based hydrostatic shallow-water equations (SWE) model to a non-hydrostatic model. A series of benchmark tests are used to validate the proposed numerical model. Our results demonstrate that the proposed model is robust and well-balanced, and it captures the wet-dry fronts accurately. A comparison between the SNH and SWE models indicates the importance of considering the wave dispersion effect in simulations when the wave amplitude to water depth ratio is large.
Two step hybrid methods of 7th and 8th order for the numerical integration of second order IVPs
NASA Astrophysics Data System (ADS)
Kalogiratou, Z.; Monovasilis, Th.; 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 two step hybrid methods with six and seven stages and seventh and eighth algebraic order. We apply the new methods on the numerical integration of several test problems.
Some remarks on the numerical computation of integrals on an unbounded interval
NASA Astrophysics Data System (ADS)
Capobianco, M.; Criscuolo, G.
2007-08-01
An account of the error and the convergence theory is given for Gauss?Laguerre and Gauss?Radau?Laguerre quadrature formulae. We develop also truncated models of the original Gauss rules to compute integrals extended over the positive real axis. Numerical examples confirming the theoretical results are given comparing these rules among themselves and with different quadrature formulae proposed by other authors (Evans, Int. J. Comput. Math. 82:721?730, 2005; Gautschi, BIT 31:438?446, 1991).
The Use of Phase-Lag Derivatives in the Numerical Integration of ODEs with Oscillating Solutions
Anastassi, Z. A.; Vlachos, D. S.; Simos, T. E.
2008-09-01
In this paper we consider the fitting of the coefficients of a numerical method, not only due to the nullification of the phase-lag, but also to its derivatives. We show that the method gains efficiency with each derivative of the phase-lag nullified for various problems with oscillating solutions. The analysis of the local truncation error analysis and the stability of the methods show the importance of zero phase-lag derivatives when integrating oscillatory differential equations.
A leaky integrate-and-fire model with adaptation for the generation of a spike train.
Buonocore, Aniello; Caputo, Luigia; Pirozzi, Enrica; Carfora, Maria Francesca
2016-06-01
A model is proposed to describe the spike-frequency adaptation observed in many neuronal systems. We assume that adaptation is mainly due to a calcium-activated potassium current, and we consider two coupled stochastic differential equations for which an analytical approach combined with simulation techniques and numerical methods allow to obtain both qualitative and quantitative results about asymptotic mean firing rate, mean calcium concentration and the firing probability density. A related algorithm, based on the Hazard Rate Method, is also devised and described. PMID:27106179
Behrens, Bernd-Arno; Nolte, Ingo; Wefstaedt, Patrick; Stukenborg-Colsman, Christina; Bouguecha, Anas
2009-01-01
Background There are several numerical investigations on bone remodelling after total hip arthroplasty (THA) on the basis of the finite element analysis (FEA). For such computations certain boundary conditions have to be defined. The authors chose a maximum of three static load situations, usually taken from the gait cycle because this is the most frequent dynamic activity of a patient after THA. Materials and methods The numerical study presented here investigates whether it is useful to consider only one static load situation of the gait cycle in the FE calculation of the bone remodelling. For this purpose, 5 different loading cases were examined in order to determine their influence on the change in the physiological load distribution within the femur and on the resulting strain-adaptive bone remodelling. First, four different static loading cases at 25%, 45%, 65% and 85% of the gait cycle, respectively, and then the whole gait cycle in a loading regime were examined in order to regard all the different loadings of the cycle in the simulation. Results The computed evolution of the apparent bone density (ABD) and the calculated mass losses in the periprosthetic femur show that the simulation results are highly dependent on the chosen boundary conditions. Conclusion These numerical investigations prove that a static load situation is insufficient for representing the whole gait cycle. This causes severe deviations in the FE calculation of the bone remodelling. However, accompanying clinical examinations are necessary to calibrate the bone adaptation law and thus to validate the FE calculations. PMID:19371424
Adaptive learning of Multi-Sensor Integration techniques with genetic algorithms
Baker, J.E.
1994-06-01
This research focuses on automating the time-consuming process of developing and optimizing multi-sensor integration techniques. Our approach is currently based on adaptively learning how to exploit low-level image detail. Although this system is specifically designed to be both sensor and application domain independent, an empirical validation with actual multi-modal sensor data is presented.
Career Adaptability: An Integrative Construct for Life-Span, Life-Space Theory.
ERIC Educational Resources Information Center
Savickas, Mark L.
1997-01-01
Examines the origin and current status of lifespan, life-space theory and proposes one way in which to integrate its three segments. Discusses a functionalist strategy for theory construction and the outcomes and consequences of this strategy. Discusses future directions for theory development, such as career adaptability and planful attitudes.…
ERIC Educational Resources Information Center
International Migration, 1979
1979-01-01
This document contains working papers prepared for a seminar on Adaptation and Integration of Permanent Immigrants, along with general and specific recommendations formulated by seminar participants. Conclusions and recommendations from each paper are presented in English, French, and Spanish; the conference papers themselves are presented only in…
Integrated and adaptive management of water resources: Tensions, legacies, and the next best thing
Engle, Nathan L.; Johns, Owen R.; Lemos, Maria Carmen; Nelson, Donald
2011-02-01
Integrated water resources management (IWRM) and adaptive management (AM) are two institutional and management paradigms designed to address shortcomings within water systems governance – the limits of hierarchical water institutional arrangements in the case of IWRM and the challenge of making water management decisions under uncertainty in the case of AM. Recently, there has been a trend to merge these paradigms to address the growing complexity of stressors shaping water management, such as globalization and climate change. However, because many of these joint approaches have received little empirical attention, questions remain about how they might work (or not) in practice. Here, we explore a few of these issues using empirical research carried out in Brazil. We focus on highlighting the potentially negative interactions, tensions, and tradeoffs between different institutions/mechanisms perceived as desirable as research and practice attempt to make water systems management simultaneously integrated and adaptive. Our examples pertain mainly on the use of techno-scientific knowledge in water management and governance in Brazil’s IWRM model and how it relates to participation, democracy, deliberation, diversity, and adaptability. We show that a legacy of technical and hierarchical management has shaped the integration of management, and subsequently, the degree to which management might also be adaptive. While integrated systems may be more legitimate and accountable than top-down command and control ones, the mechanisms of IWRM may be at odds with the flexible, experimental, and self-organizing nature of AM.
NASA Technical Reports Server (NTRS)
Spratlin, Kenneth Milton
1987-01-01
An adaptive numeric predictor-corrector guidance is developed for atmospheric entry vehicles which utilize lift to achieve maximum footprint capability. Applicability of the guidance design to vehicles with a wide range of performance capabilities is desired so as to reduce the need for algorithm redesign with each new vehicle. Adaptability is desired to minimize mission-specific analysis and planning. The guidance algorithm motivation and design are presented. Performance is assessed for application of the algorithm to the NASA Entry Research Vehicle (ERV). The dispersions the guidance must be designed to handle are presented. The achievable operational footprint for expected worst-case dispersions is presented. The algorithm performs excellently for the expected dispersions and captures most of the achievable footprint.
Climate change adaptation and Integrated Water Resource Management in the water sector
NASA Astrophysics Data System (ADS)
Ludwig, Fulco; van Slobbe, Erik; Cofino, Wim
2014-10-01
Integrated Water Resources Management (IWRM) was introduced in 1980s to better optimise water uses between different water demanding sectors. However, since it was introduced water systems have become more complicated due to changes in the global water cycle as a result of climate change. The realization that climate change will have a significant impact on water availability and flood risks has driven research and policy making on adaptation. This paper discusses the main similarities and differences between climate change adaptation and IWRM. The main difference between the two is the focus on current and historic issues of IWRM compared to the (long-term) future focus of adaptation. One of the main problems of implementing climate change adaptation is the large uncertainties in future projections. Two completely different approaches to adaptation have been developed in response to these large uncertainties. A top-down approach based on large scale biophysical impacts analyses focussing on quantifying and minimizing uncertainty by using a large range of scenarios and different climate and impact models. The main problem with this approach is the propagation of uncertainties within the modelling chain. The opposite is the bottom up approach which basically ignores uncertainty. It focusses on reducing vulnerabilities, often at local scale, by developing resilient water systems. Both these approaches however are unsuitable for integrating into water management. The bottom up approach focuses too much on socio-economic vulnerability and too little on developing (technical) solutions. The top-down approach often results in an “explosion” of uncertainty and therefore complicates decision making. A more promising direction of adaptation would be a risk based approach. Future research should further develop and test an approach which starts with developing adaptation strategies based on current and future risks. These strategies should then be evaluated using a range
Analytical Solutions Using Integral Formulations and Their Coupling with Numerical Approaches.
Morel-Seytoux, Hubert J
2015-01-01
Analytical and numerical approaches have their own distinct domains of merit and application. Unfortunately there has been a tendency to use either one or the other even when their domains overlap. Yet there is definite advantage in combining the two approaches. Being relatively new this emerging technique of combining the approaches is, at this stage, more of an art than a science. In this article we suggest approaches for the combination through simple examples. We also suggest that the integral formulation of the analytical problems may have some advantages over the differential formulation. The differential formulation limits somewhat the range of linear system descriptions that can be applied to a variety of practical problems. On the other hand the integral approach tends to focus attention to overall integrated behavior and properties of the system rather than on minute details. This is particularly useful in the coupling with a numerical model as in practice it generally deals also with only the integrated behavior of the system. The thesis of this article is illustrated with some simple stream-aquifer flow exchange examples. PMID:25213772
DE 102 - A numerically integrated ephemeris of the moon and planets spanning forty-four centuries
NASA Technical Reports Server (NTRS)
Newhall, X. X.; Standish, E. M.; Willams, J. G.
1983-01-01
It is pointed out that the 1960's were the turning point for the generation of lunar and planetary ephemerides. All previous measurements of the positions of solar system bodies were optical angular measurements. New technological improvements leading to immense changes in observational accuracy are related to developments concerning radar, Viking landers on Mars, and laser ranges to lunar corner cube retroreflectors. Suitable numerical integration techniques and more comprehensive physical models were developed to match the accuracy of the modern data types. The present investigation is concerned with the first integrated ephemeris, DE 102, which covers the entire span of the historical astronomical observations of usable accuracy which are known. The fit is made to modern data. The integration spans the time period from 1411 BC to 3002 AD.
Valentín, A.; Humphrey, J. D.; Holzapfel, G. A.
2013-01-01
We implemented a constrained mixture model of arterial growth and remodeling (G&R) in a nonlinear finite element framework to facilitate numerical analyses of diverse cases of arterial adaptation and maladaptation, including disease progression, resulting in complex evolving geometries and compositions. This model enables hypothesis testing by predicting consequences of postulated characteristics of cell and matrix turnover, including evolving quantities and orientations of fibrillar constituents and non-homogenous degradation of elastin or loss of smooth muscle function. The non-linear finite element formulation is general within the context of arterial mechanics, but we restricted our present numerical verification to cylindrical geometries to allow comparisons to prior results for two special cases: uniform transmural changes in mass and differential G&R within a two-layered cylindrical model of the human aorta. The present finite element model recovers the results of these simplified semi-inverse analyses with good agreement. PMID:23713058
NASA Astrophysics Data System (ADS)
Grenga, Temistocle
The aim of this research is to further develop a dynamically adaptive algorithm based on wavelets that is able to solve efficiently multi-dimensional compressible reactive flow problems. This work demonstrates the great potential for the method to perform direct numerical simulation (DNS) of combustion with detailed chemistry and multi-component diffusion. In particular, it addresses the performance obtained using a massive parallel implementation and demonstrates important savings in memory storage and computational time over conventional methods. In addition, fully-resolved simulations of challenging three dimensional problems involving mixing and combustion processes are performed. These problems are particularly challenging due to their strong multiscale characteristics. For these solutions, it is necessary to combine the advanced numerical techniques applied to modern computational resources.
NASA Astrophysics Data System (ADS)
Macomber, B.; Woollands, R. M.; Probe, A.; Younes, A.; Bai, X.; Junkins, J.
2013-09-01
Modified Chebyshev Picard Iteration (MCPI) is an iterative numerical method for approximating solutions of linear or non-linear Ordinary Differential Equations (ODEs) to obtain time histories of system state trajectories. Unlike other step-by-step differential equation solvers, the Runge-Kutta family of numerical integrators for example, MCPI approximates long arcs of the state trajectory with an iterative path approximation approach, and is ideally suited to parallel computation. Orthogonal Chebyshev Polynomials are used as basis functions during each path iteration; the integrations of the Picard iteration are then done analytically. Due to the orthogonality of the Chebyshev basis functions, the least square approximations are computed without matrix inversion; the coefficients are computed robustly from discrete inner products. As a consequence of discrete sampling and weighting adopted for the inner product definition, Runge phenomena errors are minimized near the ends of the approximation intervals. The MCPI algorithm utilizes a vector-matrix framework for computational efficiency. Additionally, all Chebyshev coefficients and integrand function evaluations are independent, meaning they can be simultaneously computed in parallel for further decreased computational cost. Over an order of magnitude speedup from traditional methods is achieved in serial processing, and an additional order of magnitude is achievable in parallel architectures. This paper presents a new MCPI library, a modular toolset designed to allow MCPI to be easily applied to a wide variety of ODE systems. Library users will not have to concern themselves with the underlying mathematics behind the MCPI method. Inputs are the boundary conditions of the dynamical system, the integrand function governing system behavior, and the desired time interval of integration, and the output is a time history of the system states over the interval of interest. Examples from the field of astrodynamics are
ERIC Educational Resources Information Center
Raiche, Gilles; Blais, Jean-Guy
In a computerized adaptive test (CAT), it would be desirable to obtain an acceptable precision of the proficiency level estimate using an optimal number of items. Decreasing the number of items is accompanied, however, by a certain degree of bias when the true proficiency level differs significantly from the a priori estimate. G. Raiche (2000) has…
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
Quantitative evaluation of numerical integration schemes for Lagrangian particle dispersion models
NASA Astrophysics Data System (ADS)
Ramli, Huda Mohd.; Esler, J. Gavin
2016-07-01
A rigorous methodology for the evaluation of integration schemes for Lagrangian particle dispersion models (LPDMs) is presented. A series of one-dimensional test problems are introduced, for which the Fokker-Planck equation is solved numerically using a finite-difference discretisation in physical space and a Hermite function expansion in velocity space. Numerical convergence errors in the Fokker-Planck equation solutions are shown to be much less than the statistical error associated with a practical-sized ensemble (N = 106) of LPDM solutions; hence, the former can be used to validate the latter. The test problems are then used to evaluate commonly used LPDM integration schemes. The results allow for optimal time-step selection for each scheme, given a required level of accuracy. The following recommendations are made for use in operational models. First, if computational constraints require the use of moderate to long time steps, it is more accurate to solve the random displacement model approximation to the LPDM rather than use existing schemes designed for long time steps. Second, useful gains in numerical accuracy can be obtained, at moderate additional computational cost, by using the relatively simple "small-noise" scheme of Honeycutt.
Broom, Donald M
2006-01-01
The term adaptation is used in biology in three different ways. It may refer to changes which occur at the cell and organ level, or at the individual level, or at the level of gene action and evolutionary processes. Adaptation by cells, especially nerve cells helps in: communication within the body, the distinguishing of stimuli, the avoidance of overload and the conservation of energy. The time course and complexity of these mechanisms varies. Adaptive characters of organisms, including adaptive behaviours, increase fitness so this adaptation is evolutionary. The major part of this paper concerns adaptation by individuals and its relationships to welfare. In complex animals, feed forward control is widely used. Individuals predict problems and adapt by acting before the environmental effect is substantial. Much of adaptation involves brain control and animals have a set of needs, located in the brain and acting largely via motivational mechanisms, to regulate life. Needs may be for resources but are also for actions and stimuli which are part of the mechanism which has evolved to obtain the resources. Hence pigs do not just need food but need to be able to carry out actions like rooting in earth or manipulating materials which are part of foraging behaviour. The welfare of an individual is its state as regards its attempts to cope with its environment. This state includes various adaptive mechanisms including feelings and those which cope with disease. The part of welfare which is concerned with coping with pathology is health. Disease, which implies some significant effect of pathology, always results in poor welfare. Welfare varies over a range from very good, when adaptation is effective and there are feelings of pleasure or contentment, to very poor. A key point concerning the concept of individual adaptation in relation to welfare is that welfare may be good or poor while adaptation is occurring. Some adaptation is very easy and energetically cheap and
NASA Astrophysics Data System (ADS)
Vanderlinden, J. P.; Baztan, J.
2014-12-01
The prupose of this paper is to present the "Adaptation Research a Transdisciplinary community and policy centered appoach" (ARTisticc) project. ARTisticc's goal is to apply innovative standardized transdisciplinary art and science integrative approaches to foster robust, socially, culturally and scientifically, community centred adaptation to climate change. The approach used in the project is based on the strong understanding that adaptation is: (a) still "a concept of uncertain form"; (b) a concept dealing with uncertainty; (c) a concept that calls for an analysis that goes beyond the traditional disciplinary organization of science, and; (d) an unconventional process in the realm of science and policy integration. The project is centered on case studies in France, Greenland, Russia, India, Canada, Alaska, and Senegal. In every site we jointly develop artwork while we analyzing how natural science, essentially geosciences can be used in order to better adapt in the future, how society adapt to current changes and how memories of past adaptations frames current and future processes. Artforms are mobilized in order to share scientific results with local communities and policy makers, this in a way that respects cultural specificities while empowering stakeholders, ARTISTICC translates these "real life experiments" into stories and artwork that are meaningful to those affected by climate change. The scientific results and the culturally mediated productions will thereafter be used in order to co-construct, with NGOs and policy makers, policy briefs, i.e. robust and scientifically legitimate policy recommendations regarding coastal adaptation. This co-construction process will be in itself analysed with the goal of increasing arts and science's performative functions in the universe of evidence-based policy making. The project involves scientists from natural sciences, the social sciences and the humanities, as well as artitis from the performing arts (playwriters
Rivera, Claudia
2014-01-01
This paper analyses the perceptions of disaster risk reduction (DRR) practitioners concerning the on-going integration of climate change adaptation (CCA) into their practices in urban contexts in Nicaragua. Understanding their perceptions is important as this will provide information on how this integration can be improved. Exploring the perceptions of practitioners in Nicaragua is important as the country has a long history of disasters, and practitioners have been developing the current DRR planning framework for more than a decade. The analysis is based on semi-structured interviews designed to collect information about practitioners’ understanding of: (a) CCA, (b) the current level of integration of CCA into DRR and urban planning, (c) the opportunities and constraints of this integration, and (d) the potential to adapt cities to climate change. The results revealed that practitioners’ perception is that the integration of CCA into their practice is at an early stage, and that they need to improve their understanding of CCA in terms of a development issue. Three main constraints on improved integration were identified: (a) a recognized lack of understanding of CCA, (b) insufficient guidance on how to integrate it, and (c) the limited opportunities to integrate it into urban planning due to a lack of instruments and capacity in this field. Three opportunities were also identified: (a) practitioners’ awareness of the need to integrate CCA into their practices, (b) the robust structure of the DRR planning framework in the country, which provides a suitable channel for facilitating integration, and (c) the fact that CCA is receiving more attention and financial and technical support from the international community. PMID:24475365
Modeling, mesh generation, and adaptive numerical methods for partial differential equations
Babuska, I.; Henshaw, W.D.; Oliger, J.E.; Flaherty, J.E.; Hopcroft, J.E.; Tezduyar, T.
1995-12-31
Mesh generation is one of the most time consuming aspects of computational solutions of problems involving partial differential equations. It is, furthermore, no longer acceptable to compute solutions without proper verification that specified accuracy criteria are being satisfied. Mesh generation must be related to the solution through computable estimates of discretization errors. Thus, an iterative process of alternate mesh and solution generation evolves in an adaptive manner with the end result that the solution is computed to prescribed specifications in an optimal, or at least efficient, manner. While mesh generation and adaptive strategies are becoming available, major computational challenges remain. One, in particular, involves moving boundaries and interfaces, such as free-surface flows and fluid-structure interactions. A 3-week program was held from July 5 to July 23, 1993 with 173 participants and 66 keynote, invited, and contributed presentations. This volume represents written versions of 21 of these lectures. These proceedings are organized roughly in order of their presentation at the workshop. Thus, the initial papers are concerned with geometry and mesh generation and discuss the representation of physical objects and surfaces on a computer and techniques to use this data to generate, principally, unstructured meshes of tetrahedral or hexahedral elements. The remainder of the papers cover adaptive strategies, error estimation, and applications. Several submissions deal with high-order p- and hp-refinement methods where mesh refinement/coarsening (h-refinement) is combined with local variation of method order (p-refinement). Combinations of mathematically verified and physically motivated approaches to error estimation are represented. Applications center on fluid mechanics. Selected papers are indexed separately for inclusion in the Energy Science and Technology Database.
NASA Technical Reports Server (NTRS)
Hu, Fang Q.
1994-01-01
It is known that the exact analytic solutions of wave scattering by a circular cylinder, when they exist, are not in a closed form but in infinite series which converges slowly for high frequency waves. In this paper, we present a fast number solution for the scattering problem in which the boundary integral equations, reformulated from the Helmholtz equation, are solved using a Fourier spectral method. It is shown that the special geometry considered here allows the implementation of the spectral method to be simple and very efficient. The present method differs from previous approaches in that the singularities of the integral kernels are removed and dealt with accurately. The proposed method preserves the spectral accuracy and is shown to have an exponential rate of convergence. Aspects of efficient implementation using FFT are discussed. Moreover, the boundary integral equations of combined single and double-layer representation are used in the present paper. This ensures the uniqueness of the numerical solution for the scattering problem at all frequencies. Although a strongly singular kernel is encountered for the Neumann boundary conditions, we show that the hypersingularity can be handled easily in the spectral method. Numerical examples that demonstrate the validity of the method are also presented.
NASA Technical Reports Server (NTRS)
Brislawn, Kristi D.; Brown, David L.; Chesshire, Geoffrey S.; Saltzman, Jeffrey S.
1995-01-01
Adaptive mesh refinement (AMR) in conjunction with higher-order upwind finite-difference methods have been used effectively on a variety of problems in two and three dimensions. In this paper we introduce an approach for resolving problems that involve complex geometries in which resolution of boundary geometry is important. The complex geometry is represented by using the method of overlapping grids, while local resolution is obtained by refining each component grid with the AMR algorithm, appropriately generalized for this situation. The CMPGRD algorithm introduced by Chesshire and Henshaw is used to automatically generate the overlapping grid structure for the underlying mesh.
Numerical simulation of strain-adaptive bone remodelling in the ankle joint
2011-01-01
Background The use of artificial endoprostheses has become a routine procedure for knee and hip joints while ankle arthritis has traditionally been treated by means of arthrodesis. Due to its advantages, the implantation of endoprostheses is constantly increasing. While finite element analyses (FEA) of strain-adaptive bone remodelling have been carried out for the hip joint in previous studies, to our knowledge there are no investigations that have considered remodelling processes of the ankle joint. In order to evaluate and optimise new generation implants of the ankle joint, as well as to gain additional knowledge regarding the biomechanics, strain-adaptive bone remodelling has been calculated separately for the tibia and the talus after providing them with an implant. Methods FE models of the bone-implant assembly for both the tibia and the talus have been developed. Bone characteristics such as the density distribution have been applied corresponding to CT scans. A force of 5,200 N, which corresponds to the compression force during normal walking of a person with a weight of 100 kg according to Stauffer et al., has been used in the simulation. The bone adaptation law, previously developed by our research team, has been used for the calculation of the remodelling processes. Results A total bone mass loss of 2% in the tibia and 13% in the talus was calculated. The greater decline of density in the talus is due to its smaller size compared to the relatively large implant dimensions causing remodelling processes in the whole bone tissue. In the tibia, bone remodelling processes are only calculated in areas adjacent to the implant. Thus, a smaller bone mass loss than in the talus can be expected. There is a high agreement between the simulation results in the distal tibia and the literature regarding. Conclusions In this study, strain-adaptive bone remodelling processes are simulated using the FE method. The results contribute to a better understanding of the
Numerical integration techniques for curved-element discretizations of molecule-solvent interfaces.
Bardhan, Jaydeep P; Altman, Michael D; Willis, David J; Lippow, Shaun M; Tidor, Bruce; White, Jacob K
2007-07-01
Surface formulations of biophysical modeling problems offer attractive theoretical and computational properties. Numerical simulations based on these formulations usually begin with discretization of the surface under consideration; often, the surface is curved, possessing complicated structure and possibly singularities. Numerical simulations commonly are based on approximate, rather than exact, discretizations of these surfaces. To assess the strength of the dependence of simulation accuracy on the fidelity of surface representation, here methods were developed to model several important surface formulations using exact surface discretizations. Following and refining Zauhar's work [J. Comput.-Aided Mol. Des. 9, 149 (1995)], two classes of curved elements were defined that can exactly discretize the van der Waals, solvent-accessible, and solvent-excluded (molecular) surfaces. Numerical integration techniques are presented that can accurately evaluate nonsingular and singular integrals over these curved surfaces. After validating the exactness of the surface discretizations and demonstrating the correctness of the presented integration methods, a set of calculations are presented that compare the accuracy of approximate, planar-triangle-based discretizations and exact, curved-element-based simulations of surface-generalized-Born (sGB), surface-continuum van der Waals (scvdW), and boundary-element method (BEM) electrostatics problems. Results demonstrate that continuum electrostatic calculations with BEM using curved elements, piecewise-constant basis functions, and centroid collocation are nearly ten times more accurate than planar-triangle BEM for basis sets of comparable size. The sGB and scvdW calculations give exceptional accuracy even for the coarsest obtainable discretized surfaces. The extra accuracy is attributed to the exact representation of the solute-solvent interface; in contrast, commonly used planar-triangle discretizations can only offer improved
CALL FOR PAPERS: Special Issue on `Geometric Numerical Integration of Differential Equations'
NASA Astrophysics Data System (ADS)
Quispel, G. R. W.; McLachlan, R. I.
2005-02-01
This is a call for contributions to a special issue of Journal of Physics A: Mathematical and General entitled `Geometric Numerical Integration of Differential Equations'. This issue should be a repository for high quality original work. We are interested in having the topic interpreted broadly, that is, to include contributions dealing with symplectic or multisymplectic integration; volume-preserving integration; symmetry-preserving integration; integrators that preserve first integrals, Lyapunov functions, or dissipation; exponential integrators; integrators for highly oscillatory systems; Lie-group integrators, etc. Papers on geometric integration of both ODEs and PDEs will be considered, as well as application to molecular-scale integration, celestial mechanics, particle accelerators, fluid flows, population models, epidemiological models and/or any other areas of science. We believe that this issue is timely, and hope that it will stimulate further development of this new and exciting field. The Editorial Board has invited G R W Quispel and R I McLachlan to serve as Guest Editors for the special issue. Their criteria for acceptance of contributions are the following: • The subject of the paper should relate to geometric numerical integration in the sense described above. • Contributions will be refereed and processed according to the usual procedure of the journal. • Papers should be original; reviews of a work published elsewhere will not be accepted. The guidelines for the preparation of contributions are as follows: • The DEADLINE for submission of contributions is 1 September 2005. This deadline will allow the special issue to appear in late 2005 or early 2006. • There is a strict page limit of 16 printed pages (approximately 9600 words) per contribution. For papers exceeding this limit, the Guest Editors reserve the right to request a reduction in length. Further advice on publishing your work in Journal of Physics A: Mathematical and General
Sensitivity of inelastic response to numerical integration of strain energy. [for cantilever beam
NASA Technical Reports Server (NTRS)
Kamat, M. P.
1976-01-01
The exact solution to the quasi-static, inelastic response of a cantilever beam of rectangular cross section subjected to a bending moment at the tip is obtained. The material of the beam is assumed to be linearly elastic-linearly strain-hardening. This solution is then compared with three different numerical solutions of the same problem obtained by minimizing the total potential energy using Gaussian quadratures of two different orders and a Newton-Cotes scheme for integrating the strain energy of deformation. Significant differences between the exact dissipative strain energy and its numerical counterpart are emphasized. The consequence of this on the nonlinear transient responses of a beam with solid cross section and that of a thin-walled beam on elastic supports under impulsive loads are examined.
Does Integration Help Adapt to Climate Change? Case of Increased US Corn Yield Volatility
NASA Astrophysics Data System (ADS)
Verma, M.; Diffenbaugh, N. S.; Hertel, T. W.
2012-12-01
In absence of of new crop varieties or significant shifts in the geography of corn production, US national corn yields variation could double by the year 2040 as a result of climate change and without adaptation this could lead the variability in US corn prices to quadruple (Diffenbaugh et al. 2012). In addition to climate induced price changes, analysis of recent commodity price spikes suggests that interventionist trade policies are partly to blame. Assuming we cannot much influence the future climate outcome, what policies can we undertake to adapt better? Can we use markets to blunt this edge? Diffenbaugh et al. find that sale of corn- ethanol for use in liquid fuel, when governed by quotas such as US Renewable Fuel Standard (RFS), could make US corn prices even more variable; in contrast the same food-fuel market link (we refer to it as intersectoral link) may well dampen price volatility when the sale of corn to ethanol industry is driven by higher future oil prices. The latter however comes at the cost of exposing corn prices to the greater volatility in oil markets. Similarly intervention in corn trade can make US corn prices less or more volatile by distorting international corn price transmission. A negative US corn yield shock shows that domestic corn supply falls and domestic prices to go up irrespective of whether or not markets are integrated. How much the prices go up depends on how much demand adjusts to accommodate the supply shock. Based on the forgoing analysis, one should expect that demand would adjust more readily when markets are integrated and therefore reduce the resulting price fluctuation. Simulation results confirm this response of corn markets. In terms of relative comparisons however a policy driven intersectoral integration is least effective and prices rise much more. Similarly, a positive world oil price shock makes the US oil imports expensive and with oil being used to produce gasoline blends, it increases the price of gasoline
NASA Astrophysics Data System (ADS)
Carbillet, Marcel; Riccardi, Armando; Esposito, Simone
2004-10-01
We present our latest results concerning the simulation studies performed for the first-light adaptive optics (AO) system of the Large Binocular Telescope (LBT), namely WLBT. After a brief description of the "raw" performance evaluation results, in terms of Strehl ratios attained in the various considered bands (from V to K), we focus on the "scientific" performance that will be obtained when considering the subsequent instrumentation that will benefit from the correction given by the AO system WLBT and the adaptive secondary mirrors LBT 672. In particular, we discuss the performance of the coupling with the instrument LUCIFER, working at near-infrared bands, in terms of signal-to-noise values and limiting magnitudes, and in both the cases of spectroscopy and photometric detection. We also give the encircled energies that are expected in the visible bands, result relevant in one hand for the instrument PEPSI, and in other hand for the "technical viewer" that will be on board the WLBT system itself.
Numerical and analytical tests of quasi-integrability in modified sine-Gordon models
NASA Astrophysics Data System (ADS)
Ferreira, L. A.; Zakrzewski, Wojtek J.
2014-01-01
Following our attempts to define quasi-integrability in which we related this concept to a particular symmetry of the two-soliton function we check this condition in three classes of modified sine-Gordon models in (1 + 1) dimensions. We find that the numerical results seen in various scatterings of two solitons and in the time evolution of breather-like structures support our ideas about the symmetry of the field configurations and its effects on the anomalies of the conservation laws of the charges.
Time transformations and Cowell's method. [for numerical integration of satellite motion equations
NASA Technical Reports Server (NTRS)
Velez, C. E.; Hilinski, S.
1978-01-01
The precise numerical integration of Cowell's equations of satellite motion is frequently performed with an independent variable s defined by an equation of the form dt = cr to the n-th power ds, where t represents time, r the radial distance from the center of attraction, c is a constant, and n is a parameter. This has been primarily motivated by the 'uniformizing' effects of such a transformation resulting in desirable 'analytic' stepsize control for elliptical orbits. This report discusses the 'proper' choice of the parameter n defining the independent variable s for various types of orbits and perturbation models, and develops a criterion for its selection.
Suárez-Pellicioni, Macarena; Núñez-Peña, María Isabel; Colomé, Àngels
2014-01-01
This study uses event-related brain potentials (ERPs) to investigate the electrophysiological correlates of numeric conflict monitoring in math-anxious individuals, by analyzing whether math anxiety is related to abnormal processing in early conflict detection (as shown by the N450 component) and/or in a later, response-related stage of processing (as shown by the conflict sustained potential; Conflict-SP). Conflict adaptation effects were also studied by analyzing the effect of the previous trial’s congruence in current interference. To this end, 17 low math-anxious (LMA) and 17 high math-anxious (HMA) individuals were presented with a numerical Stroop task. Groups were extreme in math anxiety but did not differ in trait or state anxiety or in simple math ability. The interference effect of the current trial (incongruent-congruent) and the interference effect preceded by congruence and by incongruity were analyzed both for behavioral measures and for ERPs. A greater interference effect was found for response times in the HMA group than in the LMA one. Regarding ERPs, the LMA group showed a greater N450 component for the interference effect preceded by congruence than when preceded by incongruity, while the HMA group showed greater Conflict-SP amplitude for the interference effect preceded by congruence than when preceded by incongruity. Our study showed that the electrophysiological correlates of numeric interference in HMA individuals comprise the absence of a conflict adaptation effect in the first stage of conflict processing (N450) and an abnormal subsequent up-regulation of cognitive control in order to overcome the conflict (Conflict-SP). More concretely, our study shows that math anxiety is related to a reactive and compensatory recruitment of control resources that is implemented only when previously exposed to a stimuli presenting conflicting information. PMID:24918584
Integration of adaptive process control with computational simulation for spin-forming
Raboin, P. J., LLNL
1998-03-10
Improvements in spin-forming capabilities through upgrades to a metrology and machine control system and advances in numerical simulation techniques were studied in a two year project funded by Laboratory Directed Research and Development (LDRD) at Lawrence Livermore National Laboratory. Numerical analyses were benchmarked with spin-forming experiments and computational speeds increased sufficiently to now permit actual part forming simulations. Extensive modeling activities examined the simulation speeds and capabilities of several metal forming computer codes for modeling flat plate and cylindrical spin-forming geometries. Shape memory research created the first numerical model to describe this highly unusual deformation behavior in Uranium alloys. A spin-forming metrology assessment led to sensor and data acquisition improvements that will facilitate future process accuracy enhancements, such as a metrology frame. Finally, software improvements (SmartCAM) to the manufacturing process numerically integrate the part models to the spin-forming process and to computational simulations.
NASA Technical Reports Server (NTRS)
Majumdar, Alok K.; LeClair, Andre C.; Hedayat, Ali
2016-01-01
This paper presents a numerical model of pressurization of a cryogenic propellant tank for the Integrated Vehicle Fluid (IVF) system using the Generalized Fluid System Simulation Program (GFSSP). The IVF propulsion system, being developed by United Launch Alliance, uses boiloff propellants to drive thrusters for the reaction control system as well as to run internal combustion engines to develop power and drive compressors to pressurize propellant tanks. NASA Marshall Space Flight Center (MSFC) has been running tests to verify the functioning of the IVF system using a flight tank. GFSSP, a finite volume based flow network analysis software developed at MSFC, has been used to develop an integrated model of the tank and the pressurization system. This paper presents an iterative algorithm for converging the interface boundary conditions between different component models of a large system model. The model results have been compared with test data.
Three-dimensional numerical modeling of photonic integration with dielectric-loaded SPP waveguides
NASA Astrophysics Data System (ADS)
Krasavin, A. V.; Zayats, A. V.
2008-07-01
Using full three-dimensional numerical modeling, we demonstrate highly efficient passive and active photonic circuit elements based on dielectric-loaded surface plasmon polariton waveguides (DLSPPWs). Highly confined surface plasmon polariton (SPP) mode having subwavelength cross section allows high level of integration of DLSPPW circuitry. We demonstrate very efficient guiding and routing of SPP signals with the passive waveguide elements such as bends, splitters, and Bragg reflectors, having a functional size of just a few microns at telecommunication wavelengths. Introducing a gain in the dielectric, we have found the requirement for lossless waveguiding and estimated the performance of DLSPPW lossless and active elements. DLSPPW based components have prospective implementation in photonic integrated chips, hybrid optical-electronic circuits, and lab-on-a-chip applications.
NASA Astrophysics Data System (ADS)
Zanotti, O.; Dumbser, M.; Fambri, F.
2016-05-01
We describe a new method for the solution of the ideal MHD equations in special relativity which adopts the following strategy: (i) the main scheme is based on Discontinuous Galerkin (DG) methods, allowing for an arbitrary accuracy of order N+1, where N is the degree of the basis polynomials; (ii) in order to cope with oscillations at discontinuities, an ”a-posteriori” sub-cell limiter is activated, which scatters the DG polynomials of the previous time-step onto a set of 2N+1 sub-cells, over which the solution is recomputed by means of a robust finite volume scheme; (iii) a local spacetime Discontinuous-Galerkin predictor is applied both on the main grid of the DG scheme and on the sub-grid of the finite volume scheme; (iv) adaptive mesh refinement (AMR) with local time-stepping is used. We validate the new scheme and comment on its potential applications in high energy astrophysics.
NASA Astrophysics Data System (ADS)
Kuraz, Michal
2016-06-01
This paper presents pseudo-deterministic catchment runoff model based on the Richards equation model [1] - the governing equation for the subsurface flow. The subsurface flow in a catchment is described here by two-dimensional variably saturated flow (unsaturated and saturated). The governing equation is the Richards equation with a slight modification of the time derivative term as considered e.g. by Neuman [2]. The nonlinear nature of this problem appears in unsaturated zone only, however the delineation of the saturated zone boundary is a nonlinear computationally expensive issue. The simple one-dimensional Boussinesq equation was used here as a rough estimator of the saturated zone boundary. With this estimate the dd-adaptivity algorithm (see Kuraz et al. [4, 5, 6]) could always start with an optimal subdomain split, so it is now possible to avoid solutions of huge systems of linear equations in the initial iteration level of our Richards equation based runoff model.
Antithetic Integral Feedback Ensures Robust Perfect Adaptation in Noisy Biomolecular Networks.
Briat, Corentin; Gupta, Ankit; Khammash, Mustafa
2016-01-27
The ability to adapt to stimuli is a defining feature of many biological systems and critical to maintaining homeostasis. While it is well appreciated that negative feedback can be used to achieve homeostasis when networks behave deterministically, the effect of noise on their regulatory function is not understood. Here, we combine probability and control theory to develop a theory of biological regulation that explicitly takes into account the noisy nature of biochemical reactions. We introduce tools for the analysis and design of robust homeostatic circuits and propose a new regulation motif, which we call antithetic integral feedback. This motif exploits stochastic noise, allowing it to achieve precise regulation in scenarios where similar deterministic regulation fails. Specifically, antithetic integral feedback preserves the stability of the overall network, steers the population of any regulated species to a desired set point, and adapts perfectly. We suggest that this motif may be prevalent in endogenous biological circuits and useful when creating synthetic circuits. PMID:27136686
NASA Technical Reports Server (NTRS)
Baer-Riedhart, Jennifer L.; Landy, Robert J.
1987-01-01
The highly integrated digital electronic control (HIDEC) program at NASA Ames Research Center, Dryden Flight Research Facility is a multiphase flight research program to quantify the benefits of promising integrated control systems. McDonnell Aircraft Company is the prime contractor, with United Technologies Pratt and Whitney Aircraft, and Lear Siegler Incorporated as major subcontractors. The NASA F-15A testbed aircraft was modified by the HIDEC program by installing a digital electronic flight control system (DEFCS) and replacing the standard F100 (Arab 3) engines with F100 engine model derivative (EMD) engines equipped with digital electronic engine controls (DEEC), and integrating the DEEC's and DEFCS. The modified aircraft provides the capability for testing many integrated control modes involving the flight controls, engine controls, and inlet controls. This paper focuses on the first two phases of the HIDEC program, which are the digital flight control system/aircraft model identification (DEFCS/AMI) phase and the adaptive engine control system (ADECS) phase.
Assessment of Disaster Risk Reduction and Climate Change Adaptation policy integration in Zambia
NASA Astrophysics Data System (ADS)
Pilli-Sihvola, K.; Väätäinen-Chimpuku, S.
2015-12-01
Integration of Disaster Risk Management (DRM) and Climate Change Adaptation (CCA) policies, their implementation measures and the contribution of these to development has been gaining attention recently. Due to the shared objectives of CCA and particularly Disaster Risk Reduction (DRR), a component of DRM, their integration provides many benefits. At the implementation level, DRR and CCA are usually integrated. Policy integration, however, is often lacking. This study presents a novel analysis of the policy integration of DRR and CCA by 1) suggesting a definition for their integration at a general and further at horizontal and vertical levels, 2) using an analysis framework for policy integration cycle, which separates the policy formulation and implementation processes, and 3) applying these to a case study in Zambia. Moreover, the study identifies the key gaps in the integration process, obtains an understanding of identified key factors for creating an enabling environment for the integration, and provides recommendations for further progress. The study is based on a document analysis of the relevant DRM, climate change (CC), agriculture, forestry, water management and meteorology policy documents and Acts, and 21 semi-structured interviews with key stakeholders. Horizontal integration has occurred both ways, as the revised DRM policy draft has incorporated CCA, and the new CC policy draft has incorporated DRR. This is not necessarily an optimal strategy and unless carefully implemented, it may create pressure on institutional structures and duplication of efforts in the implementation. Much less vertical integration takes place, and where it does, no guidance on how potential goal conflicts with sectorial and development objectives ought to be handled. The objectives of the instruments show convergence. At the programme stage, the measures are fully integrated as they can be classified as robust CCA measures, providing benefits in the current and future
Owolabi, Kolade M; Patidar, Kailash C
2016-01-01
In this paper, we consider the numerical simulations of an extended nonlinear form of Kierstead-Slobodkin reaction-transport system in one and two dimensions. We employ the popular fourth-order exponential time differencing Runge-Kutta (ETDRK4) schemes proposed by Cox and Matthew (J Comput Phys 176:430-455, 2002), that was modified by Kassam and Trefethen (SIAM J Sci Comput 26:1214-1233, 2005), for the time integration of spatially discretized partial differential equations. We demonstrate the supremacy of ETDRK4 over the existing exponential time differencing integrators that are of standard approaches and provide timings and error comparison. Numerical results obtained in this paper have granted further insight to the question 'What is the minimal size of the spatial domain so that the population persists?' posed by Kierstead and Slobodkin (J Mar Res 12:141-147, 1953), with a conclusive remark that the population size increases with the size of the domain. In attempt to examine the biological wave phenomena of the solutions, we present the numerical results in both one- and two-dimensional space, which have interesting ecological implications. Initial data and parameter values were chosen to mimic some existing patterns. PMID:27064984
NASA Astrophysics Data System (ADS)
Gong, He; Fan, Yubo; Zhang, Ming
2008-04-01
The objective of this paper is to identify the effects of mechanical disuse and basic multi-cellular unit (BMU) activation threshold on the form of trabecular bone during menopause. A bone adaptation model with mechanical- biological factors at BMU level was integrated with finite element analysis to simulate the changes of trabecular bone structure during menopause. Mechanical disuse and changes in the BMU activation threshold were applied to the model for the period from 4 years before to 4 years after menopause. The changes in bone volume fraction, trabecular thickness and fractal dimension of the trabecular structures were used to quantify the changes of trabecular bone in three different cases associated with mechanical disuse and BMU activation threshold. It was found that the changes in the simulated bone volume fraction were highly correlated and consistent with clinical data, and that the trabecular thickness reduced significantly during menopause and was highly linearly correlated with the bone volume fraction, and that the change trend of fractal dimension of the simulated trabecular structure was in correspondence with clinical observations. The numerical simulation in this paper may help to better understand the relationship between the bone morphology and the mechanical, as well as biological environment; and can provide a quantitative computational model and methodology for the numerical simulation of the bone structural morphological changes caused by the mechanical environment, and/or the biological environment.
NASA Technical Reports Server (NTRS)
Radhakrishnan, K.
1984-01-01
The efficiency and accuracy of several algorithms recently developed for the efficient numerical integration of stiff ordinary differential equations are compared. The methods examined include two general-purpose codes, EPISODE and LSODE, and three codes (CHEMEQ, CREK1D, and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes are applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code currently available for the integration of combustion kinetic rate equations. An important finding is that an interactive solution of the algebraic energy conservation equation to compute the temperature does not result in significant errors. In addition, this method is more efficient than evaluating the temperature by integrating its time derivative. Significant reductions in computational work are realized by updating the rate constants (k = at(supra N) N exp(-E/RT) only when the temperature change exceeds an amount delta T that is problem dependent. An approximate expression for the automatic evaluation of delta T is derived and is shown to result in increased efficiency.
Mie Light-Scattering Granulometer with an Adaptive Numerical Filtering Method. II. Experiments.
Hespel, L; Delfour, A; Guillame, B
2001-02-20
A nephelometer is presented that theoretically requires no absolute calibration. This instrument is used for determining the particle-size distribution of various scattering media (aerosols, fogs, rocket exhausts, engine plumes, and the like) from angular static light-scattering measurements. An inverse procedure is used, which consists of a least-squares method and a regularization scheme based on numerical filtering. To retrieve the distribution function one matches the experimental data with theoretical patterns derived from Mie theory. The main principles of the inverse method are briefly presented, and the nephelometer is then described with the associated partial calibration procedure. Finally, the whole granulometer system (inverse method and nephelometer) is validated by comparison of measurements of scattering media with calibrated monodisperse or known size distribution functions. PMID:18357082
Zhou, Qianqian; Panduro, Toke Emil; Thorsen, Bo Jellesmark; Arnbjerg-Nielsen, Karsten
2013-03-01
This paper presents a cross-disciplinary framework for assessment of climate change adaptation to increased precipitation extremes considering pluvial flood risk as well as additional environmental services provided by some of the adaptation options. The ability of adaptation alternatives to cope with extreme rainfalls is evaluated using a quantitative flood risk approach based on urban inundation modeling and socio-economic analysis of corresponding costs and benefits. A hedonic valuation model is applied to capture the local economic gains or losses from more water bodies in green areas. The framework was applied to the northern part of the city of Aarhus, Denmark. We investigated four adaptation strategies that encompassed laissez-faire, larger sewer pipes, local infiltration units, and open drainage system in the urban green structure. We found that when taking into account environmental amenity effects, an integration of open drainage basins in urban recreational areas is likely the best adaptation strategy, followed by pipe enlargement and local infiltration strategies. All three were improvements compared to the fourth strategy of no measures taken. PMID:23334752
NASA Astrophysics Data System (ADS)
Zhou, Qianqian; Panduro, Toke Emil; Thorsen, Bo Jellesmark; Arnbjerg-Nielsen, Karsten
2013-03-01
This paper presents a cross-disciplinary framework for assessment of climate change adaptation to increased precipitation extremes considering pluvial flood risk as well as additional environmental services provided by some of the adaptation options. The ability of adaptation alternatives to cope with extreme rainfalls is evaluated using a quantitative flood risk approach based on urban inundation modeling and socio-economic analysis of corresponding costs and benefits. A hedonic valuation model is applied to capture the local economic gains or losses from more water bodies in green areas. The framework was applied to the northern part of the city of Aarhus, Denmark. We investigated four adaptation strategies that encompassed laissez-faire, larger sewer pipes, local infiltration units, and open drainage system in the urban green structure. We found that when taking into account environmental amenity effects, an integration of open drainage basins in urban recreational areas is likely the best adaptation strategy, followed by pipe enlargement and local infiltration strategies. All three were improvements compared to the fourth strategy of no measures taken.
NASA Astrophysics Data System (ADS)
Banaś, Krzysztof; Krużel, Filip; Bielański, Jan
2016-06-01
The paper presents investigations on the implementation and performance of the finite element numerical integration algorithm for first order approximations and three processor architectures, popular in scientific computing, classical CPU, Intel Xeon Phi and NVIDIA Kepler GPU. A unifying programming model and portable OpenCL implementation is considered for all architectures. Variations of the algorithm due to different problems solved and different element types are investigated and several optimizations aimed at proper optimization and mapping of the algorithm to computer architectures are demonstrated. Performance models of execution are developed for different processors and tested in practical experiments. The results show the varying levels of performance for different architectures, but indicate that the algorithm can be effectively ported to all of them. The general conclusion is that the finite element numerical integration can achieve sufficient performance on different multi- and many-core architectures and should not become a performance bottleneck for finite element simulation codes. Specific observations lead to practical advises on how to optimize the kernels and what performance can be expected for the tested architectures.
NASA Astrophysics Data System (ADS)
Ding, Ye; Zhu, Limin; Zhang, Xiaojian; Ding, Han
2012-09-01
As one of the bases of gradient-based optimization algorithms, sensitivity analysis is usually required to calculate the derivatives of the system response with respect to the machining parameters. The most widely used approaches for sensitivity analysis are based on time-consuming numerical methods, such as finite difference methods. This paper presents a semi-analytical method for calculation of the sensitivity of the stability boundary in milling. After transforming the delay-differential equation with time-periodic coefficients governing the dynamic milling process into the integral form, the Floquet transition matrix is constructed by using the numerical integration method. Then, the analytical expressions of derivatives of the Floquet transition matrix with respect to the machining parameters are obtained. Thereafter, the classical analytical expression of the sensitivity of matrix eigenvalues is employed to calculate the sensitivity of the stability lobe diagram. The two-degree-of-freedom milling example illustrates the accuracy and efficiency of the proposed method. Compared with the existing methods, the unique merit of the proposed method is that it can be used for analytically computing the sensitivity of the stability boundary in milling, without employing any finite difference methods. Therefore, the high accuracy and high efficiency are both achieved. The proposed method can serve as an effective tool for machining parameter optimization and uncertainty analysis in high-speed milling.
Application of Numerical Integration and Data Fusion in Unit Vector Method
NASA Astrophysics Data System (ADS)
Zhang, J.
2012-01-01
The Unit Vector Method (UVM) is a series of orbit determination methods which are designed by Purple Mountain Observatory (PMO) and have been applied extensively. It gets the conditional equations for different kinds of data by projecting the basic equation to different unit vectors, and it suits for weighted process for different kinds of data. The high-precision data can play a major role in orbit determination, and accuracy of orbit determination is improved obviously. The improved UVM (PUVM2) promoted the UVM from initial orbit determination to orbit improvement, and unified the initial orbit determination and orbit improvement dynamically. The precision and efficiency are improved further. In this thesis, further research work has been done based on the UVM: Firstly, for the improvement of methods and techniques for observation, the types and decision of the observational data are improved substantially, it is also asked to improve the decision of orbit determination. The analytical perturbation can not meet the requirement. So, the numerical integration for calculating the perturbation has been introduced into the UVM. The accuracy of dynamical model suits for the accuracy of the real data, and the condition equations of UVM are modified accordingly. The accuracy of orbit determination is improved further. Secondly, data fusion method has been introduced into the UVM. The convergence mechanism and the defect of weighted strategy have been made clear in original UVM. The problem has been solved in this method, the calculation of approximate state transition matrix is simplified and the weighted strategy has been improved for the data with different dimension and different precision. Results of orbit determination of simulation and real data show that the work of this thesis is effective: (1) After the numerical integration has been introduced into the UVM, the accuracy of orbit determination is improved obviously, and it suits for the high-accuracy data of
Integrated Numerical Experiments (INEX) and the Free-Electron Laser Physical Process Code (FELPPC)
Thode, L.E.; Chan, K.C.D.; Schmitt, M.J.; McKee, J.; Ostic, J.; Elliott, C.J.; McVey, B.D.
1990-01-01
The strong coupling of subsystem elements, such as the accelerator, wiggler, and optics, greatly complicates the understanding and design of a free electron laser (FEL), even at the conceptual level. Given the requirements for high-performance FELs, the strong coupling between the laser subsystems must be included to obtain a realistic picture of the potential operational capability. To address the strong coupling character of the FEL the concept of an Integrated Numerical Experiment (INEX) was proposed. Unique features of the INEX approach are consistency and numerical equivalence of experimental diagnostics. The equivalent numerical diagnostics mitigates the major problem of misinterpretation that often occurs when theoretical and experimental data are compared. The INEX approach has been applied to a large number of accelerator and FEL experiments. Overall, the agreement between INEX and the experiments is very good. Despite the success of INEX, the approach is difficult to apply to trade-off and initial design studies because of the significant manpower and computational requirements. On the other hand, INEX provides a base from which realistic accelerator, wiggler, and optics models can be developed. The Free Electron Laser Physical Process Code (FELPPC) includes models developed from INEX, provides coupling between the subsystems models and incorporates application models relevant to a specific trade-off or design study.
Integrated Numerical Experiments (INEX) and the Free-Electron Laser Physical Process Code (FELPPC)
NASA Astrophysics Data System (ADS)
Thode, L. E.; Chan, K. C. D.; Schmitt, M. J.; McKee, J.; Ostic, J.; Elliott, C. J.; McVey, B. D.
The strong coupling of subsystem elements, such as the accelerator, wiggler, and optics, greatly complicates the understanding and design of a free electron laser (FEL), even at the conceptual level. Given the requirements for high-performance FELs, the strong coupling between the laser subsystems must be included to obtain a realistic picture of the potential operational capability. To address the strong coupling character of the FEL the concept of an Integrated Numerical Experiment (INEX) was proposed. Unique features of the INEX approach are consistency and numerical equivalence of experimental diagnostics. The equivalent numerical diagnostics mitigates the major problem of misinterpretation that often occurs when theoretical and experimental data are compared. The INEX approach has been applied to a large number of accelerator and FEL experiments. Overall, the agreement between INEX and the experiments is very good. Despite the success of INEX, the approach is difficult to apply to trade-off and initial design studies because of the significant manpower and computational requirements. On the other hand, INEX provides a base from which realistic accelerator, wiggler, and optics models can be developed. The Free Electron Laser Physical Process Code (FELPPC) includes models developed from INEX, provides coupling between the subsystems models and incorporates application models relevant to a specific trade-off or design study.
Droegemeier, Kelvin K
2009-03-13
Mesoscale weather, such as convective systems, intense local rainfall resulting in flash floods and lake effect snows, frequently is characterized by unpredictable rapid onset and evolution, heterogeneity and spatial and temporal intermittency. Ironically, most of the technologies used to observe the atmosphere, predict its evolution and compute, transmit or store information about it, operate in a static pre-scheduled framework that is fundamentally inconsistent with, and does not accommodate, the dynamic behaviour of mesoscale weather. As a result, today's weather technology is highly constrained and far from optimal when applied to any particular situation. This paper describes a new cyberinfrastructure framework, in which remote and in situ atmospheric sensors, data acquisition and storage systems, assimilation and prediction codes, data mining and visualization engines, and the information technology frameworks within which they operate, can change configuration automatically, in response to evolving weather. Such dynamic adaptation is designed to allow system components to achieve greater overall effectiveness, relative to their static counterparts, for any given situation. The associated service-oriented architecture, known as Linked Environments for Atmospheric Discovery (LEAD), makes advanced meteorological and cyber tools as easy to use as ordering a book on the web. LEAD has been applied in a variety of settings, including experimental forecasting by the US National Weather Service, and allows users to focus much more attention on the problem at hand and less on the nuances of data formats, communication protocols and job execution environments. PMID:19087934
Essays on agricultural adaptation to climate change and ethanol market integration in the U.S
NASA Astrophysics Data System (ADS)
Aisabokhae, Ruth Ada
Climate factors like precipitation and temperature, being closely intertwined with agriculture, make a changing climate a big concern for the entire human race and its basic survival. Adaptation to climate is a long-running characteristic of agriculture evidenced by the varying types and forms of agricultural enterprises associated with differing climatic conditions. Nevertheless climate change poses a substantial, additional adaptation challenge for agriculture. Mitigation encompasses efforts to reduce the current and future extent of climate change. Biofuels production, for instance, expands agriculture's role in climate change mitigation. This dissertation encompasses adaptation and mitigation strategies as a response to climate change in the U.S. by examining comprehensively scientific findings on agricultural adaptation to climate change; developing information on the costs and benefits of select adaptations to examine what adaptations are most desirable, for which society can further devote its resources; and studying how ethanol prices are interrelated across, and transmitted within the U.S., and the markets that play an important role in these dynamics. Quantitative analysis using the Forestry and Agricultural Sector Optimization Model (FASOM) shows adaptation to be highly beneficial to agriculture. On-farm varietal and other adaptations contributions outweigh a mix shift northwards significantly, implying progressive technical change and significant returns to adaptation research and investment focused on farm management and varietal adaptations could be quite beneficial over time. Northward shift of corn-acre weighted centroids observed indicates that substantial production potential may shift across regions with the possibility of less production in the South, and more in the North, and thereby, potential redistribution of income. Time series techniques employed to study ethanol price dynamics show that the markets studied are co-integrated and strongly