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.
Adaptive Fault-Tolerant Control of Uncertain Nonlinear Large-Scale Systems With Unknown Dead Zone.
Chen, Mou; Tao, Gang
2016-08-01
In this paper, an adaptive neural fault-tolerant control scheme is proposed and analyzed for a class of uncertain nonlinear large-scale systems with unknown dead zone and external disturbances. To tackle the unknown nonlinear interaction functions in the large-scale system, the radial basis function neural network (RBFNN) is employed to approximate them. To further handle the unknown approximation errors and the effects of the unknown dead zone and external disturbances, integrated as the compounded disturbances, the corresponding disturbance observers are developed for their estimations. Based on the outputs of the RBFNN and the disturbance observer, the adaptive neural fault-tolerant control scheme is designed for uncertain nonlinear large-scale systems by using a decentralized backstepping technique. The closed-loop stability of the adaptive control system is rigorously proved via Lyapunov analysis and the satisfactory tracking performance is achieved under the integrated effects of unknown dead zone, actuator fault, and unknown external disturbances. Simulation results of a mass-spring-damper system are given to illustrate the effectiveness of the proposed adaptive neural fault-tolerant control scheme for uncertain nonlinear large-scale systems. PMID:26340792
Liu, Yan-Jun; Gao, Ying; Tong, Shaocheng; Chen, C L Philip
2016-01-01
In this paper, an effective adaptive control approach is constructed to stabilize a class of nonlinear discrete-time systems, which contain unknown functions, unknown dead-zone input, and unknown control direction. Different from linear dead zone, the dead zone, in this paper, is a kind of nonlinear dead zone. To overcome the noncausal problem, which leads to the control scheme infeasible, the systems can be transformed into a m -step-ahead predictor. Due to nonlinear dead-zone appearance, the transformed predictor still contains the nonaffine function. In addition, it is assumed that the gain function of dead-zone input and the control direction are unknown. These conditions bring about the difficulties and the complicacy in the controller design. Thus, the implicit function theorem is applied to deal with nonaffine dead-zone appearance, the problem caused by the unknown control direction can be resolved through applying the discrete Nussbaum gain, and the neural networks are used to approximate the unknown function. Based on the Lyapunov theory, all the signals of the resulting closed-loop system are proved to be semiglobal uniformly ultimately bounded. Moreover, the tracking error is proved to be regulated to a small neighborhood around zero. The feasibility of the proposed approach is demonstrated by a simulation example. PMID:26353383
Peng, Jinzhu; Dubay, Rickey
2011-10-01
In this paper, an adaptive control approach based on the neural networks is presented to control a DC motor system with dead-zone characteristics (DZC), where two neural networks are proposed to formulate the traditional identification and control approaches. First, a Wiener-type neural network (WNN) is proposed to identify the motor DZC, which formulates the Wiener model with a linear dynamic block in cascade with a nonlinear static gain. Second, a feedforward neural network is proposed to formulate the traditional PID controller, termed as PID-type neural network (PIDNN), which is then used to control and compensate for the DZC. In this way, the DC motor system with DZC is identified by the WNN identifier, which provides model information to the PIDNN controller in order to make it adaptive. Back-propagation algorithms are used to train both neural networks. Also, stability and convergence analysis are conducted using the Lyapunov theorem. Finally, experiments on the DC motor system demonstrated accurate identification and good compensation for dead-zone with improved control performance over the conventional PID control. PMID:21788017
Observed-Based Adaptive Fuzzy Tracking Control for Switched Nonlinear Systems With Dead-Zone.
Tong, Shaocheng; Sui, Shuai; Li, Yongming
2015-12-01
In this paper, the problem of adaptive fuzzy output-feedback control is investigated for a class of uncertain switched nonlinear systems in strict-feedback form. The considered switched systems contain unknown nonlinearities, dead-zone, and immeasurable states. Fuzzy logic systems are utilized to approximate the unknown nonlinear functions, a switched fuzzy state observer is designed and thus the immeasurable states are obtained by it. By applying the adaptive backstepping design principle and the average dwell time method, an adaptive fuzzy output-feedback tracking control approach is developed. It is proved that the proposed control approach can guarantee that all the variables in the closed-loop system are bounded under a class of switching signals with average dwell time, and also that the system output can track a given reference signal as closely as possible. The simulation results are given to check the effectiveness of the proposed approach. PMID:25594991
Liu, Yan-Jun; Tong, Shaocheng
2015-03-01
In the paper, an adaptive tracking control design is studied for a class of nonlinear discrete-time systems with dead-zone input. The considered systems are of the nonaffine pure-feedback form and the dead-zone input appears nonlinearly in the systems. The contributions of the paper are that: 1) it is for the first time to investigate the control problem for this class of discrete-time systems with dead-zone; 2) there are major difficulties for stabilizing such systems and in order to overcome the difficulties, the systems are transformed into an n-step-ahead predictor but nonaffine function is still existent; and 3) an adaptive compensative term is constructed to compensate for the parameters of the dead-zone. The neural networks are used to approximate the unknown functions in the transformed systems. Based on the Lyapunov theory, it is proven that all the signals in the closed-loop system are semi-globally uniformly ultimately bounded and the tracking error converges to a small neighborhood of zero. Two simulation examples are provided to verify the effectiveness of the control approach in the paper. PMID:24968366
FRIT for Systems with Dead-Zone and Its Application to Ultrasonic Motors
NASA Astrophysics Data System (ADS)
Wakasa, Yuji; Kanagawa, Shinji; Tanaka, Kanya; Nishimura, Yuki
Ultrasonic motors (USMs) intrinsically have a dead-zone property which is sensitive to load changes. This paper proposes a fictitious reference iterative tuning (FRIT) method for systems with a dead-zone property such as USMs. The standard FRIT method is basically developed for linear systems and may not give a satisfactory control performance for noninvertible nonlinear systems including USMs. In contrast, the proposed FRIT method can achieve such a performance by introducing a right inverse of a dead-zone function as a dead-zone compensator. In the optimization process of FRIT, the so-called covariance matrix adaptation evolution strategy (CMA-ES) is used for simultaneously searching a dead-zone parameter as well as controller parameters. CMA-ES is a kind of stochastic multi-point search techniques and is effective for nondifferentiable and nonconvex optimization problems. Experimental results for a USM are given to show the effectiveness of the proposed method.
Investigating Aquatic Dead Zones
ERIC Educational Resources Information Center
Testa, Jeremy; Gurbisz, Cassie; Murray, Laura; Gray, William; Bosch, Jennifer; Burrell, Chris; Kemp, Michael
2010-01-01
This article features two engaging high school activities that include current scientific information, data, and authentic case studies. The activities address the physical, biological, and chemical processes that are associated with oxygen-depleted areas, or "dead zones," in aquatic systems. Students can explore these dead zones through both…
Bagherpoor, H M; Salmasi, Farzad R
2015-07-01
In this paper, robust model reference adaptive tracking controllers are considered for Single-Input Single-Output (SISO) and Multi-Input Multi-Output (MIMO) linear systems containing modeling uncertainties, unknown additive disturbances and actuator fault. Two new lemmas are proposed for both SISO and MIMO, under which dead-zone modification rule is improved such that the tracking error for any reference signal tends to zero in such systems. In the conventional approach, adaption of the controller parameters is ceased inside the dead-zone region which results tracking error, while preserving the system stability. In the proposed scheme, control signal is reinforced with an additive term based on tracking error inside the dead-zone which results in full reference tracking. In addition, no Fault Detection and Diagnosis (FDD) unit is needed in the proposed approach. Closed loop system stability and zero tracking error are proved by considering a suitable Lyapunov functions candidate. It is shown that the proposed control approach can assure that all the signals of the close loop system are bounded in faulty conditions. Finally, validity and performance of the new schemes have been illustrated through numerical simulations of SISO and MIMO systems in the presence of actuator faults, modeling uncertainty and output disturbance. PMID:25744053
Climate change and dead zones.
Altieri, Andrew H; Gedan, Keryn B
2015-04-01
Estuaries and coastal seas provide valuable ecosystem services but are particularly vulnerable to the co-occurring threats of climate change and oxygen-depleted dead zones. We analyzed the severity of climate change predicted for existing dead zones, and found that 94% of dead zones are in regions that will experience at least a 2 °C temperature increase by the end of the century. We then reviewed how climate change will exacerbate hypoxic conditions through oceanographic, ecological, and physiological processes. We found evidence that suggests numerous climate variables including temperature, ocean acidification, sea-level rise, precipitation, wind, and storm patterns will affect dead zones, and that each of those factors has the potential to act through multiple pathways on both oxygen availability and ecological responses to hypoxia. Given the variety and strength of the mechanisms by which climate change exacerbates hypoxia, and the rates at which climate is changing, we posit that climate change variables are contributing to the dead zone epidemic by acting synergistically with one another and with recognized anthropogenic triggers of hypoxia including eutrophication. This suggests that a multidisciplinary, integrated approach that considers the full range of climate variables is needed to track and potentially reverse the spread of dead zones. PMID:25385668
Long-Time Sustainability of Rossby Wave Instability in Protoplanetary Disks with Dead Zone
NASA Astrophysics Data System (ADS)
Li, S.; Li, H.
2015-10-01
We have run 2D simulations to investigate the generation and sustainability of Rossby wave instability (RWI) in proto-planetary disks with constant viscosity and for disks with low viscosity regions (dead zone). For the constant viscosity case, the development of RWI requires a low viscosity and life time of the RWI is short. We also find that the vortex, when it migrates, does so much faster than the disk's viscous drift rate. For disks with dead zone case, a much larger viscosity can be used and the RWI vortex can be sustained for a long time, even the life time of the disk, depending on the width and depth of the dead zone. For a narrow dead zone, the vortex depicts a periodic pattern with a period inversely proportional to the viscosity. If the dead-zone width exceeds some threshold, the periodicity of the RWI disappears.
Dead Zone Accretion Flows in Protostellar Disks
NASA Technical Reports Server (NTRS)
Turner, Neal; Sano, T.
2008-01-01
Planets form inside protostellar disks in a dead zone where the electrical resistivity of the gas is too high for magnetic forces to drive turbulence. We show that much of the dead zone nevertheless is active and flows toward the star while smooth, large-scale magnetic fields transfer the orbital angular momentum radially outward. Stellar X-ray and radionuclide ionization sustain a weak coupling of the dead zone gas to the magnetic fields, despite the rapid recombination of free charges on dust grains. Net radial magnetic fields are generated in the magnetorotational turbulence in the electrically conducting top and bottom surface layers of the disk, and reach the midplane by ohmic diffusion. A toroidal component to the fields is produced near the midplane by the orbital shear. The process is similar to the magnetization of the solar tachocline. The result is a laminar, magnetically driven accretion flow in the region where the planets form.
ON HYDRODYNAMIC MOTIONS IN DEAD ZONES
Oishi, Jeffrey S.; Mac Low, Mordecai-Mark E-mail: mordecai@amnh.or
2009-10-20
We investigate fluid motions near the midplane of vertically stratified accretion disks with highly resistive midplanes. In such disks, the magnetorotational instability drives turbulence in thin layers surrounding a resistive, stable dead zone. The turbulent layers in turn drive motions in the dead zone. We examine the properties of these motions using three-dimensional, stratified, local, shearing-box, non-ideal, magnetohydrodynamical simulations. Although the turbulence in the active zones provides a source of vorticity to the midplane, no evidence for coherent vortices is found in our simulations. It appears that this is because of strong vertical oscillations in the dead zone. By analyzing time series of azimuthally averaged flow quantities, we identify an axisymmetric wave mode particular to models with dead zones. This mode is reduced in amplitude, but not suppressed entirely, by changing the equation of state from isothermal to ideal. These waves are too low frequency to affect sedimentation of dust to the midplane, but may have significance for the gravitational stability of the resulting midplane dust layers.
Cheatgrass Dead Zones in Northern Nevada
Technology Transfer Automated Retrieval System (TEKTRAN)
Reports of areas of cheatgrass die-off are becoming more frequent. In 2009, we investigated cheatgrass die-off in north-central Nevada. Dead zones ranged from several to hundreds of acres in size and were largely unvegetated and covered by cheatgrass litter with a distinct gray cast. We collected re...
Gulf of Mexico dead zone - the last 150 years
Osterman, Lisa; Swarzenski, P.W.; Poore, R.Z.
2006-01-01
'Gulf of Mexico Dead Zone-The Last 150 Years' discusses the dead zone that forms seasonally in the northern Gulf of Mexico when subsurface waters become depleted in dissolved oxygen and cannot support most life.
Improving measurement of Chesapeake Bay's dead zone
NASA Astrophysics Data System (ADS)
Schultz, Colin
2013-09-01
In the 1930s, researchers first noticed that the Chesapeake Bay had a dead zone, an expanse of water with drastically reduced concentrations of oxygen. In the 1980s, hypoxia—low-oxygen conditions—gave way in some places to anoxia—a near-total depletion of dissolved oxygen. A lack of oxygen makes the water inhospitable for many marine organisms, and the Chesapeake Bay is the focus of major ecosystem rehabilitation efforts.
Recipe for Hypoxia: Playing the Dead Zone Game
ERIC Educational Resources Information Center
Kastler, Jessica A.
2009-01-01
Dead zones--areas experiencing low levels of dissolved oxygen--are growing in shallow ocean waters around the world. Research has shown that dead zones form as a result of a specific type of pollution, called nutrient enrichment or eutrophication, and are found in almost every coastal zone where humans have large populations. Concepts related to…
Turbulence, Transport, and Waves in Ohmic Dead Zones
NASA Astrophysics Data System (ADS)
Gole, Daniel; Simon, Jacob B.; Lubow, Stephen H.; Armitage, Philip J.
2016-07-01
We use local numerical simulations to study a vertically stratified accretion disk with a resistive mid-plane that damps magnetohydrodynamic (MHD) turbulence. This is an idealized model for the dead zones that may be present at some radii in protoplanetary and dwarf novae disks. We vary the relative thickness of the dead and active zones to quantify how forced fluid motions in the dead zone change. We find that the residual Reynolds stress near the mid-plane decreases with increasing dead zone thickness, becoming negligible in cases where the active to dead mass ratio is less than a few percent. This implies that purely Ohmic dead zones would be vulnerable to episodic accretion outbursts via the mechanism of Martin & Lubow. We show that even thick dead zones support a large amount of kinetic energy, but this energy is largely in fluid motions that are inefficient at angular momentum transport. Confirming results from Oishi & Mac Low, the perturbed velocity field in the dead zone is dominated by an oscillatory, vertically extended circulation pattern with a low frequency compared to the orbital frequency. This disturbance has the properties predicted for the lowest order r mode in a hydrodynamic disk. We suggest that in a global disk similar excitations would lead to propagating waves, whose properties would vary with the thickness of the dead zone and the nature of the perturbations (isothermal or adiabatic). Flows with similar amplitudes would buckle settled particle layers and could reduce the efficiency of pebble accretion.
Uncertainty evaluation of dead zone of diagnostic ultrasound equipment
NASA Astrophysics Data System (ADS)
Souza, R. M.; Alvarenga, A. V.; Braz, D. S.; Petrella, L. I.; Costa-Felix, R. P. B.
2016-07-01
This paper presents a model for evaluating measurement uncertainty of a feature used in the assessment of ultrasound images: dead zone. The dead zone was measured by two technicians of the INMETRO's Laboratory of Ultrasound using a phantom and following the standard IEC/TS 61390. The uncertainty model was proposed based on the Guide to the Expression of Uncertainty in Measurement. For the tested equipment, results indicate a dead zone of 1.01 mm, and based on the proposed model, the expanded uncertainty was 0.17 mm. The proposed uncertainty model contributes as a novel way for metrological evaluation of diagnostic imaging by ultrasound.
Metagenomic insights into important microbes from the Dead Zone
NASA Astrophysics Data System (ADS)
Thrash, C.; Baker, B.; Seitz, K.; Temperton, B.; Gillies, L.; Rabalais, N. N.; Mason, O. U.
2015-12-01
Coastal regions of eutrophication-driven oxygen depletion are widespread and increasing in number. Also known as dead zones, these regions take their name from the deleterious effects of hypoxia (dissolved oxygen less than 2 mg/L) on shrimp, demersal fish, and other animal life. Dead zones result from nutrient enrichment of primary production, concomitant consumption by chemoorganotrophic aerobic microorganisms, and strong stratification that prevents ventilation of bottom water. One of the largest dead zones in the world occurs seasonally in the northern Gulf of Mexico (nGOM), where hypoxia can reach up to 22,000 square kilometers. While this dead zone shares many features with more well-known marine oxygen minimum zones, it is nevertheless understudied with regards to the microbial assemblages involved in biogeochemical cycling. We performed metagenomic and metatranscriptomic sequencing on six samples from the 2013 nGOM dead zone from both hypoxic and oxic bottom waters. Assembly and binning led to the recovery of over fifty partial to nearly complete metagenomes from key microbial taxa previously determined to be numerically abundant from 16S rRNA data, such as Thaumarcheaota, Marine Group II Euryarchaeota, SAR406, SAR324, Synechococcus spp., and Planctomycetes. These results provide information about the roles of these taxa in the nGOM dead zone, and opportunities for comparing this region of low oxygen to others around the globe.
MAGNETIZED ACCRETION AND DEAD ZONES IN PROTOSTELLAR DISKS
Dzyurkevich, Natalia; Henning, Thomas; Turner, Neal J.; Kley, Wilhelm
2013-03-10
The edges of magnetically dead zones in protostellar disks have been proposed as locations where density bumps may arise, trapping planetesimals and helping form planets. Magneto-rotational turbulence in magnetically active zones provides both accretion of gas on the star and transport of mass to the dead zone. We investigate the location of the magnetically active regions in a protostellar disk around a solar-type star, varying the disk temperature, surface density profile, and dust-to-gas ratio. We also consider stellar masses between 0.4 and 2 M{sub Sun }, with corresponding adjustments in the disk mass and temperature. The dead zone's size and shape are found using the Elsasser number criterion with conductivities including the contributions from ions, electrons, and charged fractal dust aggregates. The charged species' abundances are found using the approach proposed by Okuzumi. The dead zone is in most cases defined by the ambipolar diffusion. In our maps, the dead zone takes a variety of shapes, including a fish tail pointing away from the star and islands located on and off the midplane. The corresponding accretion rates vary with radius, indicating locations where the surface density will increase over time, and others where it will decrease. We show that density bumps do not readily grow near the dead zone's outer edge, independently of the disk parameters and the dust properties. Instead, the accretion rate peaks at the radius where the gas-phase metals freeze out. This could lead to clearing a valley in the surface density, and to a trap for pebbles located just outside the metal freezeout line.
FORMATION OF CIRCUMBINARY PLANETS IN A DEAD ZONE
Martin, Rebecca G.; Armitage, Philip J.; Alexander, Richard D.
2013-08-10
Circumbinary planets have been observed at orbital radii where binary perturbations may have significant effects on the gas disk structure, on planetesimal velocity dispersion, and on the coupling between turbulence and planetesimals. Here, we note that the impact of all of these effects on planet formation is qualitatively altered if the circumbinary disk structure is layered, with a non-turbulent midplane layer (dead zone) and strongly turbulent surface layers. For close binaries, we find that the dead zone typically extends from a radius close to the inner disk edge up to a radius of around 10-20 AU from the center of mass of the binary. The peak in the surface density occurs within the dead zone, far from the inner disk edge, close to the snow line, and may act as a trap for aerodynamically coupled solids. We suggest that circumbinary planet formation may be easier near this preferential location than for disks around single stars. However, dead zones around wide binaries are less likely, and hence planet formation may be more difficult there.
Dead Zones in LX-17 and PBX 9502
Souers, P C; Andreski, H G; Batteux, J; Bratton, B; Cabacungan, C; Cook, III, C F; Fletcher, S; Garza, R; Grimsley, D; Handly, J; Hernandez, A; McMaster, P; Molitoris, J D; Palmer, R; Prindiville, J; Rodriguez, J; Schneberk, D; Wong, B; Vitello, P
2005-09-06
Pin and X-ray corner-turning data have been taken on ambient LX-17 and PBX 9052, and the results are listed in tables as an aid to future modeling. The results have been modeled at 4 zones/mm with a reactive flow approach that varies the burn rate as a function of pressure. A single rate format is used to simulate failure and detonation in different pressure regimes. A pressure cut-off must also be reached to initiate the burn. Corner-turning and failure are modeled using an intermediate pressure rate region, and detonation occurs at high pressure. The TATB booster is also modeled using reactive flow, and X-ray tomography is used to partition the ram-pressed hemisphere into five different density regions. The model reasonably fits the bare corner-turning experiment but predicts a smaller dead zone with steel confinement, in contradiction with experiment. The same model also calculates the confined and unconfined cylinder detonation velocities and predicts the failure of the unconfined cylinder at 3.75 mm radius. The PBX 9502 shows a smaller dead zone than LX-17. An old experiment that showed a large apparent dead zone in Comp B was repeated with X-ray transmission and no dead zone was seen. This confirms the idea that a variable burn rate is the key to modeling. The model also produces initiation delays, which are shorter than those found in time-to-detonation.
DEAD ZONE IN THE POLAR-CAP ACCELERATOR OF PULSARS
Chen, Alexander Y.; Beloborodov, Andrei M.
2013-01-10
We study plasma flows above pulsar polar caps using time-dependent simulations of plasma particles in the self-consistent electric field. The flow behavior is controlled by the dimensionless parameter {alpha} = j/c{rho}{sub GJ}, where j is the electric current density and {rho}{sub GJ} is the Goldreich-Julian charge density. The region of the polar cap where 0 < {alpha} < 1 is a {sup d}ead zone{sup -}in this zone, particle acceleration is inefficient and pair creation is not expected even for young, rapidly rotating pulsars. Pulsars with polar caps near the rotation axis are predicted to have a hollow-cone structure of radio emission, as the dead zone occupies the central part of the polar cap. Our results apply to charge-separated flows of electrons (j < 0) or ions (j > 0). In the latter case, we consider the possibility of a mixed flow consisting of different ion species, and observe the development of two-stream instability. The dead zone at the polar cap is essential for the development of an outer gap near the null surface {rho}{sub GJ} = 0.
Gulf of Mexico dead zone - 1000 year record
Osterman, L.E.; Poore, R.Z.; Swarzenski, P.W.
2010-01-01
An area of oxygen-depleted bottom- and subsurfacewater (hypoxia = dissolved oxygen Since systematic measurement of the extent of the dead zone was begun in 1985, the overall pattern indicates that the area of the dead zone is increasing. Several studies have concluded that the expansion of the Louisiana shelf dead zone is related to increased nutrients (primarily nitrogen, but possibly also phosphorous) in the Mississippi River drainage basin and is responsible for the degradation of Gulf of Mexico marine habitats. The goal of this research is to augment information on the recent expansion of Louisiana shelf hypoxia and to investigate the temporal and geographic extent of the lowoxygen bottom-water conditions prior to 1985 in sediment cores collected from the Louisiana shelf. We use a specific low-oxygen faunal proxy termed the PEB index based on the cumulative percentage of three foraminifers (= % Protononion atlanticum, + % Epistominella vitrea, + % Buliminella morgani) that has been shown statistically to represent the modern seasonal Louisiana hypoxia zone. Our hypothesis is that the increased relative abundance of PEB species in dated sediment cores accurately tracks past seasonal low-oxygen conditions on the Louisiana shelf.
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
The dead zones: oxygen-starved coastal waters.
Joyce, S
2000-03-01
After the great Mississippi River flood of 1993, the hypoxic (or low-oxygen) "dead zone" in the Gulf of Mexico more than doubled its size, reaching an all-time high of over 7,700 square miles in July of 1999. Scientists attribute the Gulf of Mexico dead zone largely to nutrient runoff from agriculture in the Mississippi River basin. During the warm months, these nutrients fuel eutrophication, or high organic production, causing large algal blooms. When the algae decay, the result is hypoxia. Reports of such hypoxic events around the world have been increasing since the mid 1960s. Eutrophication and hypoxia have resulted in mortality of bottom-dwelling life in dozens of marine ecosystems and have stressed fisheries worldwide. Some algal blooms can alter the function of coastal ecosystems or, potentially, threaten human health. Anthropogenic nutrient loading from sources such as agriculture, fossil fuel emissions, and climate events is believed to be related to the global increase in frequency, size, and duration of certain algal blooms. PMID:10706539
Early Jurassic black shales: Global anoxia or regional "Dead Zones"?
NASA Astrophysics Data System (ADS)
van de Schootbrugge, B.; Payne, J.; Wignall, P.
2012-12-01
The so-called "Schwarzer Jura" or "Black Jurassic" in Germany is informally used to designate a series of organic-rich sediments that roughly span the Early Jurassic (201.6 - 175.6 Myr), and which culminate in the Toarcian Oceanic Anoxic Event. Based on organic and inorganic geochemical as well as (micro)palaeontological data from several recently drilled cores, black shales deposited directly following the end-Triassic extinction (201.6 Ma) during the Hettangian are extremely similar to Toarcian black shales. Both events are characterized by laminated black shales that contain high amounts of the biomarker isorenieratane, a fossilized pigment derived from green sulphur bacteria. Furthermore, the two intervals show similar changes in phytoplankton assemblages from chromophyte (red) to chlorophyte (green) algae. Combined, the evidence suggests that photic zone euxinia developed repeatedly during the Early Jurassic, making wide swaths of shelf area inhospitable to benthic life. In the oceans today such areas are called "Dead Zones" and they are increasing in number and extent due to the combined effects of man-made eutrophication and global warming. During the Early Jurassic, regional anoxic events developed in response to flood basalt volcanism, which triggered global warming, increased run-off, and changes in ocean circulation. The patchiness of Early Jurassic anoxia allows comparisons to be made with present-day "Dead Zones", while at the same time ocean de-oxygenation in the past may serve to predict future perturbations in the Earth system.
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.
Dead zone characteristics of SQS and GM modes in Kr-, Ar- and Ne-Mixtures
Nohtomi, A.; Sakae, T.; Matoba, M. ); Koori, N. )
1992-08-01
In this paper, dead zones of a gas counter are studied systematically for various gas mixtures at gas pressures equal to and lower than 1 atm. Evaluated values of the dead zone for the self-quenching streamer (SQS) mode were 190-390 [mu]s [center dot] at 1 atm. Characteristics of the dead zone were explained by considering the effect of SQS ions. The critical gas pressure, at which the operation mode changes form the SQS mode to the GM mode, was determined for each gas mixture on the basis of the dead zone measurement. Operation mode diagrams of the counter were drawn according to those results.
Metagenome of a versatile chemolithoautotroph from expanding oceanic dead zones.
Walsh, David A; Zaikova, Elena; Howes, Charles G; Song, Young C; Wright, Jody J; Tringe, Susannah G; Tortell, Philippe D; Hallam, Steven J
2009-10-23
Oxygen minimum zones, also known as oceanic "dead zones," are widespread oceanographic features currently expanding because of global warming. Although inhospitable to metazoan life, they support a cryptic microbiota whose metabolic activities affect nutrient and trace gas cycling within the global ocean. Here, we report metagenomic analyses of a ubiquitous and abundant but uncultivated oxygen minimum zone microbe (SUP05) related to chemoautotrophic gill symbionts of deep-sea clams and mussels. The SUP05 metagenome harbors a versatile repertoire of genes mediating autotrophic carbon assimilation, sulfur oxidation, and nitrate respiration responsive to a wide range of water-column redox states. Our analysis provides a genomic foundation for understanding the ecological and biogeochemical role of pelagic SUP05 in oxygen-deficient oceanic waters and its potential sensitivity to environmental changes. PMID:19900896
Metagenome of a Versatile Chemolithoautotroph from Expanding Oceanic Dead Zones
Walsh, David A.; Zaikova, Elena; Howes, Charles L.; Song, Young; Wright, Jody; Tringe, Susannah G.; Tortell, Philippe D.; Hallam, Steven J.
2009-07-15
Oxygen minimum zones (OMZs), also known as oceanic"dead zones", are widespread oceanographic features currently expanding due to global warming and coastal eutrophication. Although inhospitable to metazoan life, OMZs support a thriving but cryptic microbiota whose combined metabolic activity is intimately connected to nutrient and trace gas cycling within the global ocean. Here we report time-resolved metagenomic analyses of a ubiquitous and abundant but uncultivated OMZ microbe (SUP05) closely related to chemoautotrophic gill symbionts of deep-sea clams and mussels. The SUP05 metagenome harbors a versatile repertoire of genes mediating autotrophic carbon assimilation, sulfur-oxidation and nitrate respiration responsive to a wide range of water column redox states. Thus, SUP05 plays integral roles in shaping nutrient and energy flow within oxygen-deficient oceanic waters via carbon sequestration, sulfide detoxification and biological nitrogen loss with important implications for marine productivity and atmospheric greenhouse control.
Open ocean dead zones in the tropical North Atlantic Ocean
NASA Astrophysics Data System (ADS)
Karstensen, J.; Fiedler, B.; Schütte, F.; Brandt, P.; Körtzinger, A.; Fischer, G.; Zantopp, R.; Hahn, J.; Visbeck, M.; Wallace, D.
2015-04-01
Here we present first observations, from instrumentation installed on moorings and a float, of unexpectedly low (<2 μmol kg-1) oxygen environments in the open waters of the tropical North Atlantic, a region where oxygen concentration does normally not fall much below 40 μmol kg-1. The low-oxygen zones are created at shallow depth, just below the mixed layer, in the euphotic zone of cyclonic eddies and anticyclonic-modewater eddies. Both types of eddies are prone to high surface productivity. Net respiration rates for the eddies are found to be 3 to 5 times higher when compared with surrounding waters. Oxygen is lowest in the centre of the eddies, in a depth range where the swirl velocity, defining the transition between eddy and surroundings, has its maximum. It is assumed that the strong velocity at the outer rim of the eddies hampers the transport of properties across the eddies boundary and as such isolates their cores. This is supported by a remarkably stable hydrographic structure of the eddies core over periods of several months. The eddies propagate westward, at about 4 to 5 km day-1, from their generation region off the West African coast into the open ocean. High productivity and accompanying respiration, paired with sluggish exchange across the eddy boundary, create the "dead zone" inside the eddies, so far only reported for coastal areas or lakes. We observe a direct impact of the open ocean dead zones on the marine ecosystem as such that the diurnal vertical migration of zooplankton is suppressed inside the eddies.
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.
Severe dead-zone eddies in the open North Atlantic Ocean
NASA Astrophysics Data System (ADS)
Karstensen, J.; Fiedler, B.; Brandt, P.; Körtzinger, A.; Zantopp, R.; Wallace, D.; Krahmann, G.; Visbeck, M.; Hahn, J.
2012-04-01
Ocean volumes with very low dissolved oxygen, so called "dead-zones", have been observed in many coastal areas of the world ocean. Dead-zones are characterized by a dissolved oxygen content below 2 mg/l (approx. 60 µmol/kg) and making them inhabitable for many marine organisms. Here we report on severe dead-zones in the open North Atlantic, several hundreds of kilometres away from the coast, where so far concentrations below about 40 µmol/L have not been reported. The severe dead-zones are contained within mesoscale eddies that originate from the West African upwelling region and propagate slowly (100km per month) westward. Local dynamics isolate the dead-zone eddy from surrounding waters and create, within the rather well oxygenated North Atlantic, a biogeochemical realm comparable to the major oxygen minimum zone (OMZ) of the Pacific and Indian Ocean. Below a well oxygenated upper mixed-layer of some 20 to 50m depth follows a drastic drop in oxygen, which is the actual dead-zone. In one the most dramatic case of a North Atlantic dead-zone eddy, the oxygen content right below the mixed layer (50m depth) was approximately 0 µmol/kg, while the 60µmol/kg dead-zone threshold was reached at about 200m depth, resulting in a dead-zone 150m deep. It was found that mobile marine organisms are unable to follow their diurnal vertical migration and are trapped in the mixed layer, above the dead-zone, instead. Our data suggest that most severe low-oxygen ocean conditions (~0 µmol/L) are created just below the surface mixed layer in anti-cyclonic Mode Water type eddies, but still significant (~15 µmol/L) concentrations were observed in a cyclonic eddy.
Dead zones enhance key fisheries species by providing predation refuge.
Altieri, Andrew H
2008-10-01
Natural stress gradients can reduce predation intensity and increase prey abundances. Whether the harsh conditions of anthropogenic habitat degradation can similarly reduce predation intensity and structure community dynamics remains largely unexplored. Oxygen depletion in coastal waters (hypoxia) is a form of degradation that has recently emerged as one of the greatest threats to coastal ecosystems worldwide due to increased rates of eutrophication and climate change. I conducted field experiments and surveys to test whether relaxed predation could explain the paradoxically high abundance of clams that have sustained a fishery in a degraded estuary with chronic hypoxic conditions. Hypoxia reduced predation on all experimental species but enhanced the long-term survivorship of only sufficiently hypoxia-tolerant prey due to periodic extreme conditions. As a consequence, only the harvested quahog clam (Mercenaria mercenaria) thrived in hypoxic areas that were otherwise rendered dead zones with depauperate diversity and low abundances of other species. This suggests that enhanced populations of some key species may be part of a predictable nonlinear community response that sustains ecosystem services and masks overall downward trends of habitat degradation. PMID:18959318
MODELING MAGNETOROTATIONAL TURBULENCE IN PROTOPLANETARY DISKS WITH DEAD ZONES
Okuzumi, Satoshi; Hirose, Shigenobu
2011-12-01
Turbulence driven by magnetorotational instability (MRI) crucially affects the evolution of solid bodies in protoplanetary disks. On the other hand, small dust particles stabilize MRI by capturing ionized gas particles needed for the coupling of the gas and magnetic fields. To provide an empirical basis for modeling the coevolution of dust and MRI, we perform three-dimensional, ohmic-resistive MHD simulations of a vertically stratified shearing box with an MRI-inactive 'dead zone' of various sizes and with a net vertical magnetic flux of various strengths. We find that the vertical structure of turbulence is well characterized by the vertical magnetic flux and three critical heights derived from the linear analysis of MRI in a stratified disk. In particular, the turbulent structure depends on the resistivity profile only through the critical heights and is insensitive to the details of the resistivity profile. We discover scaling relations between the amplitudes of various turbulent quantities (velocity dispersion, density fluctuation, vertical diffusion coefficient, and outflow mass flux) and vertically integrated accretion stresses. We also obtain empirical formulae for the integrated accretion stresses as a function of the vertical magnetic flux and the critical heights. These empirical relations allow us to predict the vertical turbulent structure of a protoplanetary disk for a given strength of the magnetic flux and a given resistivity profile.
Analysis of archaeal communities in Gulf of Mexico dead zone sediments.
Sediments may contribute significantly to Louisiana continental shelf “dead zone” hypoxia but limited information hinders comparison of sediment biogeochemistry between norm-oxic and hypoxic seasons. Dead zone sediment cores collected during hypoxia (September 2006) had higher l...
Restoring Life to the Dead Zone: Addressing Gulf Hypoxia, a National Problem
U.S. Geological Survey
2000-01-01
The hypoxic zone in the Gulf of Mexico, the so-called 'dead zone' lacking enough oxygen to support most marine life, is one of the largest environmental issues of the decade. Practical solutions, based on sound science, are needed.
Real-time dynamic simulation of angular velocity and suppression of dead zone in IFOG
NASA Astrophysics Data System (ADS)
Gu, Hong; Huan, Yunpeng; Wang, Ansu; Luan, Jinwen
2015-02-01
The mechanism of the dead zone in closed-loop IFOGs is theoretically analyzed. The digital closed-loop transfer model for fiber-optic gyroscope is modified and the electronic cross-coupling interference from the feedback channel in the output signal of the detector in the forward channel is added. Using this model, we achieve a real-time dynamic simulation for the dead zone. The simulation results are basically consistent with experiment results. Before the electronic cross-coupling has been eliminated, the dead zone of the high-precision fiber-optic gyro is about , while the medium precision fiber-optic gyro is . After modifying the modulator circuit boards, the dead zone and the noise of the IFOG are of the same order of magnitude, better than and , respectively.
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
Ben-Kish, A; Romalis, M V
2010-11-01
Atomic magnetometers have very high absolute precision and sensitivity to magnetic fields but suffer from a fundamental problem: the vectorial or tensorial interaction of light with atoms leads to "dead zones," certain orientations of the magnetic field where the magnetometer loses its sensitivity. We demonstrate a simple polarization modulation scheme that simultaneously creates coherent population trapping (CPT) in orientation and alignment, thereby eliminating dead zones. Using 87Rb in a 10 Torr buffer gas cell we measure narrow, high-contrast CPT transparency peaks for all orientations and also show the absence of systematic effects associated with nonlinear Zeeman splitting. PMID:21231166
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
"Bagels Anyone?": Pedagogy of the Confused and Hungry in the Dead Zone.
ERIC Educational Resources Information Center
Katz, Julie
One instructor's "dead zone" (her windowless classroom in the depths of the Humanities building) was the place where little exchange between teacher and students took place. When one day she overheard the students talking about how little money they had left on their meal cards, she took a few dozen bagels to that afternoon's writing lab. The…
NASA Technical Reports Server (NTRS)
Wyman, D.; Steinman, R. M.
1973-01-01
Recently Timberlake, Wyman, Skavenski, and Steinman (1972) concluded in a study of the oculomotor error signal in the fovea that 'the oculomotor dead zone is surely smaller than 10 min and may even be less than 5 min (smaller than the 0.25 to 0.5 deg dead zone reported by Rashbass (1961) with similar stimulus conditions).' The Timberlake et al. speculation is confirmed by demonstrating that the fixating eye consistently and accurately corrects target displacements as small as 3.4 min. The contact lens optical lever technique was used to study the manner in which the oculomotor system responds to small step displacements of the fixation target. Subjects did, without prior practice, use saccades to correct step displacements of the fixation target just as they correct small position errors during maintained fixation.
Still Water: Dead Zones and Collimated Ejecta from the Impact of Granular Jets
NASA Astrophysics Data System (ADS)
Ellowitz, Jake; Turlier, Hervé; Guttenberg, Nicholas; Zhang, Wendy W.; Nagel, Sidney R.
2013-10-01
When a dense granular jet hits a target, it forms a large dead zone and ejects a highly collimated conical sheet with a well-defined opening angle. Using experiments, simulations, and continuum modeling, we find that this opening angle is insensitive to the precise target shape and the dissipation mechanisms in the flow. We show that this surprising insensitivity arises because dense granular jet impact, though highly dissipative, is nonetheless controlled by the limit of perfect fluid flow.
Zombie Vortices: The Dead Zones of Protoplanetary Disks are Not Dead
NASA Astrophysics Data System (ADS)
Jiang, Chung-Hsiang; Marcus, Philip; Pei, Suyang; Barranco, Joe; Hassanzadeh, Pedram; Lecoanet, Daniel
2014-11-01
Numerical simulations, using both the anelastic and fully compressible equations of motion, show that the ``dead zones'' of protoplanetary disks (PPDs) around forming stars are unstable and filled with vortex-dominated turbulence with Mach and Rossby numbers of order 0.2 - 0.3. The dead zones are regions in which the temperature is too cool for the gas to ionize and be destabilized by instabilities associated with the magnetic field. The ``dead zones'' were thought, by most authors, to be stable to all purely-hydrodynamic instabilities because the flow has an angular momentum that increases with increasing radius in a PPD and is therefore stable by Rayleigh's theorem. However, that theorem in not applicable to stratified flows, such as those in a PPD. We summarize our simulations with emphasis on the finite-amplitude trigger of the instability and show that when the trigger is Kolmogorov noise, the Mach number of the noise that is needed to create instability is proportional to Re - 1 / 2 , where Re is the Reynolds number of the initial noise.
DEAD ZONES AS THERMAL BARRIERS TO RAPID PLANETARY MIGRATION IN PROTOPLANETARY DISKS
Hasegawa, Yasuhiro; Pudritz, Ralph E. E-mail: pudritz@physics.mcmaster.ca
2010-02-20
Planetary migration in standard models of gaseous protoplanetary disks is known to be very rapid ({approx}10{sup 5} years), jeopardizing the existence of planetary systems. We present a new mechanism for significantly slowing rapid planetary migration, discovered by means of radiative transfer calculations of the thermal structure of protoplanetary disks irradiated by their central stars. Rapid dust settling in a disk's dead zone-a region with very little turbulence-leaves a dusty wall at its outer edge. We show that the back-heating of the dead zone by this irradiated wall produces a positive gradient of the disk temperature, which acts as a thermal barrier to planetary migration which persists for the disk lifetime. Although we analyze in detail the migration of a super-Earth in a low-mass disk around an M star, our findings can apply to a wide variety of young planetary systems. We compare our findings with other potentially important stopping mechanisms and show that there are large parameter spaces for which dead zones are likely to play the most important role for reproducing the observed mass-period relation in longer planetary periods.
Okuzumi, Satoshi; Hirose, Shigenobu
2012-07-01
Turbulence driven by magnetorotational instability (MRI) affects planetesimal formation by inducing diffusion and collisional fragmentation of dust particles. We examine conditions preferred for planetesimal formation in MRI-inactive 'dead zones' using an analytic dead-zone model based on our recent resistive MHD simulations. We argue that successful planetesimal formation requires not only a sufficiently large dead zone (which can be produced by tiny dust grains) but also a sufficiently small net vertical magnetic flux (NVF). Although often ignored, the latter condition is indeed important since the NVF strength determines the saturation level of turbulence in MRI-active layers. We show that direct collisional formation of icy planetesimal across the fragmentation barrier is possible when the NVF strength is lower than 10 mG (for the minimum-mass solar nebula model). Formation of rocky planetesimals via the secular gravitational instability is also possible within a similar range of the NVF strength. Our results indicate that the fate of planet formation largely depends on how the NVF is radially transported in the initial disk formation and subsequent disk accretion processes.
Planet filtering at the inner edges of dead zones in protoplanetary disks
NASA Astrophysics Data System (ADS)
Faure, Julien; Nelson, Richard P.
2016-02-01
Context. The interface between the dead zone and the inner active zone in a protoplanetary disk provides a promising region where the inward migration of planets may be halted owing to the existence of strong corotation torques. Recent work has indicated that this region may be prone to supporting a vortex cycle, during which vortices form at the dead-active zone interface and migrate into the active region before being destroyed, after which a new vortex forms and the cycle repeats. Aims: The aim of this paper is to examine the interaction between migrating planets and this vortex cycle, and to determine the conditions under which planets are able to remain trapped at the dead-active zone interface. Methods: We use the magnetohydrodynamics (MHD) codes PLUTO and RAMSES to perform 2D viscous disk simulations and 3D MHD simulations of protoplanetary disks containing migrating planets. A temperature switch is used to control the effective viscosity at the dead-active zone interface. Results: We find that both low mass and non-gap forming higher mass planets are able to escape from the planet trap at the inner edge of the dead zone as a result of their interaction with the migrating vortices, whereas intermediate mass planets remain trapped for the duration of simulation run times. Conclusions: Our results indicate that the vortex cycle causes the dead zone inner edge to act as an effective and mass-dependent planet filter, allowing some planets to pass through this region and others to remain there over long timescales.
Nebular dead zone effects on the D/H ratio in chondrites and comets
NASA Astrophysics Data System (ADS)
Ali-Dib, Mohamad; Martin, R. G.; Petit, J.-M.; Mousis, O.; Vernazza, P.; Lunine, J. I.
2015-11-01
Comets and chondrites show non-monotonic behaviour of their Deuterium to Hydrogen (D/H) ratio as a function of their formation location from the Sun. This is difficult to explain with a classical protoplanetary disk model that has a decreasing temperature structure with radius from the Sun.We want to understand if a protoplanetary disc with a dead zone, a region of zero or low turbulence, can explain the measured D/H values in comets and chondrites. We use time snapshots of a vertically layered disk model with turbulent surface layers and a dead zone at the midplane. The disc has a non-monotonic temperature structure due to increased heating from self-gravity in the outer parts of the dead zone. We couple this to a D/H ratio evolution model in order to quantify the effect of such thermal profiles on D/H enrichment in the nebula.We find that the local temperature peak in the disk can explain the diversity in the D/H ratios of different chondritic families. This disk temperature profile leads to a non-monotonic D/H enrichment evolution, allowing these families to acquire their different D/H values while forming in close proximity. The formation order we infer for these families is compatible with that inferred from their water abundances. However, we find that even for very young disks, the thermal profile reversal is too close to the Sun to be relevant for comets.[1] Ali-Dib, M., Martin, R. G., Petit, J.-M., Mousis, O., Vernazza, P., and Lunine, J. I. (2015, in press A&A). arXiv:1508.00263.
Trapping solids at the inner edge of the dead zone: 3-D global MHD simulations
NASA Astrophysics Data System (ADS)
Dzyurkevich, N.; Flock, M.; Turner, N. J.; Klahr, H.; Henning, Th.
2010-06-01
Context. The poorly-ionized interior of the protoplanetary disk or “dead zone” is the location where dust coagulation processes may be most efficient. However even here, planetesimal formation may be limited by the loss of solid material through radial drift, and by collisional fragmentation of the particles. Both depend on the turbulent properties of the gas. Aims: Our aim here is to investigate the possibility that solid particles are trapped at local pressure maxima in the dynamically evolving disk. We perform the first 3-D global non-ideal magnetohydrodynamical (MHD) calculations of a section of the disk treating the turbulence driven by the magneto-rotational instability (MRI). Methods: We use the ZeusMP code with a fixed Ohmic resistivity distribution. The domain contains an inner MRI-active region near the young star and an outer midplane dead zone, with the transition between the two modeled by a sharp increase in the magnetic diffusivity. Results: The azimuthal magnetic fields generated in the active zone oscillate over time, changing sign about every 150 years. We thus observe the radial structure of the “butterfly pattern” seen previously in local shearing-box simulations. The mean magnetic field diffuses from the active zone into the dead zone, where the Reynolds stress nevertheless dominates, giving a residual α between 10-4 and 10-3. The greater total accretion stress in the active zone leads to a net reduction in the surface density, so that after 800 years an approximate steady state is reached in which a local radial maximum in the midplane pressure lies near the transition radius. We also observe the formation of density ridges within the active zone. Conclusions: The dead zone in our models possesses a mean magnetic field, significant Reynolds stresses and a steady local pressure maximum at the inner edge, where the outward migration of planetary embryos and the efficient trapping of solid material are possible.
How big is the Ocean Dead Zone off the Coast of California?
NASA Astrophysics Data System (ADS)
Hofmann, A. F.; Peltzer, E. T.; Walz, P. M.; Brewer, P. G.
2010-12-01
The term “Ocean Dead Zone”, generally referring to a zone that is devoid of aerobic marine life of value to humans, is now widely used in the press and scientific literature but it appears to be not universally defined. The global assessment and monitoring of ocean dead zones, however, is of high public concern due to the considerable economic value associated with impacted fisheries and with questions over the growth of these zones forced by climate change. We report on the existence of a zone at ~850m depth off Santa Monica, California where dissolved oxygen (DO) levels are 1 μmol/kg; an order of magnitude below any existing definition of an “Ocean Dead Zone”. ROV dives show the region to be visually devoid of all aerobic marine life. But how large is this dead zone, and how may its boundaries be defined? Without an accepted definition we cannot report this nor can we compare it to other dead zones reported elsewhere in the world. “Dead zones” are now assessed solely by DO levels. A multitude of values in different units are used (Fig 1), which are clearly not universally applicable. This seriously hampers an integrated global monitoring and management effort and frustrates the development of valid connections with climate change and assessment of the consequences. Furthermore, input of anthropogenic CO2 can also stress marine life. Recent work supported by classical data suggests that higher pCO2 influences the thermodynamic energy efficiency of oxic respiration (CH2O + O2 -> CO2 + H2O). The ratio pO2/pCO2, called the respiration index (RI), emerges as the critical variable, combining the impacts of warming on DO and rising CO2 levels within a single, well defined quantity. We advocate that future monitoring efforts report pO2 and pCO2 concurrently, thus making it possible to classify, monitor and manage “dead zones” within a standard reference system that may include, as with e.g, hurricanes, differing categories of intensity. Fig.1. A DO
Characterization and impact of "dead-zone" eddies in the tropical Northeast Atlantic Ocean
NASA Astrophysics Data System (ADS)
Schuette, Florian; Karstensen, Johannes; Krahmann, Gerd; Hauss, Helena; Fiedler, Björn; Brandt, Peter; Visbeck, Martin; Körtzinger, Arne
2016-04-01
Localized open-ocean low-oxygen dead-zones in the tropical Northeast Atlantic are recently discovered ocean features that can develop in dynamically isolated water masses within cyclonic eddies (CE) and anticyclonic modewater eddies (ACME). Analysis of a comprehensive oxygen dataset obtained from gliders, moorings, research vessels and Argo floats shows that eddies with low oxygen concentrations at 50-150 m depths can be found in surprisingly high numbers and in a large area (from about 5°N to 20°N, from the shelf at the eastern boundary to 30°W). Minimum oxygen concentrations of about 9 μmol/kg in CEs and close to anoxic concentrations (< 1 μmol/kg) in ACMEs were observed. In total, 495 profiles with oxygen concentrations below the minimum background concentration of 40 μmol/kg could be associated with 27 independent "dead-zone" eddies (10 CEs; 17 ACMEs). The low oxygen concentration right beneath the mixed layer has been attributed to the combination of high productivity in the surface waters of the eddies and the isolation of the eddies' cores. Indeed eddies of both types feature a cold sea surface temperature anomaly and enhanced chlorophyll concentrations in their center. The oxygen minimum is located in the eddy core beneath the mixed layer at around 80 m depth. The mean oxygen anomaly between 50 to 150 m depth for CEs (ACMEs) is -49 (-81) μmol/kg. Eddies south of 12°N carry weak hydrographic anomalies in their cores and seem to be generated in the open ocean away from the boundary. North of 12°N, eddies of both types carry anomalously low salinity water of South Atlantic Central Water origin from the eastern boundary upwelling region into the open ocean. This points to an eddy generation near the eastern boundary. A conservative estimate yields that around 5 dead-zone eddies (4 CEs; 1 ACME) per year entering the area north of 12°N between the Cap Verde Islands and 19°W. The associated contribution to the oxygen budget of the shallow oxygen minimum
GAUGE RUN-TO-DETONATION DATA AND FAILURE/DEAD ZONE MODELING
Souers, P C; Vitello, P; Vandersall, K S
2009-06-26
Previous shock initiation run-to-detonation experiments on energetic materials were plotted with distance and time to get a single distance/time to detonation. Modern shots utilize enough gauges so that the distance-time data can be differentiated, which shows not only the usual inflection pressure point before detonation, referred to here as P{sub b}, but also a second, low-pressure inflection, referred to here as P{sub a}, that marks rapid ramp-up of the initiation. An analysis of the TATB based LX-17 and PBX 9502 in addition to the LLM-105 based RX-55 data shows that both P{sub a} and P{sub b} increase linearly with the initiation pressure created by the flyer plate. This contradicts the current method in the Tarantula failure/dead zone model, which uses constant pressure boundaries between reaction regions. Modeling changes required by the new data will be considered.
Sweep-up growth at the inner edge of dead zones
NASA Astrophysics Data System (ADS)
Drazkowska, Joanna; Windmark, Fredrik; Dullemond, Cornelis P.
2013-07-01
Planetesimal formation is still not understood. Coagulation models have revealed numerous obstacles to the dust growth. One of them is the bouncing barrier. The growth of small dust grains was shown to be completely halted already for cm-sized silicate particles. This barrier can be actually beneficial to the growth. When a limited number of grains is inserted into a population halted by the bouncing, growth to planetesimal sizes is possible. This is because as long as a collision between two big particles generally leads to fragmentation, a collision involving non-equal sized aggregates can lead to growth via so-called fragmentation with mass transfer. The origin of the first seeds is a problem for this scenario. We propose a new method of providing the seeds. We find that a steep radial variation in the turbulence efficiency that takes place at the inner edge of a dead zone, promotes planetesimal formation via sweep-up in several ways. It provides a pressure trap that saves the dust from the radial drift barrier. It also causes a change in the maximum size of aggregates at which growth barriers occur. The seeds can grow in the dead zone, where the bouncing barrier occurs for larger grains, and then be delivered by radial mixing to the MRI active region, where the growth via sweep-up occurs. In the presented model, which is employing an ad hoc turbulent viscosity change near the snow line, it is possible to grow planetesimals by incremental growth on timescales relevant for planet formation.
NASA Astrophysics Data System (ADS)
Lesur, Geoffroy; Kunz, Matthew W.; Fromang, Sébastien
2014-06-01
Protoplanetary discs are poorly ionised due to their low temperatures and high column densities and are therefore subject to three "non-ideal" magnetohydrodynamic (MHD) effects: Ohmic dissipation, ambipolar diffusion, and the Hall effect. The existence of magnetically driven turbulence in these discs has been a central question since the discovery of the magnetorotational instability (MRI). Early models considered Ohmic diffusion only and led to a scenario of layered accretion, in which a magnetically "dead" zone in the disc midplane is embedded within magnetically "active" surface layers at distances of about 1-10 au from the central protostellar object. Recent work has suggested that a combination of Ohmic dissipation and ambipolar diffusion can render both the midplane and surface layers of the disc inactive and that torques due to magnetically driven outflows are required to explain the observed accretion rates. We reassess this picture by performing three-dimensional numerical simulations that include all three non-ideal MHD effects for the first time. We find that the Hall effect can generically "revive" dead zones by producing a dominant azimuthal magnetic field and a large-scale Maxwell stress throughout the midplane, provided that the angular velocity and magnetic field satisfy Ω·B > 0. The attendant large magnetic pressure modifies the vertical density profile and substantially increases the disc scale height beyond its hydrostatic value. Outflows are produced but are not necessary to explain accretion rates ≲ 10-7 M⊙ yr-1. The flow in the disc midplane is essentially laminar, suggesting that dust sedimentation may be efficient. These results demonstrate that if the MRI is relevant for driving mass accretion in protoplanetary discs, one must include the Hall effect to obtain even qualitatively correct results. Appendices are available in electronic form at http://www.aanda.org
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.
A multidisciplinary glider survey of an open ocean dead-zone eddy
NASA Astrophysics Data System (ADS)
Karstensen, Johannes; Schütte, Florian; Pietri, Alice; Krahmann, Gerd; Fiedler, Björn; Löscher, Carolin; Grundle, Damian; Hauss, Helena; Körtzinger, Arne; Testor, Pierre; Viera, Nuno
2016-04-01
The physical (temperature, salinity) and biogeochemical (oxygen, nitrate, chlorophyll fluorescence, turbidity) structure of an anticyclonic modewater eddy, hosting an open ocean dead zone, is investigated using observational data sampled in high temporal and spatial resolution with autonomous gliders in March and April 2014. The core of the eddy is identified in the glider data as a volume of fresher (on isopycnals) water in the depth range from the mixed layer base (about 70m) to about 200m depth. The width is about 80km. The core aligns well with the 40 μmolkg‑1 oxygen contour. From two surveys about 1 month apart, changes in the minimal oxygen concentrations (below 5μmolkg‑1) are observed that indicate that small scale processes are in operation. Several scales of coherent variability of physical and biogeochemical variable are identified - from a few meters to the mesoscale. One of the gliders carried an autonomous Nitrate (N) sensor and the data is used to analyse the possible nitrogen pathways within the eddy. Also the highest N is accompanied by lowest oxygen concentrations, the AOU:N ratio reveals a preferred oxygen cycling per N.
NASA Astrophysics Data System (ADS)
Miyasato, Yoshihiko
The problem of constructing model reference adaptive H∞ control for distributed parameter systems of hyperbolic type preceded by unknown input nonlinearity such as dead zone or backlash, is considered in this paper. Distributed parameter systems are infinite dimensional processes, but the proposed control scheme is constructed from finite dimensional controllers. An adaptive inverse model is introduced to estimate and compensate the input nonlinearity. The stabilizing control signal is added to regulate the effect of spill-over terms, and it is derived as a solution of certain H∞ control problem where the residual part of the inverse model and the spill-over term are considered as external disturbances to the process.
NASA Astrophysics Data System (ADS)
Marcus, Philip S.; Pei, Suyang; Jiang, Chung-Hsiang; Barranco, Joseph A.; Hassanzadeh, Pedram; Lecoanet, Daniel
2015-07-01
There is considerable interest in hydrodynamic instabilities in dead zones of protoplanetary disks as a mechanism for driving angular momentum transport and as a source of particle-trapping vortices to mix chondrules and incubate planetesimal formation. We present simulations with a pseudo-spectral anelastic code and with the compressible code Athena, showing that stably stratified flows in a shearing, rotating box are violently unstable and produce space-filling, sustained turbulence dominated by large vortices with Rossby numbers of order ˜0.2-0.3. This Zombie Vortex Instability (ZVI) is observed in both codes and is triggered by Kolmogorov turbulence with Mach numbers less than ˜0.01. It is a common view that if a given constant density flow is stable, then stable vertical stratification should make the flow even more stable. Yet, we show that sufficient vertical stratification can be unstable to ZVI. ZVI is robust and requires no special tuning of boundary conditions, or initial radial entropy or vortensity gradients (though we have studied ZVI only in the limit of infinite cooling time). The resolution of this paradox is that stable stratification allows for a new avenue to instability: baroclinic critical layers. ZVI has not been seen in previous studies of flows in rotating, shearing boxes because those calculations frequently lacked vertical density stratification and/or sufficient numerical resolution. Although we do not expect appreciable angular momentum transport from ZVI in the small domains in this study, we hypothesize that ZVI in larger domains with compressible equations may lead to angular transport via spiral density waves.
Rossby wave instability and long-term evolution of dead zones in protoplanetary discs
NASA Astrophysics Data System (ADS)
Miranda, Ryan; Lai, Dong; Méheut, Héloïse
2016-04-01
The physical mechanism of angular momentum transport in poorly ionized regions of protoplanetary discs, the dead zones (DZs), is not understood. The presence of a DZ naturally leads to conditions susceptible to the Rossby wave instability (RWI), which produces vortices and spiral density waves that may revive the DZ and be responsible for observed large-scale disc structures. We present a series of two-dimensional hydrodynamic simulations to investigate the role of the RWI in DZs, including its impact on the long-term evolution of the disc and its morphology. The non-linear RWI can generate Reynolds stresses (effective α parameter) as large as 0.01-0.05 in the DZ, helping to sustain quasi-steady accretion throughout the disc. It also produces novel disc morphologies, including azimuthal asymmetries with m = 1, 2, and atypical vortex shapes. The angular momentum transport strength and morphology are most sensitive to two parameters: the radial extent of the DZ and the disc viscosity. The largest Reynolds stresses are produced when the radial extent of the DZ is less than its distance to the central star. Such narrow DZs lead to a single vortex or two coherent antipodal vortices in the quasi-steady state. The edges of wider DZs evolve separately, resulting in two independent vortices and reduced angular momentum transport efficiency. In either case, we find that, because of the Reynolds stresses generated by the non-linear RWI, gravitational instability is unlikely to play a role in angular momentum transport across the DZ, unless the accretion rate is sufficiently high.
Thomas, Peter; Rahman, Md Saydur
2010-01-01
Recently evidence has been obtained for reproductive impairment in estuarine populations of Atlantic croaker exposed to seasonal hypoxia. However, it is not known whether a similar disruption of reproductive function occurs in croaker inhabiting a much larger hypoxic area, the extensive dead zone in coastal regions of the northern Gulf of Mexico extending from the Mississippi Delta to East Texas. Gonadal development in male Atlantic croaker collected in September 2008 at six sites in the dead zone was compared to that in male fish sampled from three reference sites east of the Mississippi Delta which do not experience persistent hypoxia. Croaker testes collected from the dead zone were at an earlier stage of spermatogenesis than those from the reference sites. Histological examination of the testes collected from the dead zone showed that their tubules had small lumens that contained very little sperm compared to the lumens of the reference fish. Overall, sperm production was 26.2% that of the control fish at the reference sites. This decrease in spermatogenesis at the dead zone sites was accompanied by an approximately 50% decrease in testicular growth compared to that in the reference fish. The results suggest that reproductive impairment can occur over regional scales in marine fish populations exposed to extensive seasonal hypoxia in dead zones with potential long-term impacts on population abundance. PMID:19931178
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
Lyra, Wladimir; Mac Low, Mordecai-Mark E-mail: mordecai@amnh.org
2012-09-01
It has been suggested that the transition between magnetorotationally active and dead zones in protoplanetary disks should be prone to the excitation of vortices via Rossby wave instability (RWI). However, the only numerical evidence for this has come from alpha disk models, where the magnetic field evolution is not followed, and the effect of turbulence is parameterized by Laplacian viscosity. We aim to establish the phenomenology of the flow in the transition in three-dimensional resistive-magnetohydrodynamical models. We model the transition by a sharp jump in resistivity, as expected in the inner dead zone boundary, using the PENCIL CODE to simulate the flow. We find that vortices are readily excited in the dead side of the transition. We measure the mass accretion rate finding similar levels of Reynolds stress at the dead and active zones, at the {alpha} Almost-Equal-To 10{sup -2} level. The vortex sits in a pressure maximum and does not migrate, surviving until the end of the simulation. A pressure maximum in the active zone also triggers the RWI. The magnetized vortex that results should be disrupted by parasitical magneto-elliptic instabilities, yet it subsists in high resolution. This suggests that either the parasitic modes are still numerically damped or that the RWI supplies vorticity faster than they can destroy it. We conclude that the resistive transition between the active and dead zones in the inner regions of protoplanetary disks, if sharp enough, can indeed excite vortices via RWI. Our results lend credence to previous works that relied on the alpha-disk approximation, and caution against the use of overly reduced azimuthal coverage on modeling this transition.
Charnoz, Sebastien; Taillifet, Esther
2012-07-10
Dust is a major component of protoplanetary and debris disks as it is the main observable signature of planetary formation. However, since dust dynamics are size-dependent (because of gas drag or radiation pressure) any attempt to understand the full dynamical evolution of circumstellar dusty disks that neglect the coupling of collisional evolution with dynamical evolution is thwarted because of the feedback between these two processes. Here, a new hybrid Lagrangian/Eulerian code is presented that overcomes some of these difficulties. The particles representing 'dust clouds' are tracked individually in a Lagrangian way. This system is then mapped on an Eulerian spatial grid, inside the cells of which the local collisional evolutions are computed. Finally, the system is remapped back in a collection of discrete Lagrangian particles, keeping their number constant. An application example of dust growth in a turbulent protoplanetary disk at 1 AU is presented. First, the growth of dust is considered in the absence of a dead zone and the vertical distribution of dust is self-consistently computed. It is found that the mass is rapidly dominated by particles about a fraction of a millimeter in size. Then the same case with an embedded dead zone is investigated and it is found that coagulation is much more efficient and produces, in a short timescale, 1-10 cm dust pebbles that dominate the mass. These pebbles may then be accumulated into embryo-sized objects inside large-scale turbulent structures as shown recently.
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
In search of the dead zone: Use of otoliths for tracking fish exposure to hypoxia
Limburg, Karin E.; Walther, Benjamin D.; Lu, Zunli; Jackman, George; Mohan, John; Walther, Yvonne; Nissling, Anders; Weber, Peter K.; Schmitt, Axel K.
2015-01-01
Otolith chemistry is often useful for tracking provenance of fishes, as well as examining migration histories. Whereas elements such as strontium and barium correlate well with salinity and temperature, experiments that examine manganese uptake as a function of these parameters have found no such correlation. Instead, dissolved manganese is available as a redox product, and as such, is indicative of low-oxygen conditions. Here we present evidence for that mechanism in a range of habitats from marine to freshwater, across species, and also present ancillary proxies that support the mechanism as well. For example, iodine is redox-sensitive and varies inversely with Mn; and sulfur stable isotope ratios provide evidence of anoxic sulfate reduction in some circumstances.
In search of the dead zone: Use of otoliths for tracking fish exposure to hypoxia
Limburg, Karin E.; Walther, Benjamin D.; Lu, Zunli; Jackman, George; Mohan, John; Walther, Yvonne; Nissling, Anders; Weber, Peter K.; Schmitt, Axel K.
2015-01-01
Otolith chemistry is often useful for tracking provenance of fishes, as well as examining migration histories. Whereas elements such as strontium and barium correlate well with salinity and temperature, experiments that examine manganese uptake as a function of these parameters have found no such correlation. Instead, dissolved manganese is available as a redox product, and as such, is indicative of low-oxygen conditions. Here we present evidence for that mechanism in a range of habitats from marine to freshwater, across species, and also present ancillary proxies that support the mechanism as well. For example, iodine is redox-sensitive and varies inversely withmore » Mn; and sulfur stable isotope ratios provide evidence of anoxic sulfate reduction in some circumstances.« less
In search of the dead zone: Use of otoliths for tracking fish exposure to hypoxia
NASA Astrophysics Data System (ADS)
Limburg, Karin E.; Walther, Benjamin D.; Lu, Zunli; Jackman, George; Mohan, John; Walther, Yvonne; Nissling, Anders; Weber, Peter K.; Schmitt, Axel K.
2015-01-01
Otolith chemistry is often useful for tracking provenance of fishes, as well as examining migration histories. Whereas elements such as strontium and barium correlate well with salinity and temperature, experiments that examine manganese uptake as a function of these parameters have found no such correlation. Instead, dissolved manganese is available as a redox product, and as such, is indicative of low-oxygen conditions. Here we present evidence for that mechanism in a range of habitats from marine to freshwater, across species, and also present ancillary proxies that support the mechanism as well. For example, iodine is redox-sensitive and varies inversely with Mn; and sulfur stable isotope ratios provide evidence of anoxic sulfate reduction in some circumstances. Further, S may be incorporated trophically whereas other elements appear to be taken up directly from water. This research suggests a potential means to identify individual fish exposure to hypoxia, over entire lifetimes. With further testing and understanding, in the future fish may be able to be used as "mobile monitors" of hypoxic conditions.
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
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.
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.
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…
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.
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
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.
Circadian clocks optimally adapt to sunlight for reliable synchronization
Hasegawa, Yoshihiko; Arita, Masanori
2014-01-01
Circadian oscillation provides selection advantages through synchronization to the daylight cycle. However, a reliable clock must be designed through two conflicting properties: entrainability to synchronize internal time with periodic stimuli such as sunlight, and regularity to oscillate with a precise period. These two aspects do not easily coexist, because better entrainability favours higher sensitivity which may sacrifice regularity. To investigate conditions for satisfying the two properties, we analytically calculated the optimal phase–response curve with a variational method. Our results indicate an existence of a dead zone, i.e. a time period during which input stimuli neither advance nor delay the clock. A dead zone appears only when input stimuli obey the time course of actual solar radiation, but a simple sine curve cannot yield a dead zone. Our calculation demonstrates that every circadian clock with a dead zone is optimally adapted to the daylight cycle. PMID:24352677
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).
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.
Gressel, O.; Nelson, R. P.; Turner, N. J.; Ziegler, U. E-mail: r.p.nelson@qmul.ac.uk E-mail: uziegler@aip.de
2013-12-10
We present global hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations with mesh refinement of accreting planets embedded in protoplanetary disks (PPDs). The magnetized disk includes Ohmic resistivity that depends on the overlying mass column, leading to turbulent surface layers and a dead zone near the midplane. The main results are: (1) the accretion flow in the Hill sphere is intrinsically three-dimensional for HD and MHD models. Net inflow toward the planet is dominated by high-latitude flows. A circumplanetary disk (CPD) forms. Its midplane flows outward in a pattern whose details differ between models. (2) The opening of a gap magnetically couples and ignites the dead zone near the planet, leading to stochastic accretion, a quasi-turbulent flow in the Hill sphere, and a CPD whose structure displays high levels of variability. (3) Advection of magnetized gas onto the rotating CPD generates helical fields that launch magnetocentrifugally driven outflows. During one specific epoch, a highly collimated, one-sided jet is observed. (4) The CPD's surface density is ∼30 g cm{sup −2}, small enough for significant ionization and turbulence to develop. (5) The accretion rate onto the planet in the MHD simulation reaches a steady value 8 × 10{sup –3} M {sub ⊕} yr{sup –1} and is similar in the viscous HD runs. Our results suggest that gas accretion onto a forming giant planet within a magnetized PPD with a dead zone allows rapid growth from Saturnian to Jovian masses. As well as being relevant for giant planet formation, these results have important implications for the formation of regular satellites around gas giant planets.
NASA Astrophysics Data System (ADS)
Carling, Paul; Kleinhans, Maarten; Leyland, Julian; Besozzi, Louison; Duranton, Pierre; Trieu, Hai; Teske, Roy
2014-05-01
Understanding of flow resistance of forested floodplains is essential for floodplain flow routing and floodplain reforestation projects. Although the flow resistance of grass-lined channels is well-known, flow retention due to flow-blocking by trees is poorly understood. Flow behaviour through tree-filled channels or over forested floodplain surfaces has largely been addressed using laboratory studies of artificial surfaces and vegetation. Herein we take advantage of a broad, shallow earthen experimental outdoor channel with headwater and tailwater controls. The channel was disused and left undisturbed for more than 20 years. During this time period, small deciduous trees and a soil cover of grass, herbs and leaf-litter established naturally. We measured flow resistance and fluid retention in fifteen controlled water discharge experiments for the following conditions: (a) natural cover of herbs and trees; (b) trees only and; (c) earthen channel only. In the b-experiments the herbaceous groundcover was first removed carefully and in the c-experiments the trees were first cut flush with the earthen channel floor. Rhodamine-B dye was used to tag the flow and the resultant fluorescence of water samples were systematically assayed through time at two stations along the length of the channel. Dilution-curve data were analysed within the Aggregated Dead Zone (ADZ) framework to yield bulk flow parameters including dispersion, fluid retention and flow resistance parameters after the procedure of Richardson & Carling (2006). The primary response of the bulk flow to vegetation removal was an increase in bulk velocity, with depth and wetted width decreasing imperceptibly at the resolution of measurement. An overall reduction in flow resistance and retention occurred as discharge increased in all experiments and flow retention. Retentiveness was more prominent during low flow and for all three experimental conditions tended to converge on a constant low value for high
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.
Patel, Ruben; Pedersen, Geir; Ona, Egil
2009-02-01
Acoustic measurement of near-bottom fish with a directional transducer is generally problematical because the powerful bottom echo interferes with weaker echoes from fish within the main lobe but at greater ranges than that of the bottom. The volume that is obscured is called the dead zone. This has already been estimated for the special case of a flat horizontal bottom when observed by an echo sounder with a stable vertical transducer beam [Ona, E., and Mitson, R. B. (1996). ICES J. Mar. Sci. 53, 677-690]. The more general case of observation by a split-beam echo sounder with a transducer mounted on a noninertial platform is addressed here. This exploits the capability of a split-beam echo sounder to measure the bottom slope relative to the beam axis and thence to allow the dead-zone volume over a flat but sloping bottom to be estimated analytically. The method is established for the Simrad EK60 scientific echo sounder, with split-beam transducers operating at 18, 38, 70, 120, and 200 kHz. It is validated by comparing their estimates of seafloor slope near the Lofoten Islands, N67-70, with simultaneous measurements made by two hydrographic multibeam sonars, the Simrad EM100295 kHz and EM30030 kHz systems working in tandem. PMID:19206847
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
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.
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
Performance-Based Adaptive Fuzzy Tracking Control for Networked Industrial Processes.
Wang, Tong; Qiu, Jianbin; Yin, Shen; Gao, Huijun; Fan, Jialu; Chai, Tianyou
2016-08-01
In this paper, the performance-based control design problem for double-layer networked industrial processes is investigated. At the device layer, the prescribed performance functions are first given to describe the output tracking performance, and then by using backstepping technique, new adaptive fuzzy controllers are designed to guarantee the tracking performance under the effects of input dead-zone and the constraint of prescribed tracking performance functions. At operation layer, by considering the stochastic disturbance, actual index value, target index value, and index prediction simultaneously, an adaptive inverse optimal controller in discrete-time form is designed to optimize the overall performance and stabilize the overall nonlinear system. Finally, a simulation example of continuous stirred tank reactor system is presented to show the effectiveness of the proposed control method. PMID:27168605
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.
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
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.
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.
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
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
Wu, T; Peng, X; Lin, Z; Guo, H
2015-10-01
We demonstrate an all-optical (4)He atomic magnetometer experimental scheme based on an original Bell-Bloom configuration. A single intensity-modulated linearly polarized laser beam is used both for generating spin polarization within a single (4)He vapor and probing the spin precessing under a static magnetic field. The transmitted light signal from the vapor is then phase-sensitively detected at the modulation frequency and its harmonics, which lead to the atomic magnetic resonance signals. Based on this structure, a liquid crystal is added in our magnetometer system and constitutes a polarization rotator. By controlling the voltage applied on the liquid crystal, the light linear polarization vector can be kept perpendicular with the ambient magnetic field direction, which in turn provides the maximum resonance signal amplitude. Moreover, the system exhibits a magnetic-field noise floor of about 2pT/√Hz, which is not degraded due to the presence of the liquid crystal and varying magnetic field direction. The experiment results prove that our method can eliminate the dead-zone effect, improve the system spatial isotropy, and thus be suitable in mobile applications. PMID:26520938
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
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
NASA Astrophysics Data System (ADS)
Briggs, M.; Gooseff, M.; Morkeski, K.; Wollheim, W.; Hopkinson, C.; Peterson, B.; Vorosmarty, C.
2007-12-01
A major enhancement to our understanding of how watersheds function would be the ability to discriminate between in-channel dead zone ( DZ) and hyporheic zone ( HZ) transient storage, and an evaluation of how these properties scale across stream orders. The nature of DZ storage is to display faster exchange rates with the main channel and less overall sediment contact time than HZ storage. These differences have great significance to many in-stream processes such as nutrient cycling. The combination of high slope, coarse bed material and fluvial structure endemic to many 1st order streams can provide greater forcing of hyporheic flow paths than occurs within the lower gradient 5th order streams. Conversely many 5th order reaches exhibit large side pool and back eddy DZ areas not common along 1st order streams. This study builds on existing methods to delineate the DZ and HZ from the integrated signal of a conservative solute's breakthrough curve ( BTC). Data for this comparison were collected over the summer of 2007 within the Ipswich River watershed, a basin which drains into Plum Island Sound on the north shore of Massachusetts, USA. The conservative solute NaCl was injected into both a 1st order medium gradient stream and a 5th order low gradient stream. The BTCs collected in thalwegs from the NaCl injections were simulated using a version of the solute transport model OTIS containing two zones of transient storage. Hydrometric measurements of stream velocity were used to estimate average main channel cross sectional area ( A) and DZ cross sectional area ( ASDZ) for each reach to constrain parameter estimates and avoid model equifinality between the storage zones. Initial values for the exchange rate between main channel flow and DZ storage ( αDZ) were estimated from DZ BTCs. Our results indicate that although the overall storage zone is much larger in proportion to the main channel for the 1st order reach than for the 5th order reach, the percentage of median
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
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
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.
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.
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.
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.
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
Robust adaptive regulation without persistent excitation
NASA Technical Reports Server (NTRS)
Lozano-Leal, Rogelio
1988-01-01
A globally convergent adaptive regulator for minimum or nonminimum phase systems subject to bounded distrubances and unmodeled dynamics is presented. The control strategy is designed for a particular input-output representation obtained from the state space representation of the system. The leading coefficient of the new representation is the product of the observability and controllability matrices of the system. The controller scheme uses a Least Squares identification algorithm with a dead zone. The dead zone is chosen to obtain convergence properties on the estimates and on the covariance matrix as well. This allows the definition of modified estimates which secure well-conditioned matrices in the adaptive control law. Explicit bounds on the plant output are given.
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…
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
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.
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.
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
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
NASA Astrophysics Data System (ADS)
Key, K.
2013-12-01
This work announces the public release of an open-source inversion code named MARE2DEM (Modeling with Adaptively Refined Elements for 2D Electromagnetics). Although initially designed for the rapid inversion of marine electromagnetic data, MARE2DEM now supports a wide variety of acquisition configurations for both offshore and onshore surveys that utilize electric and magnetic dipole transmitters or magnetotelluric plane waves. The model domain is flexibly parameterized using a grid of arbitrarily shaped polygonal regions, allowing for complicated structures such as topography or seismically imaged horizons to be easily assimilated. MARE2DEM efficiently solves the forward problem in parallel by dividing the input data parameters into smaller subsets using a parallel data decomposition algorithm. The data subsets are then solved in parallel using an automatic adaptive finite element method that iterative solves the forward problem on successively refined finite element meshes until a specified accuracy tolerance is met, thus freeing the end user from the burden of designing an accurate numerical modeling grid. Regularized non-linear inversion for isotropic or anisotropic conductivity is accomplished with a new implementation of Occam's method referred to as fast-Occam, which is able to minimize the objective function in much fewer forward evaluations than the required by the original method. This presentation will review the theoretical considerations behind MARE2DEM and use a few recent offshore EM data sets to demonstrate its capabilities and to showcase the software interface tools that streamline model building and data inversion.
NASA Technical Reports Server (NTRS)
Ianculescu, G. D.; Klop, J. J.
1992-01-01
Classical and adaptive control algorithms for the solar array pointing system of the Space Station Freedom are designed using a continuous rigid body model of the solar array gimbal assembly containing both linear and nonlinear dynamics due to various friction components. The robustness of the design solution is examined by performing a series of sensitivity analysis studies. Adaptive control strategies are examined in order to compensate for the unfavorable effect of static nonlinearities, such as dead-zone uncertainties.
Li, Zhijun; Su, Chun-Yi
2013-09-01
In this paper, adaptive neural network control is investigated for single-master-multiple-slaves teleoperation in consideration of time delays and input dead-zone uncertainties for multiple mobile manipulators carrying a common object in a cooperative manner. Firstly, concise dynamics of teleoperation systems consisting of a single master robot, multiple coordinated slave robots, and the object are developed in the task space. To handle asymmetric time-varying delays in communication channels and unknown asymmetric input dead zones, the nonlinear dynamics of the teleoperation system are transformed into two subsystems through feedback linearization: local master or slave dynamics including the unknown input dead zones and delayed dynamics for the purpose of synchronization. Then, a model reference neural network control strategy based on linear matrix inequalities (LMI) and adaptive techniques is proposed. The developed control approach ensures that the defined tracking errors converge to zero whereas the coordination internal force errors remain bounded and can be made arbitrarily small. Throughout this paper, stability analysis is performed via explicit Lyapunov techniques under specific LMI conditions. The proposed adaptive neural network control scheme is robust against motion disturbances, parametric uncertainties, time-varying delays, and input dead zones, which is validated by simulation studies. PMID:24808577
Adaptive output feedback control of flexible systems
NASA Astrophysics Data System (ADS)
Yang, Bong-Jun
Neural network-based adaptive output feedback approaches that augment a linear control design are described in this thesis, and emphasis is placed on their real-time implementation with flexible systems. Two different control architectures that are robust to parametric uncertainties and unmodelled dynamics are presented. The unmodelled effects can consist of minimum phase internal dynamics of the system together with external disturbance process. Within this context, adaptive compensation for external disturbances is addressed. In the first approach, internal model-following control, adaptive elements are designed using feedback inversion. The effect of an actuator limit is treated using control hedging, and the effect of other actuation nonlinearities, such as dead zone and backlash, is mitigated by a disturbance observer-based control design. The effectiveness of the approach is illustrated through simulation and experimental testing with a three-disk torsional system, which is subjected to control voltage limit and stiction. While the internal model-following control is limited to minimum phase systems, the second approach, external model-following control, does not involve feedback linearization and can be applied to non-minimum phase systems. The unstable zero dynamics are assumed to have been modelled in the design of the existing linear controller. The laboratory tests for this method include a three-disk torsional pendulum, an inverted pendulum, and a flexible-base robot manipulator. The external model-following control architecture is further extended in three ways. The first extension is an approach for control of multivariable nonlinear systems. The second extension is a decentralized adaptive control approach for large-scale interconnected systems. The third extension is to make use of an adaptive observer to augment a linear observer-based controller. In this extension, augmenting terms for the adaptive observer can be used to achieve adaptation in
NASA Astrophysics Data System (ADS)
Marcus, P. S.; Jiang, C.; Pei, S.; Hassanzadeh, P.
2012-12-01
The annular region of a protoplanetary disk, approximately 5 - 12 AU from the protostar and within 1 or 2 pressure scale heights of the mid-plane, has too low a temperature to significantly ionize hydrogen gas and therefore to destabilize the near-Keplerian flow via the magneto-hydrodynamic instability (MRI). Because it is assumed that Keplerian flow is linearly stable and that no hydrodynamic features can arise spontaneously within this region and thereby aid in star formation by transporting angular momentum and energy radially, this region is known as the "dead zone". Here we show that the pronouncement of this region as "dead" is premature (to paraphrase Mark Twain). Baroclinic critical layers can occur in rotating, vertically-stratified, uni-directional shear flows such as that in a protoplanetary disk. They are special cases of neutrally stable eigenmodes. Baroclinic critical layers have logarithmic singularities in density and vertical velocity. They differ from barotropic critical layers associated with Kelvin's cats-eyes in constant-density, unidirectional shear flows, which form at locations where the shear flow velocity matches the eigenmode's phase speed and have singularities only in stream-wise velocities. Baroclinic critical layers are easily excited with no special tuning of parameters by perturbations from vortices or waves. Unlike barotropic critical layers the amplitudes of baroclinic layers become large by drawing energy from the background shear. In the case of protoplanetary disk, energy is extracted from the huge reservoir of kinetic energy in the near-Keplerian shear. The large vertical velocities in the critical layers, coupled with the Coriolis parameter create large-amplitude vortex layers. These layers often roll-up into large coherent vortices. The baroclinic critical layers' growth and roll-up are robust: they form in cylindrical and Cartesian geometries, in Boussinesq fluids and ideal gases, and in flows with uniform and non
Bolea, Mario; Mora, José; Ortega, Beatriz; Capmany, José
2009-03-30
We propose theoretically and demonstrate experimentally an optical architecture for flexible Ultra-Wideband pulse generation. It is based on an N-tap reconfigurable microwave photonic filter fed by a laser array by using phase inversion in a Mach-Zehnder modulator. Since a large number of positive and negative coefficients can be easily implemented, UWB pulses fitted to the FCC mask requirements can be generated. As an example, a four tap pulse generator is experimentally demonstrated which complies with the FCC regulation. The proposed pulse generator allows different pulse modulation formats since the amplitude, polarity and time delay of generated pulse is controlled. PMID:19333263
Trinh, Quoclinh; Zhu, Pengyu; Shi, Hui; Xu, Wentao; Hao, Junran; Luo, Yunbo; Huang, Kunlun
2014-12-01
The polymerase chain reaction (PCR)-based genome walking method has been extensively used to isolate unknown flanking sequences, whereas nonspecific products are always inevitable. To resolve these problems, we developed a new strategy to isolate the unknown flanking sequences by combining A-T linker adapter PCR with inverse PCR (I-PCR) or thermal asymmetric interlaced PCR (TAIL-PCR). The result showed that this method can be efficiently achieved with the flanking sequence from the Arabidopsis mutant and papain gene. Our study provides researchers with an additional method for determining genomic DNA flanking sequences to identify the target band from bulk of bands and to eliminate the cloning step for sequencing. PMID:25086366
Broom, Donald M
2006-01-01
The term adaptation is used in biology in three different ways. It may refer to changes which occur at the cell and organ level, or at the individual level, or at the level of gene action and evolutionary processes. Adaptation by cells, especially nerve cells helps in: communication within the body, the distinguishing of stimuli, the avoidance of overload and the conservation of energy. The time course and complexity of these mechanisms varies. Adaptive characters of organisms, including adaptive behaviours, increase fitness so this adaptation is evolutionary. The major part of this paper concerns adaptation by individuals and its relationships to welfare. In complex animals, feed forward control is widely used. Individuals predict problems and adapt by acting before the environmental effect is substantial. Much of adaptation involves brain control and animals have a set of needs, located in the brain and acting largely via motivational mechanisms, to regulate life. Needs may be for resources but are also for actions and stimuli which are part of the mechanism which has evolved to obtain the resources. Hence pigs do not just need food but need to be able to carry out actions like rooting in earth or manipulating materials which are part of foraging behaviour. The welfare of an individual is its state as regards its attempts to cope with its environment. This state includes various adaptive mechanisms including feelings and those which cope with disease. The part of welfare which is concerned with coping with pathology is health. Disease, which implies some significant effect of pathology, always results in poor welfare. Welfare varies over a range from very good, when adaptation is effective and there are feelings of pleasure or contentment, to very poor. A key point concerning the concept of individual adaptation in relation to welfare is that welfare may be good or poor while adaptation is occurring. Some adaptation is very easy and energetically cheap and
NASA Astrophysics Data System (ADS)
Cheyney, S.; Fishwick, S.; Hill, I. A.; Linford, N. T.
2015-08-01
Despite the development of advanced processing and interpretation tools for magnetic data sets in the fields of mineral and hydrocarbon industries, these methods have not achieved similar levels of adoption for archaeological or very near surface surveys. Using a synthetic data set we demonstrate that certain methodologies and assumptions used to successfully invert more regional-scale data can lead to large discrepancies between the true and recovered depths when applied to archaeological-type anomalies. We propose variations to the current approach, analysing the choice of the depth-weighting function, mesh design and parameter constraints, to develop an appropriate technique for the 3-D inversion of archaeological-scale data sets. The results show a successful recovery of a synthetic scenario, as well as a case study of a Romano-Celtic temple in the UK. For the case study, the final susceptibility model is compared with two coincident ground penetrating radar surveys, showing a high correlation with the comparative depth slices. The new approach takes interpretation of archaeological data sets beyond a simple 2-D visual interpretation based on pattern recognition.
NASA Astrophysics Data System (ADS)
Bargatze, L. F.
2015-12-01
Active Data Archive Product Tracking (ADAPT) is a collection of software routines that permits one to generate XML metadata files to describe and register data products in support of the NASA Heliophysics Virtual Observatory VxO effort. ADAPT is also a philosophy. The ADAPT concept is to use any and all available metadata associated with scientific data to produce XML metadata descriptions in a consistent, uniform, and organized fashion to provide blanket access to the full complement of data stored on a targeted data server. In this poster, we present an application of ADAPT to describe all of the data products that are stored by using the Common Data File (CDF) format served out by the CDAWEB and SPDF data servers hosted at the NASA Goddard Space Flight Center. These data servers are the primary repositories for NASA Heliophysics data. For this purpose, the ADAPT routines have been used to generate data resource descriptions by using an XML schema named Space Physics Archive, Search, and Extract (SPASE). SPASE is the designated standard for documenting Heliophysics data products, as adopted by the Heliophysics Data and Model Consortium. The set of SPASE XML resource descriptions produced by ADAPT includes high-level descriptions of numerical data products, display data products, or catalogs and also includes low-level "Granule" descriptions. A SPASE Granule is effectively a universal access metadata resource; a Granule associates an individual data file (e.g. a CDF file) with a "parent" high-level data resource description, assigns a resource identifier to the file, and lists the corresponding assess URL(s). The CDAWEB and SPDF file systems were queried to provide the input required by the ADAPT software to create an initial set of SPASE metadata resource descriptions. Then, the CDAWEB and SPDF data repositories were queried subsequently on a nightly basis and the CDF file lists were checked for any changes such as the occurrence of new, modified, or deleted
Sliding mode control of wind-induced vibrations using fuzzy sliding surface and gain adaptation
NASA Astrophysics Data System (ADS)
Thenozhi, Suresh; Yu, Wen
2016-04-01
Although fuzzy/adaptive sliding mode control can reduce the chattering problem in structural vibration control applications, they require the equivalent control and the upper bounds of the system uncertainties. In this paper, we used fuzzy logic to approximate the standard sliding surface and designed a dead-zone adaptive law for tuning the switching gain of the sliding mode control. The stability of the proposed controller is established using Lyapunov stability theory. A six-storey building prototype equipped with an active mass damper has been used to demonstrate the effectiveness of the proposed controller towards the wind-induced vibrations.
Deng, Hua; Li, Han-Xiong; Wu, Yi-Hu
2008-09-01
A new feedback-linearization-based neural network (NN) adaptive control is proposed for unknown nonaffine nonlinear discrete-time systems. An equivalent model in affine-like form is first derived for the original nonaffine discrete-time systems as feedback linearization methods cannot be implemented for such systems. Then, feedback linearization adaptive control is implemented based on the affine-like equivalent model identified with neural networks. Pretraining is not required and the weights of the neural networks used in adaptive control are directly updated online based on the input-output measurement. The dead-zone technique is used to remove the requirement of persistence excitation during the adaptation. With the proposed neural network adaptive control, stability and performance of the closed-loop system are rigorously established. Illustrated examples are provided to validate the theoretical findings. PMID:18779092
NASA Astrophysics Data System (ADS)
Sergienko, Olga
2013-04-01
Since Doug MacAyeal's pioneering studies of the ice-stream basal traction optimizations by control methods, inversions for unknown parameters (e.g., basal traction, accumulation patterns, etc) have become a hallmark of the present-day ice-sheet modeling. The common feature of such inversion exercises is a direct relationship between optimized parameters and observations used in the optimization procedure. For instance, in the standard optimization for basal traction by the control method, ice-stream surface velocities constitute the control data. The optimized basal traction parameters explicitly appear in the momentum equations for the ice-stream velocities (compared to the control data). The inversion for basal traction is carried out by minimization of the cost (or objective, misfit) function that includes the momentum equations facilitated by the Lagrange multipliers. Here, we build upon this idea, and demonstrate how to optimize for parameters indirectly related to observed data using a suite of nested constraints (like Russian dolls) with additional sets of Lagrange multipliers in the cost function. This method opens the opportunity to use data from a variety of sources and types (e.g., velocities, radar layers, surface elevation changes, etc.) in the same optimization process.
Life in the end-Permian dead zone
Looy, Cindy V.; Twitchett, Richard J.; Dilcher, David L.; Van Konijnenburg-Van Cittert, Johanna H. A.; Visscher, Henk
2001-01-01
The fossil record of land plants is an obvious source of information on the dynamics of mass extinctions in the geological past. In conjunction with the end-Permian ecological crisis, ≈250 million years ago, palynological data from East Greenland reveal some unanticipated patterns. We document the significant time lag between terrestrial ecosystem collapse and selective extinction among characteristic Late Permian plants. Furthermore, ecological crisis resulted in an initial increase in plant diversity, instead of a decrease. Paradoxically, these floral patterns correspond to a “dead zone” in the end-Permian faunal record, characterized by a paucity of marine invertebrate megafossils. The time-delayed, end-Permian plant extinctions resemble modeled “extinction debt” responses of multispecies metapopulations to progressive habitat destruction. PMID:11427710
Fuzzy Adaptive Control Design and Discretization for a Class of Nonlinear Uncertain Systems.
Zhao, Xudong; Shi, Peng; Zheng, Xiaolong
2016-06-01
In this paper, tracking control problems are investigated for a class of uncertain nonlinear systems in lower triangular form. First, a state-feedback controller is designed by using adaptive backstepping technique and the universal approximation ability of fuzzy logic systems. During the design procedure, a developed method with less computation is proposed by constructing one maximum adaptive parameter. Furthermore, adaptive controllers with nonsymmetric dead-zone are also designed for the systems. Then, a sampled-data control scheme is presented to discretize the obtained continuous-time controller by using the forward Euler method. It is shown that both proposed continuous and discrete controllers can ensure that the system output tracks the target signal with a small bounded error and the other closed-loop signals remain bounded. Two simulation examples are presented to verify the effectiveness and applicability of the proposed new design techniques. PMID:26208376
Vock, David M; Wolfson, Julian; Bandyopadhyay, Sunayan; Adomavicius, Gediminas; Johnson, Paul E; Vazquez-Benitez, Gabriela; O'Connor, Patrick J
2016-06-01
Models for predicting the probability of experiencing various health outcomes or adverse events over a certain time frame (e.g., having a heart attack in the next 5years) based on individual patient characteristics are important tools for managing patient care. Electronic health data (EHD) are appealing sources of training data because they provide access to large amounts of rich individual-level data from present-day patient populations. However, because EHD are derived by extracting information from administrative and clinical databases, some fraction of subjects will not be under observation for the entire time frame over which one wants to make predictions; this loss to follow-up is often due to disenrollment from the health system. For subjects without complete follow-up, whether or not they experienced the adverse event is unknown, and in statistical terms the event time is said to be right-censored. Most machine learning approaches to the problem have been relatively ad hoc; for example, common approaches for handling observations in which the event status is unknown include (1) discarding those observations, (2) treating them as non-events, (3) splitting those observations into two observations: one where the event occurs and one where the event does not. In this paper, we present a general-purpose approach to account for right-censored outcomes using inverse probability of censoring weighting (IPCW). We illustrate how IPCW can easily be incorporated into a number of existing machine learning algorithms used to mine big health care data including Bayesian networks, k-nearest neighbors, decision trees, and generalized additive models. We then show that our approach leads to better calibrated predictions than the three ad hoc approaches when applied to predicting the 5-year risk of experiencing a cardiovascular adverse event, using EHD from a large U.S. Midwestern healthcare system. PMID:26992568
High resolution 3D nonlinear integrated inversion
NASA Astrophysics Data System (ADS)
Li, Yong; Wang, Xuben; Li, Zhirong; Li, Qiong; Li, Zhengwen
2009-06-01
The high resolution 3D nonlinear integrated inversion method is based on nonlinear theory. Under layer control, the log data from several wells (or all wells) in the study area and seismic trace data adjacent to the wells are input to a network with multiple inputs and outputs and are integratedly trained to obtain an adaptive weight function of the entire study area. Integrated nonlinear mapping relationships are built and updated by the lateral and vertical geologic variations of the reservoirs. Therefore, the inversion process and its inversion results can be constrained and controlled and a stable seismic inversion section with high resolution with velocity inversion, impedance inversion, and density inversion sections, can be gained. Good geologic effects have been obtained in model computation tests and real data processing, which verified that this method has high precision, good practicality, and can be used for quantitative reservoir analysis.
NASA Astrophysics Data System (ADS)
Nath, Saurabh; Mukherjee, Anish; Chatterjee, Souvick; Ganguly, Ranjan; Sen, Swarnendu; Mukhopadhyay, Achintya; Boreyko, Jonathan
2014-11-01
We have observed that capillarity forces may cause floatation in a few non-intuitive configurations. These may be divided into 2 categories: i) floatation of heavier liquid droplets on lighter immiscible ones and ii) fully submerged floatation of lighter liquid droplets in a heavier immiscible medium. We call these counter-intuitive because of the inverse floatation configuration. For case (i) we have identified and studied in detail the several factors affecting the shape and maximum volume of the floating drop. We used water and vegetable oil combinations as test fluids and established the relation between Bond Number and maximum volume contained in a floating drop (in the order of μL). For case (ii), we injected vegetable oil drop-wise into a pool of water. The fully submerged configuration of the drop is not stable and a slight perturbation to the system causes the droplet to burst and float in partially submerged condition. Temporal variation of a characteristic length of the droplet is analyzed using MATLAB image processing. The constraint of small Bond Number establishes the assumption of lubrication regime in the thin gap. A brief theoretical formulation also shows the temporal variation of the gap thickness. Jadavpur University, Jagadis Bose Centre of Excellence, Virginia Tech.
Darwin's "strange inversion of reasoning".
Dennett, Daniel
2009-06-16
Darwin's theory of evolution by natural selection unifies the world of physics with the world of meaning and purpose by proposing a deeply counterintuitive "inversion of reasoning" (according to a 19th century critic): "to make a perfect and beautiful machine, it is not requisite to know how to make it" [MacKenzie RB (1868) (Nisbet & Co., London)]. Turing proposed a similar inversion: to be a perfect and beautiful computing machine, it is not requisite to know what arithmetic is. Together, these ideas help to explain how we human intelligences came to be able to discern the reasons for all of the adaptations of life, including our own. PMID:19528651
Shahnazi, Reza
2015-01-01
An adaptive fuzzy output feedback controller is proposed for a class of uncertain MIMO nonlinear systems with unknown input nonlinearities. The input nonlinearities can be backlash-like hysteresis or dead-zone. Besides, the gains of unknown input nonlinearities are unknown nonlinear functions. Based on universal approximation theorem, the unknown nonlinear functions are approximated by fuzzy systems. The proposed method does not need the availability of the states and an observer based on strictly positive real (SPR) theory is designed to estimate the states. An adaptive robust structure is used to cope with fuzzy approximation error and external disturbances. The semi-global asymptotic stability of the closed-loop system is guaranteed via Lyapunov approach. The applicability of the proposed method is also shown via simulations. PMID:25104646
Precup, Radu-Emil; David, Radu-Codrut; Petriu, Emil M; Radac, Mircea-Bogdan; Preitl, Stefan
2014-11-01
This paper suggests a new generation of optimal PI controllers for a class of servo systems characterized by saturation and dead zone static nonlinearities and second-order models with an integral component. The objective functions are expressed as the integral of time multiplied by absolute error plus the weighted sum of the integrals of output sensitivity functions of the state sensitivity models with respect to two process parametric variations. The PI controller tuning conditions applied to a simplified linear process model involve a single design parameter specific to the extended symmetrical optimum (ESO) method which offers the desired tradeoff to several control system performance indices. An original back-calculation and tracking anti-windup scheme is proposed in order to prevent the integrator wind-up and to compensate for the dead zone nonlinearity of the process. The minimization of the objective functions is carried out in the framework of optimization problems with inequality constraints which guarantee the robust stability with respect to the process parametric variations and the controller robustness. An adaptive gravitational search algorithm (GSA) solves the optimization problems focused on the optimal tuning of the design parameter specific to the ESO method and of the anti-windup tracking gain. A tuning method for PI controllers is proposed as an efficient approach to the design of resilient control systems. The tuning method and the PI controllers are experimentally validated by the adaptive GSA-based tuning of PI controllers for the angular position control of a laboratory servo system. PMID:25330468
Multidimensional NMR inversion without Kronecker products: Multilinear inversion.
Medellín, David; Ravi, Vivek R; Torres-Verdín, Carlos
2016-08-01
Multidimensional NMR inversion using Kronecker products poses several challenges. First, kernel compression is only possible when the kernel matrices are separable, and in recent years, there has been an increasing interest in NMR sequences with non-separable kernels. Second, in three or more dimensions, the singular value decomposition is not unique; therefore kernel compression is not well-defined for higher dimensions. Without kernel compression, the Kronecker product yields matrices that require large amounts of memory, making the inversion intractable for personal computers. Finally, incorporating arbitrary regularization terms is not possible using the Lawson-Hanson (LH) or the Butler-Reeds-Dawson (BRD) algorithms. We develop a minimization-based inversion method that circumvents the above problems by using multilinear forms to perform multidimensional NMR inversion without using kernel compression or Kronecker products. The new method is memory efficient, requiring less than 0.1% of the memory required by the LH or BRD methods. It can also be extended to arbitrary dimensions and adapted to include non-separable kernels, linear constraints, and arbitrary regularization terms. Additionally, it is easy to implement because only a cost function and its first derivative are required to perform the inversion. PMID:27209370
Multidimensional NMR inversion without Kronecker products: Multilinear inversion
NASA Astrophysics Data System (ADS)
Medellín, David; Ravi, Vivek R.; Torres-Verdín, Carlos
2016-08-01
Multidimensional NMR inversion using Kronecker products poses several challenges. First, kernel compression is only possible when the kernel matrices are separable, and in recent years, there has been an increasing interest in NMR sequences with non-separable kernels. Second, in three or more dimensions, the singular value decomposition is not unique; therefore kernel compression is not well-defined for higher dimensions. Without kernel compression, the Kronecker product yields matrices that require large amounts of memory, making the inversion intractable for personal computers. Finally, incorporating arbitrary regularization terms is not possible using the Lawson-Hanson (LH) or the Butler-Reeds-Dawson (BRD) algorithms. We develop a minimization-based inversion method that circumvents the above problems by using multilinear forms to perform multidimensional NMR inversion without using kernel compression or Kronecker products. The new method is memory efficient, requiring less than 0.1% of the memory required by the LH or BRD methods. It can also be extended to arbitrary dimensions and adapted to include non-separable kernels, linear constraints, and arbitrary regularization terms. Additionally, it is easy to implement because only a cost function and its first derivative are required to perform the inversion.
Human inversions and their functional consequences
Puig, Marta; Casillas, Sònia; Villatoro, Sergi
2015-01-01
Polymorphic inversions are a type of structural variants that are difficult to analyze owing to their balanced nature and the location of breakpoints within complex repeated regions. So far, only a handful of inversions have been studied in detail in humans and current knowledge about their possible functional effects is still limited. However, inversions have been related to phenotypic changes and adaptation in multiple species. In this review, we summarize the evidences of the functional impact of inversions in the human genome. First, given that inversions have been shown to inhibit recombination in heterokaryotes, chromosomes displaying different orientation are expected to evolve independently and this may lead to distinct gene-expression patterns. Second, inversions have a role as disease-causing mutations both by directly affecting gene structure or regulation in different ways, and by predisposing to other secondary arrangements in the offspring of inversion carriers. Finally, several inversions show signals of being selected during human evolution. These findings illustrate the potential of inversions to have phenotypic consequences also in humans and emphasize the importance of their inclusion in genome-wide association studies. PMID:25998059
An Inversion Recovery NMR Kinetics Experiment
ERIC Educational Resources Information Center
Williams, Travis J.; Kershaw, Allan D.; Li, Vincent; Wu, Xinping
2011-01-01
A convenient laboratory experiment is described in which NMR magnetization transfer by inversion recovery is used to measure the kinetics and thermochemistry of amide bond rotation. The experiment utilizes Varian spectrometers with the VNMRJ 2.3 software, but can be easily adapted to any NMR platform. The procedures and sample data sets in this…
Bédard, Catherine; Belin, Pascal
2004-07-01
Voice is the carrier of speech but is also an "auditory face" rich in information on the speaker's identity and affective state. Three experiments explored the possibility of a "voice inversion effect," by analogy to the classical "face inversion effect," which could support the hypothesis of a voice-specific module. Experiment 1 consisted of a gender identification task on two syllables pronounced by 90 speakers (boys, girls, men, and women). Experiment 2 consisted of a speaker discrimination task on pairs of syllables (8 men and 8 women). Experiment 3 consisted of an instrument discrimination task on pairs of melodies (8 string and 8 wind instruments). In all three experiments, stimuli were presented in 4 conditions: (1) no inversion; (2) temporal inversion (e.g., backwards speech); (3) frequency inversion centered around 4000 Hz; and (4) around 2500 Hz. Results indicated a significant decrease in performance caused by sound inversion, with a much stronger effect for frequency than for temporal inversion. Interestingly, although frequency inversion markedly affected timbre for both voices and instruments, subjects' performance was still above chance. However, performance at instrument discrimination was much higher than for voices, preventing comparison of inversion effects for voices vs. non-vocal stimuli. Additional experiments will be necessary to conclude on the existence of a possible "voice inversion effect." PMID:15177788
Adaptive control of robotic manipulators
NASA Technical Reports Server (NTRS)
Seraji, H.
1987-01-01
The author presents a novel approach to adaptive control of manipulators to achieve trajectory tracking by the joint angles. The central concept in this approach is the utilization of the manipulator inverse as a feedforward controller. The desired trajectory is applied as an input to the feedforward controller which behaves as the inverse of the manipulator at any operating point; the controller output is used as the driving torque for the manipulator. The controller gains are then updated by an adaptation algorithm derived from MRAC (model reference adaptive control) theory to cope with variations in the manipulator inverse due to changes of the operating point. An adaptive feedback controller and an auxiliary signal are also used to enhance closed-loop stability and to achieve faster adaptation. The proposed control scheme is computationally fast and does not require a priori knowledge of the complex dynamic model or the parameter values of the manipulator or the payload.
ERIC Educational Resources Information Center
Bedard, Catherine; Belin, Pascal
2004-01-01
Voice is the carrier of speech but is also an ''auditory face'' rich in information on the speaker's identity and affective state. Three experiments explored the possibility of a ''voice inversion effect,'' by analogy to the classical ''face inversion effect,'' which could support the hypothesis of a voice-specific module. Experiment 1 consisted…
Jackiewicz, Jason
2009-09-16
With the rapid advances in sophisticated solar modeling and the abundance of high-quality solar pulsation data, efficient and robust inversion techniques are crucial for seismic studies. We present some aspects of an efficient Fourier Optimally Localized Averaging (OLA) inversion method with an example applied to time-distance helioseismology.
Pearson, Bruce R.; Water, Willem van de
2005-03-01
While the ordinary structure function in turbulence is concerned with the statistical moments of the velocity increment {delta}u measured over a distance r, the inverse structure function is related to the distance r where the turbulent velocity exits the interval {delta}u. We study inverse structure functions of wind-tunnel turbulence which covers a range of Reynolds numbers Re{sub {lambda}}=400-1100. We test a recently proposed relation between the scaling exponents of the ordinary structure functions and those of the inverse structure functions [S. Roux and M. H. Jensen, Phys. Rev. E 69, 16309 (2004)]. The relatively large range of Reynolds numbers in our experiment also enables us to address the scaling with Reynolds number that is expected to highlight the intermediate dissipative range. While we firmly establish the (relative) scaling of inverse structure functions, our experimental results fail both predictions. Therefore, the question of the significance of inverse structure functions remains open.
Plasma inverse transition acceleration
Xie, Ming
2001-06-18
It can be proved fundamentally from the reciprocity theorem with which the electromagnetism is endowed that corresponding to each spontaneous process of radiation by a charged particle there is an inverse process which defines a unique acceleration mechanism, from Cherenkov radiation to inverse Cherenkov acceleration (ICA) [1], from Smith-Purcell radiation to inverse Smith-Purcell acceleration (ISPA) [2], and from undulator radiation to inverse undulator acceleration (IUA) [3]. There is no exception. Yet, for nearly 30 years after each of the aforementioned inverse processes has been clarified for laser acceleration, inverse transition acceleration (ITA), despite speculation [4], has remained the least understood, and above all, no practical implementation of ITA has been found, until now. Unlike all its counterparts in which phase synchronism is established one way or the other such that a particle can continuously gain energy from an acceleration wave, the ITA to be discussed here, termed plasma inverse transition acceleration (PITA), operates under fundamentally different principle. As a result, the discovery of PITA has been delayed for decades, waiting for a conceptual breakthrough in accelerator physics: the principle of alternating gradient acceleration [5, 6, 7, 8, 9, 10]. In fact, PITA was invented [7, 8] as one of several realizations of the new principle.
Generalized emissivity inverse problem.
Ming, DengMing; Wen, Tao; Dai, XianXi; Dai, JiXin; Evenson, William E
2002-04-01
Inverse problems have recently drawn considerable attention from the physics community due to of potential widespread applications [K. Chadan and P. C. Sabatier, Inverse Problems in Quantum Scattering Theory, 2nd ed. (Springer Verlag, Berlin, 1989)]. An inverse emissivity problem that determines the emissivity g(nu) from measurements of only the total radiated power J(T) has recently been studied [Tao Wen, DengMing Ming, Xianxi Dai, Jixin Dai, and William E. Evenson, Phys. Rev. E 63, 045601(R) (2001)]. In this paper, a new type of generalized emissivity and transmissivity inverse (GETI) problem is proposed. The present problem differs from our previous work on inverse problems by allowing the unknown (emissivity) function g(nu) to be temperature dependent as well as frequency dependent. Based on published experimental information, we have developed an exact solution formula for this GETI problem. A universal function set suggested for numerical calculation is shown to be robust, making this inversion method practical and convenient for realistic calculations. PMID:12005916
The inverse electroencephalography pipeline
NASA Astrophysics Data System (ADS)
Weinstein, David Michael
The inverse electroencephalography (EEG) problem is defined as determining which regions of the brain are active based on remote measurements recorded with scalp EEG electrodes. An accurate solution to this problem would benefit both fundamental neuroscience research and clinical neuroscience applications. However, constructing accurate patient-specific inverse EEG solutions requires complex modeling, simulation, and visualization algorithms, and to date only a few systems have been developed that provide such capabilities. In this dissertation, a computational system for generating and investigating patient-specific inverse EEG solutions is introduced, and the requirements for each stage of this Inverse EEG Pipeline are defined and discussed. While the requirements of many of the stages are satisfied with existing algorithms, others have motivated research into novel modeling and simulation methods. The principal technical results of this work include novel surface-based volume modeling techniques, an efficient construction for the EEG lead field, and the Open Source release of the Inverse EEG Pipeline software for use by the bioelectric field research community. In this work, the Inverse EEG Pipeline is applied to three research problems in neurology: comparing focal and distributed source imaging algorithms; separating measurements into independent activation components for multifocal epilepsy; and localizing the cortical activity that produces the P300 effect in schizophrenia.
Direct and indirect inversions
NASA Astrophysics Data System (ADS)
Virieux, Jean; Brossier, Romain; Métivier, Ludovic; Operto, Stéphane; Ribodetti, Alessandra
2016-06-01
A bridge is highlighted between the direct inversion and the indirect inversion. They are based on fundamental different approaches: one is looking after a projection from the data space to the model space while the other one is reducing a misfit between observed data and synthetic data obtained from a given model. However, it is possible to obtain similar structures for model perturbation, and we shall focus on P-wave velocity reconstruction. This bridge is built up through the Born approximation linearizing the forward problem with respect to model perturbation and through asymptotic approximations of the Green functions of the wave propagation equation. We first describe the direct inversion and its ingredients and then we focus on a specific misfit function design leading to a indirect inversion. Finally, we shall compare this indirect inversion with more standard least-squares inversion as the FWI, enabling the focus on small weak velocity perturbations on one side and the speed-up of the velocity perturbation reconstruction on the other side. This bridge has been proposed by the group led by Raul Madariaga in the early nineties, emphasizing his leading role in efficient imaging workflows for seismic velocity reconstruction, a drastic requirement at that time.
Aquifer Structure Identification Using Stochastic Inversion
Harp, Dylan R; Dai, Zhenxue; Wolfsberg, Andrew V; Vrugt, Jasper A
2008-01-01
This study presents a stochastic inverse method for aquifer structure identification using sparse geophysical and hydraulic response data. The method is based on updating structure parameters from a transition probability model to iteratively modify the aquifer structure and parameter zonation. The method is extended to the adaptive parameterization of facies hydraulic parameters by including these parameters as optimization variables. The stochastic nature of the statistical structure parameters leads to nonconvex objective functions. A multi-method genetically adaptive evolutionary approach (AMALGAM-SO) was selected to perform the inversion given its search capabilities. Results are obtained as a probabilistic assessment of facies distribution based on indicator cokriging simulation of the optimized structural parameters. The method is illustrated by estimating the structure and facies hydraulic parameters of a synthetic example with a transient hydraulic response.
On the global dynamics of adaptive systems - A study of an elementary example
NASA Technical Reports Server (NTRS)
Espana, Martin D.; Praly, Laurent
1993-01-01
The inherent nonlinear character of adaptive systems poses serious theoretical problems for the analysis of their dynamics. On the other hand, the importance of their dynamic behavior is directly related to the practical interest in predicting such undesirable phenomena as nonlinear oscillations, abrupt transients, intermittence or a high sensitivity with respect to initial conditions. A geometrical/qualitative description of the phase portrait of a discrete-time adaptive system with unmodeled disturbances is given. For this, the motions in the phase space are referred to normally hyperbolic (structurally stable) locally invariant sets. The study is complemented with a local stability analysis of the equilibrium point and periodic solutions. The critical character of adaptive systems under rather usual working conditions is discussed. Special emphasis is put on the causes leading to intermittence. A geometric interpretation of the effects of some commonly used palliatives to this problem is given. The 'dead-zone' approach is studied in more detail. The predicted dynamics are compared with simulation results.
NASA Astrophysics Data System (ADS)
De la sen, M.
2006-03-01
This paper deals with the pole-placement type robust adaptive control of continuous linear systems in the presence of bounded noise and a common class of unmodeled dynamics with the use of multiple estimation schemes working in parallel. The multiestimation scheme consisting of the above set of various single estimation schemes is a tool used to minimize the plant identification error by building an estimate which is a convex combination of the estimates at all time. The weighting functions of the individual estimates are provided at each time by a suboptimization scheme for a quadratic loss function of a possibly filtered tracking error and/or control input. The robust stability of the overall adaptive scheme is ensured by an adaptation relative dead zone which takes into account the contribution of the unmodeled dynamics and bounded noise. The basic results are derived for two different estimation strategies which have either a shared regressor with the plant or individual regressors for the input contribution and its relevant time-derivatives. In this second case, the plant input is obtained through a similar convex combination rule as the one used for the estimators in the first approach. An extension of the basic strategies is also pointed out including a combined use of the suboptimization scheme with a supervisor of past measures for the on-line calculation of the estimator weights in the convex combination.
Electromagnetic inverse scattering
NASA Technical Reports Server (NTRS)
Bojarski, N. N.
1972-01-01
A three-dimensional electromagnetic inverse scattering identity, based on the physical optics approximation, is developed for the monostatic scattered far field cross section of perfect conductors. Uniqueness of this inverse identity is proven. This identity requires complete scattering information for all frequencies and aspect angles. A nonsingular integral equation is developed for the arbitrary case of incomplete frequence and/or aspect angle scattering information. A general closed-form solution to this integral equation is developed, which yields the shape of the scatterer from such incomplete information. A specific practical radar solution is presented. The resolution of this solution is developed, yielding short-pulse target resolution radar system parameter equations. The special cases of two- and one-dimensional inverse scattering and the special case of a priori knowledge of scatterer symmetry are treated in some detail. The merits of this solution over the conventional radar imaging technique are discussed.
Inverse temperature in Superstatistics
NASA Astrophysics Data System (ADS)
Loguercio, Humberto; Davis, Sergio
2016-05-01
In this work, it is shown that there are (at least) three alternative definitions of the inverse temperature for a non-canonical ensemble. These definitions coincide in expectation but, in general, not in their higher moments. We explore in detail the application to the recent formalism of Superstatistics (C. Beck, 2003), and, in particular, to the configurational probability distribution in the microcanonical ensemble.
Application of full-wave inversion to real crosshole data
Song, Z.; Williamson, P.R.
1994-12-31
A 2.5D acoustic frequency domain fullwave inversion method was applied to a real dataset from an open-cast coal exploration site. The only data processing required was the removal of tube waves, because no shear wave arrivals were observed. The inversion is efficient because only a few frequency components are needed. The authors encounter two site-specific problems (source inconsistency and anisotropy) which are addressed by simple adaptations of the inversion algorithm. High resolution results are achieved for both velocity and attenuation reconstructions. The fullwave inversion method combines the advantages of first-arrival travel-time tomography and reflected waves migration. To evaluate the inversion result, they model time domain traces using a source signature estimated by fitting the frequency domain response of the reconstructed model to the observed data across the spectrum. The synthetic traces match the early arrivals in the real data reasonably well.
Gravity inversion in spherical coordinates using tesseroids
NASA Astrophysics Data System (ADS)
Uieda, Leonardo; Barbosa, Valeria C. F.
2014-05-01
Satellite observations of the gravity field have provided geophysicists with exceptionally dense and uniform coverage of data over vast areas. This enables regional or global scale high resolution geophysical investigations. Techniques like forward modeling and inversion of gravity anomalies are routinely used to investigate large geologic structures, such as large igneous provinces, suture zones, intracratonic basins, and the Moho. Accurately modeling such large structures requires taking the sphericity of the Earth into account. A reasonable approximation is to assume a spherical Earth and use spherical coordinates. In recent years, efforts have been made to advance forward modeling in spherical coordinates using tesseroids, particularly with respect to speed and accuracy. Conversely, traditional space domain inverse modeling methods have not yet been adapted to use spherical coordinates and tesseroids. In the literature there are a range of inversion methods that have been developed for Cartesian coordinates and right rectangular prisms. These include methods for estimating the relief of an interface, like the Moho or the basement of a sedimentary basin. Another category includes methods to estimate the density distribution in a medium. The latter apply many algorithms to solve the inverse problem, ranging from analytic solutions to random search methods as well as systematic search methods. We present an adaptation for tesseroids of the systematic search method of "planting anomalous densities". This method can be used to estimate the geometry of geologic structures. As prior information, it requires knowledge of the approximate densities and positions of the structures. The main advantage of this method is its computational efficiency, requiring little computer memory and processing time. We demonstrate the shortcomings and capabilities of this approach using applications to synthetic and field data. Performing the inversion of gravity and gravity gradient
AVO inversion based on inverse operator estimation in trust region
NASA Astrophysics Data System (ADS)
Yin, Xing-Yao; Deng, Wei; Zong, Zhao-Yun
2016-04-01
Amplitude variation with offset (AVO) inversion is widely utilized in exploration geophysics, especially for reservoir prediction and fluid identification. Inverse operator estimation in the trust region algorithm is applied for solving AVO inversion problems in which optimization and inversion directly are integrated. The L1 norm constraint is considered on the basis of reasonable initial model in order to improve effciency and stability during the AVO inversion process. In this study, high-order Zoeppritz approximation is utilized to establish the inversion objective function in which variation of {{v}\\text{p}}/{{v}\\text{s}} with time is taken into consideration. A model test indicates that the algorithm has a relatively higher stability and accuracy than the damp least-squares algorithm. Seismic data inversion is feasible and inversion values of three parameters ({{v}\\text{p}},{{v}\\text{s}},ρ ) maintain good consistency with logging curves.
Molecular inversion probe assay.
Absalan, Farnaz; Ronaghi, Mostafa
2007-01-01
We have described molecular inversion probe technologies for large-scale genetic analyses. This technique provides a comprehensive and powerful tool for the analysis of genetic variation and enables affordable, large-scale studies that will help uncover the genetic basis of complex disease and explain the individual variation in response to therapeutics. Major applications of the molecular inversion probes (MIP) technologies include targeted genotyping from focused regions to whole-genome studies, and allele quantification of genomic rearrangements. The MIP technology (used in the HapMap project) provides an efficient, scalable, and affordable way to score polymorphisms in case/control populations for genetic studies. The MIP technology provides the highest commercially available multiplexing levels and assay conversion rates for targeted genotyping. This enables more informative, genome-wide studies with either the functional (direct detection) approach or the indirect detection approach. PMID:18025701
NASA Astrophysics Data System (ADS)
Kim, Sunghwan; Mitropoulos, Alexander N.; Spitzberg, Joshua D.; Tao, Hu; Kaplan, David L.; Omenetto, Fiorenzo G.
2012-12-01
Periodic nanostructures provide the facility to control and manipulate the flow of light through their lattices. Three-dimensional photonic crystals enable the controlled design of structural colour, which can be varied by infiltrating the structure with different (typically liquid) fillers. Here, we report three-dimensional photonic crystals composed entirely of a purified natural protein (silk fibroin). The biocompatibility of this protein, as well as its favourable material properties and ease of biological functionalization, present opportunities for otherwise unattainable device applications such as bioresorbable integration of structural colour within living tissue or lattice functionalization by means of organic and inorganic material doping. We present a silk inverse opal that demonstrates a pseudo-photonic bandgap in the visible spectrum and show its associated structural colour beneath biological tissue. We also leverage silk's facile dopability to manufacture a gold nanoparticle silk inverse opal and demonstrate patterned heating mediated by enhancement of nanoparticle absorption at the band-edge frequency of the photonic crystal.
Intersections, ideals, and inversion
Vasco, D.W.
1998-10-01
Techniques from computational algebra provide a framework for treating large classes of inverse problems. In particular, the discretization of many types of integral equations and of partial differential equations with undetermined coefficients lead to systems of polynomial equations. The structure of the solution set of such equations may be examined using algebraic techniques.. For example, the existence and dimensionality of the solution set may be determined. Furthermore, it is possible to bound the total number of solutions. The approach is illustrated by a numerical application to the inverse problem associated with the Helmholtz equation. The algebraic methods are used in the inversion of a set of transverse electric (TE) mode magnetotelluric data from Antarctica. The existence of solutions is demonstrated and the number of solutions is found to be finite, bounded from above at 50. The best fitting structure is dominantly onedimensional with a low crustal resistivity of about 2 ohm-m. Such a low value is compatible with studies suggesting lower surface wave velocities than found in typical stable cratons.
Darwin's “strange inversion of reasoning”
Dennett, Daniel
2009-01-01
Darwin's theory of evolution by natural selection unifies the world of physics with the world of meaning and purpose by proposing a deeply counterintuitive “inversion of reasoning” (according to a 19th century critic): “to make a perfect and beautiful machine, it is not requisite to know how to make it” [MacKenzie RB (1868) (Nisbet & Co., London)]. Turing proposed a similar inversion: to be a perfect and beautiful computing machine, it is not requisite to know what arithmetic is. Together, these ideas help to explain how we human intelligences came to be able to discern the reasons for all of the adaptations of life, including our own. PMID:19528651
Inverse avalanches on Abelian sandpiles
Chau, H.F. Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, Illinois 61801-3080 )
1994-11-01
A simple and computationally efficient way of finding inverse avalanches for Abelian sandpiles, called the inverse particle addition operator, is presented. In addition, the method is shown to be optimal in the sense that it requires the minimum amount of computation among methods of the same kind. The method is also conceptually succinct because avalanche and inverse avalanche are placed in the same footing.
Haber, Eldad
2014-03-17
The focus of research was: Developing adaptive mesh for the solution of Maxwell's equations; Developing a parallel framework for time dependent inverse Maxwell's equations; Developing multilevel methods for optimization problems with inequal- ity constraints; A new inversion code for inverse Maxwell's equations in the 0th frequency (DC resistivity); A new inversion code for inverse Maxwell's equations in low frequency regime. Although the research concentrated on electromagnetic forward and in- verse problems the results of the research was applied to the problem of image registration.
Multiscale full waveform inversion
NASA Astrophysics Data System (ADS)
Fichtner, Andreas; Trampert, Jeannot; Cupillard, Paul; Saygin, Erdinc; Taymaz, Tuncay; Capdeville, Yann; Villaseñor, Antonio
2013-07-01
We develop and apply a full waveform inversion method that incorporates seismic data on a wide range of spatio-temporal scales, thereby constraining the details of both crustal and upper-mantle structure. This is intended to further our understanding of crust-mantle interactions that shape the nature of plate tectonics, and to be a step towards improved tomographic models of strongly scale-dependent earth properties, such as attenuation and anisotropy. The inversion for detailed regional earth structure consistently embedded within a large-scale model requires locally refined numerical meshes that allow us to (1) model regional wave propagation at high frequencies, and (2) capture the inferred fine-scale heterogeneities. The smallest local grid spacing sets the upper bound of the largest possible time step used to iteratively advance the seismic wave field. This limitation leads to extreme computational costs in the presence of fine-scale structure, and it inhibits the construction of full waveform tomographic models that describe earth structure on multiple scales. To reduce computational requirements to a feasible level, we design a multigrid approach based on the decomposition of a multiscale earth model with widely varying grid spacings into a family of single-scale models where the grid spacing is approximately uniform. Each of the single-scale models contains a tractable number of grid points, which ensures computational efficiency. The multi-to-single-scale decomposition is the foundation of iterative, gradient-based optimization schemes that simultaneously and consistently invert data on all scales for one multi-scale model. We demonstrate the applicability of our method in a full waveform inversion for Eurasia, with a special focus on Anatolia where coverage is particularly dense. Continental-scale structure is constrained by complete seismic waveforms in the 30-200 s period range. In addition to the well-known structural elements of the Eurasian mantle
Adaptive nonlinear flight control
NASA Astrophysics Data System (ADS)
Rysdyk, Rolf Theoduor
1998-08-01
Research under supervision of Dr. Calise and Dr. Prasad at the Georgia Institute of Technology, School of Aerospace Engineering. has demonstrated the applicability of an adaptive controller architecture. The architecture successfully combines model inversion control with adaptive neural network (NN) compensation to cancel the inversion error. The tiltrotor aircraft provides a specifically interesting control design challenge. The tiltrotor aircraft is capable of converting from stable responsive fixed wing flight to unstable sluggish hover in helicopter configuration. It is desirable to provide the pilot with consistency in handling qualities through a conversion from fixed wing flight to hover. The linear model inversion architecture was adapted by providing frequency separation in the command filter and the error-dynamics, while not exiting the actuator modes. This design of the architecture provides for a model following setup with guaranteed performance. This in turn allowed for convenient implementation of guaranteed handling qualities. A rigorous proof of boundedness is presented making use of compact sets and the LaSalle-Yoshizawa theorem. The analysis allows for the addition of the e-modification which guarantees boundedness of the NN weights in the absence of persistent excitation. The controller is demonstrated on the Generic Tiltrotor Simulator of Bell-Textron and NASA Ames R.C. The model inversion implementation is robustified with respect to unmodeled input dynamics, by adding dynamic nonlinear damping. A proof of boundedness of signals in the system is included. The effectiveness of the robustification is also demonstrated on the XV-15 tiltrotor. The SHL Perceptron NN provides a more powerful application, based on the universal approximation property of this type of NN. The SHL NN based architecture is also robustified with the dynamic nonlinear damping. A proof of boundedness extends the SHL NN augmentation with robustness to unmodeled actuator
Neighborhood inverse consistency preprocessing
Freuder, E.C.; Elfe, C.D.
1996-12-31
Constraint satisfaction consistency preprocessing methods are used to reduce search effort. Time and especially space costs limit the amount of preprocessing that will be cost effective. A new form of consistency preprocessing, neighborhood inverse consistency, can achieve more problem pruning than the usual arc consistency preprocessing in a cost effective manner. There are two basic ideas: (1) Common forms of consistency enforcement basically operate by identifying and remembering solutions to subproblems for which a consistent value cannot be found for some additional problem variable. The space required for this memory can quickly become prohibitive. Inverse consistency basically operates by removing values for variables that are not consistent with any solution to some subproblem involving additional variables. The space requirement is at worst linear. (2) Typically consistency preprocessing achieves some level of consistency uniformly throughout the problem. A subproblem solution will be tested against each additional variable that constrains any subproblem variable. Neighborhood consistency focuses attention on the subproblem formed by the variables that are all constrained by the value in question. By targeting highly relevant subproblems we hope to {open_quotes}skim the cream{close_quotes}, obtaining a high payoff for a limited cost.
When Schools Become Dead Zones of the Imagination: A Critical Pedagogy Manifesto
ERIC Educational Resources Information Center
Giroux, Henry A.
2014-01-01
This article examines the so-called new school reform movement led by a host of right-wing ideologues, billionaires, and foundations. It argues that instead of being reformers, the latter are part of a counter-revolution in American education to dismantle public schools not because they are failing but because they are public and make a claim,…
Open ocean dead-zone in the tropical North Atlantic Ocean
NASA Astrophysics Data System (ADS)
Karstensen, J.; Fiedler, B.; Schütte, F.; Brandt, P.; Körtzinger, A.; Fischer, G.; Zantopp, R.; Hahn, J.; Visbeck, M.; Wallace, D.
2014-12-01
The intermittent appearances of low oxygen environments are a particular thread for marine ecosystems. Here we present first observations of unexpected low (<2 μmol kg-1) oxygen environments in the open waters of the eastern tropical North Atlantic, a region where typically oxygen concentration does not fall below 40 μmol kg-1. The low oxygen zones are created just below the mixed-layer, in the euphotic zone of high productive cyclonic and anticyclonic-modewater eddies. A dynamic boundary is created from the large swirl-velocity against the weak background flow. Hydrographic properties within the eddies are kept constant over periods of several months, while net respiration is elevated by a factor of 3 to 5 reducing the oxygen content. We repeatedly observed low oxygen eddies in the region. The direct impact on the ecosystem is evident from anomalous backscatter behaviour. Satellite derived global eddy statistics do not allow to estimate the large-scale impact of the eddies because their vertical structure (mixed-layer depth, euphotic depth) play a key role in creating the low oxygen environment.
THE FATE OF PLANETESIMALS IN TURBULENT DISKS WITH DEAD ZONES. I. THE TURBULENT STIRRING RECIPE
Okuzumi, Satoshi; Ormel, Chris W.
2013-07-01
Turbulence in protoplanetary disks affects planet formation in many ways. While small dust particles are mainly affected by the aerodynamical coupling with turbulent gas velocity fields, planetesimals and larger bodies are more affected by gravitational interaction with gas density fluctuations. For the latter process, a number of numerical simulations have been performed in recent years, but a fully parameter-independent understanding has not been yet established. In this study, we present simple scaling relations for the planetesimal stirring rate in turbulence driven by magnetorotational instability (MRI), taking into account the stabilization of MRI due to ohmic resistivity. We begin with order-of-magnitude estimates of the turbulence-induced gravitational force acting on solid bodies and associated diffusion coefficients for their orbital elements. We then test the predicted scaling relations using the results of recent ohmic-resistive MHD simulations by Gressel et al. We find that these relations successfully explain the simulation results if we properly fix order-of-unity uncertainties within the estimates. We also update the saturation predictor for the density fluctuation amplitude in MRI-driven turbulence originally proposed by Okuzumi and Hirose. Combination of the scaling relations and saturation predictor allows us to know how the turbulent stirring rate of planetesimals depends on disk parameters such as the gas column density, distance from the central star, vertical resistivity distribution, and net vertical magnetic flux. In Paper II, we apply our recipe to planetesimal accretion to discuss its viability in turbulent disks.
ERIC Educational Resources Information Center
Wyner, Yael
2010-01-01
This inquiry-based activity provides a real-world example that connects to students' everyday seafood choices. In fact, many students went home and insisted to their parents that they should only buy "green" seafood choices. It was also an effective activity because students were able to use what they learned about ocean ecosystems and apply it to…
Generalized multi-point inverse airfoil design
NASA Technical Reports Server (NTRS)
Selig, Michael S.; Maughmer, Mark D.
1991-01-01
In a rather general sense, inverse airfoil design can be taken to mean the problem of specifying a desired set of airfoil characteristics, such as the airfoil maximum thickness ratio, pitching moment, part of the velocity distribution or boundary-layer development, etc., then from this information determine the corresponding airfoil shape. This paper presents a method which approaches the design problem from this perspective. In particular, the airfoil is divided into segments along which, together with the design conditions, either the velocity distribution or boundary-layer development may be prescribed. In addition to these local desired distributions, single parameters like the airfoil thickness can be specified. The problem of finding the airfoil shape is determined by coupling an incompressible, inviscid, inverse airfoil design method with a direct integral boundary-layer analysis method and solving the resulting nonlinear equations via a multidimensional Newton iteration technique. The approach is fast and easily allows for interactive design. It is also flexible and could be adapted to solving compressible, inverse airfoil design problems.
NASA Astrophysics Data System (ADS)
Suortti, Pekka
2016-04-01
A novel concept for a high resolution Compton spectrometer is introduced. 88 keV radiation from an Inverse Compton Compact Source is focused using crossed cylindrically bent Laue-type Si perfect crystals, and dispersed on the sample with a constant energy gradient. Dispersion is compensated exactly at a Ge crystal analyzer, so that the same wavelength shift is observed for all wavelengths of the incident beam. The ThomX source is used as a concrete example. Detailed dimensions and flux estimates at successive locations of the spectrometer are given, and the performance is compared with the dispersion compensating spectrometer at ID15 of the ESRF. The momentum resolution is better than 0.1 atomic units in both cases. The intensity of scattering with the compact source is an order of magnitude smaller, but still adequate for high resolution Compton profile measurements.
Inverse magnetorheological fluids.
Rodríguez-Arco, L; López-López, M T; Zubarev, A Y; Gdula, K; Durán, J D G
2014-09-01
We report a new kind of field-responsive fluid consisting of suspensions of diamagnetic (DM) and ferromagnetic (FM) microparticles in ferrofluids. We designate them as inverse magnetorheological (IMR) fluids for analogy with inverse ferrofluids (IFFs). Observations on the particle self-assembly in IMR fluids upon magnetic field application showed that DM and FM microparticles were assembled into alternating chains oriented along the field direction. We explain such assembly on the basis of the dipolar interaction energy between particles. We also present results on the rheological properties of IMR fluids and, for comparison, those of IFFs and bidispersed magnetorheological (MR) fluids. Interestingly, we found that upon magnetic field application, the rheological properties of IMR fluids were enhanced with respect to bidispersed MR fluids with the same FM particle concentration, by an amount greater than the sum of the isolated contribution of DM particles. Furthermore, the field-induced yield stress was moderately increased when up to 30% of the total FM particle content was replaced with DM particles. Beyond this point, the dependence of the yield stress on the DM content was non-monotonic, as expected for FM concentrations decreasing to zero. We explain these synergistic results by two separate phenomena: the formation of exclusion areas for FM particles due to the perturbation of the magnetic field by DM particles and the dipole-dipole interaction between DM and FM particles, which enhances the field-induced structures. Based on the second phenomenon, we present a theoretical model for the yield stress that semi-quantitatively predicts the experimental results. PMID:25022363
Wavelet Sparse Approximate Inverse Preconditioners
NASA Technical Reports Server (NTRS)
Chan, Tony F.; Tang, W.-P.; Wan, W. L.
1996-01-01
There is an increasing interest in using sparse approximate inverses as preconditioners for Krylov subspace iterative methods. Recent studies of Grote and Huckle and Chow and Saad also show that sparse approximate inverse preconditioner can be effective for a variety of matrices, e.g. Harwell-Boeing collections. Nonetheless a drawback is that it requires rapid decay of the inverse entries so that sparse approximate inverse is possible. However, for the class of matrices that, come from elliptic PDE problems, this assumption may not necessarily hold. Our main idea is to look for a basis, other than the standard one, such that a sparse representation of the inverse is feasible. A crucial observation is that the kind of matrices we are interested in typically have a piecewise smooth inverse. We exploit this fact, by applying wavelet techniques to construct a better sparse approximate inverse in the wavelet basis. We shall justify theoretically and numerically that our approach is effective for matrices with smooth inverse. We emphasize that in this paper we have only presented the idea of wavelet approximate inverses and demonstrated its potential but have not yet developed a highly refined and efficient algorithm.
Inverse problem for Bremsstrahlung radiation
Voss, K.E.; Fisch, N.J.
1991-10-01
For certain predominantly one-dimensional distribution functions, an analytic inversion has been found which yields the velocity distribution of superthermal electrons given their Bremsstrahlung radiation. 5 refs.
Inverse Problems of Thermoelectricity
NASA Astrophysics Data System (ADS)
Anatychuk, L. I.; Luste, O. J.; Kuz, R. V.; Strutinsky, M. N.
2011-05-01
Classical thermoelectricity is based on the use of the Seebeck and Thomson effects that occur in the near-contact areas between n- and p-type materials. A conceptually different approach to thermoelectric power converter design that is based on the law of thermoelectric induction of currents is also known. The efficiency of this approach has already been demonstrated by its first applications. More than 10 basically new types of thermoelements were discovered with properties that cannot be achieved by thermocouple power converters. Therefore, further development of this concept is of practical interest. This paper provides a classification and theory for solving the inverse problems of thermoelectricity that form the basis for devising new thermoelement types. Computer methods for their solution for anisotropic and inhomogeneous media are elaborated. Regularities related to thermoelectric current excitation in anisotropic and inhomogeneous media are established. The possibility of obtaining eddy currents of a particular configuration through control of the temperature field and material parameters for the creation of new thermo- element types is demonstrated for three-dimensional (3D) models of anisotropic and inhomogeneous media.
Inverse problem in hydrogeology
NASA Astrophysics Data System (ADS)
Carrera, Jesús; Alcolea, Andrés; Medina, Agustín; Hidalgo, Juan; Slooten, Luit J.
2005-03-01
The state of the groundwater inverse problem is synthesized. Emphasis is placed on aquifer characterization, where modelers have to deal with conceptual model uncertainty (notably spatial and temporal variability), scale dependence, many types of unknown parameters (transmissivity, recharge, boundary conditions, etc.), nonlinearity, and often low sensitivity of state variables (typically heads and concentrations) to aquifer properties. Because of these difficulties, calibration cannot be separated from the modeling process, as it is sometimes done in other fields. Instead, it should be viewed as one step in the process of understanding aquifer behavior. In fact, it is shown that actual parameter estimation methods do not differ from each other in the essence, though they may differ in the computational details. It is argued that there is ample room for improvement in groundwater inversion: development of user-friendly codes, accommodation of variability through geostatistics, incorporation of geological information and different types of data (temperature, occurrence and concentration of isotopes, age, etc.), proper accounting of uncertainty, etc. Despite this, even with existing codes, automatic calibration facilitates enormously the task of modeling. Therefore, it is contended that its use should become standard practice. L'état du problème inverse des eaux souterraines est synthétisé. L'accent est placé sur la caractérisation de l'aquifère, où les modélisateurs doivent jouer avec l'incertitude des modèles conceptuels (notamment la variabilité spatiale et temporelle), les facteurs d'échelle, plusieurs inconnues sur différents paramètres (transmissivité, recharge, conditions aux limites, etc.), la non linéarité, et souvent la sensibilité de plusieurs variables d'état (charges hydrauliques, concentrations) des propriétés de l'aquifère. A cause de ces difficultés, le calibrage ne peut êtreséparé du processus de modélisation, comme c'est le
Application of Adaptive Autopilot Designs for an Unmanned Aerial Vehicle
NASA Technical Reports Server (NTRS)
Shin, Yoonghyun; Calise, Anthony J.; Motter, Mark A.
2005-01-01
This paper summarizes the application of two adaptive approaches to autopilot design, and presents an evaluation and comparison of the two approaches in simulation for an unmanned aerial vehicle. One approach employs two-stage dynamic inversion and the other employs feedback dynamic inversions based on a command augmentation system. Both are augmented with neural network based adaptive elements. The approaches permit adaptation to both parametric uncertainty and unmodeled dynamics, and incorporate a method that permits adaptation during periods of control saturation. Simulation results for an FQM-117B radio controlled miniature aerial vehicle are presented to illustrate the performance of the neural network based adaptation.
Adaptive management is an approach to natural resource management that emphasizes learning through management where knowledge is incomplete, and when, despite inherent uncertainty, managers and policymakers must act. Unlike a traditional trial and error approach, adaptive managem...
Uncertainty estimation in finite fault inversion
NASA Astrophysics Data System (ADS)
Dettmer, Jan; Cummins, Phil R.; Benavente, Roberto
2016-04-01
This work considers uncertainty estimation for kinematic rupture models in finite fault inversion by Bayesian sampling. Since the general problem of slip estimation on an unknown fault from incomplete and noisy data is highly non-linear and currently intractable, assumptions are typically made to simplify the problem. These almost always include linearization of the time dependence of rupture by considering multiple discrete time windows, and a tessellation of the fault surface into a set of 'subfaults' whose dimensions are fixed below what is subjectively thought to be resolvable by the data. Even non-linear parameterizations are based on a fixed discretization. This results in over-parametrized models which include more parameters than resolvable by the data and require regularization criteria that stabilize the inversion. While it is increasingly common to consider slip uncertainties arising from observational error, the effects of the assumptions implicit in parameterization choices are rarely if ever considered. Here, we show that linearization and discretization assumptions can strongly affect both slip and uncertainty estimates and that therefore the selection of parametrizations should be included in the inference process. We apply Bayesian model selection to study the effect of parametrization choice on inversion results. The Bayesian sampling method which produces inversion results is based on a trans-dimensional rupture discretization which adapts the spatial and temporal parametrization complexity based on data information and does not require regularization. Slip magnitude, direction and rupture velocity are unknowns across the fault and causal first rupture times are obtained by solving the Eikonal equation for a spatially variable rupture-velocity field. The method provides automated local adaptation of rupture complexity based on data information and does not assume globally constant resolution. This is an important quality since seismic data do not
Microwave inverse Cerenkov accelerator
Zhang, T.B.; Marshall, T.C.; LaPointe, M.A.; Hirshfield, J.L.
1997-03-01
A Microwave Inverse Cerenkov Accelerator (MICA) is currently under construction at the Yale Beam Physics Laboratory. The accelerating structure in MICA consists of an axisymmetric dielectrically lined waveguide. For the injection of 6 MeV microbunches from a 2.856 GHz RF gun, and subsequent acceleration by the TM{sub 01} fields, particle simulation studies predict that an acceleration gradient of 6.3 MV/m can be achieved with a traveling-wave power of 15 MW applied to the structure. Synchronous injection into a narrow phase window is shown to allow trapping of all injected particles. The RF fields of the accelerating structure are shown to provide radial focusing, so that longitudinal and transverse emittance growth during acceleration is small, and that no external magnetic fields are required for focusing. For 0.16 nC, 5 psec microbunches, the normalized emittance of the accelerated beam is predicted to be less than 5{pi}mm-mrad. Experiments on sample alumina tubes have been conducted that verify the theoretical dispersion relation for the TM{sub 01} mode over a two-to-one range in frequency. No excitation of axisymmetric or non-axisymmetric competing waveguide modes was observed. High power tests showed that tangential electric fields at the inner surface of an uncoated sample of alumina pipe could be sustained up to at least 8.4 MV/m without breakdown. These considerations suggest that a MICA test accelerator can be built to examine these predictions using an available RF power source, 6 MeV RF gun and associated beam line. {copyright} {ital 1997 American Institute of Physics.}
Integrated inversion using combined wave-equation tomography and full waveform inversion
NASA Astrophysics Data System (ADS)
Wang, Haiyang; Singh, Satish C.; Calandra, Henri
2014-07-01
Wave-equation tomography (WT) and full waveform inversion (FWI) are combined through a hybrid misfit function to estimate high-resolution subsurface structures starting from a poorly constrained initial velocity model. Both methods share the same wavefield forward modelling and inversion schemes, while they differ only on the ways to calculate misfit functions and hence the ways to sample in the model space. Aiming at minimizing the cross-correlation phase delay between synthetic and real data, WT can be used to retrieve the long- and middle-wavelength model components, which are essential to FWI. Compared to ray-based traveltime tomography that is based on asymptotic high-frequency approximation, WT provides a better resolution by exploring the band-limited feature of seismic wavefield. On the other hand, FWI is capable of resolving the short-wavelength model component, complementing the WT. In this study, we apply WT to surface first-arrival refraction data, and apply FWI to both refraction and reflection data. We assign adaptive weights to the two different misfit measurements and build a progressive inversion strategy. To illustrate the advantage of our strategy over conventional `ray tomography + FWI' approach, we show in a synthetic lens test that WT can provide extra subsurface information that is critical for a successful FWI application. To further show the efficiency, we test our strategy on the 2-D Marmousi model where satisfactory inversion results are achieved without much manual intervention. Finally, we apply the inversion strategy to a deep-water seismic data set acquired offshore Sumatra with a 12-km-long streamer. In order to alleviate several practical problems posed by the deep-water setting, we apply downward continuation (DC) to generate a virtual ocean bottom experiment data set prior to inversion. The new geometry after DC boosts up the shallow refractions, as well as avoiding cumbersome modelling through the thick water column, thus
A generalized inversion method: Simultaneous source localization and environmental inversion
NASA Astrophysics Data System (ADS)
Neilsen, Tracianne B.; Knobles, David P.
2002-05-01
The problem of localizing and tracking a source in the shallow ocean is often complicated by uncertainty in the environmental parameters. Likewise, the estimates of environmental parameters in the shallow ocean obtained by inversion methods can be degraded by incorrect information about the source location. To overcome both these common obstacles-environmental mismatch in matched field processing and incorrect source location in geoacoustic inversions-a generalized inversion scheme is developed that includes both source and environmental parameters as unknowns in the inversion. The new technique called systematic decoupling using rotated coordinates (SDRC) expands the original idea of rotated coordinates [M. D. Collins and L. Fishman, J. Acoust. Soc. Am. 98, 1637-1644 (1995)] by using multiple sets of coherent broadband rotated coordinates, each corresponding to a different set of bounds, to systematically decouple the unknowns in a series of simulated annealing inversions. The results of applying the SDRC inversion method to data from the Area Characterization Test II experiment performed on the New Jersey continental shelf are presented. [Work supported by ONR.
Temperature Inversions Have Cold Bottoms.
ERIC Educational Resources Information Center
Bohren, Craig F.; Brown, Gail M.
1982-01-01
Uses discussion and illustrations of several demonstrations on air temperature differences and atmospheric stability to explain the phenomena of temperature inversions. Relates this to the smog in Los Angeles and discusses the implications. (DC)
NASA Astrophysics Data System (ADS)
Mahan, G. D.
2014-09-01
We calculate the binding energy of an electron bound to a donor in a semiconductor inverse opal. Inverse opals have two kinds of cavities, which we call octahedral and tetrahedral, according to their group symmetry. We put the donor in the center of each of these two cavities and obtain the binding energy. The binding energies become very large when the inverse opal is made from templates with small spheres. For spheres less than 50 nm in diameter, the donor binding can increase to several times its unconfined value. Then electrons become tightly bound to the donor and are unlikely to be thermally activated to the semiconductor conduction band. This conclusion suggests that inverse opals will be poor conductors.
Mahan, G. D.
2014-09-21
We calculate the binding energy of an electron bound to a donor in a semiconductor inverse opal. Inverse opals have two kinds of cavities, which we call octahedral and tetrahedral, according to their group symmetry. We put the donor in the center of each of these two cavities and obtain the binding energy. The binding energies become very large when the inverse opal is made from templates with small spheres. For spheres less than 50 nm in diameter, the donor binding can increase to several times its unconfined value. Then electrons become tightly bound to the donor and are unlikely to be thermally activated to the semiconductor conduction band. This conclusion suggests that inverse opals will be poor conductors.
Inversion layer MOS solar cells
NASA Technical Reports Server (NTRS)
Ho, Fat Duen
1986-01-01
Inversion layer (IL) Metal Oxide Semiconductor (MOS) solar cells were fabricated. The fabrication technique and problems are discussed. A plan for modeling IL cells is presented. Future work in this area is addressed.
Uterine Inversion; A case report.
Bouchikhi, C; Saadi, H; Fakhir, B; Chaara, H; Bouguern, H; Banani, A; Melhouf, Ma
2008-01-01
The puerperal uterine inversion is a rare and severe complication occurring in the third stage of labour. The mechanisms are not completely known. However, extrinsic factors such as oxytocic arrests after a prolonged labour, umbilical cord traction or abdominal expression are pointed. Other intrinsic factors such as primiparity, uterine hypotonia, various placental localizations, fundic myoma or short umbilical cord were also reported. The diagnosis of the uterine inversion is mainly supported by clinical symptoms. It is based on three elements: haemorrhage, shock and a strong pelvic pain. The immediate treatment of the uterine inversion is required. It is based on a medical reanimation associated with firstly a manual reduction then surgical treatment using various techniques. We report an observation of a 25 years old grand multiparous patient with a subacute uterine inversion after delivery at home. PMID:21516244
Uterine Inversion; A case report
Bouchikhi, C; Saadi, H; Fakhir, B; Chaara, H; Bouguern, H; Banani, A; Melhouf, MA
2008-01-01
The puerperal uterine inversion is a rare and severe complication occurring in the third stage of labour. The mechanisms are not completely known. However, extrinsic factors such as oxytocic arrests after a prolonged labour, umbilical cord traction or abdominal expression are pointed. Other intrinsic factors such as primiparity, uterine hypotonia, various placental localizations, fundic myoma or short umbilical cord were also reported. The diagnosis of the uterine inversion is mainly supported by clinical symptoms. It is based on three elements: haemorrhage, shock and a strong pelvic pain. The immediate treatment of the uterine inversion is required. It is based on a medical reanimation associated with firstly a manual reduction then surgical treatment using various techniques. We report an observation of a 25 years old grand multiparous patient with a subacute uterine inversion after delivery at home. PMID:21516244
An inverse acoustic waveguide problem in the time domain
NASA Astrophysics Data System (ADS)
Monk, Peter; Selgas, Virginia
2016-05-01
We consider the problem of locating an obstacle in a waveguide from time domain measurements of causal waves. More precisely, we assume that we are given the scattered field due to point sources placed on a surface located inside the waveguide away from the obstacle, where the scattered field is measured on the same surface. From this multi-static scattering data we wish to determine the position and shape of an obstacle in the waveguide. To deal with this inverse problem, we adapt and analyze the time domain linear sampling method. This involves proving new time domain estimates for the forward problem, as well as analyzing several time domain operators arising in the inversion scheme. We also implement the inversion algorithm and provide numerical results in two-dimensions using synthetic data.
Computation of inverse magnetic cascades
NASA Technical Reports Server (NTRS)
Montgomery, D.
1981-01-01
Inverse cascades of magnetic quantities for turbulent incompressible magnetohydrodynamics are reviewed, for two and three dimensions. The theory is extended to the Strauss equations, a description intermediate between two and three dimensions appropriate to Tokamak magnetofluids. Consideration of the absolute equilibrium Gibbs ensemble for the system leads to a prediction of an inverse cascade of magnetic helicity, which may manifest itself as a major disruption. An agenda for computational investigation of this conjecture is proposed.
Computation of inverse magnetic cascades
Montgomery, D.
1981-10-01
Inverse cascades of magnetic quantities for turbulent incompressible magnetohydrodynamics are reviewed, for two and three dimensions. The theory is extended to the Strauss equations, a description intermediate between two and three dimensions appropriate to tokamak magnetofluids. Consideration of the absolute equilibrium Gibbs ensemble for the system leads to a prediction of an inverse cascade of magnetic helicity, which may manifest itself as a major disruption. An agenda for computational investigation of this conjecture is proposed.
Inverse Ising inference with correlated samples
NASA Astrophysics Data System (ADS)
Obermayer, Benedikt; Levine, Erel
2014-12-01
Correlations between two variables of a high-dimensional system can be indicative of an underlying interaction, but can also result from indirect effects. Inverse Ising inference is a method to distinguish one from the other. Essentially, the parameters of the least constrained statistical model are learned from the observed correlations such that direct interactions can be separated from indirect correlations. Among many other applications, this approach has been helpful for protein structure prediction, because residues which interact in the 3D structure often show correlated substitutions in a multiple sequence alignment. In this context, samples used for inference are not independent but share an evolutionary history on a phylogenetic tree. Here, we discuss the effects of correlations between samples on global inference. Such correlations could arise due to phylogeny but also via other slow dynamical processes. We present a simple analytical model to address the resulting inference biases, and develop an exact method accounting for background correlations in alignment data by combining phylogenetic modeling with an adaptive cluster expansion algorithm. We find that popular reweighting schemes are only marginally effective at removing phylogenetic bias, suggest a rescaling strategy that yields better results, and provide evidence that our conclusions carry over to the frequently used mean-field approach to the inverse Ising problem.
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.
Inverse solutions for electric and potential field imaging
NASA Astrophysics Data System (ADS)
Johnson, Christopher R.; MacLeod, Robert S.
1993-08-01
One of the fundamental problems in theoretical electrocardiography can be characterized by an inverse problem. In this paper, we present new methods for achieving better estimates of heart surface potential distributions in terms of torso potentials through an inverse procedure. First, an adaptive meshing algorithm is described which minimizes the error in the forward problem due to spatial discretization. We have found that since the inverse problem relies directly on the accuracy of the forward solution, adaptive meshing produces a more accurate inverse transfer matrix. Secondly, we introduce a new local regularization procedure. This method works by breaking the global transfer matrix into sub-matrices and performing regularization only on those sub-matrices which have large condition numbers. Furthermore, the regularization parameters are specifically 'tuned' for each sub-matrix using an a priori scheme based on the L-curve method. This local regularization method provides substantial increases in accuracy when compared to global regularization schemes. Finally, we present specific examples of the implementation of these schemes using models derived from magnetic resonance imaging data from a human subject.
Barrett, Harrison H.; Furenlid, Lars R.; Freed, Melanie; Hesterman, Jacob Y.; Kupinski, Matthew A.; Clarkson, Eric; Whitaker, Meredith K.
2008-01-01
Adaptive imaging systems alter their data-acquisition configuration or protocol in response to the image information received. An adaptive pinhole single-photon emission computed tomography (SPECT) system might acquire an initial scout image to obtain preliminary information about the radiotracer distribution and then adjust the configuration or sizes of the pinholes, the magnifications, or the projection angles in order to improve performance. This paper briefly describes two small-animal SPECT systems that allow this flexibility and then presents a framework for evaluating adaptive systems in general, and adaptive SPECT systems in particular. The evaluation is in terms of the performance of linear observers on detection or estimation tasks. Expressions are derived for the ideal linear (Hotelling) observer and the ideal linear (Wiener) estimator with adaptive imaging. Detailed expressions for the performance figures of merit are given, and possible adaptation rules are discussed. PMID:18541485
Megabase-Scale Inversion Polymorphism in the Wild Ancestor of Maize
Fang, Zhou; Pyhäjärvi, Tanja; Weber, Allison L.; Dawe, R. Kelly; Glaubitz, Jeffrey C.; González, José de Jesus Sánchez; Ross-Ibarra, Claudia; Doebley, John; Morrell, Peter L.; Ross-Ibarra, Jeffrey
2012-01-01
Chromosomal inversions are thought to play a special role in local adaptation, through dramatic suppression of recombination, which favors the maintenance of locally adapted alleles. However, relatively few inversions have been characterized in population genomic data. On the basis of single-nucleotide polymorphism (SNP) genotyping across a large panel of Zea mays, we have identified an ∼50-Mb region on the short arm of chromosome 1 where patterns of polymorphism are highly consistent with a polymorphic paracentric inversion that captures >700 genes. Comparison to other taxa in Zea and Tripsacum suggests that the derived, inverted state is present only in the wild Z. mays subspecies parviglumis and mexicana and is completely absent in domesticated maize. Patterns of polymorphism suggest that the inversion is ancient and geographically widespread in parviglumis. Cytological screens find little evidence for inversion loops, suggesting that inversion heterozygotes may suffer few crossover-induced fitness consequences. The inversion polymorphism shows evidence of adaptive evolution, including a strong altitudinal cline, a statistical association with environmental variables and phenotypic traits, and a skewed haplotype frequency spectrum for inverted alleles. PMID:22542971
Joint inversion of aquifer test, MRS, and TEM data
NASA Astrophysics Data System (ADS)
Vilhelmsen, Troels N.; Behroozmand, Ahmad A.; Christensen, Steen; Nielsen, Toke H.
2014-05-01
This paper presents two methods for joint inversion of aquifer test data, magnetic resonance sounding (MRS) data, and transient electromagnetic data acquired from a multilayer hydrogeological system. The link between the MRS model and the groundwater model is created by tying hydraulic conductivities (k) derived from MRS parameters to those of the groundwater model. Method 1 applies k estimated from MRS directly in the groundwater model, during the inversion. Method 2 on the other hand uses the petrophysical relation as a regularization constraint that only enforces k estimated for the groundwater model to be equal to MRS derived k to the extent that data can be fitted. Both methodologies can jointly calibrate parameters pertaining to the individual models as well as a parameter pertaining to the petrophysical relation. This allows the petrophysical relation to adapt to the local conditions during the inversion. The methods are tested using a synthetic data set as well as a field data set. In combination, the two case studies show that the joint methods can constrain the inversion to achieve estimates of k, decay times, and water contents for a leaky confined aquifer system. We show that the geophysical data can assist in determining otherwise insensitive k, and vice versa. Based on our experiments and results, we mainly advocate the future application of method 2 since this seems to produce the most reliable results, has a faster inversion runtime, and is applicable also for linking k of 3-D groundwater flow models to multiple MRS soundings.
Global inversion for anisotropy during full-waveform inversion
NASA Astrophysics Data System (ADS)
Debens, H. A.; Warner, M.; Umpleby, A.
2015-12-01
Full-waveform inversion (FWI) is a powerful tool for quantitative estimation of high-resolution high-fidelity models of subsurface seismic parameters, typically P-wave velocity. The solution to FWI's posed nonlinear inverse problem is obtained via an iterative series of linearized local updates to a start model, assuming this model lies within the basin of attraction to the global minimum. Thanks to many successful published applications to three-dimensional (3D) field datasets, its advance has been rapid and driven in large-part by the oil and gas industry. The consideration of seismic anisotropy during FWI is of vital importance, as it holds influence over both the kinematics and dynamics of seismic waveforms. If not appropriately taken into account then inadequacies in the anisotropy model are likely to manifest as significant error in the recovered velocity model. Conventionally, anisotropic FWI employs either an a priori anisotropy model, held fixed during FWI, or it uses a multi-parameter local inversion scheme to recover the anisotropy as part of the FWI; both of these methods can be problematic. Constructing an anisotropy model prior to FWI often involves intensive (and hence expensive) iterative procedures, such as travel-time tomography or moveout velocity analysis. On the other hand, introducing multiple parameters to FWI itself increases the complexity of what is already an underdetermined inverse problem. We propose that global rather than local FWI can be used to recover the long-wavelength acoustic anisotropy model, and that this can then be followed by more-conventional local FWI to recover the detailed model. We validate this approach using a full 3D field dataset, demonstrating that it avoids problems associated to crosstalk that can bedevil local inversion schemes, and reconciles well with in situ borehole measurements. Although our approach includes a global inversion for anisotropy, it is nonetheless affordable and practical for 3D field data.
Geophysical Inversion Through Hierarchical Scheme
NASA Astrophysics Data System (ADS)
Furman, A.; Huisman, J. A.
2010-12-01
Geophysical investigation is a powerful tool that allows non-invasive and non-destructive mapping of subsurface states and properties. However, non-uniqueness associated with the inversion process prevents the quantitative use of these methods. One major direction researchers are going is constraining the inverse problem by hydrological observations and models. An alternative to the commonly used direct inversion methods are global optimization schemes (such as genetic algorithms and Monte Carlo Markov Chain methods). However, the major limitation here is the desired high resolution of the tomographic image, which leads to a large number of parameters and an unreasonably high computational effort when using global optimization schemes. Two innovative schemes are presented here. First, a hierarchical approach is used to reduce the computational effort for the global optimization. Solution is achieved for coarse spatial resolution, and this solution is used as the starting point for finer scheme. We show that the computational effort is reduced in this way dramatically. Second, we use a direct ERT inversion as the starting point for global optimization. In this case preliminary results show that the outcome is not necessarily beneficial, probably because of spatial mismatch between the results of the direct inversion and the true resistivity field.
Optimization and geophysical inverse problems
Barhen, J.; Berryman, J.G.; Borcea, L.; Dennis, J.; de Groot-Hedlin, C.; Gilbert, F.; Gill, P.; Heinkenschloss, M.; Johnson, L.; McEvilly, T.; More, J.; Newman, G.; Oldenburg, D.; Parker, P.; Porto, B.; Sen, M.; Torczon, V.; Vasco, D.; Woodward, N.B.
2000-10-01
A fundamental part of geophysics is to make inferences about the interior of the earth on the basis of data collected at or near the surface of the earth. In almost all cases these measured data are only indirectly related to the properties of the earth that are of interest, so an inverse problem must be solved in order to obtain estimates of the physical properties within the earth. In February of 1999 the U.S. Department of Energy sponsored a workshop that was intended to examine the methods currently being used to solve geophysical inverse problems and to consider what new approaches should be explored in the future. The interdisciplinary area between inverse problems in geophysics and optimization methods in mathematics was specifically targeted as one where an interchange of ideas was likely to be fruitful. Thus about half of the participants were actively involved in solving geophysical inverse problems and about half were actively involved in research on general optimization methods. This report presents some of the topics that were explored at the workshop and the conclusions that were reached. In general, the objective of a geophysical inverse problem is to find an earth model, described by a set of physical parameters, that is consistent with the observational data. It is usually assumed that the forward problem, that of calculating simulated data for an earth model, is well enough understood so that reasonably accurate synthetic data can be generated for an arbitrary model. The inverse problem is then posed as an optimization problem, where the function to be optimized is variously called the objective function, misfit function, or fitness function. The objective function is typically some measure of the difference between observational data and synthetic data calculated for a trial model. However, because of incomplete and inaccurate data, the objective function often incorporates some additional form of regularization, such as a measure of smoothness
Complex Patterns of Local Adaptation in Teosinte
Pyhäjärvi, Tanja; Hufford, Matthew B.; Mezmouk, Sofiane; Ross-Ibarra, Jeffrey
2013-01-01
Populations of widely distributed species encounter and must adapt to local environmental conditions. However, comprehensive characterization of the genetic basis of adaptation is demanding, requiring genome-wide genotype data, multiple sampled populations, and an understanding of population structure and potential selection pressures. Here, we used single-nucleotide polymorphism genotyping and data on numerous environmental variables to describe the genetic basis of local adaptation in 21 populations of teosinte, the wild ancestor of maize. We found complex hierarchical genetic structure created by altitude, dispersal events, and admixture among subspecies, which complicated identification of locally beneficial alleles. Patterns of linkage disequilibrium revealed four large putative inversion polymorphisms showing clinal patterns of frequency. Population differentiation and environmental correlations suggest that both inversions and intergenic polymorphisms are involved in local adaptation. PMID:23902747
Thermoelectric properties of inverse opals
NASA Astrophysics Data System (ADS)
Mahan, G. D.; Poilvert, N.; Crespi, V. H.
2016-02-01
Rayleigh's method [Philos. Mag. Ser. 5 34, 481 (1892)] is used to solve for the classical thermoelectric equations in inverse opals. His theory predicts that in an inverse opal, with periodic holes, the Seebeck coefficient and the figure of merit are identical to that of the bulk material. We also provide a major revision to Rayleigh's method, in using the electrochemical potential as an important variable, instead of the electrostatic potential. We also show that in some cases, the thermal boundary resistance is important in the effective thermal conductivity.
Multiphase inverse modeling: An Overview
Finsterle, S.
1998-03-01
Inverse modeling is a technique to derive model-related parameters from a variety of observations made on hydrogeologic systems, from small-scale laboratory experiments to field tests to long-term geothermal reservoir responses. If properly chosen, these observations contain information about the system behavior that is relevant to the performance of a geothermal field. Estimating model-related parameters and reducing their uncertainty is an important step in model development, because errors in the parameters constitute a major source of prediction errors. This paper contains an overview of inverse modeling applications using the ITOUGH2 code, demonstrating the possibilities and limitations of a formalized approach to the parameter estimation problem.
NASA Astrophysics Data System (ADS)
Ghasemi-Nejhad, Mehrdad N.
2013-04-01
designed to have tip-tilt pointing and simultaneous multi-degree-of-freedom vibration isolation capability for pointing stabilization. Analytical approaches have been employed for determining the loads in the components as well as optimizing the design of the system. The different critical components such as telescope tube struts, flexure joints, and the secondary mirror mount have been designed and analyzed using finite element technique. The Simultaneous Precision Positioning and Vibration Suppression (SPPVS) smart composites platforms for the adaptive TVC and adaptive composite telescope are analogous (e.g., see work by Ghasemi-Nejhad and co-workers [1, 2]), where innovative concepts and control strategies are introduced, and experimental verifications of simultaneous thrust vector control and vibration isolation of satellites were performed. The smart composite platforms function as an active structural interface between the main thruster of a satellite and the satellite structure for the adaptive TVC application and as an active structural interface between the main smart composite telescope and the satellite structure for the adaptive laser communication application. The cascaded multiple feedback loops compensate the hysteresis (for piezoelectric stacks inside the three linear actuators that individually have simultaneous precision positioning and vibration suppression), dead-zone, back-lash, and friction nonlinearities very well, and provide precision and quick smart platform control and satisfactory thrust vector control capability. In addition, for example for the adaptive TVC, the experimental results show that the smart composite platform satisfactorily provided precision and fast smart platform control as well as the satisfactory thrust vector control capability. The vibration controller isolated 97% of the vibration energy due to the thruster firing.
Inversions and the origin of behavioral differences in cod.
Samuk, Kieran
2016-05-01
How does adaptation manage to occur in the face of overwhelming gene flow? One popular idea is that the suppression of recombination, for example the fixation of a chromosomal inversion, can maintain linkage disequilibrium between groups of locally adapted alleles that would otherwise be degraded by gene flow. This idea has captured the imagination of many geneticists and evolutionary biologists, but we still have only a basic understanding of its general importance. In this issue of Molecular Ecology, Kirubakaran et al. () examine the role of recombination suppression in a particularly fascinating example of adaptation in the face of gene flow: the evolution of migratory differences between interbreeding populations of cod. Along the north coast of Norway, two types of cod breed in the near-shore waters: a 'stationary' form that lives near the coast year round, and a 'migratory' form that lives far offshore and only returns to the coast to breed. Using a combination of approaches, Kirubakaran et al. () deftly demonstrate that the migratory form has completely fixed two adjacent inversions containing a suite of genes closely connected to migratory behaviour and feeding differences. This work provides an excellent example of how recombination suppression can facilitate adaptive divergence, and helps us understand the geographic and temporal scales over which genomic structural variation evolves. PMID:27213696
Inverse free electron laser accelerator
Fisher, A.; Gallardo, J.; Sandweiss, J.; van Steenbergen, A. )
1992-07-01
The study of the INVERSE FREE ELECTRON LASER, as a potential mode of electron acceleration, is being pursued at Brookhaven National Laboratory. Recent studies have focussed on the development of a low energy, high gradient, multi stage linear accelerator. The elementary ingredients for the IFEL interaction are the 50 MeV Linac e[sup [minus
Inverse Free Electron Laser accelerator
Fisher, A.; Gallardo, J.; van Steenbergen, A. ); Sandweiss, J. )
1992-09-01
The study of the INVERSE FREE ELECTRON LASER, as a potential mode of electron acceleration, is being pursued at Brookhaven National Laboratory. Recent studies have focussed on the development of a low energy, high gradient, multi stage linear accelerator. The elementary ingredients for the IFEL interaction are the 50 MeV Linac e[sup [minus
Action Understanding as Inverse Planning
ERIC Educational Resources Information Center
Baker, Chris L.; Saxe, Rebecca; Tenenbaum, Joshua B.
2009-01-01
Humans are adept at inferring the mental states underlying other agents' actions, such as goals, beliefs, desires, emotions and other thoughts. We propose a computational framework based on Bayesian inverse planning for modeling human action understanding. The framework represents an intuitive theory of intentional agents' behavior based on the…
High-Resolution Inverse-Based Determination of Seismic-Velocity Structure in Basins
NASA Astrophysics Data System (ADS)
Akcelik, V.; Bielak, J.; Epanomeritakis, I.; Ghattas, O.
2004-12-01
Starting with the pioneering work of Aki, Christoffersson, and Husebye in 1976, there has been an increasing interest in developing inversion techniques for determining the three-dimensional crustal velocity structure in seismic regions. In this paper we describe a methodology that capitalizes on recent advances in optimization methods to adapt, extend, and refine powerful nonlinear Newton-Krylov adjoint-based inverse wave propagation algorithms to two- and three-dimensional velocity structure and kinematic source inversion problems. We present results of high resolution models for two-dimensional sedimentary valleys undergoing antiplane motion, and three dimensional acoustic approximations of models of the San Fernando Valley using parallel scalable inversion algorithms that overcome many of the difficulties particular to inverse heterogeneous wave propagation problems.
Non-linearity in Bayesian 1-D magnetotelluric inversion
NASA Astrophysics Data System (ADS)
Guo, Rongwen; Dosso, Stan E.; Liu, Jianxin; Dettmer, Jan; Tong, Xiaozhong
2011-05-01
This paper applies a Bayesian approach to examine non-linearity for the 1-D magnetotelluric (MT) inverse problem. In a Bayesian formulation the posterior probability density (PPD), which combines data and prior information, is interpreted in terms of parameter estimates and uncertainties, which requires optimizing and integrating the PPD. Much work on 1-D MT inversion has been based on (approximate) linearized solutions, but more recently fully non-linear (numerical) approaches have been applied. This paper directly compares results of linearized and non-linear uncertainty estimation for 1-D MT inversion; to do so, advanced methods for both approaches are applied. In the non-linear formulation used here, numerical optimization is carried out using an adaptive-hybrid algorithm. Numerical integration applies Metropolis-Hastings sampling, rotated to a principal-component parameter space for efficient sampling of correlated parameters, and employing non-unity sampling temperatures to ensure global sampling. Since appropriate model parametrizations are generally not known a priori, both under- and overparametrized approaches are considered. For underparametrization, the Bayesian information criterion is applied to determine the number of layers consistent with the resolving power of the data. For overparametrization, prior information is included which favours simple structure in a manner similar to regularized inversion. The data variance and/or trade-off parameter regulating data and prior information are treated in several ways, including applying fixed optimal estimates (an empirical Bayesian approach) or including them as hyperparameters in the sampling (hierarchical Bayesian). The latter approach has the benefit of accounting for the uncertainty in the hyperparameters in estimating model parameter uncertainties. Non-linear and linearized inversion results are compared for synthetic test cases and for the measured COPROD1 MT data by considering marginal probability
Adaptive compressive sensing camera
NASA Astrophysics Data System (ADS)
Hsu, Charles; Hsu, Ming K.; Cha, Jae; Iwamura, Tomo; Landa, Joseph; Nguyen, Charles; Szu, Harold
2013-05-01
We have embedded Adaptive Compressive Sensing (ACS) algorithm on Charge-Coupled-Device (CCD) camera based on the simplest concept that each pixel is a charge bucket, and the charges comes from Einstein photoelectric conversion effect. Applying the manufactory design principle, we only allow altering each working component at a minimum one step. We then simulated what would be such a camera can do for real world persistent surveillance taking into account of diurnal, all weather, and seasonal variations. The data storage has saved immensely, and the order of magnitude of saving is inversely proportional to target angular speed. We did design two new components of CCD camera. Due to the matured CMOS (Complementary metal-oxide-semiconductor) technology, the on-chip Sample and Hold (SAH) circuitry can be designed for a dual Photon Detector (PD) analog circuitry for changedetection that predicts skipping or going forward at a sufficient sampling frame rate. For an admitted frame, there is a purely random sparse matrix [Φ] which is implemented at each bucket pixel level the charge transport bias voltage toward its neighborhood buckets or not, and if not, it goes to the ground drainage. Since the snapshot image is not a video, we could not apply the usual MPEG video compression and Hoffman entropy codec as well as powerful WaveNet Wrapper on sensor level. We shall compare (i) Pre-Processing FFT and a threshold of significant Fourier mode components and inverse FFT to check PSNR; (ii) Post-Processing image recovery will be selectively done by CDT&D adaptive version of linear programming at L1 minimization and L2 similarity. For (ii) we need to determine in new frames selection by SAH circuitry (i) the degree of information (d.o.i) K(t) dictates the purely random linear sparse combination of measurement data a la [Φ]M,N M(t) = K(t) Log N(t).
Adaptive support vector regression for UAV flight control.
Shin, Jongho; Jin Kim, H; Kim, Youdan
2011-01-01
This paper explores an application of support vector regression for adaptive control of an unmanned aerial vehicle (UAV). Unlike neural networks, support vector regression (SVR) generates global solutions, because SVR basically solves quadratic programming (QP) problems. With this advantage, the input-output feedback-linearized inverse dynamic model and the compensation term for the inversion error are identified off-line, which we call I-SVR (inversion SVR) and C-SVR (compensation SVR), respectively. In order to compensate for the inversion error and the unexpected uncertainty, an online adaptation algorithm for the C-SVR is proposed. Then, the stability of the overall error dynamics is analyzed by the uniformly ultimately bounded property in the nonlinear system theory. In order to validate the effectiveness of the proposed adaptive controller, numerical simulations are performed on the UAV model. PMID:20970303
Simultaneous inversion for velocity and attenuation by waveform tomography
NASA Astrophysics Data System (ADS)
Gao, Fengxia; Wang, Yanghua
2016-08-01
Seismic waveform tomography can invert for the velocity and attenuation (Q- 1) variations simultaneously. For this simultaneous inversion, we propose two strategies for waveform tomography. First, we analyze the contributions of the real part and the imaginary part of the gradients, associated with the velocity and attenuation parameters respectively, and determine that the combination of the real part of the gradient subvector for the velocity parameter and the imaginary part of the gradient subvector for the attenuation parameter would produce an optimal inversion result. Second, we attempt to balance the sensitivities of the objective function to the velocity and the attenuation parameters. Considering the magnitude differences between these two-type parameters in the simultaneous inversion, we apply preliminarily a normalization to both the velocity model and the attenuation model. However, for balancing their sensitivities, we further adjust the corresponding model updates using a tuning factor. We determine this tuning parameter adaptively, based on the sensitivities of these two parameters, at each iteration. Numerical tests demonstrate the feasibility and reliability of these two strategies in full waveform inversion.
Nonlinear functional approximation with networks using adaptive neurons
NASA Technical Reports Server (NTRS)
Tawel, Raoul
1992-01-01
A novel mathematical framework for the rapid learning of nonlinear mappings and topological transformations is presented. It is based on allowing the neuron's parameters to adapt as a function of learning. This fully recurrent adaptive neuron model (ANM) has been successfully applied to complex nonlinear function approximation problems such as the highly degenerate inverse kinematics problem in robotics.
Puerma, Eva; Orengo, Dorcas J.; Aguadé, Montserrat
2016-01-01
Chromosomal inversions can contribute to the adaptation of organisms to their environment by capturing particular advantageous allelic combinations of a set of genes included in the inverted fragment and also by advantageous functional changes due to the inversion process itself that might affect not only the expression of flanking genes but also their dose and structure. Of the two mechanisms originating inversions —ectopic recombination, and staggered double-strand breaks and subsequent repair— only the latter confers the inversion the potential to have dosage effects and/or to generate advantageous chimeric genes. In Drosophila subobscura, there is ample evidence for the adaptive character of its chromosomal polymorphism, with an important contribution of some warm-climate arrangements such as E1+2+9+12. Here, we have characterized the breakpoints of inversion E12 and established that it originated through the staggered-break mechanism like four of the five inversions of D. subobscura previously studied. This mechanism that also predominates in the D. melanogaster lineage might be prevalent in the Sophophora subgenus and contribute to the adaptive character of the polymorphic and fixed inversions of its species. Finally, we have shown that the D. subobscura inversion breakpoint regions have generally been disrupted by additional structural changes occurred at different time scales. PMID:27470196
Puerma, Eva; Orengo, Dorcas J; Aguadé, Montserrat
2016-01-01
Chromosomal inversions can contribute to the adaptation of organisms to their environment by capturing particular advantageous allelic combinations of a set of genes included in the inverted fragment and also by advantageous functional changes due to the inversion process itself that might affect not only the expression of flanking genes but also their dose and structure. Of the two mechanisms originating inversions -ectopic recombination, and staggered double-strand breaks and subsequent repair- only the latter confers the inversion the potential to have dosage effects and/or to generate advantageous chimeric genes. In Drosophila subobscura, there is ample evidence for the adaptive character of its chromosomal polymorphism, with an important contribution of some warm-climate arrangements such as E1+2+9+12. Here, we have characterized the breakpoints of inversion E12 and established that it originated through the staggered-break mechanism like four of the five inversions of D. subobscura previously studied. This mechanism that also predominates in the D. melanogaster lineage might be prevalent in the Sophophora subgenus and contribute to the adaptive character of the polymorphic and fixed inversions of its species. Finally, we have shown that the D. subobscura inversion breakpoint regions have generally been disrupted by additional structural changes occurred at different time scales. PMID:27470196
NASA Technical Reports Server (NTRS)
Hedland, D. A.; Degonia, P. K.
1974-01-01
The RAE-1 spacecraft inversion performed October 31, 1972 is described based upon the in-orbit dynamical data in conjunction with results obtained from previously developed computer simulation models. The computer simulations used are predictive of the satellite dynamics, including boom flexing, and are applicable during boom deployment and retraction, inter-phase coast periods, and post-deployment operations. Attitude data, as well as boom tip data, were analyzed in order to obtain a detailed description of the dynamical behavior of the spacecraft during and after the inversion. Runs were made using the computer model and the results were analyzed and compared with the real time data. Close agreement between the actual recorded spacecraft attitude and the computer simulation results was obtained.
Simplified, inverse, ejector design tool
NASA Technical Reports Server (NTRS)
Dechant, Lawrence J.
1993-01-01
A simple lumped parameter based inverse design tool has been developed which provides flow path geometry and entrainment estimates subject to operational, acoustic, and design constraints. These constraints are manifested through specification of primary mass flow rate or ejector thrust, fully-mixed exit velocity, and static pressure matching. Fundamentally, integral forms of the conservation equations coupled with the specified design constraints are combined to yield an easily invertible linear system in terms of the flow path cross-sectional areas. Entrainment is computed by back substitution. Initial comparison with experimental and analogous one-dimensional methods show good agreement. Thus, this simple inverse design code provides an analytically based, preliminary design tool with direct application to High Speed Civil Transport (HSCT) design studies.
Broadband synthetic aperture geoacoustic inversion.
Tan, Bien Aik; Gerstoft, Peter; Yardim, Caglar; Hodgkiss, William S
2013-07-01
A typical geoacoustic inversion procedure involves powerful source transmissions received on a large-aperture receiver array. A more practical approach is to use a single moving source and/or receiver in a low signal to noise ratio (SNR) setting. This paper uses single-receiver, broadband, frequency coherent matched-field inversion and exploits coherently repeated transmissions to improve estimation of the geoacoustic parameters. The long observation time creates a synthetic aperture due to relative source-receiver motion. This approach is illustrated by studying the transmission of multiple linear frequency modulated (LFM) pulses which results in a multi-tonal comb spectrum that is Doppler sensitive. To correlate well with the measured field across a receiver trajectory and to incorporate transmission from a source trajectory, waveguide Doppler and normal mode theory is applied. The method is demonstrated with low SNR, 100-900 Hz LFM pulse data from the Shallow Water 2006 experiment. PMID:23862809
Southern California Adjoint Source Inversions
NASA Astrophysics Data System (ADS)
Tromp, J.; Kim, Y.
2007-12-01
Southern California Centroid-Moment Tensor (CMT) solutions with 9 components (6 moment tensor elements, latitude, longitude, and depth) are sought to minimize a misfit function computed from waveform differences. The gradient of a misfit function is obtained based upon two numerical simulations for each earthquake: one forward calculation for the southern California model, and an adjoint calculation that uses time-reversed signals at the receivers. Conjugate gradient and square-root variable metric methods are used to iteratively improve the earthquake source model while reducing the misfit function. The square-root variable metric algorithm has the advantage of providing a direct approximation to the posterior covariance operator. We test the inversion procedure by perturbing each component of the CMT solution, and see how the algorithm converges. Finally, we demonstrate full inversion capabilities using data for real Southern California earthquakes.
Momentum resolution in inverse photoemission
Zumbülte, A.; Schmidt, A. B.; Donath, M.
2015-01-15
We present a method to determine the electron beam divergence, and thus the momentum resolution, of an inverse-photoemission setup directly from a series of spectra measured on Cu(111). Simulating these spectra with different beam divergences shows a distinct influence of the divergence on the appearance of the Shockley surface state. Upon crossing the Fermi level, its rise in intensity can be directly linked with the beam divergence. A comparison of measurement and simulation enables us to quantify the momentum resolution independent of surface quality, energy resolution, and experimental geometry. With spin resolution, a single spectrum taken around the Fermi momentum of a spin-split surface state, e.g., on Au(111), is sufficient to derive the momentum resolution of an inverse-photoemission setup.
Inverse psoriasis treated with ustekinumab.
Campos, Manuel António; Varela, Paulo; Baptista, Armando; Moreira, Ana Isabel
2016-01-01
Inverse psoriasis is characterised by the involvement of flexural skin folds. This form of psoriasis has distinct clinical and therapeutic features. This report refers to the case of a 48-year-old Caucasian man who was observed in our department, with a clinically and biopsy proven diagnosis of inverse psoriasis. For 2 years, the patient was treated with different combinations of corticosteroids, vitamin D analogues and methotrexate, with no satisfactory response. Given the lack of a clinical response and comorbidities, latent tuberculosis was excluded, and we started treatment with ustekinumab. We chose this biological agent because the patient was a long-distance truck driver and refused the possibility of autoinjections. The patient underwent three ustekinumab injections, which resulted in significant improvement of pruritus, erythaematous lesions and quality of life. PMID:27222277
Techniques in Doppler gravity inversion
NASA Technical Reports Server (NTRS)
Phillips, R. J.
1974-01-01
The types of Doppler gravity data available for local as opposed to planetwide geophysical modeling are reviewed. Those gravity fields that are determined dynamically in orbit determination programs yield a smoothed representation of the local gravity field that may be used for quantitative modeling. An estimate of the difference between smoothed and true fields can be considered as a noise limitation in generating local gravity models. A nonlinear inversion for the geometry, depth, and density of the Mare Serenitatis mascon using an ellipsoidal model yielded a global least squares minimum in horizontal dimensions, depth, and thickness-density contrast product. It was subsequently found, by using a linear model, that there were an infinite number of solutions corresponding to various combinations of depth and lateral inhomogeneity. Linear modeling was performed by means of generalized inverse theory.
Pyramidal inversion domain boundaries revisited
Remmele, T.; Albrecht, M.; Irmscher, K.; Fornari, R.; Strassburg, M.
2011-10-03
The structure of pyramidal inversion domain boundaries in GaN:Mg was investigated by aberration corrected transmission electron microscopy. The analysis shows the upper (0001) boundary to consist of a single Mg layer inserted between polarity inverted GaN layers in an abcab stacking. The Mg bound in these defects is at least one order of magnitude lower than the chemical Mg concentration. Temperature dependent Hall effect measurements show that up to 27% of the Mg acceptors is electrically compensated.
Stochastic inversion by ray continuation
Haas, A.; Viallix
1989-05-01
The conventional tomographic inversion consists in minimizing residuals between measured and modelled traveltimes. The process tends to be unstable and some additional constraints are required to stabilize it. The stochastic formulation generalizes the technique and sets it on firmer theoretical bases. The Stochastic Inversion by Ray Continuation (SIRC) is a probabilistic approach, which takes a priori geological information into account and uses probability distributions to characterize data correlations and errors. It makes it possible to tie uncertainties to the results. The estimated parameters are interval velocities and B-spline coefficients used to represent smoothed interfaces. Ray tracing is done by a continuation technique between source and receives. The ray coordinates are computed from one path to the next by solving a linear system derived from Fermat's principle. The main advantages are fast computations, accurate traveltimes and derivatives. The seismic traces are gathered in CMPs. For a particular CMP, several reflecting elements are characterized by their time gradient measured on the stacked section, and related to a mean emergence direction. The program capabilities are tested on a synthetic example as well as on a field example. The strategy consists in inverting the parameters for one layer, then for the next one down. An inversion step is divided in two parts. First the parameters for the layer concerned are inverted, while the parameters for the upper layers remain fixed. Then all the parameters are reinverted. The velocity-depth section computed by the program together with the corresponding errors can be used directly for the interpretation, as an initial model for depth migration or for the complete inversion program under development.
NASA Technical Reports Server (NTRS)
2005-01-01
The goal of this research is to develop and demonstrate innovative adaptive seal technologies that can lead to dramatic improvements in engine performance, life, range, and emissions, and enhance operability for next generation gas turbine engines. This work is concentrated on the development of self-adaptive clearance control systems for gas turbine engines. Researchers have targeted the high-pressure turbine (HPT) blade tip seal location for following reasons: Current active clearance control (ACC) systems (e.g., thermal case-cooling schemes) cannot respond to blade tip clearance changes due to mechanical, thermal, and aerodynamic loads. As such they are prone to wear due to the required tight running clearances during operation. Blade tip seal wear (increased clearances) reduces engine efficiency, performance, and service life. Adaptive sealing technology research has inherent impact on all envisioned 21st century propulsion systems (e.g. distributed vectored, hybrid and electric drive propulsion concepts).
Adaptive method of realizing natural gradient learning for multilayer perceptrons.
Amari, S; Park, H; Fukumizu, K
2000-06-01
The natural gradient learning method is known to have ideal performances for on-line training of multilayer perceptrons. It avoids plateaus, which give rise to slow convergence of the backpropagation method. It is Fisher efficient, whereas the conventional method is not. However, for implementing the method, it is necessary to calculate the Fisher information matrix and its inverse, which is practically very difficult. This article proposes an adaptive method of directly obtaining the inverse of the Fisher information matrix. It generalizes the adaptive Gauss-Newton algorithms and provides a solid theoretical justification of them. Simulations show that the proposed adaptive method works very well for realizing natural gradient learning. PMID:10935719
Online monitoring of polymerization reactions in inverse emulsions.
Alb, Alina M; Farinato, Ray; Calbick, Joe; Reed, Wayne F
2006-01-17
Automatic continuous online monitoring of polymerization reactions (ACOMP) was adapted to the monitoring of acrylamide polymerization in inverse emulsions. This is the first application of ACOMP to heterogeneous phase polymerization. The conversion and reduced viscosity were monitored by continuously inverting and diluting the emulsion phase using a small reactor sample stream and a breaker surfactant solution, followed by UV absorption and viscometric detection. This inversion into a stable portion of the polymer/surfactant phase diagram is accomplished in tens of seconds, yielding dilute solutions containing acrylamide (Aam), polyacrylamide (PA), oil droplets, and small quantities of surfactant, initiator and other debris, and low molecular weight compounds. After establishing the means of making ACOMP measurements, a first application of the method is made to resolving some of the kinetic issues involved in emulsion polymerization, including the evolution of molecular mass, and the simultaneous action of an "intrinsic" initiator and an added chemical initiator. PMID:16401138
A reduced basis Landweber method for nonlinear inverse problems
NASA Astrophysics Data System (ADS)
Garmatter, Dominik; Haasdonk, Bernard; Harrach, Bastian
2016-03-01
We consider parameter identification problems in parametrized partial differential equations (PDEs). These lead to nonlinear ill-posed inverse problems. One way of solving them is using iterative regularization methods, which typically require numerous amounts of forward solutions during the solution process. In this article we consider the nonlinear Landweber method and couple it with the reduced basis method as a model order reduction technique in order to reduce the overall computational time. In particular, we consider PDEs with a high-dimensional parameter space, which are known to pose difficulties in the context of reduced basis methods. We present a new method that is able to handle such high-dimensional parameter spaces by combining the nonlinear Landweber method with adaptive online reduced basis updates. It is then applied to the inverse problem of reconstructing the conductivity in the stationary heat equation.
Restricted Complexity Framework for Nonlinear Adaptive Control in Complex Systems
Williams, Rube B.
2004-02-04
Control law adaptation that includes implicit or explicit adaptive state estimation, can be a fundamental underpinning for the success of intelligent control in complex systems, particularly during subsystem failures, where vital system states and parameters can be impractical or impossible to measure directly. A practical algorithm is proposed for adaptive state filtering and control in nonlinear dynamic systems when the state equations are unknown or are too complex to model analytically. The state equations and inverse plant model are approximated by using neural networks. A framework for a neural network based nonlinear dynamic inversion control law is proposed, as an extrapolation of prior developed restricted complexity methodology used to formulate the adaptive state filter. Examples of adaptive filter performance are presented for an SSME simulation with high pressure turbine failure to support extrapolations to adaptive control problems.
Restricted Complexity Framework for Nonlinear Adaptive Control in Complex Systems
NASA Astrophysics Data System (ADS)
Williams, Rube B.
2004-02-01
Control law adaptation that includes implicit or explicit adaptive state estimation, can be a fundamental underpinning for the success of intelligent control in complex systems, particularly during subsystem failures, where vital system states and parameters can be impractical or impossible to measure directly. A practical algorithm is proposed for adaptive state filtering and control in nonlinear dynamic systems when the state equations are unknown or are too complex to model analytically. The state equations and inverse plant model are approximated by using neural networks. A framework for a neural network based nonlinear dynamic inversion control law is proposed, as an extrapolation of prior developed restricted complexity methodology used to formulate the adaptive state filter. Examples of adaptive filter performance are presented for an SSME simulation with high pressure turbine failure to support extrapolations to adaptive control problems.