Inverse dynamics of adaptive structures used as space cranes
NASA Technical Reports Server (NTRS)
Das, S. K.; Utku, S.; Wada, B. K.
1990-01-01
As a precursor to the real-time control of fast moving adaptive structures used as space cranes, a formulation is given for the flexibility induced motion relative to the nominal motion (i.e., the motion that assumes no flexibility) and for obtaining the open loop time varying driving forces. An algorithm is proposed for the computation of the relative motion and driving forces. The governing equations are given in matrix form with explicit functional dependencies. A simulator is developed to implement the algorithm on a digital computer. In the formulations, the distributed mass of the crane is lumped by two schemes, vz., 'trapezoidal' lumping and 'Simpson's rule' lumping. The effects of the mass lumping schemes are shown by simulator runs.
NASA Technical Reports Server (NTRS)
Campbell, Stefan F.; Kaneshige, John T.
2010-01-01
Presented here is a Predictor-Based Model Reference Adaptive Control (PMRAC) architecture for a generic transport aircraft. At its core, this architecture features a three-axis, non-linear, dynamic-inversion controller. Command inputs for this baseline controller are provided by pilot roll-rate, pitch-rate, and sideslip commands. This paper will first thoroughly present the baseline controller followed by a description of the PMRAC adaptive augmentation to this control system. Results are presented via a full-scale, nonlinear simulation of NASA s Generic Transport Model (GTM).
Adaptive dynamic inversion robust control for BTT missile based on wavelet neural network
NASA Astrophysics Data System (ADS)
Li, Chuanfeng; Wang, Yongji; Deng, Zhixiang; Wu, Hao
2009-10-01
A new nonlinear control strategy incorporated the dynamic inversion method with wavelet neural networks is presented for the nonlinear coupling system of Bank-to-Turn(BTT) missile in reentry phase. The basic control law is designed by using the dynamic inversion feedback linearization method, and the online learning wavelet neural network is used to compensate the inversion error due to aerodynamic parameter errors, modeling imprecise and external disturbance in view of the time-frequency localization properties of wavelet transform. Weights adjusting laws are derived according to Lyapunov stability theory, which can guarantee the boundedness of all signals in the whole system. Furthermore, robust stability of the closed-loop system under this tracking law is proved. Finally, the six degree-of-freedom(6DOF) simulation results have shown that the attitude angles can track the anticipant command precisely under the circumstances of existing external disturbance and in the presence of parameter uncertainty. It means that the dependence on model by dynamic inversion method is reduced and the robustness of control system is enhanced by using wavelet neural network(WNN) to reconstruct inversion error on-line.
Hybrid Adaptive Flight Control with Model Inversion Adaptation
NASA Technical Reports Server (NTRS)
Nguyen, Nhan
2011-01-01
This study investigates a hybrid adaptive flight control method as a design possibility for a flight control system that can enable an effective adaptation strategy to deal with off-nominal flight conditions. The hybrid adaptive control blends both direct and indirect adaptive control in a model inversion flight control architecture. The blending of both direct and indirect adaptive control provides a much more flexible and effective adaptive flight control architecture than that with either direct or indirect adaptive control alone. The indirect adaptive control is used to update the model inversion controller by an on-line parameter estimation of uncertain plant dynamics based on two methods. The first parameter estimation method is an indirect adaptive law based on the Lyapunov theory, and the second method is a recursive least-squares indirect adaptive law. The model inversion controller is therefore made to adapt to changes in the plant dynamics due to uncertainty. As a result, the modeling error is reduced that directly leads to a decrease in the tracking error. In conjunction with the indirect adaptive control that updates the model inversion controller, a direct adaptive control is implemented as an augmented command to further reduce any residual tracking error that is not entirely eliminated by the indirect adaptive control.
Zhang, Zhihua; Sheng, Zheng; Shi, Hanqing; Fan, Zhiqiang
2016-01-01
Using the RFC technique to estimate refractivity parameters is a complex nonlinear optimization problem. In this paper, an improved cuckoo search (CS) algorithm is proposed to deal with this problem. To enhance the performance of the CS algorithm, a parameter dynamic adaptive operation and crossover operation were integrated into the standard CS (DACS-CO). Rechenberg's 1/5 criteria combined with learning factor were used to control the parameter dynamic adaptive adjusting process. The crossover operation of genetic algorithm was utilized to guarantee the population diversity. The new hybrid algorithm has better local search ability and contributes to superior performance. To verify the ability of the DACS-CO algorithm to estimate atmospheric refractivity parameters, the simulation data and real radar clutter data are both implemented. The numerical experiments demonstrate that the DACS-CO algorithm can provide an effective method for near-real-time estimation of the atmospheric refractivity profile from radar clutter. PMID:27212938
Zhang, Zhihua; Sheng, Zheng; Shi, Hanqing; Fan, Zhiqiang
2016-01-01
Using the RFC technique to estimate refractivity parameters is a complex nonlinear optimization problem. In this paper, an improved cuckoo search (CS) algorithm is proposed to deal with this problem. To enhance the performance of the CS algorithm, a parameter dynamic adaptive operation and crossover operation were integrated into the standard CS (DACS-CO). Rechenberg's 1/5 criteria combined with learning factor were used to control the parameter dynamic adaptive adjusting process. The crossover operation of genetic algorithm was utilized to guarantee the population diversity. The new hybrid algorithm has better local search ability and contributes to superior performance. To verify the ability of the DACS-CO algorithm to estimate atmospheric refractivity parameters, the simulation data and real radar clutter data are both implemented. The numerical experiments demonstrate that the DACS-CO algorithm can provide an effective method for near-real-time estimation of the atmospheric refractivity profile from radar clutter. PMID:27212938
Adaptive Inverse Hyperbolic Tangent Algorithm for Dynamic Contrast Adjustment in Displaying Scenes
NASA Astrophysics Data System (ADS)
Yu, Cheng-Yi; Ouyang, Yen-Chieh; Wang, Chuin-Mu; Chang, Chein-I.
2010-12-01
Contrast has a great influence on the quality of an image in human visual perception. A poorly illuminated environment can significantly affect the contrast ratio, producing an unexpected image. This paper proposes an Adaptive Inverse Hyperbolic Tangent (AIHT) algorithm to improve the display quality and contrast of a scene. Because digital cameras must maintain the shadow in a middle range of luminance that includes a main object such as a face, a gamma function is generally used for this purpose. However, this function has a severe weakness in that it decreases highlight contrast. To mitigate this problem, contrast enhancement algorithms have been designed to adjust contrast to tune human visual perception. The proposed AIHT determines the contrast levels of an original image as well as parameter space for different contrast types so that not only the original histogram shape features can be preserved, but also the contrast can be enhanced effectively. Experimental results show that the proposed algorithm is capable of enhancing the global contrast of the original image adaptively while extruding the details of objects simultaneously.
Adaptation through chromosomal inversions in Anopheles.
Ayala, Diego; Ullastres, Anna; González, Josefa
2014-01-01
Chromosomal inversions have been repeatedly involved in local adaptation in a large number of animals and plants. The ecological and behavioral plasticity of Anopheles species-human malaria vectors-is mirrored by high amounts of polymorphic inversions. The adaptive significance of chromosomal inversions has been consistently attested by strong and significant correlations between their frequencies and a number of phenotypic traits. Here, we provide an extensive literature review of the different adaptive traits associated with chromosomal inversions in the genus Anopheles. Traits having important consequences for the success of present and future vector control measures, such as insecticide resistance and behavioral changes, are discussed. PMID:24904633
Adaptive regularization of earthquake slip distribution inversion
NASA Astrophysics Data System (ADS)
Wang, Chisheng; Ding, Xiaoli; Li, Qingquan; Shan, Xinjian; Zhu, Jiasong; Guo, Bo; Liu, Peng
2016-04-01
Regularization is a routine approach used in earthquake slip distribution inversion to avoid numerically abnormal solutions. To date, most slip inversion studies have imposed uniform regularization on all the fault patches. However, adaptive regularization, where each retrieved parameter is regularized differently, has exhibited better performances in other research fields such as image restoration. In this paper, we implement an investigation into adaptive regularization for earthquake slip distribution inversion. It is found that adaptive regularization can achieve a significantly smaller mean square error (MSE) than uniform regularization, if it is set properly. We propose an adaptive regularization method based on weighted total least squares (WTLS). This approach assumes that errors exist in both the regularization matrix and observation, and an iterative algorithm is used to solve the solution. A weight coefficient is used to balance the regularization matrix residual and the observation residual. An experiment using four slip patterns was carried out to validate the proposed method. The results show that the proposed regularization method can derive a smaller MSE than uniform regularization and resolution-based adaptive regularization, and the improvement in MSE is more significant for slip patterns with low-resolution slip patches. In this paper, we apply the proposed regularization method to study the slip distribution of the 2011 Mw 9.0 Tohoku earthquake. The retrieved slip distribution is less smooth and more detailed than the one retrieved with the uniform regularization method, and is closer to the existing slip model from joint inversion of the geodetic and seismic data.
Adaptive Inverse optimal neuromuscular electrical stimulation.
Wang, Qiang; Sharma, Nitin; Johnson, Marcus; Gregory, Chris M; Dixon, Warren E
2013-12-01
Neuromuscular electrical stimulation (NMES) is a prescribed treatment for various neuromuscular disorders, where an electrical stimulus is provided to elicit a muscle contraction. Barriers to the development of NMES controllers exist because the muscle response to an electrical stimulation is nonlinear and the muscle model is uncertain. Efforts in this paper focus on the development of an adaptive inverse optimal NMES controller. The controller yields desired limb trajectory tracking while simultaneously minimizing a cost functional that is positive in the error states and stimulation input. The development of this framework allows tradeoffs to be made between tracking performance and control effort by putting different penalties on error states and control input, depending on the clinical goal or functional task. The controller is examined through a Lyapunov-based analysis. Experiments on able-bodied individuals are provided to demonstrate the performance of the developed controller. PMID:23757569
Chromosome inversions, adaptive cassettes and the evolution of species' ranges.
Kirkpatrick, Mark; Barrett, Brian
2015-05-01
A chromosome inversion can spread when it captures locally adapted alleles or when it is introduced into a species by hybridization with adapted alleles that were previously absent. We present a model that shows how both processes can cause a species range to expand. Introgression of an inversion that carries novel, locally adapted alleles is a particularly powerful mechanism for range expansion. The model supports the earlier proposal that introgression of an inversion triggered a large range expansion of a malaria mosquito. These results suggest a role for inversions as cassettes of genes that can accelerate adaptation by crossing species boundaries, rather than protecting genomes from introgression. PMID:25583098
Dynamically consistent Jacobian inverse for mobile manipulators
NASA Astrophysics Data System (ADS)
Ratajczak, Joanna; Tchoń, Krzysztof
2016-06-01
By analogy to the definition of the dynamically consistent Jacobian inverse for robotic manipulators, we have designed a dynamically consistent Jacobian inverse for mobile manipulators built of a non-holonomic mobile platform and a holonomic on-board manipulator. The endogenous configuration space approach has been exploited as a source of conceptual guidelines. The new inverse guarantees a decoupling of the motion in the operational space from the forces exerted in the endogenous configuration space and annihilated by the dual Jacobian inverse. A performance study of the new Jacobian inverse as a tool for motion planning is presented.
Dynamical Adaptation in Photoreceptors
Clark, Damon A.; Benichou, Raphael; Meister, Markus; Azeredo da Silveira, Rava
2013-01-01
Adaptation is at the heart of sensation and nowhere is it more salient than in early visual processing. Light adaptation in photoreceptors is doubly dynamical: it depends upon the temporal structure of the input and it affects the temporal structure of the response. We introduce a non-linear dynamical adaptation model of photoreceptors. It is simple enough that it can be solved exactly and simulated with ease; analytical and numerical approaches combined provide both intuition on the behavior of dynamical adaptation and quantitative results to be compared with data. Yet the model is rich enough to capture intricate phenomenology. First, we show that it reproduces the known phenomenology of light response and short-term adaptation. Second, we present new recordings and demonstrate that the model reproduces cone response with great precision. Third, we derive a number of predictions on the response of photoreceptors to sophisticated stimuli such as periodic inputs, various forms of flickering inputs, and natural inputs. In particular, we demonstrate that photoreceptors undergo rapid adaptation of response gain and time scale, over ∼ 300 ms—i. e., over the time scale of the response itself—and we confirm this prediction with data. For natural inputs, this fast adaptation can modulate the response gain more than tenfold and is hence physiologically relevant. PMID:24244119
Adaptive Dynamic Bayesian Networks
Ng, B M
2007-10-26
A discrete-time Markov process can be compactly modeled as a dynamic Bayesian network (DBN)--a graphical model with nodes representing random variables and directed edges indicating causality between variables. Each node has a probability distribution, conditional on the variables represented by the parent nodes. A DBN's graphical structure encodes fixed conditional dependencies between variables. But in real-world systems, conditional dependencies between variables may be unknown a priori or may vary over time. Model errors can result if the DBN fails to capture all possible interactions between variables. Thus, we explore the representational framework of adaptive DBNs, whose structure and parameters can change from one time step to the next: a distribution's parameters and its set of conditional variables are dynamic. This work builds on recent work in nonparametric Bayesian modeling, such as hierarchical Dirichlet processes, infinite-state hidden Markov networks and structured priors for Bayes net learning. In this paper, we will explain the motivation for our interest in adaptive DBNs, show how popular nonparametric methods are combined to formulate the foundations for adaptive DBNs, and present preliminary results.
Nonlinear Dynamic Inversion Baseline Control Law: Architecture and Performance Predictions
NASA Technical Reports Server (NTRS)
Miller, Christopher J.
2011-01-01
A model reference dynamic inversion control law has been developed to provide a baseline control law for research into adaptive elements and other advanced flight control law components. This controller has been implemented and tested in a hardware-in-the-loop simulation; the simulation results show excellent handling qualities throughout the limited flight envelope. A simple angular momentum formulation was chosen because it can be included in the stability proofs for many basic adaptive theories, such as model reference adaptive control. Many design choices and implementation details reflect the requirements placed on the system by the nonlinear flight environment and the desire to keep the system as basic as possible to simplify the addition of the adaptive elements. Those design choices are explained, along with their predicted impact on the handling qualities.
Fast Parallel Computation Of Manipulator Inverse Dynamics
NASA Technical Reports Server (NTRS)
Fijany, Amir; Bejczy, Antal K.
1991-01-01
Method for fast parallel computation of inverse dynamics problem, essential for real-time dynamic control and simulation of robot manipulators, undergoing development. Enables exploitation of high degree of parallelism and, achievement of significant computational efficiency, while minimizing various communication and synchronization overheads as well as complexity of required computer architecture. Universal real-time robotic controller and simulator (URRCS) consists of internal host processor and several SIMD processors with ring topology. Architecture modular and expandable: more SIMD processors added to match size of problem. Operate asynchronously and in MIMD fashion.
Genomic Evidence for Adaptive Inversion Clines in Drosophila melanogaster.
Kapun, Martin; Fabian, Daniel K; Goudet, Jérôme; Flatt, Thomas
2016-05-01
Clines in chromosomal inversion polymorphisms-presumably driven by climatic gradients-are common but there is surprisingly little evidence for selection acting on them. Here we address this long-standing issue in Drosophila melanogaster by using diagnostic single nucleotide polymorphism (SNP) markers to estimate inversion frequencies from 28 whole-genome Pool-seq samples collected from 10 populations along the North American east coast. Inversions In(3L)P, In(3R)Mo, and In(3R)Payne showed clear latitudinal clines, and for In(2L)t, In(2R)NS, and In(3R)Payne the steepness of the clinal slopes changed between summer and fall. Consistent with an effect of seasonality on inversion frequencies, we detected small but stable seasonal fluctuations of In(2R)NS and In(3R)Payne in a temperate Pennsylvanian population over 4 years. In support of spatially varying selection, we observed that the cline in In(3R)Payne has remained stable for >40 years and that the frequencies of In(2L)t and In(3R)Payne are strongly correlated with climatic factors that vary latitudinally, independent of population structure. To test whether these patterns are adaptive, we compared the amount of genetic differentiation of inversions versus neutral SNPs and found that the clines in In(2L)t and In(3R)Payne are maintained nonneutrally and independent of admixture. We also identified numerous clinal inversion-associated SNPs, many of which exhibit parallel differentiation along the Australian cline and reside in genes known to affect fitness-related traits. Together, our results provide strong evidence that inversion clines are maintained by spatially-and perhaps also temporally-varying selection. We interpret our data in light of current hypotheses about how inversions are established and maintained. PMID:26796550
An adaptive inverse kinematics algorithm for robot manipulators
NASA Technical Reports Server (NTRS)
Colbaugh, R.; Glass, K.; Seraji, H.
1990-01-01
An adaptive algorithm for solving the inverse kinematics problem for robot manipulators is presented. The algorithm is derived using model reference adaptive control (MRAC) theory and is computationally efficient for online applications. The scheme requires no a priori knowledge of the kinematics of the robot if Cartesian end-effector sensing is available, and it requires knowledge of only the forward kinematics if joint position sensing is used. Computer simulation results are given for the redundant seven-DOF robotics research arm, demonstrating that the proposed algorithm yields accurate joint angle trajectories for a given end-effector position/orientation trajectory.
Effects of adaptive refinement on the inverse EEG solution
NASA Astrophysics Data System (ADS)
Weinstein, David M.; Johnson, Christopher R.; Schmidt, John A.
1995-10-01
One of the fundamental problems in electroencephalography can be characterized by an inverse problem. Given a subset of electrostatic potentials measured on the surface of the scalp and the geometry and conductivity properties within the head, calculate the current vectors and potential fields within the cerebrum. Mathematically the generalized EEG problem can be stated as solving Poisson's equation of electrical conduction for the primary current sources. The resulting problem is mathematically ill-posed i.e., the solution does not depend continuously on the data, such that small errors in the measurement of the voltages on the scalp can yield unbounded errors in the solution, and, for the general treatment of a solution of Poisson's equation, the solution is non-unique. However, if accurate solutions the general treatment of a solution of Poisson's equation, the solution is non-unique. However, if accurate solutions to such problems could be obtained, neurologists would gain noninvasive accesss to patient-specific cortical activity. Access to such data would ultimately increase the number of patients who could be effectively treated for pathological cortical conditions such as temporal lobe epilepsy. In this paper, we present the effects of spatial adaptive refinement on the inverse EEG problem and show that the use of adaptive methods allow for significantly better estimates of electric and potential fileds within the brain through an inverse procedure. To test these methods, we have constructed several finite element head models from magneteic resonance images of a patient. The finite element meshes ranged in size from 2724 nodes and 12,812 elements to 5224 nodes and 29,135 tetrahedral elements, depending on the level of discretization. We show that an adaptive meshing algorithm minimizes the error in the forward problem due to spatial discretization and thus increases the accuracy of the inverse solution.
Adaptive critics for dynamic optimization.
Kulkarni, Raghavendra V; Venayagamoorthy, Ganesh Kumar
2010-06-01
A novel action-dependent adaptive critic design (ACD) is developed for dynamic optimization. The proposed combination of a particle swarm optimization-based actor and a neural network critic is demonstrated through dynamic sleep scheduling of wireless sensor motes for wildlife monitoring. The objective of the sleep scheduler is to dynamically adapt the sleep duration to node's battery capacity and movement pattern of animals in its environment in order to obtain snapshots of the animal on its trajectory uniformly. Simulation results show that the sleep time of the node determined by the actor critic yields superior quality of sensory data acquisition and enhanced node longevity. PMID:20223635
Direct inversion of rigid-body rotational dynamics
NASA Technical Reports Server (NTRS)
Bach, Ralph; Paielli, Russell
1990-01-01
The global linearization (inversion) of rigid-body rotational dynamics is reviewed and representations in terms of quaternions and direction cosines are compared. Certain properties common to quaternions and direction cosines that make their use preferable to Euler angles and that simplify the inversion procedure are described. Applications of the inversion procedure for state estimation and attitude control are discussed. To avoid complexities caused by aerodynamics, an example of direct inversion for linear feedback control of spacecraft attitude is given.
Adaptive dynamics of saturated polymorphisms.
Kisdi, Éva; Geritz, Stefan A H
2016-03-01
We study the joint adaptive dynamics of n scalar-valued strategies in ecosystems where n is the maximum number of coexisting strategies permitted by the (generalized) competitive exclusion principle. The adaptive dynamics of such saturated systems exhibits special characteristics, which we first demonstrate in a simple example of a host-pathogen-predator model. The main part of the paper characterizes the adaptive dynamics of saturated polymorphisms in general. In order to investigate convergence stability, we give a new sufficient condition for absolute stability of an arbitrary (not necessarily saturated) polymorphic singularity and show that saturated evolutionarily stable polymorphisms satisfy it. For the case [Formula: see text], we also introduce a method to construct different pairwise invasibility plots of the monomorphic population without changing the selection gradients of the saturated dimorphism. PMID:26676357
Effects of Tape and Exercise on Dynamic Ankle Inversion
Ricard, Mark D.; Sherwood, Stephen M.; Schulthies, Shane S.; Knight, Kenneth L.
2000-01-01
Objective: To compare the effects of tape, with and without prewrap, on dynamic ankle inversion before and after exercise. Design and Setting: Doubly multivariate analyses of variance were used to compare the taping and exercise conditions. Subjects were randomly assigned to a fixed treatment order as determined by a balanced latin square. The independent variables were tape application (no tape, tape with prewrap, tape to skin) and exercise (before and after). The dependent variables were average inversion velocity, total inversion, maximum inversion velocity, and time to maximum inversion. Subjects: Thirty college-age male and female students (17 males, 13 females; mean age = 24.9 ± 4.3 years, range, 19 to 39 years) were tested. Subjects were excluded from the study if they exhibited a painful gait or painful range of motion or had a past history of ankle surgery or an ankle sprain within the past 4 weeks. Measurements: We collected data using electronic goniometers while subjects balanced on the right leg on an inversion platform tilted about the medial-lateral axis to produce 15° of plantar flexion. Sudden ankle inversion was induced by pulling the inversion platform support, allowing the platform support base to rotate 37°. Ten satisfactory trials were recorded on the inversion platform before and after a prescribed exercise bout. We calculated total inversion, time to maximum inversion, average inversion velocity, and maximum inversion velocity after sudden inversion. Results: We found no significant differences between taping to the skin and taping over prewrap for any of the variables measured. There were significant differences between both taping conditions and no-tape postexercise for average inversion velocity, maximum inversion, maximum inversion velocity, and time to maximum inversion. The total inversion mean for no-tape postexercise was 38.8° ± 6.3°, whereas the means for tape and skin and for tape and prewrap were 28.3° ± 4.6° and 29.1°
Parallel computation of manipulator inverse dynamics
NASA Technical Reports Server (NTRS)
Fijany, Amir; Bejczy, Antal K.
1991-01-01
In this article, parallel computation of manipulator inverse dynamics is investigated. A hierarchical graph-based mapping approach is devised to analyze the inherent parallelism in the Newton-Euler formulation at several computational levels, and to derive the features of an abstract architecture for exploitation of parallelism. At each level, a parallel algorithm represents the application of a parallel model of computation that transforms the computation into a graph whose structure defines the features of an abstract architecture, i.e., number of processors, communication structure, etc. Data-flow analysis is employed to derive the time lower bound in the computation as well as the sequencing of the abstract architecture. The features of the target architecture are defined by optimization of the abstract architecture to exploit maximum parallelism while minimizing architectural complexity. An architecture is designed and implemented that is capable of efficient exploitation of parallelism at several computational levels. The computation time of the Newton-Euler formulation for a 6-degree-of-freedom (dof) general manipulator is measured as 187 microsec. The increase in computation time for each additional dof is 23 microsec, which leads to a computation time of less than 500 microsec, even for a 12-dof redundant arm.
Criticality of Adaptive Control Dynamics
NASA Astrophysics Data System (ADS)
Patzelt, Felix; Pawelzik, Klaus
2011-12-01
We show, that stabilization of a dynamical system can annihilate observable information about its structure. This mechanism induces critical points as attractors in locally adaptive control. It also reveals, that previously reported criticality in simple controllers is caused by adaptation and not by other controller details. We apply these results to a real-system example: human balancing behavior. A model of predictive adaptive closed-loop control subject to some realistic constraints is introduced and shown to reproduce experimental observations in unprecedented detail. Our results suggests, that observed error distributions in between the Lévy and Gaussian regimes may reflect a nearly optimal compromise between the elimination of random local trends and rare large errors.
Pegueroles, Cinta; Ferrés-Coy, Albert; Martí-Solano, Maria; Aquadro, Charles F; Pascual, Marta; Mestres, Francesc
2016-01-01
Adaptation is defined as an evolutionary process allowing organisms to succeed in certain habitats or conditions. Chromosomal inversions have the potential to be key in the adaptation processes, since they can contribute to the maintenance of favoured combinations of adaptive alleles through reduced recombination between individuals carrying different inversions. We have analysed six genes (Pif1A, Abi, Sqd, Yrt, Atpα and Fmr1), located inside and outside three inversions of the O chromosome in European populations of Drosophila subobscura. Genetic differentiation was significant between inversions despite extensive recombination inside inverted regions, irrespective of gene distance to the inversion breakpoints. Surprisingly, the highest level of genetic differentiation between arrangements was found for the Atpα gene, which is located outside the O1 and O7 inversions. Two derived unrelated arrangements (O3+4+1 and O3+4+7) are nearly fixed for several amino acid substitutions at the Atpα gene that have been described to confer resistance in other species to the cardenolide ouabain, a plant toxin capable of blocking ATPases. Similarities in the Atpα variants, conferring ouabain resistance in both arrangements, may be the result of convergent substitution and be favoured in response to selective pressures presumably related to the presence of plants containing ouabain in the geographic locations where both inversions are present. PMID:27029337
A model following inverse controller with adaptive compensation for General Aviation aircraft
NASA Astrophysics Data System (ADS)
Bruner, Hugh S.
The theory for an adaptive inverse flight controller, suitable for use on General Aviation aircraft, is developed in this research. The objectives of this controller are to separate the normally coupled modes of the basic aircraft and thereby permit direct control of airspeed and flight-path angle, meet prescribed performance characteristics as defined by damping ratio and natural frequency, adapt to uncertainties in the physical plant, and be computationally efficient. The three basic elements of the controller are a linear prefilter, an inverse transfer function, and an adaptive neural network compensator. The linear prefilter shapes accelerations required of the overall system in order to achieve the desired system performance characteristics. The inverse transfer function is used to compute the aircraft control inputs required to achieve the necessary accelerations. The adaptive neural network compensator is used to compensate for modeling errors during design or real-time changes in the physical plant. This architecture is patterned after the work of Calise, but differs by not requiring dynamic feedback of the state variables. The controller is coded in ANSI C and integrated with a simulation of a typical General Aviation aircraft. Twenty-three cases are simulated to prove that the objectives for the controller are met. Among these cases are simulated stability and controllability failures in the physical plant, as well as several simulated failures of the neural network. With the exception of some bounded speed-tracking error, the controller is capable of continued flight with any foreseeable failure of the neural network. Recommendations are provided for follow-on investigations by other researchers.
Dynamic Adaption of Vascular Morphology
Okkels, Fridolin; Jacobsen, Jens Christian Brings
2012-01-01
The structure of vascular networks adapts continuously to meet changes in demand of the surrounding tissue. Most of the known vascular adaptation mechanisms are based on local reactions to local stimuli such as pressure and flow, which in turn reflects influence from the surrounding tissue. Here we present a simple two-dimensional model in which, as an alternative approach, the tissue is modeled as a porous medium with intervening sharply defined flow channels. Based on simple, physiologically realistic assumptions, flow-channel structure adapts so as to reach a configuration in which all parts of the tissue are supplied. A set of model parameters uniquely determine the model dynamics, and we have identified the region of the best-performing model parameters (a global optimum). This region is surrounded in parameter space by less optimal model parameter values, and this separation is characterized by steep gradients in the related fitness landscape. Hence it appears that the optimal set of parameters tends to localize close to critical transition zones. Consequently, while the optimal solution is stable for modest parameter perturbations, larger perturbations may cause a profound and permanent shift in systems characteristics. We suggest that the system is driven toward a critical state as a consequence of the ongoing parameter optimization, mimicking an evolutionary pressure on the system. PMID:23060814
Yan Di; Liang Jian
2013-02-15
Purpose: To construct expected treatment dose for adaptive inverse planning optimization, and evaluate it on head and neck (h and n) cancer adaptive treatment modification. Methods: Adaptive inverse planning engine was developed and integrated in our in-house adaptive treatment control system. The adaptive inverse planning engine includes an expected treatment dose constructed using the daily cone beam (CB) CT images in its objective and constrains. Feasibility of the adaptive inverse planning optimization was evaluated retrospectively using daily CBCT images obtained from the image guided IMRT treatment of 19 h and n cancer patients. Adaptive treatment modification strategies with respect to the time and the number of adaptive inverse planning optimization during the treatment course were evaluated using the cumulative treatment dose in organs of interest constructed using all daily CBCT images. Results: Expected treatment dose was constructed to include both the delivered dose, to date, and the estimated dose for the remaining treatment during the adaptive treatment course. It was used in treatment evaluation, as well as in constructing the objective and constraints for adaptive inverse planning optimization. The optimization engine is feasible to perform planning optimization based on preassigned treatment modification schedule. Compared to the conventional IMRT, the adaptive treatment for h and n cancer illustrated clear dose-volume improvement for all critical normal organs. The dose-volume reductions of right and left parotid glands, spine cord, brain stem and mandible were (17 {+-} 6)%, (14 {+-} 6)%, (11 {+-} 6)%, (12 {+-} 8)%, and (5 {+-} 3)% respectively with the single adaptive modification performed after the second treatment week; (24 {+-} 6)%, (22 {+-} 8)%, (21 {+-} 5)%, (19 {+-} 8)%, and (10 {+-} 6)% with three weekly modifications; and (28 {+-} 5)%, (25 {+-} 9)%, (26 {+-} 5)%, (24 {+-} 8)%, and (15 {+-} 9)% with five weekly modifications. Conclusions
Dynamic Inversion based Control of a Docking Mechanism
NASA Technical Reports Server (NTRS)
Kulkarni, Nilesh V.; Ippolito, Corey; Krishnakumar, Kalmanje
2006-01-01
The problem of position and attitude control of the Stewart platform based docking mechanism is considered motivated by its future application in space missions requiring the autonomous docking capability. The control design is initiated based on the framework of the intelligent flight control architecture being developed at NASA Ames Research Center. In this paper, the baseline position and attitude control system is designed using dynamic inversion with proportional-integral augmentation. The inverse dynamics uses a Newton-Euler formulation that includes the platform dynamics, the dynamics of the individual legs along with viscous friction in the joints. Simulation results are presented using forward dynamics simulated by a commercial physics engine that builds the system as individual elements with appropriate joints and uses constrained numerical integration,
A spatiotemporal dynamic distributed solution to the MEG inverse problem.
Lamus, Camilo; Hämäläinen, Matti S; Temereanca, Simona; Brown, Emery N; Purdon, Patrick L
2012-11-01
MEG/EEG are non-invasive imaging techniques that record brain activity with high temporal resolution. However, estimation of brain source currents from surface recordings requires solving an ill-conditioned inverse problem. Converging lines of evidence in neuroscience, from neuronal network models to resting-state imaging and neurophysiology, suggest that cortical activation is a distributed spatiotemporal dynamic process, supported by both local and long-distance neuroanatomic connections. Because spatiotemporal dynamics of this kind are central to brain physiology, inverse solutions could be improved by incorporating models of these dynamics. In this article, we present a model for cortical activity based on nearest-neighbor autoregression that incorporates local spatiotemporal interactions between distributed sources in a manner consistent with neurophysiology and neuroanatomy. We develop a dynamic maximum a posteriori expectation-maximization (dMAP-EM) source localization algorithm for estimation of cortical sources and model parameters based on the Kalman Filter, the Fixed Interval Smoother, and the EM algorithms. We apply the dMAP-EM algorithm to simulated experiments as well as to human experimental data. Furthermore, we derive expressions to relate our dynamic estimation formulas to those of standard static models, and show how dynamic methods optimally assimilate past and future data. Our results establish the feasibility of spatiotemporal dynamic estimation in large-scale distributed source spaces with several thousand source locations and hundreds of sensors, with resulting inverse solutions that provide substantial performance improvements over static methods. PMID:22155043
A spatiotemporal dynamic distributed solution to the MEG inverse problem
Lamus, Camilo; Hämäläinen, Matti S.; Temereanca, Simona; Brown, Emery N.; Purdon, Patrick L.
2012-01-01
MEG/EEG are non-invasive imaging techniques that record brain activity with high temporal resolution. However, estimation of brain source currents from surface recordings requires solving an ill-conditioned inverse problem. Converging lines of evidence in neuroscience, from neuronal network models to resting-state imaging and neurophysiology, suggest that cortical activation is a distributed spatiotemporal dynamic process, supported by both local and long-distance neuroanatomic connections. Because spatiotemporal dynamics of this kind are central to brain physiology, inverse solutions could be improved by incorporating models of these dynamics. In this article, we present a model for cortical activity based on nearest-neighbor autoregression that incorporates local spatiotemporal interactions between distributed sources in a manner consistent with neurophysiology and neuroanatomy. We develop a dynamic Maximum a Posteriori Expectation-Maximization (dMAP-EM) source localization algorithm for estimation of cortical sources and model parameters based on the Kalman Filter, the Fixed Interval Smoother, and the EM algorithms. We apply the dMAP-EM algorithm to simulated experiments as well as to human experimental data. Furthermore, we derive expressions to relate our dynamic estimation formulas to those of standard static models, and show how dynamic methods optimally assimilate past and future data. Our results establish the feasibility of spatiotemporal dynamic estimation in large-scale distributed source spaces with several thousand source locations and hundreds of sensors, with resulting inverse solutions that provide substantial performance improvements over static methods. PMID:22155043
Adaptation and dynamics of cat retinal ganglion cells
Enroth-Cugell, Christina; Shapley, R. M.
1973-01-01
1. The impulse/quantum (I/Q) ratio was measured as a function of background illumination for rod-dominated, pure central, linear square-wave responses of retinal ganglion cells in the cat. 2. The I/Q ratio was constant at low backgrounds (dark adapted state) and inversely proportional to the 0·9 power of the background at high backgrounds (the light adapted state). There was an abrupt transition from the dark-adapted state to the light-adapted state. 3. It was possible to define the adaptation level at a particular background as the ratio (I/Q ratio at that background)/(dark adapted I/Q ratio). 4. The time course of the square-wave response was correlated with the adaptation level. The response was sustained in the dark-adapted state, partially transient at the transition level, and progressively more transient the lower the impulse/quantum ratio of the ganglion cell became. This was true both for on-centre and off-centre cells. 5. The frequency response of the central response mechanism at different adaptation levels was measured. It was a low-pass characteristic in the dark-adapted state and became progressively more of a bandpass characteristic as the cell became more light-adapted. 6. The rapidity of onset of adaptation was measured with a time-varying adapting light. The impulse/quantum ratio is reset within 100 msec of the onset of the conditioning light, and is kept at the new value throughout the time the conditioning light is on. 7. These results can be explained by a nonlinear feedback model. In the model, it is postulated that the exponential function of the horizontal cell potential controls transmission from rods to bipolars. This model has an abrupt transition from dark- to light-adapted states, and its response dynamics are correlated with adaptation level. PMID:4747229
From seismic images to plate dynamics: Towards the full inverse
NASA Astrophysics Data System (ADS)
Gurnis, M.; Ratnaswamy, V.; Stadler, G.; Ghattas, O.; Alisic, L.
2014-12-01
Three-dimensional seismic images of slabs and other mantle structures provide a first order constraint on the forces driving plate motions. Previous attempts to invert for plate motions from seismic images have blurry slabs that do not act as stress guides. Using forward models, we describe characteristics needed to capture the coupling between mantle structures and plates. In forward models, we capitalized on advances in adaptive mesh refinement and scalable solvers to simulate global mantle flow and plate motions, with plate margins resolved down to 1 km. Cold thermal anomalies within the lower mantle are coupled into oceanic plates through narrow high-viscosity slabs, altering the velocity of oceanic plates. Back-arc extension and slab rollback are emergent consequences of slab descent in the upper mantle. The forward models require the solution of a highly ill-conditioned non-linear Stokes equation. Based on a realistic rheological model with yielding and strain rate weakening from dislocation creep, we formulate inverse problems casted as PDE-constrained optimization problems and derive adjoints of the nonlinear Stokes and incompressibility equations. An inexact-Gauss Newton method is used to infer the rheological parameters while quantifying the uncertainty using the Hessian at the maximum a posteriori (MAP) point. Through 2-D numerical experiments we demonstrate that when the temperature field is known from seismic images, we can recover all of these properties to varying levels of certainty: strength of plate boundaries, yield stress and strain rate exponent in the upper mantle. When the system becomes more unconstrained (when all three mechanical properties are unknown), there can be tradeoffs depending on how well the data approximates the realistic dynamics. As plate boundaries become weaker beyond a limiting value, the uncertainty of the inferred parameters increases due to insensitivity of plate motion to plate coupling. Using the inverse of the
Robust inverse kinematics using damped least squares with dynamic weighting
NASA Technical Reports Server (NTRS)
Schinstock, D. E.; Faddis, T. N.; Greenway, R. B.
1994-01-01
This paper presents a general method for calculating the inverse kinematics with singularity and joint limit robustness for both redundant and non-redundant serial-link manipulators. Damped least squares inverse of the Jacobian is used with dynamic weighting matrices in approximating the solution. This reduces specific joint differential vectors. The algorithm gives an exact solution away from the singularities and joint limits, and an approximate solution at or near the singularities and/or joint limits. The procedure is here implemented for a six d.o.f. teleoperator and a well behaved slave manipulator resulted under teleoperational control.
Efficient algorithms for linear dynamic inverse problems with known motion
NASA Astrophysics Data System (ADS)
Hahn, B. N.
2014-03-01
An inverse problem is called dynamic if the object changes during the data acquisition process. This occurs e.g. in medical applications when fast moving organs like the lungs or the heart are imaged. Most regularization methods are based on the assumption that the object is static during the measuring procedure. Hence, their application in the dynamic case often leads to serious motion artefacts in the reconstruction. Therefore, an algorithm has to take into account the temporal changes of the investigated object. In this paper, a reconstruction method that compensates for the motion of the object is derived for dynamic linear inverse problems. The algorithm is validated at numerical examples from computerized tomography.
Adaptive EAGLE dynamic solution adaptation and grid quality enhancement
NASA Technical Reports Server (NTRS)
Luong, Phu Vinh; Thompson, J. F.; Gatlin, B.; Mastin, C. W.; Kim, H. J.
1992-01-01
In the effort described here, the elliptic grid generation procedure in the EAGLE grid code was separated from the main code into a subroutine, and a new subroutine which evaluates several grid quality measures at each grid point was added. The elliptic grid routine can now be called, either by a computational fluid dynamics (CFD) code to generate a new adaptive grid based on flow variables and quality measures through multiple adaptation, or by the EAGLE main code to generate a grid based on quality measure variables through static adaptation. Arrays of flow variables can be read into the EAGLE grid code for use in static adaptation as well. These major changes in the EAGLE adaptive grid system make it easier to convert any CFD code that operates on a block-structured grid (or single-block grid) into a multiple adaptive code.
NASA Astrophysics Data System (ADS)
Yang, Dikun; Oldenburg, Douglas W.; Haber, Eldad
2014-03-01
Airborne electromagnetic (AEM) methods are highly efficient tools for assessing the Earth's conductivity structures in a large area at low cost. However, the configuration of AEM measurements, which typically have widely distributed transmitter-receiver pairs, makes the rigorous modelling and interpretation extremely time-consuming in 3-D. Excessive overcomputing can occur when working on a large mesh covering the entire survey area and inverting all soundings in the data set. We propose two improvements. The first is to use a locally optimized mesh for each AEM sounding for the forward modelling and calculation of sensitivity. This dedicated local mesh is small with fine cells near the sounding location and coarse cells far away in accordance with EM diffusion and the geometric decay of the signals. Once the forward problem is solved on the local meshes, the sensitivity for the inversion on the global mesh is available through quick interpolation. Using local meshes for AEM forward modelling avoids unnecessary computing on fine cells on a global mesh that are far away from the sounding location. Since local meshes are highly independent, the forward modelling can be efficiently parallelized over an array of processors. The second improvement is random and dynamic down-sampling of the soundings. Each inversion iteration only uses a random subset of the soundings, and the subset is reselected for every iteration. The number of soundings in the random subset, determined by an adaptive algorithm, is tied to the degree of model regularization. This minimizes the overcomputing caused by working with redundant soundings. Our methods are compared against conventional methods and tested with a synthetic example. We also invert a field data set that was previously considered to be too large to be practically inverted in 3-D. These examples show that our methodology can dramatically reduce the processing time of 3-D inversion to a practical level without losing resolution
Adaptive Role of Inversion Polymorphism of Drosophila subobscura in Lead Stressed Environment
Kenig, Bojan; Kurbalija Novičić, Zorana; Patenković, Aleksandra; Stamenković-Radak, Marina; Anđelković, Marko
2015-01-01
Local adaptation to environmental stress at different levels of genetic polymorphism in various plants and animals has been documented through evolution of heavy metal tolerance. We used samples of Drosophila subobscura populations from two differently polluted environments to analyze the change of chromosomal inversion polymorphism as genetic marker during laboratory exposure to lead. Exposure to environmental contamination can affect the genetic content within a particular inversion and produce targets for selection in populations from different environments. The aims were to discover whether the inversion polymorphism is shaped by the local natural environments, and if lead as a selection pressure would cause adaptive divergence of two populations during the multigenerational laboratory experiment. The results showed that populations retain signatures from past contamination events, and that heavy metal pollution can cause adaptive changes in population. Differences in inversion polymorphism between the two populations increased over generations under lead contamination in the laboratory. The inversion polymorphism of population originating from the more polluted natural environment was more stable during the experiment, both under conditions with and without lead. Therefore, results showed that inversion polymorphism as a genetic marker reflects a strong signature of adaptation to the local environment, and that historical demographic events and selection are important for both prediction of evolutionary potential and long-term viability of natural populations. PMID:26102201
Dynamic optimization and adaptive controller design
NASA Astrophysics Data System (ADS)
Inamdar, S. R.
2010-10-01
In this work I present a new type of controller which is an adaptive tracking controller which employs dynamic optimization for optimizing current value of controller action for the temperature control of nonisothermal continuously stirred tank reactor (CSTR). We begin with a two-state model of nonisothermal CSTR which are mass and heat balance equations and then add cooling system dynamics to eliminate input multiplicity. The initial design value is obtained using local stability of steady states where approach temperature for cooling action is specified as a steady state and a design specification. Later we make a correction in the dynamics where material balance is manipulated to use feed concentration as a system parameter as an adaptive control measure in order to avoid actuator saturation for the main control loop. The analysis leading to design of dynamic optimization based parameter adaptive controller is presented. The important component of this mathematical framework is reference trajectory generation to form an adaptive control measure.
GARCH modelling of covariance in dynamical estimation of inverse solutions
NASA Astrophysics Data System (ADS)
Galka, Andreas; Yamashita, Okito; Ozaki, Tohru
2004-12-01
The problem of estimating unobserved states of spatially extended dynamical systems poses an inverse problem, which can be solved approximately by a recently developed variant of Kalman filtering; in order to provide the model of the dynamics with more flexibility with respect to space and time, we suggest to combine the concept of GARCH modelling of covariance, well known in econometrics, with Kalman filtering. We formulate this algorithm for spatiotemporal systems governed by stochastic diffusion equations and demonstrate its feasibility by presenting a numerical simulation designed to imitate the situation of the generation of electroencephalographic recordings by the human cortex.
An adaptive subspace trust-region method for frequency-domain seismic full waveform inversion
NASA Astrophysics Data System (ADS)
Zhang, Huan; Li, Xiaofan; Song, Hanjie; Liu, Shaolin
2015-05-01
Full waveform inversion is currently considered as a promising seismic imaging method to obtain high-resolution and quantitative images of the subsurface. It is a nonlinear ill-posed inverse problem, the main difficulty of which that prevents the full waveform inversion from widespread applying to real data is the sensitivity to incorrect initial models and noisy data. Local optimization theories including Newton's method and gradient method always lead the convergence to local minima, while global optimization algorithms such as simulated annealing are computationally costly. To confront this issue, in this paper we investigate the possibility of applying the trust-region method to the full waveform inversion problem. Different from line search methods, trust-region methods force the new trial step within a certain neighborhood of the current iterate point. Theoretically, the trust-region methods are reliable and robust, and they have very strong convergence properties. The capability of this inversion technique is tested with the synthetic Marmousi velocity model and the SEG/EAGE Salt model. Numerical examples demonstrate that the adaptive subspace trust-region method can provide solutions closer to the global minima compared to the conventional Approximate Hessian approach and the L-BFGS method with a higher convergence rate. In addition, the match between the inverted model and the true model is still excellent even when the initial model deviates far from the true model. Inversion results with noisy data also exhibit the remarkable capability of the adaptive subspace trust-region method for low signal-to-noise data inversions. Promising numerical results suggest this adaptive subspace trust-region method is suitable for full waveform inversion, as it has stronger convergence and higher convergence rate.
Impulse radar imaging for dispersive concrete using inverse adaptive filtering techniques
Arellano, J.; Hernandez, J.M.; Brase, J.
1993-05-01
This publication addresses applications of a delayed inverse model adaptive filter for modeled data obtained from short-pulse radar reflectometry. To determine the integrity of concrete, a digital adaptive filter was used, which allows compensation of dispersion and clutter generated by the concrete. A standard set of weights produced by an adaptive filter are used on modeled data to obtain the inverse-impulse response of the concrete. The data for this report include: Multiple target, nondispersive data; single-target, variable-size dispersive data; single-target, variable-depth dispersive data; and single-target, variable transmitted-pulse-width dispersive data. Results of this simulation indicate that data generated by the weights of the adaptive filter, coupled with a two-dimensional, synthetic-aperture focusing technique, successfully generate two-dimensional images of targets within the concrete from modeled data.
Dynamic adaptivity of "smart" piezoelectric structures
NASA Astrophysics Data System (ADS)
Tzou, Horn-Sen; Zhong, Jianping P.
1990-10-01
Active smart" space and machine structures with adaptive dynamic characteristics have long been interested in a variety of high-performance systems, e.g., flexible robots, flexible space structures, "smart" machines, etc. In this paper, an active adaptive structure made of piezoelectric materials is proposed and evaluated. The structural adaptivity is achieved by a voltage feedback (open or closed loops) utilizing the converse piezoelectric effect. A mathematical model is proposed and the electrodynamic equations of motion and the generalized boundary conditions of a generic piezoelectric shell subjected to mechanical and electrical excitations are derived using Hamilton's principle and the linear piezoelectric theory. The dynamic adaptivity of the structure is introduced using a feedback control system. The theory is demonstrated in a case study in which the structural adaptivity (natural frequency) is investigated.
Adaptive Thouless-Anderson-Palmer approach to inverse Ising problems with quenched random fields.
Huang, Haiping; Kabashima, Yoshiyuki
2013-06-01
The adaptive Thouless-Anderson-Palmer equation is derived for inverse Ising problems in the presence of quenched random fields. We test the proposed scheme on Sherrington-Kirkpatrick, Hopfield, and random orthogonal models and find that the adaptive Thouless-Anderson-Palmer approach allows accurate inference of quenched random fields whose distribution can be either Gaussian or bimodal. In particular, another competitive method for inferring external fields, namely, the naive mean field method with diagonal weights, is compared and discussed. PMID:23848649
ERIC Educational Resources Information Center
Brown, Malcolm
2009-01-01
Inversions are fascinating phenomena. They are reversals of the normal or expected order. They occur across a wide variety of contexts. What do inversions have to do with learning spaces? The author suggests that they are a useful metaphor for the process that is unfolding in higher education with respect to education. On the basis of…
Solving the inverse problem of noise-driven dynamic networks.
Zhang, Zhaoyang; Zheng, Zhigang; Niu, Haijing; Mi, Yuanyuan; Wu, Si; Hu, Gang
2015-01-01
Nowadays, massive amounts of data are available for analysis in natural and social systems and the tasks to depict system structures from the data, i.e., the inverse problems, become one of the central issues in wide interdisciplinary fields. In this paper, we study the inverse problem of dynamic complex networks driven by white noise. A simple and universal inference formula of double correlation matrices and noise-decorrelation (DCMND) method is derived analytically, and numerical simulations confirm that the DCMND method can accurately depict both network structures and noise correlations by using available output data only. This inference performance has never been regarded possible by theoretical derivation, numerical computation, and experimental design. PMID:25679664
Solving the inverse problem of noise-driven dynamic networks
NASA Astrophysics Data System (ADS)
Zhang, Zhaoyang; Zheng, Zhigang; Niu, Haijing; Mi, Yuanyuan; Wu, Si; Hu, Gang
2015-01-01
Nowadays, massive amounts of data are available for analysis in natural and social systems and the tasks to depict system structures from the data, i.e., the inverse problems, become one of the central issues in wide interdisciplinary fields. In this paper, we study the inverse problem of dynamic complex networks driven by white noise. A simple and universal inference formula of double correlation matrices and noise-decorrelation (DCMND) method is derived analytically, and numerical simulations confirm that the DCMND method can accurately depict both network structures and noise correlations by using available output data only. This inference performance has never been regarded possible by theoretical derivation, numerical computation, and experimental design.
Geophysical Inversion with Adaptive Array Processing of Ambient Noise
NASA Astrophysics Data System (ADS)
Traer, James
2011-12-01
Land-based seismic observations of microseisms generated during Tropical Storms Ernesto and Florence are dominated by signals in the 0.15--0.5Hz band. Data from seafloor hydrophones in shallow water (70m depth, 130 km off the New Jersey coast) show dominant signals in the gravity-wave frequency band, 0.02--0.18Hz and low amplitudes from 0.18--0.3Hz, suggesting significant opposing wave components necessary for DF microseism generation were negligible at the site. Both storms produced similar spectra, despite differing sizes, suggesting near-coastal shallow water as the dominant region for observed microseism generation. A mathematical explanation for a sign-inversion induced to the passive fathometer response by minimum variance distortionless response (MVDR) beamforming is presented. This shows that, in the region containing the bottom reflection, the MVDR fathometer response is identical to that obtained with conventional processing multiplied by a negative factor. A model is presented for the complete passive fathometer response to ocean surface noise, interfering discrete noise sources, and locally uncorrelated noise in an ideal waveguide. The leading order term of the ocean surface noise produces the cross-correlation of vertical multipaths and yields the depth of sub-bottom reflectors. Discrete noise incident on the array via multipaths give multiple peaks in the fathometer response. These peaks may obscure the sub-bottom reflections but can be attenuated with use of Minimum Variance Distortionless Response (MVDR) steering vectors. A theory is presented for the Signal-to-Noise-Ratio (SNR) for the seabed reflection peak in the passive fathometer response as a function of seabed depth, seabed reflection coefficient, averaging time, bandwidth and spatial directivity of the noise field. The passive fathometer algorithm was applied to data from two drifting array experiments in the Mediterranean, Boundary 2003 and 2004, with 0.34s of averaging time. In the 2004
Doss, S D; Ezzedine, S; Gelinas, R; Chawathe, A
2001-06-11
A novel approach called Forward-Inverse Adaptive Techniques (FIAT) for reservoir characterization is developed and applied to three representative exploration cases. Inverse modeling refers to the determination of the entire reservoir permeability under steady state single-phase flow regime, given only field permeability, pressure and production well measurements. FIAT solves the forward and inverse partial differential equations (PDEs) simultaneously by adding a regularization term and filtering pressure gradients. An implicit adaptive-grid, Galerkin, numerical scheme is used to numerically solve the set of PDEs subject to pressure and permeability boundary conditions. Three examples are presented. Results from all three cases demonstrate attainable and reasonably accurate solutions and, more importantly, provide insights into the consequences of data undersampling.
Dynamical similarities of the direct and inverse turbulent cascades
NASA Astrophysics Data System (ADS)
Vela-Martin, Alberto; Jimenez, Javier
2015-11-01
A fully reversible homogeneous isotropic turbulent system is constructed using inviscid LES to model energy fluxes in the inertial range. It recovers energy and other turbulent quantities when reversed after being allowed to decay. During the first phase, a direct cascade transfers energy from large to small scales while, during the second, an inverse cascade does the opposite. Short-time Lyapunov (STL) analysis is used to study and compare the dynamics of both cascades. This allows us to identify a smallest length scale for the chaotic flow behavior, below which the system behaves as a unit dynamically enslaved to larger motions by the contracting effect of the model. Above it, the inertial forces become relevant and the system is fully chaotic. When the inertial scales are isolated, the leading STL exponent is similar for both cascades, suggesting that the dynamics of the inertial range is conservative and time-symmetric, and that the direct and inverse energy cascades share similar energy transfer mechanisms. The cascade would thus be a bi-directional reversible process with similar up and down mechanisms, although, because the L2 norm used in the STL analysis respects the geometry of phase space, the entropy-driven cascade directionally is not precluded. Funded by the ERC Multiflow program.
Goal Directed Model Inversion: A Study of Dynamic Behavior
NASA Technical Reports Server (NTRS)
Colombano, Silvano P.; Compton, Michael; Raghavan, Bharathi; Lum, Henry, Jr. (Technical Monitor)
1994-01-01
Goal Directed Model Inversion (GDMI) is an algorithm designed to generalize supervised learning to the case where target outputs are not available to the learning system. The output of the learning system becomes the input to some external device or transformation, and only the output of this device or transformation can be compared to a desired target. The fundamental driving mechanism of GDMI is to learn from success. Given that a wrong outcome is achieved, one notes that the action that produced that outcome 0 "would have been right if the outcome had been the desired one." The algorithm then proceeds as follows: (1) store the action that produced the wrong outcome as a "target" (2) redefine the wrong outcome as a desired goal (3) submit the new desired goal to the system (4) compare the new action with the target action and modify the system by using a suitable algorithm for credit assignment (Back propagation in our example) (5) resubmit the original goal. Prior publications by our group in this area focused on demonstrating empirical results based on the inverse kinematic problem for a simulated robotic arm. In this paper we apply the inversion process to much simpler analytic functions in order to elucidate the dynamic behavior of the system and to determine the sensitivity of the learning process to various parameters. This understanding will be necessary for the acceptance of GDMI as a practical tool.
Bridi, L C; Rafael, M S
2016-02-01
Anopheles darlingi is the main malaria vector in humans in South America. In the Amazon basin, it lives along the banks of rivers and lakes, which responds to the annual hydrological cycle (dry season and rainy season). In these breeding sites, the larvae of this mosquito feed on decomposing organic and microorganisms, which can be pathogenic and trigger the activation of innate immune system pathways, such as proteins Gram-negative binding protein (GNBP). Such environmental changes affect the occurrence of polymorphic inversions especially at the heterozygote frequency, which confer adaptative advantage compared to homozygous inversions. We mapped the GNBP probe to the An. darlingi 2Rd inversion by fluorescent in situ hybridization (FISH), which was a good indicator of the GNBP immune response related to the chromosomal polymorphic inversions and adaptative evolution. To better understand the evolutionary relations and time of divergence of the GNBP of An. darlingi, we compared it with nine other mosquito GNBPs. The results of the phylogenetic analysis of the GNBP sequence between the species of mosquitoes demonstrated three clades. Clade I and II included the GNBPB5 sequence, and clade III the sequence of GNBPB1. Most of these sequences of GNBP analyzed were homologous with that of subfamily B, including that of An. gambiae (87 %), therefore suggesting that GNBP of An. darling belongs to subfamily B. This work helps us understand the role of inversion polymorphism in evolution of An. darlingi. PMID:26767379
Inverse problem of nonlinear dynamical systems: a constructive approach
Gonzalez-Gascon, F.; Moreno-Insertis, F.; Rodriguez-Camino, E.
1980-08-01
A quite simple and practical method is developed for the construction of two dimensional nonlinear dynamical systems (plane vector fields) possessing an arbitrary number of given limit cycles. The method is applied to the construction of n-dimensional dynamical systems (R/sup n/ vector fields) possessing at least one limit cycle and, under certain circumstances, more than one, or even a numerable infinity. Interesting open problems arise when n is greater than two, or where more than one limit cycle appears. Our constructive algorithm for this type of inverse problem is also applied to the construction of second order differential equations (Newtonian differential equations) possessing a finite or infinite number of invariant speeds. This last problem is relevant for certain aspects of the special theory of relativity.
Dettmer, Jan; Dosso, Stan E
2013-05-01
This paper develops a probabilistic two-dimensional (2D) inversion for geoacoustic seabed and water-column parameters in a strongly range-dependent environment. Range-dependent environments in shelf and shelf-break regions are of increasing importance to the acoustical-oceanography community, and recent advances in nonlinear inverse theory and sampling methods are applied here for efficient probabilistic range-dependent inversion. The 2D seabed and water column are parameterized using highly efficient, self-adapting irregular grids which intrinsically match the local resolving power of the data and provide parsimonious solutions requiring few parameters to capture complex environments. The self-adapting parameterization is achieved by implementing the irregular grid as a trans-dimensional hierarchical Bayesian model with an unknown number of nodes which is sampled with the Metropolis-Hastings-Green algorithm. To improve sampling, population Monte Carlo is applied with a large number of interacting parallel Markov chains with adaptive proposal distributions. The inversion is applied to simulated data for a vertical-line array and several source locations to several kilometers range. Complex acoustic-pressure fields are computed using a parabolic equation model and results are considered in terms of 2D ensemble parameter estimates and credibility intervals. PMID:23654369
Adaptive inverse control for rotorcraft vibration reduction. Ph.D. Thesis
NASA Technical Reports Server (NTRS)
Jacklin, S. A.
1985-01-01
The Least Mean Square (LMS) algorithm is extended to solve the multiple-input, multiple-output problem of alleviating N/Rev helicopter fuselage vibration by means of adaptive inverse control. A frequency domain locally linear model is used to represent the transfer matrix relating the high harmonic pitch control inputs to the harmonic vibration outputs to be controlled. By using the inverse matrix as the controller gain matrix, an adaptive inverse regulator is formed to alleviate the N/Rev vibration. The stability and rate of convergence properties of the extended LMS algorithm are discussed. It is shown that the stability ranges for the elements of the stability gain matrix are directly related to the eigenvalues of the vibration signal information matrix for the learning phase, but not for the control phase. The overall conclusion is that the LMS adaptive inverse control method can form a robust vibration control system, but will require some tuning of the input sensor gains, the stability gain matrix, and the amount of control relaxation to be used. The learning curve of the controller during the learning phase is shown to be quantitatively close to that predicted by averaging the learning curves of the normal modes. It is shown that the best selections of the stability gain matrix elements and the amount of control relaxation is basically a compromise between slow, stable convergence and fast convergence with increased possibility of unstable identification.
Success Stories in Control: Nonlinear Dynamic Inversion Control
NASA Technical Reports Server (NTRS)
Bosworth, John T.
2010-01-01
NASA plays an important role in advancing the state of the art in flight control systems. In the case of Nonlinear Dynamic Inversion (NDI) NASA supported initial implementation of the theory in an aircraft and demonstration in a space vehicle. Dr. Dale Enns of Honeywell Aerospace Advanced Technology performed this work in cooperation with NASA and under NASA contract. Honeywell and Lockheed Martin were subsequently contracted by AFRL to create "Design Guidelines for Multivariable Control Theory". This foundational work directly contributed to the advancement of the technology and the credibility of the control law as a design option. As a result Honeywell collaborated with Lockheed Martin to produce a Nonlinear Dynamic Inversion controller for the X-35 and subsequently Lockheed Martin did the same for the production Lockheed Martin F-35 vehicle. The theory behind NDI is to use a systematic generalized approach to controlling a vehicle. Using general aircraft nonlinear equations of motion and onboard aerodynamic, mass properties, and engine models specific to the vehicle, a relationship between control effectors and desired aircraft motion can be formulated. Using this formulation a control combination is used that provides a predictable response to commanded motion. Control loops around this formulation shape the response as desired and provide robustness to modeling errors. Once the control law is designed it can be used on a similar class of vehicle with only an update to the vehicle specific onboard models.
Automated adaptive inference of phenomenological dynamical models
Daniels, Bryan C.; Nemenman, Ilya
2015-01-01
Dynamics of complex systems is often driven by large and intricate networks of microscopic interactions, whose sheer size obfuscates understanding. With limited experimental data, many parameters of such dynamics are unknown, and thus detailed, mechanistic models risk overfitting and making faulty predictions. At the other extreme, simple ad hoc models often miss defining features of the underlying systems. Here we develop an approach that instead constructs phenomenological, coarse-grained models of network dynamics that automatically adapt their complexity to the available data. Such adaptive models produce accurate predictions even when microscopic details are unknown. The approach is computationally tractable, even for a relatively large number of dynamical variables. Using simulated data, it correctly infers the phase space structure for planetary motion, avoids overfitting in a biological signalling system and produces accurate predictions for yeast glycolysis with tens of data points and over half of the interacting species unobserved. PMID:26293508
Automated adaptive inference of phenomenological dynamical models
NASA Astrophysics Data System (ADS)
Daniels, Bryan C.; Nemenman, Ilya
2015-08-01
Dynamics of complex systems is often driven by large and intricate networks of microscopic interactions, whose sheer size obfuscates understanding. With limited experimental data, many parameters of such dynamics are unknown, and thus detailed, mechanistic models risk overfitting and making faulty predictions. At the other extreme, simple ad hoc models often miss defining features of the underlying systems. Here we develop an approach that instead constructs phenomenological, coarse-grained models of network dynamics that automatically adapt their complexity to the available data. Such adaptive models produce accurate predictions even when microscopic details are unknown. The approach is computationally tractable, even for a relatively large number of dynamical variables. Using simulated data, it correctly infers the phase space structure for planetary motion, avoids overfitting in a biological signalling system and produces accurate predictions for yeast glycolysis with tens of data points and over half of the interacting species unobserved.
NASA Astrophysics Data System (ADS)
Ry, Rexha Verdhora; Nugraha, Andri Dian
2015-04-01
Observation of earthquakes is routinely used widely in tectonic activity observation, and also in local scale such as volcano tectonic and geothermal activity observation. It is necessary for determining the location of precise hypocenter which the process involves finding a hypocenter location that has minimum error between the observed and the calculated travel times. When solving this nonlinear inverse problem, simulated annealing inversion method can be applied to such global optimization problems, which the convergence of its solution is independent of the initial model. In this study, we developed own program codeby applying adaptive simulated annealing inversion in Matlab environment. We applied this method to determine earthquake hypocenter using several data cases which are regional tectonic, volcano tectonic, and geothermal field. The travel times were calculated using ray tracing shooting method. We then compared its results with the results using Geiger's method to analyze its reliability. Our results show hypocenter location has smaller RMS error compared to the Geiger's result that can be statistically associated with better solution. The hypocenter of earthquakes also well correlated with geological structure in the study area. Werecommend using adaptive simulated annealing inversion to relocate hypocenter location in purpose to get precise and accurate earthquake location.
Ry, Rexha Verdhora; Nugraha, Andri Dian
2015-04-24
Observation of earthquakes is routinely used widely in tectonic activity observation, and also in local scale such as volcano tectonic and geothermal activity observation. It is necessary for determining the location of precise hypocenter which the process involves finding a hypocenter location that has minimum error between the observed and the calculated travel times. When solving this nonlinear inverse problem, simulated annealing inversion method can be applied to such global optimization problems, which the convergence of its solution is independent of the initial model. In this study, we developed own program codeby applying adaptive simulated annealing inversion in Matlab environment. We applied this method to determine earthquake hypocenter using several data cases which are regional tectonic, volcano tectonic, and geothermal field. The travel times were calculated using ray tracing shooting method. We then compared its results with the results using Geiger’s method to analyze its reliability. Our results show hypocenter location has smaller RMS error compared to the Geiger’s result that can be statistically associated with better solution. The hypocenter of earthquakes also well correlated with geological structure in the study area. Werecommend using adaptive simulated annealing inversion to relocate hypocenter location in purpose to get precise and accurate earthquake location.
X-38 Application of Dynamic Inversion Flight Control
NASA Technical Reports Server (NTRS)
Wacker, Roger; Munday, Steve; Merkle, Scott
2001-01-01
This paper summarizes the application of a nonlinear dynamic inversion (DI) flight control system (FCS) to an autonomous flight test vehicle in NASA's X-38 Project, a predecessor to the International Space Station (ISS) Crew Return Vehicle (CRV). Honeywell's Multi-Application Control-H (MACH) is a parameterized FCS design architecture including both model-based DI rate-compensation and classical P+I command-tracking. MACH was adopted by X-38 in order to shorten the design cycle time for different vehicle shapes and flight envelopes and evolving aerodynamic databases. Specific design issues and analysis results are presented for the application of MACH to the 3rd free flight (FF3) of X-38 Vehicle 132 (V132). This B-52 drop test, occurring on March 30, 2000, represents the first flight test of MACH and one of the first few known applications of DI in the primary FCS of an autonomous flight test vehicle.
Dynamical Adaptation in Terrorist Cells/Networks
NASA Astrophysics Data System (ADS)
Hussain, D. M. Akbar; Ahmed, Zaki
Typical terrorist cells/networks have dynamical structure as they evolve or adapt to changes which may occur due to capturing or killing of a member of the cell/network. Analytical measures in graph theory like degree centrality, betweenness and closeness centralities are very common and have long history of their successful use in revealing the importance of various members of the network. However, modeling of covert, terrorist or criminal networks through social graph dose not really provide the hierarchical structure which exist in these networks as these networks are composed of leaders and followers etc. In this research we analyze and predict the most likely role a particular node can adapt once a member of the network is either killed or caught. The adaptation is based on computing Bayes posteriori probability of each node and the level of the said node in the network structure.
Dynamic Load Balancing for Adaptive Unstructured Grids
NASA Technical Reports Server (NTRS)
Biswas, Rupak; Saini, Subhash (Technical Monitor)
1998-01-01
Dynamic mesh adaptation on unstructured grids is a powerful tool for computing unsteady three-dimensional problems that require grid modifications to efficiently resolve solution features. By locally refining and coarsening the mesh to capture phenomena of interest, such procedures make standard computational methods more cost effective. Highly refined meshes are required to accurately capture shock waves, contact discontinuities, vortices, and shear layers in fluid flow problems. Adaptive meshes have also proved to be useful in several other areas of computational science and engineering like computer vision and graphics, semiconductor device modeling, and structural mechanics. Local mesh adaptation provides the opportunity to obtain solutions that are comparable to those obtained on globally-refined grids but at a much lower cost. Additional information is contained in the original extended abstract.
Target tracking with dynamically adaptive correlation
NASA Astrophysics Data System (ADS)
Gaxiola, Leopoldo N.; Diaz-Ramirez, Victor H.; Tapia, Juan J.; García-Martínez, Pascuala
2016-04-01
A reliable algorithm for target tracking based on dynamically adaptive correlation filtering is presented. The algorithm is capable of tracking with high accuracy the location of a target in an input video sequence without using an offline training process. The target is selected at the beginning of the algorithm. Afterwards, a composite correlation filter optimized for distortion tolerant pattern recognition is designed to recognize the target in the next frame. The filter is dynamically adapted to each frame using information of current and past scene observations. Results obtained with the proposed algorithm in synthetic and real-life video sequences, are analyzed and compared with those obtained with recent state-of-the-art tracking algorithms in terms of objective metrics.
Adaptive synchronization and anticipatory dynamical systems
NASA Astrophysics Data System (ADS)
Yang, Ying-Jen; Chen, Chun-Chung; Lai, Pik-Yin; Chan, C. K.
2015-09-01
Many biological systems can sense periodical variations in a stimulus input and produce well-timed, anticipatory responses after the input is removed. Such systems show memory effects for retaining timing information in the stimulus and cannot be understood from traditional synchronization consideration of passive oscillatory systems. To understand this anticipatory phenomena, we consider oscillators built from excitable systems with the addition of an adaptive dynamics. With such systems, well-timed post-stimulus responses similar to those from experiments can be obtained. Furthermore, a well-known model of working memory is shown to possess similar anticipatory dynamics when the adaptive mechanism is identified with synaptic facilitation. The last finding suggests that this type of oscillator can be common in neuronal systems with plasticity.
Adaptive Dynamic Event Tree in RAVEN code
Alfonsi, Andrea; Rabiti, Cristian; Mandelli, Diego; Cogliati, Joshua Joseph; Kinoshita, Robert Arthur
2014-11-01
RAVEN is a software tool that is focused on performing statistical analysis of stochastic dynamic systems. RAVEN has been designed in a high modular and pluggable way in order to enable easy integration of different programming languages (i.e., C++, Python) and coupling with other applications (system codes). Among the several capabilities currently present in RAVEN, there are five different sampling strategies: Monte Carlo, Latin Hyper Cube, Grid, Adaptive and Dynamic Event Tree (DET) sampling methodologies. The scope of this paper is to present a new sampling approach, currently under definition and implementation: an evolution of the DET me
Synaptic dynamics: linear model and adaptation algorithm.
Yousefi, Ali; Dibazar, Alireza A; Berger, Theodore W
2014-08-01
In this research, temporal processing in brain neural circuitries is addressed by a dynamic model of synaptic connections in which the synapse model accounts for both pre- and post-synaptic processes determining its temporal dynamics and strength. Neurons, which are excited by the post-synaptic potentials of hundred of the synapses, build the computational engine capable of processing dynamic neural stimuli. Temporal dynamics in neural models with dynamic synapses will be analyzed, and learning algorithms for synaptic adaptation of neural networks with hundreds of synaptic connections are proposed. The paper starts by introducing a linear approximate model for the temporal dynamics of synaptic transmission. The proposed linear model substantially simplifies the analysis and training of spiking neural networks. Furthermore, it is capable of replicating the synaptic response of the non-linear facilitation-depression model with an accuracy better than 92.5%. In the second part of the paper, a supervised spike-in-spike-out learning rule for synaptic adaptation in dynamic synapse neural networks (DSNN) is proposed. The proposed learning rule is a biologically plausible process, and it is capable of simultaneously adjusting both pre- and post-synaptic components of individual synapses. The last section of the paper starts with presenting the rigorous analysis of the learning algorithm in a system identification task with hundreds of synaptic connections which confirms the learning algorithm's accuracy, repeatability and scalability. The DSNN is utilized to predict the spiking activity of cortical neurons and pattern recognition tasks. The DSNN model is demonstrated to be a generative model capable of producing different cortical neuron spiking patterns and CA1 Pyramidal neurons recordings. A single-layer DSNN classifier on a benchmark pattern recognition task outperforms a 2-Layer Neural Network and GMM classifiers while having fewer numbers of free parameters and
NASA Technical Reports Server (NTRS)
Devasia, Santosh; Bayo, Eduardo
1993-01-01
This paper addresses the problem of inverse dynamics for articulated flexible structures with both lumped and distributed actuators. This problem arises, for example, in the combined vibration minimization and trajectory control of space robots and structures. A new inverse dynamics scheme for computing the nominal lumped and distributed inputs for tracking a prescribed trajectory is given.
Riemannian mean and space-time adaptive processing using projection and inversion algorithms
NASA Astrophysics Data System (ADS)
Balaji, Bhashyam; Barbaresco, Frédéric
2013-05-01
The estimation of the covariance matrix from real data is required in the application of space-time adaptive processing (STAP) to an airborne ground moving target indication (GMTI) radar. A natural approach to estimation of the covariance matrix that is based on the information geometry has been proposed. In this paper, the output of the Riemannian mean is used in inversion and projection algorithms. It is found that the projection class of algorithms can yield very significant gains, even when the gains due to inversion-based algorithms are marginal over standard algorithms. The performance of the projection class of algorithms does not appear to be overly sensitive to the projected subspace dimension.
An adaptive importance sampling algorithm for Bayesian inversion with multimodal distributions
Li, Weixuan; Lin, Guang
2015-08-01
Parametric uncertainties are encountered in the simulations of many physical systems, and may be reduced by an inverse modeling procedure that calibrates the simulation results to observations on the real system being simulated. Following Bayes' rule, a general approach for inverse modeling problems is to sample from the posterior distribution of the uncertain model parameters given the observations. However, the large number of repetitive forward simulations required in the sampling process could pose a prohibitive computational burden. This difficulty is particularly challenging when the posterior is multimodal. We present in this paper an adaptive importance sampling algorithm to tackle these challenges. Two essential ingredients of the algorithm are: 1) a Gaussian mixture (GM) model adaptively constructed as the proposal distribution to approximate the possibly multimodal target posterior, and 2) a mixture of polynomial chaos (PC) expansions, built according to the GM proposal, as a surrogate model to alleviate the computational burden caused by computational-demanding forward model evaluations. In three illustrative examples, the proposed adaptive importance sampling algorithm demonstrates its capabilities of automatically finding a GM proposal with an appropriate number of modes for the specific problem under study, and obtaining a sample accurately and efficiently representing the posterior with limited number of forward simulations.
An adaptive importance sampling algorithm for Bayesian inversion with multimodal distributions
Li, Weixuan; Lin, Guang
2015-03-21
Parametric uncertainties are encountered in the simulations of many physical systems, and may be reduced by an inverse modeling procedure that calibrates the simulation results to observations on the real system being simulated. Following Bayes’ rule, a general approach for inverse modeling problems is to sample from the posterior distribution of the uncertain model parameters given the observations. However, the large number of repetitive forward simulations required in the sampling process could pose a prohibitive computational burden. This difficulty is particularly challenging when the posterior is multimodal. We present in this paper an adaptive importance sampling algorithm to tackle these challenges. Two essential ingredients of the algorithm are: 1) a Gaussian mixture (GM) model adaptively constructed as the proposal distribution to approximate the possibly multimodal target posterior, and 2) a mixture of polynomial chaos (PC) expansions, built according to the GM proposal, as a surrogate model to alleviate the computational burden caused by computational-demanding forward model evaluations. In three illustrative examples, the proposed adaptive importance sampling algorithm demonstrates its capabilities of automatically finding a GM proposal with an appropriate number of modes for the specific problem under study, and obtaining a sample accurately and efficiently representing the posterior with limited number of forward simulations.
An adaptive importance sampling algorithm for Bayesian inversion with multimodal distributions
Li, Weixuan; Lin, Guang
2015-03-21
Parametric uncertainties are encountered in the simulations of many physical systems, and may be reduced by an inverse modeling procedure that calibrates the simulation results to observations on the real system being simulated. Following Bayes’ rule, a general approach for inverse modeling problems is to sample from the posterior distribution of the uncertain model parameters given the observations. However, the large number of repetitive forward simulations required in the sampling process could pose a prohibitive computational burden. This difficulty is particularly challenging when the posterior is multimodal. We present in this paper an adaptive importance sampling algorithm to tackle thesemore » challenges. Two essential ingredients of the algorithm are: 1) a Gaussian mixture (GM) model adaptively constructed as the proposal distribution to approximate the possibly multimodal target posterior, and 2) a mixture of polynomial chaos (PC) expansions, built according to the GM proposal, as a surrogate model to alleviate the computational burden caused by computational-demanding forward model evaluations. In three illustrative examples, the proposed adaptive importance sampling algorithm demonstrates its capabilities of automatically finding a GM proposal with an appropriate number of modes for the specific problem under study, and obtaining a sample accurately and efficiently representing the posterior with limited number of forward simulations.« less
Reentry Vehicle Flight Controls Design Guidelines: Dynamic Inversion
NASA Technical Reports Server (NTRS)
Ito, Daigoro; Georgie, Jennifer; Valasek, John; Ward, Donald T.
2002-01-01
This report addresses issues in developing a flight control design for vehicles operating across a broad flight regime and with highly nonlinear physical descriptions of motion. Specifically it addresses the need for reentry vehicles that could operate through reentry from space to controlled touchdown on Earth. The latter part of controlled descent is achieved by parachute or paraglider - or by all automatic or a human-controlled landing similar to that of the Orbiter. Since this report addresses the specific needs of human-carrying (not necessarily piloted) reentry vehicles, it deals with highly nonlinear equations of motion, and then-generated control systems must be robust across a very wide range of physics. Thus, this report deals almost exclusively with some form of dynamic inversion (DI). Two vital aspects of control theory - noninteracting control laws and the transformation of nonlinear systems into equivalent linear systems - are embodied in DI. Though there is no doubt that the mathematical tools and underlying theory are widely available, there are open issues as to the practicality of using DI as the only or primary design approach for reentry articles. This report provides a set of guidelines that can be used to determine the practical usefulness of the technique.
Molecular dynamics simulations of ring inversion in RDX
NASA Astrophysics Data System (ADS)
Wallis, Eric P.; Thompson, Donald L.
1993-08-01
Molecular dynamics simulations, using the finite volume method of Murrell and co-workers [J. Chem. Phys. 94, 3908 (1991)], have been carried out to study conformational changes in hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in isolation and in dense Xe gas. The configurational distributions for RDX in a Xe bath and in the gas-phase are markedly different. The results show that as the solvent concentration increases, the concentration of RDX molecules in the boat conformation increases by a factor of about 4. The rate constant for the chair→boat ring inversion was calculated as a function of the xenon concentration [Xe]. The rate constant obeys Lindemann behavior at low concentrations, i.e., it increases with increasing solvent density. At [Xe]˜6.2 mol dm-3, the rate constant reaches a maximum (Kramer's turnover) and becomes a decreasing function of the solvent concentration. For [Xe] above 16.2 mol dm-3, the rate constant again increases as a function of the solvent density.
A covariance-adaptive approach for regularized inversion in linear models
NASA Astrophysics Data System (ADS)
Kotsakis, Christopher
2007-11-01
The optimal inversion of a linear model under the presence of additive random noise in the input data is a typical problem in many geodetic and geophysical applications. Various methods have been developed and applied for the solution of this problem, ranging from the classic principle of least-squares (LS) estimation to other more complex inversion techniques such as the Tikhonov-Philips regularization, truncated singular value decomposition, generalized ridge regression, numerical iterative methods (Landweber, conjugate gradient) and others. In this paper, a new type of optimal parameter estimator for the inversion of a linear model is presented. The proposed methodology is based on a linear transformation of the classic LS estimator and it satisfies two basic criteria. First, it provides a solution for the model parameters that is optimally fitted (in an average quadratic sense) to the classic LS parameter solution. Second, it complies with an external user-dependent constraint that specifies a priori the error covariance (CV) matrix of the estimated model parameters. The formulation of this constrained estimator offers a unified framework for the description of many regularization techniques that are systematically used in geodetic inverse problems, particularly for those methods that correspond to an eigenvalue filtering of the ill-conditioned normal matrix in the underlying linear model. Our study lies on the fact that it adds an alternative perspective on the statistical properties and the regularization mechanism of many inversion techniques commonly used in geodesy and geophysics, by interpreting them as a family of `CV-adaptive' parameter estimators that obey a common optimal criterion and differ only on the pre-selected form of their error CV matrix under a fixed model design.
Use of reduced basis technique in the inverse dynamics of large space cranes
NASA Technical Reports Server (NTRS)
Das, S. K.; Utku, S.; Wada, B. K.
1990-01-01
The inverse dynamics of adaptive structures used as space cranes can prove computationally expensive in the case of large structures, due to the large number of degrees of freedom involved. Consequently, reduced basis techniques (reduction techniques) are frequently used to reduce the problem size to a time manageable level (for possible use in real time control). A reduced basis technique is proposed which is different from, but related to, the path-derivatives reduction technique. A linearly independent set of deflection n-tuples is used, chosen at the beginning of the time range in which it is wished to reduce the equations, in whose subspace it is assumed that the deflection vectors of the unreduced problem will lie (approximately).
Emerging hierarchies in dynamically adapting webs
NASA Astrophysics Data System (ADS)
Katifori, Eleni; Graewer, Johannes; Magnasco, Marcelo; Modes, Carl
Transport networks play a key role across four realms of eukaryotic life: slime molds, fungi, plants, and animals. In addition to the developmental algorithms that build them, many also employ adaptive strategies to respond to stimuli, damage, and other environmental changes. We model these adapting network architectures using a generic dynamical system on weighted graphs and find in simulation that these networks ultimately develop a hierarchical organization of the final weighted architecture accompanied by the formation of a system-spanning backbone. We quantify the hierarchical organization of the networks by developing an algorithm that decomposes the architecture to multiple scales and analyzes how the organization in each scale relates to that of the scale above and below it. The methodologies developed in this work are applicable to a wide range of systems including the slime mold physarum polycephalum, human microvasculature, and force chains in granular media.
Cardiac fluid dynamics anticipates heart adaptation.
Pedrizzetti, Gianni; Martiniello, Alfonso R; Bianchi, Valter; D'Onofrio, Antonio; Caso, Pio; Tonti, Giovanni
2015-01-21
Hemodynamic forces represent an epigenetic factor during heart development and are supposed to influence the pathology of the grown heart. Cardiac blood motion is characterized by a vortical dynamics, and it is common belief that the cardiac vortex has a role in disease progressions or regression. Here we provide a preliminary demonstration about the relevance of maladaptive intra-cardiac vortex dynamics in the geometrical adaptation of the dysfunctional heart. We employed an in vivo model of patients who present a stable normal heart function in virtue of the cardiac resynchronization therapy (CRT, bi-ventricular pace-maker) and who are expected to develop left ventricle remodeling if pace-maker was switched off. Intra-ventricular fluid dynamics is analyzed by echocardiography (Echo-PIV). Under normal conditions, the flow presents a longitudinal alignment of the intraventricular hemodynamic forces. When pacing is temporarily switched off, flow forces develop a misalignment hammering onto lateral walls, despite no other electro-mechanical change is noticed. Hemodynamic forces result to be the first event that evokes a physiological activity anticipating cardiac changes and could help in the prediction of longer term heart adaptations. PMID:25529139
Mental workload dynamics in adaptive interface design
NASA Technical Reports Server (NTRS)
Hancock, Peter A.; Chignell, Mark H.
1988-01-01
In examining the role of time in mental workload, the authors present a different perspective from which to view the problem of assessment. Mental workload is plotted in three dimensions, whose axes represent effective time for action, perceived distance from desired goal state, level of effort required to achieve the time-constrained goal. This representation allows the generation of isodynamic workload contours that incorporate the factors of operator skill and equifinality of effort. An adaptive interface for dynamic task reallocation is described that uses this form of assessment to reconcile the joint aims of stable operator loading and acceptable primary task performance by the total system.
An Adaptive ANOVA-based PCKF for High-Dimensional Nonlinear Inverse Modeling
LI, Weixuan; Lin, Guang; Zhang, Dongxiao
2014-02-01
The probabilistic collocation-based Kalman filter (PCKF) is a recently developed approach for solving inverse problems. It resembles the ensemble Kalman filter (EnKF) in every aspect—except that it represents and propagates model uncertainty by polynomial chaos expansion (PCE) instead of an ensemble of model realizations. Previous studies have shown PCKF is a more efficient alternative to EnKF for many data assimilation problems. However, the accuracy and efficiency of PCKF depends on an appropriate truncation of the PCE series. Having more polynomial chaos bases in the expansion helps to capture uncertainty more accurately but increases computational cost. Bases selection is particularly important for high-dimensional stochastic problems because the number of polynomial chaos bases required to represent model uncertainty grows dramatically as the number of input parameters (random dimensions) increases. In classic PCKF algorithms, the PCE bases are pre-set based on users’ experience. Also, for sequential data assimilation problems, the bases kept in PCE expression remain unchanged in different Kalman filter loops, which could limit the accuracy and computational efficiency of classic PCKF algorithms. To address this issue, we present a new algorithm that adaptively selects PCE bases for different problems and automatically adjusts the number of bases in different Kalman filter loops. The algorithm is based on adaptive functional ANOVA (analysis of variance) decomposition, which approximates a high-dimensional function with the summation of a set of low-dimensional functions. Thus, instead of expanding the original model into PCE, we implement the PCE expansion on these low-dimensional functions, which is much less costly. We also propose a new adaptive criterion for ANOVA that is more suited for solving inverse problems. The new algorithm is tested with different examples and demonstrated great effectiveness in comparison with non-adaptive PCKF and En
An adaptive ANOVA-based PCKF for high-dimensional nonlinear inverse modeling
Li, Weixuan; Lin, Guang; Zhang, Dongxiao
2014-02-01
The probabilistic collocation-based Kalman filter (PCKF) is a recently developed approach for solving inverse problems. It resembles the ensemble Kalman filter (EnKF) in every aspect—except that it represents and propagates model uncertainty by polynomial chaos expansion (PCE) instead of an ensemble of model realizations. Previous studies have shown PCKF is a more efficient alternative to EnKF for many data assimilation problems. However, the accuracy and efficiency of PCKF depends on an appropriate truncation of the PCE series. Having more polynomial chaos basis functions in the expansion helps to capture uncertainty more accurately but increases computational cost. Selection of basis functions is particularly important for high-dimensional stochastic problems because the number of polynomial chaos basis functions required to represent model uncertainty grows dramatically as the number of input parameters (random dimensions) increases. In classic PCKF algorithms, the PCE basis functions are pre-set based on users' experience. Also, for sequential data assimilation problems, the basis functions kept in PCE expression remain unchanged in different Kalman filter loops, which could limit the accuracy and computational efficiency of classic PCKF algorithms. To address this issue, we present a new algorithm that adaptively selects PCE basis functions for different problems and automatically adjusts the number of basis functions in different Kalman filter loops. The algorithm is based on adaptive functional ANOVA (analysis of variance) decomposition, which approximates a high-dimensional function with the summation of a set of low-dimensional functions. Thus, instead of expanding the original model into PCE, we implement the PCE expansion on these low-dimensional functions, which is much less costly. We also propose a new adaptive criterion for ANOVA that is more suited for solving inverse problems. The new algorithm was tested with different examples and demonstrated
Direct adaptive control for nonlinear uncertain dynamical systems
NASA Astrophysics Data System (ADS)
Hayakawa, Tomohisa
In light of the complex and highly uncertain nature of dynamical systems requiring controls, it is not surprising that reliable system models for many high performance engineering and life science applications are unavailable. In the face of such high levels of system uncertainty, robust controllers may unnecessarily sacrifice system performance whereas adaptive controllers are clearly appropriate since they can tolerate far greater system uncertainty levels to improve system performance. In this dissertation, we develop a Lyapunov-based direct adaptive and neural adaptive control framework that addresses parametric uncertainty, unstructured uncertainty, disturbance rejection, amplitude and rate saturation constraints, and digital implementation issues. Specifically, we consider the following research topics; direct adaptive control for nonlinear uncertain systems with exogenous disturbances; robust adaptive control for nonlinear uncertain systems; adaptive control for nonlinear uncertain systems with actuator amplitude and rate saturation constraints; adaptive reduced-order dynamic compensation for nonlinear uncertain systems; direct adaptive control for nonlinear matrix second-order dynamical systems with state-dependent uncertainty; adaptive control for nonnegative and compartmental dynamical systems with applications to general anesthesia; direct adaptive control of nonnegative and compartmental dynamical systems with time delay; adaptive control for nonlinear nonnegative and compartmental dynamical systems with applications to clinical pharmacology; neural network adaptive control for nonlinear nonnegative dynamical systems; passivity-based neural network adaptive output feedback control for nonlinear nonnegative dynamical systems; neural network adaptive dynamic output feedback control for nonlinear nonnegative systems using tapped delay memory units; Lyapunov-based adaptive control framework for discrete-time nonlinear systems with exogenous disturbances
NASA Astrophysics Data System (ADS)
Ma, Xiang; Zabaras, Nicholas
2009-03-01
A new approach to modeling inverse problems using a Bayesian inference method is introduced. The Bayesian approach considers the unknown parameters as random variables and seeks the probabilistic distribution of the unknowns. By introducing the concept of the stochastic prior state space to the Bayesian formulation, we reformulate the deterministic forward problem as a stochastic one. The adaptive hierarchical sparse grid collocation (ASGC) method is used for constructing an interpolant to the solution of the forward model in this prior space which is large enough to capture all the variability/uncertainty in the posterior distribution of the unknown parameters. This solution can be considered as a function of the random unknowns and serves as a stochastic surrogate model for the likelihood calculation. Hierarchical Bayesian formulation is used to derive the posterior probability density function (PPDF). The spatial model is represented as a convolution of a smooth kernel and a Markov random field. The state space of the PPDF is explored using Markov chain Monte Carlo algorithms to obtain statistics of the unknowns. The likelihood calculation is performed by directly sampling the approximate stochastic solution obtained through the ASGC method. The technique is assessed on two nonlinear inverse problems: source inversion and permeability estimation in flow through porous media.
ADAPTIVE MULTILEVEL SPLITTING IN MOLECULAR DYNAMICS SIMULATIONS*
Aristoff, David; Lelièvre, Tony; Mayne, Christopher G.; Teo, Ivan
2014-01-01
Adaptive Multilevel Splitting (AMS) is a replica-based rare event sampling method that has been used successfully in high-dimensional stochastic simulations to identify trajectories across a high potential barrier separating one metastable state from another, and to estimate the probability of observing such a trajectory. An attractive feature of AMS is that, in the limit of a large number of replicas, it remains valid regardless of the choice of reaction coordinate used to characterize the trajectories. Previous studies have shown AMS to be accurate in Monte Carlo simulations. In this study, we extend the application of AMS to molecular dynamics simulations and demonstrate its effectiveness using a simple test system. Our conclusion paves the way for useful applications, such as molecular dynamics calculations of the characteristic time of drug dissociation from a protein target. PMID:26005670
Circuit dynamics of adaptive and maladaptive behaviour
Deisseroth, Karl
2014-01-01
The recent development of technologies for investigating specific components of intact biological systems has allowed elucidation of the neural circuitry underlying adaptive and maladaptive behaviours. Investigators are now able to observe and control, with high spatio-temporal resolution, structurally defined intact pathways along which electrical activity flows during and after the performance of complex behaviours. These investigations have revealed that control of projection-specific dynamics is well suited to modulating behavioural patterns that are relevant to a broad range of psychiatric diseases. Structural dynamics principles have emerged to provide diverse, unexpected and causal insights into the operation of intact and diseased nervous systems, linking form and function in the brain. PMID:24429629
Approximated Stable Inversion for Nonlinear Systems with Nonhyperbolic Internal Dynamics. Revised
NASA Technical Reports Server (NTRS)
Devasia, Santosh
1999-01-01
A technique to achieve output tracking for nonminimum phase nonlinear systems with non- hyperbolic internal dynamics is presented. The present paper integrates stable inversion techniques (that achieve exact-tracking) with approximation techniques (that modify the internal dynamics) to circumvent the nonhyperbolicity of the internal dynamics - this nonhyperbolicity is an obstruction to applying presently available stable inversion techniques. The theory is developed for nonlinear systems and the method is applied to a two-cart with inverted-pendulum example.
Force and Moment Approach for Achievable Dynamics Using Nonlinear Dynamic Inversion
NASA Technical Reports Server (NTRS)
Ostroff, Aaron J.; Bacon, Barton J.
1999-01-01
This paper describes a general form of nonlinear dynamic inversion control for use in a generic nonlinear simulation to evaluate candidate augmented aircraft dynamics. The implementation is specifically tailored to the task of quickly assessing an aircraft's control power requirements and defining the achievable dynamic set. The achievable set is evaluated while undergoing complex mission maneuvers, and perfect tracking will be accomplished when the desired dynamics are achievable. Variables are extracted directly from the simulation model each iteration, so robustness is not an issue. Included in this paper is a description of the implementation of the forces and moments from simulation variables, the calculation of control effectiveness coefficients, methods for implementing different types of aerodynamic and thrust vectoring controls, adjustments for control effector failures, and the allocation approach used. A few examples illustrate the perfect tracking results obtained.
Adaptive dynamics for physiologically structured population models.
Durinx, Michel; Metz, J A J Hans; Meszéna, Géza
2008-05-01
We develop a systematic toolbox for analyzing the adaptive dynamics of multidimensional traits in physiologically structured population models with point equilibria (sensu Dieckmann et al. in Theor. Popul. Biol. 63:309-338, 2003). Firstly, we show how the canonical equation of adaptive dynamics (Dieckmann and Law in J. Math. Biol. 34:579-612, 1996), an approximation for the rate of evolutionary change in characters under directional selection, can be extended so as to apply to general physiologically structured population models with multiple birth states. Secondly, we show that the invasion fitness function (up to and including second order terms, in the distances of the trait vectors to the singularity) for a community of N coexisting types near an evolutionarily singular point has a rational form, which is model-independent in the following sense: the form depends on the strategies of the residents and the invader, and on the second order partial derivatives of the one-resident fitness function at the singular point. This normal form holds for Lotka-Volterra models as well as for physiologically structured population models with multiple birth states, in discrete as well as continuous time and can thus be considered universal for the evolutionary dynamics in the neighbourhood of singular points. Only in the case of one-dimensional trait spaces or when N = 1 can the normal form be reduced to a Taylor polynomial. Lastly we show, in the form of a stylized recipe, how these results can be combined into a systematic approach for the analysis of the (large) class of evolutionary models that satisfy the above restrictions. PMID:17943289
NASA Technical Reports Server (NTRS)
Bayo, Eduardo; Ledesma, Ragnar
1993-01-01
A technique is presented for solving the inverse dynamics of flexible planar multibody systems. This technique yields the non-causal joint efforts (inverse dynamics) as well as the internal states (inverse kinematics) that produce a prescribed nominal trajectory of the end effector. A non-recursive global Lagrangian approach is used in formulating the equations for motion as well as in solving the inverse dynamics equations. Contrary to the recursive method previously presented, the proposed method solves the inverse problem in a systematic and direct manner for both open-chain as well as closed-chain configurations. Numerical simulation shows that the proposed procedure provides an excellent tracking of the desired end effector trajectory.
Opinion dynamics on an adaptive random network
NASA Astrophysics Data System (ADS)
Benczik, I. J.; Benczik, S. Z.; Schmittmann, B.; Zia, R. K. P.
2009-04-01
We revisit the classical model for voter dynamics in a two-party system with two basic modifications. In contrast to the original voter model studied in regular lattices, we implement the opinion formation process in a random network of agents in which interactions are no longer restricted by geographical distance. In addition, we incorporate the rapidly changing nature of the interpersonal relations in the model. At each time step, agents can update their relationships. This update is determined by their own opinion, and by their preference to make connections with individuals sharing the same opinion, or rather with opponents. In this way, the network is built in an adaptive manner, in the sense that its structure is correlated and evolves with the dynamics of the agents. The simplicity of the model allows us to examine several issues analytically. We establish criteria to determine whether consensus or polarization will be the outcome of the dynamics and on what time scales these states will be reached. In finite systems consensus is typical, while in infinite systems a disordered metastable state can emerge and persist for infinitely long time before consensus is reached.
Modified Dynamic Inversion to Control Large Flexible Aircraft: What's Going On?
NASA Technical Reports Server (NTRS)
Gregory, Irene M.
1999-01-01
High performance aircraft of the future will be designed lighter, more maneuverable, and operate over an ever expanding flight envelope. One of the largest differences from the flight control perspective between current and future advanced aircraft is elasticity. Over the last decade, dynamic inversion methodology has gained considerable popularity in application to highly maneuverable fighter aircraft, which were treated as rigid vehicles. This paper explores dynamic inversion application to an advanced highly flexible aircraft. An initial application has been made to a large flexible supersonic aircraft. In the course of controller design for this advanced vehicle, modifications were made to the standard dynamic inversion methodology. The results of this application were deemed rather promising. An analytical study has been undertaken to better understand the nature of the made modifications and to determine its general applicability. This paper presents the results of this initial analytical look at the modifications to dynamic inversion to control large flexible aircraft.
Nitric oxide regulates vascular adaptive mitochondrial dynamics.
Miller, Matthew W; Knaub, Leslie A; Olivera-Fragoso, Luis F; Keller, Amy C; Balasubramaniam, Vivek; Watson, Peter A; Reusch, Jane E B
2013-06-15
Cardiovascular disease risk factors, such as diabetes, hypertension, dyslipidemia, obesity, and physical inactivity, are all correlated with impaired endothelial nitric oxide synthase (eNOS) function and decreased nitric oxide (NO) production. NO-mediated regulation of mitochondrial biogenesis has been established in many tissues, yet the role of eNOS in vascular mitochondrial biogenesis and dynamics is unclear. We hypothesized that genetic eNOS deletion and 3-day nitric oxide synthase (NOS) inhibition in rodents would result in impaired mitochondrial biogenesis and defunct fission/fusion and autophagy profiles within the aorta. We observed a significant, eNOS expression-dependent decrease in mitochondrial electron transport chain (ETC) protein subunits from complexes I, II, III, and V in eNOS heterozygotes and eNOS null mice compared with age-matched controls. In response to NOS inhibition with NG-nitro-L-arginine methyl ester (L-NAME) treatment in Sprague Dawley rats, significant decreases were observed in ETC protein subunits from complexes I, III, and IV as well as voltage-dependent anion channel 1. Decreased protein content of upstream regulators of mitochondrial biogenesis, cAMP response element-binding protein and peroxisome proliferator-activated receptor-γ coactivator-1α, were observed in response to 3-day L-NAME treatment. Both genetic eNOS deletion and NOS inhibition resulted in decreased manganese superoxide dismutase protein. L-NAME treatment resulted in significant changes to mitochondrial dynamic protein profiles with decreased fusion, increased fission, and minimally perturbed autophagy. In addition, L-NAME treatment blocked mitochondrial adaptation to an exercise intervention in the aorta. These results suggest that eNOS/NO play a role in basal and adaptive mitochondrial biogenesis in the vasculature and regulation of mitochondrial turnover. PMID:23585138
Mission to Mars: Adaptive Identifier for the Solution of Inverse Optical Metrology Tasks
NASA Astrophysics Data System (ADS)
Krapivin, Vladimir F.; Varotsos, Costas A.; Christodoulakis, John
2016-06-01
A human mission to Mars requires the solution of many problems that mainly linked to the safety of life, the reliable operational control of drinking water as well as health care. The availability of liquid fuels is also an important issue since the existing tools cannot fully provide the required liquid fuels quantities for the mission return journey. This paper presents the development of new methods and technology for reliable, operational, and with high availability chemical analysis of liquid solutions of various types. This technology is based on the employment of optical sensors (such as the multi-channel spectrophotometers or spectroellipsometers and microwave radiometers) and the development of a database of spectral images for typical liquid solutions that could be the objects of life on Mars. This database exploits the adaptive recognition of optical images of liquids using specific algorithms that are based on spectral analysis, cluster analysis and methods for solving the inverse optical metrology tasks.
Mission to Mars: Adaptive Identifier for the Solution of Inverse Optical Metrology Tasks
NASA Astrophysics Data System (ADS)
Krapivin, Vladimir F.; Varotsos, Costas A.; Christodoulakis, John
2016-04-01
A human mission to Mars requires the solution of many problems that mainly linked to the safety of life, the reliable operational control of drinking water as well as health care. The availability of liquid fuels is also an important issue since the existing tools cannot fully provide the required liquid fuels quantities for the mission return journey. This paper presents the development of new methods and technology for reliable, operational, and with high availability chemical analysis of liquid solutions of various types. This technology is based on the employment of optical sensors (such as the multi-channel spectrophotometers or spectroellipsometers and microwave radiometers) and the development of a database of spectral images for typical liquid solutions that could be the objects of life on Mars. This database exploits the adaptive recognition of optical images of liquids using specific algorithms that are based on spectral analysis, cluster analysis and methods for solving the inverse optical metrology tasks.
High resolution imaging of the Earth with adaptive full-waveform inversion
NASA Astrophysics Data System (ADS)
Morgan, J. V.; Warner, M.; Guasch, L.; Umpleby, A.; Yao, G.; Herrmann, F. J.
2014-12-01
Three-dimensional full-waveform inversion (FWI) is a high-resolution, high-fidelity, quantitative, seismic imaging technique that has advanced rapidly within the oil and gas industry. The method involves the iterative improvement of a starting model using a series of local linearized updates to solve the full non-linear inversion problem. During the inversion, forward modeling employs the full two-way three-dimensional heterogeneous anisotropic acoustic or elastic wave equation to predict the observed raw field data, wiggle-for-wiggle, trace-by-trace. The method is computationally demanding; it is highly parallelized, and runs on large multi-core multi-node clusters. A recently developed adaptive version of FWI is able to overcome the requirement for a good starting model and low frequencies in the data, and this opens up the range of datasets and problems to which FWI can be applied. Here, we demonstrate what can be achieved by applying this newly practical technique to high-density 3D seismic datasets acquired to image petroleum targets. We show that the resulting anisotropic p-wave velocity models match in situ measurements in boreholes, reproduce detailed structure observed independently on high-resolution seismic reflection sections, accurately predict the raw seismic data, and simplify and sharpen reverse-time-migrated reflection images of deeper horizons. The velocity models image individual faults, gas clouds, channels, and other geological features with previously unobtainable resolution and clarity. These same benefits can be obtained when this technique is applied to scientific targets provided that the data coverage is adequate in three-dimensions, and that an appropriate range of offsets and azimuths are available. Possible targets range from the water column, ice sheets, and Holocene deposits, through active faults, spreading centres, collision zones, rifted margins, magma plumbing, lower-continental crust, and deep crustal hot zones, to whole
Nie Xiaobo; Liang Jian; Yan Di
2012-12-15
Purpose: To create an organ sample generator (OSG) for expected treatment dose construction and adaptive inverse planning optimization. The OSG generates random samples of organs of interest from a distribution obeying the patient specific organ variation probability density function (PDF) during the course of adaptive radiotherapy. Methods: Principle component analysis (PCA) and a time-varying least-squares regression (LSR) method were used on patient specific geometric variations of organs of interest manifested on multiple daily volumetric images obtained during the treatment course. The construction of the OSG includes the determination of eigenvectors of the organ variation using PCA, and the determination of the corresponding coefficients using time-varying LSR. The coefficients can be either random variables or random functions of the elapsed treatment days depending on the characteristics of organ variation as a stationary or a nonstationary random process. The LSR method with time-varying weighting parameters was applied to the precollected daily volumetric images to determine the function form of the coefficients. Eleven h and n cancer patients with 30 daily cone beam CT images each were included in the evaluation of the OSG. The evaluation was performed using a total of 18 organs of interest, including 15 organs at risk and 3 targets. Results: Geometric variations of organs of interest during h and n cancer radiotherapy can be represented using the first 3 {approx} 4 eigenvectors. These eigenvectors were variable during treatment, and need to be updated using new daily images obtained during the treatment course. The OSG generates random samples of organs of interest from the estimated organ variation PDF of the individual. The accuracy of the estimated PDF can be improved recursively using extra daily image feedback during the treatment course. The average deviations in the estimation of the mean and standard deviation of the organ variation PDF for h
NASA Astrophysics Data System (ADS)
Ryerson, F. J.; Ezzedine, S. M.; Antoun, T.
2013-12-01
equation for the distribution of k is solved, provided that Cauchy data are appropriately assigned. In the next stage, only a limited number of passive measurements are provided. In this case, the forward and inverse PDEs are solved simultaneously. This is accomplished by adding regularization terms and filtering the pressure gradients in the inverse problem. Both the forward and the inverse problem are either simultaneously or sequentially coupled and solved using implicit schemes, adaptive mesh refinement, Galerkin finite elements. The final case arises when P, k, and Q data only exist at producing wells. This exceedingly ill posed problem calls for additional constraints on the forward-inverse coupling to insure that the production rates are satisfied at the desired locations. Results from all three cases are presented demonstrating stability and accuracy of the proposed approach and, more importantly, providing some insights into the consequences of data under sampling, uncertainty propagation and quantification. We illustrate the advantages of this novel approach over the common UQ forward drivers on several subsurface energy problems in either porous or fractured or/and faulted reservoirs. This work performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.
Inversion of dynamically repositioned multi-axis electromagnetic data for ordnance characterization
NASA Astrophysics Data System (ADS)
Keranen, Joe; Shubitidze, Fridon; Besaw, Lance; Casari, Matthew J.; Miller, Jonathan; Schultz, Gregory
2011-06-01
The challenges associated with removing UXO and explosive remnants of war have led to a variety of methods for detection and discrimination of buried metallic objects using time-domain electromagnetic induction (EMI). Recent work has shown that parameters recovered from physics-based inversions can discriminate and classify buried ordnance from non-ordnance. We present results of applying advanced processing to data from a dynamically repositioned multiaxis EMI instrument. Data are collected using an adaptive sampling process to find the center of the anomaly and collect minimal data while maintaining model fidelity. An ortho-normalized volume magnetic source (ONVMS) model is used to resolve various targets at different depths. The ONVMS model is a generalized volume dipole model, with the single dipole model being a special limiting case. Using the ONVMS model, an object's response to a sensor's primary magnetic field is modeled mathematically by a set of equivalent magnetic dipoles distributed inside a volume containing the object. We assess the utility and veracity of the dynamic sampling strategy coupled with the ONVMS model on data acquired over a set of calibration and simulant targets. Rapid target characterization codes are aggregated into a software package with particular focus on ease of use for non-expert users.
Introducing interactive inverse FEM simulation and its application for adaptive radiotherapy.
Coevoet, Eulalie; Reynaert, Nick; Lartigau, Eric; Schiappacasse, Luis; Dequidt, Jérémie; Duriez, Christian
2014-01-01
We introduce a new methodology for semi-automatic deformable registration of anatomical structures, using interactive inverse simulations. The method relies on non-linear real-time Finite Element Method (FEM) within a constraint-based framework. Given a set of few registered points provided by the user, a real-time optimization adapts the boundary conditions and(/or) some parameters of the FEM in order to obtain the adequate geometrical deformations. To dramatically fasten the process, the method relies on a projection of the model in the space of the optimization variables. In this reduced space, a quadratic programming problem is formulated and solved very quickly. The method is validated with numerical examples for retrieving Young's modulus and some pressures on the boundaries. Then, we apply the approach for the registration of the parotid glands during the radiotherapy of the head and neck cancer. Radiotherapy treatment induces weight loss that modifies the shape and the positions of these structures and they eventually intersect the target volume. We show how we could adapt the planning to limit the radiation of these glands. PMID:25485365
Efficient solution of an inverse problem in cell population dynamics
NASA Astrophysics Data System (ADS)
Groh, Andreas; Krebs, Jochen; Wagner, Mathias
2011-06-01
In this paper, a size-structured model for cell division is examined and the question of determining the division (birth) rate from a measurable stable size distribution of the population is addressed. This inverse problem can be formulated as a differential-dilation equation. We propose a novel solution scheme based on mollification. The method of approximate inverse allows us to shift the derivative from the data to a precomputable reconstruction kernel. To comprise all available a priori information, a presmoothing step based on regression in reproducing kernel Hilbert spaces is introduced. We establish an error theory for the emerging algorithm, prove convergence and deduce a parameter strategy. The results are substantiated with extensive numerical tests both for artificial and real data based on proliferating tumor cells.
Adaptive fuzzy control with smooth inverse for nonlinear systems preceded by non-symmetric dead-zone
NASA Astrophysics Data System (ADS)
Wang, Xingjian; Wang, Shaoping
2016-07-01
In this study, the adaptive output feedback control problem of a class of nonlinear systems preceded by non-symmetric dead-zone is considered. To cope with the possible control signal chattering phenomenon which is caused by non-smooth dead-zone inverse, a new smooth inverse is proposed for non-symmetric dead-zone compensation. For the systematic design procedure of the adaptive fuzzy control algorithm, we combine the backstepping technique and small-gain approach. The Takagi-Sugeno fuzzy logic systems are used to approximate unknown system nonlinearities. The closed-loop stability is studied by using small gain theorem and the closed-loop system is proved to be semi-globally uniformly ultimately bounded. Simulation results indicate that, compared to the algorithm with the non-smooth inverse, the proposed control strategy can achieve better tracking performance and the chattering phenomenon can be avoided effectively.
Investigating the reliability of kinematic source inversion with dynamic rupture models
NASA Astrophysics Data System (ADS)
Zhang, Y.; Song, S.; Dalguer, L. A.; Clinton, J. F.
2011-12-01
An essential element of understanding the earthquake source processes is obtaining a reliable source model via geophysical data inversion. However, the epistemic uncertainties in the kinematic source inversion produce a variety of source model estimates for any given event. Thus, as done in the Source Inversion Validation (SIV) project, it is important to validate our inversion methods with synthetic data by testing forward Green's function calculation and comparing various inversion methods. Spontaneous dynamic rupture modeling, which incorporates the conservation laws of continuum mechanics and the constitutive behavior of rocks under frictional sliding, is capable of producing physically self-consistent kinematic description of the fault and its associated seismic wave propagation resulting in ground motions on the surface. Here we develop accurate dynamic rupture simulation of a vertical strike slip fault. Our source model is composed of well-defined asperities (patches of large stress drop) and we assume that fault rupture is governed by the linear slip weakening friction model. The resulting near-source ground motions dominated by low frequency (up to 1Hz) are used for testing our inversion method. We performed various inversion tests and compared estimated solutions with true solutions obtained by the forward dynamic rupture modeling. Our preliminary results show that estimated model spaces could be significantly perturbed, depending on data and modeling schemes used in the inversion, not only in terms of spatial distribution of model parameters, but also in terms of their auto- and cross-correlation structure. The Bayesian approach in source inversion is becoming increasingly popular because of the recent common availability of high performance computing capabilities. We adopted the Bayesian approach in our source inversion test, so that we can more effectively analyze the uncertainty of estimated models and also implement physically guided regularization
A fast inverse dynamics model of walking for use in optimisation studies.
Salehi, Hadi; Ren, Lei; Howard, David
2016-08-01
Computer simulation of human gait, based on measured motion data, is a well-established technique in biomechanics. However, optimisation studies requiring many iterative gait cycle simulations have not yet found widespread application because of their high computational cost. Therefore, a computationally efficient inverse dynamics model of 3D human gait has been designed and compared with an equivalent model, created using a commercial multi-body dynamics package. The fast inverse dynamics model described in this paper led to an eight fold increase in execution speed. Sufficient detail is provided to allow readers to implement the model themselves. PMID:26745213
Analysis of forward and inverse problems in chemical dynamics and spectroscopy
Rabitz, H.
1993-12-01
The overall scope of this research concerns the development and application of forward and inverse analysis tools for problems in chemical dynamics and chemical kinetics. The chemical dynamics work is specifically associated with relating features in potential surfaces and resultant dynamical behavior. The analogous inverse research aims to provide stable algorithms for extracting potential surfaces from laboratory data. In the case of chemical kinetics, the focus is on the development of systematic means to reduce the complexity of chemical kinetic models. Recent progress in these directions is summarized below.
On trajectory generation for flexible space crane: Inverse dynamics analysis by LATDYN
NASA Technical Reports Server (NTRS)
Chen, G.-S.; Housner, J. M.; Wu, S.-C.; Chang, C.-W.
1989-01-01
For future in-space construction facility, one or more space cranes capable of manipulating and positioning large and massive spacecraft components will be needed. Inverse dynamics was extensively studied as a basis for trajectory generation and control of robot manipulators. The focus here is on trajectory generation in the gross-motion phase of space crane operation. Inverse dynamics of the flexible crane body is much more complex and intricate as compared with rigid robot link. To model and solve the space crane's inverse dynamics problem, LATDYN program which employs a three-dimensional finite element formulation for the multibody truss-type structures will be used. The formulation is oriented toward a joint dominated structure which is suitable for the proposed space crane concept. To track a planned trajectory, procedures will be developed to obtain the actuation profile and dynamics envelope which are pertinent to the design and performance requirements of the space crane concept.
Adaptive forward-inverse modeling of reservoir fluids away from wellbores
Ziagos, J P; Gelinas, R J; Doss, S K; Nelson, R G
1999-07-30
This Final Report contains the deliverables of the DeepLook Phase I project entitled, ''Adaptive Forward-Inverse Modeling of Reservoir Fluids Away from Wellbores''. The deliverables are: (i) a description of 2-D test problem results, analyses, and technical descriptions of the techniques used, (ii) a listing of program setup commands that construct and execute the codes for selected test problems (these commands are in mathematical terminology, which reinforces technical descriptions in the text), and (iii) an evaluation and recommendation regarding continuance of this project, including considerations of possible extensions to 3-D codes, additional technical scope, and budget for the out-years. The far-market objective in this project is to develop advanced technologies that can help locate and enhance the recovery of oil from heterogeneous rock formations. The specific technical objective in Phase I was to develop proof-of-concept of new forward and inverse (F-I) modeling techniques [Gelinas et al, 1998] that seek to enhance estimates (images) of formation permeability distributions and fluid motion away from wellbore volumes. This goes to the heart of improving industry's ability to jointly image reservoir permeability and flow predictions of trapped and recovered oil versus time. The estimation of formation permeability away from borehole measurements is an ''inverse'' problem. It is an inseparable part of modeling fluid flows throughout the reservoir in efforts to increase the efficiency of oil recovery at minimum cost. Classic issues of non-uniqueness, mathematical instability, noise effects, and inadequate numerical solution techniques have historically impeded progress in reservoir parameter estimations. Because information pertaining to fluid and rock properties is always sampled sparsely by wellbore measurements, a successful method for interpolating permeability and fluid data between the measurements must be: (i) physics-based, (ii) conditioned by
Dynamics of the inverse MAPLE nanoparticle deposition process
NASA Astrophysics Data System (ADS)
Steiner, Matthew A.; Fitz-Gerald, James M.
2015-05-01
Matrix-assisted pulsed laser evaporation (MAPLE) is a processing technique by which laser-sensitive materials are dissolved or placed into colloidal solution with a strongly absorbing sacrificial solvent, which when frozen into a solid target and irradiated under vacuum disperses the undamaged solute material onto a desired substrate. We present an inversion of the original MAPLE process, where the irradiation of metal-based acetate precursors in solution with UV transparent water results in the deposition of inorganic nanoparticles. A theory is forwarded to explain the underlying multiscale sequence of events that control the inverse MAPLE process from acetate decomposition to nanoparticle formation and subsequent ejection. Support for this theory is provided through the analysis of deposited nanoparticles and by novel characterization of MAPLE targets post-irradiation via cryostage scanning electron microscopy. Ejection is shown to proceed through the same phase-explosion mechanism that drives conventional MAPLE, relating the two techniques and advancing the broader understanding of MAPLE deposition processes.
An inverse dynamics approach to trajectory optimization for an aerospace plane
NASA Technical Reports Server (NTRS)
Lu, Ping
1992-01-01
An inverse dynamics approach for trajectory optimization is proposed. This technique can be useful in many difficult trajectory optimization and control problems. The application of the approach is exemplified by ascent trajectory optimization for an aerospace plane. Both minimum-fuel and minimax types of performance indices are considered. When rocket augmentation is available for ascent, it is shown that accurate orbital insertion can be achieved through the inverse control of the rocket in the presence of disturbances.
A 2-D dynamical model of mesospheric temperature inversions in winter
Hauchecorne, A.; Maillard, A. )
1990-11-01
A 2-D stratospheric and mesospheric dynamical model including drag and diffusion due to gravity wave breaking is used to simulate winter mesospheric temperature inversions similar to those observed by Rayleigh lidar. It is shown that adiabatic heating associated to descending velocities in the mesosphere is the main mechanism involved in the formation of such inversions. Sensitivity tests are performed with the model and confirm this assumption. It is also explained why other previous similar studies with 2-D models did not show mesospheric inversion layers.
Neural network with dynamically adaptable neurons
NASA Technical Reports Server (NTRS)
Tawel, Raoul (Inventor)
1994-01-01
This invention is an adaptive neuron for use in neural network processors. The adaptive neuron participates in the supervised learning phase of operation on a co-equal basis with the synapse matrix elements by adaptively changing its gain in a similar manner to the change of weights in the synapse IO elements. In this manner, training time is decreased by as much as three orders of magnitude.
Can a pseudo-dynamic source inversion approach improve earthquake source imaging?
NASA Astrophysics Data System (ADS)
Zhang, Youbing; Song, Seok Goo; Dalguer, Luis; Clinton, John
2014-05-01
Imaging a high resolution spatio-temporal slip distribution of an earthquake rupture is a core research goal in seismology. In general we expect to obtain a higher quality source image by improving the observational input data (e.g. using more, higher quality, near field stations). However, recent studies show that increasing the surface station density alone does not significantly improve source inversion results (Custodio et al. 2005; Zhang et al. in review). Song et al. (2009) and Song and Dalguer (2013) found interesting correlation structures between kinematic source parameters (e.g. slip, peak slip velocity and rupture velocity) obtained both from kinematic inversion and dynamic modeling. These correlation structures that effectively regularize the model space may improve source imaging more than by simply improving the observational data. In this 'pseudo-dynamic' source inversion, source images are constrained by both physical constraints derived from rupture dynamics as well all the observational data, without compromising the computational efficiency of kinematic inversion. We investigate the efficiency of the pseudo-dynamic source inversion using synthetic dynamic rupture models. Our target model is a buried vertical strike-slip event (Mw 7.3) in a homogeneous half space. In the inversion, we model low frequency (below 1Hz) waveforms using a genetic algorithm in a Bayesian framework (Moneli et al. 2008). A dynamically consistent regularized Yoffe function (Tinti, et al. 2005) was applied as a single-window slip velocity function. We have first implemented the autocorrelation of slip in the prior distribution in the Bayesian inversion - preliminary results show that estimated kinematic source models closely match the target dynamic model. The prior information describing the auto-correlation of source parameters (e.g. slip) improves the imaging of spatial distribution of source parameters. By implementing both auto- and cross-correlation of kinematic
Llacer, J; Solberg, T D; Promberger, C
2001-10-01
This paper presents a description of tests carried out to compare the behaviour of five algorithms in inverse radiation therapy planning: (1) The Dynamically Penalized Likelihood (DPL), an algorithm based on statistical estimation theory; (2) an accelerated version of the same algorithm: (3) a new fast adaptive simulated annealing (ASA) algorithm; (4) a conjugate gradient method; and (5) a Newton gradient method. A three-dimensional mathematical phantom and two clinical cases have been studied in detail. The phantom consisted of a U-shaped tumour with a partially enclosed 'spinal cord'. The clinical examples were a cavernous sinus meningioma and a prostate case. The algorithms have been tested in carefully selected and controlled conditions so as to ensure fairness in the assessment of results. It has been found that all five methods can yield relatively similar optimizations, except when a very demanding optimization is carried out. For the easier cases. the differences are principally in robustness, ease of use and optimization speed. In the more demanding case, there are significant differences in the resulting dose distributions. The accelerated DPL emerges as possibly the algorithm of choice for clinical practice. An appendix describes the differences in behaviour between the new ASA method and the one based on a patent by the Nomos Corporation. PMID:11686280
The Branching Bifurcation of Adaptive Dynamics
NASA Astrophysics Data System (ADS)
Della Rossa, Fabio; Dercole, Fabio; Landi, Pietro
2015-06-01
We unfold the bifurcation involving the loss of evolutionary stability of an equilibrium of the canonical equation of Adaptive Dynamics (AD). The equation deterministically describes the expected long-term evolution of inheritable traits — phenotypes or strategies — of coevolving populations, in the limit of rare and small mutations. In the vicinity of a stable equilibrium of the AD canonical equation, a mutant type can invade and coexist with the present — resident — types, whereas the fittest always win far from equilibrium. After coexistence, residents and mutants effectively diversify, according to the enlarged canonical equation, only if natural selection favors outer rather than intermediate traits — the equilibrium being evolutionarily unstable, rather than stable. Though the conditions for evolutionary branching — the joint effect of resident-mutant coexistence and evolutionary instability — have been known for long, the unfolding of the bifurcation has remained a missing tile of AD, the reason being related to the nonsmoothness of the mutant invasion fitness after branching. In this paper, we develop a methodology that allows the approximation of the invasion fitness after branching in terms of the expansion of the (smooth) fitness before branching. We then derive a canonical model for the branching bifurcation and perform its unfolding around the loss of evolutionary stability. We cast our analysis in the simplest (but classical) setting of asexual, unstructured populations living in an isolated, homogeneous, and constant abiotic environment; individual traits are one-dimensional; intra- as well as inter-specific ecological interactions are described in the vicinity of a stationary regime.
Inversion Of Dynamical Equations For Control Of Attitude
NASA Technical Reports Server (NTRS)
Bach, Ralph; Paielli, Russell
1995-01-01
Method of inverting nonlinear equations of rotational dynamics of rigid body used to design feedback control of orientation of body. Applicable to both direction-cosine and quaternion formulations suitable for large-angle maneuvers. Exploiting some apparently little-known properties of direction cosine and quaternion formulations, method leads to equations for model-follower control system that exhibits exactly linear attitude-error dynamics. Quarternion system more robust in responding to large roll-angle commands.
Classical and quantum dynamics in an inverse square potential
Guillaumín-España, Elisa; Núñez-Yépez, H. N.; Salas-Brito, A. L.
2014-10-15
The classical motion of a particle in a 3D inverse square potential with negative energy, E, is shown to be geodesic, i.e., equivalent to the particle's free motion on a non-compact phase space manifold irrespective of the sign of the coupling constant. We thus establish that all its classical orbits with E < 0 are unbounded. To analyse the corresponding quantum problem, the Schrödinger equation is solved in momentum space. No discrete energy levels exist in the unrenormalized case and the system shows a complete “fall-to-the-center” with an energy spectrum unbounded by below. Such behavior corresponds to the non-existence of bound classical orbits. The symmetry of the problem is SO(3) × SO(2, 1) corroborating previously obtained results.
Static and dynamic responses of an ultrathin adaptive secondary mirror
NASA Astrophysics Data System (ADS)
del Vecchio, Ciro; Brusa, Guido; Gallieni, Daniele; Lloyd-Hart, Michael; Davison, Warren B.
1999-09-01
We present the results of a compete set of static and dynamic runs of the FEA model of the MMT adaptive secondary. The thin mirror is the most delicate component of the MMT adaptive secondary unit, as it provides the deformable optical surface able to correct the incoming wavefront. The static performances are evaluated as a function of the various load cases arising form gravitational loads and from the forces deriving from the magnetic interactions between actuators. In addition, computations were performed to assess the dynamic response to the high bandwidth, adaptive correcting force.s In both cases, the performances of the adaptive mirror design are able to accommodate the severe specifications.
Kalman filtering, smoothing and recursive robot arm forward and inverse dynamics
NASA Technical Reports Server (NTRS)
Rodriguez, G.
1986-01-01
The inverse and forward dynamics problems for multi-link serial manipulators are solved by using recursive techniques from linear filtering and smoothing theory. The pivotal step is to cast the system dynamics and kinematics as a two-point boundary-value problem. Solution of this problem leads to filtering and smoothing techniques identical to the equations of Kalman filtering and Bryson-Frazier fixed time-interval smoothing. The solutions prescribe an inward filtering recursion to compute a sequence of constraint moments and forces followed by an outward recursion to determine a corresponding sequence of angular and linear accelerations. In addition to providing techniques to compute joint accelerations from applied joint moments (and vice versa), the report provides an approach to evaluate recursively the composite multi-link system inertia matrix and its inverse. The report lays the foundation for the potential use of filtering and smoothing techniques in robot inverse and forward dynamics and in robot control design.
Recruitment dynamics in adaptive social networks.
Shkarayev, Maxim S; Schwartz, Ira B; Shaw, Leah B
2013-01-01
We model recruitment in adaptive social networks in the presence of birth and death processes. Recruitment is characterized by nodes changing their status to that of the recruiting class as a result of contact with recruiting nodes. Only a susceptible subset of nodes can be recruited. The recruiting individuals may adapt their connections in order to improve recruitment capabilities, thus changing the network structure adaptively. We derive a mean field theory to predict the dependence of the growth threshold of the recruiting class on the adaptation parameter. Furthermore, we investigate the effect of adaptation on the recruitment level, as well as on network topology. The theoretical predictions are compared with direct simulations of the full system. We identify two parameter regimes with qualitatively different bifurcation diagrams depending on whether nodes become susceptible frequently (multiple times in their lifetime) or rarely (much less than once per lifetime). PMID:25395989
Using a pseudo-dynamic source inversion approach to improve earthquake source imaging
NASA Astrophysics Data System (ADS)
Zhang, Y.; Song, S. G.; Dalguer, L. A.; Clinton, J. F.
2014-12-01
Imaging a high-resolution spatio-temporal slip distribution of an earthquake rupture is a core research goal in seismology. In general we expect to obtain a higher quality source image by improving the observational input data (e.g. using more higher quality near-source stations). However, recent studies show that increasing the surface station density alone does not significantly improve source inversion results (Custodio et al. 2005; Zhang et al. 2014). We introduce correlation structures between the kinematic source parameters: slip, rupture velocity, and peak slip velocity (Song et al. 2009; Song and Dalguer 2013) in the non-linear source inversion. The correlation structures are physical constraints derived from rupture dynamics that effectively regularize the model space and may improve source imaging. We name this approach pseudo-dynamic source inversion. We investigate the effectiveness of this pseudo-dynamic source inversion method by inverting low frequency velocity waveforms from a synthetic dynamic rupture model of a buried vertical strike-slip event (Mw 6.5) in a homogeneous half space. In the inversion, we use a genetic algorithm in a Bayesian framework (Moneli et al. 2008), and a dynamically consistent regularized Yoffe function (Tinti, et al. 2005) was used for a single-window slip velocity function. We search for local rupture velocity directly in the inversion, and calculate the rupture time using a ray-tracing technique. We implement both auto- and cross-correlation of slip, rupture velocity, and peak slip velocity in the prior distribution. Our results suggest that kinematic source model estimates capture the major features of the target dynamic model. The estimated rupture velocity closely matches the target distribution from the dynamic rupture model, and the derived rupture time is smoother than the one we searched directly. By implementing both auto- and cross-correlation of kinematic source parameters, in comparison to traditional smoothing
Identification of dynamic characteristics of flexible rotors as dynamic inverse problem
NASA Technical Reports Server (NTRS)
Roisman, W. P.; Vajingortin, L. D.
1991-01-01
The problem of dynamic and balancing of flexible rotors were considered, which were set and solved as the problem of the identification of flexible rotor systems, which is the same as the inverse problem of the oscillation theory dealing with the task of the identifying the outside influences and system parameters on the basis of the known laws of motion. This approach to the problem allows the disclosure the picture of disbalances throughout the rotor-under-test (which traditional methods of flexible rotor balancing, based on natural oscillations, could not provide), and identify dynamic characteristics of the system, which correspond to a selected mathematical model. Eventually, various methods of balancing were developed depending on the special features of the machines as to their design, technology, and operation specifications. Also, theoretical and practical methods are given for the flexible rotor balancing at far from critical rotation frequencies, which does not necessarily require the knowledge forms of oscillation, dissipation, and elasticity and inertia characteristics, and to use testing masses.
Non-negative constraint research of Tikhonov regularization inversion for dynamic light scattering
NASA Astrophysics Data System (ADS)
Wang, Y. J.; Shen, J.; Liu, W.; Sun, X. M.; Dou, Z. H.
2013-08-01
In dynamic light scattering (DLS) technology, a non-negative constraint on the solution can improve the inversion accuracy of the particle size distribution (PSD). Different non-negative constraint methods have different effects on the inversion results. Combined with the Tikhonov regularization inversion method, the following non-negativity constraint methods: negative to zero (N-to-Z), multi-negative to zero (Multi-N-to-Z), Lin-projected gradient (LPG), oblique projected Landweber (OPL), projected sequential subspace optimization (PSESOP), interior point Newton (IPN), gradient projection conjugate gradient (GPCG) and trust-region method based on the interior reflective Newton (TR-IRN) method are studied in DLS inversion. In different inversion ranges and noise levels, autocorrelation functions of unimodal and bimodal particle distributions were inverted using different non-negativity constraint methods. From the inversion results, the characteristics of the various methods were obtained, which can be treated as a reference for the implementation of non-negative constraints in Tikhonov regularization inversion of DLS.
Persistent inversion dynamics and wintertime PM10 air pollution in Alpine valleys
NASA Astrophysics Data System (ADS)
Largeron, Yann; Staquet, Chantal
2016-06-01
The present study investigates persistent inversions dynamics during a whole winter in Alpine valleys of the area of Grenoble (French Alps), and their relationship to PM10 air pollution episodes and synoptic scale meteorology. For this purpose, hourly time series from November to March of PM10 concentration measurements at the bottom of the valleys and of ground-based temperature data at different altitudes are used. A methodology is developed to quantify a simple estimate of the inversion strength from temperature profiles deduced from the ground-based observations. This estimate is shown to be equivalent to the boundary layer heat deficit. A criterion based on this estimate is proposed to identify persistent (more than 3 days) inversions. Persistent inversions are found to occur from November to February and span 35% of the time. It is shown that they are closely related to PM10 pollution episodes, the PM10 concentration increasing with the boundary layer stability as the inversion develops. Polluted episodes are primarily driven by persistent inversions and consequently, pollution is of fully local origin from November to February. In March local dynamics become less important and long-range transport can dominate. Persistent inversions occur systematically during a high-pressure regime, which first triggers a synoptic scale elevated inversion due to the advection of warm air masses in the mid-troposphere. In valleys, the sheltered boundary layer becomes decoupled from the free troposphere, which allows a ground-based inversion to intensify in the following days. An inversion layer of quasi-constant temperature gradient, greater than 5 K km-1, then forms up to an altitude of about 1600 m, close to the average elevation of the summits. If the episode is sufficiently long, a stagnation stage is reached during which daytime insolation produces a shallow convective surface layer which does not destroy the persistent inversion. The inversion break-up occurs rapidly
Dynamic mesh adaption for triangular and tetrahedral grids
NASA Technical Reports Server (NTRS)
Biswas, Rupak; Strawn, Roger
1993-01-01
The following topics are discussed: requirements for dynamic mesh adaption; linked-list data structure; edge-based data structure; adaptive-grid data structure; three types of element subdivision; mesh refinement; mesh coarsening; additional constraints for coarsening; anisotropic error indicator for edges; unstructured-grid Euler solver; inviscid 3-D wing; and mesh quality for solution-adaptive grids. The discussion is presented in viewgraph form.
NASA Astrophysics Data System (ADS)
Grayver, Alexander V.
2015-07-01
This paper presents a distributed magnetotelluric inversion scheme based on adaptive finite-element method (FEM). The key novel aspect of the introduced algorithm is the use of automatic mesh refinement techniques for both forward and inverse modelling. These techniques alleviate tedious and subjective procedure of choosing a suitable model parametrization. To avoid overparametrization, meshes for forward and inverse problems were decoupled. For calculation of accurate electromagnetic (EM) responses, automatic mesh refinement algorithm based on a goal-oriented error estimator has been adopted. For further efficiency gain, EM fields for each frequency were calculated using independent meshes in order to account for substantially different spatial behaviour of the fields over a wide range of frequencies. An automatic approach for efficient initial mesh design in inverse problems based on linearized model resolution matrix was developed. To make this algorithm suitable for large-scale problems, it was proposed to use a low-rank approximation of the linearized model resolution matrix. In order to fill a gap between initial and true model complexities and resolve emerging 3-D structures better, an algorithm for adaptive inverse mesh refinement was derived. Within this algorithm, spatial variations of the imaged parameter are calculated and mesh is refined in the neighborhoods of points with the largest variations. A series of numerical tests were performed to demonstrate the utility of the presented algorithms. Adaptive mesh refinement based on the model resolution estimates provides an efficient tool to derive initial meshes which account for arbitrary survey layouts, data types, frequency content and measurement uncertainties. Furthermore, the algorithm is capable to deliver meshes suitable to resolve features on multiple scales while keeping number of unknowns low. However, such meshes exhibit dependency on an initial model guess. Additionally, it is demonstrated
Method and system for training dynamic nonlinear adaptive filters which have embedded memory
NASA Technical Reports Server (NTRS)
Rabinowitz, Matthew (Inventor)
2002-01-01
Described herein is a method and system for training nonlinear adaptive filters (or neural networks) which have embedded memory. Such memory can arise in a multi-layer finite impulse response (FIR) architecture, or an infinite impulse response (IIR) architecture. We focus on filter architectures with separate linear dynamic components and static nonlinear components. Such filters can be structured so as to restrict their degrees of computational freedom based on a priori knowledge about the dynamic operation to be emulated. The method is detailed for an FIR architecture which consists of linear FIR filters together with nonlinear generalized single layer subnets. For the IIR case, we extend the methodology to a general nonlinear architecture which uses feedback. For these dynamic architectures, we describe how one can apply optimization techniques which make updates closer to the Newton direction than those of a steepest descent method, such as backpropagation. We detail a novel adaptive modified Gauss-Newton optimization technique, which uses an adaptive learning rate to determine both the magnitude and direction of update steps. For a wide range of adaptive filtering applications, the new training algorithm converges faster and to a smaller value of cost than both steepest-descent methods such as backpropagation-through-time, and standard quasi-Newton methods. We apply the algorithm to modeling the inverse of a nonlinear dynamic tracking system 5, as well as a nonlinear amplifier 6.
Analytical simulation and inversion of dynamic urban land surface effects
NASA Astrophysics Data System (ADS)
Bayer, P.; Rivera, J.; Blum, P.; Schweizer, D.; Rybach, L.
2015-12-01
Long-term thermal changes at the land surface can be backtracked from borehole temperature profiles. The main focus so far has been on past climate changes, assuming perfect coupling of surface air and ground temperature. In many urbanized areas, however, temperature profiles are heavily perturbed. We find a characteristic bending of urban profiles towards shallow depth, which indicates strong heating from the ground surface during recent decades. This phenomenon is generally described as subsurface urban heat island (UHI) effect, which exists beneath many cities worldwide. Major drivers are land use changes and urban structures that act as long-term heat sources that artificially load the top 100 m of the ground. While variability in land use and coverage are critical factors for reliable borehole climatology, temperature profiles can also be inverted to trace back the combined effect of past urbanization and climate. We present an analytical framework based on the superposition of specific Green's functions for simulating transient land use changes and their effects on borehole temperature profiles. By inversion in a Bayesian framework, flexible calibration of unknown spatially distributed parameter values and their correlation is feasible. The procedure is applied to four temperature logs which are around 200-400 m deep from the city and suburbs of Zurich, Switzerland. These were recorded recently by a temperature sensor and data logger introduced in closed borehole heat exchangers before the start of geothermal operation. At the sites, long-term land use changes are well documented for more than the last century. This facilitated focusing on a few unknown parameters, and we selected the contribution by asphalt and by basements of buildings. It is revealed that for three of the four sites, these two factors dominate the subsurface UHI evolution. At one site, additional factors such as buried district heating networks may play a role. It is demonstrated that site
Reinbolt, Jeffrey A.; Haftka, Raphael T.; Chmielewski, Terese L.; Fregly, Benjamin J.
2013-01-01
Variations in joint parameter values (axis positions and orientations in body segments) and inertial parameter values (segment masses, mass centers, and moments of inertia) as well as kinematic noise alter the results of inverse dynamics analyses of gait. Three-dimensional linkage models with joint constraints have been proposed as one way to minimize the effects of noisy kinematic data. Such models can also be used to perform gait optimizations to predict post-treatment function given pre-treatment gait data. This study evaluates whether accurate patient-specific joint and inertial parameter values are needed in three-dimensional linkage models to produce accurate inverse dynamics results for gait. The study was performed in two stages. First, we used optimization analyses to evaluate whether patient-specific joint and inertial parameter values can be calibrated accurately from noisy kinematic data, and second, we used Monte Carlo analyses to evaluate how errors in joint and inertial parameter values affect inverse dynamics calculations. Both stages were performed using a dynamic, 27 degree-of-freedom, full-body linkage model and synthetic (i.e., computer generated) gait data corresponding to a nominal experimental gait motion. In general, joint but not inertial parameter values could be found accurately from noisy kinematic data. Root-mean-square (RMS) errors were 3° and 4 mm for joint parameter values and 1 kg, 22 mm, and 74,500 kg*mm2 for inertial parameter values. Furthermore, errors in joint but not inertial parameter values had a significant effect on calculated lower-extremity inverse dynamics joint torques. The worst RMS torque error averaged 4% bodyweight*height (BW*H) due to joint parameter variations but less than 0.25% BW*H due to inertial parameter variations. These results suggest that inverse dynamics analyses of gait utilizing linkage models with joint constraints should calibrate the model’s joint parameter values to obtain accurate joint
Inverse Dynamics Control of Constrained Robots in the Presence of Joint Flexibility
NASA Astrophysics Data System (ADS)
IDER, S. KEMAL
1999-07-01
An inverse dynamics control algorithm for constrained flexible-joint robots is developed. It is shown that in a flexible-joint robot, the acceleration level inverse dynamic equations are singular because of the elastic media. Implicit numerical integration methods that account for the higher order derivative information are utilized for solving the singular set of differential equations. The control law proposed linearizes and decouples the system and achieves simultaneous and asymptotically stable trajectory tracking control of the end-effector motion and contact forces. Together with the integrators for improving robustness due to modelling errors and disturbances, a fifth order position error dynamics and a third order contact force error dynamics are obtained. A 3R spatial robot with all joints flexible is simulated to illustrate the performance of the method.
Cheng, Ching-An; Huang, Han-Pang; Hsu, Huan-Kun; Lai, Wei-Zh; Cheng, Chih-Chun
2016-07-01
We investigate the modeling of inverse dynamics without prior kinematic information for holonomic rigid-body robots. Despite success in compensating robot dynamics and friction, general inverse dynamics models are nontrivial. Rigid-body models are restrictive or inefficient; learning-based models are generalizable yet require large training data. The structured kernels address the dilemma by embedding the robot dynamics in reproducing kernel Hilbert space. The proposed kernels autonomously converge to rigid-body models but require fewer samples; with a semi-parametric framework that incorporates additional parametric basis for friction, the structured kernels can efficiently model general rigid-body robots. We tested the proposed scheme in simulations and experiments; the models that consider the structure of function space are more accurate. PMID:26316286
Control of a high beta maneuvering reentry vehicle using dynamic inversion.
Watts, Alfred Chapman
2005-05-01
The design of flight control systems for high performance maneuvering reentry vehicles presents a significant challenge to the control systems designer. These vehicles typically have a much higher ballistic coefficient than crewed vehicles like as the Space Shuttle or proposed crew return vehicles such as the X-38. Moreover, the missions of high performance vehicles usually require a steeper reentry flight path angle, followed by a pull-out into level flight. These vehicles then must transit the entire atmosphere and robustly perform the maneuvers required for the mission. The vehicles must also be flown with small static margins in order to perform the required maneuvers, which can result in highly nonlinear aerodynamic characteristics that frequently transition from being aerodynamically stable to unstable as angle of attack increases. The control system design technique of dynamic inversion has been applied successfully to both high performance aircraft and low beta reentry vehicles. The objective of this study was to explore the application of this technique to high performance maneuvering reentry vehicles, including the basic derivation of the dynamic inversion technique, followed by the extension of that technique to the use of tabular trim aerodynamic models in the controller. The dynamic inversion equations are developed for high performance vehicles and augmented to allow the selection of a desired response for the control system. A six degree of freedom simulation is used to evaluate the performance of the dynamic inversion approach, and results for both nominal and off nominal aerodynamic characteristics are presented.
Multithreaded Model for Dynamic Load Balancing Parallel Adaptive PDE Computations
NASA Technical Reports Server (NTRS)
Chrisochoides, Nikos
1995-01-01
We present a multithreaded model for the dynamic load-balancing of numerical, adaptive computations required for the solution of Partial Differential Equations (PDE's) on multiprocessors. Multithreading is used as a means of exploring concurrency in the processor level in order to tolerate synchronization costs inherent to traditional (non-threaded) parallel adaptive PDE solvers. Our preliminary analysis for parallel, adaptive PDE solvers indicates that multithreading can be used an a mechanism to mask overheads required for the dynamic balancing of processor workloads with computations required for the actual numerical solution of the PDE's. Also, multithreading can simplify the implementation of dynamic load-balancing algorithms, a task that is very difficult for traditional data parallel adaptive PDE computations. Unfortunately, multithreading does not always simplify program complexity, often makes code re-usability not an easy task, and increases software complexity.
Light-Directed Dynamic Chirality Inversion in Functional Self-Organized Helical Superstructures.
Bisoyi, Hari Krishna; Li, Quan
2016-02-24
Helical superstructures are widely observed in nature, in synthetic polymers, and in supramolecular assemblies. Controlling the chirality (the handedness) of dynamic helical superstructures of molecular and macromolecular systems by external stimuli is a challenging task, but is of great fundamental significance with appealing morphology-dependent applications. Light-driven chirality inversion in self-organized helical superstructures (i.e. cholesteric, chiral nematic liquid crystals) is currently in the limelight because inversion of the handedness alters the chirality of the circularly polarized light that they selectively reflect, which has wide potential for application. Here we discuss the recent developments toward inversion of the handedness of cholesteric liquid crystals enabled by photoisomerizable chiral molecular switches or motors. Different classes of chiral photoresponsive dopants (guests) capable of conferring light-driven reversible chirality inversion of helical superstructures fabricated from different nematic hosts are discussed. Rational molecular designs of chiral molecular switches toward endowing handedness inversion to the induced helical superstructures of cholesteric liquid crystals are highlighted. This Review is concluded by throwing light on the challenges and opportunities in this emerging frontier, and it is expected to provide useful guidelines toward the development of self-organized soft materials with stimuli-directed chirality inversion capability and multifunctional host-guest systems. PMID:26764018
Quantitative Adaptation Analytics for Assessing Dynamic Systems of Systems.
Gauthier, John H.; Miner, Nadine E.; Wilson, Michael L.; Le, Hai D.; Kao, Gio K; Melander, Darryl J.; Longsine, Dennis Earl; Vander Meer, Robert Charles,
2015-01-01
Our society is increasingly reliant on systems and interoperating collections of systems, known as systems of systems (SoS). These SoS are often subject to changing missions (e.g., nation- building, arms-control treaties), threats (e.g., asymmetric warfare, terrorism), natural environments (e.g., climate, weather, natural disasters) and budgets. How well can SoS adapt to these types of dynamic conditions? This report details the results of a three year Laboratory Directed Research and Development (LDRD) project aimed at developing metrics and methodologies for quantifying the adaptability of systems and SoS. Work products include: derivation of a set of adaptability metrics, a method for combining the metrics into a system of systems adaptability index (SoSAI) used to compare adaptability of SoS designs, development of a prototype dynamic SoS (proto-dSoS) simulation environment which provides the ability to investigate the validity of the adaptability metric set, and two test cases that evaluate the usefulness of a subset of the adaptability metrics and SoSAI for distinguishing good from poor adaptability in a SoS. Intellectual property results include three patents pending: A Method For Quantifying Relative System Adaptability, Method for Evaluating System Performance, and A Method for Determining Systems Re-Tasking.
Dynamics of adaptive agents with asymmetric information
NASA Astrophysics Data System (ADS)
DeMartino, Andrea; Galla, Tobias
2005-08-01
We apply path integral techniques to study the dynamics of agent-based models with asymmetric information structures. In particular, we devise a batch version of a model proposed originally by Berg et al (2001 Quantitative Finance 1 203), and convert the coupled multi-agent processes into an effective-agent problem from which the dynamical order parameters in ergodic regimes can be derived self-consistently together with the corresponding phase structure. Our dynamical study complements and extends the available static theory. Results are confirmed by numerical simulations.
Bayesian inversions of a dynamic vegetation model at four European grassland sites
NASA Astrophysics Data System (ADS)
Minet, J.; Laloy, E.; Tychon, B.; Francois, L.
2015-05-01
Eddy covariance data from four European grassland sites are used to probabilistically invert the CARAIB (CARbon Assimilation In the Biosphere) dynamic vegetation model (DVM) with 10 unknown parameters, using the DREAM(ZS) (DiffeRential Evolution Adaptive Metropolis) Markov chain Monte Carlo (MCMC) sampler. We focus on comparing model inversions, considering both homoscedastic and heteroscedastic eddy covariance residual errors, with variances either fixed a priori or jointly inferred together with the model parameters. Agreements between measured and simulated data during calibration are comparable with previous studies, with root mean square errors (RMSEs) of simulated daily gross primary productivity (GPP), ecosystem respiration (RECO) and evapotranspiration (ET) ranging from 1.73 to 2.19, 1.04 to 1.56 g C m-2 day-1 and 0.50 to 1.28 mm day-1, respectively. For the calibration period, using a homoscedastic eddy covariance residual error model resulted in a better agreement between measured and modelled data than using a heteroscedastic residual error model. However, a model validation experiment showed that CARAIB models calibrated considering heteroscedastic residual errors perform better. Posterior parameter distributions derived from using a heteroscedastic model of the residuals thus appear to be more robust. This is the case even though the classical linear heteroscedastic error model assumed herein did not fully remove heteroscedasticity of the GPP residuals. Despite the fact that the calibrated model is generally capable of fitting the data within measurement errors, systematic bias in the model simulations are observed. These are likely due to model inadequacies such as shortcomings in the photosynthesis modelling. Besides the residual error treatment, differences between model parameter posterior distributions among the four grassland sites are also investigated. It is shown that the marginal distributions of the specific leaf area and characteristic
Complexity and network dynamics in physiological adaptation: an integrated view.
Baffy, György; Loscalzo, Joseph
2014-05-28
Living organisms constantly interact with their surroundings and sustain internal stability against perturbations. This dynamic process follows three fundamental strategies (restore, explore, and abandon) articulated in historical concepts of physiological adaptation such as homeostasis, allostasis, and the general adaptation syndrome. These strategies correspond to elementary forms of behavior (ordered, chaotic, and static) in complex adaptive systems and invite a network-based analysis of the operational characteristics, allowing us to propose an integrated framework of physiological adaptation from a complex network perspective. Applicability of this concept is illustrated by analyzing molecular and cellular mechanisms of adaptation in response to the pervasive challenge of obesity, a chronic condition resulting from sustained nutrient excess that prompts chaotic exploration for system stability associated with tradeoffs and a risk of adverse outcomes such as diabetes, cardiovascular disease, and cancer. Deconstruction of this complexity holds the promise of gaining novel insights into physiological adaptation in health and disease. PMID:24751342
Dynamics of adaptive structures: Design through simulations
NASA Technical Reports Server (NTRS)
Park, K. C.; Alexander, S.
1993-01-01
The use of a helical bi-morph actuator/sensor concept by mimicking the change of helical waveform in bacterial flagella is perhaps the first application of bacterial motions (living species) to longitudinal deployment of space structures. However, no dynamical considerations were analyzed to explain the waveform change mechanisms. The objective is to review various deployment concepts from the dynamics point of view and introduce the dynamical considerations from the outset as part of design considerations. Specifically, the impact of the incorporation of the combined static mechanisms and dynamic design considerations on the deployment performance during the reconfiguration stage is studied in terms of improved controllability, maneuvering duration, and joint singularity index. It is shown that intermediate configurations during articulations play an important role for improved joint mechanisms design and overall structural deployability.
Stability Result For Dynamic Inversion Devised to Control Large Flexible Aircraft
NASA Technical Reports Server (NTRS)
Gregory, Irene M.
2001-01-01
High performance aircraft of the future will be designed lighter, more maneuverable, and operate over an ever expanding flight envelope. One of the largest differences from the flight control perspective between current and future advanced aircraft is elasticity. Over the last decade, dynamic inversion methodology has gained considerable popularity in application to highly maneuverable fighter aircraft, which were treated as rigid vehicles. This paper is an initial attempt to establish global stability results for dynamic inversion methodology as applied to a large, flexible aircraft. This work builds on a previous result for rigid fighter aircraft and adds a new level of complexity that is the flexible aircraft dynamics, which cannot be ignored even in the most basic flight control. The results arise from observations of the control laws designed for a new generation of the High-Speed Civil Transport aircraft.
Lai, Guanyu; Liu, Zhi; Zhang, Yun; Philip Chen, C L
2016-06-01
This paper is concentrated on the problem of adaptive fuzzy tracking control for an uncertain nonlinear system whose actuator is encountered by the asymmetric backlash behavior. First, we propose a new smooth inverse model which can approximate the asymmetric actuator backlash arbitrarily. By applying it, two adaptive fuzzy control scenarios, namely, the compensation-based control scheme and nonlinear decomposition-based control scheme, are then developed successively. It is worth noticing that the first fuzzy controller exhibits a better tracking control performance, although it recourses to a known slope ratio of backlash nonlinearity. The second one further removes the restriction, and also gets a desirable control performance. By the strict Lyapunov argument, both adaptive fuzzy controllers guarantee that the output tracking error is convergent to an adjustable region of zero asymptotically, while all the signals remain semiglobally uniformly ultimately bounded. Lastly, two comparative simulations are conducted to verify the effectiveness of the proposed fuzzy controllers. PMID:27187937
Analog forecasting with dynamics-adapted kernels
NASA Astrophysics Data System (ADS)
Zhao, Zhizhen; Giannakis, Dimitrios
2016-09-01
Analog forecasting is a nonparametric technique introduced by Lorenz in 1969 which predicts the evolution of states of a dynamical system (or observables defined on the states) by following the evolution of the sample in a historical record of observations which most closely resembles the current initial data. Here, we introduce a suite of forecasting methods which improve traditional analog forecasting by combining ideas from kernel methods developed in harmonic analysis and machine learning and state-space reconstruction for dynamical systems. A key ingredient of our approach is to replace single-analog forecasting with weighted ensembles of analogs constructed using local similarity kernels. The kernels used here employ a number of dynamics-dependent features designed to improve forecast skill, including Takens’ delay-coordinate maps (to recover information in the initial data lost through partial observations) and a directional dependence on the dynamical vector field generating the data. Mathematically, our approach is closely related to kernel methods for out-of-sample extension of functions, and we discuss alternative strategies based on the Nyström method and the multiscale Laplacian pyramids technique. We illustrate these techniques in applications to forecasting in a low-order deterministic model for atmospheric dynamics with chaotic metastability, and interannual-scale forecasting in the North Pacific sector of a comprehensive climate model. We find that forecasts based on kernel-weighted ensembles have significantly higher skill than the conventional approach following a single analog.
Yavari, Fatemeh; Mahdavi, Shirin; Towhidkhah, Farzad; Ahmadi-Pajouh, Mohammad-Ali; Ekhtiari, Hamed; Darainy, Mohammad
2016-04-01
Despite several pieces of evidence, which suggest that the human brain employs internal models for motor control and learning, the location of these models in the brain is not yet clear. In this study, we used transcranial direct current stimulation (tDCS) to manipulate right cerebellar function, while subjects adapt to a visuomotor task. We investigated the effect of this manipulation on the internal forward and inverse models by measuring two kinds of behavior: generalization of training in one direction to neighboring directions (as a proxy for inverse models) and localization of the hand position after movement without visual feedback (as a proxy for forward model). The experimental results showed no effect of cerebellar tDCS on generalization, but significant effect on localization. These observations support the idea that the cerebellum is a possible brain region for internal forward, but not inverse model formation. We also used a realistic human head model to calculate current density distribution in the brain. The result of this model confirmed the passage of current through the cerebellum. Moreover, to further explain some observed experimental results, we modeled the visuomotor adaptation process with the help of a biologically inspired method known as population coding. The effect of tDCS was also incorporated in the model. The results of this modeling study closely match our experimental data and provide further evidence in line with the idea that tDCS manipulates FM's function in the cerebellum. PMID:26706039
Adaptive dynamic FBG interrogation utilising erbium-doped fibre
NASA Astrophysics Data System (ADS)
John, R. N.; Read, I.; MacPherson, W. N.
2013-04-01
A dynamic fibre Bragg grating interrogation scheme is investigated using two-wave mixing in erbium-doped fibre, capable of adapting to quasistatic strain and temperature drifts. An interference pattern set up in the erbium-doped fibre creates, due to the photorefractive effect, a dynamic grating capable of wavelength demodulating the FBG signal. The presence of a dynamic grating was verified and then dynamic strain signals from a fibre stretcher were measured. The adaptive nature of the technique was successfully demonstrated by heating the FBG while it underwent dynamic straining leading to detection unlike an alternative arrayed waveguide grating system which simultaneously failed detection. Two gratings were then wavelength division multiplexed with the signal grating receiving approximately 30dB greater signal showing that there was little cross talk in the system.
Tran, Anh Phuong; Dafflon, Baptiste; Hubbard, Susan S.; Kowalsky, Michael B.; Long, Philip; Tokunaga, Tetsu K.; Williams, Kenneth H.
2016-08-31
Improving our ability to estimate the parameters that control water and heat fluxes in the shallow subsurface is particularly important due to their strong control on recharge, evaporation and biogeochemical processes. The objectives of this study are to develop and test a new inversion scheme to simultaneously estimate subsurface hydrological, thermal and petrophysical parameters using hydrological, thermal and electrical resistivity tomography (ERT) data. The inversion scheme – which is based on a nonisothermal, multiphase hydrological model – provides the desired subsurface property estimates in high spatiotemporal resolution. A particularly novel aspect of the inversion scheme is the explicit incorporation of themore » dependence of the subsurface electrical resistivity on both moisture and temperature. The scheme was applied to synthetic case studies, as well as to real datasets that were autonomously collected at a biogeochemical field study site in Rifle, Colorado. At the Rifle site, the coupled hydrological-thermal-geophysical inversion approach well predicted the matric potential, temperature and apparent resistivity with the Nash–Sutcliffe efficiency criterion greater than 0.92. Synthetic studies found that neglecting the subsurface temperature variability, and its effect on the electrical resistivity in the hydrogeophysical inversion, may lead to an incorrect estimation of the hydrological parameters. The approach is expected to be especially useful for the increasing number of studies that are taking advantage of autonomously collected ERT and soil measurements to explore complex terrestrial system dynamics.« less
NASA Astrophysics Data System (ADS)
Goldman, S. P.; Turnbull, D.; Johnson, C.; Chen, J. Z.; Battista, J. J.
2009-05-01
A fast, accurate and stable optimization algorithm is very important for inverse planning of intensity-modulated radiation therapy (IMRT), and for implementing dose-adaptive radiotherapy in the future. Conventional numerical search algorithms with positive beam weight constraints generally require numerous iterations and may produce suboptimal dose results due to trapping in local minima regions of the objective function landscape. A direct solution of the inverse problem using conventional quadratic objective functions without positive beam constraints is more efficient but it will result in unrealistic negative beam weights. We review here a direct solution of the inverse problem that is efficient and does not yield unphysical negative beam weights. In fast inverse dose optimization (FIDO) method the objective function for the optimization of a large number of beamlets is reformulated such that the optimization problem is reducible to a linear set of equations. The optimal set of intensities is then found through a matrix inversion, and negative beamlet intensities are avoided without the need for externally imposed ad hoc conditions. In its original version [S. P. Goldman, J. Z. Chen, and J. J. Battista, in Proceedings of the XIVth International Conference on the Use of Computers in Radiation Therapy, 2004, pp. 112-115; S. P. Goldman, J. Z. Chen, and J. J. Battista, Med. Phys. 32, 3007 (2005)], FIDO was tested on single two-dimensional computed tomography (CT) slices with sharp KERMA beams without scatter, in order to establish a proof of concept which demonstrated that FIDO could be a viable method for the optimization of cancer treatment plans. In this paper we introduce the latest advancements in FIDO that now include not only its application to three-dimensional volumes irradiated by beams with full scatter but include as well a complete implementation of clinical dose-volume constraints including maximum and minimum dose as well as equivalent uniform dose
Robust adaptive dynamic programming with an application to power systems.
Jiang, Yu; Jiang, Zhong-Ping
2013-07-01
This brief presents a novel framework of robust adaptive dynamic programming (robust-ADP) aimed at computing globally stabilizing and suboptimal control policies in the presence of dynamic uncertainties. A key strategy is to integrate ADP theory with techniques in modern nonlinear control with a unique objective of filling up a gap in the past literature of ADP without taking into account dynamic uncertainties. Neither the system dynamics nor the system order are required to be precisely known. As an illustrative example, the computational algorithm is applied to the controller design of a two-machine power system. PMID:24808528
NASA Astrophysics Data System (ADS)
Foks, Nathan Leon
The interpretation of geophysical data plays an important role in the analysis of potential field data in resource exploration industries. Two categories of interpretation techniques are discussed in this thesis; boundary detection and geophysical inversion. Fault or boundary detection is a method to interpret the locations of subsurface boundaries from measured data, while inversion is a computationally intensive method that provides 3D information about subsurface structure. My research focuses on these two aspects of interpretation techniques. First, I develop a method to aid in the interpretation of faults and boundaries from magnetic data. These processes are traditionally carried out using raster grid and image processing techniques. Instead, I use unstructured meshes of triangular facets that can extract inferred boundaries using mesh edges. Next, to address the computational issues of geophysical inversion, I develop an approach to reduce the number of data in a data set. The approach selects the data points according to a user specified proxy for its signal content. The approach is performed in the data domain and requires no modification to existing inversion codes. This technique adds to the existing suite of compressive inversion algorithms. Finally, I develop an algorithm to invert gravity data for an interfacing surface using an unstructured mesh of triangular facets. A pertinent property of unstructured meshes is their flexibility at representing oblique, or arbitrarily oriented structures. This flexibility makes unstructured meshes an ideal candidate for geometry based interface inversions. The approaches I have developed provide a suite of algorithms geared towards large-scale interpretation of potential field data, by using an unstructured representation of both the data and model parameters.
Performance evaluation of the inverse dynamics method for optimal spacecraft reorientation
NASA Astrophysics Data System (ADS)
Ventura, Jacopo; Romano, Marcello; Walter, Ulrich
2015-05-01
This paper investigates the application of the inverse dynamics in the virtual domain method to Euler angles, quaternions, and modified Rodrigues parameters for rapid optimal attitude trajectory generation for spacecraft reorientation maneuvers. The impact of the virtual domain and attitude representation is numerically investigated for both minimum time and minimum energy problems. Owing to the nature of the inverse dynamics method, it yields sub-optimal solutions for minimum time problems. Furthermore, the virtual domain improves the optimality of the solution, but at the cost of more computational time. The attitude representation also affects solution quality and computational speed. For minimum energy problems, the optimal solution can be obtained without the virtual domain with any considered attitude representation.
A forward-muscular inverse-skeletal dynamics framework for human musculoskeletal simulations.
S Shourijeh, Mohammad; Smale, Kenneth B; Potvin, Brigitte M; Benoit, Daniel L
2016-06-14
This study provides a forward-muscular inverse-skeletal dynamics framework for musculoskeletal simulations. The simulation framework works based on solving the muscle redundancy problem forward in time parallel to a torque tracking between the musculotendon net torques and joint moments from inverse dynamics. The proposed framework can be used by any musculoskeletal modeling software package; however, just to exemplify, here in this study it is wrapped around OpenSim and the optimization is done in MATLAB. The novel simulation framework was highly robust for repeated runs and produced relatively high correlations between predicted muscle excitations and experimental EMGs for level gait trials. This simulation framework represents an efficient and robust approach to predict muscle excitation, musculotendon unit force, and to estimate net joint torque. PMID:27106173
On-line identification of forward/inverse systems for adaptive control applications
NASA Technical Reports Server (NTRS)
Horta, Lucas G.; Sandridge, Chris A.
1992-01-01
The paper is concerned with the on-line identification of system Markov parameters using observers. Two types of parameters are sought: the forward system parameters, which provide estimates of the outputs for given inputs, and inverse system parameters, which provide estimates of inputs for given outputs. A procedure is proposed which uses the inverse system to generate commands that produce a desired system response. The feasibility of the procedure is demonstrated using test results from an implementation of a truss structure fitted with piezoelectric actuators and collocated strain gauges.
Saha, Debajyoti Kumar Shaw, Pankaj; Janaki, M. S.; Sekar Iyengar, A. N.; Ghosh, Sabuj; Mitra, Vramori Michael Wharton, Alpha
2014-03-15
Order-chaos-order was observed in the relaxation oscillations of a glow discharge plasma with variation in the discharge voltage. The first transition exhibits an inverse homoclinic bifurcation followed by a homoclinic bifurcation in the second transition. For the two regimes of observations, a detailed analysis of correlation dimension, Lyapunov exponent, and Renyi entropy was carried out to explore the complex dynamics of the system.
Inverse kinematic and forward dynamic models of the 2002 Denali fault earthquake, Alaska
Oglesby, D.D.; Dreger, Douglas S.; Harris, R.A.; Ratchkovski, N.; Hansen, R.
2004-01-01
We perform inverse kinematic and forward dynamic models of the M 7.9 2002 Denali fault, Alaska, earthquake to shed light on the rupture process and dynamics of this event, which took place on a geometrically complex fault system in central Alaska. We use a combination of local seismic and Global Positioning System (GPS) data for our kinematic inversion and find that the slip distribution of this event is characterized by three major asperities on the Denali fault. The rupture nucleated on the Susitna Glacier thrust fault, and after a pause, propagated onto the strike-slip Denali fault. Approximately 216 km to the east, the rupture abandoned the Denali fault in favor of the more southwesterly directed Totschunda fault. Three-dimensional dynamic models of this event indicate that the abandonment of the Denali fault for the Totschunda fault can be explained by the Totschunda fault's more favorable orientation with respect to the local stress field. However, a uniform tectonic stress field cannot explain the complex slip pattern in this event. We also find that our dynamic models predict discontinuous rupture from the Denali to Totschunda fault segments. Such discontinuous rupture helps to qualitatively improve our kinematic inverse models. Two principal implications of our study are (1) a combination of inverse and forward modeling can bring insight into earthquake processes that are not possible with either technique alone, and (2) the stress field on geometrically complex fault systems is most likely not due to a uniform tectonic stress field that is resolved onto fault segments of different orientations; rather, other forms of stress heterogeneity must be invoked to explain the observed slip patterns.
Adaptive planning for applications with dynamic objectives
NASA Technical Reports Server (NTRS)
Hadavi, Khosrow; Hsu, Wen-Ling; Pinedo, Michael
1992-01-01
We devise a qualitative control layer to be integrated into a real-time multi-agent reactive planner. The reactive planning system consists of distributed planning agents attending to various perspectives of the task environment. Each perspective corresponds to an objective. The set of objectives considered are sometimes in conflict with each other. Each agent receives information about events as they occur, and a set of actions based on heuristics can be taken by the agents. Within the qualitative control scheme, we use a set of qualitative feature vectors to describe the effects of applying actions. A qualitative transition vector is used to denote the qualitative distance between the current state and the target state. We will then apply on-line learning at the qualitative control level to achieve adaptive planning. Our goal is to design a mechanism to refine the heuristics used by the reactive planner every time an action is taken toward achieving the objectives, using feedback from the results of the actions. When the outcome is compared with expectations, our prior objectives may be modified and a new set of objectives (or a new assessment of the relative importance of the different objectives) can be introduced. Because we are able to obtain better estimates of the time-varying objectives, the reactive strategies can be improved and better prediction can be achieved.
Feasible muscle activation ranges based on inverse dynamics analyses of human walking.
Simpson, Cole S; Sohn, M Hongchul; Allen, Jessica L; Ting, Lena H
2015-09-18
Although it is possible to produce the same movement using an infinite number of different muscle activation patterns owing to musculoskeletal redundancy, the degree to which observed variations in muscle activity can deviate from optimal solutions computed from biomechanical models is not known. Here, we examined the range of biomechanically permitted activation levels in individual muscles during human walking using a detailed musculoskeletal model and experimentally-measured kinetics and kinematics. Feasible muscle activation ranges define the minimum and maximum possible level of each muscle's activation that satisfy inverse dynamics joint torques assuming that all other muscles can vary their activation as needed. During walking, 73% of the muscles had feasible muscle activation ranges that were greater than 95% of the total muscle activation range over more than 95% of the gait cycle, indicating that, individually, most muscles could be fully active or fully inactive while still satisfying inverse dynamics joint torques. Moreover, the shapes of the feasible muscle activation ranges did not resemble previously-reported muscle activation patterns nor optimal solutions, i.e. static optimization and computed muscle control, that are based on the same biomechanical constraints. Our results demonstrate that joint torque requirements from standard inverse dynamics calculations are insufficient to define the activation of individual muscles during walking in healthy individuals. Identifying feasible muscle activation ranges may be an effective way to evaluate the impact of additional biomechanical and/or neural constraints on possible versus actual muscle activity in both normal and impaired movements. PMID:26300401
Robust adaptive dynamic programming and feedback stabilization of nonlinear systems.
Jiang, Yu; Jiang, Zhong-Ping
2014-05-01
This paper studies the robust optimal control design for a class of uncertain nonlinear systems from a perspective of robust adaptive dynamic programming (RADP). The objective is to fill up a gap in the past literature of adaptive dynamic programming (ADP) where dynamic uncertainties or unmodeled dynamics are not addressed. A key strategy is to integrate tools from modern nonlinear control theory, such as the robust redesign and the backstepping techniques as well as the nonlinear small-gain theorem, with the theory of ADP. The proposed RADP methodology can be viewed as an extension of ADP to uncertain nonlinear systems. Practical learning algorithms are developed in this paper, and have been applied to the controller design problems for a jet engine and a one-machine power system. PMID:24808035
Dynamics and Adaptive Control for Stability Recovery of Damaged Aircraft
NASA Technical Reports Server (NTRS)
Nguyen, Nhan; Krishnakumar, Kalmanje; Kaneshige, John; Nespeca, Pascal
2006-01-01
This paper presents a recent study of a damaged generic transport model as part of a NASA research project to investigate adaptive control methods for stability recovery of damaged aircraft operating in off-nominal flight conditions under damage and or failures. Aerodynamic modeling of damage effects is performed using an aerodynamic code to assess changes in the stability and control derivatives of a generic transport aircraft. Certain types of damage such as damage to one of the wings or horizontal stabilizers can cause the aircraft to become asymmetric, thus resulting in a coupling between the longitudinal and lateral motions. Flight dynamics for a general asymmetric aircraft is derived to account for changes in the center of gravity that can compromise the stability of the damaged aircraft. An iterative trim analysis for the translational motion is developed to refine the trim procedure by accounting for the effects of the control surface deflection. A hybrid direct-indirect neural network, adaptive flight control is proposed as an adaptive law for stabilizing the rotational motion of the damaged aircraft. The indirect adaptation is designed to estimate the plant dynamics of the damaged aircraft in conjunction with the direct adaptation that computes the control augmentation. Two approaches are presented 1) an adaptive law derived from the Lyapunov stability theory to ensure that the signals are bounded, and 2) a recursive least-square method for parameter identification. A hardware-in-the-loop simulation is conducted and demonstrates the effectiveness of the direct neural network adaptive flight control in the stability recovery of the damaged aircraft. A preliminary simulation of the hybrid adaptive flight control has been performed and initial data have shown the effectiveness of the proposed hybrid approach. Future work will include further investigations and high-fidelity simulations of the proposed hybrid adaptive Bight control approach.
Kalman filtering, smoothing, and recursive robot arm forward and inverse dynamics
NASA Technical Reports Server (NTRS)
Rodriguez, Guillermo
1987-01-01
The recursive difference equations of Kalman filtering and Bryson-Frazier fixed time-interval smoothing, arising in the state estimation theory for linear state space systems, are used here to solve problems of serial manipulator inverse and forward dynamics. The configuration analyzed is that of a joint connected N-link serial manipulator attached to an immobile base. The equivalence between the filtering and smoothing techniques from state estimation theory and recursive robot dynamics methods is demonstrated. Several areas for future research are suggested.
Inverse Dynamics Model for the Ankle Joint with Applications in Tibia Malleolus Fracture
NASA Astrophysics Data System (ADS)
Budescu, E.; Merticaru, E.; Chirazi, M.
The paper presents a biomechanical model of the ankle joint, in order to determine the force and the torque of reaction into the articulation, through inverse dynamic analysis, in various stages of the gait. Thus, knowing the acceleration of the foot and the reaction force between foot and ground during the gait, determined by experimental measurement, there was calculated, for five different positions of the foot, the joint reaction forces, on the basis of dynamic balance equations. The values numerically determined were compared with the admissible forces appearing in the technical systems of osteosynthesis of tibia malleolus fracture, in order to emphasize the motion restrictions during bone healing.
Sex speeds adaptation by altering the dynamics of molecular evolution.
McDonald, Michael J; Rice, Daniel P; Desai, Michael M
2016-03-10
Sex and recombination are pervasive throughout nature despite their substantial costs. Understanding the evolutionary forces that maintain these phenomena is a central challenge in biology. One longstanding hypothesis argues that sex is beneficial because recombination speeds adaptation. Theory has proposed several distinct population genetic mechanisms that could underlie this advantage. For example, sex can promote the fixation of beneficial mutations either by alleviating interference competition (the Fisher-Muller effect) or by separating them from deleterious load (the ruby in the rubbish effect). Previous experiments confirm that sex can increase the rate of adaptation, but these studies did not observe the evolutionary dynamics that drive this effect at the genomic level. Here we present the first, to our knowledge, comparison between the sequence-level dynamics of adaptation in experimental sexual and asexual Saccharomyces cerevisiae populations, which allows us to identify the specific mechanisms by which sex speeds adaptation. We find that sex alters the molecular signatures of evolution by changing the spectrum of mutations that fix, and confirm theoretical predictions that it does so by alleviating clonal interference. We also show that substantially deleterious mutations hitchhike to fixation in adapting asexual populations. In contrast, recombination prevents such mutations from fixing. Our results demonstrate that sex both speeds adaptation and alters its molecular signature by allowing natural selection to more efficiently sort beneficial from deleterious mutations. PMID:26909573
Gradient-based adaptation of continuous dynamic model structures
NASA Astrophysics Data System (ADS)
La Cava, William G.; Danai, Kourosh
2016-01-01
A gradient-based method of symbolic adaptation is introduced for a class of continuous dynamic models. The proposed model structure adaptation method starts with the first-principles model of the system and adapts its structure after adjusting its individual components in symbolic form. A key contribution of this work is its introduction of the model's parameter sensitivity as the measure of symbolic changes to the model. This measure, which is essential to defining the structural sensitivity of the model, not only accommodates algebraic evaluation of candidate models in lieu of more computationally expensive simulation-based evaluation, but also makes possible the implementation of gradient-based optimisation in symbolic adaptation. The proposed method is applied to models of several virtual and real-world systems that demonstrate its potential utility.
NASA Astrophysics Data System (ADS)
Muta, Osamu; Akaiwa, Yoshihiko
In this paper, we propose a simple peak power reduction (PPR) method based on adaptive inversion of parity-check block of codeword in BCH-coded OFDM system. In the proposed method, the entire parity-check block of the codeword is adaptively inversed by multiplying weighting factors (WFs) so as to minimize PAPR of the OFDM signal, symbol-by-symbol. At the receiver, these WFs are estimated based on the property of BCH decoding. When the primitive BCH code with single error correction such as (31,26) code is used, to estimate the WFs, the proposed method employs a significant bit protection method which assigns a significant bit to the best subcarrier selected among all possible subcarriers. With computer simulation, when (31,26), (31,21) and (32,21) BCH codes are employed, PAPR of the OFDM signal at the CCDF (Complementary Cumulative Distribution Function) of 10-4 is reduced by about 1.9, 2.5 and 2.5dB by applying the PPR method, while achieving the BER performance comparable to the case with the perfect WF estimation in exponentially decaying 12-path Rayleigh fading condition.
NASA Astrophysics Data System (ADS)
Yamada, Masumi; Mangeney, Anne; Matsushi, Yuki; Moretti, Laurent
2016-06-01
We performed numerical simulations of the 2011 deep-seated Akatani landslide in central Japan to understand the dynamic evolution of friction of the landslide. By comparing the forces obtained from numerical simulation to those resolved from seismic waveform inversion, the coefficient of the friction during sliding was investigated in the range of 0.1 to 0.4. The simulation assuming standard Coulomb friction shows that the forces obtained by the seismic waveform inversion are well explained using a constant friction of μ = 0.3. A small difference between the residuals of Coulomb simulation and a velocity-dependent simulation suggests that the coefficient of friction over the volume is well constrained as 0.3 most of time during sliding. It suggests the sudden loss of shearing resistance at the onset of sliding, i.e., sudden drop of the initial coefficient of friction in our model, which accelerates the deep-seated landslide. Our numerical simulation calibrated by seismic data provides the evolution of dynamic friction with a reasonable resolution in time, which is difficult to obtain from a conventional runout simulation, or seismic waveform inversion alone.
Kemppainen, Petri; Knight, Christopher G; Sarma, Devojit K; Hlaing, Thaung; Prakash, Anil; Maung Maung, Yan Naung; Somboon, Pradya; Mahanta, Jagadish; Walton, Catherine
2015-09-01
Recent advances in sequencing allow population-genomic data to be generated for virtually any species. However, approaches to analyse such data lag behind the ability to generate it, particularly in nonmodel species. Linkage disequilibrium (LD, the nonrandom association of alleles from different loci) is a highly sensitive indicator of many evolutionary phenomena including chromosomal inversions, local adaptation and geographical structure. Here, we present linkage disequilibrium network analysis (LDna), which accesses information on LD shared between multiple loci genomewide. In LD networks, vertices represent loci, and connections between vertices represent the LD between them. We analysed such networks in two test cases: a new restriction-site-associated DNA sequence (RAD-seq) data set for Anopheles baimaii, a Southeast Asian malaria vector; and a well-characterized single nucleotide polymorphism (SNP) data set from 21 three-spined stickleback individuals. In each case, we readily identified five distinct LD network clusters (single-outlier clusters, SOCs), each comprising many loci connected by high LD. In A. baimaii, further population-genetic analyses supported the inference that each SOC corresponds to a large inversion, consistent with previous cytological studies. For sticklebacks, we inferred that each SOC was associated with a distinct evolutionary phenomenon: two chromosomal inversions, local adaptation, population-demographic history and geographic structure. LDna is thus a useful exploratory tool, able to give a global overview of LD associated with diverse evolutionary phenomena and identify loci potentially involved. LDna does not require a linkage map or reference genome, so it is applicable to any population-genomic data set, making it especially valuable for nonmodel species. PMID:25573196
Kemppainen, Petri; Knight, Christopher G; Sarma, Devojit K; Hlaing, Thaung; Prakash, Anil; Maung Maung, Yan Naung; Somboon, Pradya; Mahanta, Jagadish; Walton, Catherine
2015-01-01
Recent advances in sequencing allow population-genomic data to be generated for virtually any species. However, approaches to analyse such data lag behind the ability to generate it, particularly in nonmodel species. Linkage disequilibrium (LD, the nonrandom association of alleles from different loci) is a highly sensitive indicator of many evolutionary phenomena including chromosomal inversions, local adaptation and geographical structure. Here, we present linkage disequilibrium network analysis (LDna), which accesses information on LD shared between multiple loci genomewide. In LD networks, vertices represent loci, and connections between vertices represent the LD between them. We analysed such networks in two test cases: a new restriction-site-associated DNA sequence (RAD-seq) data set for Anopheles baimaii, a Southeast Asian malaria vector; and a well-characterized single nucleotide polymorphism (SNP) data set from 21 three-spined stickleback individuals. In each case, we readily identified five distinct LD network clusters (single-outlier clusters, SOCs), each comprising many loci connected by high LD. In A. baimaii, further population-genetic analyses supported the inference that each SOC corresponds to a large inversion, consistent with previous cytological studies. For sticklebacks, we inferred that each SOC was associated with a distinct evolutionary phenomenon: two chromosomal inversions, local adaptation, population-demographic history and geographic structure. LDna is thus a useful exploratory tool, able to give a global overview of LD associated with diverse evolutionary phenomena and identify loci potentially involved. LDna does not require a linkage map or reference genome, so it is applicable to any population-genomic data set, making it especially valuable for nonmodel species. PMID:25573196
Content-adaptive ghost imaging of dynamic scenes.
Li, Ziwei; Suo, Jinli; Hu, Xuemei; Dai, Qionghai
2016-04-01
Limited by long acquisition time of 2D ghost imaging, current ghost imaging systems are so far inapplicable for dynamic scenes. However, it's been demonstrated that nature images are spatiotemporally redundant and the redundancy is scene dependent. Inspired by that, we propose a content-adaptive computational ghost imaging approach to achieve high reconstruction quality under a small number of measurements, and thus achieve ghost imaging of dynamic scenes. To utilize content-adaptive inter-frame redundancy, we put the reconstruction under an iterative reweighted optimization, with non-uniform weight computed from temporal-correlated frame sequences. The proposed approach can achieve dynamic imaging at 16fps with 64×64-pixel resolution. PMID:27137022
Analysis and inverse substructuring computation on dynamic quality of mechanical assembly
NASA Astrophysics Data System (ADS)
Lü, Guangqing; Yi, Chuijie; Fang, Ke
2016-05-01
Mechanical assembly has its own dynamic quality directly affecting the dynamic quality of whole product and should be considered in quality inspection and estimation of mechanical assembly. Based on functional relations between dynamic characteristics involved in mechanical assembly, the effects of assembling process on dynamic characteristics of substructural components of an assembly system are investigated by substructuring analysis. Assembly-coupling dynamic stiffness is clarified as the dominant factor of the effects and can be used as a quantitative measure of assembly dynamic quality. Two computational schemes using frequency response functions(FRFs) to determine the stiffness are provided and discussed by inverse substructuring analysis, including their applicable conditions and implementation procedure in application. Eigenvalue analysis on matrix-ratios of FRFs before and after assembling is employed and well validates the analytical outcomes and the schemes via both a lumped-parameter model and its analogic experimental counterpart. Applying the two schemes to inspect the dynamic quality provides the message of dynamic performance of the assembly system, and therefore improves conventional quality inspection and estimation of mechanical assembly in completeness.
Analysis and inverse substructuring computation on dynamic quality of mechanical assembly
NASA Astrophysics Data System (ADS)
Lü, Guangqing; Yi, Chuijie; Fang, Ke
2016-04-01
Mechanical assembly has its own dynamic quality directly affecting the dynamic quality of whole product and should be considered in quality inspection and estimation of mechanical assembly. Based on functional relations between dynamic characteristics involved in mechanical assembly, the effects of assembling process on dynamic characteristics of substructural components of an assembly system are investigated by substructuring analysis. Assembly-coupling dynamic stiffness is clarified as the dominant factor of the effects and can be used as a quantitative measure of assembly dynamic quality. Two computational schemes using frequency response functions(FRFs) to determine the stiffness are provided and discussed by inverse substructuring analysis, including their applicable conditions and implementation procedure in application. Eigenvalue analysis on matrix-ratios of FRFs before and after assembling is employed and well validates the analytical outcomes and the schemes via both a lumped-parameter model and its analogic experimental counterpart. Applying the two schemes to inspect the dynamic quality provides the message of dynamic performance of the assembly system, and therefore improves conventional quality inspection and estimation of mechanical assembly in completeness.
Serial and parallel dynamic adaptation of general hybrid meshes
NASA Astrophysics Data System (ADS)
Kavouklis, Christos
The Navier-Stokes equations are a standard mathematical representation of viscous fluid flow. Their numerical solution in three dimensions remains a computationally intensive and challenging task, despite recent advances in computer speed and memory. A strategy to increase accuracy of Navier-Stokes simulations, while maintaining computing resources to a minimum, is local refinement of the associated computational mesh in regions of large solution gradients and coarsening in regions where the solution does not vary appreciably. In this work we consider adaptation of general hybrid meshes for Computational Fluid Dynamics (CFD) applications. Hybrid meshes are composed of four types of elements; hexahedra, prisms, pyramids and tetrahedra, and have been proven a promising technology in accurately resolving fluid flow for complex geometries. The first part of this dissertation is concerned with the design and implementation of a serial scheme for the adaptation of general three dimensional hybrid meshes. We have defined 29 refinement types, for all four kinds of elements. The core of the present adaptation scheme is an iterative algorithm that flags mesh edges for refinement, so that the adapted mesh is conformal. Of primary importance is considered the design of a suitable dynamic data structure that facilitates refinement and coarsening operations and furthermore minimizes memory requirements. A special dynamic list is defined for mesh elements, in contrast with the usual tree structures. It contains only elements of the current adaptation step and minimal information that is utilized to reconstruct parent elements when the mesh is coarsened. In the second part of this work, a new parallel dynamic mesh adaptation and load balancing algorithm for general hybrid meshes is presented. Partitioning of a hybrid mesh reduces to partitioning of the corresponding dual graph. Communication among processors is based on the faces of the interpartition boundary. The distributed
Coherent source imaging and dynamic support tracking for inverse scattering using compressive MUSIC
NASA Astrophysics Data System (ADS)
Lee, Okkyun; Kim, Jong Min; Yoo, Jaejoon; Jin, Kyunghwan; Ye, Jong Chul
2011-09-01
The goal of this paper is to develop novel algorithms for inverse scattering problems such as EEG/MEG, microwave imaging, and/or diffuse optical tomograpahy, and etc. One of the main contributions of this paper is a class of novel non-iterative exact nonlinear inverse scattering theory for coherent source imaging and moving targets. Specifically, the new algorithms guarantee the exact recovery under a very relaxed constraint on the number of source and receivers, under which the conventional methods fail. Such breakthrough was possible thanks to the recent theory of compressive MUSIC and its extension using support correction criterion, where partial support are estimated using the conventional compressed sensing approaches, then the remaining supports are estimated using a novel generalized MUSIC criterion. Numerical results using coherent sources in EEG/MEG and dynamic targets confirm that the new algorithms outperform the conventional ones.
Dynamic analysis of neural encoding by point process adaptive filtering.
Eden, Uri T; Frank, Loren M; Barbieri, Riccardo; Solo, Victor; Brown, Emery N
2004-05-01
Neural receptive fields are dynamic in that with experience, neurons change their spiking responses to relevant stimuli. To understand how neural systems adapt their representations of biological information, analyses of receptive field plasticity from experimental measurements are crucial. Adaptive signal processing, the well-established engineering discipline for characterizing the temporal evolution of system parameters, suggests a framework for studying the plasticity of receptive fields. We use the Bayes' rule Chapman-Kolmogorov paradigm with a linear state equation and point process observation models to derive adaptive filters appropriate for estimation from neural spike trains. We derive point process filter analogues of the Kalman filter, recursive least squares, and steepest-descent algorithms and describe the properties of these new filters. We illustrate our algorithms in two simulated data examples. The first is a study of slow and rapid evolution of spatial receptive fields in hippocampal neurons. The second is an adaptive decoding study in which a signal is decoded from ensemble neural spiking activity as the receptive fields of the neurons in the ensemble evolve. Our results provide a paradigm for adaptive estimation for point process observations and suggest a practical approach for constructing filtering algorithms to track neural receptive field dynamics on a millisecond timescale. PMID:15070506
Parallel tetrahedral mesh adaptation with dynamic load balancing
Oliker, Leonid; Biswas, Rupak; Gabow, Harold N.
2000-06-28
The ability to dynamically adapt an unstructured grid is a powerful tool for efficiently solving computational problems with evolving physical features. In this paper, we report on our experience parallelizing an edge-based adaptation scheme, called 3D-TAG, using message passing. Results show excellent speedup when a realistic helicopter rotor mesh is randomly refined. However, performance deteriorates when the mesh is refined using a solution-based error indicator since mesh adaptation for practical problems occurs in a localized region, creating a severe load imbalance. To address this problem, we have developed PLUM, a global dynamic load balancing framework for adaptive numerical computations. Even though PLUM primarily balances processor workloads for the solution phase, it reduces the load imbalance problem within mesh adaptation by repartitioning the mesh after targeting edges for refinement but before the actual subdivision. This dramatically improves the performance of parallel 3D-TAG since refinement occurs in a more load balanced fashion. We also present optimal and heuristic algorithms that, when applied to the default mapping of a parallel repartitioner, significantly reduce the data redistribution overhead. Finally, portability is examined by comparing performance on three state-of-the-art parallel machines.
Model-adaptive hybrid dynamic control for robotic assembly tasks
Austin, D.J.; McCarragher, B.J.
1999-10-01
A new task-level adaptive controller is presented for the hybrid dynamic control of robotic assembly tasks. Using a hybrid dynamic model of the assembly task, velocity constraints are derived from which satisfactory velocity commands are obtained. Due to modeling errors and parametric uncertainties, the velocity commands may be erroneous and may result in suboptimal performance. Task-level adaptive control schemes, based on the occurrence of discrete events, are used to change the model parameters from which the velocity commands are determined. Two adaptive schemes are presented: the first is based on intuitive reasoning about the vector spaces involved whereas the second uses a search region that is reduced with each iteration. For the first adaptation law, asymptotic convergence to the correct model parameters is proven except for one case. This weakness motivated the development of the second adaptation law, for which asymptotic convergence is proven in all cases. Automated control of a peg-in-hole assembly task is given as an example, and simulations and experiments for this task are presented. These results demonstrate the success of the method and also indicate properties for rapid convergence.
Parallel Tetrahedral Mesh Adaptation with Dynamic Load Balancing
NASA Technical Reports Server (NTRS)
Oliker, Leonid; Biswas, Rupak; Gabow, Harold N.
1999-01-01
The ability to dynamically adapt an unstructured grid is a powerful tool for efficiently solving computational problems with evolving physical features. In this paper, we report on our experience parallelizing an edge-based adaptation scheme, called 3D_TAG. using message passing. Results show excellent speedup when a realistic helicopter rotor mesh is randomly refined. However. performance deteriorates when the mesh is refined using a solution-based error indicator since mesh adaptation for practical problems occurs in a localized region., creating a severe load imbalance. To address this problem, we have developed PLUM, a global dynamic load balancing framework for adaptive numerical computations. Even though PLUM primarily balances processor workloads for the solution phase, it reduces the load imbalance problem within mesh adaptation by repartitioning the mesh after targeting edges for refinement but before the actual subdivision. This dramatically improves the performance of parallel 3D_TAG since refinement occurs in a more load balanced fashion. We also present optimal and heuristic algorithms that, when applied to the default mapping of a parallel repartitioner, significantly reduce the data redistribution overhead. Finally, portability is examined by comparing performance on three state-of-the-art parallel machines.
Eco-evolutionary feedbacks, adaptive dynamics and evolutionary rescue theory
Ferriere, Regis; Legendre, Stéphane
2013-01-01
Adaptive dynamics theory has been devised to account for feedbacks between ecological and evolutionary processes. Doing so opens new dimensions to and raises new challenges about evolutionary rescue. Adaptive dynamics theory predicts that successive trait substitutions driven by eco-evolutionary feedbacks can gradually erode population size or growth rate, thus potentially raising the extinction risk. Even a single trait substitution can suffice to degrade population viability drastically at once and cause ‘evolutionary suicide’. In a changing environment, a population may track a viable evolutionary attractor that leads to evolutionary suicide, a phenomenon called ‘evolutionary trapping’. Evolutionary trapping and suicide are commonly observed in adaptive dynamics models in which the smooth variation of traits causes catastrophic changes in ecological state. In the face of trapping and suicide, evolutionary rescue requires that the population overcome evolutionary threats generated by the adaptive process itself. Evolutionary repellors play an important role in determining how variation in environmental conditions correlates with the occurrence of evolutionary trapping and suicide, and what evolutionary pathways rescue may follow. In contrast with standard predictions of evolutionary rescue theory, low genetic variation may attenuate the threat of evolutionary suicide and small population sizes may facilitate escape from evolutionary traps. PMID:23209163
Dynamical singularities in adaptive delayed-feedback control.
Saito, Asaki; Konishi, Keiji
2011-09-01
We demonstrate the dynamical characteristics of adaptive delayed-feedback control systems, exploiting a discrete-time adaptive control method derived for carrying out detailed analysis. In particular, the systems exhibit singularities such as power-law decay of the distribution of transient times and almost zero finite-time Lyapunov exponents. We can explain these results by characterizing such systems as having (1) a Jacobian matrix with unity eigenvalue in the whole phase space, and (2) parameters approaching a stability boundary proven to be identical with that of (nonadaptive) delayed-feedback control. PMID:22060398
Kinematic and dynamic inversion of the 16 December earthquake in Northern Chile
NASA Astrophysics Data System (ADS)
Ruiz, S.; Lancieri, M.; Madariaga, R. I.; Sobiesiak, M.; Campos, J. A.
2009-12-01
We study the kinematic and dynamic rupture propagation of the M 6.7, intraplate, intermediate depth, slab push earthquake that occurred 16 December 2007, a month after the large interplate thrust event of Tocopilla, Chile (M 7.7). The occurrence of a slab push event after a large subduction earthquake is well explained by Coulomb stress transfer and crack dynamics. A dense seismic network, equipped with short period and accelerometers was deployed after the event of 14 November 2007 by the Task Force of GFZ Potsdam and the University of Chile in Santiago. This network was in place on December 16 providing the best seismic data set ever recorded for a Chilean earthquake. We have used it to do a detailed study of rupture processes. We localized the main event of December 16 and the aftershocks that occurred within 24 h of the main event. The main event was located at 43 km depth, while the aftershocks distribution covered a circular zone of 5 to 8 km of radius centered on the main shock epicenter and with depth ranging between [39 - 49] km. The aftershocks are distributed on an almost vertical plane that agrees with the almost vertical plane of the fault mechanism (86° dip) and all the aftershock have the same mechanism as the main event. We used eight of the nearest accelerometric records low pass filtered at 1 Hz, two of which were situated right above the hypocenter. We performed a non-linear kinematic inversion based on the neighborhood algorithm (NA) with an L2 norm. The velocity model was derived from previous work by GFZ. The earthquake is very well modeled by a circular rupture of radius between 5 and 8 km that propagated with a very low rupture velocity, that varies between 1 and 2 km/s. We need only a few non-linear parameters to model this event, parameter space has a dimension close to 6. The kinematic solution was validated using a full dynamic inversion method in which the rupture process is modeled using finite differences on a coarse grid with a slip
Dynamics of a many-particle Landau-Zener model: Inverse sweep
Itin, A. P.
2009-05-15
We consider dynamics of a slowly time-dependent Dicke model, which represents a many-body generalization of the Landau-Zener model. In particular, the model describes narrow Feshbach resonance passage in an ultracold gas of Fermi atoms. Adiabaticity is destroyed when a parameter crosses a critical value, even at very slow sweeping rates of a parameter. The dynamics crucially depends on direction of the sweep. We apply our recent analysis (A. P. Itin and P. Toermae, e-print arXiv:0901.4778) to the 'inverse' sweep through the resonance, corresponding (in a context of Feshbach resonance passage) to dissociation of molecules. On a level of the mean-field approximation, the dynamics is equivalent to a molecular condensate formation from Bose atoms within a two-mode model. Mapping the system to a Painleve equation allows us to calculate deviation from adiabaticity at very slow sweeps analytically.
Adaptive network dynamics and evolution of leadership in collective migration
NASA Astrophysics Data System (ADS)
Pais, Darren; Leonard, Naomi E.
2014-01-01
The evolution of leadership in migratory populations depends not only on costs and benefits of leadership investments but also on the opportunities for individuals to rely on cues from others through social interactions. We derive an analytically tractable adaptive dynamic network model of collective migration with fast timescale migration dynamics and slow timescale adaptive dynamics of individual leadership investment and social interaction. For large populations, our analysis of bifurcations with respect to investment cost explains the observed hysteretic effect associated with recovery of migration in fragmented environments. Further, we show a minimum connectivity threshold above which there is evolutionary branching into leader and follower populations. For small populations, we show how the topology of the underlying social interaction network influences the emergence and location of leaders in the adaptive system. Our model and analysis can be extended to study the dynamics of collective tracking or collective learning more generally. Thus, this work may inform the design of robotic networks where agents use decentralized strategies that balance direct environmental measurements with agent interactions.
Direct and Inverse Problems of Item Pool Design for Computerized Adaptive Testing
ERIC Educational Resources Information Center
Belov, Dmitry I.; Armstrong, Ronald D.
2009-01-01
The recent literature on computerized adaptive testing (CAT) has developed methods for creating CAT item pools from a large master pool. Each CAT pool is designed as a set of nonoverlapping forms reflecting the skill levels of an assumed population of test takers. This article presents a Monte Carlo method to obtain these CAT pools and discusses…
Dynamic inversion of a Slab-push earthquake in Northern Chile
NASA Astrophysics Data System (ADS)
Ruiz, Sergio; Madariaga, Raul; Lancieri, Maria; Sobesiak, Monika
2010-05-01
We study the dynamic rupture propagation of a M 6.7 intraplate earthquake that occurred 16 December 2007, a month after a large thrust event of Tocopilla, Chile (M 7.7). The occurrence of a slab push event after a large subduction earthquake is well explained by Coulomb stress transfer models and crack dynamics. A dense seismic network, equipped with short period and accelerometers was deployed after the event of 14 November 2007 by the Task Force of GFZ Potsdam and the University of Chile in Santiago. This network was in place on December 16 providing an excellent data set for this earthquake. We used these data to make a detailed study of rupture processes. We localized the main event of December 16 and the aftershocks that occurred within 24 h of the main event. The main event was located at 43 km depth, while the aftershocks distribution covered a circular zone of 5 to 8 km of radius centred on the main shock epicentre. The aftershocks are distributed on an almost vertical plane that agrees with one of the fault planes of the mechanism (86° dip) and all the aftershock have the same mechanism as the main event. We used nearest accelerometric records in order to do dynamic inversion, two of these accelerometers were situated right above the hypocentre. We performed a non-linear dynamic inversion based on the neighbourhood algorithm (NA) and MonteCarlo methods with an L2 norm. The data was initially filtered in the 0.05-1 Hz. The velocity model was derived from previous work by GFZ. The earthquake was modelled using finite differences on a grid of variable size. Friction was modelled by the standard Ida slip weakening friction law. At each step of the inversion more than 32 full numerical simulations are carried in parallel. These simulations have been optimized in order to reduce the computer time to a minimum. The best models that result from dynamic inversion reduced the variance by more than 30 %, these models ruptured a relatively small zone of the fault
Shkvarko, Yuriy; Tuxpan, José; Santos, Stewart
2011-01-01
We consider a problem of high-resolution array radar/SAR imaging formalized in terms of a nonlinear ill-posed inverse problem of nonparametric estimation of the power spatial spectrum pattern (SSP) of the random wavefield scattered from a remotely sensed scene observed through a kernel signal formation operator and contaminated with random Gaussian noise. First, the Sobolev-type solution space is constructed to specify the class of consistent kernel SSP estimators with the reproducing kernel structures adapted to the metrics in such the solution space. Next, the "model-free" variational analysis (VA)-based image enhancement approach and the "model-based" descriptive experiment design (DEED) regularization paradigm are unified into a new dynamic experiment design (DYED) regularization framework. Application of the proposed DYED framework to the adaptive array radar/SAR imaging problem leads to a class of two-level (DEED-VA) regularized SSP reconstruction techniques that aggregate the kernel adaptive anisotropic windowing with the projections onto convex sets to enforce the consistency and robustness of the overall iterative SSP estimators. We also show how the proposed DYED regularization method may be considered as a generalization of the MVDR, APES and other high-resolution nonparametric adaptive radar sensing techniques. A family of the DYED-related algorithms is constructed and their effectiveness is finally illustrated via numerical simulations. PMID:22163859
Combined dynamic inversion and QFT flight control of an unstable high performance aircraft
NASA Astrophysics Data System (ADS)
Stout, Perry Walter
Quantitative Feedback Theory (QFT) is a control system synthesis, technique that directly considers system uncertainties and disturbance magnitudes when formulating closed-loop control algorithms. Dynamic Inversion is a nonlinear control system design technique that relies on accurate mathematical models to compute control inputs producing arbitrary system responses. Both techniques have been applied to unstable high performance aircraft flight control, and produced effective aircraft controllers. Both techniques have certain drawbacks: Nonlinear QFT controllers tend to be unnecessarily conservative (the computed controllers have excessive bandwidth) because known system properties are treated as "unknown" disturbances during loop synthesis. Meanwhile Dynamic Inversion control is sensitive to differences between assumed mathematical models and actual system dynamic properties. Combining the two control techniques provides the benefit of both while suffering the drawbacks of neither, as demonstrated by Single Input, Single Output (SISO) control of a constant airspeed, no roll, no yaw nonlinear model of the F-16 aircraft, and by Multi-Input, Multi-Output (MIMO) control of a full six-degree-of-freedom version. Design performance of the combined controllers is verified by reduced actuator efforts and by reduced sensor noise to actuator input (U( s)/n(s)) transfer function magnitudes compared to standard QFT versions.
Lehikoinen, A.; Huttunen, J.M.J.; Finsterle, S.; Kowalsky, M.B.; Kaipio, J.P.
2009-08-01
We propose an approach for imaging the dynamics of complex hydrological processes. The evolution of electrically conductive fluids in porous media is imaged using time-lapse electrical resistance tomography. The related dynamic inversion problem is solved using Bayesian filtering techniques, that is, it is formulated as a sequential state estimation problem in which the target is an evolving posterior probability density of the system state. The dynamical inversion framework is based on the state space representation of the system, which involves the construction of a stochastic evolution model and an observation model. The observation model used in this paper consists of the complete electrode model for ERT, with Archie's law relating saturations to electrical conductivity. The evolution model is an approximate model for simulating flow through partially saturated porous media. Unavoidable modeling and approximation errors in both the observation and evolution models are considered by computing approximate statistics for these errors. These models are then included in the construction of the posterior probability density of the estimated system state. This approximation error method allows the use of approximate - and therefore computationally efficient - observation and evolution models in the Bayesian filtering. We consider a synthetic example and show that the incorporation of an explicit model for the model uncertainties in the state space representation can yield better estimates than a frame-by-frame imaging approach.
Negre, Christian F A; Mniszewski, Susan M; Cawkwell, Marc J; Bock, Nicolas; Wall, Michael E; Niklasson, Anders M N
2016-07-12
We present a reduced complexity algorithm to compute the inverse overlap factors required to solve the generalized eigenvalue problem in a quantum-based molecular dynamics (MD) simulation. Our method is based on the recursive, iterative refinement of an initial guess of Z (inverse square root of the overlap matrix S). The initial guess of Z is obtained beforehand by using either an approximate divide-and-conquer technique or dynamical methods, propagated within an extended Lagrangian dynamics from previous MD time steps. With this formulation, we achieve long-term stability and energy conservation even under the incomplete, approximate, iterative refinement of Z. Linear-scaling performance is obtained using numerically thresholded sparse matrix algebra based on the ELLPACK-R sparse matrix data format, which also enables efficient shared-memory parallelization. As we show in this article using self-consistent density-functional-based tight-binding MD, our approach is faster than conventional methods based on the diagonalization of overlap matrix S for systems as small as a few hundred atoms, substantially accelerating quantum-based simulations even for molecular structures of intermediate size. For a 4158-atom water-solvated polyalanine system, we find an average speedup factor of 122 for the computation of Z in each MD step. PMID:27267207
Inversion of Dynamical Scattering from Large-Angle Rocking-Beam Electron Diffraction Patterns
NASA Astrophysics Data System (ADS)
Wang, Feng; Pennington, Robert S.; Koch, Christoph T.
2016-07-01
A method for ab initio structure factor retrieval from large-angle rocking-beam electron diffraction data of thin crystals is described and tested with experimental and simulated data. No additional information, such as atomicity or information about chemical composition, has been made use of. Our numerical experiments show that the inversion of dynamical scattering works best, if the beam tilt range is large and the specimen not too thick, because for moderate multiple scattering, the large tilt amplitude effectively removes local minima in this global optimization problem.
Conn, Charlotte E.; Ces, Oscar; Mulet, Xavier; Seddon, John M.; Templer, Richard H.; Finet, Stephanie; Winter, Roland
2006-03-17
The liquid crystalline lamellar (L{sub {alpha}}) to double-diamond inverse bicontinuous cubic (Q{sub II}{sup D}) phase transition for the amphiphile monoelaidin in excess water exhibits a remarkable sequence of structural transformations for pressure or temperature jumps. Our data imply that the transition dynamics depends on a coupling between changes in molecular shape and the geometrical and topological constraints of domain size. We propose a qualitative model for this coupling based on theories of membrane fusion via stalks and existing knowledge of the structure and energetics of bicontinuous cubic phases.
Inversion of Dynamical Scattering from Large-Angle Rocking-Beam Electron Diffraction Patterns.
Wang, Feng; Pennington, Robert S; Koch, Christoph T
2016-07-01
A method for ab initio structure factor retrieval from large-angle rocking-beam electron diffraction data of thin crystals is described and tested with experimental and simulated data. No additional information, such as atomicity or information about chemical composition, has been made use of. Our numerical experiments show that the inversion of dynamical scattering works best, if the beam tilt range is large and the specimen not too thick, because for moderate multiple scattering, the large tilt amplitude effectively removes local minima in this global optimization problem. PMID:27419576
NASA Technical Reports Server (NTRS)
Bacon, Barton J.; Ostroff, Aaron J.
2000-01-01
This paper presents an approach to on-line control design for aircraft that have suffered either actuator failure, missing effector surfaces, surface damage, or any combination. The approach is based on a modified version of nonlinear dynamic inversion. The approach does not require a model of the baseline vehicle (effectors at zero deflection), but does require feedback of accelerations and effector positions. Implementation issues are addressed and the method is demonstrated on an advanced tailless aircraft. An experimental simulation analysis tool is used to directly evaluate the nonlinear system's stability robustness.
NASA Technical Reports Server (NTRS)
Gherlone, Marco; Cerracchio, Priscilla; Mattone, Massimiliano; Di Sciuva, Marco; Tessler, Alexander
2011-01-01
A robust and efficient computational method for reconstructing the three-dimensional displacement field of truss, beam, and frame structures, using measured surface-strain data, is presented. Known as shape sensing , this inverse problem has important implications for real-time actuation and control of smart structures, and for monitoring of structural integrity. The present formulation, based on the inverse Finite Element Method (iFEM), uses a least-squares variational principle involving strain measures of Timoshenko theory for stretching, torsion, bending, and transverse shear. Two inverse-frame finite elements are derived using interdependent interpolations whose interior degrees-of-freedom are condensed out at the element level. In addition, relationships between the order of kinematic-element interpolations and the number of required strain gauges are established. As an example problem, a thin-walled, circular cross-section cantilevered beam subjected to harmonic excitations in the presence of structural damping is modeled using iFEM; where, to simulate strain-gauge values and to provide reference displacements, a high-fidelity MSC/NASTRAN shell finite element model is used. Examples of low and high-frequency dynamic motion are analyzed and the solution accuracy examined with respect to various levels of discretization and the number of strain gauges.
Function-valued adaptive dynamics and the calculus of variations.
Parvinen, Kalle; Dieckmann, Ulf; Heino, Mikko
2006-01-01
Adaptive dynamics has been widely used to study the evolution of scalar-valued, and occasionally vector-valued, strategies in ecologically realistic models. In many ecological situations, however, evolving strategies are best described as function-valued, and thus infinite-dimensional, traits. So far, such evolution has only been studied sporadically, mostly based on quantitative genetics models with limited ecological realism. In this article we show how to apply the calculus of variations to find evolutionarily singular strategies of function-valued adaptive dynamics: such a strategy has to satisfy Euler's equation with environmental feedback. We also demonstrate how second-order derivatives can be used to investigate whether or not a function-valued singular strategy is evolutionarily stable. We illustrate our approach by presenting several worked examples. PMID:16012801
Dynamic Load Balancing for Adaptive Meshes using Symmetric Broadcast Networks
NASA Technical Reports Server (NTRS)
Das, Sajal K.; Harvey, Daniel J.; Biswas, Rupak; Saini, Subhash (Technical Monitor)
1998-01-01
Many scientific applications involve grids that lack a uniform underlying structure. These applications are often dynamic in the sense that the grid structure significantly changes between successive phases of execution. In parallel computing environments, mesh adaptation of grids through selective refinement/coarsening has proven to be an effective approach. However, achieving load balance while minimizing inter-processor communication and redistribution costs is a difficult problem. Traditional dynamic load balancers are mostly inadequate because they lack a global view across processors. In this paper, we compare a novel load balancer that utilizes symmetric broadcast networks (SBN) to a successful global load balancing environment (PLUM) created to handle adaptive unstructured applications. Our experimental results on the IBM SP2 demonstrate that performance of the proposed SBN load balancer is comparable to results achieved under PLUM.
Adaptive neural information processing with dynamical electrical synapses
Xiao, Lei; Zhang, Dan-ke; Li, Yuan-qing; Liang, Pei-ji; Wu, Si
2013-01-01
The present study investigates a potential computational role of dynamical electrical synapses in neural information process. Compared with chemical synapses, electrical synapses are more efficient in modulating the concerted activity of neurons. Based on the experimental data, we propose a phenomenological model for short-term facilitation of electrical synapses. The model satisfactorily reproduces the phenomenon that the neuronal correlation increases although the neuronal firing rates attenuate during the luminance adaptation. We explore how the stimulus information is encoded in parallel by firing rates and correlated activity of neurons, and find that dynamical electrical synapses mediate a transition from the firing rate code to the correlation one during the luminance adaptation. The latter encodes the stimulus information by using the concerted, but lower neuronal firing rate, and hence is economically more efficient. PMID:23596413
Effects of adaptive dynamical linking in networked games.
Yang, Zhihu; Li, Zhi; Wu, Te; Wang, Long
2013-10-01
The role of dynamical topologies in the evolution of cooperation has received considerable attention, as some studies have demonstrated that dynamical networks are much better than static networks in terms of boosting cooperation. Here we study a dynamical model of evolution of cooperation on stochastic dynamical networks in which there are no permanent partners to each agent. Whenever a new link is created, its duration is randomly assigned without any bias or preference. We allow the agent to adaptively adjust the duration of each link during the evolution in accordance with the feedback from game interactions. By Monte Carlo simulations, we find that cooperation can be remarkably promoted by this adaptive dynamical linking mechanism both for the game of pairwise interactions, such as the Prisoner's Dilemma game (PDG), and for the game of group interactions, illustrated by the public goods game (PGG). And the faster the adjusting rate, the more successful the evolution of cooperation. We also show that in this context weak selection favors cooperation much more than strong selection does. What is particularly meaningful is that the prosperity of cooperation in this study indicates that the rationality and selfishness of a single agent in adjusting social ties can lead to the progress of altruism of the whole population. PMID:24229137
Effects of adaptive dynamical linking in networked games
NASA Astrophysics Data System (ADS)
Yang, Zhihu; Li, Zhi; Wu, Te; Wang, Long
2013-10-01
The role of dynamical topologies in the evolution of cooperation has received considerable attention, as some studies have demonstrated that dynamical networks are much better than static networks in terms of boosting cooperation. Here we study a dynamical model of evolution of cooperation on stochastic dynamical networks in which there are no permanent partners to each agent. Whenever a new link is created, its duration is randomly assigned without any bias or preference. We allow the agent to adaptively adjust the duration of each link during the evolution in accordance with the feedback from game interactions. By Monte Carlo simulations, we find that cooperation can be remarkably promoted by this adaptive dynamical linking mechanism both for the game of pairwise interactions, such as the Prisoner's Dilemma game (PDG), and for the game of group interactions, illustrated by the public goods game (PGG). And the faster the adjusting rate, the more successful the evolution of cooperation. We also show that in this context weak selection favors cooperation much more than strong selection does. What is particularly meaningful is that the prosperity of cooperation in this study indicates that the rationality and selfishness of a single agent in adjusting social ties can lead to the progress of altruism of the whole population.
NASA Astrophysics Data System (ADS)
Grayver, Alexander V.; Kuvshinov, Alexey V.
2016-05-01
This paper presents a methodology to sample equivalence domain (ED) in nonlinear partial differential equation (PDE)-constrained inverse problems. For this purpose, we first applied state-of-the-art stochastic optimization algorithm called Covariance Matrix Adaptation Evolution Strategy (CMAES) to identify low-misfit regions of the model space. These regions were then randomly sampled to create an ensemble of equivalent models and quantify uncertainty. CMAES is aimed at exploring model space globally and is robust on very ill-conditioned problems. We show that the number of iterations required to converge grows at a moderate rate with respect to number of unknowns and the algorithm is embarrassingly parallel. We formulated the problem by using the generalized Gaussian distribution. This enabled us to seamlessly use arbitrary norms for residual and regularization terms. We show that various regularization norms facilitate studying different classes of equivalent solutions. We further show how performance of the standard Metropolis-Hastings Markov chain Monte Carlo algorithm can be substantially improved by using information CMAES provides. This methodology was tested by using individual and joint inversions of magneotelluric, controlled-source electromagnetic (EM) and global EM induction data.
NASA Astrophysics Data System (ADS)
Grayver, Alexander V.; Kuvshinov, Alexey V.
2016-02-01
This paper presents a methodology to sample equivalence domain (ED) in non-linear PDE-constrained inverse problems. For this purpose, we first applied state-of-the-art stochastic optimization algorithm called Covariance Matrix Adaptation Evolution Strategy (CMAES) to identify low misfit regions of the model space. These regions were then randomly sampled to create an ensemble of equivalent models and quantify uncertainty. CMAES is aimed at exploring model space globally and is robust on very ill-conditioned problems. We show that the number of iterations required to converge grows at a moderate rate with respect to number of unknowns and the algorithm is embarrassingly parallel. We formulated the problem by using the generalized Gaussian distribution. This enabled us to seamlessly use arbitrary norms for residual and regularization terms. We show that various regularization norms facilitate studying different classes of equivalent solutions. We further show how performance of the standard Metropolis-Hastings Markov chain Monte Carlo (MCMC) algorithm can be substantially improved by using information CMAES provides. This methodology was tested by using individual and joint inversions of Magneotelluric, Controlled-source Electromagnetic (EM) and Global EM induction data.
LASERS: Emission dynamics of coupled Nd3+ : YAG lasers with a shared population inversion source
NASA Astrophysics Data System (ADS)
Kaptsov, L. N.; Yatskiv, A. M.
1995-08-01
A calculation is made of the frequency spectrum of relaxation oscillations of a cw Nd3+ : YAG laser with several lasing channels intersecting in the same active element. It is shown that the highest frequency of relaxation oscillations of isolated channels is virtually retained when these channels are coupled by population inversion, but the other frequencies are reduced. The results of these calculations are in agreement with measurements carried out on a double-beam laser. Experiments are reported on the transition of such a system to dynamic chaos. As in the case of a solid-state laser which emits multimode (in respect of the longitudinal index) radiation, near the frequencies of relaxation oscillations the transition of the investigated system to dynamic chaos follows the Ruelle—Takens—Newhouse scenario.
The Dynamics of Tachyon Field with AN Inverse Square Potential in Loop Quantum Cosmology
NASA Astrophysics Data System (ADS)
Huang, Fei; Zhu, Jian-Yang; Xiao, Kui
2013-05-01
The dynamical behavior of tachyon field with an inverse potential is investigated in loop quantum cosmology. It reveals that the late-time behavior of tachyon field with this potential leads to a power-law expansion. In addition, an additional barotropic perfect fluid with the adiabatic index 0 < γ < 2 is added and the dynamical system is shown to be an autonomous one. The stability of this autonomous system is discussed using phase plane analysis. There exist up to five fixed points with only two of them possibly stable. The two stable node (attractor) solutions are specified and their cosmological indications are discussed. For the tachyon dominated solution, the further discussion is stretched to the possibility of considering tachyon field as a combination of two parts which respectively behave like dark matter and dark energy.
Mates, Steven P; Forster, Aaron M; Hunston, Donald; Rhorer, Richard; Everett, Richard K; Simmonds, Kirth E; Bagchi, Amit
2012-10-01
Soft elastomeric materials that mimic real soft human tissues are sought to provide realistic experimental devices to simulate the human body's response to blast loading to aid the development of more effective protective equipment. The dynamic mechanical behavior of these materials is often measured using a Kolsky bar because it can achieve both the high strain rates (>100s(-1)) and the large strains (>20%) that prevail in blast scenarios. Obtaining valid results is challenging, however, due to poor dynamic equilibrium, friction, and inertial effects. To avoid these difficulties, an inverse method was employed to determine the dynamic response of a soft, prospective biomimetic elastomer using Kolsky bar tests coupled with high-speed 3D digital image correlation. Individual tests were modeled using finite elements, and the dynamic stiffness of the elastomer was identified by matching the simulation results with test data using numerical optimization. Using this method, the average dynamic response was found to be nearly equivalent to the quasi-static response measured with stress-strain curves at compressive strains up to 60%, with an uncertainty of ±18%. Moreover, the behavior was consistent with the results in stress relaxation experiments and oscillatory tests although the latter were performed at lower strain levels. PMID:22982958
Yamasaki, Taiga; Idehara, Katsutoshi; Xin, Xin
2016-07-01
We propose a new method to estimate muscle activity in a straightforward manner with high accuracy and relatively small computational costs by using the external input of the joint angle and its first to fourth derivatives with respect to time. The method solves the inverse dynamics problem of the skeletal system, the forward dynamics problem of the muscular system, and the load-sharing problem of muscles as a static optimization of neural excitation signals. The external input including the higher-order derivatives is required for a calculation of constraints imposed on the load-sharing problem. The feasibility of the method is demonstrated by the simulation of a simple musculoskeletal model with a single joint. Moreover, the influences of the muscular dynamics, and the higher-order derivatives on the estimation of the muscle activity are demonstrated, showing the results when the time constants of the activation dynamics are very small, and the third and fourth derivatives of the external input are ignored, respectively. It is concluded that the method can have the potential to improve estimation accuracy of muscle activity of highly dynamic motions. PMID:27211782
Finnveden, Svante; Hörlin, Nils-Erik; Barbagallo, Mathias
2014-04-01
Viscoelastic properties of porous materials, typical of those used in vehicles for noise insulation and absorption, are estimated from measurements and inverse finite element procedures. The measurements are taken in a near vacuum and cover a broad frequency range: 20 Hz to 1 kHz. The almost cubic test samples were made of 25 mm foam covered by a "heavy layer" of rubber. They were mounted in a vacuum chamber on an aluminum table, which was excited in the vertical and horizontal directions with a shaker. Three kinds of response are measured allowing complete estimates of the viscoelastic moduli for isotropic materials and also providing some information on the degree of material anisotropicity. First, frequency independent properties are estimated, where dissipation is described by constant loss factors. Then, fractional derivative models that capture the variation with frequency of the stiffness and damping are adapted. The measurement setup is essentially two-dimensional and calculations are three-dimensional and for a state of plane strain. The good agreement between measured and calculated response provides some confidence in the presented procedures. If, however, the material model cannot fit the measurements well, the inverse procedure yields a certain degree of arbitrariness to the parameter estimation. PMID:25234982
NASA Astrophysics Data System (ADS)
Hong, M.; Zhang, R.; Li, J. X.; Ge, J. J.; Liu, K. F.
2013-02-01
Based on time series data of 500 hPa potential field from NCEP/NCAR (National Center for Environmental Forecast of American/National Center for Atmospheric Research), a novel consideration of empirical orthogonal function (EOF) time-space separation and dynamic system reconstruction for time series is introduced. This method consists of two parts: first, the dynamical model inversion and model parameter optimization are carried out on the EOF time coefficient series using the genetic algorithm (GA), and, second, a nonlinear dynamic model representing the subtropical high (SH) activity and its abnormality is established. The SH activity and its abnormal mechanism is studied using the developed dynamical model. Results show that the configuration and diversification of the SH equilibriums have good correspondence with the actual short-medium term abnormal activity of the SH. Change of SH potential field brought by the combination of equilibriums is more complex than that by mutation, and their exhibition patterns are different. The mutation behavior from high-value to low-value equilibriums of the SH in summer corresponds with the southward drop of the SH in the observed weather process. The combination behavior of the two steady equilibriums corresponds with disappearance of the "double-ridge" phenomenon of the SH. Dynamical mechanisms of these phenomena are explained.
Glassy Dynamics in the Adaptive Immune Response Prevents Autoimmune Disease
NASA Astrophysics Data System (ADS)
Sun, Jun; Earl, David J.; Deem, Michael W.
2005-09-01
The immune system normally protects the human host against death by infection. However, when an immune response is mistakenly directed at self-antigens, autoimmune disease can occur. We describe a model of protein evolution to simulate the dynamics of the adaptive immune response to antigens. Computer simulations of the dynamics of antibody evolution show that different evolutionary mechanisms, namely, gene segment swapping and point mutation, lead to different evolved antibody binding affinities. Although a combination of gene segment swapping and point mutation can yield a greater affinity to a specific antigen than point mutation alone, the antibodies so evolved are highly cross reactive and would cause autoimmune disease, and this is not the chosen dynamics of the immune system. We suggest that in the immune system’s search for antibodies, a balance has evolved between binding affinity and specificity.
Adaptive dynamic programming for auto-resilient video streaming
NASA Astrophysics Data System (ADS)
Zhao, Juan; Li, Xingmei; Wang, Wei; Wu, Guoping
2007-11-01
Wireless video transmission encounters higher error rate than in wired network, which introduces distortion into the error-sensitive compressed data, reducing the quality of the playback video. Therefore, to ensure the end-to-end quality, wireless video needs a transmission system including both efficient source coding scheme and transmission technology against the influence of the channel error. This paper tackles a dynamic programming algorithm for robust video streaming over error-prone channels. An auto-resilient multiple-description coding with optimized transmission strategy has been proposed. Further study is done on the computational complexity of rate-distortion optimized video streaming and a dynamic programming algorithm is considered. Experiment results show that video streaming with adaptive dynamic programming gains better playback video quality at the receiver when transmitted through error-prone mobile channel.
Generalization in Adaptation to Stable and Unstable Dynamics
Kadiallah, Abdelhamid; Franklin, David W.; Burdet, Etienne
2012-01-01
Humans skillfully manipulate objects and tools despite the inherent instability. In order to succeed at these tasks, the sensorimotor control system must build an internal representation of both the force and mechanical impedance. As it is not practical to either learn or store motor commands for every possible future action, the sensorimotor control system generalizes a control strategy for a range of movements based on learning performed over a set of movements. Here, we introduce a computational model for this learning and generalization, which specifies how to learn feedforward muscle activity in a function of the state space. Specifically, by incorporating co-activation as a function of error into the feedback command, we are able to derive an algorithm from a gradient descent minimization of motion error and effort, subject to maintaining a stability margin. This algorithm can be used to learn to coordinate any of a variety of motor primitives such as force fields, muscle synergies, physical models or artificial neural networks. This model for human learning and generalization is able to adapt to both stable and unstable dynamics, and provides a controller for generating efficient adaptive motor behavior in robots. Simulation results exhibit predictions consistent with all experiments on learning of novel dynamics requiring adaptation of force and impedance, and enable us to re-examine some of the previous interpretations of experiments on generalization. PMID:23056191
Generalization in adaptation to stable and unstable dynamics.
Kadiallah, Abdelhamid; Franklin, David W; Burdet, Etienne
2012-01-01
Humans skillfully manipulate objects and tools despite the inherent instability. In order to succeed at these tasks, the sensorimotor control system must build an internal representation of both the force and mechanical impedance. As it is not practical to either learn or store motor commands for every possible future action, the sensorimotor control system generalizes a control strategy for a range of movements based on learning performed over a set of movements. Here, we introduce a computational model for this learning and generalization, which specifies how to learn feedforward muscle activity in a function of the state space. Specifically, by incorporating co-activation as a function of error into the feedback command, we are able to derive an algorithm from a gradient descent minimization of motion error and effort, subject to maintaining a stability margin. This algorithm can be used to learn to coordinate any of a variety of motor primitives such as force fields, muscle synergies, physical models or artificial neural networks. This model for human learning and generalization is able to adapt to both stable and unstable dynamics, and provides a controller for generating efficient adaptive motor behavior in robots. Simulation results exhibit predictions consistent with all experiments on learning of novel dynamics requiring adaptation of force and impedance, and enable us to re-examine some of the previous interpretations of experiments on generalization. PMID:23056191
NASA Technical Reports Server (NTRS)
Lee, C. S. G.; Chen, C. L.
1989-01-01
Two efficient mapping algorithms for scheduling the robot inverse dynamics computation consisting of m computational modules with precedence relationship to be executed on a multiprocessor system consisting of p identical homogeneous processors with processor and communication costs to achieve minimum computation time are presented. An objective function is defined in terms of the sum of the processor finishing time and the interprocessor communication time. The minimax optimization is performed on the objective function to obtain the best mapping. This mapping problem can be formulated as a combination of the graph partitioning and the scheduling problems; both have been known to be NP-complete. Thus, to speed up the searching for a solution, two heuristic algorithms were proposed to obtain fast but suboptimal mapping solutions. The first algorithm utilizes the level and the communication intensity of the task modules to construct an ordered priority list of ready modules and the module assignment is performed by a weighted bipartite matching algorithm. For a near-optimal mapping solution, the problem can be solved by the heuristic algorithm with simulated annealing. These proposed optimization algorithms can solve various large-scale problems within a reasonable time. Computer simulations were performed to evaluate and verify the performance and the validity of the proposed mapping algorithms. Finally, experiments for computing the inverse dynamics of a six-jointed PUMA-like manipulator based on the Newton-Euler dynamic equations were implemented on an NCUBE/ten hypercube computer to verify the proposed mapping algorithms. Computer simulation and experimental results are compared and discussed.
Dynamic data-driven sensor network adaptation for border control
NASA Astrophysics Data System (ADS)
Bein, Doina; Madan, Bharat B.; Phoha, Shashi; Rajtmajer, Sarah; Rish, Anna
2013-06-01
Given a specific scenario for the border control problem, we propose a dynamic data-driven adaptation of the associated sensor network via embedded software agents which make sensor network control, adaptation and collaboration decisions based on the contextual information value of competing data provided by different multi-modal sensors. We further propose the use of influence diagrams to guide data-driven decision making in selecting the appropriate action or course of actions which maximize a given utility function by designing a sensor embedded software agent that uses an influence diagram to make decisions about whether to engage or not engage higher level sensors for accurately detecting human presence in the region. The overarching goal of the sensor system is to increase the probability of target detection and classification and reduce the rate of false alarms. The proposed decision support software agent is validated experimentally on a laboratory testbed for multiple border control scenarios.
Elucidating Microbial Adaptation Dynamics via Autonomous Exposure and Sampling
NASA Astrophysics Data System (ADS)
Grace, J. M.; Verseux, C.; Gentry, D.; Moffet, A.; Thayabaran, R.; Wong, N.; Rothschild, L.
2013-12-01
The adaptation of micro-organisms to their environments is a complex process of interaction between the pressures of the environment and of competition. Reducing this multifactorial process to environmental exposure in the laboratory is a common tool for elucidating individual mechanisms of evolution, such as mutation rates[Wielgoss et al., 2013]. Although such studies inform fundamental questions about the way adaptation and even speciation occur, they are often limited by labor-intensive manual techniques[Wassmann et al., 2010]. Current methods for controlled study of microbial adaptation limit the length of time, the depth of collected data, and the breadth of applied environmental conditions. Small idiosyncrasies in manual techniques can have large effects on outcomes; for example, there are significant variations in induced radiation resistances following similar repeated exposure protocols[Alcántara-Díaz et al., 2004; Goldman and Travisano, 2011]. We describe here a project under development to allow rapid cycling of multiple types of microbial environmental exposure. The system allows continuous autonomous monitoring and data collection of both single species and sampled communities, independently and concurrently providing multiple types of controlled environmental pressure (temperature, radiation, chemical presence or absence, and so on) to a microbial community in dynamic response to the ecosystem's current status. When combined with DNA sequencing and extraction, such a controlled environment can cast light on microbial functional development, population dynamics, inter- and intra-species competition, and microbe-environment interaction. The project's goal is to allow rapid, repeatable iteration of studies of both natural and artificial microbial adaptation. As an example, the same system can be used both to increase the pH of a wet soil aliquot over time while periodically sampling it for genetic activity analysis, or to repeatedly expose a culture of
Development of a dynamically adaptive grid method for multidimensional problems
NASA Astrophysics Data System (ADS)
Holcomb, J. E.; Hindman, R. G.
1984-06-01
An approach to solution adaptive grid generation for use with finite difference techniques, previously demonstrated on model problems in one space dimension, has been extended to multidimensional problems. The method is based on the popular elliptic steady grid generators, but is 'dynamically' adaptive in the sense that a grid is maintained at all times satisfying the steady grid law driven by a solution-dependent source term. Testing has been carried out on Burgers' equation in one and two space dimensions. Results appear encouraging both for inviscid wave propagation cases and viscous boundary layer cases, suggesting that application to practical flow problems is now possible. In the course of the work, obstacles relating to grid correction, smoothing of the solution, and elliptic equation solvers have been largely overcome. Concern remains, however, about grid skewness, boundary layer resolution and the need for implicit integration methods. Also, the method in 3-D is expected to be very demanding of computer resources.
Direct Adaptive Aircraft Control Using Dynamic Cell Structure Neural Networks
NASA Technical Reports Server (NTRS)
Jorgensen, Charles C.
1997-01-01
A Dynamic Cell Structure (DCS) Neural Network was developed which learns topology representing networks (TRNS) of F-15 aircraft aerodynamic stability and control derivatives. The network is integrated into a direct adaptive tracking controller. The combination produces a robust adaptive architecture capable of handling multiple accident and off- nominal flight scenarios. This paper describes the DCS network and modifications to the parameter estimation procedure. The work represents one step towards an integrated real-time reconfiguration control architecture for rapid prototyping of new aircraft designs. Performance was evaluated using three off-line benchmarks and on-line nonlinear Virtual Reality simulation. Flight control was evaluated under scenarios including differential stabilator lock, soft sensor failure, control and stability derivative variations, and air turbulence.
NASA Astrophysics Data System (ADS)
Piacentino, Michael R.; Berends, David C.; Zhang, David C.; Gudis, Eduardo
2013-05-01
Two of the biggest challenges in designing U×V vision systems are properly representing high dynamic range scene content using low dynamic range components and reducing camera motion blur. SRI's MASI-HDR (Motion Adaptive Signal Integration-High Dynamic Range) is a novel technique for generating blur-reduced video using multiple captures for each displayed frame while increasing the effective camera dynamic range by four bits or more. MASI-HDR processing thus provides high performance video from rapidly moving platforms in real-world conditions in low latency real time, enabling even the most demanding applications on air, ground and water.
An error function minimization approach for the inverse problem of adaptive mirrors tuning
NASA Astrophysics Data System (ADS)
Vannoni, Maurizio; Yang, Fan; Siewert, Frank; Sinn, Harald
2014-09-01
Adaptive x-ray optics are more and more used in synchrotron beamlines, and it is probable that they will be considered for the future high-power free-electron laser sources, as the European XFEL now under construction in Hamburg, or similar projects now in discussion. These facilities will deliver a high power x-ray beam, with an expected high heat load delivered on the optics. For this reason, bendable mirrors are required to actively compensate the resulting wavefront distortion. On top of that, the mirror could have also intrinsic surface defects, as polishing errors or mounting stresses. In order to be able to correct the mirror surface with a high precision to maintain its challenging requirements, the mirror surface is usually characterized with a high accuracy metrology to calculate the actuators pulse functions and to assess its initial shape. After that, singular value decomposition (SVD) is used to find the signals to be applied into the actuators, to reach the desired surface deformation or correction. But in some cases this approach could be not robust enough for the needed performance. We present here a comparison between the classical SVD method and an error function minimization based on root-mean-square calculation. Some examples are provided, using a simulation of the European XFEL mirrors design as a case of study, and performances of the algorithms are evaluated in order to reach the ultimate quality in different scenarios. The approach could be easily generalized to other situations as well.
Adaptive and neuroadaptive control for nonnegative and compartmental dynamical systems
NASA Astrophysics Data System (ADS)
Volyanskyy, Kostyantyn Y.
Neural networks have been extensively used for adaptive system identification as well as adaptive and neuroadaptive control of highly uncertain systems. The goal of adaptive and neuroadaptive control is to achieve system performance without excessive reliance on system models. To improve robustness and the speed of adaptation of adaptive and neuroadaptive controllers several controller architectures have been proposed in the literature. In this dissertation, we develop a new neuroadaptive control architecture for nonlinear uncertain dynamical systems. The proposed framework involves a novel controller architecture with additional terms in the update laws that are constructed using a moving window of the integrated system uncertainty. These terms can be used to identify the ideal system weights of the neural network as well as effectively suppress system uncertainty. Linear and nonlinear parameterizations of the system uncertainty are considered and state and output feedback neuroadaptive controllers are developed. Furthermore, we extend the developed framework to discrete-time dynamical systems. To illustrate the efficacy of the proposed approach we apply our results to an aircraft model with wing rock dynamics, a spacecraft model with unknown moment of inertia, and an unmanned combat aerial vehicle undergoing actuator failures, and compare our results with standard neuroadaptive control methods. Nonnegative systems are essential in capturing the behavior of a wide range of dynamical systems involving dynamic states whose values are nonnegative. A sub-class of nonnegative dynamical systems are compartmental systems. These systems are derived from mass and energy balance considerations and are comprised of homogeneous interconnected microscopic subsystems or compartments which exchange variable quantities of material via intercompartmental flow laws. In this dissertation, we develop direct adaptive and neuroadaptive control framework for stabilization, disturbance
Adaptive mesh refinement for 1-dimensional gas dynamics
Hedstrom, G.; Rodrigue, G.; Berger, M.; Oliger, J.
1982-01-01
We consider the solution of the one-dimensional equation of gas-dynamics. Accurate numerical solutions are difficult to obtain on a given spatial mesh because of the existence of physical regions where components of the exact solution are either discontinuous or have large gradient changes. Numerical methods treat these phenomena in a variety of ways. In this paper, the method of adaptive mesh refinement is used. A thorough description of this method for general hyperbolic systems is given elsewhere and only properties of the method pertinent to the system are elaborated.
Neural Network Assisted Inverse Dynamic Guidance for Terminally Constrained Entry Flight
Chen, Wanchun
2014-01-01
This paper presents a neural network assisted entry guidance law that is designed by applying Bézier approximation. It is shown that a fully constrained approximation of a reference trajectory can be made by using the Bézier curve. Applying this approximation, an inverse dynamic system for an entry flight is solved to generate guidance command. The guidance solution thus gotten ensures terminal constraints for position, flight path, and azimuth angle. In order to ensure terminal velocity constraint, a prediction of the terminal velocity is required, based on which, the approximated Bézier curve is adjusted. An artificial neural network is used for this prediction of the terminal velocity. The method enables faster implementation in achieving fully constrained entry flight. Results from simulations indicate improved performance of the neural network assisted method. The scheme is expected to have prospect for further research on automated onboard control of terminal velocity for both reentry and terminal guidance laws. PMID:24723821
Neural network assisted inverse dynamic guidance for terminally constrained entry flight.
Zhou, Hao; Rahman, Tawfiqur; Chen, Wanchun
2014-01-01
This paper presents a neural network assisted entry guidance law that is designed by applying Bézier approximation. It is shown that a fully constrained approximation of a reference trajectory can be made by using the Bézier curve. Applying this approximation, an inverse dynamic system for an entry flight is solved to generate guidance command. The guidance solution thus gotten ensures terminal constraints for position, flight path, and azimuth angle. In order to ensure terminal velocity constraint, a prediction of the terminal velocity is required, based on which, the approximated Bézier curve is adjusted. An artificial neural network is used for this prediction of the terminal velocity. The method enables faster implementation in achieving fully constrained entry flight. Results from simulations indicate improved performance of the neural network assisted method. The scheme is expected to have prospect for further research on automated onboard control of terminal velocity for both reentry and terminal guidance laws. PMID:24723821
Adaptation to shape switching by component selection in a constitutional dynamic system.
Ulrich, Sébastien; Lehn, Jean-Marie
2009-04-22
Molecules having different accessible shape states, which can be addressed in an effector-controlled manner, may be termed morphological switches. A dynamic covalent system can undergo adaptation to each state of a two-state morphological switch by generation of an optimal constitution through component selection. We have studied such a component selection in the dynamic covalent constituents generated by metal cation-induced shape switching of a core component between two states of W and U shape, characterized by both different geometries and different coordination features. The system performs shape-dependent self-sorting of metal ions and components. The origin of the selectivity was investigated through competition experiments, in solution and by analysis of solid state structures, which reveal the role of the molecular shape in the formation of a particular self-assembled architecture. The coordination features of each state as well as phase change also play an important role, in addition to the shape plasticity, in steering the covalent dynamic system toward the formation of a given entity by the selection of the most appropriate components. Different examples are described which show that the morphological switching of one component of a given self-assembled entity can lead to the exchange of the complementary one, which is no longer the best partner, for a new partner, able to form a more stable new assembly. Thus, the constitutional evolution of these dynamic systems is steered by the shape of a given state via both its geometry and its coordination features toward metal ions, leading to incorporation/decorporation of the most appropriate components. The controlled interconversion of the shape states of the morphological switches, induced by addition/removal of metal ions, results in a constitutional adaptation behavior through inversion of the selection preferences. PMID:19206535
Inter-limb interference during bimanual adaptation to dynamic environments
Casadio, Maura; Sanguineti, Vittorio; Squeri, Valentina; Masia, Lorenzo; Morasso, Pietro
2015-01-01
Skillful manipulation of objects often requires the spatio-temporal coordination of both hands and, at the same time, the compensation of environmental forces. In bimanual coordination, movements of the two hands may be coupled because each hand needs to compensate the forces generated by the other hand or by an object operated by both hands (dynamic coupling), or because the two hands share the same workspace (spatial coupling). We examined how spatial coupling influences bimanual coordination, by looking at the adaptation of velocity-dependent force fields during a task in which the two hands simultaneously perform center-out reaching movements with the same initial position and the same targets, equally spaced on a circle. Subjects were randomly allocated to two groups, which differed in terms of the force fields they were exposed to: in one group (CW-CW), force fields had equal clockwise orientations in both hands; in the other group (CCW-CW), they had opposite orientations. In both groups, in randomly selected trials (catch trials) of the adaptation phase, the force fields were unexpectedly removed. Adaptation was quantified in terms of the changes of directional error for both hand trajectories. Bimanual coordination was quantified in terms of inter-limb longitudinal and sideways displacements, in force field and in catch trials. Experimental results indicate that both arms could simultaneously adapt to the two force fields. However, in the CCW-CW group, adaptation was incomplete for the movements from the central position to the more distant targets with respect to the body. In addition, in this group the left hand systematically leads in the movements toward targets on the left of the starting position, whereas the right hand leads in the movements to targets on the right. We show that these effects are due to a gradual sideways shift of the hands, so that during movements the left hand tends to consistently remain at the left of the right hand. These
Hydrodynamics in adaptive resolution particle simulations: Multiparticle collision dynamics
NASA Astrophysics Data System (ADS)
Alekseeva, Uliana; Winkler, Roland G.; Sutmann, Godehard
2016-06-01
A new adaptive resolution technique for particle-based multi-level simulations of fluids is presented. In the approach, the representation of fluid and solvent particles is changed on the fly between an atomistic and a coarse-grained description. The present approach is based on a hybrid coupling of the multiparticle collision dynamics (MPC) method and molecular dynamics (MD), thereby coupling stochastic and deterministic particle-based methods. Hydrodynamics is examined by calculating velocity and current correlation functions for various mixed and coupled systems. We demonstrate that hydrodynamic properties of the mixed fluid are conserved by a suitable coupling of the two particle methods, and that the simulation results agree well with theoretical expectations.
Experimental Dynamic Characterization of a Reconfigurable Adaptive Precision Truss
NASA Technical Reports Server (NTRS)
Hinkle, J. D.; Peterson, L. D.
1994-01-01
The dynamic behavior of a reconfigurable adaptive truss structure with non-linear joints is investigated. The objective is to experimentally examine the effects of the local non-linearities on the global dynamics of the structure. Amplitude changes in the frequency response functions are measured at micron levels of motion. The amplitude and frequency variations of a number of modes indicate a non-linear Coulomb friction response. Hysteretic bifurcation behavior is also measured at an amplitude approximately equal to the specified free-play in the joint. Under the 1 g pre-load, however, the non-linearity was dominantly characteristic of Coulomb friction with little evidence of free-play stiffening.
Adaptive integral dynamic surface control of a hypersonic flight vehicle
NASA Astrophysics Data System (ADS)
Aslam Butt, Waseem; Yan, Lin; Amezquita S., Kendrick
2015-07-01
In this article, non-linear adaptive dynamic surface air speed and flight path angle control designs are presented for the longitudinal dynamics of a flexible hypersonic flight vehicle. The tracking performance of the control design is enhanced by introducing a novel integral term that caters to avoiding a large initial control signal. To ensure feasibility, the design scheme incorporates magnitude and rate constraints on the actuator commands. The uncertain non-linear functions are approximated by an efficient use of the neural networks to reduce the computational load. A detailed stability analysis shows that all closed-loop signals are uniformly ultimately bounded and the ? tracking performance is guaranteed. The robustness of the design scheme is verified through numerical simulations of the flexible flight vehicle model.
NASA Astrophysics Data System (ADS)
Shubitidze, Fridon; Miller, Jonathan S.; Schultz, Gregory M.; Marble, Jay A.
2010-04-01
This paper reports vehicle based electromagnetic induction (EMI) array sensor data inversion and discrimination results. Recent field studies show that EMI arrays, such as the Minelab Single Transmitter Multiple Receiver (STMR), and the Geophex GEM-5 EMI array, provide a fast and safe way to detect subsurface metallic targets such as landmines, unexploded ordnance (UXO) and buried explosives. The array sensors are flexible and easily adaptable for a variety of ground vehicles and mobile platforms, which makes them very attractive for safe and cost effective detection operations in many applications, including but not limited to explosive ordnance disposal and humanitarian UXO and demining missions. Most state-of-the-art EMI arrays measure the vertical or full vector field, or gradient tensor fields and utilize them for real-time threat detection based on threshold analysis. Real field practice shows that the threshold-level detection has high false alarms. One way to reduce these false alarms is to use EMI numerical techniques that are capable of inverting EMI array data in real time. In this work a physically complete model, known as the normalized volume/surface magnetic sources (NV/SMS) model is adapted to the vehicle-based EMI array, such as STMR and GEM-5, data. The NV/SMS model can be considered as a generalized volume or surface dipole model, which in a special limited case coincides with an infinitesimal dipole model approach. According to the NV/SMS model, an object's response to a sensor's primary field is modeled mathematically by a set of equivalent magnetic dipoles, distributed inside the object (i.e. NVMS) or over a surface surrounding the object (i.e. NSMS). The scattered magnetic field of the NSMS is identical to that produced by a set of interacting magnetic dipoles. The amplitudes of the magnetic dipoles are normalized to the primary magnetic field, relating induced magnetic dipole polarizability and the primary magnetic field. The magnitudes of
Cunha, J Adam; Hsu, I-Chow; Pouliot, Jean; Roach Iii, Mack; Shinohara, Katsuto; Kurhanewicz, John; Reed, Galen; Stoianovici, Dan
2010-08-01
To translate any robot into a clinical environment, it is critical that the robot can seamlessly integrate with all the technology of a modern clinic. MRBot, an MR-stealth brachytherapy delivery device, was used in a closed-bore 3T MRI and a clinical brachytherapy cone beam CT suite. Targets included ceramic dummy seeds, MR-Spectroscopy-sensitive metabolite, and a prostate phantom. Acquired DICOM images were exported to planning software to register the robot coordinates in the imager's frame, contour and verify target locations, create dose plans, and export needle and seed positions to the robot. The coordination of each system element (imaging device, brachytherapy planning system, robot control, robot) was validated with a seed delivery accuracy of within 2 mm in both a phantom and soft tissue. An adaptive workflow was demonstrated by acquiring images after needle insertion and prior to seed deposition. This allows for adjustment if the needle is in the wrong position. Inverse planning (IPSA) was used to generate a seed placement plan and coordinates for ten needles and 29 seeds were transferred to the robot. After every two needles placed, an image was acquired. The placed seeds were identified and validated prior to placing the seeds in the next two needles. The ability to robotically deliver seeds to locations determined by IPSA and the ability of the system to incorporate novel needle patterns were demonstrated. Shown here is the ability to overcome this critical step. An adaptive brachytherapy workflow is demonstrated which integrates a clinical anatomy-based seed location optimization engine and a robotic brachytherapy device. Demonstration of this workflow is a key element of a successful translation to the clinic of the MRI stealth robotic delivery system, MRBot. PMID:20642386
CUNHA, J. ADAM; HSU, I-CHOW; POULIOT, JEAN; ROACH, MACK; SHINOHARA, KATSUTO; KURHANEWICZ, JOHN; REED, GALEN; STOIANOVICI, DAN
2011-01-01
To translate any robot into a clinical environment, it is critical that the robot can seamlessly integrate with all the technology of a modern clinic. MRBot, an MR-stealth brachytherapy delivery device, was used in a closed-bore 3T MRI and a clinical brachytherapy cone beam CT suite. Targets included ceramic dummy seeds, MR-Spectroscopy-sensitive metabolite, and a prostate phantom. Acquired DICOM images were exported to planning software to register the robot coordinates in the imager’s frame, contour and verify target locations, create dose plans, and export needle and seed positions to the robot. The coordination of each system element (imaging device, brachytherapy planning system, robot control, robot) was validated with a seed delivery accuracy of within 2 mm in both a phantom and soft tissue. An adaptive workflow was demonstrated by acquiring images after needle insertion and prior to seed deposition. This allows for adjustment if the needle is in the wrong position. Inverse planning (IPSA) was used to generate a seed placement plan and coordinates for ten needles and 29 seeds were transferred to the robot. After every two needles placed, an image was acquired. The placed seeds were identified and validated prior to placing the seeds in the next two needles. The ability to robotically deliver seeds to locations determined by IPSA and the ability of the system to incorporate novel needle patterns were demonstrated. Shown here is the ability to overcome this critical step. An adaptive brachytherapy workflow is demonstrated which integrates a clinical anatomy-based seed location optimization engine and a robotic brachytherapy device. Demonstration of this workflow is a key element of a successful translation to the clinic of the MRI stealth robotic delivery system, MRBot. PMID:20642386
Ultra-Low Power Dynamic Knob in Adaptive Compressed Sensing Towards Biosignal Dynamics.
Wang, Aosen; Lin, Feng; Jin, Zhanpeng; Xu, Wenyao
2016-06-01
Compressed sensing (CS) is an emerging sampling paradigm in data acquisition. Its integrated analog-to-information structure can perform simultaneous data sensing and compression with low-complexity hardware. To date, most of the existing CS implementations have a fixed architectural setup, which lacks flexibility and adaptivity for efficient dynamic data sensing. In this paper, we propose a dynamic knob (DK) design to effectively reconfigure the CS architecture by recognizing the biosignals. Specifically, the dynamic knob design is a template-based structure that comprises a supervised learning module and a look-up table module. We model the DK performance in a closed analytic form and optimize the design via a dynamic programming formulation. We present the design on a 130 nm process, with a 0.058 mm (2) fingerprint and a 187.88 nJ/event energy-consumption. Furthermore, we benchmark the design performance using a publicly available dataset. Given the energy constraint in wireless sensing, the adaptive CS architecture can consistently improve the signal reconstruction quality by more than 70%, compared with the traditional CS. The experimental results indicate that the ultra-low power dynamic knob can provide an effective adaptivity and improve the signal quality in compressed sensing towards biosignal dynamics. PMID:26800548
Structural dynamic health monitoring of adaptive CFRP structures
NASA Astrophysics Data System (ADS)
Kaiser, Stephan; Melcher, Joerg; Breitbach, Elmar J.; Sachau, Delf
1999-07-01
The DLR Institute of Structural Mechanics is engaged in the construction and optimization of adaptive structures for aerospace and terrestrial applications. Due to the FFS- Project, one of the recent works of the Institute is the reduction of buffet induced vibration loads at a fin. The construction of modern aircrafts is influenced b the increasing use of fiber composites. They have more specific stiffness and strength properties than metals. On the other hand the layered structure leads to new kinds of damages like delaminations. In the fin interface there are actuators and sensors integrated. Therefore the fin is connected with a controller. For the extension of this adaptive system towards an on-line tool for health monitoring this controller can be used as an identifier of the structure's modal parameters. The most promising procedure is based on MX filters. These filters constitute the filter coefficients from which a fast transformation procedure extracts the modal parameters. The changes of these parameters are related to the location and extent of the damage. So when using the already integrate controller for system identification, one can have a low-cost on-line damage detection for dynamic adaptive structures. First off-line test at CFRP plates have shown the ability to detect delaminations.
Optimal spectral tracking--adapting to dynamic regime change.
Brittain, John-Stuart; Halliday, David M
2011-01-30
Real world data do not always obey the statistical restraints imposed upon them by sophisticated analysis techniques. In spectral analysis for instance, an ergodic process--the interchangeability of temporal for spatial averaging--is assumed for a repeat-trial design. Many evolutionary scenarios, such as learning and motor consolidation, do not conform to such linear behaviour and should be approached from a more flexible perspective. To this end we previously introduced the method of optimal spectral tracking (OST) in the study of trial-varying parameters. In this extension to our work we modify the OST routines to provide an adaptive implementation capable of reacting to dynamic transitions in the underlying system state. In so doing, we generalise our approach to characterise both slow-varying and rapid fluctuations in time-series, simultaneously providing a metric of system stability. The approach is first applied to a surrogate dataset and compared to both our original non-adaptive solution and spectrogram approaches. The adaptive OST is seen to display fast convergence and desirable statistical properties. All three approaches are then applied to a neurophysiological recording obtained during a study on anaesthetic monitoring. Local field potentials acquired from the posterior hypothalamic region of a deep brain stimulation patient undergoing anaesthesia were analysed. The characterisation of features such as response delay, time-to-peak and modulation brevity are considered. PMID:21115043
Autonomous Path-Following by Approximate Inverse Dynamics and Vector Field Prediction
NASA Astrophysics Data System (ADS)
Gerlach, Adam R.
In this dissertation, we develop two general frameworks for the navigation and control of autonomous vehicles that must follow predefined paths. These frameworks are designed such that they inherently provide accurate navigation and control of a wide class of systems directly from a model of the vehicle's dynamics. The first framework introduced is the inverse dynamics by radial basis function (IDRBF) algorithm, which exploits the best approximation property of radial basis functions to accurately approximate the inverse dynamics of non-linear systems. This approximation is then used with the known, desired state of the system at a future time point to generate the system input that must be applied to reach the desired state in the specified time interval. The IDRBF algorithm is then tested on two non-linear dynamic systems, and accurate path-following is demonstrated. The second framework introduced is the predictive vector field (PVF) algorithm. The PVF algorithm uses the equations of motion and constraints of the system to predict a set of reachable states by sampling the system's configuration space. By finding and minimizing a continuous mapping between the system's configuration space and a cost space relating the reachable states of the system with a vector field (VF), one can determine the system inputs required to follow the VF. The PVF algorithm is then tested on the Dubin's vehicle and aircraft models, and accurate path-following is demonstrated. As the PVF algorithm's performance is dependent on the quality of the underlying system model and VF, algorithms are introduced for automatically generating VFs for constant altitude paths defined by a series of waypoints and for handling modeling uncertainties. Additionally, we provide a mathematical proof showing that this method can automatically produce VFs of the desired form. To handle modeling uncertainties, we enhance the PVF algorithm with the Gaussian process machine learning framework, enabling the
Adaptive dynamic programming as a theory of sensorimotor control.
Jiang, Yu; Jiang, Zhong-Ping
2014-08-01
Many characteristics of sensorimotor control can be explained by models based on optimization and optimal control theories. However, most of the previous models assume that the central nervous system has access to the precise knowledge of the sensorimotor system and its interacting environment. This viewpoint is difficult to be justified theoretically and has not been convincingly validated by experiments. To address this problem, this paper presents a new computational mechanism for sensorimotor control from a perspective of adaptive dynamic programming (ADP), which shares some features of reinforcement learning. The ADP-based model for sensorimotor control suggests that a command signal for the human movement is derived directly from the real-time sensory data, without the need to identify the system dynamics. An iterative learning scheme based on the proposed ADP theory is developed, along with rigorous convergence analysis. Interestingly, the computational model as advocated here is able to reproduce the motor learning behavior observed in experiments where a divergent force field or velocity-dependent force field was present. In addition, this modeling strategy provides a clear way to perform stability analysis of the overall system. Hence, we conjecture that human sensorimotor systems use an ADP-type mechanism to control movements and to achieve successful adaptation to uncertainties present in the environment. PMID:24962078
NASA Astrophysics Data System (ADS)
Cheung, Mark C. M.; Boerner, P.; Schrijver, C. J.; Testa, P.; Chen, F.; Peter, H.; Malanushenko, A.
2015-07-01
We present a new method for performing differential emission measure (DEM) inversions on narrow-band EUV images from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory. The method yields positive definite DEM solutions by solving a linear program. This method has been validated against a diverse set of thermal models of varying complexity and realism. These include (1) idealized Gaussian DEM distributions, (2) 3D models of NOAA Active Region 11158 comprising quasi-steady loop atmospheres in a nonlinear force-free field, and (3) thermodynamic models from a fully compressible, 3D MHD simulation of active region (AR) corona formation following magnetic flux emergence. We then present results from the application of the method to AIA observations of Active Region 11158, comparing the region's thermal structure on two successive solar rotations. Additionally, we show how the DEM inversion method can be adapted to simultaneously invert AIA and Hinode X-ray Telescope data, and how supplementing AIA data with the latter improves the inversion result. The speed of the method allows for routine production of DEM maps, thus facilitating science studies that require tracking of the thermal structure of the solar corona in time and space.
Analysis of dynamic deformation processes with adaptive KALMAN-filtering
NASA Astrophysics Data System (ADS)
Eichhorn, Andreas
2007-05-01
In this paper the approach of a full system analysis is shown quantifying a dynamic structural ("white-box"-) model for the calculation of thermal deformations of bar-shaped machine elements. The task was motivated from mechanical engineering searching new methods for the precise prediction and computational compensation of thermal influences in the heating and cooling phases of machine tools (i.e. robot arms, etc.). The quantification of thermal deformations under variable dynamic loads requires the modelling of the non-stationary spatial temperature distribution inside the object. Based upon FOURIERS law of heat flow the high-grade non-linear temperature gradient is represented by a system of partial differential equations within the framework of a dynamic Finite Element topology. It is shown that adaptive KALMAN-filtering is suitable to quantify relevant disturbance influences and to identify thermal parameters (i.e. thermal diffusivity) with a deviation of only 0,2%. As result an identified (and verified) parametric model for the realistic prediction respectively simulation of dynamic temperature processes is presented. Classifying the thermal bend as the main deformation quantity of bar-shaped machine tools, the temperature model is extended to a temperature deformation model. In lab tests thermal load steps are applied to an aluminum column. Independent control measurements show that the identified model can be used to predict the columns bend with a mean deviation (
Costa, L; Mantha, V R; Silva, A J; Fernandes, R J; Marinho, D A; Vilas-Boas, J P; Machado, L; Rouboa, A
2015-07-16
Computational fluid dynamics (CFD) plays an important role to quantify, understand and "observe" the water movements around the human body and its effects on drag (D). We aimed to investigate the flow effects around the swimmer and to compare the drag and drag coefficient (CD) values obtained from experiments (using cable velocimetry in a swimming pool) with those of CFD simulations for the two ventral gliding positions assumed during the breaststroke underwater cycle (with shoulders flexed and upper limbs extended above the head-GP1; with shoulders in neutral position and upper limbs extended along the trunk-GP2). Six well-trained breaststroke male swimmers (with reasonable homogeneity of body characteristics) participated in the experimental tests; afterwards a 3D swimmer model was created to fit within the limits of the sample body size profile. The standard k-ε turbulent model was used to simulate the fluid flow around the swimmer model. Velocity ranged from 1.30 to 1.70 m/s for GP1 and 1.10 to 1.50 m/s for GP2. Values found for GP1 and GP2 were lower for CFD than experimental ones. Nevertheless, both CFD and experimental drag/drag coefficient values displayed a tendency to jointly increase/decrease with velocity, except for GP2 CD where CFD and experimental values display opposite tendencies. Results suggest that CFD values obtained by single model approaches should be considered with caution due to small body shape and dimension differences to real swimmers. For better accuracy of CFD studies, realistic individual 3D models of swimmers are required, and specific kinematics respected. PMID:26087879
NASA Astrophysics Data System (ADS)
Peevey, T. R.; Gille, J. C.; Homeyer, C. R.; Manney, G. L.
2014-09-01
Using High Resolution Dynamic Limb Sounder observations and ERA-Interim reanalysis this study demonstrates that the warm conveyor belt (WCB) is a mechanism responsible for the relationship between the double tropopause (DT) and the tropopause inversion layer (TIL), a relationship recently suggested in the literature based on idealized model simulations of baroclinic disturbances. Using these data sets, spatial and temporal characteristics of the DT-TIL relationship are examined over a 3 year period, 2005-2008. In the extratropics, results from satellite data show that as the TIL increases in strength, so does the frequency of the DT, regardless of season or hemisphere. The inverse relationship is found in the tropics. Using only DT profiles, zonal composites of wind, relative vorticity, and temperature from reanalysis data show that as the TIL increases in strength, the upper tropospheric circulation switches from cyclonic to anticyclonic, and the upward vertical motion increases. This result suggests the WCB as a mechanism since it is on the anticyclonic side of the jet and is characterized by the movement of tropical air poleward and upward from the surface. To verify this relationship, the vertical and horizontal development of a synoptic-scale baroclinic system is analyzed over a 4 day period. Results show the equatorward extension of the polar tropopause, and thus the formation of the DT, due to the strengthening of the TIL in the region of vertical motion associated with the WCB. Moreover, this result suggests that air movement within the DT could originate from high latitudes when associated with a baroclinic disturbance.
Time-lapse AVO fluid inversion for dynamic reservoir characterization in Delhi Field, Louisiana
NASA Astrophysics Data System (ADS)
Putri, Indah Hermansyah
In the development stage, CO2 injection is becoming more widely used in enhanced oil recovery (EOR). Delhi Oil Field is part of Phases XIII and XIV of the Reservoir Characterization Project (RCP) Colorado School of Mines. The focus of these phases is to monitor the effectiveness of the CO 2 injection in Delhi Field by using multicomponent time-lapse seismic data. In this study, I analyze the amplitude versus offset (AVO) response of the time-lapse P-wave seismic data in order to quantify the fluid probability in the field. RCP acquired four square miles of multicomponent time-lapse seismic in Delhi Field to characterize the field dynamically. RCP's two surveys, monitor 1 and monitor 2, were shot in 2010 and 2011 after the start of CO2 injection in November 2009. Time-lapse AVO modeling was performed. The modeling results show that both the top Tuscaloosa and Paluxy Formations are class III AVO, and change toward class IV AVO by increasing the CO2 saturation in the reservoir. In addition, the Paluxy Formation shows a consistent result between the synthetic and real data, however, the Tuscaloosa Formation is not consistent as it is affected by tuning. AVO fluid inversion (AFI) was performed on both the Tuscaloosa and Paluxy Formations in order to quantify the fluid probability in these formations. The inversion results are confirmed by the pseudo gamma ray model, the porosity model, the permeability model, the pressure model, and the production data. In the Tuscaloosa and Paluxy Formations, oil and CO2 are located in the good quality, high porosity, and high permeability sandstones. The presence of CO2 is also confirmed by the pressure interpretation. Furthermore, production data from both Tuscaloosa and Paluxy Formations confirm the fluid presence in the reservoir.
Bayesian inversions of a dynamic vegetation model in four European grassland sites
NASA Astrophysics Data System (ADS)
Minet, J.; Laloy, E.; Tychon, B.; François, L.
2015-01-01
Eddy covariance data from four European grassland sites are used to probabilistically invert the CARAIB dynamic vegetation model (DVM) with ten unknown parameters, using the DREAM(ZS) Markov chain Monte Carlo (MCMC) sampler. We compare model inversions considering both homoscedastic and heteroscedastic eddy covariance residual errors, with variances either fixed a~priori or jointly inferred with the model parameters. Agreements between measured and simulated data during calibration are comparable with previous studies, with root-mean-square error (RMSE) of simulated daily gross primary productivity (GPP), ecosystem respiration (RECO) and evapotranspiration (ET) ranging from 1.73 to 2.19 g C m-2 day-1, 1.04 to 1.56 g C m-2 day-1, and 0.50 to 1.28 mm day-1, respectively. In validation, mismatches between measured and simulated data are larger, but still with Nash-Sutcliffe efficiency scores above 0.5 for three out of the four sites. Although measurement errors associated with eddy covariance data are known to be heteroscedastic, we showed that assuming a classical linear heteroscedastic model of the residual errors in the inversion do not fully remove heteroscedasticity. Since the employed heteroscedastic error model allows for larger deviations between simulated and measured data as the magnitude of the measured data increases, this error model expectedly lead to poorer data fitting compared to inversions considering a constant variance of the residual errors. Furthermore, sampling the residual error variances along with model parameters results in overall similar model parameter posterior distributions as those obtained by fixing these variances beforehand, while slightly improving model performance. Despite the fact that the calibrated model is generally capable of fitting the data within measurement errors, systematic bias in the model simulations are observed. These are likely due to model inadequacies such as shortcomings in the photosynthesis modelling
Glassy Dynamics in the Adaptive Immune Response Prevents Autoimmune Disease
NASA Astrophysics Data System (ADS)
Sun, Jun; Deem, Michael
2006-03-01
The immune system normally protects the human host against death by infection. However, when an immune response is mistakenly directed at self antigens, autoimmune disease can occur. We describe a model of protein evolution to simulate the dynamics of the adaptive immune response to antigens. Computer simulations of the dynamics of antibody evolution show that different evolutionary mechanisms, namely gene segment swapping and point mutation, lead to different evolved antibody binding affinities. Although a combination of gene segment swapping and point mutation can yield a greater affinity to a specific antigen than point mutation alone, the antibodies so evolved are highly cross-reactive and would cause autoimmune disease, and this is not the chosen dynamics of the immune system. We suggest that in the immune system a balance has evolved between binding affinity and specificity in the mechanism for searching the amino acid sequence space of antibodies. Our model predicts that chronic infection may lead to autoimmune disease as well due to cross-reactivity and suggests a broad distribution for the time of onset of autoimmune disease due to chronic exposure. The slow search of antibody sequence space by point mutation leads to the broad of distribution times.
Patient-adaptive lesion metabolism analysis by dynamic PET images.
Gao, Fei; Liu, Huafeng; Shi, Pengcheng
2012-01-01
Dynamic PET imaging provides important spatial-temporal information for metabolism analysis of organs and tissues, and generates a great reference for clinical diagnosis and pharmacokinetic analysis. Due to poor statistical properties of the measurement data in low count dynamic PET acquisition and disturbances from surrounding tissues, identifying small lesions inside the human body is still a challenging issue. The uncertainties in estimating the arterial input function will also limit the accuracy and reliability of the metabolism analysis of lesions. Furthermore, the sizes of the patients and the motions during PET acquisition will yield mismatch against general purpose reconstruction system matrix, this will also affect the quantitative accuracy of metabolism analyses of lesions. In this paper, we present a dynamic PET metabolism analysis framework by defining a patient adaptive system matrix to improve the lesion metabolism analysis. Both patient size information and potential small lesions are incorporated by simulations of phantoms of different sizes and individual point source responses. The new framework improves the quantitative accuracy of lesion metabolism analysis, and makes the lesion identification more precisely. The requirement of accurate input functions is also reduced. Experiments are conducted on Monte Carlo simulated data set for quantitative analysis and validation, and on real patient scans for assessment of clinical potential. PMID:23286175
Dynamic modeling and adaptive control for space stations
NASA Technical Reports Server (NTRS)
Ih, C. H. C.; Wang, S. J.
1985-01-01
Of all large space structural systems, space stations present a unique challenge and requirement to advanced control technology. Their operations require control system stability over an extremely broad range of parameter changes and high level of disturbances. During shuttle docking the system mass may suddenly increase by more than 100% and during station assembly the mass may vary even more drastically. These coupled with the inherent dynamic model uncertainties associated with large space structural systems require highly sophisticated control systems that can grow as the stations evolve and cope with the uncertainties and time-varying elements to maintain the stability and pointing of the space stations. The aspects of space station operational properties are first examined, including configurations, dynamic models, shuttle docking contact dynamics, solar panel interaction, and load reduction to yield a set of system models and conditions. A model reference adaptive control algorithm along with the inner-loop plant augmentation design for controlling the space stations under severe operational conditions of shuttle docking, excessive model parameter errors, and model truncation are then investigated. The instability problem caused by the zero-frequency rigid body modes and a proposed solution using plant augmentation are addressed. Two sets of sufficient conditions which guarantee the globablly asymptotic stability for the space station systems are obtained.
Takamuku, Shinya; Gomi, Hiroaki
2015-07-22
How our central nervous system (CNS) learns and exploits relationships between force and motion is a fundamental issue in computational neuroscience. While several lines of evidence have suggested that the CNS predicts motion states and signals from motor commands for control and perception (forward dynamics), it remains controversial whether it also performs the 'inverse' computation, i.e. the estimation of force from motion (inverse dynamics). Here, we show that the resistive sensation we experience while moving a delayed cursor, perceived purely from the change in visual motion, provides evidence of the inverse computation. To clearly specify the computational process underlying the sensation, we systematically varied the visual feedback and examined its effect on the strength of the sensation. In contrast to the prevailing theory that sensory prediction errors modulate our perception, the sensation did not correlate with errors in cursor motion due to the delay. Instead, it correlated with the amount of exposure to the forward acceleration of the cursor. This indicates that the delayed cursor is interpreted as a mechanical load, and the sensation represents its visually implied reaction force. Namely, the CNS automatically computes inverse dynamics, using visually detected motions, to monitor the dynamic forces involved in our actions. PMID:26156766
Dynamic adaptive learning for decision-making supporting systems
NASA Astrophysics Data System (ADS)
He, Haibo; Cao, Yuan; Chen, Sheng; Desai, Sachi; Hohil, Myron E.
2008-03-01
This paper proposes a novel adaptive learning method for data mining in support of decision-making systems. Due to the inherent characteristics of information ambiguity/uncertainty, high dimensionality and noisy in many homeland security and defense applications, such as surveillances, monitoring, net-centric battlefield, and others, it is critical to develop autonomous learning methods to efficiently learn useful information from raw data to help the decision making process. The proposed method is based on a dynamic learning principle in the feature spaces. Generally speaking, conventional approaches of learning from high dimensional data sets include various feature extraction (principal component analysis, wavelet transform, and others) and feature selection (embedded approach, wrapper approach, filter approach, and others) methods. However, very limited understandings of adaptive learning from different feature spaces have been achieved. We propose an integrative approach that takes advantages of feature selection and hypothesis ensemble techniques to achieve our goal. Based on the training data distributions, a feature score function is used to provide a measurement of the importance of different features for learning purpose. Then multiple hypotheses are iteratively developed in different feature spaces according to their learning capabilities. Unlike the pre-set iteration steps in many of the existing ensemble learning approaches, such as adaptive boosting (AdaBoost) method, the iterative learning process will automatically stop when the intelligent system can not provide a better understanding than a random guess in that particular subset of feature spaces. Finally, a voting algorithm is used to combine all the decisions from different hypotheses to provide the final prediction results. Simulation analyses of the proposed method on classification of different US military aircraft databases show the effectiveness of this method.
Community dynamics of cellulose-adapted thermophilic bacterial consortia.
Eichorst, Stephanie A; Varanasi, Patanjali; Stavila, Vatalie; Zemla, Marcin; Auer, Manfred; Singh, Seema; Simmons, Blake A; Singer, Steven W
2013-09-01
Enzymatic hydrolysis of cellulose is a key process in the global carbon cycle and the industrial conversion of biomass to biofuels. In natural environments, cellulose hydrolysis is predominately performed by microbial communities. However, detailed understanding of bacterial cellulose hydrolysis is primarily confined to a few model isolates. Developing models for cellulose hydrolysis by mixed microbial consortia will complement these isolate studies and may reveal new mechanisms for cellulose deconstruction. Microbial communities were adapted to microcrystalline cellulose under aerobic, thermophilic conditions using green waste compost as the inoculum to study cellulose hydrolysis in a microbial consortium. This adaptation selected for three dominant taxa--the Firmicutes, Bacteroidetes and Thermus. A high-resolution profile of community development during the enrichment demonstrated a community transition from Firmicutes to a novel Bacteroidetes population that clusters in the Chitinophagaceae family. A representative strain of this population, strain NYFB, was successfully isolated, and sequencing of a nearly full-length 16S rRNA gene demonstrated that it was only 86% identical compared with other validated strains in the phylum Bacteroidetes. Strain NYFB grew well on soluble polysaccharide substrates, but grew poorly on insoluble polysaccharide substrates. Similar communities were observed in companion thermophilic enrichments on insoluble wheat arabinoxylan, a hemicellulosic substrate, suggesting a common model for deconstruction of plant polysaccharides. Combining observations of community dynamics and the physiology of strain NYFB, a cooperative successional model for polysaccharide hydrolysis by the Firmicutes and Bacteroidetes in the thermophilic cellulolytic consortia is proposed. PMID:23763762
PAQ: Persistent Adaptive Query Middleware for Dynamic Environments
NASA Astrophysics Data System (ADS)
Rajamani, Vasanth; Julien, Christine; Payton, Jamie; Roman, Gruia-Catalin
Pervasive computing applications often entail continuous monitoring tasks, issuing persistent queries that return continuously updated views of the operational environment. We present PAQ, a middleware that supports applications' needs by approximating a persistent query as a sequence of one-time queries. PAQ introduces an integration strategy abstraction that allows composition of one-time query responses into streams representing sophisticated spatio-temporal phenomena of interest. A distinguishing feature of our middleware is the realization that the suitability of a persistent query's result is a function of the application's tolerance for accuracy weighed against the associated overhead costs. In PAQ, programmers can specify an inquiry strategy that dictates how information is gathered. Since network dynamics impact the suitability of a particular inquiry strategy, PAQ associates an introspection strategy with a persistent query, that evaluates the quality of the query's results. The result of introspection can trigger application-defined adaptation strategies that alter the nature of the query. PAQ's simple API makes developing adaptive querying systems easily realizable. We present the key abstractions, describe their implementations, and demonstrate the middleware's usefulness through application examples and evaluation.
Fully Threaded Tree for Adaptive Refinement Fluid Dynamics Simulations
NASA Technical Reports Server (NTRS)
Khokhlov, A. M.
1997-01-01
A fully threaded tree (FTT) for adaptive refinement of regular meshes is described. By using a tree threaded at all levels, tree traversals for finding nearest neighbors are avoided. All operations on a tree including tree modifications are O(N), where N is a number of cells, and are performed in parallel. An efficient implementation of the tree is described that requires 2N words of memory. A filtering algorithm for removing high frequency noise during mesh refinement is described. A FTT can be used in various numerical applications. In this paper, it is applied to the integration of the Euler equations of fluid dynamics. An adaptive mesh time stepping algorithm is described in which different time steps are used at different l evels of the tree. Time stepping and mesh refinement are interleaved to avoid extensive buffer layers of fine mesh which were otherwise required ahead of moving shocks. Test examples are presented, and the FTT performance is evaluated. The three dimensional simulation of the interaction of a shock wave and a spherical bubble is carried out that shows the development of azimuthal perturbations on the bubble surface.
NASA Astrophysics Data System (ADS)
Rivest-Hénault, David; Dowson, Nicholas; Greer, Peter; Dowling, Jason
2014-03-01
MRI-alone treatment planning and adaptive MRI-based prostate radiation therapy are two promising techniques that could significantly increase the accuracy of the curative dose delivery processes while reducing the total radiation dose. State-of-the-art methods rely on the registration of a patient MRI with a MR-CT atlas for the estimation of pseudo-CT [5]. This atlas itself is generally created by registering many CT and MRI pairs. Most registration methods are not symmetric, but the order of the images influences the result [8]. The computed transformation is therefore biased, introducing unwanted variability. This work examines how much a symmetric algorithm improves the registration. Methods: A robust symmetric registration algorithm is proposed that simultaneously optimises a half space transform and its inverse. During the registration process, the two input volumetric images are transformed to a common position in space, therefore minimising any computational bias. An asymmetrical implementation of the same algorithm was used for comparison purposes. Results: Whole pelvis MRI and CT scans from 15 prostate patients were registered, as in the creation of MR-CT atlases. In each case, two registrations were performed, with different input image orders, and the transformation error quantified. Mean residuals of 0.63±0.26 mm (translation) and (8.7±7.3) × 10--3 rad (rotation) were found for the asymmetrical implementation with corresponding values of 0.038±0.039 mm and (1.6 ± 1.3) × 10--3 rad for the proposed symmetric algorithm, a substantial improvement. Conclusions: The increased registration precision will enhance the generation of pseudo-CT from MRI for atlas based MR planning methods.
He, Qixin; Knowles, L Lacey
2016-05-01
Chromosomal inversions are important structural changes that may facilitate divergent selection when they capture co-adaptive loci in the face of gene flow. However, identifying selection targets within inversions can be challenging. The high degrees of differentiation between heterokaryotypes, as well as the differences in demographic histories of collinear regions compared with inverted ones, reduce the power of traditional outlier analyses for detecting selected loci. Here, we develop a new approach that uses discriminant functions informed from inversion-specific expectations to classify loci that are under selection (or drift). Analysis of RAD sequencing data we collected in a classic dipteran species with polymorphic inversion clines-Anopheles gambiae, a malaria vector species from sub-Saharan Africa-demonstrates the benefits of the approach compared with traditional outlier analyses. We focus specifically on two polymorphic inversions, the 2La and 2Rb arrangements that predominate in dry habitats and the 2L+(a) and 2R+(b) arrangements in wet habitats, which contrast with the minimal geographic structure of SNPs from collinear regions. With our approach, we identify two strongly selected regions within 2La associated with dry habitat. Moreover, we also show that the prevalence of selection is greater in the arrangement 2L+(a) that is associated with wet habitat (unlike presumed importance of selective divergence associated with the shift of the mosquitoes into dry habitats). We discuss the implications of these results with respect to studies of rapid adaptation in these malaria vectors, and in particular, the insights our newly developed approach offers for identifying not only potential targets of selection, but also the population that has undergone adaptive change. PMID:26994406
Jiang, Xiaoming; Van den Broek, Wouter; Koch, Christoph T
2016-04-01
Inverse dynamical photon scattering (IDPS), an artificial neural network based algorithm for three-dimensional quantitative imaging in optical microscopy, is introduced. Because the inverse problem entails numerical minimization of an explicit error metric, it becomes possible to freely choose a more robust metric, to introduce regularization of the solution, and to retrieve unknown experimental settings or microscope values, while the starting guess is simply set to zero. The regularization is accomplished through an alternate directions augmented Lagrangian approach, implemented on a graphics processing unit. These improvements are demonstrated on open source experimental data, retrieving three-dimensional amplitude and phase for a thick specimen. PMID:27136994
Wiebeler, Christian; Schumacher, Stefan
2014-09-11
Photochromism is a light-induced molecular process that is likely to find its way into future optoelectronic devices. In further optimization of photochromic materials, light-induced conversion efficiencies as well as reaction times can usually only be determined once a new molecule was synthesized. Here we use nonadiabatic ab initio molecular dynamics to study the electrocyclic reaction of diarylethenes, comparing normal- and inverse-type systems. Our study highlights that reaction quantum yields can be successfully predicted in accord with experimental findings. In particular, we find that inverse-type diarylethenes show a significantly higher reaction quantum yield and cycloreversion on times typically as short as 100 fs. PMID:25140609
Dynamic skeletal muscle stimulation and its potential in bone adaptation
Qin, Y-X.; Lam, H.; Ferreri, S.; Rubin, C.
2016-01-01
To identify mechanotransductive signals for combating musculoskeletal deterioration, it is essential to determine the components and mechanisms critical to the anabolic processes of musculoskeletal tissues. It is hypothesized that the interaction between bone and muscle may depend on fluid exchange in these tissues by mechanical loading. It has been shown that intramedullary pressure (ImP) and low-level bone strain induced by muscle stimulation (MS) has the potential to mitigate bone loss induced by disuse osteopenia. Optimized MS signals, i.e., low-intensity and high frequency, may be critical in maintaining bone mass and mitigating muscle atrophy. The objectives for this review are to discuss the potential for MS to induce ImP and strains on bone, to regulate bone adaptation, and to identify optimized stimulation frequency in the loading regimen. The potential for MS to regulate blood and fluid flow will also be discussed. The results suggest that oscillatory MS regulates fluid dynamics with minimal mechanical strain in bone. The response was shown to be dependent on loading frequency, serving as a critical mediator in mitigating bone loss. A specific regimen of dynamic MS may be optimized in vivo to attenuate disuse osteopenia and serve as a biomechanical intervention in the clinical setting. PMID:20190376
On Cognition, Structured Sequence Processing, and Adaptive Dynamical Systems
NASA Astrophysics Data System (ADS)
Petersson, Karl Magnus
2008-11-01
Cognitive neuroscience approaches the brain as a cognitive system: a system that functionally is conceptualized in terms of information processing. We outline some aspects of this concept and consider a physical system to be an information processing device when a subclass of its physical states can be viewed as representational/cognitive and transitions between these can be conceptualized as a process operating on these states by implementing operations on the corresponding representational structures. We identify a generic and fundamental problem in cognition: sequentially organized structured processing. Structured sequence processing provides the brain, in an essential sense, with its processing logic. In an approach addressing this problem, we illustrate how to integrate levels of analysis within a framework of adaptive dynamical systems. We note that the dynamical system framework lends itself to a description of asynchronous event-driven devices, which is likely to be important in cognition because the brain appears to be an asynchronous processing system. We use the human language faculty and natural language processing as a concrete example through out.