Production facility layout by comparing moment displacement using BLOCPLAN and ALDEP Algorithms

NASA Astrophysics Data System (ADS)

Tambunan, M.; Ginting, E.; Sari, R. M.

2018-02-01

Production floor layout settings include the organizing of machinery, materials, and all the equipments used in the production process in the available area. PT. XYZ is a company that manufactures rubber and rubber compounds for retreading tire threaded with hot and cold cooking system. In the production of PT. XYZ is divided into three interrelated parts, namely Masterbatch Department, Department Compound, and Procured Thread Line Department. PT. XYZ has a production process with material flow is irregular and the arrangement of machine is complicated and need to be redesigned. The purpose of this study is comparing movement displacement using BLOCPLAN and ALDEP algorithm in order to redesign existing layout. Redesigning the layout of the production floor is done by applying algorithms of BLOCPLAN and ALDEP. The algorithm used to find the best layout design by comparing the moment displacement and the flow pattern. Moment displacement on the floor layout of the company’s production currently amounts to 2,090,578.5 meters per year and material flow pattern is irregular. Based on the calculation, the moment displacement for the BLOCPLAN is 1,551,344.82 meter per year and ALDEP is 1,600,179 meter per year. Flow Material resulted is in the form of straight the line.

Implementation and Initial Testing of Advanced Processing and Analysis Algorithms for Correlated Neutron Counting

DOE Office of Scientific and Technical Information (OSTI.GOV)

Santi, Peter Angelo; Cutler, Theresa Elizabeth; Favalli, Andrea

In order to improve the accuracy and capabilities of neutron multiplicity counting, additional quantifiable information is needed in order to address the assumptions that are present in the point model. Extracting and utilizing higher order moments (Quads and Pents) from the neutron pulse train represents the most direct way of extracting additional information from the measurement data to allow for an improved determination of the physical properties of the item of interest. The extraction of higher order moments from a neutron pulse train required the development of advanced dead time correction algorithms which could correct for dead time effects inmore » all of the measurement moments in a self-consistent manner. In addition, advanced analysis algorithms have been developed to address specific assumptions that are made within the current analysis model, namely that all neutrons are created at a single point within the item of interest, and that all neutrons that are produced within an item are created with the same energy distribution. This report will discuss the current status of implementation and initial testing of the advanced dead time correction and analysis algorithms that have been developed in an attempt to utilize higher order moments to improve the capabilities of correlated neutron measurement techniques.« less

Regional skew for California, and flood frequency for selected sites in the Sacramento-San Joaquin River Basin, based on data through water year 2006

USGS Publications Warehouse

Parrett, Charles; Veilleux, Andrea; Stedinger, J.R.; Barth, N.A.; Knifong, Donna L.; Ferris, J.C.

2011-01-01

Improved flood-frequency information is important throughout California in general and in the Sacramento-San Joaquin River Basin in particular, because of an extensive network of flood-control levees and the risk of catastrophic flooding. A key first step in updating flood-frequency information is determining regional skew. A Bayesian generalized least squares (GLS) regression method was used to derive a regional-skew model based on annual peak-discharge data for 158 long-term (30 or more years of record) stations throughout most of California. The desert areas in southeastern California had too few long-term stations to reliably determine regional skew for that hydrologically distinct region; therefore, the desert areas were excluded from the regional skew analysis for California. Of the 158 long-term stations used to determine regional skew, 145 have minimally regulated annual-peak discharges, and 13 stations are dam sites for which unregulated peak discharges were estimated from unregulated daily maximum discharge data furnished by the U.S. Army Corp of Engineers. Station skew was determined by using an expected moments algorithm (EMA) program for fitting the Pearson Type 3 flood-frequency distribution to the logarithms of annual peak-discharge data. The Bayesian GLS regression method previously developed was modified because of the large cross correlations among concurrent recorded peak discharges in California and the use of censored data and historical flood information with the new expected moments algorithm. In particular, to properly account for these cross-correlation problems and develop a suitable regression model and regression diagnostics, a combination of Bayesian weighted least squares and generalized least squares regression was adopted. This new methodology identified a nonlinear function relating regional skew to mean basin elevation. The regional skew values ranged from -0.62 for a mean basin elevation of zero to 0.61 for a mean basin elevation of 11,000 feet. This relation between skew and elevation reflects the interaction of snow with rain, which increases with increased elevation. The equivalent record length for the new regional skew ranges from 52 to 65 years of record, depending upon mean basin elevation. The old regional skew map in Bulletin 17B, published by the Hydrology Subcommittee of the Interagency Advisory Committee on Water Data (1982), reported an equivalent record length of only 17 years. The newly developed regional skew relation for California was used to update flood frequency for the 158 sites used in the regional skew analysis as well as 206 selected sites in the Sacramento-San Joaquin River Basin. For these sites, annual-peak discharges having recurrence intervals of 2, 5, 10, 25, 50, 100, 200, and 500 years were determined on the basis of data through water year 2006. The expected moments algorithm was used for determining the magnitude and frequency of floods at gaged sites by using regional skew values and using the basic approach outlined in Bulletin

Relationships between expected, online and remembered enjoyment for food products.

PubMed

Robinson, Eric

2014-03-01

How enjoyable a food product is remembered to be is likely to shape future choice. The present study tested the influence that expectations and specific moments during consumption experiences have on remembered enjoyment for food products. Sixty-four participants consumed three snack foods (savoury, sweet and savoury-sweet) and rated expected and online enjoyment for each product. Twenty-four hours later participants rated remembered enjoyment and future expected enjoyment for each product. Remembered enjoyment differed to online enjoyment for two of the three products, resulting in the foods being remembered as less enjoyable than they actually were. Both expected enjoyment and specific moments during the consumption experience (e.g. the least enjoyable mouthful) influenced remembered enjoyment. However, the factors that shaped remembered enjoyment were not consistent across the different food products. Remembered enjoyment was also shown to be a better predictor of future expected enjoyment than online enjoyment. Remembered enjoyment is likely to influence choice behaviour and can be discrepant to actual enjoyment. Specific moments during a consumption experience can have disproportionately large influence on remembered enjoyment (whilst others are neglected), but the factors that determine which moments influence remembered enjoyment are unclear. Copyright © 2013 Elsevier Ltd. All rights reserved.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Pebay, Philippe; Terriberry, Timothy B.; Kolla, Hemanth

Formulas for incremental or parallel computation of second order central moments have long been known, and recent extensions of these formulas to univariate and multivariate moments of arbitrary order have been developed. Formulas such as these, are of key importance in scenarios where incremental results are required and in parallel and distributed systems where communication costs are high. We survey these recent results, and improve them with arbitrary-order, numerically stable one-pass formulas which we further extend with weighted and compound variants. We also develop a generalized correction factor for standard two-pass algorithms that enables the maintenance of accuracy over nearlymore » the full representable range of the input, avoiding the need for extended-precision arithmetic. We then empirically examine algorithm correctness for pairwise update formulas up to order four as well as condition number and relative error bounds for eight different central moment formulas, each up to degree six, to address the trade-offs between numerical accuracy and speed of the various algorithms. Finally, we demonstrate the use of the most elaborate among the above mentioned formulas, with the utilization of the compound moments for a practical large-scale scientific application.« less

DOE Office of Scientific and Technical Information (OSTI.GOV)

Pebay, Philippe; Terriberry, Timothy B.; Kolla, Hemanth

Formulas for incremental or parallel computation of second order central moments have long been known, and recent extensions of these formulas to univariate and multivariate moments of arbitrary order have been developed. Such formulas are of key importance in scenarios where incremental results are required and in parallel and distributed systems where communication costs are high. We survey these recent results, and improve them with arbitrary-order, numerically stable one-pass formulas which we further extend with weighted and compound variants. We also develop a generalized correction factor for standard two-pass algorithms that enables the maintenance of accuracy over nearly the fullmore » representable range of the input, avoiding the need for extended-precision arithmetic. We then empirically examine algorithm correctness for pairwise update formulas up to order four as well as condition number and relative error bounds for eight different central moment formulas, each up to degree six, to address the trade-offs between numerical accuracy and speed of the various algorithms. Finally, we demonstrate the use of the most elaborate among the above mentioned formulas, with the utilization of the compound moments for a practical large-scale scientific application.« less

Application of temporal moments and other signal processing algorithms to analysis of ultrasonic signals through melting wax

DOE PAGES

Lau, Sarah J.; Moore, David G.; Stair, Sarah L.; ...

2016-01-01

Ultrasonic analysis is being explored as a way to capture events during melting of highly dispersive wax. Typical events include temperature changes in the material, phase transition of the material, surface flows and reformations, and void filling as the material melts. Melt tests are performed with wax to evaluate the usefulness of different signal processing algorithms in capturing event data. Several algorithm paths are being pursued. The first looks at changes in the velocity of the signal through the material. This is only appropriate when the changes from one ultrasonic signal to the next can be represented by a linearmore » relationship, which is not always the case. The second tracks changes in the frequency content of the signal. The third algorithm tracks changes in the temporal moments of a signal over a full test. This method does not require that the changes in the signal be represented by a linear relationship, but attaching changes in the temporal moments to physical events can be difficult. This study describes the algorithm paths applied to experimental data from ultrasonic signals as wax melts and explores different ways to display the results.« less

Optimum Actuator Selection with a Genetic Algorithm for Aircraft Control

NASA Technical Reports Server (NTRS)

Rogers, James L.

2004-01-01

The placement of actuators on a wing determines the control effectiveness of the airplane. One approach to placement maximizes the moments about the pitch, roll, and yaw axes, while minimizing the coupling. For example, the desired actuators produce a pure roll moment without at the same time causing much pitch or yaw. For a typical wing, there is a large set of candidate locations for placing actuators, resulting in a substantially larger number of combinations to examine in order to find an optimum placement satisfying the mission requirements and mission constraints. A genetic algorithm has been developed for finding the best placement for four actuators to produce an uncoupled pitch moment. The genetic algorithm has been extended to find the minimum number of actuators required to provide uncoupled pitch, roll, and yaw control. A simplified, untapered, unswept wing is the model for each application.

Estimation of full moment tensors, including uncertainties, for earthquakes, volcanic events, and nuclear explosions

NASA Astrophysics Data System (ADS)

Alvizuri, Celso; Silwal, Vipul; Krischer, Lion; Tape, Carl

2017-04-01

A seismic moment tensor is a 3 × 3 symmetric matrix that provides a compact representation of seismic events within Earth's crust. We develop an algorithm to estimate moment tensors and their uncertainties from observed seismic data. For a given event, the algorithm performs a grid search over the six-dimensional space of moment tensors by generating synthetic waveforms at each grid point and then evaluating a misfit function between the observed and synthetic waveforms. 'The' moment tensor M for the event is then the moment tensor with minimum misfit. To describe the uncertainty associated with M, we first convert the misfit function to a probability function. The uncertainty, or rather the confidence, is then given by the 'confidence curve' P(V ), where P(V ) is the probability that the true moment tensor for the event lies within the neighborhood of M that has fractional volume V . The area under the confidence curve provides a single, abbreviated 'confidence parameter' for M. We apply the method to data from events in different regions and tectonic settings: small (Mw < 2.5) events at Uturuncu volcano in Bolivia, moderate (Mw > 4) earthquakes in the southern Alaska subduction zone, and natural and man-made events at the Nevada Test Site. Moment tensor uncertainties allow us to better discriminate among moment tensor source types and to assign physical processes to the events.

Toward an autonomous brain machine interface: integrating sensorimotor reward modulation and reinforcement learning.

PubMed

Marsh, Brandi T; Tarigoppula, Venkata S Aditya; Chen, Chen; Francis, Joseph T

2015-05-13

For decades, neurophysiologists have worked on elucidating the function of the cortical sensorimotor control system from the standpoint of kinematics or dynamics. Recently, computational neuroscientists have developed models that can emulate changes seen in the primary motor cortex during learning. However, these simulations rely on the existence of a reward-like signal in the primary sensorimotor cortex. Reward modulation of the primary sensorimotor cortex has yet to be characterized at the level of neural units. Here we demonstrate that single units/multiunits and local field potentials in the primary motor (M1) cortex of nonhuman primates (Macaca radiata) are modulated by reward expectation during reaching movements and that this modulation is present even while subjects passively view cursor motions that are predictive of either reward or nonreward. After establishing this reward modulation, we set out to determine whether we could correctly classify rewarding versus nonrewarding trials, on a moment-to-moment basis. This reward information could then be used in collaboration with reinforcement learning principles toward an autonomous brain-machine interface. The autonomous brain-machine interface would use M1 for both decoding movement intention and extraction of reward expectation information as evaluative feedback, which would then update the decoding algorithm as necessary. In the work presented here, we show that this, in theory, is possible. Copyright © 2015 the authors 0270-6474/15/357374-14$15.00/0.

Sensitivity Analysis of Linear Programming and Quadratic Programming Algorithms for Control Allocation

NASA Technical Reports Server (NTRS)

Frost, Susan A.; Bodson, Marc; Acosta, Diana M.

2009-01-01

The Next Generation (NextGen) transport aircraft configurations being investigated as part of the NASA Aeronautics Subsonic Fixed Wing Project have more control surfaces, or control effectors, than existing transport aircraft configurations. Conventional flight control is achieved through two symmetric elevators, two antisymmetric ailerons, and a rudder. The five effectors, reduced to three command variables, produce moments along the three main axes of the aircraft and enable the pilot to control the attitude and flight path of the aircraft. The NextGen aircraft will have additional redundant control effectors to control the three moments, creating a situation where the aircraft is over-actuated and where a simple relationship does not exist anymore between the required effector deflections and the desired moments. NextGen flight controllers will incorporate control allocation algorithms to determine the optimal effector commands and attain the desired moments, taking into account the effector limits. Approaches to solving the problem using linear programming and quadratic programming algorithms have been proposed and tested. It is of great interest to understand their relative advantages and disadvantages and how design parameters may affect their properties. In this paper, we investigate the sensitivity of the effector commands with respect to the desired moments and show on some examples that the solutions provided using the l2 norm of quadratic programming are less sensitive than those using the l1 norm of linear programming.

A tyre slip-based integrated chassis control of front/rear traction distribution and four-wheel independent brake from moderate driving to limit handling

NASA Astrophysics Data System (ADS)

Joa, Eunhyek; Park, Kwanwoo; Koh, Youngil; Yi, Kyongsu; Kim, Kilsoo

2018-04-01

This paper presents a tyre slip-based integrated chassis control of front/rear traction distribution and four-wheel braking for enhanced performance from moderate driving to limit handling. The proposed algorithm adopted hierarchical structure: supervisor - desired motion tracking controller - optimisation-based control allocation. In the supervisor, by considering transient cornering characteristics, desired vehicle motion is calculated. In the desired motion tracking controller, in order to track desired vehicle motion, virtual control input is determined in the manner of sliding mode control. In the control allocation, virtual control input is allocated to minimise cost function. The cost function consists of two major parts. First part is a slip-based tyre friction utilisation quantification, which does not need a tyre force estimation. Second part is an allocation guideline, which guides optimally allocated inputs to predefined solution. The proposed algorithm has been investigated via simulation from moderate driving to limit handling scenario. Compared to Base and direct yaw moment control system, the proposed algorithm can effectively reduce tyre dissipation energy in the moderate driving situation. Moreover, the proposed algorithm enhances limit handling performance compared to Base and direct yaw moment control system. In addition to comparison with Base and direct yaw moment control, comparison the proposed algorithm with the control algorithm based on the known tyre force information has been conducted. The results show that the performance of the proposed algorithm is similar with that of the control algorithm with the known tyre force information.

Performance of the split-symbol moments SNR estimator in the presence of inter-symbol interference

NASA Technical Reports Server (NTRS)

Shah, B.; Hinedi, S.

1989-01-01

The Split-Symbol Moments Estimator (SSME) is an algorithm that is designed to estimate symbol signal-to-noise ratio (SNR) in the presence of additive white Gaussian noise (AWGN). The performance of the SSME algorithm in band-limited channels is examined. The effects of the resulting inter-symbol interference (ISI) are quantified. All results obtained are in closed form and can be easily evaluated numerically for performance prediction purposes. Furthermore, they are validated through digital simulations.

Multi-fidelity stochastic collocation method for computation of statistical moments

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhu, Xueyu, E-mail: xueyu-zhu@uiowa.edu; Linebarger, Erin M., E-mail: aerinline@sci.utah.edu; Xiu, Dongbin, E-mail: xiu.16@osu.edu

We present an efficient numerical algorithm to approximate the statistical moments of stochastic problems, in the presence of models with different fidelities. The method extends the multi-fidelity approximation method developed in . By combining the efficiency of low-fidelity models and the accuracy of high-fidelity models, our method exhibits fast convergence with a limited number of high-fidelity simulations. We establish an error bound of the method and present several numerical examples to demonstrate the efficiency and applicability of the multi-fidelity algorithm.

Elementary quantum mechanics of the neutron with an electric dipole moment

PubMed Central

Baym, Gordon; Beck, D. H.

2016-01-01

