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…
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.
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.
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
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.
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.
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.
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.
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
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.
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
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.}
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
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.
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.
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
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.
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
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…
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.
On the Magic Square and Inverse
ERIC Educational Resources Information Center
Elzaidi, S. M.
2005-01-01
In this note, we give a method for finding the inverse of a three by three magic square matrix without using the usual methods for finding the inverse of a matrix. Also we give a method for finding the inverse of a three by three magic square matrix whose entries are also matrices. By using these ideas, we can construct large matrices whose…
Fast Computation of the Inverse CMH Model
NASA Technical Reports Server (NTRS)
Patel, Umesh D.; Torre, Edward Della; Day, John H. (Technical Monitor)
2001-01-01
A fast computational method based on differential equation approach for inverse DOK model has been extended for the inverse CMH model. Also, a cobweb technique for calculating the inverse CMH model is also presented. The two techniques are differed from the point of view of flexibility and computation time.
Dynamically consistent Jacobian inverse for mobile manipulators
NASA Astrophysics Data System (ADS)
Ratajczak, Joanna; Tchoń, Krzysztof
2016-06-01
By analogy to the definition of the dynamically consistent Jacobian inverse for robotic manipulators, we have designed a dynamically consistent Jacobian inverse for mobile manipulators built of a non-holonomic mobile platform and a holonomic on-board manipulator. The endogenous configuration space approach has been exploited as a source of conceptual guidelines. The new inverse guarantees a decoupling of the motion in the operational space from the forces exerted in the endogenous configuration space and annihilated by the dual Jacobian inverse. A performance study of the new Jacobian inverse as a tool for motion planning is presented.
The continuation inverse problem revisited
NASA Astrophysics Data System (ADS)
Huestis, Stephen P.
1998-06-01
The non-uniqueness of the continuation of a finite collection of harmonic potential field data to a level surface in the source-free region forces its treatment as an inverse problem. A formalism is proposed for the construction of continuation functions which are extremal by various measures. The problem is cast in such a form that the inverse problem solution is the potential function on the lowest horizontal surface above all sources, serving as the boundary function for the Dirichlet problem in the upper half-plane. The desired continuation, at the higher level of interest, must then be in the range of the upward continuation operator acting on this boundary function, rather than being allowed the full freedom of itself being part of a Dirichlet problem boundary function. Extremal solutions minimize non-linear functionals of the continuation function, which are re-expressed as different functionals of the boundary function. A crux of the method is that there is no essential distinction between the upward and downward continuation inverse problems to levels above or below data locations. Casting the optimization as a Lagrange multiplier problem leads to an integral equation for the boundary function, which is readily solved in the Fourier domain for a certain class of functionals. The desired extremal continuation is then given by upward continuation. It is found that for some functionals, application of the Lagrange multiplier theorem requires a further restriction on the set of allowable boundary functions: bandlimitedness is a natural choice for the continuation problem. With this imposition, the theory is developed in detail for semi-norm functionals penalizing departure from a constant potential, in the 2-norm and Sobelev norm senses, and illustrated by application for a small synthetic Deep Tow magnetic field data set.
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.
Parallel inverse iteration with reorthogonalization
Fann, G.I.; Littlefield, R.J.
1993-03-01
A parallel method for finding orthogonal eigenvectors of real symmetric tridiagonal is described. The method uses inverse iteration with repeated Modified Gram-Schmidt (MGS) reorthogonalization of the unconverged iterates for clustered eigenvalues. This approach is more parallelizable than reorthogonalizing against fully converged eigenvectors, as is done by LAPACK's current DSTEIN routine. The new method is found to provide accuracy and speed comparable to DSTEIN's and to have good parallel scalability even for matrices with large clusters of eigenvalues. We present al results for residual and orthogonality tests, plus timings on IBM RS/6000 (sequential) and Intel Touchstone DELTA (parallel) computers.
Parallel inverse iteration with reorthogonalization
Fann, G.I.; Littlefield, R.J.
1993-03-01
A parallel method for finding orthogonal eigenvectors of real symmetric tridiagonal is described. The method uses inverse iteration with repeated Modified Gram-Schmidt (MGS) reorthogonalization of the unconverged iterates for clustered eigenvalues. This approach is more parallelizable than reorthogonalizing against fully converged eigenvectors, as is done by LAPACK`s current DSTEIN routine. The new method is found to provide accuracy and speed comparable to DSTEIN`s and to have good parallel scalability even for matrices with large clusters of eigenvalues. We present al results for residual and orthogonality tests, plus timings on IBM RS/6000 (sequential) and Intel Touchstone DELTA (parallel) computers.
Rotation and inversion in nitrosamines
NASA Astrophysics Data System (ADS)
Kirste, Karl; Rademacher, Paul
1981-04-01
Geometry optimizations of the ground states as well as of the transition states for internal rotation and inversion have been performed by the semiempirical MNDO method for dimethyl nitrosamine (1), perfluordimethyl nitrosamine (2), N-nitroso aziridine (3), and N-nitroso azetidine (4). It was found that the potential barrier to internal rotation about the N-N bond is always of lower energy than that to inversion on the nitroso nitrogen. While the ground states tend to adopt structures which enable mesomerism, the lowest transition state is characterized by a pyramidal sp3-hybridized amino nitrogen. In accordance with experimental results the low barriers to rotation of 2 (7.96 kcal mol -1), 3 (3.38 kcal mol -1) and 4 (9.97 kcal mol -1) in comparison with 1 (12.54 kcal mol -1) indicate that in donor-acceptor molecules the transfer of charge can be limited by electronic and stereochemical effects. In particular, the equivalence of the α-methylene hydrogens which was observed in the NMR-spectrum of 3 is due to unhindered rotation and ring inveirsion.
MOSES Inversions using Multiresolution SMART
NASA Astrophysics Data System (ADS)
Rust, Thomas; Fox, Lewis; Kankelborg, Charles; Courrier, Hans; Plovanic, Jacob
2014-06-01
We present improvements to the SMART inversion algorithm for the MOSES imaging spectrograph. MOSES, the Multi-Order Solar EUV Spectrograph, is a slitless extreme ultraviolet spectrograph designed to measure cotemporal narrowband spectra over a wide field of view via tomographic inversion of images taken at three orders of a concave diffraction grating. SMART, the Smooth Multiplicative Algebraic Reconstruction Technique, relies on a global chi squared goodness of fit criterion, which enables overfit and underfit regions to "balance out" when judging fit quality. "Good" reconstructions show poor fits at some positions and length scales. Here we take a multiresolution approach to SMART, applying corrections to the reconstruction at positions and scales where correction is warranted based on the noise. The result is improved fit residuals that more closely resemble the expected noise in the images. Within the multiresolution framework it is also easy to include a regularized deconvolution of the instrument point spread functions, which we do. Different point spread functions among MOSES spectral orders results in spurious doppler shifts in the reconstructions, most notable near bright compact emission. We estimate the point spread funtions from the data. Deconvolution is done using the Richardson-Lucy method, which is algorithmically similar to SMART. Regularization results from only correcting the reconstruction at positions and scales where correction is warranted based on the noise. We expect the point spread function deconvolution to increase signal to noise and reduce systematic error in MOSES reconstructions.
Multi-Skip Tomographic Inversion
NASA Astrophysics Data System (ADS)
Volker, Arno; Bloom, Joost; Lorenz, Maarten
2011-06-01
Inspection of corrosion at pipe support locations is difficult because of accessibility limitations. Recently a screening technique has been developed called Multi-Skip ultrasonics. The method utilizes a pitch-catch set-up. Shear waves are transmitted that reflect multiple times in the pipe wall, from which integral wall thickness information is obtained. The method turns out to be very sensitive in detecting the presence of wall loss, but it turns out to be difficult to determine the extent of the wall loss. If the extent is not known, only a conservative estimate of depth can be derived from the Multi-Skip signals, because of the accumulative nature of the change in arrival time due to wall loss. Multi-Skip tomography appears to be a promising method in addition to Multi-Skip screening as a follow-up inspection technique. It uses full wave field inversion to determine a wall thickness profile at a particular location of the pipe on the support. As with the Multi-Skip screening method, Multi-Skip tomography is applied with the transmitter and receiver on both sides of the pipe support location and waves traveling in the axial pipe direction. The wave field inversion consists of a forward modeling step that predicts the measured wave field after which an iterative comparison process with the actually measured wave field results in an estimate of the wall thickness profile under the support.
Recent Advances in Geostatistical Inversion
NASA Astrophysics Data System (ADS)
Kitanidis, P. K.
2011-12-01
Inverse problems are common in hydrologic applications, such as in subsurface imaging which is the identification of parameters characterizing geologic formations from hydrologic and geophysical observations. In such problems, the data do not suffice to constrain the solution to be unique. The geostatistical approach utilizes probability theory and statistical inference to assimilate data and information about structure and to explore the range of possible solutions in a systematic way. This is a progress report on recent advances in terms of formulation and computational methods. The standard implementation of the geostatistical approach to the inverse method is computationally very expensive when there are millions of unknowns and hundreds of thousands of observations, as is the case in fusing data from many sources in hydrogeology. However, depending on the specific problem, alternative formulations and numerical methods can reduce the computational problem quite dramatically. One approach can utilize formulations that involve matrices with a very high degree of sparsity combined with indirect methods of solution and strategies that minimize the number of required forward runs. The potential for this method is illustrated with an application to transient hydraulic tomography. Another approach speeds up matrix-vector multiplications by utilizing hierarchical sparsity in commonly encountered matrices, particularly prior covariance matrices. A couple of examples show how the computational cost scales with the size of the problem (number of observations and unknowns) in conventional and newer methods. Yet another fruitful approach is to rely on a large number of realizations to represent ensembles of solutions and we illustrate this approach in fusing data from two different geophysical methods. In all of these approaches, utilizing parallel processors and mixed CPU/GPU programming can significantly reduce the computational cost and make it possible to solve very
Wake Vortex Inverse Model User's Guide
NASA Technical Reports Server (NTRS)
Lai, David; Delisi, Donald
2008-01-01
NorthWest Research Associates (NWRA) has developed an inverse model for inverting landing aircraft vortex data. The data used for the inversion are the time evolution of the lateral transport position and vertical position of both the port and starboard vortices. The inverse model performs iterative forward model runs using various estimates of vortex parameters, vertical crosswind profiles, and vortex circulation as a function of wake age. Forward model predictions of lateral transport and altitude are then compared with the observed data. Differences between the data and model predictions guide the choice of vortex parameter values, crosswind profile and circulation evolution in the next iteration. Iterations are performed until a user-defined criterion is satisfied. Currently, the inverse model is set to stop when the improvement in the rms deviation between the data and model predictions is less than 1 percent for two consecutive iterations. The forward model used in this inverse model is a modified version of the Shear-APA model. A detailed description of this forward model, the inverse model, and its validation are presented in a different report (Lai, Mellman, Robins, and Delisi, 2007). This document is a User's Guide for the Wake Vortex Inverse Model. Section 2 presents an overview of the inverse model program. Execution of the inverse model is described in Section 3. When executing the inverse model, a user is requested to provide the name of an input file which contains the inverse model parameters, the various datasets, and directories needed for the inversion. A detailed description of the list of parameters in the inversion input file is presented in Section 4. A user has an option to save the inversion results of each lidar track in a mat-file (a condensed data file in Matlab format). These saved mat-files can be used for post-inversion analysis. A description of the contents of the saved files is given in Section 5. An example of an inversion input
Dark radiative inverse seesaw mechanism
NASA Astrophysics Data System (ADS)
Ahriche, Amine; Boucenna, Sofiane M.; Nasri, Salah
2016-04-01
We present a minimal model that simultaneously accounts for neutrino masses and the origin of dark matter (DM) and where the electroweak phase transition is strong enough to allow for electroweak baryogenesis. The Standard Model is enlarged with a Majorana fermion, three generations of chiral fermion pairs, and a single complex scalar that plays a central role in DM production and phenomenology, neutrino masses, and the strength of the phase transition. All the new fields are singlets under the SM gauge group. Neutrino masses are generated via a new variant of the radiative inverse seesaw mechanism, where the required small mass term is generated via loops involving DM and no large hierarchy is assumed among the mass scales. The model offers all the advantage of low-scale neutrino mass models as well as a viable dark matter candidate that is testable with direct detection experiments.
Hierarchical Bayesian inverse reinforcement learning.
Choi, Jaedeug; Kim, Kee-Eung
2015-04-01
Inverse reinforcement learning (IRL) is the problem of inferring the underlying reward function from the expert's behavior data. The difficulty in IRL mainly arises in choosing the best reward function since there are typically an infinite number of reward functions that yield the given behavior data as optimal. Another difficulty comes from the noisy behavior data due to sub-optimal experts. We propose a hierarchical Bayesian framework, which subsumes most of the previous IRL algorithms as well as models the sub-optimality of the expert's behavior. Using a number of experiments on a synthetic problem, we demonstrate the effectiveness of our approach including the robustness of our hierarchical Bayesian framework to the sub-optimal expert behavior data. Using a real dataset from taxi GPS traces, we additionally show that our approach predicts the driving behavior with a high accuracy. PMID:25291805
Applications of matrix inversion tomosynthesis
NASA Astrophysics Data System (ADS)
Warp, Richard J.; Godfrey, Devon J.; Dobbins, James T., III
2000-04-01
The improved image quality and characteristics of new flat- panel x-ray detectors have renewed interest in advanced algorithms such as tomosynthesis. Digital tomosynthesis is a method of acquiring and reconstructing a three-dimensional data set with limited-angle tube movement. Historically, conventional tomosynthesis reconstruction has suffered contamination of the planes of interest by blurred out-of- plane structures. This paper focuses on a Matrix Inversion Tomosynthesis (MITS) algorithm to remove unwanted blur from adjacent planes. The algorithm uses a set of coupled equations to solve for the blurring function in each reconstructed plane. This paper demonstrates the use of the MITS algorithm in three imaging applications: small animal microscopy, chest radiography, and orthopedics. The results of the MITS reconstruction process demonstrate an improved reduction of blur from out-of-plane structures when compared to conventional tomosynthesis. We conclude that the MITS algorithm holds potential in a variety of applications to improve three-dimensional image reconstruction.
Forward model nonlinearity versus inverse model nonlinearity
Mehl, S.
2007-01-01
The issue of concern is the impact of forward model nonlinearity on the nonlinearity of the inverse model. The question posed is, "Does increased nonlinearity in the head solution (forward model) always result in increased nonlinearity in the inverse solution (estimation of hydraulic conductivity)?" It is shown that the two nonlinearities are separate, and it is not universally true that increased forward model nonlinearity increases inverse model nonlinearity. ?? 2007 National Ground Water Association.
NICOLE: NLTE Stokes Synthesis/Inversion Code
NASA Astrophysics Data System (ADS)
Socas-Navarro, H.
2015-08-01
NICOLE, written in Fortran 90, seeks the model atmosphere that provides the best fit to the Stokes profiles (in a least-squares sense) of an arbitrary number of simultaneously-observes spectral lines from solar/stellar atmospheres. The inversion core used for the development of NICOLE is the LORIEN engine (the Lovely Reusable Inversion ENgine), which combines the SVD technique with the Levenberg-Marquardt minimization method to solve the inverse problem.
Superresolution imaging from nonlinear inverse scattering
NASA Astrophysics Data System (ADS)
Ritter, R. Shane; Fiddy, M. A.
2015-09-01
Inverse scattering algorithms typically rely on weak scattering approximations and the inversion of far field data on an Ewald sphere. This, in turn, fixes the resolution of the computed image. However, it has long been observed that when multiple scattering occurs in a strongly interacting object, and a nonlinear inversion method is employed to image it, the resulting image can reveal subwavelength resolution. We have observed this phenomenon using a cepstral filtering approach and characterize it more fully here.
Inversion methods for interpretation of asteroid lightcurves
NASA Technical Reports Server (NTRS)
Kaasalainen, Mikko; Lamberg, L.; Lumme, K.
1992-01-01
We have developed methods of inversion that can be used in the determination of the three-dimensional shape or the albedo distribution of the surface of a body from disk-integrated photometry, assuming the shape to be strictly convex. In addition to the theory of inversion methods, we have studied the practical aspects of the inversion problem and applied our methods to lightcurve data of 39 Laetitia and 16 Psyche.
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
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
Laterally constrained inversion for CSAMT data interpretation
NASA Astrophysics Data System (ADS)
Wang, Ruo; Yin, Changchun; Wang, Miaoyue; Di, Qingyun
2015-10-01
Laterally constrained inversion (LCI) has been successfully applied to the inversion of dc resistivity, TEM and airborne EM data. However, it hasn't been yet applied to the interpretation of controlled-source audio-frequency magnetotelluric (CSAMT) data. In this paper, we apply the LCI method for CSAMT data inversion by preconditioning the Jacobian matrix. We apply a weighting matrix to Jacobian to balance the sensitivity of model parameters, so that the resolution with respect to different model parameters becomes more uniform. Numerical experiments confirm that this can improve the convergence of the inversion. We first invert a synthetic dataset with and without noise to investigate the effect of LCI applications to CSAMT data, for the noise free data, the results show that the LCI method can recover the true model better compared to the traditional single-station inversion; and for the noisy data, the true model is recovered even with a noise level of 8%, indicating that LCI inversions are to some extent noise insensitive. Then, we re-invert two CSAMT datasets collected respectively in a watershed and a coal mine area in Northern China and compare our results with those from previous inversions. The comparison with the previous inversion in a coal mine shows that LCI method delivers smoother layer interfaces that well correlate to seismic data, while comparison with a global searching algorithm of simulated annealing (SA) in a watershed shows that though both methods deliver very similar good results, however, LCI algorithm presented in this paper runs much faster. The inversion results for the coal mine CSAMT survey show that a conductive water-bearing zone that was not revealed by the previous inversions has been identified by the LCI. This further demonstrates that the method presented in this paper works for CSAMT data inversion.
NASA Astrophysics Data System (ADS)
Nielsen, O. F.; Ploug, C.; Mendoza, J. A.; Martínez, K.
2009-05-01
The need for increaseding accuracy and reduced ambiguities in the inversion results has resulted in focus on the development of more advanced inversion methods of geophysical data. Over the past few years more advanced inversion techniques have been developed to improve the results. Real 3D-inversion is time consuming and therefore often not the best solution in a cost-efficient perspective. This has motivated the development of 3D constrained inversions, where 1D-models are constrained in 3D, also known as a Spatial Constrained Inversion (SCI). Moreover, inversion of several different data types in one inversion has been developed, known as Mutually Constrained Inversion (MCI). In this paper a presentation of a Spatial Mutually Constrained Inversion method (SMCI) is given. This method allows 1D-inversion applied to different geophysical datasets and geological information constrained in 3D. Application of two or more types of geophysical methods in the inversion has proved to reduce the equivalence problem and to increase the resolution in the inversion results. The use of geological information from borehole data or digital geological models can be integrated in the inversion. In the SMCI, a 1D inversion code is used to model soundings that are constrained in three dimensions according to their relative position in space. This solution enhances the accuracy of the inversion and produces distinct layers thicknesses and resistivities. It is very efficient in the mapping of a layered geology but still also capable of mapping layer discontinuities that are, in many cases, related to fracturing and faulting or due to valley fills. Geological information may be included in the inversion directly or used only to form a starting model for the individual soundings in the inversion. In order to show the effectiveness of the method, examples are presented from both synthetic data and real data. The examples include DC-soundings as well as land-based and airborne TEM
An exact inverse method for subsonic flows
NASA Technical Reports Server (NTRS)
Daripa, Prabir
1988-01-01
A new inverse method for the aerodynamic design of airfoils is presented for subcritical flows. The pressure distribution in this method can be prescribed as a function of the arclength of the still unknown body. It is shown that this inverse problem is mathematically equivalent to solving only one nonlinear boundary value problem subject to known Dirichlet data on the boundary.
Inversion and approximation of Laplace transforms
NASA Technical Reports Server (NTRS)
Lear, W. M.
1980-01-01
A method of inverting Laplace transforms by using a set of orthonormal functions is reported. As a byproduct of the inversion, approximation of complicated Laplace transforms by a transform with a series of simple poles along the left half plane real axis is shown. The inversion and approximation process is simple enough to be put on a programmable hand calculator.
Package inspection using inverse diffraction
NASA Astrophysics Data System (ADS)
McAulay, Alastair D.
2008-08-01
More efficient cost-effective hand-held methods of inspecting packages without opening them are in demand for security. Recent new work in TeraHertz sources,1 millimeter waves, presents new possibilities. Millimeter waves pass through cardboard and styrofoam, common packing materials, and also pass through most materials except those with high conductivity like metals which block light and are easily spotted. Estimating refractive index along the path of the beam through the package from observations of the beam passing out of the package provides the necessary information to inspect the package and is a nonlinear problem. So we use a generalized linear inverse technique that we first developed for finding oil by reflection in geophysics.2 The computation assumes parallel slices in the packet of homogeneous material for which the refractive index is estimated. A beam is propagated through this model in a forward computation. The output is compared with the actual observations for the package and an update computed for the refractive indices. The loop is repeated until convergence. The approach can be modified for a reflection system or to include estimation of absorption.
Inverse problem in radionuclide transport
Yu, C.
1988-01-01
The disposal of radioactive waste must comply with the performance objectives set forth in 10 CFR 61 for low-level waste (LLW) and 10 CFR 60 for high-level waste (HLW). To determine probable compliance, the proposed disposal system can be modeled to predict its performance. One of the difficulties encountered in such a study is modeling the migration of radionuclides through a complex geologic medium for the long term. Although many radionuclide transport models exist in the literature, the accuracy of the model prediction is highly dependent on the model parameters used. The problem of using known parameters in a radionuclide transport model to predict radionuclide concentrations is a direct problem (DP); whereas the reverse of DP, i.e., the parameter identification problem of determining model parameters from known radionuclide concentrations, is called the inverse problem (IP). In this study, a procedure to solve IP is tested, using the regression technique. Several nonlinear regression programs are examined, and the best one is recommended. 13 refs., 1 tab.
Inverse magnetic/shear catalysis
NASA Astrophysics Data System (ADS)
McInnes, Brett
2016-05-01
It is well known that very large magnetic fields are generated when the Quark-Gluon Plasma is formed during peripheral heavy-ion collisions. Lattice, holographic, and other studies strongly suggest that these fields may, for observationally relevant field values, induce "inverse magnetic catalysis", signalled by a lowering of the critical temperature for the chiral/deconfinement transition. The theoretical basis of this effect has recently attracted much attention; yet so far these investigations have not included another, equally dramatic consequence of the peripheral collision geometry: the QGP acquires a large angular momentum vector, parallel to the magnetic field. Here we use holographic techniques to argue that the angular momentum can also, independently, have an effect on transition temperatures, and we obtain a rough estimate of the relative effects of the presence of both a magnetic field and an angular momentum density. We find that the shearing angular momentum reinforces the effect of the magnetic field at low values of the baryonic chemical potential, but that it can actually decrease that effect at high chemical potentials.
Inverse Common-Reflection-Surface
NASA Astrophysics Data System (ADS)
Perroud, H.; Tygel, M.; Freitas, L.
2010-12-01
The Common-Reflection-Surface (CRS) stack method is a powerful tool to produce high-quality stacked images of multicoverage seismic data. As a result of the CRS stack, not only a stacked section, but also a number of attributes defined at each point of that section, are produced. In this way, one can think of the CRS stack method as a transformation from data space to attribute space. Being a purely kinematic method, the CRS stack lacks amplitude information that can be useful for many purposes. Here we propose to fill this gap by means of a combined use of a zero-offset section (that could be a short-offset or amplitude-corrected stacked section) and common midpoint gather. We present an algorithm for an inverse CRS transformation, namely one that (approximately) transforms the CRS attributes back to data space. First synthetic tests provide satisfying results for the two simple cases of single dipping-plane and single circular reflectors with a homogeneous overburden, and provide estimates of the range of applicability, in both midpoint and offset directions. We further present an application for interpolating missing traces in a near-surface, high-resolution seismic experiment, conducted in the alluvial plain of the river Gave de Pau, near Assat, southern France, showing its ability to build coherent signals, where recording was not available. A somewhat unexpected good feature of the algorithm, is that it seems capable to reconstruct signals even in muted parts of the section.
Inverse Free Electron Laser accelerator
NASA Astrophysics Data System (ADS)
Fisher, A.; Gallardo, J.; Vansteenbergen, 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(-) beam and the 10(exp 11) Watt CO2 laser beam of BNL's Accelerator Test Facility (ATF), Center for Accelerator Physics (CAP) and a wiggler. The latter element is designed as a fast excitation unit making use of alternating stacks of Vanadium Permendur (VaP) ferromagnetic laminations, periodically interspersed with conductive, nonmagnetic laminations, which act as eddy current induced field reflectors. Wiggler parameters and field distribution data will be presented for a prototype wiggler in a constant period and in a approximately 1.5 percent/cm tapered period configuration. The CO2 laser beam will be transported through the IFEL interaction region by means of a low loss, dielectric coated, rectangular waveguide. Short waveguide test sections have been constructed and have been tested using a low power CW CO2 laser. Preliminary results of guide attenuation and mode selectivity will be given, together with a discussion of the optical issues for the IFEL accelerator. The IFEL design is supported by the development and use of 1D and 3D simulation programs. The results of simulation computations, including also wiggler errors, for a single module accelerator and for a multi-module accelerator will be presented.
Inverse free electron laser accelerator
NASA Astrophysics Data System (ADS)
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- beam and the 1011 Watt CO2 laser beam of BNL's Accelerator Test Facility (ATF), Center for Accelerator Physics (CAP), and a wiggler. The latter element is designed as a fast excitation unit making use of alternating stacks of Vanadium Permendur (VaP) ferromagnetic laminations, periodically interspersed with conductive, nonmagnetic laminations, which act as eddy current induced field reflectors. Wiggler parameters and field distribution data will be presented for a prototype wiggler in a constant period and in a ≊1.5%/cm tapered period configuration. The CO2 laser beam will be transported through the IFEL interaction region by means of a low loss, dielectric coated, rectangular waveguide. Short waveguide test sections have been constructed and have been tested using a low power cw CO2 laser. Preliminary results of guide attenuation and mode selectivity will be given, together with a discussion of the optical issues for the IFEL accelerator. The IFEL design is supported by the development and use of 1D and 3D simulation programs. The results of simulation computations, including also wiggler errors, for a single module accelerator and for a multi-module accelerator will be presented.
Inversion based on computational simulations
Hanson, K.M.; Cunningham, G.S.; Saquib, S.S.
1998-09-01
A standard approach to solving inversion problems that involve many parameters uses gradient-based optimization to find the parameters that best match the data. The authors discuss enabling techniques that facilitate application of this approach to large-scale computational simulations, which are the only way to investigate many complex physical phenomena. Such simulations may not seem to lend themselves to calculation of the gradient with respect to numerous parameters. However, adjoint differentiation allows one to efficiently compute the gradient of an objective function with respect to all the variables of a simulation. When combined with advanced gradient-based optimization algorithms, adjoint differentiation permits one to solve very large problems of optimization or parameter estimation. These techniques will be illustrated through the simulation of the time-dependent diffusion of infrared light through tissue, which has been used to perform optical tomography. The techniques discussed have a wide range of applicability to modeling including the optimization of models to achieve a desired design goal.
Fast wavelet based sparse approximate inverse preconditioner
Wan, W.L.
1996-12-31
Incomplete LU factorization is a robust preconditioner for both general and PDE problems but unfortunately not easy to parallelize. Recent study of Huckle and Grote and Chow and Saad showed that sparse approximate inverse could be a potential alternative while readily parallelizable. However, for special class of matrix A that comes from elliptic PDE problems, their preconditioners are not optimal in the sense that independent of mesh size. A reason may be that no good sparse approximate inverse exists for the dense inverse matrix. Our observation is that for this kind of matrices, its inverse entries typically have piecewise smooth changes. We can take advantage of this fact and use wavelet compression techniques to construct a better sparse approximate inverse preconditioner. We shall show numerically that our approach is effective for this kind of matrices.
Phase inversion emulsification: Current understanding and applications.
Perazzo, A; Preziosi, V; Guido, S
2015-08-01
This review is addressed to the phase inversion process, which is not only a common, low-energy route to make stable emulsions for a variety of industrial products spanning from food to pharmaceuticals, but can also be an undesired effect in some applications, such as crude oil transportation in pipelines. Two main ways to induce phase inversion are described in the literature, i.e., phase inversion composition (PIC or catastrophic) and phase inversion temperature (PIT or transitional). In the former, starting from one phase (oil or water) with surfactants, the other phase is more or less gradually added until it reverts to the continuous phase. In PIT, phase inversion is driven by a temperature change without varying system composition. Given its industrial relevance and scientific challenge, phase inversion has been the subject of a number of papers in the literature, including extensive reviews. Due to the variety of applications and the complexity of the problem, most of the publications have been focused either on the phase behavior or the interfacial properties or the mixing process of the two phases. Although all these aspects are quite important in studying phase inversion and much progress has been done on this topic, a comprehensive picture is still lacking. In particular, the general mechanisms governing the inversion phenomenon have not been completely elucidated and quantitative predictions of the phase inversion point are limited to specific systems and experimental conditions. Here, we review the different approaches on phase inversion and highlight some related applications, including future and emerging perspectives. PMID:25632889
Rapid approximate inversion of airborne TEM
NASA Astrophysics Data System (ADS)
Fullagar, Peter K.; Pears, Glenn A.; Reid, James E.; Schaa, Ralf
2015-11-01
Rapid interpretation of large airborne transient electromagnetic (ATEM) datasets is highly desirable for timely decision-making in exploration. Full solution 3D inversion of entire airborne electromagnetic (AEM) surveys is often still not feasible on current day PCs. Therefore, two algorithms to perform rapid approximate 3D interpretation of AEM have been developed. The loss of rigour may be of little consequence if the objective of the AEM survey is regional reconnaissance. Data coverage is often quasi-2D rather than truly 3D in such cases, belying the need for `exact' 3D inversion. Incorporation of geological constraints reduces the non-uniqueness of 3D AEM inversion. Integrated interpretation can be achieved most readily when inversion is applied to a geological model, attributed with lithology as well as conductivity. Geological models also offer several practical advantages over pure property models during inversion. In particular, they permit adjustment of geological boundaries. In addition, optimal conductivities can be determined for homogeneous units. Both algorithms described here can operate on geological models; however, they can also perform `unconstrained' inversion if the geological context is unknown. VPem1D performs 1D inversion at each ATEM data location above a 3D model. Interpretation of cover thickness is a natural application; this is illustrated via application to Spectrem data from central Australia. VPem3D performs 3D inversion on time-integrated (resistive limit) data. Conversion to resistive limits delivers a massive increase in speed since the TEM inverse problem reduces to a quasi-magnetic problem. The time evolution of the decay is lost during the conversion, but the information can be largely recovered by constructing a starting model from conductivity depth images (CDIs) or 1D inversions combined with geological constraints if available. The efficacy of the approach is demonstrated on Spectrem data from Brazil. Both separately and in
Neurotensin inversely modulates maternal aggression.
Gammie, S C; D'Anna, K L; Gerstein, H; Stevenson, S A
2009-02-18
Neurotensin (NT) is a versatile neuropeptide involved in analgesia, hypothermia, and schizophrenia. Although NT is released from and acts upon brain regions involved in social behaviors, it has not been linked to a social behavior. We previously selected mice for high maternal aggression (maternal defense), an important social behavior that protects offspring, and found significantly lower NT expression in the CNS of highly protective females. Our current study directly tested NT's role in maternal defense. Intracerebroventricular (i.c.v.) injections of NT significantly impaired defense in terms of time aggressive and number of attacks at all doses tested (0.05, 0.1, 1.0, and 3.0 microg). Other maternal behaviors, including pup retrieval, were unaltered following NT injections (0.05 microg) relative to vehicle, suggesting specificity of NT action on defense. Further, i.c.v. injections of the NT receptor 1 (NT1) antagonist, SR 48692 (30 microg), significantly elevated maternal aggression in terms of time aggressive and attack number. To understand where NT may regulate aggression, we examined Fos following injection of either 0.1 microg NT or vehicle. Thirteen of 26 brain regions examined exhibited significant Fos increases with NT, including regions expressing NT1 and previously implicated in maternal aggression, such as lateral septum, bed nucleus of stria terminalis, paraventricular nucleus, and central amygdala. Together, our results indicate that NT inversely regulates maternal aggression and provide the first direct evidence that lowering of NT signaling can be a mechanism for maternal aggression. To our knowledge, this is the first study to directly link NT to a social behavior. PMID:19118604
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}} beam and the 10{sup 11} Watt CO{sub 2} laser beam of BNL`s Accelerator Test Facility (ATF), Center for Accelerator Physics (CAP) and a wiggler. The latter element is designed as a fast excitation unit making use of alternating stacks of Vanadium Permendur (VaP) ferromagnetic laminations, periodically interspersed with conductive, nonmagnetic laminations, which act as eddy current induced field reflectors. Wiggler parameters and field distribution data will be presented for a prototype wiggler in a constant period and in a {approximately} 1.5 %/cm tapered period configuration. The CO{sub 2} laser beam will be transported through the IFEL interaction region by means of a low loss, dielectric coated, rectangular waveguide. Short waveguide test sections have been constructed and have been tested using a low power cw CO{sub 2} laser. Preliminary results of guide attenuation and mode selectivity will be given, together with a discussion of the optical issues for the IFEL accelerator. The IFEL design is supported by the development and use of 1D and 3D simulation programs. The results of simulation computations, including also wiggler errors, for a single module accelerator and for a multi-module accelerator will be presented.
Support minimized inversion of acoustic and elastic wave scattering
Safaeinili, A.
1994-04-24
This report discusses the following topics on support minimized inversion of acoustic and elastic wave scattering: Minimum support inversion; forward modelling of elastodynamic wave scattering; minimum support linearized acoustic inversion; support minimized nonlinear acoustic inversion without absolute phase; and support minimized nonlinear elastic inversion.
Methodology for determining multilayered temperature inversions
NASA Astrophysics Data System (ADS)
Fochesatto, G. J.
2015-05-01
Temperature sounding of the atmospheric boundary layer (ABL) and lower troposphere exhibits multilayered temperature inversions specially in high latitudes during extreme winters. These temperature inversion layers are originated based on the combined forcing of local- and large-scale synoptic meteorology. At the local scale, the thermal inversion layer forms near the surface and plays a central role in controlling the surface radiative cooling and air pollution dispersion; however, depending upon the large-scale synoptic meteorological forcing, an upper level thermal inversion can also exist topping the local ABL. In this article a numerical methodology is reported to determine thermal inversion layers present in a given temperature profile and deduce some of their thermodynamic properties. The algorithm extracts from the temperature profile the most important temperature variations defining thermal inversion layers. This is accomplished by a linear interpolation function of variable length that minimizes an error function. The algorithm functionality is demonstrated on actual radiosonde profiles to deduce the multilayered temperature inversion structure with an error fraction set independently.
On the inversion-indel distance
2013-01-01
Background The inversion distance, that is the distance between two unichromosomal genomes with the same content allowing only inversions of DNA segments, can be computed thanks to a pioneering approach of Hannenhalli and Pevzner in 1995. In 2000, El-Mabrouk extended the inversion model to allow the comparison of unichromosomal genomes with unequal contents, thus insertions and deletions of DNA segments besides inversions. However, an exact algorithm was presented only for the case in which we have insertions alone and no deletion (or vice versa), while a heuristic was provided for the symmetric case, that allows both insertions and deletions and is called the inversion-indel distance. In 2005, Yancopoulos, Attie and Friedberg started a new branch of research by introducing the generic double cut and join (DCJ) operation, that can represent several genome rearrangements (including inversions). Among others, the DCJ model gave rise to two important results. First, it has been shown that the inversion distance can be computed in a simpler way with the help of the DCJ operation. Second, the DCJ operation originated the DCJ-indel distance, that allows the comparison of genomes with unequal contents, considering DCJ, insertions and deletions, and can be computed in linear time. Results In the present work we put these two results together to solve an open problem, showing that, when the graph that represents the relation between the two compared genomes has no bad components, the inversion-indel distance is equal to the DCJ-indel distance. We also give a lower and an upper bound for the inversion-indel distance in the presence of bad components. PMID:24564182
3-D inversion of magnetotelluric Phase Tensor
NASA Astrophysics Data System (ADS)
Patro, Prasanta; Uyeshima, Makoto
2010-05-01
Three-dimensional (3-D) inversion of the magnetotelluric (MT) has become a routine practice among the MT community due to progress of algorithms for 3-D inverse problems (e.g. Mackie and Madden, 1993; Siripunvaraporn et al., 2005). While availability of such 3-D inversion codes have increased the resolving power of the MT data and improved the interpretation, on the other hand, still the galvanic effects poses difficulties in interpretation of resistivity structure obtained from the MT data. In order to tackle the galvanic distortion of MT data, Caldwell et al., (2004) introduced the concept of phase tensor. They demonstrated how the regional phase information can be retrieved from the observed impedance tensor without any assumptions for structural dimension, where both the near surface inhomogeneity and the regional conductivity structures can be 3-D. We made an attempt to modify a 3-D inversion code (Siripunvaraporn et al., 2005) to directly invert the phase tensor elements. We present here the main modification done in the sensitivity calculation and then show a few synthetic studies and its application to the real data. The synthetic model study suggests that the prior model (m_0) setting is important in retrieving the true model. This is because estimation of correct induction scale length lacks in the phase tensor inversion process. Comparison between results from conventional impedance inversion and new phase tensor inversion suggests that, in spite of presence of the galvanic distortion (due to near surface checkerboard anomalies in our case), the new inverion algorithm retrieves the regional conductivitity structure reliably. We applied the new inversion to the real data from the Indian sub continent and compared with the results from conventional impedance inversion.
Support vector machines for geophysical inversion
NASA Astrophysics Data System (ADS)
Kuzma, Heidi Anderson
This thesis explores what it means to replace classical non-linear geophysical inversion with computer learning via Support Vector Machines. Geophysical inverse problems are almost always ill-posed which means that many different models (i.e. descriptions of the earth) can be found to explain a given noisy or incomplete data set. Regularization and constraints encourage inversions to find physically realistic models. The set of preferred models needs to be defined a priori using as much geologic knowledge as is available. In inversion, it is assumed that data and a forward modeling process is known. The goal is to solve for a model. In the SVM paradigm, a series of models and associated data are known. The goal is to solve for a reverse modeling process. Starting with a series of initial models assembled using all available geologic information, synthetic data is created using the most realistic forward modeling program available. With the synthetic data as inputs and the known models as outputs, a Support Vector Machine is trained to approximate a local inverse to the forward modeling program. The advantages of this approach are that it is honest about the need to establish, a priori, the kinds of models that are reasonable in a particular field situation. There is no need to adjust the forward process to accommodate inversion, because SVMs can be easily modified to capture complicated, non-linear relationships. SVMs are transparent and require very little programming. If an SVM is trained using model/data pairs that are drawn from the same probability distribution that is implicit in the regularization of an inversion, then it will get very similar results to the inversion. Because SVMs can interpret as much data as desired so long as the conditions of an experiment do not change, they can be used to perform otherwise computationally expensive procedures. The SVMs in this paper are trained to emulate linear and non-linear seismic Amplitude Variation with Offset
Inverse Raman effect: applications and detection techniques
Hughes, L.J. Jr.
1980-08-01
The processes underlying the inverse Raman effect are qualitatively described by comparing it to the more familiar phenomena of conventional and stimulated Raman scattering. An experession is derived for the inverse Raman absorption coefficient, and its relationship to the stimulated Raman gain is obtained. The power requirements of the two fields are examined qualitatively and quantitatively. The assumption that the inverse Raman absorption coefficient is constant over the interaction length is examined. Advantages of the technique are discussed and a brief survey of reported studies is presented.
An inverse problem in thermal imaging
NASA Technical Reports Server (NTRS)
Bryan, Kurt; Caudill, Lester F., Jr.
1994-01-01
This paper examines uniqueness and stability results for an inverse problem in thermal imaging. The goal is to identify an unknown boundary of an object by applying a heat flux and measuring the induced temperature on the boundary of the sample. The problem is studied both in the case in which one has data at every point on the boundary of the region and the case in which only finitely many measurements are available. An inversion procedure is developed and used to study the stability of the inverse problem for various experimental configurations.
Inversion symmetry protected topological insulators and superconductors
NASA Astrophysics Data System (ADS)
Lee, Dung-Hai; Lu, Yuan-Ming
2015-03-01
Three dimensional topological insulator represents a class of novel quantum phases hosting robust gapless boundary excitations, which is protected by global symmetries such as time reversal, charge conservation and spin rotational symmetry. In this work we systematically study another class of topological phases of weakly interacting electrons protected by spatial inversion symmetry, which generally don't support stable gapless boundary states. We classify these inversion-symmetric topological insulators and superconductors in the framework of K-theory, and construct their lattice models. We also discuss quantized response functions of these inversion-protected topological phases, which serve as their experimental signatures.
Population inversion in a stationary recombining plasma
Otsuka, M.
1980-12-01
Population inversion, which occurs in a recombining plasma when a stationary He plasma is brought into contact with a neutral gas, is examined. With hydrogen as a contact gas, noticeable inversion between low-lying levels of H as been found. The overpopulation density is of the order of 10/sup 8/ cm/sup -3/, which is much higher then that (approx. =10/sup 5/ cm/sup -3/) obtained previously with He as a contact gas. Relations between these experimental results and the conditions for population inversion are discussed with the CR model.
Analysis of Temperature Distributions in Nighttime Inversions
NASA Astrophysics Data System (ADS)
Telyak, Oksana; Krasouski, Aliaksandr; Svetashev, Alexander; Turishev, Leonid; Barodka, Siarhei
2015-04-01
Adequate prediction of temperature inversion in the atmospheric boundary layer is one of prerequisites for successful forecasting of meteorological parameters and severe weather events. Examples include surface air temperature and precipitation forecasting as well as prediction of fog, frosts and smog with hazardous levels of atmospheric pollution. At the same time, reliable forecasting of temperature inversions remains an unsolved problem. For prediction of nighttime inversions over some specific territory, it is important to study characteristic features of local circulation cells formation and to properly take local factors into account to develop custom modeling techniques for operational use. The present study aims to investigate and analyze vertical temperature distributions in tropospheric inversions (isotherms) over the territory of Belarus. We study several specific cases of formation, evolution and decay of deep nighttime temperature inversions in Belarus by means of mesoscale numerical simulations with WRF model, considering basic mechanisms of isothermal and inverse temperature layers formation in the troposphere and impact of these layers on local circulation cells. Our primary goal is to assess the feasibility of advance prediction of inversions formation with WRF. Modeling results reveal that all cases under consideration have characteristic features of radiative inversions (e.g., their formation times, development phases, inversion intensities, etc). Regions of "blocking" layers formation are extensive and often spread over the entire territory of Belarus. Inversions decay starts from the lowermost (near surface) layer (altitudes of 5 to 50 m). In all cases, one can observe formation of temperature gradients that substantially differ from the basic inversion gradient, i.e. the layer splits into smaller layers, each having a different temperature stratification (isothermal, adiabatic, etc). As opposed to various empirical techniques as well as
BOOK REVIEW: Inverse Problems. Activities for Undergraduates
NASA Astrophysics Data System (ADS)
Yamamoto, Masahiro
2003-06-01
This book is a valuable introduction to inverse problems. In particular, from the educational point of view, the author addresses the questions of what constitutes an inverse problem and how and why we should study them. Such an approach has been eagerly awaited for a long time. Professor Groetsch, of the University of Cincinnati, is a world-renowned specialist in inverse problems, in particular the theory of regularization. Moreover, he has made a remarkable contribution to educational activities in the field of inverse problems, which was the subject of his previous book (Groetsch C W 1993 Inverse Problems in the Mathematical Sciences (Braunschweig: Vieweg)). For this reason, he is one of the most qualified to write an introductory book on inverse problems. Without question, inverse problems are important, necessary and appear in various aspects. So it is crucial to introduce students to exercises in inverse problems. However, there are not many introductory books which are directly accessible by students in the first two undergraduate years. As a consequence, students often encounter diverse concrete inverse problems before becoming aware of their general principles. The main purpose of this book is to present activities to allow first-year undergraduates to learn inverse theory. To my knowledge, this book is a rare attempt to do this and, in my opinion, a great success. The author emphasizes that it is very important to teach inverse theory in the early years. He writes; `If students consider only the direct problem, they are not looking at the problem from all sides .... The habit of always looking at problems from the direct point of view is intellectually limiting ...' (page 21). The book is very carefully organized so that teachers will be able to use it as a textbook. After an introduction in chapter 1, sucessive chapters deal with inverse problems in precalculus, calculus, differential equations and linear algebra. In order to let one gain some insight
Stochastic inverse consistency in medical image registration.
Yeung, Sai Kit; Shi, Pengcheng
2005-01-01
An essential goal in medical image registration is, the forward and reverse mapping matrices should be inverse to each other, i.e., inverse consistency. Conventional approaches enforce consistency in deterministic fashions, through incorporation of sub-objective cost function to impose source-destination symmetric property during the registration process. Assuming that the initial forward and reverse matching matrices have been computed and used as the inputs to our system, this paper presents a stochastic framework which yields perfect inverse consistency with the simultaneous considerations of the errors underneath the registration matrices and the imperfectness of the consistent constraint. An iterative generalized total least square (GTLS) strategy has been developed such that the inverse consistency is optimally imposed. PMID:16685959
Zinc oxide inverse opal enzymatic biosensor
NASA Astrophysics Data System (ADS)
You, Xueqiu; Pikul, James H.; King, William P.; Pak, James J.
2013-06-01
We report ZnO inverse opal- and nanowire (NW)-based enzymatic glucose biosensors with extended linear detection ranges. The ZnO inverse opal sensors have 0.01-18 mM linear detection range, which is 2.5 times greater than that of ZnO NW sensors and 1.5 times greater than that of other reported ZnO sensors. This larger range is because of reduced glucose diffusivity through the inverse opal geometry. The ZnO inverse opal sensors have an average sensitivity of 22.5 μA/(mM cm2), which diminished by 10% after 35 days, are more stable than ZnO NW sensors whose sensitivity decreased by 10% after 7 days.
Inverse Doppler Effects in Broadband Acoustic Metamaterials.
Zhai, S L; Zhao, X P; Liu, S; Shen, F L; Li, L L; Luo, C R
2016-01-01
The Doppler effect refers to the change in frequency of a wave source as a consequence of the relative motion between the source and an observer. Veselago theoretically predicted that materials with negative refractions can induce inverse Doppler effects. With the development of metamaterials, inverse Doppler effects have been extensively investigated. However, the ideal material parameters prescribed by these metamaterial design approaches are complex and also challenging to obtain experimentally. Here, we demonstrated a method of designing and experimentally characterising arbitrary broadband acoustic metamaterials. These omni-directional, double-negative, acoustic metamaterials are constructed with 'flute-like' acoustic meta-cluster sets with seven double meta-molecules; these metamaterials also overcome the limitations of broadband negative bulk modulus and mass density to provide a region of negative refraction and inverse Doppler effects. It was also shown that inverse Doppler effects can be detected in a flute, which has been popular for thousands of years in Asia and Europe. PMID:27578317
Inverse problem of electro-seismic conversion
NASA Astrophysics Data System (ADS)
Chen, Jie; Yang, Yang
2013-11-01
When a porous rock is saturated with an electrolyte, electrical fields are coupled with seismic waves via the electro-seismic conversion. Pride (1994 Phys. Rev. B 50 15678-96) derived the governing models, in which Maxwell equations are coupled with Biot's equations through the electro-kinetic mobility parameter. The inverse problem of the linearized electro-seismic conversion consists in two steps, namely the inversion of Biot's equations and the inversion of Maxwell equations. We analyze the reconstruction of conductivity and electro-kinetic mobility parameter in Maxwell equations with internal measurements, while the internal measurements are provided by the results of the inversion of Biot's equations. We show that knowledge of two internal data based on well-chosen boundary conditions uniquely determines these two parameters. Moreover, a Lipschitz-type stability is proved based on the same sets of well-chosen boundary conditions.
Electromagnetic inverse applications for functional brain imaging
Wood, C.C.
1997-10-01
This is the final report of a one-year, Laboratory Directed Research and Development (LDRD) project at Los Alamos National Laboratory (LANL). This project addresses an important mathematical and computational problem in functional brain imaging, namely the electromagnetic {open_quotes}inverse problem.{close_quotes} Electromagnetic brain imaging techniques, magnetoencephalography (MEG) and electroencephalography (EEG), are based on measurements of electrical potentials and magnetic fields at hundreds of locations outside the human head. The inverse problem is the estimation of the locations, magnitudes, and time-sources of electrical currents in the brain from surface measurements. This project extends recent progress on the inverse problem by combining the use of anatomical constraints derived from magnetic resonance imaging (MRI) with Bayesian and other novel algorithmic approaches. The results suggest that we can achieve significant improvements in the accuracy and robustness of inverse solutions by these two approaches.
Statistical versus nonstatistical temperature inversion methods
NASA Technical Reports Server (NTRS)
Smith, W. L.; Fleming, H. E.
1972-01-01
Vertical temperature profiles are derived from radiation measurements by inverting the integral equation of radiative transfer. Because of the nonuniqueness of the solution, the particular temperature profile obtained depends on the numerical inversion technique used and the type of auxiliary information incorporated in the solution. The choice of an inversion algorithm depends on many factors; including the speed and size of computer, the availability of representative statistics, and the accuracy of initial data. Results are presented for a numerical study comparing two contrasting inversion methods: the statistical-matrix inversion method and the nonstatistical-iterative method. These were found to be the most applicable to the problem of determining atmospheric temperature profiles. Tradeoffs between the two methods are discussed.
Inverse agonism and its therapeutic significance
Khilnani, Gurudas; Khilnani, Ajeet Kumar
2011-01-01
A large number of G-protein-coupled receptors (GPCRs) show varying degrees of basal or constitutive activity. This constitutive activity is usually minimal in natural receptors but is markedly observed in wild type and mutated (naturally or induced) receptors. According to conventional two-state drug receptor interaction model, binding of a ligand may initiate activity (agonist with varying degrees of positive intrinsic activity) or prevent the effect of an agonist (antagonist with zero intrinsic activity). Inverse agonists bind with the constitutively active receptors, stabilize them, and thus reduce the activity (negative intrinsic activity). Receptors of many classes (α-and β-adrenergic, histaminergic, GABAergic, serotoninergic, opiate, and angiotensin receptors) have shown basal activity in suitable in vitro models. Several drugs that have been conventionally classified as antagonists (β-blockers, antihistaminics) have shown inverse agonist effects on corresponding constitutively active receptors. Nearly all H1 and H2 antihistaminics (antagonists) have been shown to be inverse agonists. Among the β-blockers, carvedilol and bucindolol demonstrate low level of inverse agonism as compared to propranolol and nadolol. Several antipsychotic drugs (D2 receptors antagonist), antihypertensive (AT1 receptor antagonists), antiserotoninergic drugs and opioid antagonists have significant inverse agonistic activity that contributes partly or wholly to their therapeutic value. Inverse agonism may also help explain the underlying mechanism of beneficial effects of carvedilol in congestive failure, naloxone-induced withdrawal syndrome in opioid dependence, clozapine in psychosis, and candesartan in cardiac hypertrophy. Understanding inverse agonisms has paved a way for newer drug development. It is now possible to develop agents, which have only desired therapeutic value and are devoid of unwanted adverse effect. Pimavanserin (ACP-103), a highly selective 5-HT2A inverse
Inversion of elastic impedance for unconsolidated sediments
Lee, Myung W.
2006-01-01
Elastic properties of gas-hydrate-bearing sediments are important for quantifying gas hydrate amounts as well as discriminating the gas hydrate effect on velocity from free gas or pore pressure. This paper presents an elastic inversion method for estimating elastic properties of gas-hydrate-bearing sediments from angle stacks using sequential inversion of P-wave impedance from the zero-offset stack and S-wave impedance from the far-offset stack without assuming velocity ratio.
A fluorophosphate-based inverse Keggin structure
Fielden, John; Quasdorf, Kyle; Cronin, Leroy; Kogerler, Paul
2012-07-17
An unusual PFO(3)(2-)-templated "inverse Keggin" polyanion, [Mo(12)O(46)(PF)(4)](4-), has been isolated from the degradation reaction of an {Mo(132)}-type Keplerate to [PMo(12)O(40)](3-) by [Cu(MeCN)(4)](PF(6)) in acetonitrile. (31)P-NMR studies suggest a structure-directing role for [Cu(MeCN)(4)](+) in the formation of the highly unusual all-inorganic inverse Keggin structure.
Inverse scattering problems with multi-frequencies
NASA Astrophysics Data System (ADS)
Bao, Gang; Li, Peijun; Lin, Junshan; Triki, Faouzi
2015-09-01
This paper is concerned with computational approaches and mathematical analysis for solving inverse scattering problems in the frequency domain. The problems arise in a diverse set of scientific areas with significant industrial, medical, and military applications. In addition to nonlinearity, there are two common difficulties associated with the inverse problems: ill-posedness and limited resolution (diffraction limit). Due to the diffraction limit, for a given frequency, only a low spatial frequency part of the desired parameter can be observed from measurements in the far field. The main idea developed here is that if the reconstruction is restricted to only the observable part, then the inversion will become stable. The challenging task is how to design stable numerical methods for solving these inverse scattering problems inspired by the diffraction limit. Recently, novel recursive linearization based algorithms have been presented in an attempt to answer the above question. These methods require multi-frequency scattering data and proceed via a continuation procedure with respect to the frequency from low to high. The objective of this paper is to give a brief review of these methods, their error estimates, and the related mathematical analysis. More attention is paid to the inverse medium and inverse source problems. Numerical experiments are included to illustrate the effectiveness of these methods.
On the edge of an inverse cascade.
Seshasayanan, Kannabiran; Benavides, Santiago Jose; Alexakis, Alexandros
2014-11-01
We demonstrate that systems with a parameter-controlled inverse cascade can exhibit critical behavior for which at the critical value of the control parameter the inverse cascade stops. In the vicinity of such a critical point, standard phenomenological estimates for the energy balance will fail since the energy flux towards large length scales becomes zero. We demonstrate this using the computationally tractable model of two-dimensional (2D) magnetohydrodynamics in a periodic box. In the absence of any external magnetic forcing, the system reduces to hydrodynamic fluid turbulence with an inverse energy cascade. In the presence of strong magnetic forcing, the system behaves as 2D magnetohydrodynamic turbulence with forward energy cascade. As the amplitude of the magnetic forcing is varied, a critical value is met for which the energy flux towards the large scales becomes zero. Close to this point, the energy flux scales as a power law with the departure from the critical point and the normalized amplitude of the fluctuations diverges. Similar behavior is observed for the flux of the square vector potential for which no inverse flux is observed for weak magnetic forcing, while a finite inverse flux is observed for magnetic forcing above the critical point. We conjecture that this behavior is generic for systems of variable inverse cascade. PMID:25493730
Inverse Bremsstrahlung in Shocked Astrophysical Plasmas
NASA Technical Reports Server (NTRS)
Baring, Matthew G.; Jones, Frank C.; Ellison, Donald C.
2000-01-01
There has recently been interest in the role of inverse bremsstrahlung, the emission of photons by fast suprathermal ions in collisions with ambient electrons possessing relatively low velocities, in tenuous plasmas in various astrophysical contexts. This follows a long hiatus in the application of suprathermal ion bremsstrahlung to astrophysical models since the early 1970s. The potential importance of inverse bremsstrahlung relative to normal bremsstrahlung, i.e. where ions are at rest, hinges upon the underlying velocity distributions of the interacting species. In this paper, we identify the conditions under which the inverse bremsstrahlung emissivity is significant relative to that for normal bremsstrahlung in shocked astrophysical plasmas. We determine that, since both observational and theoretical evidence favors electron temperatures almost comparable to, and certainly not very deficient relative to proton temperatures in shocked plasmas, these environments generally render inverse bremsstrahlung at best a minor contributor to the overall emission. Hence inverse bremsstrahlung can be safely neglected in most models invoking shock acceleration in discrete sources such as supernova remnants. However, on scales approximately > 100 pc distant from these sources, Coulomb collisional losses can deplete the cosmic ray electrons, rendering inverse bremsstrahlung, and perhaps bremsstrahlung from knock-on electrons, possibly detectable.
Iterative inverse kinematics with manipulator configuration control
Grudic, G.Z.; Lawrence, P.D.
1993-08-01
A new method, termed the offset modification method (OM method), for solving the manipulator inverse kinematics problem is presented. The OM method works by modifying the link offset values of a manipulator until it is possible to derive closed-form inverse kinematics equations for the resulting manipulator (termed the model manipulator). This procedure allows one to derive a set of three nonlinear equations in three unknowns that, when numerically solved, give an inverse kinematics solution for the original manipulator. The OM method can be applied to manipulators with any number of degrees of freedom, as long as the manipulator satisfies a given set of conditions (Theorem 1). The OM method is tested on a 6-degree-of-freedom manipulator that has no known closed-form inverse kinematics equations. It is shown that the OM method is applicable to real-time manipulator control, can be used to guarantee convergence to a desired endpoint position and orientation (if it exists), and allows one to directly choose which inverse kinematics solution the algorithm will converge to (as specified in the model manipulator closed-form inverse kinematics equations). Applications of the method to other 6-DOF manipulator geometries and to redundant manipulators (i.e. greater than 6 DOF geometries) are discussed.
Geoacoustic inversion with ships as sources
NASA Astrophysics Data System (ADS)
Koch, Robert A.; Knobles, David P.
2005-02-01
Estimation of geoacoustic parameters using acoustic data from a surface ship was performed for a shallow water region in the Gulf of Mexico. The data were recorded from hydrophones in a bottom mounted, horizontal line array (HLA). The techniques developed to produce the geoacoustic inversion are described, and an efficient method for geoacoustic inversion with broadband beam cross-spectral data is demonstrated. The performance of cost functions that involve coherent or incoherent sums over frequency and one or multiple time segments is discussed. Successful inversions for the first sediment layer sound speed and thickness and some of the parameters for the deeper layers were obtained with the surface ship at nominal ranges of 20, 30, or 50 water depths. The data for these inversions were beam cross-spectra from four subapertures of the HLA spanning a little more than two water depths. The subaperture beams included ten frequencies equally spaced in the 120-200 Hz band. The values of the geoacoustic parameters from the inversions are validated by comparisons with geophysical observations and with the parameter values from previous inversions by other invesigators, and by comparing transmission loss (TL) measured in the experiment with modeled TL based on the inverted geoacoustic parameters. .
Sol-gel co-assembly of hollow cylindrical inverse opals and inverse opal columns.
Haibin, Ni; Ming, Wang; Wei, Chen
2011-12-19
A facile approach of fabricating hollow cylindrical inverse opals and inverse opal columns by sol-gel co-assembly method was proposed. Polystyrene (PS) colloidal suspension added with hydrolyzed silicate precursor solution was used to self-assemble composite colloidal crystals which consist of PS colloidal crystal template and infiltrated silica gel in the interstitial of microspheres. Continuous hollow cylindrical composite colloidal crystal films have been produced on capillaries' outside and internal surface. Composite colloidal crystal columns which filling up the interior of a capillary were fabricated by pressure assisted sol-gel co-assembly method. Hollow cylindrical inverse opals and inverse opal columns were obtained after removing PS colloidal crystal from the composite colloidal crystal. Optical properties of the silica hollow cylindrical inverse opals were characterized by transmission spectrum and a stop band was observed. Structure and optical properties of the inverse opal columns were investigated. PMID:22274178
NASA Technical Reports Server (NTRS)
Sabatier, P. C.
1972-01-01
The progressive realization of the consequences of nonuniqueness imply an evolution of both the methods and the centers of interest in inverse problems. This evolution is schematically described together with the various mathematical methods used. A comparative description is given of inverse methods in scientific research, with examples taken from mathematics, quantum and classical physics, seismology, transport theory, radiative transfer, electromagnetic scattering, electrocardiology, etc. It is hoped that this paper will pave the way for an interdisciplinary study of inverse problems.
Application of the least-squares inversion method: Fourier series versus waveform inversion
NASA Astrophysics Data System (ADS)
Min, Dong-Joo; Shin, Jungkyun; Shin, Changsoo
2015-11-01
We describe an implicit link between waveform inversion and Fourier series based on inversion methods such as gradient, Gauss-Newton, and full Newton methods. Fourier series have been widely used as a basic concept in studies on seismic data interpretation, and their coefficients are obtained in the classical Fourier analysis. We show that Fourier coefficients can also be obtained by inversion algorithms, and compare the method to seismic waveform inversion algorithms. In that case, Fourier coefficients correspond to model parameters (velocities, density or elastic constants), whereas cosine and sine functions correspond to components of the Jacobian matrix, that is, partial derivative wavefields in seismic inversion. In the classical Fourier analysis, optimal coefficients are determined by the sensitivity of a given function to sine and cosine functions. In the inversion method for Fourier series, Fourier coefficients are obtained by measuring the sensitivity of residuals between given functions and test functions (defined as the sum of weighted cosine and sine functions) to cosine and sine functions. The orthogonal property of cosine and sine functions makes the full or approximate Hessian matrix become a diagonal matrix in the inversion for Fourier series. In seismic waveform inversion, the Hessian matrix may or may not be a diagonal matrix, because partial derivative wavefields correlate with each other to some extent, making them semi-orthogonal. At the high-frequency limits, however, the Hessian matrix can be approximated by either a diagonal matrix or a diagonally-dominant matrix. Since we usually deal with relatively low frequencies in seismic waveform inversion, it is not diagonally dominant and thus it is prohibitively expensive to compute the full or approximate Hessian matrix. By interpreting Fourier series with the inversion algorithms, we note that the Fourier series can be computed at an iteration step using any inversion algorithms such as the
Anisotropy effects on 3D waveform inversion
NASA Astrophysics Data System (ADS)
Stekl, I.; Warner, M.; Umpleby, A.
2010-12-01
In the recent years 3D waveform inversion has become achievable procedure for seismic data processing. A number of datasets has been inverted and presented (Warner el al 2008, Ben Hadj at all, Sirgue et all 2010) using isotropic 3D waveform inversion. However the question arises will the results be affected by isotropic assumption. Full-wavefield inversion techniques seek to match field data, wiggle-for-wiggle, to synthetic data generated by a high-resolution model of the sub-surface. In this endeavour, correctly matching the travel times of the principal arrivals is a necessary minimal requirement. In many, perhaps most, long-offset and wide-azimuth datasets, it is necessary to introduce some form of p-wave velocity anisotropy to match the travel times successfully. If this anisotropy is not also incorporated into the wavefield inversion, then results from the inversion will necessarily be compromised. We have incorporated anisotropy into our 3D wavefield tomography codes, characterised as spatially varying transverse isotropy with a tilted axis of symmetry - TTI anisotropy. This enhancement approximately doubles both the run time and the memory requirements of the code. We show that neglect of anisotropy can lead to significant artefacts in the recovered velocity models. We will present inversion results of inverting anisotropic 3D dataset by assuming isotropic earth and compare them with anisotropic inversion result. As a test case Marmousi model extended to 3D with no velocity variation in third direction and with added spatially varying anisotropy is used. Acquisition geometry is assumed as OBC with sources and receivers everywhere at the surface. We attempted inversion using both 2D and full 3D acquisition for this dataset. Results show that if no anisotropy is taken into account although image looks plausible most features are miss positioned in depth and space, even for relatively low anisotropy, which leads to incorrect result. This may lead to
Havla, Lukas; Basha, Tamer; Rayatzadeh, Hussein; Shaw, Jaime L.; Manning, Warren J.; Reeder, Scott B.; Kozerke, Sebastian; Nezafat, Reza
2012-01-01
Purpose To develop an improved chemical shift-based water-fat separation sequence using a water-selective inversion pulse for inversion-recovery 3D contrast-enhanced cardiac MR. Materials and Methods In inversion-recovery sequences, the fat signal is substantially reduced due to the application of a non-selective inversion pulse. Therefore, for simultaneous visualization of water, fat, and myocardial enhancement in inversion-recovery based sequences such as late Gadolinium enhancement imaging, two separate scans are used. To overcome this, the non-selective inversion pulse is replaced with a water-selective inversion pulse. Imaging was performed in phantoms, 9 healthy subjects and 9 patients with suspected arrhythmogenic right ventricular cardiomyopathy plus 1 patient for tumor/mass imaging. In patients, images with conventional turbo-spin echo (TSE) with and without fat saturation were acquired prior to contrast injection for fat assessment. Subjective image scores (1=poor, 4=excellent) were used for image assessment. Results Phantom experiments showed a fat SNR increase between 1.7 to 5.9 times for inversion times of 150 and 300ms, respectively. The water-selective inversion pulse retains the fat signal in contrast-enhanced cardiac MR, allowing improved visualization of fat in the water-fat separated images of healthy subjects with a score of 3.7 ± 0.6. Patient images acquired with the proposed sequence were scored higher when compared with TSE sequence (3.5 ± 0.7 vs. 2.2 ± 0.5, p<0.05). Conclusion The water-selective inversion pulse retains the fat signal in inversion-recovery based contrast-enhanced cardiac MR, allowing simultaneous visualization of water and fat. PMID:22927327
The rotation-inversion spectrum of cyanamide
NASA Astrophysics Data System (ADS)
Read, W. G.; Cohen, E. A.; Pickett, H. M.
1986-02-01
The microwave, millimeter, and submillimeter spectra of cyanamide were studied to better determine the inversion-rotation parameters of the ground and first excited states. A total of 146 transitions including 64 rotation-inversion frequencies between 7 and 500 GHz have been measured at this laboratory. An additional 118 a-type R-branch transitions between 139 and 262 GHz measured by Möller and Winnewisser at Justus Liebig University, Giessen, have also been included in the analysis. The data were fitted to a Hamiltonian which contains a rotation-inversion interaction as an off-diagonal inertial term. The interpretation of the interaction term in terms of the molecular structure and inversion motion is in good agreement with experiment. Higher order rotational effects were handled with a Watson " S" centrifugal distortion Hamiltonian. The inversion splitting, rotational constants, centrifugal distortion constants, nitrogen nuclear quadrupole coupling tensors including χac for the amino nitrogen, and the a and c components of the electric dipole moment are reported.
A semisimultaneous inversion algorithm for SAGE III
NASA Astrophysics Data System (ADS)
Ward, Dale M.
2002-12-01
The Stratospheric Aerosol and Gas Experiment (SAGE) III instrument was successfully launched into orbit on 10 December 2001. The planned operational species separation inversion algorithm will utilize a stepwise retrieval strategy. This paper presents an alternative, semisimultaneous species separation inversion that simultaneously retrieves all species over user-specified vertical intervals or blocks. By overlapping these vertical blocks, retrieved species profiles over the entire vertical range of the measurements are obtained. The semisimultaneous retrieval approach provides a more straightforward method for evaluating the error coupling that occurs among the retrieved profiles due to various types of input uncertainty. Simulation results are presented to show how the semisimultaneous inversion can enhance understanding of the SAGE III retrieval process. In the future, the semisimultaneous inversion algorithm will be used to help evaluate the results and performance of the operational inversion. Compared to SAGE II, SAGE III will provide expanded and more precise spectral measurements. This alone is shown to significantly reduce the uncertainties in the retrieved ozone, nitrogen dioxide, and aerosol extinction profiles for SAGE III. Additionally, the well-documented concern that SAGE II retrievals are biased by the level of volcanic aerosol is greatly alleviated for SAGE III.
Inverse kinematic-based robot control
NASA Technical Reports Server (NTRS)
Wolovich, W. A.; Flueckiger, K. F.
1987-01-01
A fundamental problem which must be resolved in virtually all non-trivial robotic operations is the well-known inverse kinematic question. More specifically, most of the tasks which robots are called upon to perform are specified in Cartesian (x,y,z) space, such as simple tracking along one or more straight line paths or following a specified surfacer with compliant force sensors and/or visual feedback. In all cases, control is actually implemented through coordinated motion of the various links which comprise the manipulator; i.e., in link space. As a consequence, the control computer of every sophisticated anthropomorphic robot must contain provisions for solving the inverse kinematic problem which, in the case of simple, non-redundant position control, involves the determination of the first three link angles, theta sub 1, theta sub 2, and theta sub 3, which produce a desired wrist origin position P sub xw, P sub yw, and P sub zw at the end of link 3 relative to some fixed base frame. Researchers outline a new inverse kinematic solution and demonstrate its potential via some recent computer simulations. They also compare it to current inverse kinematic methods and outline some of the remaining problems which will be addressed in order to render it fully operational. Also discussed are a number of practical consequences of this technique beyond its obvious use in solving the inverse kinematic question.
Improved SOLA Inversions of MDI Data
NASA Astrophysics Data System (ADS)
Larsen, R. M.; Christensen-Dalsgaard, J.; Kosovichev, A. G.; Schou, J.
We present a new version of 2d-SOLA, where the target functions have been modified to match the behavior of the mode kernels near the rotation axis and to minimize near-surface contributions. Inversion of artificial data show that these modifications significantly improve the effective resolution near the pole, which allows us to assess the reliability of the high-latitude features seen by other inversion methods. Most importantly, our new inversions seem to confirm the detection of a submerged polar jet previously seen in the 2d-RLS inversions reported by Schou et al. 1998. A test of the robustness of the improved method is carried out by inverting artificial data from the MDI Hare and Hounds exercise. We analyze the averaging kernels and error propagation of the method, and also describe the error-correlation between different points in the solution, the latter being a potential source of spurious features in the solutions as pointed out by Howe and Thompson, 1996. So far, helioseismic datasets given in the form of a-coefficients have been inverted under the assumption that the errors in different a-coefficients are uncorrelated. The MDI peak-bagging procedure, however, does produce estimates of the error-correlation between a-coefficients within the same multiplet. Here we investigate the effect of including this knowledge in the inversions.
On the nonuniqueness of receiver function inversions
Ammon, C.J. ); Randall, G.E. ); Zandt, G. )
1990-09-10
To study the resolving power of teleseismic P waveforms for receiver structure, the authors model synthetic waveforms using a time domain waveform inversion scheme beginning with a range of initial models to estimate the range of acceptable velocity structures. To speed up the waveform inversions, they implement Randall's (1989) efficient algorithms for calculating differential seismograms and include a smoothness constraint on all the resulting velocity models utilizing the jumping inversion technique of Shaw and Orcutt (1985). They present the results of more than 235 waveform inversions for one-dimensional velocity structures that indicate that the primary sensitivity of a receiver function is to high wavenumber velocity changes, and a depth-velocity product, not simply velocity. The range of slownesses in a typical receiver function study does not appear to be broad enough to remove the depth-velocity ambiguity; the inclusion of a priori information is necessary. They also present inversion results for station RSCP, located in the Cumberland Plateau, Tennessee. The results are similar to those from a previous study by Owens et al. (1984) and demonstrate the uncertainties in the resulting velocity estimate more clearly.
Yao, Jie; Lesage, Anne-Cécile; Hussain, Fazle; Bodmann, Bernhard G.; Kouri, Donald J.
2014-12-15
The reversion of the Born-Neumann series of the Lippmann-Schwinger equation is one of the standard ways to solve the inverse acoustic scattering problem. One limitation of the current inversion methods based on the reversion of the Born-Neumann series is that the velocity potential should have compact support. However, this assumption cannot be satisfied in certain cases, especially in seismic inversion. Based on the idea of distorted wave scattering, we explore an inverse scattering method for velocity potentials without compact support. The strategy is to decompose the actual medium as a known single interface reference medium, which has the same asymptotic form as the actual medium and a perturbative scattering potential with compact support. After introducing the method to calculate the Green’s function for the known reference potential, the inverse scattering series and Volterra inverse scattering series are derived for the perturbative potential. Analytical and numerical examples demonstrate the feasibility and effectiveness of this method. Besides, to ensure stability of the numerical computation, the Lanczos averaging method is employed as a filter to reduce the Gibbs oscillations for the truncated discrete inverse Fourier transform of each order. Our method provides a rigorous mathematical framework for inverse acoustic scattering with a non-compact support velocity potential.
Inversion for seismic anisotropy using genetic algorithms
Horne, S. Univ. of Edinburgh . Dept. of Geology and Geophysics); MacBeth, C. . Dept. of Geology and Geophysics)
1994-11-01
A general inversion scheme based on a genetic algorithm is developed to invert seismic observations for anisotropic parameters. The technique is applied to the inversion of shear-wave observations from two azimuthal VSP data sets from the Conoco test site in Oklahoma. Horizontal polarizations and time-delays are inverted for hexagonal and orthorhombic symmetries. The model solutions are consistent with previous studies using trial and error matching of full waveform synthetics. The shear-wave splitting observations suggest the presence of a shear-wave line singularity and are consistent with a dipping fracture system which is known to exist at the test site. Application of the inversion scheme prior to full waveform modeling demonstrates that a considerable saving in time is possible while retaining the same degree of accuracy.
Optimization and Inverse Design of Pump Impeller
NASA Astrophysics Data System (ADS)
Miyauchi, S.; Zhu, B.; Luo, X.; Piao, B.; Matsumoto, H.; Sano, M.; Kassai, N.
2012-11-01
As for pump impellers, the meridional flow channel and blade-to-blade flow channel, which are relatively independent of each other but greatly affect performance, are designed in parallel. And the optimization design is used for the former and the inverse design is used for the latter. To verify this new design method, a mixed-flow impeller was made. Next, we use Tani's inverse design method for the blade loading of inverse design. It is useful enough to change a deceleration rate freely and greatly. And it can integrally express the rear blade loading of various methods by NACA, Zangeneh and Stratford. We controlled the deceleration rate by shape parameter m, and its value became almost same with Tani's recommended value of the laminar airfoil.
Inverse scattering approach to improving pattern recognition
NASA Astrophysics Data System (ADS)
Chapline, George; Fu, Chi-Yung
2005-05-01
The Helmholtz machine provides what may be the best existing model for how the mammalian brain recognizes patterns. Based on the observation that the "wake-sleep" algorithm for training a Helmholtz machine is similar to the problem of finding the potential for a multi-channel Schrodinger equation, we propose that the construction of a Schrodinger potential using inverse scattering methods can serve as a model for how the mammalian brain learns to extract essential information from sensory data. In particular, inverse scattering theory provides a conceptual framework for imagining how one might use EEG and MEG observations of brain-waves together with sensory feedback to improve human learning and pattern recognition. Longer term, implementation of inverse scattering algorithms on a digital or optical computer could be a step towards mimicking the seamless information fusion of the mammalian brain.
Tabu optimization for matched field inversion
NASA Astrophysics Data System (ADS)
Michalopoulou, Zoi-Heleni; Ghosh-Dastidar, Urmi
2002-11-01
Matched field processing is a powerful tool for source localization and geoacoustic inversion. Because of significant environmental and geometry uncertainties, however, matched field processing usually involves multiparameter searches. To facilitate these searches, global optimization techniques such as genetic algorithms and simulated annealing have been successfully employed. In this work, a different approach, tabu, is implemented for optimization in matched field inversion. Tabu is a technique relying on the use of memory; it searches for the global maximum of the objective function through a navigation process that avoids already revisited models, also making use of aspiration criteria and diversification for faster convergence. The tabu performance in localization and geoacoustic inversion is demonstrated through experimentation with both synthetic and real (SWellEX 96) data. The approach is shown to provide reliable estimates in an efficient manner. [Work supported by ONR.
Fast Inversion of Solar Ca II Spectra
NASA Astrophysics Data System (ADS)
Beck, C.; Choudhary, D. P.; Rezaei, R.; Louis, R. E.
2015-01-01
We present a fast (Lt1 s per profile) inversion code for solar Ca II lines. The code uses an archive of spectra that are synthesized prior to the inversion under the assumption of local thermodynamic equilibrium (LTE). We show that it can be successfully applied to spectrograph data or more sparsely sampled spectra from two-dimensional spectrometers. From a comparison to a non-LTE inversion of the same set of spectra, we derive a first-order non-LTE correction to the temperature stratifications derived in the LTE approach. The correction factor is close to unity up to log τ ~ -3 and increases to values of 2.5 and 4 at log τ = -6 in the quiet Sun and the umbra, respectively.
Inverse Scattering Approach to Improving Pattern Recognition
Chapline, G; Fu, C
2005-02-15
The Helmholtz machine provides what may be the best existing model for how the mammalian brain recognizes patterns. Based on the observation that the ''wake-sleep'' algorithm for training a Helmholtz machine is similar to the problem of finding the potential for a multi-channel Schrodinger equation, we propose that the construction of a Schrodinger potential using inverse scattering methods can serve as a model for how the mammalian brain learns to extract essential information from sensory data. In particular, inverse scattering theory provides a conceptual framework for imagining how one might use EEG and MEG observations of brain-waves together with sensory feedback to improve human learning and pattern recognition. Longer term, implementation of inverse scattering algorithms on a digital or optical computer could be a step towards mimicking the seamless information fusion of the mammalian brain.
Psycholinguistic Evidence for Inverse Scope in Korean.
Lee, Sunyoung; O'Grady, William
2016-08-01
We use experimental data to shed light on the ongoing question of whether Korean allows inverse scope interpretation in sentences containing an indefinite subject and a universally quantified direct object (e.g., 'Someone bought each loaf of bread at the bakery'). The results of an off-line acceptability judgment task (n = 38) and an online self-paced reading task (n [Formula: see text] 22) indicate that inverse scope interpretations are in fact permitted in Korean as a secondary option, as is also the case in English. We argue that the dispreference for the inverse scope reading reflects processing considerations related to burden on working memory. PMID:26022290
Reservoir parameter inversion based on weighted statistics
NASA Astrophysics Data System (ADS)
Gui, Jin-Yong; Gao, Jian-Hu; Yong, Xue-Shan; Li, Sheng-Jun; Liu, Bin-Yang; Zhao, Wan-Jin
2015-12-01
Variation of reservoir physical properties can cause changes in its elastic parameters. However, this is not a simple linear relation. Furthermore, the lack of observations, data overlap, noise interference, and idealized models increases the uncertainties of the inversion result. Thus, we propose an inversion method that is different from traditional statistical rock physics modeling. First, we use deterministic and stochastic rock physics models considering the uncertainties of elastic parameters obtained by prestack seismic inversion and introduce weighting coefficients to establish a weighted statistical relation between reservoir and elastic parameters. Second, based on the weighted statistical relation, we use Markov chain Monte Carlo simulations to generate the random joint distribution space of reservoir and elastic parameters that serves as a sample solution space of an objective function. Finally, we propose a fast solution criterion to maximize the posterior probability density and obtain reservoir parameters. The method has high efficiency and application potential.
Error handling strategies in multiphase inverse modeling
Finsterle, S.; Zhang, Y.
2010-12-01
Parameter estimation by inverse modeling involves the repeated evaluation of a function of residuals. These residuals represent both errors in the model and errors in the data. In practical applications of inverse modeling of multiphase flow and transport, the error structure of the final residuals often significantly deviates from the statistical assumptions that underlie standard maximum likelihood estimation using the least-squares method. Large random or systematic errors are likely to lead to convergence problems, biased parameter estimates, misleading uncertainty measures, or poor predictive capabilities of the calibrated model. The multiphase inverse modeling code iTOUGH2 supports strategies that identify and mitigate the impact of systematic or non-normal error structures. We discuss these approaches and provide an overview of the error handling features implemented in iTOUGH2.
A new strategy for helioseismic inversions
NASA Astrophysics Data System (ADS)
Eff-Darwich, A.; Perez Hernandez, F.
1997-10-01
Helioseismic inversion techniques have been revealed as powerful tools for inferring the internal structure and dynamics of the Sun. One of the most popular techniques is Regularized Least Squares. When it is used, it is necessary to define an inversion mesh and a penalty function, without an a priori knowledge of the behaviour of the solution. In addition, this penalty function is weighted by a trade-off parameter that must be fixed in order to obtain the solution. We present here a new technique, developed in order to find the optimal mesh and smoothing function by means of a deep analysis of the basis functions of the inversion problem. We have found that the method is suitable in particular for obtaining the sound speed and density profiles simultaneously, without any reference to the equation of state.
FAST INVERSION OF SOLAR Ca II SPECTRA
Beck, C.; Choudhary, D. P.; Rezaei, R.; Louis, R. E.
2015-01-10
We present a fast (<<1 s per profile) inversion code for solar Ca II lines. The code uses an archive of spectra that are synthesized prior to the inversion under the assumption of local thermodynamic equilibrium (LTE). We show that it can be successfully applied to spectrograph data or more sparsely sampled spectra from two-dimensional spectrometers. From a comparison to a non-LTE inversion of the same set of spectra, we derive a first-order non-LTE correction to the temperature stratifications derived in the LTE approach. The correction factor is close to unity up to log τ ∼ –3 and increases to values of 2.5 and 4 at log τ = –6 in the quiet Sun and the umbra, respectively.
Acoustic Inversion in Optoacoustic Tomography: A Review
Rosenthal, Amir; Ntziachristos, Vasilis; Razansky, Daniel
2013-01-01
Optoacoustic tomography enables volumetric imaging with optical contrast in biological tissue at depths beyond the optical mean free path by the use of optical excitation and acoustic detection. The hybrid nature of optoacoustic tomography gives rise to two distinct inverse problems: The optical inverse problem, related to the propagation of the excitation light in tissue, and the acoustic inverse problem, which deals with the propagation and detection of the generated acoustic waves. Since the two inverse problems have different physical underpinnings and are governed by different types of equations, they are often treated independently as unrelated problems. From an imaging standpoint, the acoustic inverse problem relates to forming an image from the measured acoustic data, whereas the optical inverse problem relates to quantifying the formed image. This review focuses on the acoustic aspects of optoacoustic tomography, specifically acoustic reconstruction algorithms and imaging-system practicalities. As these two aspects are intimately linked, and no silver bullet exists in the path towards high-performance imaging, we adopt a holistic approach in our review and discuss the many links between the two aspects. Four classes of reconstruction algorithms are reviewed: time-domain (so called back-projection) formulae, frequency-domain formulae, time-reversal algorithms, and model-based algorithms. These algorithms are discussed in the context of the various acoustic detectors and detection surfaces which are commonly used in experimental studies. We further discuss the effects of non-ideal imaging scenarios on the quality of reconstruction and review methods that can mitigate these effects. Namely, we consider the cases of finite detector aperture, limited-view tomography, spatial under-sampling of the acoustic signals, and acoustic heterogeneities and losses. PMID:24772060
Geometry and kinematic evolution of inversion structures
Mitra, S. )
1993-07-01
Positive inversion structures form by the compressional reactivation of preexisting extensional structures. Experimental models and observations of natural structures are used to develop quantitative models for the geometry and kinematic evolution of inversion structures. In this paper, I analyze two main formation mechanisms of inversion structures: (1) fault-propagation folding on planar faults, and (2) fault-bend folding on listric faults. Inversion structures formed by fault-propagation folding occur in the southern North Sea, the Central Montana platform, and the Kangean Basin. During extension, a broad fault-propagation (or drape) fold develops above the master fault, with the fault subsequently breaking through the fold. Synextensional growth units deposited in the hanging wall typically thicken into the basin. Compressional reactivation results in slip reversal on the master and secondary faults, their rotation to shallower dips, and the development of a compressional fault-propagation fold. Inversion structures formed by fault-bend folding on listric faults occur in the Taranaki Basin, and possibly in the southern North Sea. Rollover folding in the hanging wall occurs during extension, possibly accompanied by a small component of fault-propagation folding in the vicinity of the fault tip. Deformation is primarily along a system of antithetic faults. Synextensional growth sediments typically thicken into the fault, but also show thinning in the immediate vicinity of the fault. During compression, the extensional fold is first unfolded and then folded into a compressional fault-bend fold. The characteristic variations in bed geometry and thickness provide predictive models for interpreting the subsurface geometries of these two classes of inversion structures in areas with poor seismic data. These models are particularly useful in exploring for structural traps in the complex and relatively unexplored synextensional growth units. 31 refs., 29 figs.
Trimming and procrastination as inversion techniques
NASA Astrophysics Data System (ADS)
Backus, George E.
1996-12-01
By examining the processes of truncating and approximating the model space (trimming it), and by committing to neither the objectivist nor the subjectivist interpretation of probability (procrastinating), we construct a formal scheme for solving linear and non-linear geophysical inverse problems. The necessary prior information about the correct model xE can be either a collection of inequalities or a probability measure describing where xE was likely to be in the model space X before the data vector y0 was measured. The results of the inversion are (1) a vector z0 that estimates some numerical properties zE of xE; (2) an estimate of the error δz = z0 - zE. As y0 is finite dimensional, so is z0, and hence in principle inversion cannot describe all of xE. The error δz is studied under successively more specialized assumptions about the inverse problem, culminating in a complete analysis of the linear inverse problem with a prior quadratic bound on xE. Our formalism appears to encompass and provide error estimates for many of the inversion schemes current in geomagnetism, and would be equally applicable in geodesy and seismology if adequate prior information were available there. As an idealized example we study the magnetic field at the core-mantle boundary, using satellite measurements of field elements at sites assumed to be almost uniformly distributed on a single spherical surface. Magnetospheric currents are neglected and the crustal field is idealized as a random process with rotationally invariant statistics. We find that an appropriate data compression diagonalizes the variance matrix of the crustal signal and permits an analytic trimming of the idealized problem.
Computational 3-D inversion for seismic exploration
Gavrilov, E.M.; Forslund, D.W.; Fehler, M.C.
1997-10-01
This is the final report of a four-month, Laboratory Directed Research and Development (LDRD) project carried out at the Los Alamos National Laboratory (LANL). There is a great need for a new and effective technology with a wide scope of industrial applications to investigate media internal properties of which can be explored only from the backscattered data. The project was dedicated to the development of a three-dimensional computational inversion tool for seismic exploration. The new computational concept of the inversion algorithm was suggested. The goal of the project was to prove the concept and the practical validity of the algorithm for petroleum exploration.
Direct and Inverse problems in Electrocardiography
NASA Astrophysics Data System (ADS)
Boulakia, M.; Fernández, M. A.; Gerbeau, J. F.; Zemzemi, N.
2008-09-01
We present numerical results related to the direct and the inverse problems in electrocardiography. The electrical activity of the heart is described by the bidomain equations. The electrocardiograms (ECGs) recorded in different points on the body surface are obtained by coupling the bidomain equation to a Laplace equation in the torso. The simulated ECGs are quite satisfactory. As regards the inverse problem, our goal is to estimate the parameters of the bidomain-torso model. Here we present some preliminary results of a parameter estimation for the torso model.
Bayesian Inference in Satellite Gravity Inversion
NASA Technical Reports Server (NTRS)
Kis, K. I.; Taylor, Patrick T.; Wittmann, G.; Kim, Hyung Rae; Torony, B.; Mayer-Guerr, T.
2005-01-01
To solve a geophysical inverse problem means applying measurements to determine the parameters of the selected model. The inverse problem is formulated as the Bayesian inference. The Gaussian probability density functions are applied in the Bayes's equation. The CHAMP satellite gravity data are determined at the altitude of 400 kilometer altitude over the South part of the Pannonian basin. The model of interpretation is the right vertical cylinder. The parameters of the model are obtained from the minimum problem solved by the Simplex method.
Dispersion analysis with inverse dielectric function modelling.
Mayerhöfer, Thomas G; Ivanovski, Vladimir; Popp, Jürgen
2016-11-01
We investigate how dispersion analysis can profit from the use of a Lorentz-type description of the inverse dielectric function. In particular at higher angles of incidence, reflectance spectra using p-polarized light are dominated by bands from modes that have their transition moments perpendicular to the surface. Accordingly, the spectra increasingly resemble inverse dielectric functions. A corresponding description can therefore eliminate the complex dependencies of the dispersion parameters, allow their determination and facilitate a more accurate description of the optical properties of single crystals. PMID:27294550
Directional wetting in anisotropic inverse opals.
Phillips, Katherine R; Vogel, Nicolas; Burgess, Ian B; Perry, Carole C; Aizenberg, Joanna
2014-07-01
Porous materials display interesting transport phenomena due to restricted motion of fluids within the nano- to microscale voids. Here, we investigate how liquid wetting in highly ordered inverse opals is affected by anisotropy in pore geometry. We compare samples with different degrees of pore asphericity and find different wetting patterns depending on the pore shape. Highly anisotropic structures are infiltrated more easily than their isotropic counterparts. Further, the wetting of anisotropic inverse opals is directional, with liquids filling from the side more easily. This effect is supported by percolation simulations as well as direct observations of wetting using time-resolved optical microscopy. PMID:24941308
Inverse melting in stressed fused silica
NASA Astrophysics Data System (ADS)
Bouchut, Philippe
2012-12-01
The emissive properties of proton implanted fused silica surfaces have been studied by laser beam annealing. When submitted to a high thermal step from a focused CO2 laser, an intense near infra-red thermoluminescence peak rises at a heating rate threshold. The in plane tensile stress relaxes and silica melts. We show that in the irreversible inverse melting of stressed fused silica, the protons exo-diffuse through internal modes coupling. The heat and mass transfer is one entropy flux whose dynamics are regulated by the mass transport. Inverse melting is the thermodynamic process that initiates the glass transition when heating.
3D Electromagnetic inversion using conjugate gradients
Newman, G.A.; Alumbaugh, D.L.
1997-06-01
In large scale 3D EM inverse problems it may not be possible to directly invert a full least-squares system matrix involving model sensitivity elements. Thus iterative methods must be employed. For the inverse problem, we favor either a linear or non-linear (NL) CG scheme, depending on the application. In a NL CG scheme, the gradient of the objective function is required at each relaxation step along with a univariate line search needed to determine the optimum model update. Solution examples based on both approaches will be presented.
Aerosol physical properties from satellite horizon inversion
NASA Technical Reports Server (NTRS)
Gray, C. R.; Malchow, H. L.; Merritt, D. C.; Var, R. E.; Whitney, C. K.
1973-01-01
The feasibility is investigated of determining the physical properties of aerosols globally in the altitude region of 10 to 100 km from a satellite horizon scanning experiment. The investigation utilizes a horizon inversion technique previously developed and extended. Aerosol physical properties such as number density, size distribution, and the real and imaginary components of the index of refraction are demonstrated to be invertible in the aerosol size ranges (0.01-0.1 microns), (0.1-1.0 microns), (1.0-10 microns). Extensions of previously developed radiative transfer models and recursive inversion algorithms are displayed.
Antarctic Crustal Thickness from Gravity Inversion
NASA Astrophysics Data System (ADS)
Vaughan, A. P.; Kusznir, N. J.; Ferraccioli, F.; Jordan, T. A.
2013-12-01
Using gravity anomaly inversion, we have produced the first comprehensive regional maps of crustal thickness and oceanic lithosphere distribution for Antarctica and the Southern Ocean. We determine Moho depth, crustal basement thickness, continental lithosphere thinning (1-1/β) and ocean-continent transition location using a 3D spectral domain gravity inversion method, which incorporates a lithosphere thermal gravity anomaly correction. The continental lithosphere thinning distribution, used to define the initial thermal model temperature perturbation is derived from the gravity inversion and uses no a priori isochron information; as a consequence the gravity inversion method provides a prediction of ocean-continent transition location, which is independent of ocean isochron information. The gravity anomaly contribution from ice thickness is included in the gravity inversion, as is the contribution from sediments which assumes a compaction controlled sediment density increase with depth. Data used in the gravity inversion are elevation and bathymetry, free-air gravity anomaly, the most recent Bedmap2 ice thickness and bedrock topography compilation south of 60 degrees south (Fretwell et al., 2013) and relatively sparse constraints on sediment thickness. Our gravity inversion study predicts thick crust (> 45 km) under interior East Antarctica penetrated by narrow continental rifts that feature relatively thinner crust. The East Antarctic Rift System (EARS) is a major Permian to Cretaceous age rift system that appears to extend from the continental margin at the Lambert Rift to the South Pole region, a distance of 2500 km. This is comparable in scale to the well-studied East African rift system. Intermediate crustal thickness with an inferred linear rift fabric is predicted under Coates Land. An extensive region of either thick oceanic crust or highly thinned continental crust is predicted offshore Oates Land and north Victoria Land, and also off West Antarctica
Uncertainty estimations for seismic source inversions
NASA Astrophysics Data System (ADS)
Duputel, Z.; Rivera, L. A.; Fukahata, Y.; Kanamori, H.
2011-12-01
Source inversion is a very widely used practice in seismology. Magnitudes, moment tensors, slip distributions are now routinely calculated and disseminated by several agencies and research groups whenever an earthquake occurs. The estimated source models can be used as inputs for various algorithms such as ShakeMap computation, tsunami modeling, stress transfer calculation or waveform modeling for tomography studies. Despite the importance of these applications, the source inversion algorithms often do not include proper error analyses, and the results are often given without any estimates of uncertainties. In centroid moment tensor (CMT) inversion studies, we often estimate the uncertainty on the model parameters by using various resampling techniques such as bootstrap or jacknife. The strength of these computer-based methods lies in their simplicity. We can implement them considering the inversion procedure as a "black-box" without any knowledge about the model and data statistical properties. However, these methods can suffer from too simplistic assumptions (such as the independence of data samples) and provide the first order error estimates only without the possibility of improving the source model itself. We explore here an alternative approach by taking errors explicitly into account in source inversion problems. In this perspective we use the W-phase source inversion algorithm recently developed to provide fast and robust CMT estimations for moderate to large earthquakes. We assume that the initial probability densities can be modeled by Gaussian distributions. Formally, we can separate two sources of error which generally contribute to the model parameter uncertainties. On one side we consider the error induced by the more or less imperfect data. This information is carried by the covariance matrix for the data Cd. A key point which is practically always ignored is the possibility of having non-diagonal elements in Cd; such non-diagonal elements are due to
Effects of Tape and Exercise on Dynamic Ankle Inversion
Ricard, Mark D.; Sherwood, Stephen M.; Schulthies, Shane S.; Knight, Kenneth L.
2000-01-01
Objective: To compare the effects of tape, with and without prewrap, on dynamic ankle inversion before and after exercise. Design and Setting: Doubly multivariate analyses of variance were used to compare the taping and exercise conditions. Subjects were randomly assigned to a fixed treatment order as determined by a balanced latin square. The independent variables were tape application (no tape, tape with prewrap, tape to skin) and exercise (before and after). The dependent variables were average inversion velocity, total inversion, maximum inversion velocity, and time to maximum inversion. Subjects: Thirty college-age male and female students (17 males, 13 females; mean age = 24.9 ± 4.3 years, range, 19 to 39 years) were tested. Subjects were excluded from the study if they exhibited a painful gait or painful range of motion or had a past history of ankle surgery or an ankle sprain within the past 4 weeks. Measurements: We collected data using electronic goniometers while subjects balanced on the right leg on an inversion platform tilted about the medial-lateral axis to produce 15° of plantar flexion. Sudden ankle inversion was induced by pulling the inversion platform support, allowing the platform support base to rotate 37°. Ten satisfactory trials were recorded on the inversion platform before and after a prescribed exercise bout. We calculated total inversion, time to maximum inversion, average inversion velocity, and maximum inversion velocity after sudden inversion. Results: We found no significant differences between taping to the skin and taping over prewrap for any of the variables measured. There were significant differences between both taping conditions and no-tape postexercise for average inversion velocity, maximum inversion, maximum inversion velocity, and time to maximum inversion. The total inversion mean for no-tape postexercise was 38.8° ± 6.3°, whereas the means for tape and skin and for tape and prewrap were 28.3° ± 4.6° and 29.1°
Inversion Symmetry Breaking in Endohedral C_60
NASA Astrophysics Data System (ADS)
Clougherty, Dennis; Anderson, Frederick
1998-03-01
A pseudo--Jahn--Teller model describing central atom distortions is proposed for endohedral fullerenes of the form A@C_60 where A is either a rare gas or a metal atom. A critical (dimensionless) coupling gc is found, at or below which the symmetric configuration is stable and above which inversion symmetry is broken. Vibronic parameters are given for selected endohedral fullerenes.
Inverse-Square Orbits: A Geometric Approach.
ERIC Educational Resources Information Center
Rainwater, James C.; Weinstock, Robert
1979-01-01
Presents a derivation of Kepler's first law of planetary motion from Newtonian principles. Analogus derivations of the hyperbolic and parabolic orbits of nonreturning comets and the hyperbolic orbit for a particle in a repulsive inverse-square field are also presented. (HM)
Towards Lasing Without Inversion in Neutral Mercury
NASA Astrophysics Data System (ADS)
Rein, Benjamin; Sturm, Martin R.; Walser, Reinhold; Walther, Thomas
2015-03-01
Currently, we are implementing a lasing without inversion (LWI) scheme in mercury based on calculations by Fry et al. [1]. A recent detailed analysis predicting exact experimental parameters shows the feasibility of our LWI scheme [2]. In this paper we report on the progress of our experiment and compare its parameters with the prior analysis.
Inverse Doppler Effects in Broadband Acoustic Metamaterials
Zhai, S. L.; Zhao, X. P.; Liu, S.; Shen, F. L.; Li, L. L.; Luo, C. R.
2016-01-01
The Doppler effect refers to the change in frequency of a wave source as a consequence of the relative motion between the source and an observer. Veselago theoretically predicted that materials with negative refractions can induce inverse Doppler effects. With the development of metamaterials, inverse Doppler effects have been extensively investigated. However, the ideal material parameters prescribed by these metamaterial design approaches are complex and also challenging to obtain experimentally. Here, we demonstrated a method of designing and experimentally characterising arbitrary broadband acoustic metamaterials. These omni-directional, double-negative, acoustic metamaterials are constructed with ‘flute-like’ acoustic meta-cluster sets with seven double meta-molecules; these metamaterials also overcome the limitations of broadband negative bulk modulus and mass density to provide a region of negative refraction and inverse Doppler effects. It was also shown that inverse Doppler effects can be detected in a flute, which has been popular for thousands of years in Asia and Europe. PMID:27578317
NASA Inverse Methods/Data Assimilation
NASA Technical Reports Server (NTRS)
Bennett, Andrew
2003-01-01
An overview of NASA's Third International Summer School on Inverse Methods and Data Assimilation which was conducted at Oregon State University from July 22 to August 2, 2002, is presented. Items listed include: a roster of attendees, a description of course content and talks given.
Inverse halftoning with context driven prediction.
Jing-Ming Guo; Yun-Fu Liu; Jen-Ho Chen; Jiann-Der Lee
2014-04-01
A prior work proposed by Chung-Wu considered an edge-based lookup table to obtain good inversed image quality, yet it suffers from some drawbacks in terms of image quality, memory consumption, and complexity. In this correspondence, an improved scheme is proposed to deal with these issues. PMID:24808357
A Face Inversion Effect without a Face
ERIC Educational Resources Information Center
Brandman, Talia; Yovel, Galit
2012-01-01
Numerous studies have attributed the face inversion effect (FIE) to configural processing of internal facial features in upright but not inverted faces. Recent findings suggest that face mechanisms can be activated by faceless stimuli presented in the context of a body. Here we asked whether faceless stimuli with or without body context may induce…
Sparse matrix orderings for factorized inverse preconditioners
Benzi, M.; Tuama, M.
1998-09-01
The effect of reorderings on the performance of factorized sparse approximate inverse preconditioners is considered. It is shown that certain reorderings can be very beneficial both in the preconditioner construction phase and in terms of the rate of convergence of the preconditioned iteration.
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.
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…
An Inverse Model for TETRAD: Preliminary Results
Shook, George Michael; Renner, Joel Lawrence
2002-09-01
A model-independent parameter estimation model known as PEST has been linked to the reservoir simulator TETRAD. The method of inverse modeling is briefly reviewed, and the link between PEST and TETRAD is discussed. A single example is presented that illustrates the power of parameter estimation from well observations.
An approximation for inverse Laplace transforms
NASA Technical Reports Server (NTRS)
Lear, W. M.
1981-01-01
Programmable calculator runs simple finite-series approximation for Laplace transform inversions. Utilizing family of orthonormal functions, approximation is used for wide range of transforms, including those encountered in feedback control problems. Method works well as long as F(t) decays to zero as it approaches infinity and so is appliable to most physical systems.
Applications of high resolution inverse Raman spectroscopy
Owyoung, A.; Esherick, P.
1980-01-01
The use of high-power, narrow-band lasers has significantly improved the resolving power and sensitivity of inverse Raman spectroscopy of gases. In this paper we shall describe this technique, illustrate its capabilities by showing some Q-branch spectra of heavy spherical tops, and survey some possible future applications.
Uncertainty estimations for seismic source inversions
NASA Astrophysics Data System (ADS)
Duputel, Zacharie; Rivera, Luis; Fukahata, Yukitoshi; Kanamori, Hiroo
2012-08-01
Source inversion is a widely used practice in seismology. Magnitudes, moment tensors, slip distributions are now routinely calculated and disseminated whenever an earthquake occurs. The accuracy of such models depends on many aspects like the event magnitude, the data coverage and the data quality (instrument response, isolation, timing, etc.). Here, like in any observational problem, the error estimation should be part of the solution. It is however very rare to find a source inversion algorithm which includes realistic error analyses, and the solutions are often given without any estimates of uncertainties. Our goal here is to stress the importance of such estimation and to explore different techniques aimed at achieving such analyses. In this perspective, we use the W phase source inversion algorithm recently developed to provide fast CMT estimations for large earthquakes. We focus in particular on the linear-inverse problem of estimating the moment tensor components at a given source location. We assume that the initial probability densities can be modelled by Gaussian distributions. Formally, we can separate two sources of error which generally contribute to the model parameter uncertainties. The first source of uncertainty is the error introduced by the more or less imperfect data. This is carried by the covariance matrix for the data (Cd). The second source of uncertainty, often overlooked, is associated with modelling error or mismodelling. This is represented by the covariance matrix on the theory, CT. Among the different sources of mismodelling, we focus here on the modelling error associated with the mislocation of the centroid position. Both Cd and CT describe probability densities in the data space and it is well known that it is in fact CD=Cd+CT that should be included into the error propagation process. In source inversion problems, like in many other fields of geophysics, the data covariance (CD) is often considered as diagonal or even proportional
Optical properties of silicon inverse opals
NASA Astrophysics Data System (ADS)
Wei, Hong
Silicon inverse opals are artificial structures in which nearly monodisperse, close-packed air bubbles are embedded in a silicon matrix. If properly tailored, this structure can exhibit a photonic band gap (PBG) in the near infrared spectral region. The PBG can block light propagation in any direction, allowing the control of light flow in the material. Silicon inverse opals can be fabricated by infiltrating amorphous silicon into silica colloidal crystals and then etching away the silica. In this thesis, the structural defects of silica colloidal crystals and the optical properties of silicon inverse opals are studied. First, by using laser-scanning confocal microscopy, the concentration and distribution of stacking faults and vacancies were quantified in silica colloidal crystals. It's shown that silica colloidal crystals show strong tendency toward face-center-cubic structure with the vacancy density as small as 5 x 10-4. Second, by combining optical microscopy and Fourier Transform Infrared (FTIR) spectroscopy, the transmission and reflection spectra of silicon inverse opals along the [111] direction were measured. Combined with the calculation of transmission and reflection spectra by Transfer Matrix Methods, it is concluded that the strong light attenuation in silicon inverse opals is due to the enhanced absorption (>600%) in silicon materials. Third, by using optical pump-probe techniques, the photo-induced ultra-fast reflection changes in silicon inverse opals were examined. The pump-generated free carriers cause the reflection in the band gap region to change after ˜0.5 ps. For the first few ps, the main effect is a decrease in reflectivity due to nonlinear absorption. After ˜5 ps, this effect disappears and an unexpected blue spectral shift is seen in the photonic band gap. The refractive index decreases due to optically-induced strain born the thermal expansion mismatch between silicon and its native oxide. Finally, by infiltrating silicon inverse
Clinical knowledge-based inverse treatment planning
NASA Astrophysics Data System (ADS)
Yang, Yong; Xing, Lei
2004-11-01
Clinical IMRT treatment plans are currently made using dose-based optimization algorithms, which do not consider the nonlinear dose-volume effects for tumours and normal structures. The choice of structure specific importance factors represents an additional degree of freedom of the system and makes rigorous optimization intractable. The purpose of this work is to circumvent the two problems by developing a biologically more sensible yet clinically practical inverse planning framework. To implement this, the dose-volume status of a structure was characterized by using the effective volume in the voxel domain. A new objective function was constructed with the incorporation of the volumetric information of the system so that the figure of merit of a given IMRT plan depends not only on the dose deviation from the desired distribution but also the dose-volume status of the involved organs. The conventional importance factor of an organ was written into a product of two components: (i) a generic importance that parametrizes the relative importance of the organs in the ideal situation when the goals for all the organs are met; (ii) a dose-dependent factor that quantifies our level of clinical/dosimetric satisfaction for a given plan. The generic importance can be determined a priori, and in most circumstances, does not need adjustment, whereas the second one, which is responsible for the intractable behaviour of the trade-off seen in conventional inverse planning, was determined automatically. An inverse planning module based on the proposed formalism was implemented and applied to a prostate case and a head-neck case. A comparison with the conventional inverse planning technique indicated that, for the same target dose coverage, the critical structure sparing was substantially improved for both cases. The incorporation of clinical knowledge allows us to obtain better IMRT plans and makes it possible to auto-select the importance factors, greatly facilitating the inverse
Structural information in the inverse problem
NASA Astrophysics Data System (ADS)
Karaoulis, M.; Larson, T. H.; Ahmed, I.; Revil, A.; Thomason, J.
2013-12-01
Data integration in geophysics provides additional information to elucidate subsurface structure and reduce non-uniqueness of inverted models. There are several strategies for incorporating data integration into numerical models. A traditional approach is to use this information as a-priori knowledge, as an initial model or layer boundary specified on the mesh before the inversion. Although in some cases this approach has proven effective, the information is not directly incorporated into the inverse problem and might be lost in the final model. Another strategy is through joint inversion, where data and models are inverted simultaneous. The data integration comes from the joint operator as structural similarity or petrophysical equations. There are some limitations with this approach. In particular, structural similarity doesn't take into account different sensitivity patterns which differ for different geophysical methods, e.g. in a crosswell configuration electrical resistivity tomography is sensitive close to the borehole region while seismic waves are sensitive towards the center part. Moreover different methods have differing resolution. Therefore, a single joint operator might not be effective in all cases. In this work we demonstrate the use of image-guided inversion, where structural information is taken directly from a high resolution geophysical image (e.g. ground penetrating radar or seismic reflection) or from a geological cross-section. This structural information is introduced into the inverse problem through a weighted smoothing matrix, where it correlates and favors formations related to a specific structural feature and not just uniformly across the entire model. Both sharp and smooth features can be imaged and the recovered models can have a more realistic distribution of values. As an example of the method we use migrated seismic reflection images to extract the structural information and resistivity imaging to recover the resistivity
Estimating uncertainties in complex joint inverse problems
NASA Astrophysics Data System (ADS)
Afonso, Juan Carlos
2016-04-01
Sources of uncertainty affecting geophysical inversions can be classified either as reflective (i.e. the practitioner is aware of her/his ignorance) or non-reflective (i.e. the practitioner does not know that she/he does not know!). Although we should be always conscious of the latter, the former are the ones that, in principle, can be estimated either empirically (by making measurements or collecting data) or subjectively (based on the experience of the researchers). For complex parameter estimation problems in geophysics, subjective estimation of uncertainty is the most common type. In this context, probabilistic (aka Bayesian) methods are commonly claimed to offer a natural and realistic platform from which to estimate model uncertainties. This is because in the Bayesian approach, errors (whatever their nature) can be naturally included as part of the global statistical model, the solution of which represents the actual solution to the inverse problem. However, although we agree that probabilistic inversion methods are the most powerful tool for uncertainty estimation, the common claim that they produce "realistic" or "representative" uncertainties is not always justified. Typically, ALL UNCERTAINTY ESTIMATES ARE MODEL DEPENDENT, and therefore, besides a thorough characterization of experimental uncertainties, particular care must be paid to the uncertainty arising from model errors and input uncertainties. We recall here two quotes by G. Box and M. Gunzburger, respectively, of special significance for inversion practitioners and for this session: "…all models are wrong, but some are useful" and "computational results are believed by no one, except the person who wrote the code". In this presentation I will discuss and present examples of some problems associated with the estimation and quantification of uncertainties in complex multi-observable probabilistic inversions, and how to address them. Although the emphasis will be on sources of uncertainty related
PUBLISHER'S ANNOUNCEMENT: New developments for Inverse Problems
NASA Astrophysics Data System (ADS)
2006-12-01
2006 has proved to be a very successful year for Inverse Problems. After an increase for the fourth successive year, we achieved our highest impact factor to date, 1.541 (Source: 2005 ISI® Journal Citation Report), and the Editorial Board is keen to build on this success by continuing to improve the service we offer to our readers and authors. The Board has observed that Inverse Problems receives very few Letters to the Editor submissions, and that moreover those that we do receive rarely conform to the requirements for Letters to the Editor set out in the journal's editorial policy. The Board has therefore decided to merge the current Letters to the Editor section into our regular Papers section, which will now accommodate all research articles that meet the journal's high quality standards. Any submissions that would previously have been Letters to the Editor are still very welcome as Papers, and can be submitted by e-mail to ip@iop.org or online using our online submissions form at authors.iop.org/submit. Inverse Problems' processing times are already among the fastest in the field—on average, authors receive our decision on their paper in less than three months. Thanks to our easy-to-use online refereeing system, publishing a Paper is now just as fast as publishing a Letter to the Editor, and we are striving to ensure that the journal's high standards are applied consistently to all our Papers, maintaining Inverse Problems' position as the leading journal in the field. Our highly acclaimed Topical Review section will also continue and grow; providing timely insights into the development of all topical fields within Inverse Problems. We have many exciting Topical Reviews currently in preparation for 2007 and will continue to commission articles at the cutting edge of research. We look forward to receiving your contributions and to continuing to provide the best publication service available.
SISYPHUS: A high performance seismic inversion factory
NASA Astrophysics Data System (ADS)
Gokhberg, Alexey; Simutė, Saulė; Boehm, Christian; Fichtner, Andreas
2016-04-01
In the recent years the massively parallel high performance computers became the standard instruments for solving the forward and inverse problems in seismology. The respective software packages dedicated to forward and inverse waveform modelling specially designed for such computers (SPECFEM3D, SES3D) became mature and widely available. These packages achieve significant computational performance and provide researchers with an opportunity to solve problems of bigger size at higher resolution within a shorter time. However, a typical seismic inversion process contains various activities that are beyond the common solver functionality. They include management of information on seismic events and stations, 3D models, observed and synthetic seismograms, pre-processing of the observed signals, computation of misfits and adjoint sources, minimization of misfits, and process workflow management. These activities are time consuming, seldom sufficiently automated, and therefore represent a bottleneck that can substantially offset performance benefits provided by even the most powerful modern supercomputers. Furthermore, a typical system architecture of modern supercomputing platforms is oriented towards the maximum computational performance and provides limited standard facilities for automation of the supporting activities. We present a prototype solution that automates all aspects of the seismic inversion process and is tuned for the modern massively parallel high performance computing systems. We address several major aspects of the solution architecture, which include (1) design of an inversion state database for tracing all relevant aspects of the entire solution process, (2) design of an extensible workflow management framework, (3) integration with wave propagation solvers, (4) integration with optimization packages, (5) computation of misfits and adjoint sources, and (6) process monitoring. The inversion state database represents a hierarchical structure with
GENERATING FRACTAL PATTERNS BY USING p-CIRCLE INVERSION
NASA Astrophysics Data System (ADS)
Ramírez, José L.; Rubiano, Gustavo N.; Zlobec, Borut Jurčič
2015-10-01
In this paper, we introduce the p-circle inversion which generalizes the classical inversion with respect to a circle (p = 2) and the taxicab inversion (p = 1). We study some basic properties and we also show the inversive images of some basic curves. We apply this new transformation to well-known fractals such as Sierpinski triangle, Koch curve, dragon curve, Fibonacci fractal, among others. Then we obtain new fractal patterns. Moreover, we generalize the method called circle inversion fractal be means of the p-circle inversion.
A Bayesian method for microseismic source inversion
NASA Astrophysics Data System (ADS)
Pugh, D. J.; White, R. S.; Christie, P. A. F.
2016-05-01
Earthquake source inversion is highly dependent on location determination and velocity models. Uncertainties in both the model parameters and the observations need to be rigorously incorporated into an inversion approach. Here, we show a probabilistic Bayesian method that allows formal inclusion of the uncertainties in the moment tensor inversion. This method allows the combination of different sets of far-field observations, such as P-wave and S-wave polarities and amplitude ratios, into one inversion. Additional observations can be included by deriving a suitable likelihood function from the uncertainties. This inversion produces samples from the source posterior probability distribution, including a best-fitting solution for the source mechanism and associated probability. The inversion can be constrained to the double-couple space or allowed to explore the gamut of moment tensor solutions, allowing volumetric and other non-double-couple components. The posterior probability of the double-couple and full moment tensor source models can be evaluated from the Bayesian evidence, using samples from the likelihood distributions for the two source models, producing an estimate of whether or not a source is double-couple. Such an approach is ideally suited to microseismic studies where there are many sources of uncertainty and it is often difficult to produce reliability estimates of the source mechanism, although this can be true of many other cases. Using full-waveform synthetic seismograms, we also show the effects of noise, location, network distribution and velocity model uncertainty on the source probability density function. The noise has the largest effect on the results, especially as it can affect other parts of the event processing. This uncertainty can lead to erroneous non-double-couple source probability distributions, even when no other uncertainties exist. Although including amplitude ratios can improve the constraint on the source probability
Uncertainty quantification for ice sheet inverse problems
NASA Astrophysics Data System (ADS)
Petra, N.; Ghattas, O.; Stadler, G.; Zhu, H.
2011-12-01
Modeling the dynamics of polar ice sheets is critical for projections of future sea level rise. Yet, there remain large uncertainties in the basal boundary conditions and in the non-Newtonian constitutive relations employed within ice sheet models. In this presentation, we consider the problem of estimating uncertainty in the solution of (large-scale) ice sheet inverse problems within the framework of Bayesian inference. Computing the general solution of the inverse problem-i.e., the posterior probability density-is intractable with current methods on today's computers, due to the expense of solving the forward model (3D full Stokes flow with nonlinear rheology) and the high dimensionality of the uncertain parameters (which are discretizations of the basal slipperiness field and the Glen's law exponent field). However, under the assumption of Gaussian noise and prior probability densities, and after linearizing the parameter-to-observable map, the posterior density becomes Gaussian, and can therefore be characterized by its mean and covariance. The mean is given by the solution of a nonlinear least squares optimization problem, which is equivalent to a deterministic inverse problem with appropriate interpretation and weighting of the data misfit and regularization terms. To obtain this mean, we solve a deterministic ice sheet inverse problem; here, we infer parameters arising from discretizations of basal slipperiness and rheological exponent fields. For this purpose, we minimize a regularized misfit functional between observed and modeled surface flow velocities. The resulting least squares minimization problem is solved using an adjoint-based inexact Newton method, which uses first and second derivative information. The posterior covariance matrix is given (in the linear-Gaussian case) by the inverse of the Hessian of the least squares cost functional of the deterministic inverse problem. Direct computation of the Hessian matrix is prohibitive, since it would
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
A Bayesian method for microseismic source inversion
NASA Astrophysics Data System (ADS)
Pugh, D. J.; White, R. S.; Christie, P. A. F.
2016-08-01
Earthquake source inversion is highly dependent on location determination and velocity models. Uncertainties in both the model parameters and the observations need to be rigorously incorporated into an inversion approach. Here, we show a probabilistic Bayesian method that allows formal inclusion of the uncertainties in the moment tensor inversion. This method allows the combination of different sets of far-field observations, such as P-wave and S-wave polarities and amplitude ratios, into one inversion. Additional observations can be included by deriving a suitable likelihood function from the uncertainties. This inversion produces samples from the source posterior probability distribution, including a best-fitting solution for the source mechanism and associated probability. The inversion can be constrained to the double-couple space or allowed to explore the gamut of moment tensor solutions, allowing volumetric and other non-double-couple components. The posterior probability of the double-couple and full moment tensor source models can be evaluated from the Bayesian evidence, using samples from the likelihood distributions for the two source models, producing an estimate of whether or not a source is double-couple. Such an approach is ideally suited to microseismic studies where there are many sources of uncertainty and it is often difficult to produce reliability estimates of the source mechanism, although this can be true of many other cases. Using full-waveform synthetic seismograms, we also show the effects of noise, location, network distribution and velocity model uncertainty on the source probability density function. The noise has the largest effect on the results, especially as it can affect other parts of the event processing. This uncertainty can lead to erroneous non-double-couple source probability distributions, even when no other uncertainties exist. Although including amplitude ratios can improve the constraint on the source probability
Linear functional minimization for inverse modeling
NASA Astrophysics Data System (ADS)
Barajas-Solano, D. A.; Wohlberg, B. E.; Vesselinov, V. V.; Tartakovsky, D. M.
2015-06-01
We present a novel inverse modeling strategy to estimate spatially distributed parameters of nonlinear models. The maximum a posteriori (MAP) estimators of these parameters are based on a likelihood functional, which contains spatially discrete measurements of the system parameters and spatiotemporally discrete measurements of the transient system states. The piecewise continuity prior for the parameters is expressed via Total Variation (TV) regularization. The MAP estimator is computed by minimizing a nonquadratic objective equipped with the TV operator. We apply this inversion algorithm to estimate hydraulic conductivity of a synthetic confined aquifer from measurements of conductivity and hydraulic head. The synthetic conductivity field is composed of a low-conductivity heterogeneous intrusion into a high-conductivity heterogeneous medium. Our algorithm accurately reconstructs the location, orientation, and extent of the intrusion from the steady-state data only. Addition of transient measurements of hydraulic head improves the parameter estimation, accurately reconstructing the conductivity field in the vicinity of observation locations.
Bayesian inversion for optical diffraction tomography
NASA Astrophysics Data System (ADS)
Ayasso, H.; Duchêne, B.; Mohammad-Djafari, A.
2010-05-01
In this paper, optical diffraction tomography is considered as a non-linear inverse scattering problem and tackled within the Bayesian estimation framework. The object under test is a man-made object known to be composed of compact regions made of a finite number of different homogeneous materials. This a priori knowledge is appropriately translated by a Gauss-Markov-Potts prior. Hence, a Gauss-Markov random field is used to model the contrast distribution whereas a hidden Potts-Markov field accounts for the compactness of the regions. First, we express the a posteriori distributions of all the unknowns and then a Gibbs sampling algorithm is used to generate samples and estimate the posterior mean of the unknowns. Some preliminary results, obtained by applying the inversion algorithm to laboratory controlled data, are presented.
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
Regularity of mappings inverse to Sobolev mappings
Vodop'yanov, Sergei K
2012-10-31
For homeomorphisms {phi}:{Omega}{yields}{Omega}' on Euclidean domains in R{sup n}, n{>=}2, necessary and sufficient conditions ensuring that the inverse mapping belongs to a Sobolev class are investigated. The result obtained is used to describe a new two-index scale of homeomorphisms in some Sobolev class such that their inverses also form a two-index scale of mappings, in another Sobolev class. This scale involves quasiconformal mappings and also homeomorphisms in the Sobolev class W{sup 1}{sub n-1} such that rankD{phi}(x){<=}n-2 almost everywhere on the zero set of the Jacobian det D{phi}(x). Bibliography: 65 titles.
Earthquake source inversion with dense networks
NASA Astrophysics Data System (ADS)
Somala, S.; Ampuero, J. P.; Lapusta, N.
2012-12-01
Inversions of earthquake source slip from the recorded ground motions typically impose a number of restrictions on the source parameterization, which are needed to stabilize the inverse problem with sparse data. Such restrictions may include smoothing, causality considerations, predetermined shapes of the local source-time function, and constant rupture speed. The goal of our work is to understand whether the inversion results could be substantially improved by the availability of much denser sensor networks than currently available. The best regional networks have sensor spacing in the tens of kilometers range, much larger than the wavelengths relevant to key aspects of earthquake physics. Novel approaches to providing orders-of-magnitude denser sensing include low-cost sensors (Community Seismic Network) and space-based optical imaging (Geostationary Optical Seismometer). However, in both cases, the density of sensors comes at the expense of accuracy. Inversions that involve large number of sensors are intractable with the current source inversion codes. Hence we are developing a new approach that can handle thousands of sensors. It employs iterative conjugate gradient optimization based on an adjoint method and involves iterative time-reversed 3D wave propagation simulations using the spectral element method (SPECFEM3D). To test the developed method, and to investigate the effect of sensor density and quality on the inversion results, we have been considering kinematic and dynamic synthetic sources of several types: one or more Haskell pulses with various widths and spacings; scenarios with local rupture propagation in the opposite direction (as observed during the 2010 El Mayor-Cucapah earthquake); dynamic crack-like rupture, both subshear and supershear; and rupture that mimics supershear propagation by jumping along the fault. In each case, we produce the data by a forward SPECFEM3D calculation, choose the desired density of stations, filter the data to 1 Hz
Regeneration of stochastic processes: an inverse method
NASA Astrophysics Data System (ADS)
Ghasemi, F.; Peinke, J.; Sahimi, M.; Rahimi Tabar, M. R.
2005-10-01
We propose a novel inverse method that utilizes a set of data to construct a simple equation that governs the stochastic process for which the data have been measured, hence enabling us to reconstruct the stochastic process. As an example, we analyze the stochasticity in the beat-to-beat fluctuations in the heart rates of healthy subjects as well as those with congestive heart failure. The inverse method provides a novel technique for distinguishing the two classes of subjects in terms of a drift and a diffusion coefficients which behave completely differently for the two classes of subjects, hence potentially providing a novel diagnostic tool for distinguishing healthy subjects from those with congestive heart failure, even at the early stages of the disease development.
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.
An efficient method for inverse problems
NASA Technical Reports Server (NTRS)
Daripa, Prabir
1987-01-01
A new inverse method for aerodynamic design of subcritical airfoils is presented. The pressure distribution in this method can be prescribed in a natural way, i.e. as a function of arclength of the as yet unknown body. This inverse problem is shown to be mathematically equivalent to solving a single nonlinear boundary value problem subject to known Dirichlet data on the boundary. The solution to this problem determines the airfoil, the free stream Mach number M(sub x) and the upstream flow direction theta(sub x). The existence of a solution for any given pressure distribution is discussed. The method is easy to implement and extremely efficient. We present a series of results for which comparisons are made with the known airfoils.
Space-time inversion and its consequences
NASA Astrophysics Data System (ADS)
Chelnokov, M.
2016-07-01
The article discusses some new aspects of both the inversion of space and the inversion of time. It is shown that behind the mirror is not symmetric to in front of the mirror. It, in its turn, leads to nonconservation of spatial parity. The same situation takes place in the combined CP-parity. Further, the article shows that from the point of view of different reference systems of the Universe (from the point of view of different galaxies or accumulation of galaxies) time flows not just differently, and, in some cases, in the opposite directions. It leads to major changes in the picture of the Universe. In particular, the concept of the age of the Universe loses its meaning, serious doubts about the idea of the Big Bang and so on.
New advances in Inverse Cerenkov acceleration
Kimura, W.D.; Babzien, M.; Cline, D.B.; Fiorito, R.B.; Fontana, J.R.; Gallardo, J.C.; Gottschalk, S.C.; Kusche, K.P. |; Liu, Y.; Pogorelsky, I.V.; Quimby, D.C.; Pantell, R.H.; Rule, D.W.; Skaritka, J.; Sandweiss, J.; van Steenbergen, A.; Yakimenko, V.
1997-02-01
Inverse Cerenkov acceleration (ICA) has entered a new phase in its development. The issue of staging and rephasing the optical wave with a microbunched electron beam is now being examined. This ability to accelerate over multiple stages is important for scaling laser accelerator devices to higher energies. An inverse free electron laser (IFEL) will be positioned upstream from the ICA experiment and used to prebunch the electrons. These electrons will then be focused into the ICA interaction region for rephasing and acceleration by the laser beam. Issues that will be examined during these combined ICA/IFEL experiments include rephasing the laser beam with the microbunches, minimizing bunch smearing, and trapping the electrons in an acceleration bucket. {copyright} {ital 1997 American Institute of Physics.}
Iterative image restoration using approximate inverse preconditioning.
Nagy, J G; Plemmons, R J; Torgersen, T C
1996-01-01
Removing a linear shift-invariant blur from a signal or image can be accomplished by inverse or Wiener filtering, or by an iterative least-squares deblurring procedure. Because of the ill-posed characteristics of the deconvolution problem, in the presence of noise, filtering methods often yield poor results. On the other hand, iterative methods often suffer from slow convergence at high spatial frequencies. This paper concerns solving deconvolution problems for atmospherically blurred images by the preconditioned conjugate gradient algorithm, where a new approximate inverse preconditioner is used to increase the rate of convergence. Theoretical results are established to show that fast convergence can be expected, and test results are reported for a ground-based astronomical imaging problem. PMID:18285203
Approximate inverse preconditioners for general sparse matrices
Chow, E.; Saad, Y.
1994-12-31
Preconditioned Krylov subspace methods are often very efficient in solving sparse linear matrices that arise from the discretization of elliptic partial differential equations. However, for general sparse indifinite matrices, the usual ILU preconditioners fail, often because of the fact that the resulting factors L and U give rise to unstable forward and backward sweeps. In such cases, alternative preconditioners based on approximate inverses may be attractive. We are currently developing a number of such preconditioners based on iterating on each column to get the approximate inverse. For this approach to be efficient, the iteration must be done in sparse mode, i.e., we must use sparse-matrix by sparse-vector type operatoins. We will discuss a few options and compare their performance on standard problems from the Harwell-Boeing collection.
Study of modular inversion in RNS
NASA Astrophysics Data System (ADS)
Bajard, Jean Claude; Meloni, Nicolas; Plantard, Thomas
2005-08-01
Residue Numbers System have some features which are fine for some implementations of cryptographic protocols. The main property of RNS is the distribution of the evaluation on large values on its small residues, allowing parallelization. This last property implies that we can randomize the distribution of the bases elements. Hence, the obtained arithmetic is leak resistant, it is robust against side channel attacks. But one drawback of RNS is that modular inversion is not obvious. Thus, RNS is well suited for RSA but not really for ECC. We analyze in this paper the features of the modular inversion in RNS over GF(P). We propose a RNS Extended Euclidean Algorithm which uses a quotient approximation module.
Optical inverse-square displacement sensor
Howe, Robert D.; Kychakoff, George
1989-01-01
This invention comprises an optical displacement sensor that uses the inverse-square attenuation of light reflected from a diffused surface to calculate the distance from the sensor to the reflecting surface. Light emerging from an optical fiber or the like is directed onto the surface whose distance is to be measured. The intensity I of reflected light is angle dependent, but within a sufficiently small solid angle it falls off as the inverse square of the distance from the surface. At least a pair of optical detectors are mounted to detect the reflected light within the small solid angle, their ends being at different distances R and R+.DELTA.R from the surface. The distance R can then be found in terms of the ratio of the intensity measurements and the separation length as ##EQU1##
Optical inverse-square displacement sensor
Howe, R.D.; Kychakoff, G.
1989-09-12
This invention comprises an optical displacement sensor that uses the inverse-square attenuation of light reflected from a diffused surface to calculate the distance from the sensor to the reflecting surface. Light emerging from an optical fiber or the like is directed onto the surface whose distance is to be measured. The intensity I of reflected light is angle dependent, but within a sufficiently small solid angle it falls off as the inverse square of the distance from the surface. At least a pair of optical detectors are mounted to detect the reflected light within the small solid angle, their ends being at different distances R and R + [Delta]R from the surface. The distance R can then be found in terms of the ratio of the intensity measurements and the separation length as given in an equation. 10 figs.
Inverse free-electron laser accelerator
Pellegrini, C.; Campisi, R.
1982-01-01
We first describe the basic physical properties of an inverse free-electron laser and make an estimate of the order of magnitude of the accelerating field obtainable with such a system; then apply the general ideas to the design of an actual device and through this example we give a more accurate evaluation of the fundamental as well as the technical limitations that this acceleration scheme imposes.
Inversion identities for inhomogeneous face models
NASA Astrophysics Data System (ADS)
Frahm, Holger; Karaiskos, Nikos
2014-10-01
We derive exact inversion identities satisfied by the transfer matrix of inhomogeneous interaction-round-a-face (IRF) models with arbitrary boundary conditions using the underlying integrable structure and crossing properties of the local Boltzmann weights. For the critical restricted solid-on-solid (RSOS) models these identities together with some information on the analytical properties of the transfer matrix determine the spectrum completely and allow to derive the Bethe equations for both periodic and general open boundary conditions.
Inverse Design of a Thick Supercritical Airfoil
NASA Astrophysics Data System (ADS)
Pambagjo, Tjoetjoek Eko; Nakahashi, Kazuhiro; Obayashi, Shigeru
In this paper, a study on designing a thick supercritical airfoil by utilizing Takanashi’s inverse design method is discussed. One of the problems to design a thick supercritical airfoil by Takanashi’s method is that an oscillation of the geometry may occur during the iteration process. To reduce the oscillation, an airfoil parameterization method is utilized as the smoothing procedure. A guideline to determine the target pressure distribution to realize the thick airfoil is also discussed.
Inverse scattering problem for quantum graph vertices
Cheon, Taksu; Turek, Ondrej; Exner, Pavel
2011-06-15
We demonstrate how the inverse scattering problem of a quantum star graph can be solved by means of diagonalization of the Hermitian unitary matrix when the vertex coupling is of the scale-invariant (or Fueloep-Tsutsui) form. This enables the construction of quantum graphs with desired properties in a tailor-made fashion. The procedure is illustrated on the example of quantum vertices with equal transmission probabilities.
3D Inverse problem: Seawater intrusions
NASA Astrophysics Data System (ADS)
Steklova, K.; Haber, E.
2013-12-01
Modeling of seawater intrusions (SWI) is challenging as it involves solving the governing equations for variable density flow, multiple time scales and varying boundary conditions. Due to the nonlinearity of the equations and the large aquifer domains, 3D computations are a costly process, particularly when solving the inverse SWI problem. In addition the heads and concentration measurements are difficult to obtain due to mixing, saline wedge location is sensitive to aquifer topography, and there is general uncertainty in initial and boundary conditions and parameters. Some of these complications can be overcome by using indirect geophysical data next to standard groundwater measurements, however, the inverse problem is usually simplified, e.g. by zonation for the parameters based on geological information, steady state substitution of the unknown initial conditions, decoupling the equations or reducing the amount of unknown parameters by covariance analysis. In our work we present a discretization of the flow and solute mass balance equations for variable groundwater (GW) flow. A finite difference scheme is to solve pressure equation and a Semi - Lagrangian method for solute transport equation. In this way we are able to choose an arbitrarily large time step without losing stability up to an accuracy requirement coming from the coupled character of the variable density flow equations. We derive analytical sensitivities of the GW model for parameters related to the porous media properties and also the initial solute distribution. Analytically derived sensitivities reduce the computational cost of inverse problem, but also give insight for maximizing information in collected data. If the geophysical data are available it also enables simultaneous calibration in a coupled hydrogeophysical framework. The 3D inverse problem was tested on artificial time dependent data for pressure and solute content coming from a GW forward model and/or geophysical forward model. Two
Numerical linear algebra for reconstruction inverse problems
NASA Astrophysics Data System (ADS)
Nachaoui, Abdeljalil
2004-01-01
Our goal in this paper is to discuss various issues we have encountered in trying to find and implement efficient solvers for a boundary integral equation (BIE) formulation of an iterative method for solving a reconstruction problem. We survey some methods from numerical linear algebra, which are relevant for the solution of this class of inverse problems. We motivate the use of our constructing algorithm, discuss its implementation and mention the use of preconditioned Krylov methods.
Forecast Variance Estimates Using Dart Inversion
NASA Astrophysics Data System (ADS)
Gica, E.
2014-12-01
The tsunami forecast tool developed by the NOAA Center for Tsunami Research (NCTR) provides real-time tsunami forecast and is composed of the following major components: a pre-computed tsunami propagation database, an inversion algorithm that utilizes real-time tsunami data recorded at DART stations to define the tsunami source, and inundation models that predict tsunami wave characteristics at specific coastal locations. The propagation database is a collection of basin-wide tsunami model runs generated from 50x100 km "unit sources" with a slip of 1 meter. Linear combination and scaling of unit sources is possible since the nonlinearity in the deep ocean is negligible. To define the tsunami source using the unit sources, real-time DART data is ingested into an inversion algorithm. Based on the selected DART and length of tsunami time series, the inversion algorithm will select the best combination of unit sources and scaling factors that best fit the observed data at the selected locations. This combined source then serves as boundary condition for the inundation models. Different combinations of DARTs and length of tsunami time series used in the inversion algorithm will result in different selection of unit sources and scaling factors. Since the combined unit sources are used as boundary condition for inundation modeling, different sources will produce variations in the tsunami wave characteristics. As part of the testing procedures for the tsunami forecast tool, staff at NCTR and both National and Pacific Tsunami Warning Centers, performed post-event forecasts for several historical tsunamis. The extent of variation due to different source definitions obtained from the testing is analyzed by comparing the simulated maximum tsunami wave amplitude with recorded data at tide gauge locations. Results of the analysis will provide an error estimate defining the possible range of the simulated maximum tsunami wave amplitude for each specific inundation model.
Pole EXpansion and Selected Inversion (PEXSI)
Energy Science and Technology Software Center (ESTSC)
2014-03-01
The Pole EXpansion and Selected Inversion method (PEXSI) is a fast method for evaluating certain selected elements of a matrix function. PEXSI is highly scalable on distributed memory parallel machines. For sparse matrices, the PEXSI method can be more efficient than the widely used diagonalization method for evaluating matrix functions, especially when a relatively large number of eigenpairs are needed to be computed in the diagonalization methond
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.
Parametric optimization of inverse trapezoid oleophobic surfaces.
Cavalli, Andrea; Bøggild, Peter; Okkels, Fridolin
2012-12-18
In this paper, we introduce a comprehensive and versatile approach to the parametric shape optimization of oleophobic surfaces. We evaluate the performance of inverse trapezoid microstructures in terms of three objective parameters: apparent contact angle, maximum sustainable hydrostatic pressure, and mechanical robustness (Im, M.; Im, H:; Lee, J.H.; Yoon, J.B.; Choi, Y.K. A robust superhydrophobic and superoleophobic surface with inverse-trapezoidal microstructures on a large transparent flexible substrate. Soft Matter 2010, 6, 1401-1404; Im, M.; Im, H:; Lee, J.H.; Yoon, J.B.; Choi, Y.K. Analytical Modeling and Thermodynamic Analysis of Robust Superhydrophobic Surfaces with Inverse-Trapezoidal Microstructures. Langmuir 2010, 26, 17389-17397). We find that each of these parameters, if considered alone, would give trivial optima, while their interplay provides a well-defined optimal shape and aspect ratio. The inclusion of mechanical robustness in combination with conventional performance characteristics favors solutions relevant for practical applications, as mechanical stability is a critical issue not often addressed in idealized models. PMID:23078017
Constraining inverse curvature gravity with supernovae
Mena, Olga; Santiago, Jose; Weller, Jochen; /University Coll., London /Fermilab
2005-10-01
We show that the current accelerated expansion of the Universe can be explained without resorting to dark energy. Models of generalized modified gravity, with inverse powers of the curvature can have late time accelerating attractors without conflicting with solar system experiments. We have solved the Friedman equations for the full dynamical range of the evolution of the Universe. This allows us to perform a detailed analysis of Supernovae data in the context of such models that results in an excellent fit. Hence, inverse curvature gravity models represent an example of phenomenologically viable models in which the current acceleration of the Universe is driven by curvature instead of dark energy. If we further include constraints on the current expansion rate of the Universe from the Hubble Space Telescope and on the age of the Universe from globular clusters, we obtain that the matter content of the Universe is 0.07 {le} {omega}{sub m} {le} 0.21 (95% Confidence). Hence the inverse curvature gravity models considered can not explain the dynamics of the Universe just with a baryonic matter component.
Preparatiion of metal colloids in inverse micelles
Wilcoxon, J.P.
1990-11-23
A method is provided for preparing catalytic elemental metal colloidal particles (e.g., gold, palladium, silver, rhodium, nickel, iron, platinum, molybdenum) or colloidal alloy particles (silver/iridium or platinum/gold). A homogenous inverse micelle solution of a metal salt is first formed in a metal-salt solvent comprised of a surfactant (e.g. a nonionic or cationic surfactant) and an organic solvent. The size and number of inverse micelles is controlled by the proportions of the surfactant and the solvent. Then, the metal salt is reduced (by chemical reduction or by a pulsed or continuous wave UV laser) to colloidal particles of elemental metal. After their formation, the colloidal metal particles can be stabilized by reaction with materials that permanently add surface stabilizing groups to the surface of the colloidal metal particles. The sizes of the colloidal elemental metal particles and their size distribution is determined by the size and number of the inverse micelles. A second salt can be added with further reduction to form the colloidal alloy particles. After the colloidal elemental metal particles are formed, the homogeneous solution distributes to two phases, one phase rich in colloidal elemental metal particles and the other phase rich in surfactant. The colloidal elemental metal particles from one phase can be dried to form a powder useful as a catalyst.
Forward and Inverse Cascades in EMHD Turbulence
NASA Astrophysics Data System (ADS)
Cho, Jungyeon
2016-05-01
Electron magnetohydrodynamics (EMHD) provides a simple fluid-like description of physics below the proton gyro-scale in collisionless plasmas, such as the solar wind. In this paper, we discuss forward and inverse cascades in EMHD turbulence in the presence of a strong mean magnetic field. Similar to Alfvén waves, EMHD waves, or EMHD perturbations, propagate along magnetic field lines. Therefore, two types of EMHD waves can exist: waves moving parallel to and waves moving anti-parallel to the the magnetic field lines. For energy cascade in EMHD turbulence, the relative amplitudes of opposite-traveling waves are important. When the amplitudes are balanced, we will see fully-developed forward cascade with a k -7/3 energy spectrum and a scale-dependent anisotropy. On the other hand, when the amplitudes are imbalanced, we will see inverse cascade, as well as (presumably not fully developed) forward cascade. The underlying physics for the inverse cascade is magnetic helicity conservation.
Inverse sequential simulation: Performance and implementation details
NASA Astrophysics Data System (ADS)
Xu, Teng; Gómez-Hernández, J. Jaime
2015-12-01
For good groundwater flow and solute transport numerical modeling, it is important to characterize the formation properties. In this paper, we analyze the performance and important implementation details of a new approach for stochastic inverse modeling called inverse sequential simulation (iSS). This approach is capable of characterizing conductivity fields with heterogeneity patterns difficult to capture by standard multiGaussian-based inverse approaches. The method is based on the multivariate sequential simulation principle, but the covariances and cross-covariances used to compute the local conditional probability distributions are computed by simple co-kriging which are derived from an ensemble of conductivity and piezometric head fields, in a similar manner as the experimental covariances are computed in an ensemble Kalman filtering. A sensitivity analysis is performed on a synthetic aquifer regarding the number of members of the ensemble of realizations, the number of conditioning data, the number of piezometers at which piezometric heads are observed, and the number of nodes retained within the search neighborhood at the moment of computing the local conditional probabilities. The results show the importance of having a sufficiently large number of all of the mentioned parameters for the algorithm to characterize properly hydraulic conductivity fields with clear non-multiGaussian features.
Analysis of Raman lasing without inversion
NASA Astrophysics Data System (ADS)
Sheldon, Paul Martin
1999-12-01
Properties of lasing without inversion were studied analytically and numerically using Maple computer assisted algebra software. Gain for probe electromagnetic field without population inversion in detuned three level atomic schemes has been found. Matter density matrix dynamics and coherence is explored using Pauli matrices in 2-level systems and Gell-Mann matrices in 3-level systems. It is shown that extreme inversion produces no coherence and hence no lasing. Unitary transformation from the strict field-matter Hamiltonian to an effective two-photon Raman Hamiltonian for multilevel systems has been derived. Feynman diagrams inherent in the derivation show interesting physics. An additional picture change was achieved and showed cw gain possible. Properties of a Raman-like laser based on injection of 3- level coherently driven Λ-type atoms whose Hamiltonian contains the Raman Hamiltonian and microwave coupling the two bottom states have been studied in the limits of small and big photon numbers in the drive field. Another picture change removed the microwave coupler to all orders and simplified analysis. New possibilities of inversionless generation were found.
Combinatorial approaches for inverse metabolic engineering applications
Skretas, Georgios; Kolisis, Fragiskos N.
2013-01-01
Traditional metabolic engineering analyzes biosynthetic and physiological pathways, identifies bottlenecks, and makes targeted genetic modifications with the ultimate goal of increasing the production of high-value products in living cells. Such efforts have led to the development of a variety of organisms with industrially relevant properties. However, there are a number of cellular phenotypes important for research and the industry for which the rational selection of cellular targets for modification is not easy or possible. In these cases, strain engineering can be alternatively carried out using “inverse metabolic engineering”, an approach that first generates genetic diversity by subjecting a population of cells to a particular mutagenic process, and then utilizes genetic screens or selections to identify the clones exhibiting the desired phenotype. Given the availability of an appropriate screen for a particular property, the success of inverse metabolic engineering efforts usually depends on the level and quality of genetic diversity which can be generated. Here, we review classic and recently developed combinatorial approaches for creating such genetic diversity and discuss the use of these methodologies in inverse metabolic engineering applications. PMID:24688681
Computationally efficient Bayesian inference for inverse problems.
Marzouk, Youssef M.; Najm, Habib N.; Rahn, Larry A.
2007-10-01
Bayesian statistics provides a foundation for inference from noisy and incomplete data, a natural mechanism for regularization in the form of prior information, and a quantitative assessment of uncertainty in the inferred results. Inverse problems - representing indirect estimation of model parameters, inputs, or structural components - can be fruitfully cast in this framework. Complex and computationally intensive forward models arising in physical applications, however, can render a Bayesian approach prohibitive. This difficulty is compounded by high-dimensional model spaces, as when the unknown is a spatiotemporal field. We present new algorithmic developments for Bayesian inference in this context, showing strong connections with the forward propagation of uncertainty. In particular, we introduce a stochastic spectral formulation that dramatically accelerates the Bayesian solution of inverse problems via rapid evaluation of a surrogate posterior. We also explore dimensionality reduction for the inference of spatiotemporal fields, using truncated spectral representations of Gaussian process priors. These new approaches are demonstrated on scalar transport problems arising in contaminant source inversion and in the inference of inhomogeneous material or transport properties. We also present a Bayesian framework for parameter estimation in stochastic models, where intrinsic stochasticity may be intermingled with observational noise. Evaluation of a likelihood function may not be analytically tractable in these cases, and thus several alternative Markov chain Monte Carlo (MCMC) schemes, operating on the product space of the observations and the parameters, are introduced.
Radiation-Insensitive Inverse Majority Gates
NASA Technical Reports Server (NTRS)
Manohara, Harish; Mojarradi, Mohammad
2008-01-01
To help satisfy a need for high-density logic circuits insensitive to radiation, it has been proposed to realize inverse majority gates as microscopic vacuum electronic devices. In comparison with solid-state electronic devices ordinarily used in logic circuits, vacuum electronic devices are inherently much less adversely affected by radiation and extreme temperatures. The proposed development would involve state-of-the-art micromachining and recent advances in the fabrication of carbon-nanotube-based field emitters. A representative three-input inverse majority gate would be a monolithic, integrated structure that would include three gate electrodes, six bundles of carbon nanotubes (serving as electron emitters) at suitable positions between the gate electrodes, and an overhanging anode. The bundles of carbon nanotubes would be grown on degenerately doped silicon substrates that would be parts of the monolithic structure. The gate electrodes would be fabricated as parts of the monolithic structure by means of a double-silicon-on-insulator process developed at NASA's Jet Propulsion Laboratory. The tops of the bundles of carbon nanotubes would lie below the plane of the tops of the gate electrodes. The particular choice of shapes, dimensions, and relative positions of the electrodes and bundles of carbon nanotubes would provide for both field emission of electrons from the bundles of carbon nanotubes and control of the electron current to obtain the inverse majority function, which is described in the paper.
Inverse hydrochemical models of aqueous extracts tests
Zheng, L.; Samper, J.; Montenegro, L.
2008-10-10
Aqueous extract test is a laboratory technique commonly used to measure the amount of soluble salts of a soil sample after adding a known mass of distilled water. Measured aqueous extract data have to be re-interpreted in order to infer porewater chemical composition of the sample because porewater chemistry changes significantly due to dilution and chemical reactions which take place during extraction. Here we present an inverse hydrochemical model to estimate porewater chemical composition from measured water content, aqueous extract, and mineralogical data. The model accounts for acid-base, redox, aqueous complexation, mineral dissolution/precipitation, gas dissolution/ex-solution, cation exchange and surface complexation reactions, of which are assumed to take place at local equilibrium. It has been solved with INVERSE-CORE{sup 2D} and been tested with bentonite samples taken from FEBEX (Full-scale Engineered Barrier EXperiment) in situ test. The inverse model reproduces most of the measured aqueous data except bicarbonate and provides an effective, flexible and comprehensive method to estimate porewater chemical composition of clays. Main uncertainties are related to kinetic calcite dissolution and variations in CO2(g) pressure.
Inverse problems biomechanical imaging (Conference Presentation)
NASA Astrophysics Data System (ADS)
Oberai, Assad A.
2016-03-01
It is now well recognized that a host of imaging modalities (a list that includes Ultrasound, MRI, Optical Coherence Tomography, and optical microscopy) can be used to "watch" tissue as it deforms in response to an internal or external excitation. The result is a detailed map of the deformation field in the interior of the tissue. This deformation field can be used in conjunction with a material mechanical response to determine the spatial distribution of material properties of the tissue by solving an inverse problem. Images of material properties thus obtained can be used to quantify the health of the tissue. Recently, they have been used to detect, diagnose and monitor cancerous lesions, detect vulnerable plaque in arteries, diagnose liver cirrhosis, and possibly detect the onset of Alzheimer's disease. In this talk I will describe the mathematical and computational aspects of solving this class of inverse problems, and their applications in biology and medicine. In particular, I will discuss the well-posedness of these problems and quantify the amount of displacement data necessary to obtain a unique property distribution. I will describe an efficient algorithm for solving the resulting inverse problem. I will also describe some recent developments based on Bayesian inference in estimating the variance in the estimates of material properties. I will conclude with the applications of these techniques in diagnosing breast cancer and in characterizing the mechanical properties of cells with sub-cellular resolution.
Magnetic Resonance Elastography: Inversions in Bounded Media
Kolipaka, Arunark; McGee, Kiaran P.; Manduca, Armando; Romano, Anthony J.; Glaser, Kevin J.; Araoz, Philip A.; Ehman, Richard L.
2009-01-01
Magnetic resonance elastography (MRE) is a noninvasive imaging technique capable of quantifying and spatially resolving the shear stiffness of soft tissues by visualization of synchronized mechanical wave displacement fields. However, MRE inversions generally assume that the measured tissue motion consists primarily of shear waves propagating in a uniform, infinite medium. This assumption is not valid in organs such as the heart, eye, bladder, skin, fascia, bone and spinal cord in which the shear wavelength approaches the geometric dimensions of the object. The aim of this study was to develop and test mathematical inversion algorithms capable of resolving shear stiffness from displacement maps of flexural waves propagating in bounded media such as beams, plates and spherical shells using geometry-specific equations of motion. MRE and finite element modeling (FEM) of beam, plate, and spherical shell phantoms of various geometries were performed. Mechanical testing of the phantoms agreed with the stiffness values obtained from FEM and MRE data and a linear correlation of r2 ≥ 0.99 was observed between the stiffness values obtained using MRE and FEM data. In conclusion, we have demonstrated new inversion methods for calculating shear stiffness that may be more appropriate for waves propagating in bounded media. PMID:19780146
An Inverse Approach for Elucidating Dendritic Function
Torben-Nielsen, Benjamin; Stiefel, Klaus M.
2010-01-01
We outline an inverse approach for investigating dendritic function–structure relationships by optimizing dendritic trees for a priori chosen computational functions. The inverse approach can be applied in two different ways. First, we can use it as a “hypothesis generator” in which we optimize dendrites for a function of general interest. The optimization yields an artificial dendrite that is subsequently compared to real neurons. This comparison potentially allows us to propose hypotheses about the function of real neurons. In this way, we investigated dendrites that optimally perform input-order detection. Second, we can use it as a “function confirmation” by optimizing dendrites for functions hypothesized to be performed by classes of neurons. If the optimized, artificial, dendrites resemble the dendrites of real neurons the artificial dendrites corroborate the hypothesized function of the real neuron. Moreover, properties of the artificial dendrites can lead to predictions about yet unmeasured properties. In this way, we investigated wide-field motion integration performed by the VS cells of the fly visual system. In outlining the inverse approach and two applications, we also elaborate on the nature of dendritic function. We furthermore discuss the role of optimality in assigning functions to dendrites and point out interesting future directions. PMID:21258425
Chemical inversion in the subsurface hydrosphere
Yezhov, Yu.A.
1980-09-01
A quite common nature of chemical inversion in subsurface hydrosphere is shown in examples of several oil- and gas-bearing regions of the USSR. In particular, when the data of sampling from deep wells of the Volgo-Urals, Mangyshlak, and Western Turkmenian regions were compared, it became obvious that the composite chemical profile of subsurface hydrosphere consists of a vertical alternation of three zones: of increasing (I-II-IIIa genetic types of subsurface waters), maximum (IIIb), and decreasing water mineralization (III'a-II'-I'). The depth of occurrence of the lower inversion branch of zonality depends on the geotectonic activity at depth. It is closer to the Earth's surface in regions of Alpine tectogenesis, whereas in regions of ancient folding it lies at great depths which have not yet been reached by most deep wells. The formation of the inversion zone in the Earth's crust is connected with penetration from below ascending demineralized fluids of sodium bicarbonate type (I'). The latter is due to the presence at great depths of large quantities of free carbonic acid which is involved in hydrolytic processes of decomposition of sodium-containing minerals and produces sodium-type waters.
Geomechanical paleostress inversion using fracture data
NASA Astrophysics Data System (ADS)
Maerten, Laurent; Maerten, Frantz; Lejri, Mostfa; Gillespie, Paul
2016-08-01
We describe a fast geomechanically-based paleostress inversion technique that uses observed fracture data to constrain stress through multiple simulations. The method assumes that the local stress field around individual fractures is heterogeneous and derives the far field tectonic stress, that we also call the far field boundary conditions. We show how such far field tectonic stress can be recovered through a mechanical stress inversion technique using local observations of natural fractures (i.e. mechanical type, orientation and location). We test the paleostress inversion against outcrop analogues of fractured carbonates from both Nash Point, U.K., where there are well exposed faults and joints and the Matelles, France, where there are well exposed faults, veins and stylolites. We demonstrate through these case studies how the method can be efficiently applied to natural examples and we highlight its advantages and limitations. We discuss how such method could be applied to subsurface problems and how it can provide complementary constraints to drive discrete fracture models for better fractured reservoir characterization and modelling.
An overview of joint inversion in earthquake source imaging
NASA Astrophysics Data System (ADS)
Koketsu, Kazuki
2016-06-01
We reviewed joint inversion studies of the rupture processes of significant earthquakes, using the definition of a joint inversion in earthquake source imaging as a source inversion of multiple kinds of datasets (waveform, geodetic, or tsunami). Yoshida and Koketsu (Geophys J Int 103:355-362, 1990), and Wald and Heaton (Bull Seismol Soc Am 84:668-691, 1994) independently initiated joint inversion methods, finding that joint inversion provides more reliable rupture process models than single-dataset inversion, leading to an increase of joint inversion studies. A list of these studies was made using the finite-source rupture model database (Mai and Thingbaijam in Seismol Res Lett 85:1348-1357, 2014). Outstanding issues regarding joint inversion were also discussed.
Approximate inverse preconditioning of iterative methods for nonsymmetric linear systems
Benzi, M.; Tuma, M.
1996-12-31
A method for computing an incomplete factorization of the inverse of a nonsymmetric matrix A is presented. The resulting factorized sparse approximate inverse is used as a preconditioner in the iterative solution of Ax = b by Krylov subspace methods.
Magnetotelluric inversion via reverse time migration algorithm of seismic data
Ha, Taeyoung . E-mail: tyha@math.snu.ac.kr; Shin, Changsoo . E-mail: css@model.snu.ac.kr
2007-07-01
We propose a new algorithm for two-dimensional magnetotelluric (MT) inversion. Our algorithm is an MT inversion based on the steepest descent method, borrowed from the backpropagation technique of seismic inversion or reverse time migration, introduced in the middle 1980s by Lailly and Tarantola. The steepest descent direction can be calculated efficiently by using the symmetry of numerical Green's function derived from a mixed finite element method proposed by Nedelec for Maxwell's equation, without calculating the Jacobian matrix explicitly. We construct three different objective functions by taking the logarithm of the complex apparent resistivity as introduced in the recent waveform inversion algorithm by Shin and Min. These objective functions can be naturally separated into amplitude inversion, phase inversion and simultaneous inversion. We demonstrate our algorithm by showing three inversion results for synthetic data.
Liquid-impermeable inverse opals with invariant photonic bandgap.
Kang, Hyelim; Lee, Joon-Seok; Chang, Won Seok; Kim, Shin-Hyun
2015-02-18
Omniphobic inverse opals are created by structurally and chemically modifying the surface of inverse opals through reactive ion etching. During the etching, void arrays of the inverse opal surface evolves to a triangular post array with re-entrant geometry. The elaborate structure can efficiently pin the air-liquid interface and retain air cavities against water and oil, thereby providing liquid-impermeable inverse opals with invariant photonic bandgap. PMID:25492694
Direct inversion of rigid-body rotational dynamics
NASA Technical Reports Server (NTRS)
Bach, Ralph; Paielli, Russell
1990-01-01
The global linearization (inversion) of rigid-body rotational dynamics is reviewed and representations in terms of quaternions and direction cosines are compared. Certain properties common to quaternions and direction cosines that make their use preferable to Euler angles and that simplify the inversion procedure are described. Applications of the inversion procedure for state estimation and attitude control are discussed. To avoid complexities caused by aerodynamics, an example of direct inversion for linear feedback control of spacecraft attitude is given.
Full Waveform Inversion Using Waveform Sensitivity Kernels
NASA Astrophysics Data System (ADS)
Schumacher, Florian; Friederich, Wolfgang
2013-04-01
We present a full waveform inversion concept for applications ranging from seismological to enineering contexts, in which the steps of forward simulation, computation of sensitivity kernels, and the actual inversion are kept separate of each other. We derive waveform sensitivity kernels from Born scattering theory, which for unit material perturbations are identical to the Born integrand for the considered path between source and receiver. The evaluation of such a kernel requires the calculation of Green functions and their strains for single forces at the receiver position, as well as displacement fields and strains originating at the seismic source. We compute these quantities in the frequency domain using the 3D spectral element code SPECFEM3D (Tromp, Komatitsch and Liu, 2008) and the 1D semi-analytical code GEMINI (Friederich and Dalkolmo, 1995) in both, Cartesian and spherical framework. We developed and implemented the modularized software package ASKI (Analysis of Sensitivity and Kernel Inversion) to compute waveform sensitivity kernels from wavefields generated by any of the above methods (support for more methods is planned), where some examples will be shown. As the kernels can be computed independently from any data values, this approach allows to do a sensitivity and resolution analysis first without inverting any data. In the context of active seismic experiments, this property may be used to investigate optimal acquisition geometry and expectable resolution before actually collecting any data, assuming the background model is known sufficiently well. The actual inversion step then, can be repeated at relatively low costs with different (sub)sets of data, adding different smoothing conditions. Using the sensitivity kernels, we expect the waveform inversion to have better convergence properties compared with strategies that use gradients of a misfit function. Also the propagation of the forward wavefield and the backward propagation from the receiver
Inverse Functions: What Our Teachers Didn't Tell Us
ERIC Educational Resources Information Center
Wilson, Frank C.; Adamson, Scott; Cox, Trey; O'Bryan, Alan
2011-01-01
The mathematical topic of inverse functions is an important element of algebra courses at the high school and college levels. The inverse function concept is best understood by students when it is presented in a familiar, real-world context. In this article, the authors discuss some misconceptions about inverse functions and suggest some…
Radiation-induced chromosomal inversions in mice. Technical progress report
Roderick, T.H.
1986-01-01
Chromosomal inversions are being produced for the purpose of establishing efficient systems for assessing induced and spontaneous heritable mutations. The inversions and other chromosomal aberrations produced are used to ask basic questions about meiosis and reproductive performance. Chromosomal structure is being studied by identifying the cytological location of genes and break points related to the inversions. 2 tabs.
NASA Astrophysics Data System (ADS)
Martínez, M. D.; Lana, X.
1991-03-01
The total inversion algorithm and some elements of Mathematical Information Theory are used in the treatment of travel-time data belonging to a seismic refraction experiment from the southern segment (Sardinia Channel) of the European Geotraverse Project. The inversion algorithm allows us to improve a preliminary propagating model obtained by means of usual trial and error procedure and to quantify the resolution degree of parameters defining the crust and upper mantle of such a model. Concepts related to Mathematical Information Theory detect some seismic profiles of the refraction experiment which give the most homogeneous coverage of the model in terms of number of trajectories crossing it. Finally, the efficiency of the inversion procedure is quantified and the uncertainties regarding knowledge of different parts of the model are also evaluated.
Microwave spectrum of the HD2O+ ion: inversion-rotation transitions and inversion splitting.
Furuya, Takashi; Saito, Shuji
2008-01-21
Inversion-rotation spectral lines of the dideuterated hydronium ion, HD2O+, have been observed by a source-modulation millimeter- to submillimeter-wave spectrometer. The ion was generated by a hollow-cathode discharge in a gas mixture of D2O and H2O in a free-space cell. Ten inversion-rotation lines were measured precisely for the lowest pair of inversion doublets in the frequency region from 380 to 730 GHz. The observed lines include the most astronomically important transitions, 0(00) (-)-1(10)+ for the para species at 380 538.031(32) MHz and 1(01) (-)-1(11)+ for the ortho species at 728 420.189(34) MHz, which could be used as a radio astronomical probe investigating interstellar chemistry of deuterium fractionation. An analysis of the measured lines has yielded the rotational constants in the ground doublet states and the inversion splitting. The inversion splitting in the ground state was determined to be 808 866(34) MHz, that is, 26.980 87(113) cm(-1), where the numbers in parentheses give uncertainties estimated from the Jacobian matrix of the assumed centrifugal distortion constants. The determined inversion splitting is off by -0.51 cm(-1) from the predicted value of 27.49 cm(-1) by Rajamaki et al. using high-order coupled cluster ab initio calculation [J. Chem. Phys. 118, 10929 (2003)], and by -0.0510 cm(-1) from the observed value of 27.0318(72) cm(-1) by Dong et al. using high-resolution jet-cooled infrared spectroscopy [J. Chem. Phys. 122, 224301 (2005)] beyond the quoted uncertainty. PMID:18205502
Toward Joint Inversion of Gravity and Dyanamics
NASA Astrophysics Data System (ADS)
Jacoby, W. R.
To better understand geodynamic processes as seafloor spreading, plumes, subduction, and isostatic adjustment, gravity is inverted with "a prioriinformation from topography/bathymetry, seismic structure and dynamic models. Examples are subduction of the Juan de Fuca plate below Vancouver Island, the passive Black SeaTurkey margin and Iceland ridge-plume interaction. Gravity and other data are averaged 50 km wide strips. Mass balances are estimated (showing also that the free air anomaly is misleading for narrow structures). The mass balances represent plate forces and plate bending, affecting the gravity signals and the isostatic state of continental margins and ridge-plume effects, which are highly correlated in space and cannot be separated without a priori information from modelling. The examples from widely different tectonic situations demonstrate that the art of regional-scale gravity inversion requires extensive background knowledge and inclusion of dynamic processes. It is difficult to conceive any formal, globally applicable procedure taking care of this; it is even a question, what is data, what a priori information? They are not distinguishable if all are included as foreward routines. The "accuracy" of models cannot be perfectly determined, if the "real" mass distribution is not known if known, gravity inversion would be unnecessary. In reality only guesses are possible on the basis of observations and physical laws governing geodynamics. A priori information and gravity data limit the resolution of gravity inversion. Different model types are indistinguishable because adjustments within their parameter uncertainties permit a good fit. But gravity excludes wrong models (Karl Popper: science evolves by falsification of wrong models), and precise gravity guides and defines aims, targets and strategies for new observations.
Inversion of electromagnetic data: An overview of new techniques
NASA Astrophysics Data System (ADS)
Oldenburg, Doug
1990-09-01
This paper explores some of the newer techniques for acquiring and inverting electromagnetic data. Attention is confined primarily to the 2d magnetotelluric (MT) problem but the inverse methods are applicable to all areas of EM induction. The basis of the EMAP technique of Bostick is presented along with examples to illustrate the efficacy of that method in structural imaging and in overcoming the deleterious effects of near-surface distortions of the electric field. Reflectivity imaging methods and the application of seismic migration techniques to EM problems are also explored as imaging tools. Two new approaches to the solution of the inverse problem are presented. The AIM (Approximate Inverse Mapping) inversion of Oldenburg and Ellis uses a new way to estimate a perturbation in an iterative solution which does not involve linearization of the equations. The RRI (Rapid Relaxation Inverse) of Smith and Booker shows how approximate Fréchet derivatives and sequences of 1d inversions can be used to develop a practical inversion algorithm. The overview is structured to provide insight about the latest inversion techniques and also to touch upon most areas of the inverse problem that must be considered to carry out a practical inversion. These include model parameterization, methods of calculating first order sensitivities, and methods for setting up a linearized inversion.
EDITORIAL: Inverse Problems' 25th year of publication Inverse Problems' 25th year of publication
NASA Astrophysics Data System (ADS)
2008-01-01
2009 is Inverse Problems' 25th year of publication. In this quarter-century, the journal has established itself as the premier publication venue for inverse problems research. It has matured from its beginnings as a niche journal serving the emerging field of inverse and ill-posed problems to a monthly publication in 2009 covering all aspects of a well-established, vibrant and still-expanding subject. Along with its core readership of pure and applied mathematicians and physicists, Inverse Problems has become widely known across a broad range of researchers in areas such as geophysics, optics, radar, acoustics, communication theory, signal processing and medical imaging, amongst others. The journal's appeal to the inverse problems community and those researchers from the varied fields that encounter such problems can be attributed to our commitment to publishing only the very best papers, and to offering unique services to the community. Besides our regular research papers, which average a remarkably short five months from submission to electronic publication, we regularly publish heavily cited topical review papers and topic-specific special sections, which first appeared in 2004. These highly-downloaded invited articles focus on the latest developments and hot topics in all areas of inverse problems. No other journal in the field offers these features. I am very pleased to take Inverse Problems into its 25th year as Editor-in-Chief. The journal has an impressive tradition of scholarship, established at its inception by the founder and first Editor-in-Chief, Professor Pierre Sabatier. Professor Sabatier envisioned the journal in 1985 as providing a medium for publication of exemplary research in our intrinsically interdisciplinary field. I am glad to say that the support of our authors, readers, referees, Editors-in-Chief, Editorial Boards and Advisory Panels over the years, has resulted in Inverse Problems becoming the top publication in this field, publishing
Frequency-domain direct waveform inversion based on perturbation theory
NASA Astrophysics Data System (ADS)
Kwak, Sangmin; Kim, Youngseo; Shin, Changsoo
2014-05-01
A direct waveform inversion based on perturbation theory is proposed to delineate a subsurface velocity structure from seismic data. This technique can directly compute the difference between the actual subsurface velocity and an initial guess of the velocity, while full waveform inversion updates the velocity model in the directions of reducing the data residual. Unlike full waveform inversion using the steepest descent method, the direct waveform inversion does not require a proper step length to iteratively update the velocity model. We present an algorithm for the waveform inversion method in the frequency domain and numerical examples demonstrating how the inversion method can reconstruct subsurface velocity structures using surface seismic data. The time-domain seismograms synthesized in the inversion procedure match the corresponding shot-gather seismograms of field data.
Postpartum Prolapsed Leiomyoma with Uterine Inversion Managed by Vaginal Hysterectomy
Pieh-Holder, Kelly L.; DeVente, James E.
2014-01-01
Background. Uterine inversion is a rare, but life threatening, obstetrical emergency which occurs when the uterine fundus collapses into the endometrial cavity. Various conservative and surgical therapies have been outlined in the literature for the management of uterine inversions. Case. We present a case of a chronic, recurrent uterine inversion, which was diagnosed following spontaneous vaginal delivery and recurred seven weeks later. The uterine inversion was likely due to a leiomyoma. This late-presenting, chronic, recurring uterine inversion was treated with a vaginal hysterectomy. Conclusion. Uterine inversions can occur in both acute and chronic phases. Persistent vaginal bleeding with the appearance of a prolapsing fibroid should prompt further investigation for uterine inversion and may require surgical therapy. A vaginal hysterectomy may be an appropriate management option in select populations and may be considered in women who do not desire to maintain reproductive function. PMID:25379314
Inverse spin Hall effect by spin injection
NASA Astrophysics Data System (ADS)
Liu, S. Y.; Horing, Norman J. M.; Lei, X. L.
2007-09-01
Motivated by a recent experiment [S. O. Valenzuela and M. Tinkham, Nature (London) 442, 176 (2006)], the authors present a quantitative microscopic theory to investigate the inverse spin-Hall effect with spin injection into aluminum considering both intrinsic and extrinsic spin-orbit couplings using the orthogonalized-plane-wave method. Their theoretical results are in good agreement with the experimental data. It is also clear that the magnitude of the anomalous Hall resistivity is mainly due to contributions from extrinsic skew scattering.
The inverse problem of the optimal regulator.
NASA Technical Reports Server (NTRS)
Yokoyama, R.; Kinnen, E.
1972-01-01
The inverse problem of the optimal regulator is considered for a general class of multi-input systems with integral-type performance indices. A new phase variable canonical form is shown to be convenient for this analysis. The advantage of the canonical form is to separate the state variables into subvectors of directly controlled, indirectly controlled, and uncontrollable components. Necessary and sufficient conditions for optimized performance indices are given. With the nonlinearities of the system restricted to functions of the directly controlled state variables, additional results are developed about the nonnegative property of optimized loss functions.
Tubular inverse opal scaffolds for biomimetic vessels
NASA Astrophysics Data System (ADS)
Zhao, Ze; Wang, Jie; Lu, Jie; Yu, Yunru; Fu, Fanfan; Wang, Huan; Liu, Yuxiao; Zhao, Yuanjin; Gu, Zhongze
2016-07-01
There is a clinical need for tissue-engineered blood vessels that can be used to replace or bypass damaged arteries. The success of such grafts depends strongly on their ability to mimic native arteries; however, currently available artificial vessels are restricted by their complex processing, controversial integrity, or uncontrollable cell location and orientation. Here, we present new tubular scaffolds with specific surface microstructures for structural vessel mimicry. The tubular scaffolds are fabricated by rotationally expanding three-dimensional tubular inverse opals that are replicated from colloidal crystal templates in capillaries. Because of the ordered porous structure of the inverse opals, the expanded tubular scaffolds are imparted with circumferentially oriented elliptical pattern microstructures on their surfaces. It is demonstrated that these tailored tubular scaffolds can effectively make endothelial cells to form an integrated hollow tubular structure on their inner surface and induce smooth muscle cells to form a circumferential orientation on their outer surface. These features of our tubular scaffolds make them highly promising for the construction of biomimetic blood vessels.There is a clinical need for tissue-engineered blood vessels that can be used to replace or bypass damaged arteries. The success of such grafts depends strongly on their ability to mimic native arteries; however, currently available artificial vessels are restricted by their complex processing, controversial integrity, or uncontrollable cell location and orientation. Here, we present new tubular scaffolds with specific surface microstructures for structural vessel mimicry. The tubular scaffolds are fabricated by rotationally expanding three-dimensional tubular inverse opals that are replicated from colloidal crystal templates in capillaries. Because of the ordered porous structure of the inverse opals, the expanded tubular scaffolds are imparted with circumferentially
Full Waveform Inversion with Optimal Basis Functions
NASA Astrophysics Data System (ADS)
Sun, Gang; Chang, Qianshun; Sheng, Ping
2003-03-01
Based on the approach suggested by Tarantola, and Gauthier etal., we show that the alternate use of the step (linear) function basis and the block function (quasi-δ function) basis can give accurate full waveform inversion results for the layered acoustic systems, starting from a uniform background. Our method is robust against additive white noise (up to 20% of the signal) and can resolve layers that are comparable to or smaller than a wavelength in thickness. The physical reason for the success of our approach is illustrated through a simple example.
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
Inverse Cerenkov laser acceleration experiment at ATF
NASA Astrophysics Data System (ADS)
Wang, X. J.; Pogorelsky, I.; Fernow, R.; Kusche, K. P.; Liu, Y.; Kimura, W. D.; Kim, G. H.; Romea, R. D.; Steinhauer, L. C.
Inverse Cerenkov laser acceleration was demonstrated using an axicon optical system at the Brookhaven Accelerator Test Facility (ATF). The ATF S-band linac and a high power 10.6 MICROMETERSCO2 laser were used for the experiment. Experimental arrangement and the laser and the electron beams synchronization are discussed. The electrons were accelerated more than 0.7 MeV for a 34 MW CO2 laser power. More than 3.7 MeV acceleration was measured with 0.7 GW CO2 laser power, which is more than 20 times of the previous ICA experiment. The experimental results are compared with computer program TRANSPORT simulations.
A vector inverse algorithm for electromagnetic scattering
NASA Astrophysics Data System (ADS)
Borden, B.
1984-06-01
Investigated is an inverse electromagnetic scattering technique that uses the polarization characteristics of the scattered wave to form an image of the convex portions of the scattering body. The depolarization of an electromagnetic signal by scattering surface is related to the local principal curvatures through the measurable leading edge of the impulse response. A classic problem in differential geometry (Christoffel-Hurwitz) deals with the reconstruction of such a surface from a knowledge of this kind of information, and a differential equation relating these local measurements to the surface has long been established. A Fortran code employing a 'finite-element' solution to this equation has been constructed and tested on synthetic data.
Fast Parallel Computation Of Manipulator Inverse Dynamics
NASA Technical Reports Server (NTRS)
Fijany, Amir; Bejczy, Antal K.
1991-01-01
Method for fast parallel computation of inverse dynamics problem, essential for real-time dynamic control and simulation of robot manipulators, undergoing development. Enables exploitation of high degree of parallelism and, achievement of significant computational efficiency, while minimizing various communication and synchronization overheads as well as complexity of required computer architecture. Universal real-time robotic controller and simulator (URRCS) consists of internal host processor and several SIMD processors with ring topology. Architecture modular and expandable: more SIMD processors added to match size of problem. Operate asynchronously and in MIMD fashion.
The NYU inverse swept wing code
NASA Technical Reports Server (NTRS)
Bauer, F.; Garabedian, P.; Mcfadden, G.
1983-01-01
An inverse swept wing code is described that is based on the widely used transonic flow program FLO22. The new code incorporates a free boundary algorithm permitting the pressure distribution to be prescribed over a portion of the wing surface. A special routine is included to calculate the wave drag, which can be minimized in its dependence on the pressure distribution. An alternate formulation of the boundary condition at infinity was introduced to enhance the speed and accuracy of the code. A FORTRAN listing of the code and a listing of a sample run are presented. There is also a user's manual as well as glossaries of input and output parameters.
Solving inversion problems with neural networks
NASA Technical Reports Server (NTRS)
Kamgar-Parsi, Behzad; Gualtieri, J. A.
1990-01-01
A class of inverse problems in remote sensing can be characterized by Q = F(x), where F is a nonlinear and noninvertible (or hard to invert) operator, and the objective is to infer the unknowns, x, from the observed quantities, Q. Since the number of observations is usually greater than the number of unknowns, these problems are formulated as optimization problems, which can be solved by a variety of techniques. The feasibility of neural networks for solving such problems is presently investigated. As an example, the problem of finding the atmospheric ozone profile from measured ultraviolet radiances is studied.
Inverse Cerenkov acceleration using an IFEL prebuncher
Kimura, W.D.; Pogorelsky, I.V.; Liu, Y.; Kusche, K.P.; van Steenbergen, A.; Gallardo, J.C.; Sandweiss, J.; Cline, D.B.; Quimby, D.C.; Babzien, M.
1997-03-01
The BNL IFEL will be used to optically prebunch the e-beam before sending it into an inverse Cerenkov acceleration (ICA) stage. Prebunching the beam will greatly improve the efficiency of the ICA process. The basic experimental design and preliminary model predictions for the combined ICA/IFEL experiment are discussed. Near-term goals are to demonstrate optical prebunching, rephasing of the prebunched beam with the optical field, and more efficient acceleration. Long-term goals are to demonstrate 100 MeV net acceleration using an ICA accelerator. {copyright} {ital 1997 American Institute of Physics.}
Inverse Faraday effect driven by radiation friction
NASA Astrophysics Data System (ADS)
Liseykina, T. V.; Popruzhenko, S. V.; Macchi, A.
2016-07-01
A collective, macroscopic signature to detect radiation friction in laser–plasma experiments is proposed. In the interaction of superintense circularly polarized laser pulses with high density targets, the effective dissipation due to radiative losses allows the absorption of electromagnetic angular momentum, which in turn leads to the generation of a quasistatic axial magnetic field. This peculiar ‘inverse Faraday effect’ is investigated by analytical modeling and three-dimensional simulations, showing that multi-gigagauss magnetic fields may be generated at laser intensities \\gt {10}23 {{{W}}{{cm}}}-2.
Inversion layer solar cell fabrication and evaluation
NASA Technical Reports Server (NTRS)
Call, R. L.
1974-01-01
Inversion layer solar cells have been fabricated by etching through the diffused layer on p-type silicon wafers in a comb-like contact pattern. The charge separation comes from an induced p-n junction at the surface. This inverted surface is caused by a layer of transparent material applied to the surface that either contains free positive ions or that creates donor states at the interface. Cells have increased from 3 ma Isc to 100 ma by application of sodium silicate. The action is unstable, however, and decays with time.
A climatology of tropospheric humidity inversions in five reanalyses
NASA Astrophysics Data System (ADS)
Brunke, Michael A.; Stegall, Steve T.; Zeng, Xubin
2015-02-01
Specific humidity is generally thought to decrease with height in the troposphere. However, here we document the existence of specific humidity inversions in five reanalyses: the National Centers for Environmental Prediction (NCEP) second reanalysis (NCEP-2), the European Centre for Medium-Range Forecasts (ECMWF) 40-year reanalysis (ERA-40), the Modern Era Retrospective Analysis for Research Applications (MERRA), NCEP's Climate Forecast System Reanalysis (CFSR), and the ECMWF interim reanalysis (ERA-Interim). These inversions are most frequent in the polar regions. Inversions do occur elsewhere, most notably over the subtropical stratus regions, but are less frequent and likely overproduced depending on the location. Polar inversions are the most persistent in winter and the strongest (as defined by the humidity difference divided by the pressure difference across the inversion) in summer or autumn with low bases (at pressures > 900 hPa). Winter humidity inversions are lower, being near-surface, due to the persistence of low-level temperature inversions associated with these humidity inversions, while summer humidity inversions tend to be located near cloud top providing moisture to prevent the melt season stratus from evaporating. The most important contributions to affect humidity inversions in MERRA are dynamics, turbulence, and moist physics. However, local advection may not play as much of a role as regional humidity convergence. The subtropical stratus inversions are as thick as polar humidity inversions but with higher bases generally at pressures < 900 hPa. These inversions are confirmed by rawinsonde data, but there are discrepancies between the observed annual and diurnal cycles in inversion frequency and those portrayed in the reanalyses.
Earthquake source inversion of tsunami runup prediction
NASA Astrophysics Data System (ADS)
Sekar, Anusha
Our goal is to study two inverse problems: using seismic data to invert for earthquake parameters and using tide gauge data to invert for earthquake parameters. We focus on the feasibility of using a combination of these inverse problems to improve tsunami runup prediction. A considerable part of the thesis is devoted to studying the seismic forward operator and its modeling using immersed interface methods. We develop an immersed interface method for solving the variable coefficient advection equation in one dimension with a propagating singularity and prove a convergence result for this method. We also prove a convergence result for the one-dimensional acoustic system of partial differential equations solved using immersed interface methods with internal boundary conditions. Such systems form the building blocks of the numerical model for the earthquake. For a simple earthquake-tsunami model, we observe a variety of possibilities in the recovery of the earthquake parameters and tsunami runup prediction. In some cases the data are insufficient either to invert for the earthquake parameters or to predict the runup. When more data are added, we are able to resolve the earthquake parameters with enough accuracy to predict the runup. We expect that this variety will be true in a real world three dimensional geometry as well.
Approaching the Island of Inversion: 34P
Bender, P.C.; Hoffman, C.R.; Wiedeking, M.; Allmond, J.M.; Bernstein, L.A.; Burke, J.T.; Bleuel, D.L.; Clark, R.M.; Fallon, P.; Goldblum, B.L.; Hinners, T.A.; Jeppesen, H.B.; Lee, Sangjin; Lee, I.Y.; Lesher, S.R.; Machiavelli, A.O.; McMahan, M.A.; Morris, D.; Perry, M.; Phair, L.; Scielzo, N.D.; Tabor, S.L.; Tripathi, Vandana; Volya, A.
2011-06-14
Yrast states in 34P were investigated using the 18O(18O,pn) reaction at energies of 20, 24, 25, 30, and 44 MeV at Florida State University and at Lawrence Berkeley National Laboratory. The level scheme was expanded, ray angular distributions were measured, and lifetimes were inferred with the Doppler-shift attenuation method by detecting decay protons in coincidence with one or more rays. The results provide a clearer picture of the evolution of structure approaching the 'Island of Inversion', particularly how the 1 and 2 particle-hole (ph) states fall in energy with increasing neutro number approaching inversion. However, the agreement of the lowest few states with pure sd shell model predictions shows that the level scheme of 34P is not itself inverted. Rather, the accumulated evidence indicates that the 1-ph states start at 2.3 MeV. A good candidate for the lowest 2-ph state lies at 6236 keV, just below the neutron separation energy of 6291 keV. Shell model calculations made using a small modification of the WBP interaction reproduce the negative-parity, 1-ph states rather well.
Spontaneous sense inversion in helical mesophases
NASA Astrophysics Data System (ADS)
Wensink, H. H.
2014-08-01
We investigate the pitch sensitivity of chiral nematic phases of helicoidal patchy cylinders as a generic model for chiral (bio-)polymers and helix-shaped colloidal rods. The behaviour of the macroscopic helical pitch is studied from microscopic principles by invoking a simple density functional theory generalised to accommodate weakly twisted director fields. Upon changing the degree of alignment along the local helicoidal director we find that chiral nematic phases exhibit a sudden sense inversion whereby the helical symmetry changes from left- to right-handed and vice versa. Since the local alignment is governed by thermodynamic variables such as density, temperature or the amplitude of an external directional field, such pitch sense inversions can be expected in systems of helical mesogens of both thermotropic and lyotropic origin. We show that the spontaneous change of helical symmetry is a direct consequence of an antagonistic effective torque between helical particles with a certain prescribed internal helicity. The results may help opening up new routes towards precise control of the helical handedness of chiral assemblies by a judicious choice of external control parameters.