The neutron, in addition to possibly having a permanent electric dipole moment as a consequence of violation of time-reversal invariance, develops an induced electric dipole moment in the presence of an external electric field. We present here a unified nonrelativistic description of these two phenomena, in which the dipole moment operator, D→, is not constrained to lie along the spin operator. Although the expectation value of D→ in the neutron is less than 10−13 of the neutron radius, rn, the expectation value of D→ 2 is of order rn2. We determine the spin motion in external electric and magnetic fields, as used in past and future searches for a permanent dipole moment, and show that the neutron electric polarizability, although entering the neutron energy in an external electric field, does not affect the spin motion. In a simple nonrelativistic model we show that the expectation value of the permanent dipole is, to lowest order, proportional to the product of the time-reversal-violating coupling strength and the electric polarizability of the neutron. PMID:27325765

Elementary quantum mechanics of the neutron with an electric dipole moment.

PubMed

Baym, Gordon; Beck, D H

2016-07-05

The neutron, in addition to possibly having a permanent electric dipole moment as a consequence of violation of time-reversal invariance, develops an induced electric dipole moment in the presence of an external electric field. We present here a unified nonrelativistic description of these two phenomena, in which the dipole moment operator, [Formula: see text], is not constrained to lie along the spin operator. Although the expectation value of [Formula: see text] in the neutron is less than [Formula: see text] of the neutron radius, [Formula: see text], the expectation value of [Formula: see text] is of order [Formula: see text] We determine the spin motion in external electric and magnetic fields, as used in past and future searches for a permanent dipole moment, and show that the neutron electric polarizability, although entering the neutron energy in an external electric field, does not affect the spin motion. In a simple nonrelativistic model we show that the expectation value of the permanent dipole is, to lowest order, proportional to the product of the time-reversal-violating coupling strength and the electric polarizability of the neutron.

On the pth moment estimates of solutions to stochastic functional differential equations in the G-framework.

PubMed

Faizullah, Faiz

2016-01-01

The aim of the current paper is to present the path-wise and moment estimates for solutions to stochastic functional differential equations with non-linear growth condition in the framework of G-expectation and G-Brownian motion. Under the nonlinear growth condition, the pth moment estimates for solutions to SFDEs driven by G-Brownian motion are proved. The properties of G-expectations, Hölder's inequality, Bihari's inequality, Gronwall's inequality and Burkholder-Davis-Gundy inequalities are used to develop the above mentioned theory. In addition, the path-wise asymptotic estimates and continuity of pth moment for the solutions to SFDEs in the G-framework, with non-linear growth condition are shown.

Novel structures for Discrete Hartley Transform based on first-order moments

NASA Astrophysics Data System (ADS)

Xiong, Jun; Zheng, Wenjuan; Wang, Hao; Liu, Jianguo

2018-03-01

Discrete Hartley Transform (DHT) is an important tool in digital signal processing. In the present paper, the DHT is firstly transformed into the first-order moments-based form, then a new fast algorithm is proposed to calculate the first-order moments without multiplication. Based on the algorithm theory, the corresponding hardware architecture for DHT is proposed, which only contains shift operations and additions with no need for multipliers and large memory. To verify the availability and effectiveness, the proposed design is implemented with hardware description language and synthesized by Synopsys Design Compiler with 0.18-μm SMIC library. A series of experiments have proved that the proposed architecture has better performance in terms of the product of the hardware consumption and computation time.

Multifractal detrending moving-average cross-correlation analysis

NASA Astrophysics Data System (ADS)

Jiang, Zhi-Qiang; Zhou, Wei-Xing

2011-07-01

There are a number of situations in which several signals are simultaneously recorded in complex systems, which exhibit long-term power-law cross correlations. The multifractal detrended cross-correlation analysis (MFDCCA) approaches can be used to quantify such cross correlations, such as the MFDCCA based on the detrended fluctuation analysis (MFXDFA) method. We develop in this work a class of MFDCCA algorithms based on the detrending moving-average analysis, called MFXDMA. The performances of the proposed MFXDMA algorithms are compared with the MFXDFA method by extensive numerical experiments on pairs of time series generated from bivariate fractional Brownian motions, two-component autoregressive fractionally integrated moving-average processes, and binomial measures, which have theoretical expressions of the multifractal nature. In all cases, the scaling exponents hxy extracted from the MFXDMA and MFXDFA algorithms are very close to the theoretical values. For bivariate fractional Brownian motions, the scaling exponent of the cross correlation is independent of the cross-correlation coefficient between two time series, and the MFXDFA and centered MFXDMA algorithms have comparative performances, which outperform the forward and backward MFXDMA algorithms. For two-component autoregressive fractionally integrated moving-average processes, we also find that the MFXDFA and centered MFXDMA algorithms have comparative performances, while the forward and backward MFXDMA algorithms perform slightly worse. For binomial measures, the forward MFXDMA algorithm exhibits the best performance, the centered MFXDMA algorithms performs worst, and the backward MFXDMA algorithm outperforms the MFXDFA algorithm when the moment order q<0 and underperforms when q>0. We apply these algorithms to the return time series of two stock market indexes and to their volatilities. For the returns, the centered MFXDMA algorithm gives the best estimates of hxy(q) since its hxy(2) is closest to 0.5, as expected, and the MFXDFA algorithm has the second best performance. For the volatilities, the forward and backward MFXDMA algorithms give similar results, while the centered MFXDMA and the MFXDFA algorithms fail to extract rational multifractal nature.

Damping Rotor Nutation Oscillations in a Gyroscope with Magnetic Suspension

NASA Technical Reports Server (NTRS)

Komarov, Valentine N.

1996-01-01

A possibility of an effective damping of rotor nutations by modulating the field of the moment transducers in synchronism with the nutation frequency is considered. The algorithms for forming the control moments are proposed and their application is discussed.

Steering law design for redundant single-gimbal control moment gyroscopes. [for spacecraft attitude control

NASA Technical Reports Server (NTRS)

Bedrossian, Nazareth S.; Paradiso, Joseph; Bergmann, Edward V.; Rowell, Derek

1990-01-01

Two steering laws are presented for single-gimbal control moment gyroscopes. An approach using the Moore-Penrose pseudoinverse with a nondirectional null-motion algorithm is shown by example to avoid internal singularities for unidirectional torque commands, for which existing algorithms fail. Because this is still a tangent-based approach, however, singularity avoidance cannot be guaranteed. The singularity robust inverse is introduced as an alternative to the pseudoinverse for computing torque-producing gimbal rates near singular states. This approach, coupled with the nondirectional null algorithm, is shown by example to provide better steering law performance by allowing torque errors to be produced in the vicinity of singular states.

Calculating Higher-Order Moments of Phylogenetic Stochastic Mapping Summaries in Linear Time.

PubMed

Dhar, Amrit; Minin, Vladimir N

2017-05-01

Stochastic mapping is a simulation-based method for probabilistically mapping substitution histories onto phylogenies according to continuous-time Markov models of evolution. This technique can be used to infer properties of the evolutionary process on the phylogeny and, unlike parsimony-based mapping, conditions on the observed data to randomly draw substitution mappings that do not necessarily require the minimum number of events on a tree. Most stochastic mapping applications simulate substitution mappings only to estimate the mean and/or variance of two commonly used mapping summaries: the number of particular types of substitutions (labeled substitution counts) and the time spent in a particular group of states (labeled dwelling times) on the tree. Fast, simulation-free algorithms for calculating the mean of stochastic mapping summaries exist. Importantly, these algorithms scale linearly in the number of tips/leaves of the phylogenetic tree. However, to our knowledge, no such algorithm exists for calculating higher-order moments of stochastic mapping summaries. We present one such simulation-free dynamic programming algorithm that calculates prior and posterior mapping variances and scales linearly in the number of phylogeny tips. Our procedure suggests a general framework that can be used to efficiently compute higher-order moments of stochastic mapping summaries without simulations. We demonstrate the usefulness of our algorithm by extending previously developed statistical tests for rate variation across sites and for detecting evolutionarily conserved regions in genomic sequences.

Determining residual reduction algorithm kinematic tracking weights for a sidestep cut via numerical optimization.

PubMed

Samaan, Michael A; Weinhandl, Joshua T; Bawab, Sebastian Y; Ringleb, Stacie I

2016-12-01

Musculoskeletal modeling allows for the determination of various parameters during dynamic maneuvers by using in vivo kinematic and ground reaction force (GRF) data as inputs. Differences between experimental and model marker data and inconsistencies in the GRFs applied to these musculoskeletal models may not produce accurate simulations. Therefore, residual forces and moments are applied to these models in order to reduce these differences. Numerical optimization techniques can be used to determine optimal tracking weights of each degree of freedom of a musculoskeletal model in order to reduce differences between the experimental and model marker data as well as residual forces and moments. In this study, the particle swarm optimization (PSO) and simplex simulated annealing (SIMPSA) algorithms were used to determine optimal tracking weights for the simulation of a sidestep cut. The PSO and SIMPSA algorithms were able to produce model kinematics that were within 1.4° of experimental kinematics with residual forces and moments of less than 10 N and 18 Nm, respectively. The PSO algorithm was able to replicate the experimental kinematic data more closely and produce more dynamically consistent kinematic data for a sidestep cut compared to the SIMPSA algorithm. Future studies should use external optimization routines to determine dynamically consistent kinematic data and report the differences between experimental and model data for these musculoskeletal simulations.

Calculating Higher-Order Moments of Phylogenetic Stochastic Mapping Summaries in Linear Time

PubMed Central

Dhar, Amrit

2017-01-01

Abstract Stochastic mapping is a simulation-based method for probabilistically mapping substitution histories onto phylogenies according to continuous-time Markov models of evolution. This technique can be used to infer properties of the evolutionary process on the phylogeny and, unlike parsimony-based mapping, conditions on the observed data to randomly draw substitution mappings that do not necessarily require the minimum number of events on a tree. Most stochastic mapping applications simulate substitution mappings only to estimate the mean and/or variance of two commonly used mapping summaries: the number of particular types of substitutions (labeled substitution counts) and the time spent in a particular group of states (labeled dwelling times) on the tree. Fast, simulation-free algorithms for calculating the mean of stochastic mapping summaries exist. Importantly, these algorithms scale linearly in the number of tips/leaves of the phylogenetic tree. However, to our knowledge, no such algorithm exists for calculating higher-order moments of stochastic mapping summaries. We present one such simulation-free dynamic programming algorithm that calculates prior and posterior mapping variances and scales linearly in the number of phylogeny tips. Our procedure suggests a general framework that can be used to efficiently compute higher-order moments of stochastic mapping summaries without simulations. We demonstrate the usefulness of our algorithm by extending previously developed statistical tests for rate variation across sites and for detecting evolutionarily conserved regions in genomic sequences. PMID:28177780

Analysis of Sting Balance Calibration Data Using Optimized Regression Models

NASA Technical Reports Server (NTRS)

Ulbrich, N.; Bader, Jon B.

2010-01-01

Calibration data of a wind tunnel sting balance was processed using a candidate math model search algorithm that recommends an optimized regression model for the data analysis. During the calibration the normal force and the moment at the balance moment center were selected as independent calibration variables. The sting balance itself had two moment gages. Therefore, after analyzing the connection between calibration loads and gage outputs, it was decided to choose the difference and the sum of the gage outputs as the two responses that best describe the behavior of the balance. The math model search algorithm was applied to these two responses. An optimized regression model was obtained for each response. Classical strain gage balance load transformations and the equations of the deflection of a cantilever beam under load are used to show that the search algorithm s two optimized regression models are supported by a theoretical analysis of the relationship between the applied calibration loads and the measured gage outputs. The analysis of the sting balance calibration data set is a rare example of a situation when terms of a regression model of a balance can directly be derived from first principles of physics. In addition, it is interesting to note that the search algorithm recommended the correct regression model term combinations using only a set of statistical quality metrics that were applied to the experimental data during the algorithm s term selection process.

Moment-tensor solutions for the 24 November 1987 Superstition Hills, California, earthquakes

USGS Publications Warehouse

Sipkin, S.A.

1989-01-01

The teleseismic long-period waveforms recorded by the Global Digital Seismograph Network from the two largest Superstition Hills earthquakes are inverted using an algorithm based on optimal filter theory. These solutions differ slightly from those published in the Preliminary Determination of Epicenters Monthly Listing because a somewhat different, improved data set was used in the inversions and a time-dependent moment-tensor algorithm was used to investigate the complexity of the main shock. The foreshock (origin time 01:54:14.5, mb 5.7, Ms6.2) had a scalar moment of 2.3 ?? 1025 dyne-cm, a depth of 8km, and a mechanism of strike 217??, dip 79??, rake 4??. The main shock (origin time 13:15:56.4, mb 6.0, Ms6.6) was a complex event, consisting of at least two subevents, with a combined scalar moment of 1.0 ?? 1026 dyne-cm, a depth of 10km, and a mechanism of strike 303??, dip 89??, rake -180??. -Authors

Small sample estimation of the reliability function for technical products

NASA Astrophysics Data System (ADS)

Lyamets, L. L.; Yakimenko, I. V.; Kanishchev, O. A.; Bliznyuk, O. A.

2017-12-01

It is demonstrated that, in the absence of big statistic samples obtained as a result of testing complex technical products for failure, statistic estimation of the reliability function of initial elements can be made by the moments method. A formal description of the moments method is given and its advantages in the analysis of small censored samples are discussed. A modified algorithm is proposed for the implementation of the moments method with the use of only the moments at which the failures of initial elements occur.

Extension of moment projection method to the fragmentation process

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wu, Shaohua; Yapp, Edward K.Y.; Akroyd, Jethro

2017-04-15

The method of moments is a simple but efficient method of solving the population balance equation which describes particle dynamics. Recently, the moment projection method (MPM) was proposed and validated for particle inception, coagulation, growth and, more importantly, shrinkage; here the method is extended to include the fragmentation process. The performance of MPM is tested for 13 different test cases for different fragmentation kernels, fragment distribution functions and initial conditions. Comparisons are made with the quadrature method of moments (QMOM), hybrid method of moments (HMOM) and a high-precision stochastic solution calculated using the established direct simulation algorithm (DSA) and advantagesmore » of MPM are drawn.« less

Biologically Inspired Strategies, Algorithms and Hardware for Visual Guidance of Autonomous Helicopters

DTIC Science & Technology

2011-05-02

Dacke, J. Reinhard and M.V. Srinivasan (2010) The moment before touchdown: Landing manoeuvres of the honeybee Apis mellifera . Journal of Experimental...moment before touchdown: Landing manoeuvres of the honeybee Apis mellifera . Journal of Experimental Biology 213, 262-270. M.V. Srinivasan (2010

Fluid preconditioning for Newton–Krylov-based, fully implicit, electrostatic particle-in-cell simulations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chen, G., E-mail: gchen@lanl.gov; Chacón, L.; Leibs, C.A.

2014-02-01

A recent proof-of-principle study proposes an energy- and charge-conserving, nonlinearly implicit electrostatic particle-in-cell (PIC) algorithm in one dimension [9]. The algorithm in the reference employs an unpreconditioned Jacobian-free Newton–Krylov method, which ensures nonlinear convergence at every timestep (resolving the dynamical timescale of interest). Kinetic enslavement, which is one key component of the algorithm, not only enables fully implicit PIC as a practical approach, but also allows preconditioning the kinetic solver with a fluid approximation. This study proposes such a preconditioner, in which the linearized moment equations are closed with moments computed from particles. Effective acceleration of the linear GMRES solvemore » is demonstrated, on both uniform and non-uniform meshes. The algorithm performance is largely insensitive to the electron–ion mass ratio. Numerical experiments are performed on a 1D multi-scale ion acoustic wave test problem.« less

Radiation and scattering from bodies of translation. Volume 2: User's manual, computer program documentation

NASA Astrophysics Data System (ADS)

Medgyesi-Mitschang, L. N.; Putnam, J. M.

1980-04-01

A hierarchy of computer programs implementing the method of moments for bodies of translation (MM/BOT) is described. The algorithm treats the far-field radiation and scattering from finite-length open cylinders of arbitrary cross section as well as the near fields and aperture-coupled fields for rectangular apertures on such bodies. The theoretical development underlying the algorithm is described in Volume 1. The structure of the computer algorithm is such that no a priori knowledge of the method of moments technique or detailed FORTRAN experience are presupposed for the user. A set of carefully drawn example problems illustrates all the options of the algorithm. For more detailed understanding of the workings of the codes, special cross referencing to the equations in Volume 1 is provided. For additional clarity, comment statements are liberally interspersed in the code listings, summarized in the present volume.

Optimization and large scale computation of an entropy-based moment closure

NASA Astrophysics Data System (ADS)

Kristopher Garrett, C.; Hauck, Cory; Hill, Judith

2015-12-01

We present computational advances and results in the implementation of an entropy-based moment closure, MN, in the context of linear kinetic equations, with an emphasis on heterogeneous and large-scale computing platforms. Entropy-based closures are known in several cases to yield more accurate results than closures based on standard spectral approximations, such as PN, but the computational cost is generally much higher and often prohibitive. Several optimizations are introduced to improve the performance of entropy-based algorithms over previous implementations. These optimizations include the use of GPU acceleration and the exploitation of the mathematical properties of spherical harmonics, which are used as test functions in the moment formulation. To test the emerging high-performance computing paradigm of communication bound simulations, we present timing results at the largest computational scales currently available. These results show, in particular, load balancing issues in scaling the MN algorithm that do not appear for the PN algorithm. We also observe that in weak scaling tests, the ratio in time to solution of MN to PN decreases.

Extended AIC model based on high order moments and its application in the financial market

NASA Astrophysics Data System (ADS)

Mao, Xuegeng; Shang, Pengjian

2018-07-01

In this paper, an extended method of traditional Akaike Information Criteria(AIC) is proposed to detect the volatility of time series by combining it with higher order moments, such as skewness and kurtosis. Since measures considering higher order moments are powerful in many aspects, the properties of asymmetry and flatness can be observed. Furthermore, in order to reduce the effect of noise and other incoherent features, we combine the extended AIC algorithm with multiscale wavelet analysis, in which the newly extended AIC algorithm is applied to wavelet coefficients at several scales and the time series are reconstructed by wavelet transform. After that, we create AIC planes to derive the relationship among AIC values using variance, skewness and kurtosis respectively. When we test this technique on the financial market, the aim is to analyze the trend and volatility of the closing price of stock indices and classify them. And we also adapt multiscale analysis to measure complexity of time series over a range of scales. Empirical results show that the singularity of time series in stock market can be detected via extended AIC algorithm.

Optimization and large scale computation of an entropy-based moment closure

DOE PAGES

Hauck, Cory D.; Hill, Judith C.; Garrett, C. Kristopher

2015-09-10

We present computational advances and results in the implementation of an entropy-based moment closure, M N, in the context of linear kinetic equations, with an emphasis on heterogeneous and large-scale computing platforms. Entropy-based closures are known in several cases to yield more accurate results than closures based on standard spectral approximations, such as P N, but the computational cost is generally much higher and often prohibitive. Several optimizations are introduced to improve the performance of entropy-based algorithms over previous implementations. These optimizations include the use of GPU acceleration and the exploitation of the mathematical properties of spherical harmonics, which aremore » used as test functions in the moment formulation. To test the emerging high-performance computing paradigm of communication bound simulations, we present timing results at the largest computational scales currently available. Lastly, these results show, in particular, load balancing issues in scaling the M N algorithm that do not appear for the P N algorithm. We also observe that in weak scaling tests, the ratio in time to solution of M N to P N decreases.« less

Algorithms for Determining Physical Responses of Structures Under Load

NASA Technical Reports Server (NTRS)

Richards, W. Lance; Ko, William L.

2012-01-01

Ultra-efficient real-time structural monitoring algorithms have been developed to provide extensive information about the physical response of structures under load. These algorithms are driven by actual strain data to measure accurately local strains at multiple locations on the surface of a structure. Through a single point load calibration test, these structural strains are then used to calculate key physical properties of the structure at each measurement location. Such properties include the structure s flexural rigidity (the product of the structure's modulus of elasticity, and its moment of inertia) and the section modulus (the moment of inertia divided by the structure s half-depth). The resulting structural properties at each location can be used to determine the structure s bending moment, shear, and structural loads in real time while the structure is in service. The amount of structural information can be maximized through the use of highly multiplexed fiber Bragg grating technology using optical time domain reflectometry and optical frequency domain reflectometry, which can provide a local strain measurement every 10 mm on a single hair-sized optical fiber. Since local strain is used as input to the algorithms, this system serves multiple purposes of measuring strains and displacements, as well as determining structural bending moment, shear, and loads for assessing real-time structural health. The first step is to install a series of strain sensors on the structure s surface in such a way as to measure bending strains at desired locations. The next step is to perform a simple ground test calibration. For a beam of length l (see example), discretized into n sections and subjected to a tip load of P that places the beam in bending, the flexural rigidity of the beam can be experimentally determined at each measurement location x. The bending moment at each station can then be determined for any general set of loads applied during operation.

Time-frequency analysis-based time-windowing algorithm for the inverse synthetic aperture radar imaging of ships

NASA Astrophysics Data System (ADS)

Zhou, Peng; Zhang, Xi; Sun, Weifeng; Dai, Yongshou; Wan, Yong

2018-01-01

An algorithm based on time-frequency analysis is proposed to select an imaging time window for the inverse synthetic aperture radar imaging of ships. An appropriate range bin is selected to perform the time-frequency analysis after radial motion compensation. The selected range bin is that with the maximum mean amplitude among the range bins whose echoes are confirmed to be contributed by a dominant scatter. The criterion for judging whether the echoes of a range bin are contributed by a dominant scatter is key to the proposed algorithm and is therefore described in detail. When the first range bin that satisfies the judgment criterion is found, a sequence composed of the frequencies that have the largest amplitudes in every moment's time-frequency spectrum corresponding to this range bin is employed to calculate the length and the center moment of the optimal imaging time window. Experiments performed with simulation data and real data show the effectiveness of the proposed algorithm, and comparisons between the proposed algorithm and the image contrast-based algorithm (ICBA) are provided. Similar image contrast and lower entropy are acquired using the proposed algorithm as compared with those values when using the ICBA.

Infrared Ship Classification Using A New Moment Pattern Recognition Concept

NASA Astrophysics Data System (ADS)

Casasent, David; Pauly, John; Fetterly, Donald

1982-03-01

An analysis of the statistics of the moments and the conventional invariant moments shows that the variance of the latter become quite large as the order of the moments and the degree of invariance increases. Moreso, the need to whiten the error volume increases with the order and degree, but so does the computational load associated with computing the whitening operator. We thus advance a new estimation approach to the use of moments in pattern recog-nition that overcomes these problems. This work is supported by experimental verification and demonstration on an infrared ship pattern recognition problem. The computational load associated with our new algorithm is also shown to be very low.

Control of systematic uncertainties in the storage ring search for an electric dipole moment by measuring the electric quadrupole moment

NASA Astrophysics Data System (ADS)

Magiera, Andrzej

2017-09-01

Measurements of electric dipole moment (EDM) for light hadrons with use of a storage ring have been proposed. The expected effect is very small, therefore various subtle effects need to be considered. In particular, interaction of particle's magnetic dipole moment and electric quadrupole moment with electromagnetic field gradients can produce an effect of a similar order of magnitude as that expected for EDM. This paper describes a very promising method employing an rf Wien filter, allowing to disentangle that contribution from the genuine EDM effect. It is shown that both these effects could be separated by the proper setting of the rf Wien filter frequency and phase. In the EDM measurement the magnitude of systematic uncertainties plays a key role and they should be under strict control. It is shown that particles' interaction with field gradients offers also the possibility to estimate global systematic uncertainties with the precision necessary for an EDM measurement with the planned accuracy.

Parameterization of cloud lidar backscattering profiles by means of asymmetrical Gaussians

NASA Astrophysics Data System (ADS)

del Guasta, Massimo; Morandi, Marco; Stefanutti, Leopoldo

1995-06-01

A fitting procedure for cloud lidar data processing is shown that is based on the computation of the first three moments of the vertical-backscattering (or -extinction) profile. Single-peak clouds or single cloud layers are approximated to asymmetrical Gaussians. The algorithm is particularly stable with respect to noise and processing errors, and it is much faster than the equivalent least-squares approach. Multilayer clouds can easily be treated as a sum of single asymmetrical Gaussian peaks. The method is suitable for cloud-shape parametrization in noisy lidar signatures (like those expected from satellite lidars). It also permits an improvement of cloud radiative-property computations that are based on huge lidar data sets for which storage and careful examination of single lidar profiles can't be carried out.

Data-based comparisons of moments estimators using historical and paleoflood data

USGS Publications Warehouse

England, J.F.; Jarrett, R.D.; Salas, J.D.

2003-01-01

This paper presents the first systematic comparison, using historical and paleoflood data, of moments-based flood frequency methods. Peak flow estimates were compiled from streamflow-gaging stations with historical and/or paleoflood data at 36 sites located in the United States, Argentina, United Kingdom and China, covering a diverse range of hydrologic conditions. The Expected Moments Algorithm (EMA) and the Bulletin 17B historical weighting procedure (B17H) were compared in terms of goodness of fit using 25 of the data sets. Results from this comparison indicate that EMA is a viable alternative to current B17H procedures from an operational perspective, and performed equal to or better than B17H for the data analyzed. We demonstrate satisfactory EMA performance for the remaining 11 sites with multiple thresholds and binomial censoring, which B17H cannot accommodate. It is shown that the EMA estimator readily incorporates these types of information and the LP-III distribution provided an adequate fit to the data in most cases. The results shown here are consistent with Monte Carlo simulation studies, and demonstrate that EMA is preferred overall to B17H. The Bulletin 17B document could be revised to include an option for EMA as an alternative to the existing historical weighting approach. These results are of practical relevance to hydrologists and water resources managers for applications in floodplain management, design of hydraulic structures, and risk analysis for dams. ?? 2003 Elsevier Science B.V. All rights reserved.

Data-based comparisons of moments estimators using historical and paleoflood data

NASA Astrophysics Data System (ADS)

England, John F.; Jarrett, Robert D.; Salas, José D.

2003-07-01

This paper presents the first systematic comparison, using historical and paleoflood data, of moments-based flood frequency methods. Peak flow estimates were compiled from streamflow-gaging stations with historical and/or paleoflood data at 36 sites located in the United States, Argentina, United Kingdom and China, covering a diverse range of hydrologic conditions. The Expected Moments Algorithm (EMA) and the Bulletin 17B historical weighting procedure (B17H) were compared in terms of goodness of fit using 25 of the data sets. Results from this comparison indicate that EMA is a viable alternative to current B17H procedures from an operational perspective, and performed equal to or better than B17H for the data analyzed. We demonstrate satisfactory EMA performance for the remaining 11 sites with multiple thresholds and binomial censoring, which B17H cannot accommodate. It is shown that the EMA estimator readily incorporates these types of information and the LP-III distribution provided an adequate fit to the data in most cases. The results shown here are consistent with Monte Carlo simulation studies, and demonstrate that EMA is preferred overall to B17H. The Bulletin 17B document could be revised to include an option for EMA as an alternative to the existing historical weighting approach. These results are of practical relevance to hydrologists and water resources managers for applications in floodplain management, design of hydraulic structures, and risk analysis for dams.

Hardware-efficient implementation of digital FIR filter using fast first-order moment algorithm

NASA Astrophysics Data System (ADS)

Cao, Li; Liu, Jianguo; Xiong, Jun; Zhang, Jing

2018-03-01

As the digital finite impulse response (FIR) filter can be transformed into the shift-add form of multiple small-sized firstorder moments, based on the existing fast first-order moment algorithm, this paper presents a novel multiplier-less structure to calculate any number of sequential filtering results in parallel. The theoretical analysis on its hardware and time-complexities reveals that by appropriately setting the degree of parallelism and the decomposition factor of a fixed word width, the proposed structure may achieve better area-time efficiency than the existing two-dimensional (2-D) memoryless-based filter. To evaluate the performance concretely, the proposed designs for different taps along with the existing 2-D memoryless-based filters, are synthesized by Synopsys Design Compiler with 0.18-μm SMIC library. The comparisons show that the proposed design has less area-time complexity and power consumption when the number of filter taps is larger than 48.

Quadratures with multiple nodes, power orthogonality, and moment-preserving spline approximation

NASA Astrophysics Data System (ADS)

Milovanovic, Gradimir V.

2001-01-01

Quadrature formulas with multiple nodes, power orthogonality, and some applications of such quadratures to moment-preserving approximation by defective splines are considered. An account on power orthogonality (s- and [sigma]-orthogonal polynomials) and generalized Gaussian quadratures with multiple nodes, including stable algorithms for numerical construction of the corresponding polynomials and Cotes numbers, are given. In particular, the important case of Chebyshev weight is analyzed. Finally, some applications in moment-preserving approximation of functions by defective splines are discussed.

Comparison of recent S-wave indicating methods

NASA Astrophysics Data System (ADS)

Hubicka, Katarzyna; Sokolowski, Jakub

2018-01-01

Seismic event consists of surface waves and body waves. Due to the fact that the body waves are faster (P-waves) and more energetic (S-waves) in literature the problem of their analysis is taken more often. The most universal information that is received from the recorded wave is its moment of arrival. When this information is obtained from at least four seismometers in different locations, the epicentre of the particular event can be estimated [1]. Since the recorded body waves may overlap in signal, the problem of wave onset moment is considered more often for faster P-wave than S-wave. This however does not mean that the issue of S-wave arrival time is not taken at all. As the process of manual picking is time-consuming, methods of automatic detection are recommended (these however may be less accurate). In this paper four recently developed methods estimating S-wave arrival are compared: the method operating on empirical mode decomposition and Teager-Kaiser operator [2], the modification of STA/LTA algorithm [3], the method using a nearest neighbour-based approach [4] and the algorithm operating on characteristic of signals' second moments. The methods will be also compared to wellknown algorithm based on the autoregressive model [5]. The algorithms will be tested in terms of their S-wave arrival identification accuracy on real data originating from International Research Institutions for Seismology (IRIS) database.

Target recognition of ladar range images using slice image: comparison of four improved algorithms

NASA Astrophysics Data System (ADS)

Xia, Wenze; Han, Shaokun; Cao, Jingya; Wang, Liang; Zhai, Yu; Cheng, Yang

2017-07-01

Compared with traditional 3-D shape data, ladar range images possess properties of strong noise, shape degeneracy, and sparsity, which make feature extraction and representation difficult. The slice image is an effective feature descriptor to resolve this problem. We propose four improved algorithms on target recognition of ladar range images using slice image. In order to improve resolution invariance of the slice image, mean value detection instead of maximum value detection is applied in these four improved algorithms. In order to improve rotation invariance of the slice image, three new improved feature descriptors-which are feature slice image, slice-Zernike moments, and slice-Fourier moments-are applied to the last three improved algorithms, respectively. Backpropagation neural networks are used as feature classifiers in the last two improved algorithms. The performance of these four improved recognition systems is analyzed comprehensively in the aspects of the three invariances, recognition rate, and execution time. The final experiment results show that the improvements for these four algorithms reach the desired effect, the three invariances of feature descriptors are not directly related to the final recognition performance of recognition systems, and these four improved recognition systems have different performances under different conditions.

The algorithm of fast image stitching based on multi-feature extraction

NASA Astrophysics Data System (ADS)

Yang, Chunde; Wu, Ge; Shi, Jing

2018-05-01

This paper proposed an improved image registration method combining Hu-based invariant moment contour information and feature points detection, aiming to solve the problems in traditional image stitching algorithm, such as time-consuming feature points extraction process, redundant invalid information overload and inefficiency. First, use the neighborhood of pixels to extract the contour information, employing the Hu invariant moment as similarity measure to extract SIFT feature points in those similar regions. Then replace the Euclidean distance with Hellinger kernel function to improve the initial matching efficiency and get less mismatching points, further, estimate affine transformation matrix between the images. Finally, local color mapping method is adopted to solve uneven exposure, using the improved multiresolution fusion algorithm to fuse the mosaic images and realize seamless stitching. Experimental results confirm high accuracy and efficiency of method proposed in this paper.

The calculation of the mass moment of inertia of a fluid in a rotating rectangular tank

NASA Technical Reports Server (NTRS)

1977-01-01

This analysis calculated the mass moment of inertia of a nonviscous fluid in a slowly rotating rectangular tank. Given the dimensions of the tank in the x, y, and z coordinates, the axis of rotation, the percentage of the tank occupied by the fluid, and angle of rotation, an algorithm was written that could calculate the mass moment of inertia of the fluid. While not included in this paper, the change in the mass moment of inertia of the fluid could then be used to calculate the force exerted by the fluid on the container wall.

Compressible, multiphase semi-implicit method with moment of fluid interface representation

DOE PAGES

Jemison, Matthew; Sussman, Mark; Arienti, Marco

2014-09-16

A unified method for simulating multiphase flows using an exactly mass, momentum, and energy conserving Cell-Integrated Semi-Lagrangian advection algorithm is presented. The deforming material boundaries are represented using the moment-of-fluid method. Our new algorithm uses a semi-implicit pressure update scheme that asymptotically preserves the standard incompressible pressure projection method in the limit of infinite sound speed. The asymptotically preserving attribute makes the new method applicable to compressible and incompressible flows including stiff materials; enabling large time steps characteristic of incompressible flow algorithms rather than the small time steps required by explicit methods. Moreover, shocks are captured and material discontinuities aremore » tracked, without the aid of any approximate or exact Riemann solvers. As a result, wimulations of underwater explosions and fluid jetting in one, two, and three dimensions are presented which illustrate the effectiveness of the new algorithm at efficiently computing multiphase flows containing shock waves and material discontinuities with large “impedance mismatch.”« less

Advancements of in-flight mass moment of inertia and structural deflection algorithms for satellite attitude simulators

NASA Astrophysics Data System (ADS)

Wright, Jonathan W.

Experimental satellite attitude simulators have long been used to test and analyze control algorithms in order to drive down risk before implementation on an operational satellite. Ideally, the dynamic response of a terrestrial-based experimental satellite attitude simulator would be similar to that of an on-orbit satellite. Unfortunately, gravitational disturbance torques and poorly characterized moments of inertia introduce uncertainty into the system dynamics leading to questionable attitude control algorithm experimental results. This research consists of three distinct, but related contributions to the field of developing robust satellite attitude simulators. In the first part of this research, existing approaches to estimate mass moments and products of inertia are evaluated followed by a proposition and evaluation of a new approach that increases both the accuracy and precision of these estimates using typical on-board satellite sensors. Next, in order to better simulate the micro-torque environment of space, a new approach to mass balancing satellite attitude simulator is presented, experimentally evaluated, and verified. Finally, in the third area of research, we capitalize on the platform improvements to analyze a control moment gyroscope (CMG) singularity avoidance steering law. Several successful experiments were conducted with the CMG array at near-singular configurations. An evaluation process was implemented to verify that the platform remained near the desired test momentum, showing that the first two components of this research were effective in allowing us to conduct singularity avoidance experiments in a representative space-like test environment.

Monte Carlo closure for moment-based transport schemes in general relativistic radiation hydrodynamic simulations

NASA Astrophysics Data System (ADS)

Foucart, Francois

2018-04-01

General relativistic radiation hydrodynamic simulations are necessary to accurately model a number of astrophysical systems involving black holes and neutron stars. Photon transport plays a crucial role in radiatively dominated accretion discs, while neutrino transport is critical to core-collapse supernovae and to the modelling of electromagnetic transients and nucleosynthesis in neutron star mergers. However, evolving the full Boltzmann equations of radiative transport is extremely expensive. Here, we describe the implementation in the general relativistic SPEC code of a cheaper radiation hydrodynamic method that theoretically converges to a solution of Boltzmann's equation in the limit of infinite numerical resources. The algorithm is based on a grey two-moment scheme, in which we evolve the energy density and momentum density of the radiation. Two-moment schemes require a closure that fills in missing information about the energy spectrum and higher order moments of the radiation. Instead of the approximate analytical closure currently used in core-collapse and merger simulations, we complement the two-moment scheme with a low-accuracy Monte Carlo evolution. The Monte Carlo results can provide any or all of the missing information in the evolution of the moments, as desired by the user. As a first test of our methods, we study a set of idealized problems demonstrating that our algorithm performs significantly better than existing analytical closures. We also discuss the current limitations of our method, in particular open questions regarding the stability of the fully coupled scheme.

L-moments and TL-moments of the generalized lambda distribution

USGS Publications Warehouse

Asquith, W.H.

2007-01-01

The 4-parameter generalized lambda distribution (GLD) is a flexible distribution capable of mimicking the shapes of many distributions and data samples including those with heavy tails. The method of L-moments and the recently developed method of trimmed L-moments (TL-moments) are attractive techniques for parameter estimation for heavy-tailed distributions for which the L- and TL-moments have been defined. Analytical solutions for the first five L- and TL-moments in terms of GLD parameters are derived. Unfortunately, numerical methods are needed to compute the parameters from the L- or TL-moments. Algorithms are suggested for parameter estimation. Application of the GLD using both L- and TL-moment parameter estimates from example data is demonstrated, and comparison of the L-moment fit of the 4-parameter kappa distribution is made. A small simulation study of the 98th percentile (far-right tail) is conducted for a heavy-tail GLD with high-outlier contamination. The simulations show, with respect to estimation of the 98th-percent quantile, that TL-moments are less biased (more robost) in the presence of high-outlier contamination. However, the robustness comes at the expense of considerably more sampling variability. ?? 2006 Elsevier B.V. All rights reserved.

Development of an algorithm for controlling a multilevel three-phase converter

NASA Astrophysics Data System (ADS)

Taissariyeva, Kyrmyzy; Ilipbaeva, Lyazzat

2017-08-01

This work is devoted to the development of an algorithm for controlling transistors in a three-phase multilevel conversion system. The developed algorithm allows to organize a correct operation and describes the state of transistors at each moment of time when constructing a computer model of a three-phase multilevel converter. The developed algorithm of operation of transistors provides in-phase of a three-phase converter and obtaining a sinusoidal voltage curve at the converter output.

Computing moment to moment BOLD activation for real-time neurofeedback

PubMed Central

Hinds, Oliver; Ghosh, Satrajit; Thompson, Todd W.; Yoo, Julie J.; Whitfield-Gabrieli, Susan; Triantafyllou, Christina; Gabrieli, John D.E.

2013-01-01

Estimating moment to moment changes in blood oxygenation level dependent (BOLD) activation levels from functional magnetic resonance imaging (fMRI) data has applications for learned regulation of regional activation, brain state monitoring, and brain-machine interfaces. In each of these contexts, accurate estimation of the BOLD signal in as little time as possible is desired. This is a challenging problem due to the low signal-to-noise ratio of fMRI data. Previous methods for real-time fMRI analysis have either sacrificed the ability to compute moment to moment activation changes by averaging several acquisitions into a single activation estimate or have sacrificed accuracy by failing to account for prominent sources of noise in the fMRI signal. Here we present a new method for computing the amount of activation present in a single fMRI acquisition that separates moment to moment changes in the fMRI signal intensity attributable to neural sources from those due to noise, resulting in a feedback signal more reflective of neural activation. This method computes an incremental general linear model fit to the fMRI timeseries, which is used to calculate the expected signal intensity at each new acquisition. The difference between the measured intensity and the expected intensity is scaled by the variance of the estimator in order to transform this residual difference into a statistic. Both synthetic and real data were used to validate this method and compare it to the only other published real-time fMRI method. PMID:20682350

A robust two-node, 13 moment quadrature method of moments for dilute particle flows including wall bouncing

NASA Astrophysics Data System (ADS)

Sun, Dan; Garmory, Andrew; Page, Gary J.

2017-02-01

For flows where the particle number density is low and the Stokes number is relatively high, as found when sand or ice is ingested into aircraft gas turbine engines, streams of particles can cross each other's path or bounce from a solid surface without being influenced by inter-particle collisions. The aim of this work is to develop an Eulerian method to simulate these types of flow. To this end, a two-node quadrature-based moment method using 13 moments is proposed. In the proposed algorithm thirteen moments of particle velocity, including cross-moments of second order, are used to determine the weights and abscissas of the two nodes and to set up the association between the velocity components in each node. Previous Quadrature Method of Moments (QMOM) algorithms either use more than two nodes, leading to increased computational expense, or are shown here to give incorrect results under some circumstances. This method gives the computational efficiency advantages of only needing two particle phase velocity fields whilst ensuring that a correct combination of weights and abscissas is returned for any arbitrary combination of particle trajectories without the need for any further assumptions. Particle crossing and wall bouncing with arbitrary combinations of angles are demonstrated using the method in a two-dimensional scheme. The ability of the scheme to include the presence of drag from a carrier phase is also demonstrated, as is bouncing off surfaces with inelastic collisions. The method is also applied to the Taylor-Green vortex flow test case and is found to give results superior to the existing two-node QMOM method and is in good agreement with results from Lagrangian modelling of this case.

Study of photon correlation techniques for processing of laser velocimeter signals

NASA Technical Reports Server (NTRS)

Mayo, W. T., Jr.

1977-01-01

The objective was to provide the theory and a system design for a new type of photon counting processor for low level dual scatter laser velocimeter (LV) signals which would be capable of both the first order measurements of mean flow and turbulence intensity and also the second order time statistics: cross correlation auto correlation, and related spectra. A general Poisson process model for low level LV signals and noise which is valid from the photon-resolved regime all the way to the limiting case of nonstationary Gaussian noise was used. Computer simulation algorithms and higher order statistical moment analysis of Poisson processes were derived and applied to the analysis of photon correlation techniques. A system design using a unique dual correlate and subtract frequency discriminator technique is postulated and analyzed. Expectation analysis indicates that the objective measurements are feasible.

A hybrid multi-objective evolutionary algorithm for wind-turbine blade optimization

NASA Astrophysics Data System (ADS)

Sessarego, M.; Dixon, K. R.; Rival, D. E.; Wood, D. H.

2015-08-01

A concurrent-hybrid non-dominated sorting genetic algorithm (hybrid NSGA-II) has been developed and applied to the simultaneous optimization of the annual energy production, flapwise root-bending moment and mass of the NREL 5 MW wind-turbine blade. By hybridizing a multi-objective evolutionary algorithm (MOEA) with gradient-based local search, it is believed that the optimal set of blade designs could be achieved in lower computational cost than for a conventional MOEA. To measure the convergence between the hybrid and non-hybrid NSGA-II on a wind-turbine blade optimization problem, a computationally intensive case was performed using the non-hybrid NSGA-II. From this particular case, a three-dimensional surface representing the optimal trade-off between the annual energy production, flapwise root-bending moment and blade mass was achieved. The inclusion of local gradients in the blade optimization, however, shows no improvement in the convergence for this three-objective problem.

X-29A Lateral-Directional Stability and Control Derivatives Extracted From High-Angle-of-Attack Flight Data

NASA Technical Reports Server (NTRS)

Iliff, Kenneth W.; Wang, Kon-Sheng Charles Wang

1996-01-01

The lateral-directional stability and control derivatives of the X-29A number 2 are extracted from flight data over an angle-of-attack range of 4 degrees to 53 degrees using a parameter identification algorithm. The algorithm uses the linearized aircraft equations of motion and a maximum likelihood estimator in the presence of state and measurement noise. State noise is used to model the uncommanded forcing function caused by unsteady aerodynamics over the aircraft at angles of attack above 15 degrees. The results supported the flight-envelope-expansion phase of the X-29A number 2 by helping to update the aerodynamic mathematical model, to improve the real-time simulator, and to revise flight control system laws. Effects of the aircraft high gain flight control system on maneuver quality and the estimated derivatives are also discussed. The derivatives are plotted as functions of angle of attack and compared with the predicted aerodynamic database. Agreement between predicted and flight values is quite good for some derivatives such as the lateral force due to sideslip, the lateral force due to rudder deflection, and the rolling moment due to roll rate. The results also show significant differences in several important derivatives such as the rolling moment due to sideslip, the yawing moment due to sideslip, the yawing moment due to aileron deflection, and the yawing moment due to rudder deflection.

Comparison of PDF and Moment Closure Methods in the Modeling of Turbulent Reacting Flows

NASA Technical Reports Server (NTRS)

Norris, Andrew T.; Hsu, Andrew T.

1994-01-01

In modeling turbulent reactive flows, Probability Density Function (PDF) methods have an advantage over the more traditional moment closure schemes in that the PDF formulation treats the chemical reaction source terms exactly, while moment closure methods are required to model the mean reaction rate. The common model used is the laminar chemistry approximation, where the effects of turbulence on the reaction are assumed negligible. For flows with low turbulence levels and fast chemistry, the difference between the two methods can be expected to be small. However for flows with finite rate chemistry and high turbulence levels, significant errors can be expected in the moment closure method. In this paper, the ability of the PDF method and the moment closure scheme to accurately model a turbulent reacting flow is tested. To accomplish this, both schemes were used to model a CO/H2/N2- air piloted diffusion flame near extinction. Identical thermochemistry, turbulence models, initial conditions and boundary conditions are employed to ensure a consistent comparison can be made. The results of the two methods are compared to experimental data as well as to each other. The comparison reveals that the PDF method provides good agreement with the experimental data, while the moment closure scheme incorrectly shows a broad, laminar-like flame structure.

Gauge-free cluster variational method by maximal messages and moment matching.

PubMed

Domínguez, Eduardo; Lage-Castellanos, Alejandro; Mulet, Roberto; Ricci-Tersenghi, Federico

2017-04-01

We present an implementation of the cluster variational method (CVM) as a message passing algorithm. The kind of message passing algorithm used for CVM, usually named generalized belief propagation (GBP), is a generalization of the belief propagation algorithm in the same way that CVM is a generalization of the Bethe approximation for estimating the partition function. However, the connection between fixed points of GBP and the extremal points of the CVM free energy is usually not a one-to-one correspondence because of the existence of a gauge transformation involving the GBP messages. Our contribution is twofold. First, we propose a way of defining messages (fields) in a generic CVM approximation, such that messages arrive on a given region from all its ancestors, and not only from its direct parents, as in the standard parent-to-child GBP. We call this approach maximal messages. Second, we focus on the case of binary variables, reinterpreting the messages as fields enforcing the consistency between the moments of the local (marginal) probability distributions. We provide a precise rule to enforce all consistencies, avoiding any redundancy, that would otherwise lead to a gauge transformation on the messages. This moment matching method is gauge free, i.e., it guarantees that the resulting GBP is not gauge invariant. We apply our maximal messages and moment matching GBP to obtain an analytical expression for the critical temperature of the Ising model in general dimensions at the level of plaquette CVM. The values obtained outperform Bethe estimates, and are comparable with loop corrected belief propagation equations. The method allows for a straightforward generalization to disordered systems.

Gauge-free cluster variational method by maximal messages and moment matching

NASA Astrophysics Data System (ADS)

Domínguez, Eduardo; Lage-Castellanos, Alejandro; Mulet, Roberto; Ricci-Tersenghi, Federico

2017-04-01

We present an implementation of the cluster variational method (CVM) as a message passing algorithm. The kind of message passing algorithm used for CVM, usually named generalized belief propagation (GBP), is a generalization of the belief propagation algorithm in the same way that CVM is a generalization of the Bethe approximation for estimating the partition function. However, the connection between fixed points of GBP and the extremal points of the CVM free energy is usually not a one-to-one correspondence because of the existence of a gauge transformation involving the GBP messages. Our contribution is twofold. First, we propose a way of defining messages (fields) in a generic CVM approximation, such that messages arrive on a given region from all its ancestors, and not only from its direct parents, as in the standard parent-to-child GBP. We call this approach maximal messages. Second, we focus on the case of binary variables, reinterpreting the messages as fields enforcing the consistency between the moments of the local (marginal) probability distributions. We provide a precise rule to enforce all consistencies, avoiding any redundancy, that would otherwise lead to a gauge transformation on the messages. This moment matching method is gauge free, i.e., it guarantees that the resulting GBP is not gauge invariant. We apply our maximal messages and moment matching GBP to obtain an analytical expression for the critical temperature of the Ising model in general dimensions at the level of plaquette CVM. The values obtained outperform Bethe estimates, and are comparable with loop corrected belief propagation equations. The method allows for a straightforward generalization to disordered systems.

Three Dimensional Cross-Sectional Properties From Bone Densitometry

NASA Technical Reports Server (NTRS)

Cleek, Tammy M.; Whalen, Robert T.; Dalton, Bonnie P. (Technical Monitor)

2001-01-01

Bone densitometry has previously been used to obtain cross-sectional properties of bone in a single scan plane. Using three non-coplanar scans, we have extended the method to obtain the principal area Moments of inertia and orientations of the principal axes at each cross-section along the length of the scan. Various 5 aluminum phantoms were used to examine scanner characteristics to develop the highest accuracy possible for in vitro non-invasive analysis of mass distribution. Factors considered included X-ray photon energy, initial scan orientation, the included angle of the 3 scans, and Imin/Imax ratios. Principal moments of inertia were accurate to within 3.1% and principal angles were within 1 deg. of the expected value for phantoms scanned with included angles of 60 deg. and 90 deg. at the higher X-ray photon energy. Low standard deviations in error also 10 indicate high precision of calculated measurements with these included angles. Accuracy and precision decreased slightly when the included angle was reduced to 30 deg. The method was then successfully applied to a pair of excised cadaveric tibiae. The accuracy and insensitivity of the algorithms to cross-sectional shape and changing isotropy (Imin/Imax) values when various included angles are used make this technique viable for future in vivo studies.

Upper limb joint forces and moments during underwater cyclical movements.

PubMed

Lauer, Jessy; Rouard, Annie Hélène; Vilas-Boas, João Paulo

2016-10-03

Sound inverse dynamics modeling is lacking in aquatic locomotion research because of the difficulty in measuring hydrodynamic forces in dynamic conditions. Here we report the successful implementation and validation of an innovative methodology crossing new computational fluid dynamics and inverse dynamics techniques to quantify upper limb joint forces and moments while moving in water. Upper limb kinematics of seven male swimmers sculling while ballasted with 4kg was recorded through underwater motion capture. Together with body scans, segment inertial properties, and hydrodynamic resistances computed from a unique dynamic mesh algorithm capable to handle large body deformations, these data were fed into an inverse dynamics model to solve for joint kinetics. Simulation validity was assessed by comparing the impulse produced by the arms, calculated by integrating vertical forces over a stroke period, to the net theoretical impulse of buoyancy and ballast forces. A resulting gap of 1.2±3.5% provided confidence in the results. Upper limb joint load was within 5% of swimmer׳s body weight, which tends to supports the use of low-load aquatic exercises to reduce joint stress. We expect this significant methodological improvement to pave the way towards deeper insights into the mechanics of aquatic movement and the establishment of practice guidelines in rehabilitation, fitness or swimming performance. Copyright © 2016 Elsevier Ltd. All rights reserved.

SAMBA: Sparse Approximation of Moment-Based Arbitrary Polynomial Chaos

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ahlfeld, R., E-mail: r.ahlfeld14@imperial.ac.uk; Belkouchi, B.; Montomoli, F.

2016-09-01

A new arbitrary Polynomial Chaos (aPC) method is presented for moderately high-dimensional problems characterised by limited input data availability. The proposed methodology improves the algorithm of aPC and extends the method, that was previously only introduced as tensor product expansion, to moderately high-dimensional stochastic problems. The fundamental idea of aPC is to use the statistical moments of the input random variables to develop the polynomial chaos expansion. This approach provides the possibility to propagate continuous or discrete probability density functions and also histograms (data sets) as long as their moments exist, are finite and the determinant of the moment matrixmore » is strictly positive. For cases with limited data availability, this approach avoids bias and fitting errors caused by wrong assumptions. In this work, an alternative way to calculate the aPC is suggested, which provides the optimal polynomials, Gaussian quadrature collocation points and weights from the moments using only a handful of matrix operations on the Hankel matrix of moments. It can therefore be implemented without requiring prior knowledge about statistical data analysis or a detailed understanding of the mathematics of polynomial chaos expansions. The extension to more input variables suggested in this work, is an anisotropic and adaptive version of Smolyak's algorithm that is solely based on the moments of the input probability distributions. It is referred to as SAMBA (PC), which is short for Sparse Approximation of Moment-Based Arbitrary Polynomial Chaos. It is illustrated that for moderately high-dimensional problems (up to 20 different input variables or histograms) SAMBA can significantly simplify the calculation of sparse Gaussian quadrature rules. SAMBA's efficiency for multivariate functions with regard to data availability is further demonstrated by analysing higher order convergence and accuracy for a set of nonlinear test functions with 2, 5 and 10 different input distributions or histograms.« less

SAMBA: Sparse Approximation of Moment-Based Arbitrary Polynomial Chaos

NASA Astrophysics Data System (ADS)

Ahlfeld, R.; Belkouchi, B.; Montomoli, F.

2016-09-01

A new arbitrary Polynomial Chaos (aPC) method is presented for moderately high-dimensional problems characterised by limited input data availability. The proposed methodology improves the algorithm of aPC and extends the method, that was previously only introduced as tensor product expansion, to moderately high-dimensional stochastic problems. The fundamental idea of aPC is to use the statistical moments of the input random variables to develop the polynomial chaos expansion. This approach provides the possibility to propagate continuous or discrete probability density functions and also histograms (data sets) as long as their moments exist, are finite and the determinant of the moment matrix is strictly positive. For cases with limited data availability, this approach avoids bias and fitting errors caused by wrong assumptions. In this work, an alternative way to calculate the aPC is suggested, which provides the optimal polynomials, Gaussian quadrature collocation points and weights from the moments using only a handful of matrix operations on the Hankel matrix of moments. It can therefore be implemented without requiring prior knowledge about statistical data analysis or a detailed understanding of the mathematics of polynomial chaos expansions. The extension to more input variables suggested in this work, is an anisotropic and adaptive version of Smolyak's algorithm that is solely based on the moments of the input probability distributions. It is referred to as SAMBA (PC), which is short for Sparse Approximation of Moment-Based Arbitrary Polynomial Chaos. It is illustrated that for moderately high-dimensional problems (up to 20 different input variables or histograms) SAMBA can significantly simplify the calculation of sparse Gaussian quadrature rules. SAMBA's efficiency for multivariate functions with regard to data availability is further demonstrated by analysing higher order convergence and accuracy for a set of nonlinear test functions with 2, 5 and 10 different input distributions or histograms.

Magnetic moment of solar plasma and the Kelvin force: -The driving force of plasma up-flow -

NASA Astrophysics Data System (ADS)

Shibasaki, Kiyoto

2017-04-01

Thermal plasma in the solar atmosphere is magnetized (diamagnetic). The magnetic moment does not disappear by collisions because complete gyration is not a necessary condition to have magnetic moment. Magnetized fluid is subjected to Kelvin force in non-uniform magnetic field. Generally, magnetic field strength decreases upwards in the solar atmosphere, hence the Kelvin force is directed upwards along the field. This force is not included in the fluid treatment of MHD. By adding the Kelvin force to the MHD equation of motion, we can expect temperature dependent plasma flows along the field which are reported by many observations. The temperature dependence of the flow speed is explained by temperature dependence of magnetic moment. From the observed parameters, we can infer physical parameters in the solar atmosphere such as scale length of the magnetic field strength and the friction force acting on the flowing plasma. In case of closed magnetic field lines, loop-top concentration of hot plasma is expected which is frequently observed.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lau, Sarah J.; Moore, David G.; Stair, Sarah L.

Ultrasonic analysis is being explored as a way to capture events during melting of highly dispersive wax. Typical events include temperature changes in the material, phase transition of the material, surface flows and reformations, and void filling as the material melts. Melt tests are performed with wax to evaluate the usefulness of different signal processing algorithms in capturing event data. Several algorithm paths are being pursued. The first looks at changes in the velocity of the signal through the material. This is only appropriate when the changes from one ultrasonic signal to the next can be represented by a linearmore » relationship, which is not always the case. The second tracks changes in the frequency content of the signal. The third algorithm tracks changes in the temporal moments of a signal over a full test. This method does not require that the changes in the signal be represented by a linear relationship, but attaching changes in the temporal moments to physical events can be difficult. This study describes the algorithm paths applied to experimental data from ultrasonic signals as wax melts and explores different ways to display the results.« less

Algorithm for Surface of Translation Attached Radiators (A-STAR). Volume 2. Users manual

NASA Astrophysics Data System (ADS)

Medgyesimitschang, L. N.; Putnam, J. M.

1982-05-01

A hierarchy of computer programs implementing the method of moments for bodies of translation (MM/BOT) is described. The algorithm treats the far-field radiation from off-surface and aperture antennas on finite-length open or closed bodies of arbitrary cross section. The near fields and antenna coupling on such bodies are computed. The theoretical development underlying the algorithm is described in Volume 1 of this report.

Accelerating image recognition on mobile devices using GPGPU

NASA Astrophysics Data System (ADS)

Bordallo López, Miguel; Nykänen, Henri; Hannuksela, Jari; Silvén, Olli; Vehviläinen, Markku

2011-01-01

The future multi-modal user interfaces of battery-powered mobile devices are expected to require computationally costly image analysis techniques. The use of Graphic Processing Units for computing is very well suited for parallel processing and the addition of programmable stages and high precision arithmetic provide for opportunities to implement energy-efficient complete algorithms. At the moment the first mobile graphics accelerators with programmable pipelines are available, enabling the GPGPU implementation of several image processing algorithms. In this context, we consider a face tracking approach that uses efficient gray-scale invariant texture features and boosting. The solution is based on the Local Binary Pattern (LBP) features and makes use of the GPU on the pre-processing and feature extraction phase. We have implemented a series of image processing techniques in the shader language of OpenGL ES 2.0, compiled them for a mobile graphics processing unit and performed tests on a mobile application processor platform (OMAP3530). In our contribution, we describe the challenges of designing on a mobile platform, present the performance achieved and provide measurement results for the actual power consumption in comparison to using the CPU (ARM) on the same platform.

Pattern Recognition Control Design

NASA Technical Reports Server (NTRS)

Gambone, Elisabeth A.

2018-01-01

Spacecraft control algorithms must know the expected vehicle response to any command to the available control effectors, such as reaction thrusters or torque devices. Spacecraft control system design approaches have traditionally relied on the estimated vehicle mass properties to determine the desired force and moment, as well as knowledge of the effector performance to efficiently control the spacecraft. A pattern recognition approach was used to investigate the relationship between the control effector commands and spacecraft responses. Instead of supplying the approximated vehicle properties and the thruster performance characteristics, a database of information relating the thruster ring commands and the desired vehicle response was used for closed-loop control. A Monte Carlo simulation data set of the spacecraft dynamic response to effector commands was analyzed to establish the influence a command has on the behavior of the spacecraft. A tool developed at NASA Johnson Space Center to analyze flight dynamics Monte Carlo data sets through pattern recognition methods was used to perform this analysis. Once a comprehensive data set relating spacecraft responses with commands was established, it was used in place of traditional control methods and gains set. This pattern recognition approach was compared with traditional control algorithms to determine the potential benefits and uses.

Limiting neutrino magnetic moments with Borexino Phase-II solar neutrino data

NASA Astrophysics Data System (ADS)

Agostini, M.; Altenmüller, K.; Appel, S.; Atroshchenko, V.; Bagdasarian, Z.; Basilico, D.; Bellini, G.; Benziger, J.; Bick, D.; Bonfini, G.; Bravo, D.; Caccianiga, B.; Calaprice, F.; Caminata, A.; Caprioli, S.; Carlini, M.; Cavalcante, P.; Chepurnov, A.; Choi, K.; Collica, L.; D'Angelo, D.; Davini, S.; Derbin, A.; Ding, X. F.; Di Ludovico, A.; Di Noto, L.; Drachnev, I.; Fomenko, K.; Formozov, A.; Franco, D.; Froborg, F.; Gabriele, F.; Galbiati, C.; Ghiano, C.; Giammarchi, M.; Goretti, A.; Gromov, M.; Guffanti, D.; Hagner, C.; Houdy, T.; Hungerford, E.; Ianni, Aldo; Ianni, Andrea; Jany, A.; Jeschke, D.; Kobychev, V.; Korablev, D.; Korga, G.; Kryn, D.; Laubenstein, M.; Litvinovich, E.; Lombardi, F.; Lombardi, P.; Ludhova, L.; Lukyanchenko, G.; Lukyanchenko, L.; Machulin, I.; Manuzio, G.; Marcocci, S.; Martyn, J.; Meroni, E.; Meyer, M.; Miramonti, L.; Misiaszek, M.; Muratova, V.; Neumair, B.; Oberauer, L.; Opitz, B.; Orekhov, V.; Ortica, F.; Pallavicini, M.; Papp, L.; Penek, Ã.-.; Pilipenko, N.; Pocar, A.; Porcelli, A.; Ranucci, G.; Razeto, A.; Re, A.; Redchuk, M.; Romani, A.; Roncin, R.; Rossi, N.; Schönert, S.; Semenov, D.; Skorokhvatov, M.; Smirnov, O.; Sotnikov, A.; Stokes, L. F. F.; Suvorov, Y.; Tartaglia, R.; Testera, G.; Thurn, J.; Toropova, M.; Unzhakov, E.; Vishneva, A.; Vogelaar, R. B.; von Feilitzsch, F.; Wang, H.; Weinz, S.; Wojcik, M.; Wurm, M.; Yokley, Z.; Zaimidoroga, O.; Zavatarelli, S.; Zuber, K.; Zuzel, G.; Borexino Collaboration

2017-11-01

A search for the solar neutrino effective magnetic moment has been performed using data from 1291.5 days exposure during the second phase of the Borexino experiment. No significant deviations from the expected shape of the electron recoil spectrum from solar neutrinos have been found, and a new upper limit on the effective neutrino magnetic moment of μνeff<2.8×10 -11 μB at 90% C.L. has been set using constraints on the sum of the solar neutrino fluxes implied by the radiochemical gallium experiments. Using the limit for the effective neutrino moment, new limits for the magnetic moments of the neutrino flavor states, and for the elements of the neutrino magnetic moments matrix for Dirac and Majorana neutrinos, are derived.

Assessment of macroseismic intensity in the Nile basin, Egypt

NASA Astrophysics Data System (ADS)

Fergany, Elsayed

2018-01-01

This work intends to assess deterministic seismic hazard and risk analysis in terms of the maximum expected intensity map of the Egyptian Nile basin sector. Seismic source zone model of Egypt was delineated based on updated compatible earthquake catalog in 2015, focal mechanisms, and the common tectonic elements. Four effective seismic source zones were identified along the Nile basin. The observed macroseismic intensity data along the basin was used to develop intensity prediction equation defined in terms of moment magnitude. Expected maximum intensity map was proven based on the developed intensity prediction equation, identified effective seismic source zones, and maximum expected magnitude for each zone along the basin. The earthquake hazard and risk analysis was discussed and analyzed in view of the maximum expected moment magnitude and the maximum expected intensity values for each effective source zone. Moderate expected magnitudes are expected to put high risk at Cairo and Aswan regions. The results of this study could be a recommendation for the planners in charge to mitigate the seismic risk at these strategic zones of Egypt.

Magnetic Moment of Proton Drip-Line Nucleus (9)C

NASA Technical Reports Server (NTRS)

Matsuta, K.; Fukuda, M.; Tanigaki, M.; Minamisono, T.; Nojiri, Y.; Mihara, M.; Onishi, T.; Yamaguchi, T.; Harada, A.; Sasaki, M.

1994-01-01

The magnetic moment of the proton drip-line nucleus C-9(I(sup (pi)) = 3/2, T(sub 1/2) = 126 ms) has been measured for the first time, using the beta-NMR detection technique with polarized radioactive beams. The measure value for the magnetic moment is 1mu(C-9)! = 1.3914 +/- 0.0005 (mu)N. The deduced spin expectation value of 1.44 is unusually larger than any other ones of even-odd nuclei.

A generalized Grubbs-Beck test statistic for detecting multiple potentially influential low outliers in flood series

USGS Publications Warehouse

Cohn, T.A.; England, J.F.; Berenbrock, C.E.; Mason, R.R.; Stedinger, J.R.; Lamontagne, J.R.

2013-01-01

he Grubbs-Beck test is recommended by the federal guidelines for detection of low outliers in flood flow frequency computation in the United States. This paper presents a generalization of the Grubbs-Beck test for normal data (similar to the Rosner (1983) test; see also Spencer and McCuen (1996)) that can provide a consistent standard for identifying multiple potentially influential low flows. In cases where low outliers have been identified, they can be represented as “less-than” values, and a frequency distribution can be developed using censored-data statistical techniques, such as the Expected Moments Algorithm. This approach can improve the fit of the right-hand tail of a frequency distribution and provide protection from lack-of-fit due to unimportant but potentially influential low flows (PILFs) in a flood series, thus making the flood frequency analysis procedure more robust.

A generalized Grubbs-Beck test statistic for detecting multiple potentially influential low outliers in flood series

NASA Astrophysics Data System (ADS)

Cohn, T. A.; England, J. F.; Berenbrock, C. E.; Mason, R. R.; Stedinger, J. R.; Lamontagne, J. R.

2013-08-01

The Grubbs-Beck test is recommended by the federal guidelines for detection of low outliers in flood flow frequency computation in the United States. This paper presents a generalization of the Grubbs-Beck test for normal data (similar to the Rosner (1983) test; see also Spencer and McCuen (1996)) that can provide a consistent standard for identifying multiple potentially influential low flows. In cases where low outliers have been identified, they can be represented as "less-than" values, and a frequency distribution can be developed using censored-data statistical techniques, such as the Expected Moments Algorithm. This approach can improve the fit of the right-hand tail of a frequency distribution and provide protection from lack-of-fit due to unimportant but potentially influential low flows (PILFs) in a flood series, thus making the flood frequency analysis procedure more robust.

The Fermilab Muon g-2 experiment: laser calibration system

DOE PAGES

Karuza, M.; Anastasi, A.; Basti, A.; ...

2017-08-17

The anomalous muon dipole magnetic moment can be measured (and calculated) with great precision thus providing insight on the Standard Model and new physics. Currently an experiment is under construction at Fermilab (U.S.A.) which is expected to measure the anomalous muon dipole magnetic moment with unprecedented precision. One of the improvements with respect to the previous experiments is expected to come from the laser calibration system which has been designed and constructed by the Italian part of the collaboration (INFN). Furthermore, an emphasis of this paper will be on the calibration system that is in the final stages of constructionmore » as well as the experiment which is expected to start data taking this year.« less

Convergence of moment expansions for expectation values with embedded random matrix ensembles and quantum chaos

NASA Astrophysics Data System (ADS)

Kota, V. K. B.

2003-07-01

Smoothed forms for expectation values < K> E of positive definite operators K follow from the K-density moments either directly or in many other ways each giving a series expansion (involving polynomials in E). In large spectroscopic spaces one has to partition the many particle spaces into subspaces. Partitioning leads to new expansions for expectation values. It is shown that all the expansions converge to compact forms depending on the nature of the operator K and the operation of embedded random matrix ensembles and quantum chaos in many particle spaces. Explicit results are given for occupancies < ni> E, spin-cutoff factors < JZ2> E and strength sums < O†O> E, where O is a one-body transition operator.

Moment measurements in dynamic and quasi-static spine segment testing using eccentric compression are susceptible to artifacts based on loading configuration.

PubMed

Van Toen, Carolyn; Carter, Jarrod W; Oxland, Thomas R; Cripton, Peter A

2014-12-01

The tolerance of the spine to bending moments, used for evaluation of injury prevention devices, is often determined through eccentric axial compression experiments using segments of the cadaver spine. Preliminary experiments in our laboratory demonstrated that eccentric axial compression resulted in "unexpected" (artifact) moments. The aim of this study was to evaluate the static and dynamic effects of test configuration on bending moments during eccentric axial compression typical in cadaver spine segment testing. Specific objectives were to create dynamic equilibrium equations for the loads measured inferior to the specimen, experimentally verify these equations, and compare moment responses from various test configurations using synthetic (rubber) and human cadaver specimens. The equilibrium equations were verified by performing quasi-static (5 mm/s) and dynamic experiments (0.4 m/s) on a rubber specimen and comparing calculated shear forces and bending moments to those measured using a six-axis load cell. Moment responses were compared for hinge joint, linear slider and hinge joint, and roller joint configurations tested at quasi-static and dynamic rates. Calculated shear force and bending moment curves had similar shapes to those measured. Calculated values in the first local minima differed from those measured by 3% and 15%, respectively, in the dynamic test, and these occurred within 1.5 ms of those measured. In the rubber specimen experiments, for the hinge joint (translation constrained), quasi-static and dynamic posterior eccentric compression resulted in flexion (unexpected) moments. For the slider and hinge joints and the roller joints (translation unconstrained), extension ("expected") moments were measured quasi-statically and initial flexion (unexpected) moments were measured dynamically. In the cadaver experiments with roller joints, anterior and posterior eccentricities resulted in extension moments, which were unexpected and expected, for those configurations, respectively. The unexpected moments were due to the inertia of the superior mounting structures. This study has shown that eccentric axial compression produces unexpected moments due to translation constraints at all loading rates and due to the inertia of the superior mounting structures in dynamic experiments. It may be incorrect to assume that bending moments are equal to the product of compression force and eccentricity, particularly where the test configuration involves translational constraints and where the experiments are dynamic. In order to reduce inertial moment artifacts, the mass, and moment of inertia of any loading jig structures that rotate with the specimen should be minimized. Also, the distance between these structures and the load cell should be reduced.

Estimation of full moment tensors, including uncertainties, for earthquakes, volcanic events, and nuclear explosions

NASA Astrophysics Data System (ADS)

Alvizuri, Celso R.

We present a catalog of full seismic moment tensors for 63 events from Uturuncu volcano in Bolivia. The events were recorded during 2011-2012 in the PLUTONS seismic array of 24 broadband stations. Most events had magnitudes between 0.5 and 2.0 and did not generate discernible surface waves; the largest event was Mw 2.8. For each event we computed the misfit between observed and synthetic waveforms, and we used first-motion polarity measurements to reduce the number of possible solutions. Each moment tensor solution was obtained using a grid search over the six-dimensional space of moment tensors. For each event we show the misfit function in eigenvalue space, represented by a lune. We identify three subsets of the catalog: (1) 6 isotropic events, (2) 5 tensional crack events, and (3) a swarm of 14 events southeast of the volcanic center that appear to be double couples. The occurrence of positively isotropic events is consistent with other published results from volcanic and geothermal regions. Several of these previous results, as well as our results, cannot be interpreted within the context of either an oblique opening crack or a crack-plus-double-couple model. Proper characterization of uncertainties for full moment tensors is critical for distinguishing among physical models of source processes. A seismic moment tensor is a 3x3 symmetric matrix that provides a compact representation of a seismic source. We develop an algorithm to estimate moment tensors and their uncertainties from observed seismic data. For a given event, the algorithm performs a grid search over the six-dimensional space of moment tensors by generating synthetic waveforms for each moment tensor and then evaluating a misfit function between the observed and synthetic waveforms. 'The' moment tensor M0 for the event is then the moment tensor with minimum misfit. To describe the uncertainty associated with M0, we first convert the misfit function to a probability function. The uncertainty, or rather the confidence, is then given by the 'confidence curve' P( V), where P(V) is the probability that the true moment tensor for the event lies within the neighborhood of M that has fractional volume V. The area under the confidence curve provides a single, abbreviated 'confidence parameter' for M0. We apply the method to data from events in different regions and tectonic settings: 63 small (M w 4) earthquakes in the southern Alaska subduction zone, and 12 earthquakes and 17 nuclear explosions at the Nevada Test Site. Characterization of moment tensor uncertainties puts us in better position to discriminate among moment tensor source types and to assign physical processes to the events.

Prediction of seismic collapse risk of steel moment frame mid-rise structures by meta-heuristic algorithms

NASA Astrophysics Data System (ADS)

Jough, Fooad Karimi Ghaleh; Şensoy, Serhan

2016-12-01

Different performance levels may be obtained for sideway collapse evaluation of steel moment frames depending on the evaluation procedure used to handle uncertainties. In this article, the process of representing modelling uncertainties, record to record (RTR) variations and cognitive uncertainties for moment resisting steel frames of various heights is discussed in detail. RTR uncertainty is used by incremental dynamic analysis (IDA), modelling uncertainties are considered through backbone curves and hysteresis loops of component, and cognitive uncertainty is presented in three levels of material quality. IDA is used to evaluate RTR uncertainty based on strong ground motion records selected by the k-means algorithm, which is favoured over Monte Carlo selection due to its time saving appeal. Analytical equations of the Response Surface Method are obtained through IDA results by the Cuckoo algorithm, which predicts the mean and standard deviation of the collapse fragility curve. The Takagi-Sugeno-Kang model is used to represent material quality based on the response surface coefficients. Finally, collapse fragility curves with the various sources of uncertainties mentioned are derived through a large number of material quality values and meta variables inferred by the Takagi-Sugeno-Kang fuzzy model based on response surface method coefficients. It is concluded that a better risk management strategy in countries where material quality control is weak, is to account for cognitive uncertainties in fragility curves and the mean annual frequency.

Determination of the diffusivity, dispersion, skewness and kurtosis in heterogeneous porous flow. Part I: Analytical solutions with the extended method of moments.

NASA Astrophysics Data System (ADS)

Ginzburg, Irina; Vikhansky, Alexander

2018-05-01

The extended method of moments (EMM) is elaborated in recursive algorithmic form for the prediction of the effective diffusivity, the Taylor dispersion dyadic and the associated longitudinal high-order coefficients in mean-concentration profiles and residence-time distributions. The method applies in any streamwise-periodic stationary d-dimensional velocity field resolved in the piecewise continuous heterogeneous porosity field. It is demonstrated that EMM reduces to the method of moments and the volume-averaging formulation in microscopic velocity field and homogeneous soil, respectively. The EMM simultaneously constructs two systems of moments, the spatial and the temporal, without resorting to solving of the high-order upscaled PDE. At the same time, the EMM is supported with the reconstruction of distribution from its moments, allowing to visualize the deviation from the classical ADE solution. The EMM can be handled by any linear advection-diffusion solver with explicit mass-source and diffusive-flux jump condition on the solid boundary and permeable interface. The prediction of the first four moments is decisive in the optimization of the dispersion, asymmetry, peakedness and heavy-tails of the solute distributions, through an adequate design of the composite materials, wetlands, chemical devices or oil recovery. The symbolic solutions for dispersion, skewness and kurtosis are constructed in basic configurations: diffusion process and Darcy flow through two porous blocks in "series", straight and radial Poiseuille flow, porous flow governed by the Stokes-Brinkman-Darcy channel equation and a fracture surrounded by penetrable diffusive matrix or embedded in porous flow. We examine the moments dependency upon porosity contrast, aspect ratio, Péclet and Darcy numbers, but also for their response on the effective Brinkman viscosity applied in flow modeling. Two numerical Lattice Boltzmann algorithms, a direct solver of the microscopic ADE in heterogeneous structure and a novel scheme for EMM numerical formulation, are called for validation of the constructed analytical predictions.

Rapid Characterization of Magnetic Moment of Cells for Magnetic Separation

PubMed Central

Ooi, Chinchun; Earhart, Christopher M.; Wilson, Robert J.; Wang, Shan X.

2014-01-01

NCI-H1650 lung cancer cell lines labeled with magnetic nanoparticles via the Epithelial Cell Adhesion Molecule (EpCAM) antigen were previously shown to be captured at high efficiencies by a microfabricated magnetic sifter. If fine control and optimization of the magnetic separation process is to be achieved, it is vital to be able to characterize the labeled cells’ magnetic moment rapidly. We have thus adapted a rapid prototyping method to obtain the saturation magnetic moment of these cells. This method utilizes a cross-correlation algorithm to analyze the cells’ motion in a simple fluidic channel to obtain their magnetophoretic velocity, and is effective even when the magnetic moments of cells are small. This rapid characterization is proven useful in optimizing our microfabricated magnetic sifter procedures for magnetic cell capture. PMID:24771946

Hamiltonian approach to Ehrenfest expectation values and Gaussian quantum states

PubMed Central

Bonet-Luz, Esther

2016-01-01

The dynamics of quantum expectation values is considered in a geometric setting. First, expectation values of the canonical observables are shown to be equivariant momentum maps for the action of the Heisenberg group on quantum states. Then, the Hamiltonian structure of Ehrenfest’s theorem is shown to be Lie–Poisson for a semidirect-product Lie group, named the Ehrenfest group. The underlying Poisson structure produces classical and quantum mechanics as special limit cases. In addition, quantum dynamics is expressed in the frame of the expectation values, in which the latter undergo canonical Hamiltonian motion. In the case of Gaussian states, expectation values dynamics couples to second-order moments, which also enjoy a momentum map structure. Eventually, Gaussian states are shown to possess a Lie–Poisson structure associated with another semidirect-product group, which is called the Jacobi group. This structure produces the energy-conserving variant of a class of Gaussian moment models that have previously appeared in the chemical physics literature. PMID:27279764

Implementation of parallel moment equations in NIMROD

NASA Astrophysics Data System (ADS)

Lee, Hankyu Q.; Held, Eric D.; Ji, Jeong-Young

2017-10-01

As collisionality is low (the Knudsen number is large) in many plasma applications, kinetic effects become important, particularly in parallel dynamics for magnetized plasmas. Fluid models can capture some kinetic effects when integral parallel closures are adopted. The adiabatic and linear approximations are used in solving general moment equations to obtain the integral closures. In this work, we present an effort to incorporate non-adiabatic (time-dependent) and nonlinear effects into parallel closures. Instead of analytically solving the approximate moment system, we implement exact parallel moment equations in the NIMROD fluid code. The moment code is expected to provide a natural convergence scheme by increasing the number of moments. Work in collaboration with the PSI Center and supported by the U.S. DOE under Grant Nos. DE-SC0014033, DE-SC0016256, and DE-FG02-04ER54746.

Improved limit on the Ra 225 electric dipole moment

DOE PAGES

Bishof, Michael; Parker, Richard H.; Bailey, Kevin G.; ...

2016-08-03

In this study, octupole-deformed nuclei, such as that of 225Ra, are expected to amplify observable atomic electric dipole moments (EDMs) that arise from time-reversal and parity-violating interactions in the nuclear medium. In 2015 we reported the first “proof-of-principle” measurement of the 225Ra atomic EDM.

Improved limit on the Ra 225 electric dipole moment

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bishof, Michael; Parker, Richard H.; Bailey, Kevin G.

In this study, octupole-deformed nuclei, such as that of 225Ra, are expected to amplify observable atomic electric dipole moments (EDMs) that arise from time-reversal and parity-violating interactions in the nuclear medium. In 2015 we reported the first “proof-of-principle” measurement of the 225Ra atomic EDM.

High-frequency, transient magnetic susceptibility of ferroelectrics

NASA Astrophysics Data System (ADS)

Grimes, Craig A.

1996-10-01

A significant high-frequency magnetic susceptibility was measured both in weakly polarized and nonpolarized samples of barium titanate, lead zirconate titanate, and carnauba wax. Magnetic susceptibility measurements were made from 10 to 500 MHz using a thin film permeameter at room temperature; initial susceptibilities ranged from 0.1 to 2.5. These values are larger than expected for paramagnets and smaller than expected for ferromagnets. It was found that the magnetic susceptibility decreases rapidly with exposure to the exciting field. The origin of the magnetic susceptibility is thought to originate with the applied time varying electric field associated with the susceptibility measurements. An electric field acts to rotate an electric dipole, creating a magnetic quadrupole if the two moments are balanced, and a net magnetic dipole moment if imbalanced. It is thought that local electrostatic fields created at ferroelectric domain discontinuities associated with grain boundaries create an imbalance in the anion rotation that results in a net, measurable, magnetic moment. The origin of the magnetic aftereffect may be due to the local heating of the material through the moving charges associated with the magnetic moment.

Cross-sectional structural parameters from densitometry

NASA Technical Reports Server (NTRS)

Cleek, Tammy M.; Whalen, Robert T.

2002-01-01

Bone densitometry has previously been used to obtain cross-sectional properties of bone from a single X-ray projection across the bone width. Using three unique projections, we have extended the method to obtain the principal area moments of inertia and orientations of the principal axes at each scan cross-section along the length of the scan. Various aluminum phantoms were used to examine scanner characteristics to develop the highest accuracy possible for in vitro non-invasive analysis of cross-sectional properties. Factors considered included X-ray photon energy, initial scan orientation, the angle spanned by the three scans (included angle), and I(min)/I(max) ratios. Principal moments of inertia were accurate to within +/-3.1% and principal angles were within +/-1 degrees of the expected value for phantoms scanned with included angles of 60 degrees and 90 degrees at the higher X-ray photon energy (140 kVp). Low standard deviations in the error (0.68-1.84%) also indicate high precision of calculated measurements with these included angles. Accuracy and precision decreased slightly when the included angle was reduced to 30 degrees. The method was then successfully applied to a pair of excised cadaveric tibiae. The accuracy and insensitivity of the algorithms to cross-sectional shape and changing isotropy (I(min)/I(max)) values when various included angles are used make this technique viable for future in vivo studies.

Studies of the 4-JET Rate and of Moments of Event Shape Observables Using Jade Data

NASA Astrophysics Data System (ADS)

Kluth, S.

2005-04-01

Data from e+e- annihilation into hadrons collected by the JADE experiment at centre-of-mass energies between 14 and 44 GeV were used to study the 4-jet rate using the Durham algorithm as well as the first five moments of event shape observables. The data were compared with NLO QCD predictions, augmented by resummed NLLA calculations for the 4-jet rate, in order to extract values of the strong coupling constant αS. The preliminary results are αS(MZ0) = 0.1169 ± 0.0026 (4-jet rate) and αS(MZ0) = 0.1286 ± 0.0072 (moments) consistent with the world average value. For some of the higher moments systematic deficiencies of the QCD predictions are observed.

Noise-enhanced clustering and competitive learning algorithms.

PubMed

Osoba, Osonde; Kosko, Bart

2013-01-01

Noise can provably speed up convergence in many centroid-based clustering algorithms. This includes the popular k-means clustering algorithm. The clustering noise benefit follows from the general noise benefit for the expectation-maximization algorithm because many clustering algorithms are special cases of the expectation-maximization algorithm. Simulations show that noise also speeds up convergence in stochastic unsupervised competitive learning, supervised competitive learning, and differential competitive learning. Copyright © 2012 Elsevier Ltd. All rights reserved.

Hierarchical trie packet classification algorithm based on expectation-maximization clustering.

PubMed

Bi, Xia-An; Zhao, Junxia

2017-01-01

With the development of computer network bandwidth, packet classification algorithms which are able to deal with large-scale rule sets are in urgent need. Among the existing algorithms, researches on packet classification algorithms based on hierarchical trie have become an important packet classification research branch because of their widely practical use. Although hierarchical trie is beneficial to save large storage space, it has several shortcomings such as the existence of backtracking and empty nodes. This paper proposes a new packet classification algorithm, Hierarchical Trie Algorithm Based on Expectation-Maximization Clustering (HTEMC). Firstly, this paper uses the formalization method to deal with the packet classification problem by means of mapping the rules and data packets into a two-dimensional space. Secondly, this paper uses expectation-maximization algorithm to cluster the rules based on their aggregate characteristics, and thereby diversified clusters are formed. Thirdly, this paper proposes a hierarchical trie based on the results of expectation-maximization clustering. Finally, this paper respectively conducts simulation experiments and real-environment experiments to compare the performances of our algorithm with other typical algorithms, and analyzes the results of the experiments. The hierarchical trie structure in our algorithm not only adopts trie path compression to eliminate backtracking, but also solves the problem of low efficiency of trie updates, which greatly improves the performance of the algorithm.

Vehicle handling and stability control by the cooperative control of 4WS and DYC

NASA Astrophysics Data System (ADS)

Shen, Huan; Tan, Yun-Sheng

2017-07-01

This paper proposes an integrated control system that cooperates with the four-wheel steering (4WS) and direct yaw moment control (DYC) to improve the vehicle handling and stability. The design works of the four-wheel steering and DYC control are based on sliding mode control. The integration control system produces the suitable 4WS angle and corrective yaw moment so that the vehicle tracks the desired yaw rate and sideslip angle. Considering the change of the vehicle longitudinal velocity that means the comfort of driving conditions, both the driving torque and braking torque are used to generate the corrective yaw moment. Simulation results show the effectiveness of the proposed control algorithm.

Determination of ground and excited state dipole moments via electronic Stark spectroscopy: 5-methoxyindole

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wilke, Josefin; Wilke, Martin; Schmitt, Michael, E-mail: mschmitt@uni-duesseldorf.de

2016-01-28

The dipole moments of the ground and lowest electronically excited singlet state of 5-methoxyindole have been determined by means of optical Stark spectroscopy in a molecular beam. The resulting spectra arise from a superposition of different field configurations, one with the static electric field almost parallel to the polarization of the exciting laser radiation, the other nearly perpendicular. Each field configuration leads to different intensities in the rovibronic spectrum. With an automated evolutionary algorithm approach, the spectra can be fit and the ratio of both field configurations can be determined. A simultaneous fit of two spectra with both field configurationsmore » improved the precision of the dipole moment determination by a factor of two. We find a reduction of the absolute dipole moment from 1.59(3) D to 1.14(6) D upon electronic excitation to the lowest electronically excited singlet state. At the same time, the dipole moment orientation rotates by 54{sup ∘} showing the importance of the determination of the dipole moment components. The dipole moment in the electronic ground state can approximately be obtained from a vector addition of the indole and the methoxy group dipole moments. However, in the electronically excited state, vector addition completely fails to describe the observed dipole moment. Several reasons for this behavior are discussed.« less

A novel vehicle dynamics stability control algorithm based on the hierarchical strategy with constrain of nonlinear tyre forces

NASA Astrophysics Data System (ADS)

Li, Liang; Jia, Gang; Chen, Jie; Zhu, Hongjun; Cao, Dongpu; Song, Jian

2015-08-01

Direct yaw moment control (DYC), which differentially brakes the wheels to produce a yaw moment for the vehicle stability in a steering process, is an important part of electric stability control system. In this field, most control methods utilise the active brake pressure with a feedback controller to adjust the braked wheel. However, the method might lead to a control delay or overshoot because of the lack of a quantitative project relationship between target values from the upper stability controller to the lower pressure controller. Meanwhile, the stability controller usually ignores the implementing ability of the tyre forces, which might be restrained by the combined-slip dynamics of the tyre. Therefore, a novel control algorithm of DYC based on the hierarchical control strategy is brought forward in this paper. As for the upper controller, a correctional linear quadratic regulator, which not only contains feedback control but also contains feed forward control, is introduced to deduce the object of the stability yaw moment in order to guarantee the yaw rate and side-slip angle stability. As for the medium and lower controller, the quantitative relationship between the vehicle stability object and the target tyre forces of controlled wheels is proposed to achieve smooth control performance based on a combined-slip tyre model. The simulations with the hardware-in-the-loop platform validate that the proposed algorithm can improve the stability of the vehicle effectively.

Dynamic imaging in electrical impedance tomography of the human chest with online transition matrix identification.

PubMed

Moura, Fernando Silva; Aya, Julio Cesar Ceballos; Fleury, Agenor Toledo; Amato, Marcelo Britto Passos; Lima, Raul Gonzalez

2010-02-01

One of the electrical impedance tomography objectives is to estimate the electrical resistivity distribution in a domain based only on electrical potential measurements at its boundary generated by an imposed electrical current distribution into the boundary. One of the methods used in dynamic estimation is the Kalman filter. In biomedical applications, the random walk model is frequently used as evolution model and, under this conditions, poor tracking ability of the extended Kalman filter (EKF) is achieved. An analytically developed evolution model is not feasible at this moment. The paper investigates the identification of the evolution model in parallel to the EKF and updating the evolution model with certain periodicity. The evolution model transition matrix is identified using the history of the estimated resistivity distribution obtained by a sensitivity matrix based algorithm and a Newton-Raphson algorithm. To numerically identify the linear evolution model, the Ibrahim time-domain method is used. The investigation is performed by numerical simulations of a domain with time-varying resistivity and by experimental data collected from the boundary of a human chest during normal breathing. The obtained dynamic resistivity values lie within the expected values for the tissues of a human chest. The EKF results suggest that the tracking ability is significantly improved with this approach.

CT to Cone-beam CT Deformable Registration With Simultaneous Intensity Correction

PubMed Central

Zhen, Xin; Gu, Xuejun; Yan, Hao; Zhou, Linghong; Jia, Xun; Jiang, Steve B.

2012-01-01

Computed tomography (CT) to cone-beam computed tomography (CBCT) deformable image registration (DIR) is a crucial step in adaptive radiation therapy. Current intensity-based registration algorithms, such as demons, may fail in the context of CT-CBCT DIR because of inconsistent intensities between the two modalities. In this paper, we propose a variant of demons, called Deformation with Intensity Simultaneously Corrected (DISC), to deal with CT-CBCT DIR. DISC distinguishes itself from the original demons algorithm by performing an adaptive intensity correction step on the CBCT image at every iteration step of the demons registration. Specifically, the intensity correction of a voxel in CBCT is achieved by matching the first and the second moments of the voxel intensities inside a patch around the voxel with those on the CT image. It is expected that such a strategy can remove artifacts in the CBCT image, as well as ensuring the intensity consistency between the two modalities. DISC is implemented on computer graphics processing units (GPUs) in compute unified device architecture (CUDA) programming environment. The performance of DISC is evaluated on a simulated patient case and six clinical head-and-neck cancer patient data. It is found that DISC is robust against the CBCT artifacts and intensity inconsistency and significantly improves the registration accuracy when compared with the original demons. PMID:23032638

Pattern Recognition Control Design

NASA Technical Reports Server (NTRS)

Gambone, Elisabeth

2016-01-01

Spacecraft control algorithms must know the expected spacecraft response to any command to the available control effectors, such as reaction thrusters or torque devices. Spacecraft control system design approaches have traditionally relied on the estimated vehicle mass properties to determine the desired force and moment, as well as knowledge of the effector performance to efficiently control the spacecraft. A pattern recognition approach can be used to investigate the relationship between the control effector commands and the spacecraft responses. Instead of supplying the approximated vehicle properties and the effector performance characteristics, a database of information relating the effector commands and the desired vehicle response can be used for closed-loop control. A Monte Carlo simulation data set of the spacecraft dynamic response to effector commands can be analyzed to establish the influence a command has on the behavior of the spacecraft. A tool developed at NASA Johnson Space Center (Ref. 1) to analyze flight dynamics Monte Carlo data sets through pattern recognition methods can be used to perform this analysis. Once a comprehensive data set relating spacecraft responses with commands is established, it can be used in place of traditional control laws and gains set. This pattern recognition approach can be compared with traditional control algorithms to determine the potential benefits and uses.

Cognitive programs: software for attention's executive

PubMed Central

Tsotsos, John K.; Kruijne, Wouter

2014-01-01

What are the computational tasks that an executive controller for visual attention must solve? This question is posed in the context of the Selective Tuning model of attention. The range of required computations go beyond top-down bias signals or region-of-interest determinations, and must deal with overt and covert fixations, process timing and synchronization, information routing, memory, matching control to task, spatial localization, priming, and coordination of bottom-up with top-down information. During task execution, results must be monitored to ensure the expected results. This description includes the kinds of elements that are common in the control of any kind of complex machine or system. We seek a mechanistic integration of the above, in other words, algorithms that accomplish control. Such algorithms operate on representations, transforming a representation of one kind into another, which then forms the input to yet another algorithm. Cognitive Programs (CPs) are hypothesized to capture exactly such representational transformations via stepwise sequences of operations. CPs, an updated and modernized offspring of Ullman's Visual Routines, impose an algorithmic structure to the set of attentional functions and play a role in the overall shaping of attentional modulation of the visual system so that it provides its best performance. This requires that we consider the visual system as a dynamic, yet general-purpose processor tuned to the task and input of the moment. This differs dramatically from the almost universal cognitive and computational views, which regard vision as a passively observing module to which simple questions about percepts can be posed, regardless of task. Differing from Visual Routines, CPs explicitly involve the critical elements of Visual Task Executive (vTE), Visual Attention Executive (vAE), and Visual Working Memory (vWM). Cognitive Programs provide the software that directs the actions of the Selective Tuning model of visual attention. PMID:25505430

Flight Deck Surface Trajectory-based Operations (STBO): Results of Piloted Simulations and Implications for Concepts of Operation (ConOps)

NASA Technical Reports Server (NTRS)

Foyle, David C.; Hooey, Becky L.; Bakowski, Deborah L.

2013-01-01

The results offour piloted medium-fidelity simulations investigating flight deck surface trajectory-based operations (STBO) will be reviewed. In these flight deck STBO simulations, commercial transport pilots were given taxi clearances with time and/or speed components and required to taxi to the departing runway or an intermediate traffic intersection. Under a variety of concept of operations (ConOps) and flight deck information conditions, pilots' ability to taxi in compliance with the required time of arrival (RTA) at the designated airport location was measured. ConOps and flight deck information conditions explored included: Availability of taxi clearance speed and elapsed time information; Intermediate RTAs at intermediate time constraint points (e.g., intersection traffic flow points); STBO taxi clearances via ATC voice speed commands or datal ink; and, Availability of flight deck display algorithms to reduce STBO RTA error. Flight Deck Implications. Pilot RTA conformance for STBO clearances, in the form of ATC taxi clearances with associated speed requirements, was found to be relatively poor, unless the pilot is required to follow a precise speed and acceleration/deceleration profile. However, following such a precise speed profile results in inordinate head-down tracking of current ground speed, leading to potentially unsafe operations. Mitigating these results, and providing good taxi RTA performance without the associated safety issues, is a flight deck avionics or electronic flight bag (EFB) solution. Such a solution enables pilots to meet the taxi route RTA without moment-by-moment tracking of ground speed. An avionics or EFB "error-nulling" algorithm allows the pilot to view the STBO information when the pilot determines it is necessary and when workload alloys, thus enabling the pilot to spread his/her attention appropriately and strategically on aircraft separation airport navigation, and the many other flight deck tasks concurrently required. Surface Traffic Management (STM) System Implications. The data indicate a number of implications regarding specific parameters for ATC/STM algorithm development. Pilots have a tendency to arrive at RTA points early with slow required speeds, on time for moderate speeds, and late with faster required speeds. This implies that ATC/STM algorithms should operate with middle-range speeds, similar to that of non-STBO taxi performance. Route length has a related effect: Long taxi routes increase the earliness with slow speeds and the lateness with faster speeds. This is likely due to the" open-loop" nature of the task in which the speed error compounds over a longer time with longer routes. Results showed that this may be mitigated by imposing a small number oftime constraint points each with their own RTAs effectively tuming a long route into a series of shorter routes - and thus improving RTA performance. STBO ConOps Implications. Most important is the impact that these data have for NextGen STM system ConOps development. The results of these experiments imply that it is not reasonable to expect pilots to taxi under a "Full STBO" ConOps in which pilots are expected to be at a predictable (x,y) airport location for every time (t). An STBO ConOps with a small number of intermediate time constraint points and the departing runway, however, is feasible, but only with flight deck equipage enabling the use of a display similar to the "error-nulling algorithm/display" tested.

Algorithms for the Reduction of Wind-Tunnel Data Derived from Strain Gauge Force Balances.

DTIC Science & Technology

1984-05-01

summed. Where hinge moments are measured on a model, it is customary to express them by coefficients of the form C11 h (4.23) q Si dH where hi is the...measured hinge moment and Sit and dH are a characteristic area and length associated with the control surface. 4.6 Transformation to Body Axes...Pty. Ltd. Mr D. Pilkington Mr R. D. Bullen Commonwealth Aircraft Corporation, Libra Hawker de Havilland Aust. Pty. Ltd., Bankstown. L.ibrar

Handling of computational in vitro/in vivo correlation problems by Microsoft Excel II. Distribution functions and moments.

PubMed

Langenbucher, Frieder

2003-01-01

MS Excel is a useful tool to handle in vitro/in vivo correlation (IVIVC) distribution functions, with emphasis on the Weibull and the biexponential distribution, which are most useful for the presentation of cumulative profiles, e.g. release in vitro or urinary excretion in vivo, and differential profiles such as the plasma response in vivo. The discussion includes moments (AUC and mean) as summarizing statistics, and data-fitting algorithms for parameter estimation.

Effective equations for the quantum pendulum from momentous quantum mechanics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hernandez, Hector H.; Chacon-Acosta, Guillermo; Departamento de Matematicas Aplicadas y Sistemas, Universidad Autonoma Metropolitana-Cuajimalpa, Artificios 40, Mexico D. F. 01120

In this work we study the quantum pendulum within the framework of momentous quantum mechanics. This description replaces the Schroedinger equation for the quantum evolution of the system with an infinite set of classical equations for expectation values of configuration variables, and quantum dispersions. We solve numerically the effective equations up to the second order, and describe its evolution.

Geomagnetic dipole strength and reversal rate over the past two million years.

PubMed

Valet, Jean-Pierre; Meynadier, Laure; Guyodo, Yohan

2005-06-09

Independent records of relative magnetic palaeointensity from sediment cores in different areas of the world can be stacked together to extract the evolution of the geomagnetic dipole moment and thus provide information regarding the processes governing the geodynamo. So far, this procedure has been limited to the past 800,000 years (800 kyr; ref. 3), which does not include any geomagnetic reversals. Here we present a composite curve that shows the evolution of the dipole moment during the past two million years. This reconstruction is in good agreement with the absolute dipole moments derived from volcanic lavas, which were used for calibration. We show that, at least during this period, the time-averaged field was higher during periods without reversals but the amplitude of the short-term oscillations remained the same. As a consequence, few intervals of very low intensity, and thus fewer instabilities, are expected during periods with a strong average dipole moment, whereas more excursions and reversals are expected during periods of weak field intensity. We also observe that the axial dipole begins to decay 60-80 kyr before reversals, but rebuilds itself in the opposite direction in only a few thousand years.

Exact Dissipative Moment Closures for Simulation of Magnetospheric Plasmas

NASA Astrophysics Data System (ADS)

Newman, D. L.; Sen, N.; Goldman, M. V.

2004-11-01

Dissipative fluid closures produce a kinetic-like plasma response in simulations based on the evolution of moments of the Vlasov equation. Such methods were previously shown to approximate the kinetic susceptibility of a Maxwellian plasma.(G. W. Hammett and F. W. Perkins Phys. Rev. Lett.) 64, 3019 (1990). We show here that dissipative closures can yield the exact linear response for kappa velocity distributions (i.e., f(v)∝(v^2+w^2)^-κ in 1-D, where w∝ v_th), provided κ is an integer and κ+1 moments are retained in the closure. This finding is particularly relevant to the simulation of collisionless space plasmas, which frequently exhibit power-law tails characteristic of kappa distributions. Such dissipative algorithms can be made energy conserving by evolving the thermal parameter w. Dominant nonlinearities (e.g., ponderomotive effects) can also be incorporated into the algorithm. These methods have proven especially valuable in the context of reduced 2-D Vlasov simulations,(N. Sen, et al., Reduced 2-D Vlasov Simulationsldots), this meeting. where they have been used to model perpendicular ion dynamics in the evolution of nonlinear structures (e.g., double layers) in the auroral ionosphere.

Leonid predictions for the period 2001-2100

NASA Astrophysics Data System (ADS)

Maslov, Mikhail

2007-02-01

This article provides a set of summaries of what to expect from the Leonid meteor shower for each year of the period 2001-2100. Each summary contains the moments of maximum/maxima, their expected intensity and some comments about average meteor brightness during them. Special attention was paid to background (traditional) maxima, which are characterized with their expected times and intensities.

Recurrent procedure for constructing nonisotropic matrix elements of the collision integral of the nonlinear Boltzmann equation

NASA Astrophysics Data System (ADS)

Ender, I. A.; Bakaleinikov, L. A.; Flegontova, E. Yu.; Gerasimenko, A. B.

2017-08-01

We have proposed an algorithm for the sequential construction of nonisotropic matrix elements of the collision integral, which are required to solve the nonlinear Boltzmann equation using the moments method. The starting elements of the matrix are isotropic and assumed to be known. The algorithm can be used for an arbitrary law of interactions for any ratio of the masses of colliding particles.

Hierarchical trie packet classification algorithm based on expectation-maximization clustering

PubMed Central

Bi, Xia-an; Zhao, Junxia

2017-01-01

With the development of computer network bandwidth, packet classification algorithms which are able to deal with large-scale rule sets are in urgent need. Among the existing algorithms, researches on packet classification algorithms based on hierarchical trie have become an important packet classification research branch because of their widely practical use. Although hierarchical trie is beneficial to save large storage space, it has several shortcomings such as the existence of backtracking and empty nodes. This paper proposes a new packet classification algorithm, Hierarchical Trie Algorithm Based on Expectation-Maximization Clustering (HTEMC). Firstly, this paper uses the formalization method to deal with the packet classification problem by means of mapping the rules and data packets into a two-dimensional space. Secondly, this paper uses expectation-maximization algorithm to cluster the rules based on their aggregate characteristics, and thereby diversified clusters are formed. Thirdly, this paper proposes a hierarchical trie based on the results of expectation-maximization clustering. Finally, this paper respectively conducts simulation experiments and real-environment experiments to compare the performances of our algorithm with other typical algorithms, and analyzes the results of the experiments. The hierarchical trie structure in our algorithm not only adopts trie path compression to eliminate backtracking, but also solves the problem of low efficiency of trie updates, which greatly improves the performance of the algorithm. PMID:28704476

A solution algorithm for fluid–particle flows across all flow regimes

DOE PAGES

Kong, Bo; Fox, Rodney O.

2017-05-12

Many fluid–particle flows occurring in nature and in technological applications exhibit large variations in the local particle volume fraction. For example, in circulating fluidized beds there are regions where the particles are closepacked as well as very dilute regions where particle–particle collisions are rare. Thus, in order to simulate such fluid–particle systems, it is necessary to design a flow solver that can accurately treat all flow regimes occurring simultaneously in the same flow domain. In this work, a solution algorithm is proposed for this purpose. The algorithm is based on splitting the free-transport flux solver dynamically and locally in themore » flow. In close-packed to moderately dense regions, a hydrodynamic solver is employed, while in dilute to very dilute regions a kinetic-based finite-volume solver is used in conjunction with quadrature-based moment methods. To illustrate the accuracy and robustness of the proposed solution algorithm, it is implemented in OpenFOAM for particle velocity moments up to second order, and applied to simulate gravity-driven, gas–particle flows exhibiting cluster-induced turbulence. By varying the average particle volume fraction in the flow domain, it is demonstrated that the flow solver can handle seamlessly all flow regimes present in fluid–particle flows.« less

A solution algorithm for fluid-particle flows across all flow regimes

NASA Astrophysics Data System (ADS)

Kong, Bo; Fox, Rodney O.

2017-09-01

Many fluid-particle flows occurring in nature and in technological applications exhibit large variations in the local particle volume fraction. For example, in circulating fluidized beds there are regions where the particles are close-packed as well as very dilute regions where particle-particle collisions are rare. Thus, in order to simulate such fluid-particle systems, it is necessary to design a flow solver that can accurately treat all flow regimes occurring simultaneously in the same flow domain. In this work, a solution algorithm is proposed for this purpose. The algorithm is based on splitting the free-transport flux solver dynamically and locally in the flow. In close-packed to moderately dense regions, a hydrodynamic solver is employed, while in dilute to very dilute regions a kinetic-based finite-volume solver is used in conjunction with quadrature-based moment methods. To illustrate the accuracy and robustness of the proposed solution algorithm, it is implemented in OpenFOAM for particle velocity moments up to second order, and applied to simulate gravity-driven, gas-particle flows exhibiting cluster-induced turbulence. By varying the average particle volume fraction in the flow domain, it is demonstrated that the flow solver can handle seamlessly all flow regimes present in fluid-particle flows.

A solution algorithm for fluid–particle flows across all flow regimes

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kong, Bo; Fox, Rodney O.

Many fluid–particle flows occurring in nature and in technological applications exhibit large variations in the local particle volume fraction. For example, in circulating fluidized beds there are regions where the particles are closepacked as well as very dilute regions where particle–particle collisions are rare. Thus, in order to simulate such fluid–particle systems, it is necessary to design a flow solver that can accurately treat all flow regimes occurring simultaneously in the same flow domain. In this work, a solution algorithm is proposed for this purpose. The algorithm is based on splitting the free-transport flux solver dynamically and locally in themore » flow. In close-packed to moderately dense regions, a hydrodynamic solver is employed, while in dilute to very dilute regions a kinetic-based finite-volume solver is used in conjunction with quadrature-based moment methods. To illustrate the accuracy and robustness of the proposed solution algorithm, it is implemented in OpenFOAM for particle velocity moments up to second order, and applied to simulate gravity-driven, gas–particle flows exhibiting cluster-induced turbulence. By varying the average particle volume fraction in the flow domain, it is demonstrated that the flow solver can handle seamlessly all flow regimes present in fluid–particle flows.« less

Filament capturing with the multimaterial moment-of-fluid method*

DOE PAGES

Jemison, Matthew; Sussman, Mark; Shashkov, Mikhail

2015-01-15

A novel method for capturing two-dimensional, thin, under-resolved material configurations, known as “filaments,” is presented in the context of interface reconstruction. This technique uses a partitioning procedure to detect disconnected regions of material in the advective preimage of a cell (indicative of a filament) and makes use of the existing functionality of the Multimaterial Moment-of-Fluid interface reconstruction method to accurately capture the under-resolved feature, while exactly conserving volume. An algorithm for Adaptive Mesh Refinement in the presence of filaments is developed so that refinement is introduced only near the tips of filaments and where the Moment-of-Fluid reconstruction error is stillmore » large. Comparison to the standard Moment-of-Fluid method is made. As a result, it is demonstrated that using filament capturing at a given resolution yields gains in accuracy comparable to introducing an additional level of mesh refinement at significantly lower cost.« less

Pattern Discovery and Change Detection of Online Music Query Streams

NASA Astrophysics Data System (ADS)

Li, Hua-Fu

In this paper, an efficient stream mining algorithm, called FTP-stream (Frequent Temporal Pattern mining of streams), is proposed to find the frequent temporal patterns over melody sequence streams. In the framework of our proposed algorithm, an effective bit-sequence representation is used to reduce the time and memory needed to slide the windows. The FTP-stream algorithm can calculate the support threshold in only a single pass based on the concept of bit-sequence representation. It takes the advantage of "left" and "and" operations of the representation. Experiments show that the proposed algorithm only scans the music query stream once, and runs significant faster and consumes less memory than existing algorithms, such as SWFI-stream and Moment.

Lorentz-covariant coordinate-space representation of the leading hadronic contribution to the anomalous magnetic moment of the muon

NASA Astrophysics Data System (ADS)

Meyer, Harvey B.

2017-09-01

We present a Lorentz-covariant, Euclidean coordinate-space expression for the hadronic vacuum polarisation, the Adler function and the leading hadronic contribution to the anomalous magnetic moment of the muon. The representation offers a high degree of flexibility for an implementation in lattice QCD. We expect it to be particularly helpful for the quark-line disconnected contributions.

Theoretical electric dipole moments of SiH, GeH and SnH

NASA Technical Reports Server (NTRS)

Pettersson, L. G. M.; Langhoff, S. R.

1986-01-01

Accurate theoretical dipole moments have been computed for the X2Pi ground states of Si(-)H(+) (0.118 D), Ge(+)H(-) (0.085 D), and Sn(+)H(-) (0.357 D). The trend down the periodic table is regular and follows that expected from the electronegativities of the group IV atoms. The dipole moment of 1.24 + or - 0.1 D for GeH recently derived by Brown, Evenson and Sears (1985) from the relative intensities of electric and magnetic dipole transitions in the 10-micron spectrum of the X2Pi state is seriously questioned.

Theoretical Electric Dipole Moments of SiH, GeH and SnH

NASA Technical Reports Server (NTRS)

Pettersson, Lars G. M.; Langhoff, Stephen R.

1986-01-01

Accurate theoretical dipole moments (mu(sub c) have been computed for the X(exp 2)Pi ground states of Si(-)H(+)(0.118 D), Ge(+)H(-)(0.085 D) and Sn(+)H(-)(0.357 D). The trend down the periodic table is regular and follows that expected from the electronegativities of the group IV atoms. The dipole moment of 1.24 +/- 0.1 D for GeH recently derived by Brown, Evenson and Sears from the relative intensities of electric and magnetic dipole transitions in the 10 microns spectrum of the X(exp 2)Pi state is seriously questioned.

Formulation of the relativistic moment implicit particle-in-cell method

DOE Office of Scientific and Technical Information (OSTI.GOV)

Noguchi, Koichi; Tronci, Cesare; Zuccaro, Gianluca

2007-04-15

A new formulation is presented for the implicit moment method applied to the time-dependent relativistic Vlasov-Maxwell system. The new approach is based on a specific formulation of the implicit moment method that allows us to retain the same formalism that is valid in the classical case despite the formidable complication introduced by the nonlinear nature of the relativistic equations of motion. To demonstrate the validity of the new formulation, an implicit finite difference algorithm is developed to solve the Maxwell's equations and equations of motion. A number of benchmark problems are run: two stream instability, ion acoustic wave damping, Weibelmore » instability, and Poynting flux acceleration. The numerical results are all in agreement with analytical solutions.« less

Recent National Transonic Facility Test Process Improvements (Invited)

NASA Technical Reports Server (NTRS)

Kilgore, W. A.; Balakrishna, S.; Bobbitt, C. W., Jr.; Adcock, J. B.

2001-01-01

This paper describes the results of two recent process improvements; drag feed-forward Mach number control and simultaneous force/moment and pressure testing, at the National Transonic Facility. These improvements have reduced the duration and cost of testing. The drag feed-forward Mach number control reduces the Mach number settling time by using measured model drag in the Mach number control algorithm. Simultaneous force/moment and pressure testing allows simultaneous collection of force/moment and pressure data without sacrificing data quality thereby reducing the overall testing time. Both improvements can be implemented at any wind tunnel. Additionally the NTF is working to develop and implement continuous pitch as a testing option as an additional method to reduce costs and maintain data quality.

Recent National Transonic Facility Test Process Improvements (Invited)

NASA Technical Reports Server (NTRS)

Kilgore, W. A.; Balakrishna, S.; Bobbitt, C. W., Jr.; Adcock, J. B.

2001-01-01

This paper describes the results of two recent process improvements; drag feed-forward Mach number control and simultaneous force/moment and pressure testing, at the National Transonic Facility. These improvements have reduced the duration and cost of testing. The drag feedforward Mach number control reduces the Mach number settling time by using measured model drag in the Mach number control algorithm. Simultaneous force/moment and pressure testing allows simultaneous collection of force/moment and pressure data without sacrificing data quality thereby reducing the overall testing time. Both improvements can be implemented at any wind tunnel. Additionally the NTF is working to develop and implement continuous pitch as a testing option as an additional method to reduce costs and maintain data quality.

Speed-up of the volumetric method of moments for the approximate RCS of large arbitrary-shaped dielectric targets

NASA Astrophysics Data System (ADS)

Moreno, Javier; Somolinos, Álvaro; Romero, Gustavo; González, Iván; Cátedra, Felipe

2017-08-01

A method for the rigorous computation of the electromagnetic scattering of large dielectric volumes is presented. One goal is to simplify the analysis of large dielectric targets with translational symmetries taken advantage of their Toeplitz symmetry. Then, the matrix-fill stage of the Method of Moments is efficiently obtained because the number of coupling terms to compute is reduced. The Multilevel Fast Multipole Method is applied to solve the problem. Structured meshes are obtained efficiently to approximate the dielectric volumes. The regular mesh grid is achieved by using parallelepipeds whose centres have been identified as internal to the target. The ray casting algorithm is used to classify the parallelepiped centres. It may become a bottleneck when too many points are evaluated in volumes defined by parametric surfaces, so a hierarchical algorithm is proposed to minimize the number of evaluations. Measurements and analytical results are included for validation purposes.

Intra-pulse modulation recognition using short-time ramanujan Fourier transform spectrogram

NASA Astrophysics Data System (ADS)

Ma, Xiurong; Liu, Dan; Shan, Yunlong

2017-12-01

Intra-pulse modulation recognition under negative signal-to-noise ratio (SNR) environment is a research challenge. This article presents a robust algorithm for the recognition of 5 types of radar signals with large variation range in the signal parameters in low SNR using the combination of the Short-time Ramanujan Fourier transform (ST-RFT) and pseudo-Zernike moments invariant features. The ST-RFT provides the time-frequency distribution features for 5 modulations. The pseudo-Zernike moments provide invariance properties that are able to recognize different modulation schemes on different parameter variation conditions from the ST-RFT spectrograms. Simulation results demonstrate that the proposed algorithm achieves the probability of successful recognition (PSR) of over 90% when SNR is above -5 dB with large variation range in the signal parameters: carrier frequency (CF) for all considered signals, hop size (HS) for frequency shift keying (FSK) signals, and the time-bandwidth product for Linear Frequency Modulation (LFM) signals.

FPGA-Based Multimodal Embedded Sensor System Integrating Low- and Mid-Level Vision

PubMed Central

Botella, Guillermo; Martín H., José Antonio; Santos, Matilde; Meyer-Baese, Uwe

2011-01-01

Motion estimation is a low-level vision task that is especially relevant due to its wide range of applications in the real world. Many of the best motion estimation algorithms include some of the features that are found in mammalians, which would demand huge computational resources and therefore are not usually available in real-time. In this paper we present a novel bioinspired sensor based on the synergy between optical flow and orthogonal variant moments. The bioinspired sensor has been designed for Very Large Scale Integration (VLSI) using properties of the mammalian cortical motion pathway. This sensor combines low-level primitives (optical flow and image moments) in order to produce a mid-level vision abstraction layer. The results are described trough experiments showing the validity of the proposed system and an analysis of the computational resources and performance of the applied algorithms. PMID:22164069

A Framework for Optimal Control Allocation with Structural Load Constraints

NASA Technical Reports Server (NTRS)

Frost, Susan A.; Taylor, Brian R.; Jutte, Christine V.; Burken, John J.; Trinh, Khanh V.; Bodson, Marc

2010-01-01

Conventional aircraft generally employ mixing algorithms or lookup tables to determine control surface deflections needed to achieve moments commanded by the flight control system. Control allocation is the problem of converting desired moments into control effector commands. Next generation aircraft may have many multipurpose, redundant control surfaces, adding considerable complexity to the control allocation problem. These issues can be addressed with optimal control allocation. Most optimal control allocation algorithms have control surface position and rate constraints. However, these constraints are insufficient to ensure that the aircraft's structural load limits will not be exceeded by commanded surface deflections. In this paper, a framework is proposed to enable a flight control system with optimal control allocation to incorporate real-time structural load feedback and structural load constraints. A proof of concept simulation that demonstrates the framework in a simulation of a generic transport aircraft is presented.

Bias correction of daily satellite precipitation data using genetic algorithm

NASA Astrophysics Data System (ADS)

Pratama, A. W.; Buono, A.; Hidayat, R.; Harsa, H.

2018-05-01

Climate Hazards Group InfraRed Precipitation with Stations (CHIRPS) was producted by blending Satellite-only Climate Hazards Group InfraRed Precipitation (CHIRP) with Stasion observations data. The blending process was aimed to reduce bias of CHIRP. However, Biases of CHIRPS on statistical moment and quantil values were high during wet season over Java Island. This paper presented a bias correction scheme to adjust statistical moment of CHIRP using observation precipitation data. The scheme combined Genetic Algorithm and Nonlinear Power Transformation, the results was evaluated based on different season and different elevation level. The experiment results revealed that the scheme robustly reduced bias on variance around 100% reduction and leaded to reduction of first, and second quantile biases. However, bias on third quantile only reduced during dry months. Based on different level of elevation, the performance of bias correction process is only significantly different on skewness indicators.

FPGA-based multimodal embedded sensor system integrating low- and mid-level vision.

PubMed

Botella, Guillermo; Martín H, José Antonio; Santos, Matilde; Meyer-Baese, Uwe

2011-01-01

Motion estimation is a low-level vision task that is especially relevant due to its wide range of applications in the real world. Many of the best motion estimation algorithms include some of the features that are found in mammalians, which would demand huge computational resources and therefore are not usually available in real-time. In this paper we present a novel bioinspired sensor based on the synergy between optical flow and orthogonal variant moments. The bioinspired sensor has been designed for Very Large Scale Integration (VLSI) using properties of the mammalian cortical motion pathway. This sensor combines low-level primitives (optical flow and image moments) in order to produce a mid-level vision abstraction layer. The results are described trough experiments showing the validity of the proposed system and an analysis of the computational resources and performance of the applied algorithms.

Anomalous thermal hysteresis in the high-field magnetic moments of magnetic nanoparticles embedded in multi-walled carbon nanotubes

NASA Astrophysics Data System (ADS)

Zhao, Guo-Meng; Wang, Jun; Ren, Yang; Beeli, Pieder

2012-02-01

We report high-temperature (300-1120 K) magnetic properties of Fe and Fe3O4 nanoparticles embedded in multi-walled carbon nanotubes. We unambiguously show that the magnetic moments of Fe and Fe3O4 nanoparticles are seemingly enhanced by a factor of about 3 compared with what they would be expected to have for free (unembedded) magnetic nanoparticles. What is more intriguing is that the enhanced moments were completely lost when the sample was heated up to 1120 K and the lost moments at 1120 K were completely recovered through several thermal cycles below 1020 K. The anomalous thermal hysteresis of the high-field magnetic moments is unlikely to be explained by existing physical models except for the high-field paramagnetic Meissner effect due to the existence of ultrahigh temperature superconductivity in the multi-walled carbon nanotubes.

Dike propagation energy balance from deformation modeling and seismic release

NASA Astrophysics Data System (ADS)

Bonaccorso, Alessandro; Aoki, Yosuke; Rivalta, Eleonora

2017-06-01

Magma is transported in the crust mainly by dike intrusions. In volcanic areas, dikes can ascend toward the free surface and also move by lateral propagation, eventually feeding flank eruptions. Understanding dike mechanics is a key to forecasting the expected propagation and associated hazard. Several studies have been conducted on dike mechanisms and propagation; however, a less in-depth investigated aspect is the relation between measured dike-induced deformation and the seismicity released during its propagation. We individuated a simple x that can be used as a proxy of the expected mechanical energy released by a propagating dike and is related to its average thickness. For several intrusions around the world (Afar, Japan, and Mount Etna), we correlate such mechanical energy to the seismic moment released by the induced earthquakes. We obtain an empirical law that quantifies the expected seismic energy released before arrest. The proposed approach may be helpful to predict the total seismic moment that will be released by an intrusion and thus to control the energy status during its propagation and the time of dike arrest.