Measurement of adjoint flux at the RB reactor
The adjoint flux is of the great importance for determination of kinetic parameters of nuclear reactor (ρ, l and βeff) and for the interpretation of experiments with reactivity perturbations. In experimental reactor physics there are a few methods for the adjoint flux measurements. The method of reactivity perturbations with adequate samples is used for thermal reactors. According to the theory of reactivity perturbations the reactivity change due to sample of thermal neutrons absorbing material is proportional to product of flux and adjoint flux of thermal neutrons (Φ2(r)Φ=2(r)). The reactivity change due to fissionable nuclide is proportional to product of thermal neutron flux and adjoint flux of fast neutrons (Φ2(r)Φ=2(r)). The axial distribution of adjoint flux is determined by reactivity measurements and measurements of axial distribution of thermal neutron flux. Thi results of this measurement will be used for interpretation of other experiments with reactivity perturbations at the RB reactor
Khuat, Quang Huy; Kim, Song Hyun; Kim, Do Hyun; Shin, Chang Ho [Hanyang University, Seoul (Korea, Republic of)
2015-05-15
This technique is known as Consistent Adjoint Driven Importance Sampling (CADIS) method and it is implemented in SCALE code system. In the CADIS method, adjoint transport equation has to be solved to determine deterministic importance functions. Using the CADIS method, a problem was noted that the biased adjoint flux estimated by deterministic methods can affect the calculation efficiency and error. The biases of adjoint function are caused by the methodology, calculation strategy, tolerance of result calculated by the deterministic method and inaccurate multi-group cross section libraries. In this paper, a study to analyze the influence of the biased adjoint functions into Monte Carlo computational efficiency is pursued. In this study, a method to estimate the calculation efficiency was proposed for applying the biased adjoint fluxes in the CADIS approach. For a benchmark problem, the responses and FOMs using SCALE code system were evaluated as applying the adjoint fluxes. The results show that the biased adjoint fluxes significantly affects the calculation efficiencies.
Adjoint Monte Carlo calculation of charged plasma particle flux to wall
Äkäslompolo, Simppa
2015-01-01
This manuscript describes an adjoint/reverse Monte Carlo method to calculate the flux of charged plasma particles to the wall of e.g. a tokamak. Two applications are described: a fusion product activation probe and a neutral beam injection prompt loss measurement with a fast ion loss diagnostic. In both cases, the collisions of the particles with the background plasma can be omitted.
Humbird, Kelli D
2016-01-01
Uncertainty quantification and sensitivity analyses are a vital component for predictive modeling in the sciences and engineering. The adjoint approach to sensitivity analysis requires solving a primary system of equations and a mathematically related set of adjoint equations. The information contained in the equations can be combined to produce sensitivity information in a computationally efficient manner. In this work, sensitivity analyses are performed on systems described by flux-limited radiative diffusion using the adjoint approach. The sensitivities computed are shown to agree with standard perturbation theory, and can be obtained in significantly less computational time.
It is noted that the analog Monte Carlo method has low calculation efficiency at deep penetration problems such as radiation shielding analysis. In order to increase the calculation efficiency, variance reduction techniques have been introduced and applied for the shielding calculation. To optimize the variance reduction technique, the hybrid Monte Carlo method was introduced. For the determination of the parameters using the hybrid Monte Carlo method, the adjoint flux should be calculated by the deterministic methods. In this study, the collision probability method is applied to calculate adjoint flux. The solution of integration transport equation in the collision probability method is modified to calculate the adjoint flux approximately even for complex and arbitrary geometries. For the calculation, C++ program is developed. By using the calculated adjoint flux, importance parameters of each cell in shielding material are determined and used for variance reduction of transport calculation. In order to evaluate calculation efficiency with the proposed method, shielding calculations are performed with MCNPX 2.7. In this study, a method to calculate the adjoint flux in using the Monte Carlo variance reduction was proposed to improve Monte Carlo calculation efficiency of thick shielding problem. The importance parameter for each cell of shielding material is determined by calculating adjoint flux with the modified collision probability method. In order to calculate adjoint flux with the proposed method, C++ program is developed. The results show that the proposed method can efficiently increase the FOM of transport calculation. It is expected that the proposed method can be utilize for the calculation efficiency in thick shielding calculation
Methods of Monte Carlo biasing using two-dimensional discrete ordinates adjoint flux
Tang, J.S.; Stevens, P.N.; Hoffman, T.J.
1976-06-01
Methods of biasing three-dimensional deep penetration Monte Carlo calculations using importance functions obtained from a two-dimensional discrete ordinates adjoint calculation have been developed. The important distinction was made between the applications of the point value and the event value to alter the random walk in Monte Carlo analysis of radiation transport. The biasing techniques developed are the angular probability biasing which alters the collision kernel using the point value as the importance function and the path length biasing which alters the transport kernel using the event value as the importance function. Source location biasings using the step importance function and the scalar adjoint flux obtained from the two-dimensional discrete ordinates adjoint calculation were also investigated. The effects of the biasing techniques to Monte Carlo calculations have been investigated for neutron transport through a thick concrete shield with a penetrating duct. Source location biasing, angular probability biasing, and path length biasing were employed individually and in various combinations. Results of the biased Monte Carlo calculations were compared with the standard Monte Carlo and discrete ordinates calculations.
Stück, Arthur
2015-11-01
Inconsistent discrete expressions in the boundary treatment of Navier-Stokes solvers and in the definition of force objective functionals can lead to discrete-adjoint boundary treatments that are not a valid representation of the boundary conditions to the corresponding adjoint partial differential equations. The underlying problem is studied for an elementary 1D advection-diffusion problem first using a node-centred finite-volume discretisation. The defect of the boundary operators in the inconsistently defined discrete-adjoint problem leads to oscillations and becomes evident with the additional insight of the continuous-adjoint approach. A homogenisation of the discretisations for the primal boundary treatment and the force objective functional yields second-order functional accuracy and eliminates the defect in the discrete-adjoint boundary treatment. Subsequently, the issue is studied for aerodynamic Reynolds-averaged Navier-Stokes problems in conjunction with a standard finite-volume discretisation on median-dual grids and a strong implementation of noslip walls, found in many unstructured general-purpose flow solvers. Going out from a base-line discretisation of force objective functionals which is independent of the boundary treatment in the flow solver, two improved flux-consistent schemes are presented; based on either body wall-defined or farfield-defined control-volumes they resolve the dual inconsistency. The behaviour of the schemes is investigated on a sequence of grids in 2D and 3D.
Laboure, Vincent M; Wang, Yaqi
2016-01-01
In this paper, we derive a method for the second-order form of the transport equation that is both globally conservative and compatible with voids, using Continuous Finite Element Methods (CFEM). The main idea is to use the Least-Squares (LS) form of the transport equation in the void regions and the Self-Adjoint Angular Flux (SAAF) form elsewhere. While the SAAF formulation is globally conservative, the LS formulation need a correction in void. The price to pay for this fix is the loss of symmetry of the bilinear form. We first derive this Conservative LS (CLS) formulation in void. Second we combine the SAAF and CLS forms and end up with an hybrid SAAF-CLS method, having the desired properties. We show that extending the theory to near-void regions is a minor complication and can be done without affecting the global conservation of the scheme. Being angular discretization agnostic, this method can be applied to both discrete ordinates (SN) and spherical harmonics (PN) methods. However, since a globally conse...
Nonlinear acceleration of a continuous finite element (CFE) discretization of the transport equation requires a modification of the transport solution in order to achieve local conservation, a condition used in nonlinear acceleration to define the stopping criterion. In this work we implement a coarse-mesh finite difference acceleration for a CFE discretization of the second-order self-adjoint angular flux (SAAF) form of the transport equation and use a postprocessing to enforce local conservation. Numerical results are given for one-group source calculations of one-dimensional slabs. We also give a novel formal derivation of the boundary conditions for the SAAF. (authors)
Truchet, G.; Leconte, P.; Peneliau, Y.; Santamarina, A.; Malvagi, F.
2014-06-01
Pile-oscillation experiments are performed in the MINERVE reactor at the CEA Cadarache to improve nuclear data accuracy. In order to precisely calculate small reactivity variations (TRIPOLI-4® by using the eigenvalue difference method. This "direct" method has shown limitations in the evaluation of very small reactivity effects because it needs to reach a very small variance associated to the reactivity in both states. To answer this problem, it has been decided to implement the exact perturbation theory in TRIPOLI-4® and, consequently, to calculate a continuous-energy adjoint flux. The Iterated Fission Probability (IFP) method was chosen because it has shown great results in some other Monte Carlo codes. The IFP method uses a forward calculation to compute the adjoint flux, and consequently, it does not rely on complex code modifications but on the physical definition of the adjoint flux as a phase-space neutron importance. In the first part of this paper, the IFP method implemented in TRIPOLI-4® is described. To illustrate the effciency of the method, several adjoint fluxes are calculated and compared with their equivalent obtained by the deterministic code APOLLO-2. The new implementation can calculate angular adjoint flux. In the second part, a procedure to carry out an exact perturbation calculation is described. A single cell benchmark has been used to test the accuracy of the method, compared with the "direct" estimation of the perturbation. Once again the method based on the IFP shows good agreement for a calculation time far more inferior to the "direct" method. The main advantage of the method is that the relative accuracy of the reactivity variation does not depend on the magnitude of the variation itself, which allows us to calculate very small reactivity perturbations with high precision. Other applications of this perturbation method are presented and tested like the calculation of exact kinetic parameters (βeff, Λeff) or sensitivity parameters.
A novel approach is proposed for charged particle transport calculations using a recently developed second-order, self-adjoint angular flux (SAAF) form of the Boltzmann transport equation with continuous slowing-down. A finite element discretization that is linear continuous in space and linear discontinuous (LD) in energy is described and implemented in a one-dimensional, planar geometry, multigroup, discrete ordinates code for charged particle transport. The cross-section generating code CEPXS is used to generate the electron and photon transport cross sections employed in this code. The discrete ordinates SAAF transport equation is solved using source iteration in conjunction with an inner iteration acceleration scheme and an outer iteration acceleration scheme. Outer iterations are required with the LD energy discretization scheme because the two angular flux unknowns within each group are coupled, which gives rise to effective upscattering. The inner iteration convergence is accelerated using diffusion synthetic acceleration, and the outer iteration convergence is accelerated using a diamond difference approximation to the LD energy discretization. Computational results are given that demonstrate the effectiveness of our convergence acceleration schemes and the accuracy of our discretized SAAF equation
Metabolic Flux Analysis in Isotope Labeling Experiments using the Adjoint Approach
Mottelet, Stéphane; Gaullier, Gil; Sadaka, Georges
2016-01-01
Comprehension of metabolic pathways is considerably enhanced by metabolic flux analysis (MFA-ILE) in isotope labeling experiments. The balance equations are given by hundreds of algebraic (stationary MFA) or ordinary differential equations (nonstationary MFA), and reducing the number of operations is therefore a crucial part of reducing the computation cost. The main bottleneck for deterministic algorithms is the computation of derivatives, particularly for nonstationary MFA. In this article ...
System of adjoint P1 equations for neutron moderation
In some applications of perturbation theory, it is necessary know the adjoint neutron flux, which is obtained by the solution of adjoint neutron diffusion equation. However, the multigroup constants used for this are weighted in only the direct neutron flux, from the solution of direct P1 equations. In this work, this procedure is questioned and the adjoint P1 equations are derived by the neutron transport equation, the reversion operators rules and analogies between direct and adjoint parameters. (author)
Adjoint P1 equations solution for neutron slowing down
In some applications of perturbation theory, it is necessary know the adjoint neutron flux, which is obtained by the solution of adjoint neutron diffusion equation. However, the multigroup constants used for this are weighted in only the direct neutron flux, from the solution of direct P1 equations. In this work, the adjoint P1 equations are derived by the neutron transport equation, the reversion operators rules and analogies between direct and adjoint parameters. The direct and adjoint neutron fluxes resulting from the solution of P1 equations were used to three different weighting processes, to obtain the macrogroup macroscopic cross sections. It was found out noticeable differences among them. (author)
Fereydooni, Abolhassan; Safapour, Ahmad; Rahimi , Asghar
2012-01-01
The concept of (p,q)-pair frames is generalized to (l,l^*)-pair frames. Adjoint (conjugate) of a pair frames for dual space of a Banach space is introduced and some conditions for the existence of adjoint (conjugate) of pair frames are presented.
Martinez, Aquilino Senra; Silva, Fernando Carvalho da; Cardoso, Carlos Eduardo Santos [Universidade Federal, Rio de Janeiro, RJ (Brazil). Coordenacao dos Programas de Pos-graduacao de Engenharia. Programa de Engenharia Nuclear
2000-07-01
In some applications of perturbation theory, it is necessary know the adjoint neutron flux, which is obtained by the solution of adjoint neutron diffusion equation. However, the multigroup constants used for this are weighted in only the direct neutron flux, from the solution of direct P1 equations. In this work, this procedure is questioned and the adjoint P1 equations are derived by the neutron transport equation, the reversion operators rules and analogies between direct and adjoint parameters. (author)
Adjoint functors in graph theory
Foniok, Jan; Tardif, Claude
2013-01-01
We survey some uses of adjoint functors in graph theory pertaining to colourings, complexity reductions, multiplicativity, circular colourings and tree duality. The exposition of these applications through adjoint functors unifies the presentation to some extent, and also raises interesting questions.
Adjoint affine fusion and tadpoles
Urichuk, Andrew; Walton, Mark A.
2016-06-01
We study affine fusion with the adjoint representation. For simple Lie algebras, elementary and universal formulas determine the decomposition of a tensor product of an integrable highest-weight representation with the adjoint representation. Using the (refined) affine depth rule, we prove that equally striking results apply to adjoint affine fusion. For diagonal fusion, a coefficient equals the number of nonzero Dynkin labels of the relevant affine highest weight, minus 1. A nice lattice-polytope interpretation follows and allows the straightforward calculation of the genus-1 1-point adjoint Verlinde dimension, the adjoint affine fusion tadpole. Explicit formulas, (piecewise) polynomial in the level, are written for the adjoint tadpoles of all classical Lie algebras. We show that off-diagonal adjoint affine fusion is obtained from the corresponding tensor product by simply dropping non-dominant representations.
Adjoint affine fusion and tadpoles
Urichuk, Andrew
2016-01-01
We study affine fusion with the adjoint representation. For simple Lie algebras, elementary and universal formulas determine the decomposition of a tensor product of an integrable highest-weight representation with the adjoint representation. Using the (refined) affine depth rule, we prove that equally striking results apply to adjoint affine fusion. For diagonal fusion, a coefficient equals the number of nonzero Dynkin labels of the relevant affine highest weight, minus 1. A nice lattice-polytope interpretation follows, and allows the straightforward calculation of the genus-1 1-point adjoint Verlinde dimension, the adjoint affine fusion tadpole. Explicit formulas, (piecewise) polynomial in the level, are written for the adjoint tadpoles of all classical Lie algebras. We show that off-diagonal adjoint affine fusion is obtained from the corresponding tensor product by simply dropping non-dominant representations.
Sensitivity analysis via reduced order adjoint method
Notwithstanding the voluminous literature on adjoint sensitivity analysis, it has been generally dismissed by practitioners as cumbersome with limited value in realistic engineering models. This perception reflects two limitations about adjoint sensitivity analysis: a) its most effective application is limited to calculation of first-order variations; when higher order derivatives are required, it quickly becomes computationally inefficient; and b) the number of adjoint model evaluations depends on the number of responses, which renders it ineffective for multi-physics model where entire distributions, such as flux and power distribution, are often transferred between the various physics models. To overcome these challenges, this manuscript employs recent advances in reduced order modeling to re-cast the adjoint model equations into a form that renders its application to real reactor models practical. Past work applied reduced order modeling techniques to render reduction for general nonlinear high dimensional models by identifying mathematical subspaces, called active subspaces, that capture all dominant features of the model, including both linear and nonlinear variations. We demonstrate the application of these techniques to the calculation of first-order derivatives, or as commonly known sensitivity coefficients, for a fuel assembly model with many responses. We show that the computational cost becomes dependent on the physics model itself, via the so-called rank of the active subspace, rather than the number of responses or parameters. (author)
CHENG Qiang; CAO JianWen; WANG Bin; ZHANG HaiBin
2009-01-01
The adjoint code generator (ADG) is developed to produce the adjoint codes, which are used to analytically calculate gradients and the Hessian-vector products with the costs independent of the number of the independent variables. Different from other automatic differentiation tools, the implementation of ADG has advantages of using the least program behavior decomposition method and several static dependence analysis techniques. In this paper we first address the concerned concepts and fundamentals, and then introduce the functionality and the features of ADG. In particular, we also discuss the design architecture of ADG and implementation details including the recomputation and storing strategy and several techniques for code optimization. Some experimental results in several applications are presented at the end.
Adjoint affine fusion and tadpoles
Urichuk, Andrew; Walton, Mark A.
2016-01-01
We study affine fusion with the adjoint representation. For simple Lie algebras, elementary and universal formulas determine the decomposition of a tensor product of an integrable highest-weight representation with the adjoint representation. Using the (refined) affine depth rule, we prove that equally striking results apply to adjoint affine fusion. For diagonal fusion, a coefficient equals the number of nonzero Dynkin labels of the relevant affine highest weight, minus 1. A nice lattice-pol...
Hussein, A.; Krejčiřík, David; Siegl, P.
2015-01-01
Roč. 367, č. 4 (2015), s. 2921-2957. ISSN 0002-9947 R&D Projects: GA ČR GAP203/11/0701 Institutional support: RVO:61389005 Keywords : Laplacians on metric graphs * non-self-adjoint boundary conditions * similarity transforms to self-adjoint operators * Riesz basis Subject RIV: BE - Theoretical Physics Impact factor: 1.122, year: 2014
Implementation of Generalized Adjoint Equation Solver for DeCART
In this paper, the generalized adjoint solver based on the generalized perturbation theory is implemented on DeCART and the verification calculations were carried out. As the results, the adjoint flux for the general response coincides with the reference solution and it is expected that the solver could produce the parameters for the sensitivity and uncertainty analysis. Recently, MUSAD (Modules of Uncertainty and Sensitivity Analysis for DeCART) was developed for the uncertainty analysis of PMR200 core and the fundamental adjoint solver was implemented into DeCART. However, the application of the code was limited to the uncertainty to the multiplication factor, keff, because it was based on the classical perturbation theory. For the uncertainty analysis to the general response as like the power density, it is necessary to develop the analysis module based on the generalized perturbation theory and it needs the generalized adjoint solutions from DeCART. In this paper, the generalized adjoint solver is implemented on DeCART and the calculation results are compared with the results by TSUNAMI of SCALE 6.1
Implementation of Generalized Adjoint Equation Solver for DeCART
Han, Tae Young; Cho, Jin Young; Lee, Hyun Chul; Noh, Jae Man [Korea Atomic Energy Research Institute, Daejeon (Korea, Republic of)
2013-10-15
In this paper, the generalized adjoint solver based on the generalized perturbation theory is implemented on DeCART and the verification calculations were carried out. As the results, the adjoint flux for the general response coincides with the reference solution and it is expected that the solver could produce the parameters for the sensitivity and uncertainty analysis. Recently, MUSAD (Modules of Uncertainty and Sensitivity Analysis for DeCART) was developed for the uncertainty analysis of PMR200 core and the fundamental adjoint solver was implemented into DeCART. However, the application of the code was limited to the uncertainty to the multiplication factor, k{sub eff}, because it was based on the classical perturbation theory. For the uncertainty analysis to the general response as like the power density, it is necessary to develop the analysis module based on the generalized perturbation theory and it needs the generalized adjoint solutions from DeCART. In this paper, the generalized adjoint solver is implemented on DeCART and the calculation results are compared with the results by TSUNAMI of SCALE 6.1.
Chen, Han
2016-01-01
Many control applications can be formulated as optimization constrained by conservation laws. Such optimization can be efficiently solved by gradient-based methods, where the gradient is obtained through the adjoint method. Traditionally, the adjoint method has not been able to be implemented in "gray-box" conservation law simulations. In gray-box simulations, the analytical and numerical form of the conservation law is unknown, but the space-time solution of relevant flow quantities is available. Without the adjoint gradient, optimization can be challenging for problems with many control variables. However, much information about the gray-box simulation is contained in its space-time solution, which motivates us to estimate the adjoint gradient by leveraging the space-time solution. This article considers a type of gray-box simulations where the flux function is partially unknown. A method is introduced to estimate the adjoint gradient at a cost independent of the number of control variables. The method firs...
Estimation of Adjoint-Weighted Kinetics Parameters in Monte Carlo Wieland Calculations
Choi, Sung Hoon; Shim, Hyung Jin [Seoul National Univ., Seoul (Korea, Republic of)
2013-07-01
The effective delayed neutron fraction, β{sub eff}, and the prompt neutron generation time, Λ, in the point kinetics equation are weighted by the adjoint flux to improve the accuracy of the reactivity estimate. Recently the Monte Carlo (MC) kinetics parameter estimation methods by using the self-consistent adjoint flux calculated in the MC forward simulations have been developed and successfully applied for the research reactor analyses. However these adjoint estimation methods based on the cycle-by-cycle genealogical table require a huge memory size to store the pedigree hierarchy. In this paper, we present a new adjoint estimation in which the pedigree of a single history is utilized by applying the MC Wielandt method. The effectiveness of the new method is demonstrated in the kinetics parameter estimations for infinite homogeneous two-group problems and the Godiva critical facility.
The effective delayed neutron fraction, βeff, and the prompt neutron generation time, Λ, in the point kinetics equation are weighted by the adjoint flux to improve the accuracy of the reactivity estimate. Recently the Monte Carlo (MC) kinetics parameter estimation methods by using the adjoint flux calculated in the MC forward simulations have been developed and successfully applied for reactor analyses. However these adjoint estimation methods based on the cycle-by-cycle genealogical table require a huge memory size to store the pedigree hierarchy. In this paper, we present a new adjoint estimation method in which the pedigree of a single history is utilized by applying the MC Wielandt method. The algorithm of the new method is derived and its effectiveness is demonstrated in the kinetics parameter estimations for infinite homogeneous two-group problems and critical facilities. (author)
Estimation of Adjoint-Weighted Kinetics Parameters in Monte Carlo Wieland Calculations
The effective delayed neutron fraction, βeff, and the prompt neutron generation time, Λ, in the point kinetics equation are weighted by the adjoint flux to improve the accuracy of the reactivity estimate. Recently the Monte Carlo (MC) kinetics parameter estimation methods by using the self-consistent adjoint flux calculated in the MC forward simulations have been developed and successfully applied for the research reactor analyses. However these adjoint estimation methods based on the cycle-by-cycle genealogical table require a huge memory size to store the pedigree hierarchy. In this paper, we present a new adjoint estimation in which the pedigree of a single history is utilized by applying the MC Wielandt method. The effectiveness of the new method is demonstrated in the kinetics parameter estimations for infinite homogeneous two-group problems and the Godiva critical facility
Self-adjointness and polarization of the fermionic vacuum in the background of nontrivial topology
Singular configuration of an external static magnetic field in the form of a string polarizes vacuum in the secondly quantized theory on a plane which is orthogonal to the string axis. We consider the most general boundary conditions at the punctured singular point, which are compatible with the self-adjointness of the two-dimensional Dirac Hamiltonian. The dependence of the induced vacuum quantum numbers on the self-adjoint extension parameter and the flux of the string is determined
Adjoint Functors and Representation Dimensions
Chang Chang XI
2006-01-01
We study the global dimensions of the coherent functors over two categories that are linked by a pair of adjoint functors. This idea is then exploited to compare the representation dimensions of two algebras. In particular, we show that if an Artin algebra is switched from the other, then they have the same representation dimension.
Adjoint Weighted Kinetics Parameter Estimation in the Monte Carlo Wielandt Calculations
Choi, Sung Hoon; Shim, Hyung Jin [Seoul National Univ., Seoul (Korea, Republic of)
2013-10-15
In order to eliminate this huge memory consumption in the current adjoint estimation method, we have developed a new method in which the pedigree of a single history is utilized by applying the MC Wielandt method. The Wielandt method allows the estimations of the adjoint flux and adjoint-weighted parameters within a single cycle neutron simulations. The effectiveness of the new method is demonstrated in the kinetics parameter estimations for an infinite homogeneous two-group problem and the Godiva facility. We have developed an efficient algorithm for the adjoint-weighted kinetics parameter estimation in the MC Wielandt calculations which can significantly reduce the memory usage. From the numerical applications, it is demonstrated that the new method can predict the kinetics parameters with great accuracy.
Self-adjointness of deformed unbounded operators
Much, Albert [Instituto de Ciencias Nucleares, UNAM, México D.F. 04510 (Mexico)
2015-09-15
We consider deformations of unbounded operators by using the novel construction tool of warped convolutions. By using the Kato-Rellich theorem, we show that unbounded self-adjoint deformed operators are self-adjoint if they satisfy a certain condition. This condition proves itself to be necessary for the oscillatory integral to be well-defined. Moreover, different proofs are given for self-adjointness of deformed unbounded operators in the context of quantum mechanics and quantum field theory.
Self-Adjoint Extension of Symmetric Maps
Friedel, H. N.
2011-01-01
A densely-defined symmetric linear map from/to a real Hilbert space extends to a self-adjoint map. Extension is expressed via Riesz representation. For a case including Friedrichs extension of a strongly monotone map, self-adjoint extension is unique, and equals closure of the given map.
Curbelo, Jesus P.; Alves Filho, Hermes; Barros, Ricardo C., E-mail: jperez@iprj.uerj.br, E-mail: halves@iprj.uerj.br, E-mail: rcbarros@pq.cnpq.br [Universidade do Estado do Rio de Janeiro (UERJ), Nova Friburgo, RJ (Brazil). Instituto Politecnico. Programa de Pos-Graduacao em Modelagem Computacional; Hernandez, Carlos R.G., E-mail: cgh@instec.cu [Instituto Superior de Tecnologias y Ciencias Aplicadas (InSTEC), La Habana (Cuba)
2015-07-01
The spectral Green's function (SGF) method is a numerical method that is free of spatial truncation errors for slab-geometry fixed-source discrete ordinates (S{sub N}) adjoint problems. The method is based on the standard spatially discretized adjoint S{sub N} balance equations and a nonstandard adjoint auxiliary equation expressing the node-average adjoint angular flux, in each discretization node, as a weighted combination of the node-edge outgoing adjoint fluxes. The auxiliary equation contains parameters which act as Green's functions for the cell-average adjoint angular flux. These parameters are determined by means of a spectral analysis which yields the local general solution of the S{sub N} equations within each node of the discretization grid. In this work a number of advances in the SGF adjoint method are presented: the method is extended to adjoint S{sub N} problems considering linearly anisotropic scattering and non-zero prescribed boundary conditions for the forward source-detector problem. Numerical results to typical model problems are considered to illustrate the efficiency and accuracy of the o offered method. (author)
The spectral Green's function (SGF) method is a numerical method that is free of spatial truncation errors for slab-geometry fixed-source discrete ordinates (SN) adjoint problems. The method is based on the standard spatially discretized adjoint SN balance equations and a nonstandard adjoint auxiliary equation expressing the node-average adjoint angular flux, in each discretization node, as a weighted combination of the node-edge outgoing adjoint fluxes. The auxiliary equation contains parameters which act as Green's functions for the cell-average adjoint angular flux. These parameters are determined by means of a spectral analysis which yields the local general solution of the SN equations within each node of the discretization grid. In this work a number of advances in the SGF adjoint method are presented: the method is extended to adjoint SN problems considering linearly anisotropic scattering and non-zero prescribed boundary conditions for the forward source-detector problem. Numerical results to typical model problems are considered to illustrate the efficiency and accuracy of the o offered method. (author)
Southern California Adjoint Source Inversions
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.
Experience with Monte Carlo variance reduction using adjoint solutions in HYPER neutronics analysis
The variance reduction techniques using adjoint solutions are applied to the Monte Carlo calculation of the HYPER(HYbrid Power Extraction Reactor) core neutronics. The applied variance reduction techniques are the geometry splitting and the weight windows. The weight bounds and the cell importance needed for these techniques are generated from an adjoint discrete ordinate calculation by the two-dimensional TWODANT code. The flux distribution variances of the Monte Carlo calculations by these variance reduction techniques are compared with the results of the standard Monte Carlo calculations. It is shown that the variance reduction techniques using adjoint solutions to the HYPER core neutronics result in a decrease in the efficiency of the Monte Carlo calculation
Generalized adjoint consistent treatment of wall boundary conditions for compressible flows
Hartmann, Ralf; Leicht, Tobias
2015-11-01
In this article, we revisit the adjoint consistency analysis of Discontinuous Galerkin discretizations of the compressible Euler and Navier-Stokes equations with application to the Reynolds-averaged Navier-Stokes and k- ω turbulence equations. Here, particular emphasis is laid on the discretization of wall boundary conditions. While previously only one specific combination of discretizations of wall boundary conditions and of aerodynamic force coefficients has been shown to give an adjoint consistent discretization, in this article we generalize this analysis and provide a discretization of the force coefficients for any consistent discretization of wall boundary conditions. Furthermore, we demonstrate that a related evaluation of the cp- and cf-distributions is required. The freedom gained in choosing the discretization of boundary conditions without loosing adjoint consistency is used to devise a new adjoint consistent discretization including numerical fluxes on the wall boundary which is more robust than the adjoint consistent discretization known up to now. While this work is presented in the framework of Discontinuous Galerkin discretizations, the insight gained is also applicable to (and thus valuable for) other discretization schemes. In particular, the discretization of integral quantities, like the drag, lift and moment coefficients, as well as the discretization of local quantities at the wall like surface pressure and skin friction should follow as closely as possible the discretization of the flow equations and boundary conditions at the wall boundary.
Double-difference adjoint seismic tomography
Yuan, Yanhua O.; Simons, Frederik J.; Tromp, Jeroen
2016-06-01
We introduce a `double-difference' method for the inversion for seismic wavespeed structure based on adjoint tomography. Differences between seismic observations and model predictions at individual stations may arise from factors other than structural heterogeneity, such as errors in the assumed source-time function, inaccurate timings, and systematic uncertainties. To alleviate the corresponding nonuniqueness in the inverse problem, we construct differential measurements between stations, thereby reducing the influence of the source signature and systematic errors. We minimize the discrepancy between observations and simulations in terms of the differential measurements made on station pairs. We show how to implement the double-difference concept in adjoint tomography, both theoretically and in practice. We compare the sensitivities of absolute and differential measurements. The former provide absolute information on structure along the ray paths between stations and sources, whereas the latter explain relative (and thus higher-resolution) structural variations in areas close to the stations. Whereas in conventional tomography a measurement made on a single earthquake-station pair provides very limited structural information, in double-difference tomography one earthquake can actually resolve significant details of the structure. The double-difference methodology can be incorporated into the usual adjoint tomography workflow by simply pairing up all conventional measurements; the computational cost of the necessary adjoint simulations is largely unaffected. Rather than adding to the computational burden, the inversion of double-difference measurements merely modifies the construction of the adjoint sources for data assimilation.
ADGEN: ADjoint GENerator for computer models
This paper presents the development of a FORTRAN compiler and an associated supporting software library called ADGEN. ADGEN reads FORTRAN models as input and produces and enhanced version of the input model. The enhanced version reproduces the original model calculations but also has the capability to calculate derivatives of model results of interest with respect to any and all of the model data and input parameters. The method for calculating the derivatives and sensitivities is the adjoint method. Partial derivatives are calculated analytically using computer calculus and saved as elements of an adjoint matrix on direct assess storage. The total derivatives are calculated by solving an appropriate adjoint equation. ADGEN is applied to a major computer model of interest to the Low-Level Waste Community, the PRESTO-II model. PRESTO-II sample problem results reveal that ADGEN correctly calculates derivatives of response of interest with respect to 300 parameters. The execution time to create the adjoint matrix is a factor of 45 times the execution time of the reference sample problem. Once this matrix is determined, the derivatives with respect to 3000 parameters are calculated in a factor of 6.8 that of the reference model for each response of interest. For a single 3000 for determining these derivatives by parameter perturbations. The automation of the implementation of the adjoint technique for calculating derivatives and sensitivities eliminates the costly and manpower-intensive task of direct hand-implementation by reprogramming and thus makes the powerful adjoint technique more amenable for use in sensitivity analysis of existing models. 20 refs., 1 fig., 5 tabs
Application of adjoint operators to neural learning
Barhen, J.; Toomarian, N.; Gulati, S.
1990-01-01
A technique for the efficient analytical computation of such parameters of the neural architecture as synaptic weights and neural gain is presented as a single solution of a set of adjoint equations. The learning model discussed concentrates on the adiabatic approximation only. A problem of interest is represented by a system of N coupled equations, and then adjoint operators are introduced. A neural network is formalized as an adaptive dynamical system whose temporal evolution is governed by a set of coupled nonlinear differential equations. An approach based on the minimization of a constrained neuromorphic energylike function is applied, and the complete learning dynamics are obtained as a result of the calculations.
Adjoint-Based Uncertainty Quantification with MCNP
Seifried, Jeffrey E. [Univ. of California, Berkeley, CA (United States)
2011-09-01
This work serves to quantify the instantaneous uncertainties in neutron transport simulations born from nuclear data and statistical counting uncertainties. Perturbation and adjoint theories are used to derive implicit sensitivity expressions. These expressions are transformed into forms that are convenient for construction with MCNP6, creating the ability to perform adjoint-based uncertainty quantification with MCNP6. These new tools are exercised on the depleted-uranium hybrid LIFE blanket, quantifying its sensitivities and uncertainties to important figures of merit. Overall, these uncertainty estimates are small (< 2%). Having quantified the sensitivities and uncertainties, physical understanding of the system is gained and some confidence in the simulation is acquired.
Dual of QCD with One Adjoint Fermion
Mojaza, Matin; Nardecchia, Marco; Pica, Claudio;
2011-01-01
We construct the magnetic dual of QCD with one adjoint Weyl fermion. The dual is a consistent solution of the 't Hooft anomaly matching conditions, allows for flavor decoupling and remarkably constitutes the first nonsupersymmetric dual valid for any number of colors. The dual allows to bound the...
Aerosols Processes in the CMAQ Adjoint
Turner, M.; Henze, D.; Hakami, A.; Zhao, S.; Resler, Jaroslav; Carmichael, G.; Stanier, C.; Baek, J.; Saide, P.; Sandu, A.; Russel, A.; Jeong, G.; Nenes, A.; Capps, S.; Percell, P.; Pinder, R.; Napelenok, S.; Pye, H.; Bash, J.; Chai, T.; Byun, D
Davis: Air Quality Research Center, 2011. [IAMA 2011. International Aerosol Modeling Algorithms Conference /3./. Davis, 30.11.2011-02.12.2011] Institutional research plan: CEZ:AV0Z10300504 Keywords : air pollution * adjoint * aerosols Subject RIV: DG - Athmosphere Sciences, Meteorology http://dl.dropbox.com/u/41967626/IAMApresent/TURNER.pdf
Fine resolution modeling with CMAQ-adjoint
Resler, Jaroslav; Eben, Kryštof; Juruš, Pavel
Chapel Hill : CMAS, 2010. [Annual CMAS Conference /9./. 11.10.2010-13.10.2010, Chapel Hill] R&D Projects: GA MŽP SP/1A4/107/07 Institutional research plan: CEZ:AV0Z10300504 Keywords : CMAQ * adjoint * MPI * parallel efficiency Subject RIV: DG - Athmosphere Sciences, Meteorology http://www.cmascenter.org/conference/2010/agenda.cfm
Development of CO2 inversion system based on the adjoint of the global coupled transport model
Belikov, Dmitry; Maksyutov, Shamil; Chevallier, Frederic; Kaminski, Thomas; Ganshin, Alexander; Blessing, Simon
2014-05-01
We present the development of an inverse modeling system employing an adjoint of the global coupled transport model consisting of the National Institute for Environmental Studies (NIES) Eulerian transport model (TM) and the Lagrangian plume diffusion model (LPDM) FLEXPART. NIES TM is a three-dimensional atmospheric transport model, which solves the continuity equation for a number of atmospheric tracers on a grid spanning the entire globe. Spatial discretization is based on a reduced latitude-longitude grid and a hybrid sigma-isentropic coordinate in the vertical. NIES TM uses a horizontal resolution of 2.5°×2.5°. However, to resolve synoptic-scale tracer distributions and to have the ability to optimize fluxes at resolutions of 0.5° and higher we coupled NIES TM with the Lagrangian model FLEXPART. The Lagrangian component of the forward and adjoint models uses precalculated responses of the observed concentration to the surface fluxes and 3-D concentrations field simulated with the FLEXPART model. NIES TM and FLEXPART are driven by JRA-25/JCDAS reanalysis dataset. Construction of the adjoint of the Lagrangian part is less complicated, as LPDMs calculate the sensitivity of measurements to the surrounding emissions field by tracking a large number of "particles" backwards in time. Developing of the adjoint to Eulerian part was performed with automatic differentiation tool the Transformation of Algorithms in Fortran (TAF) software (http://www.FastOpt.com). This method leads to the discrete adjoint of NIES TM. The main advantage of the discrete adjoint is that the resulting gradients of the numerical cost function are exact, even for nonlinear algorithms. The overall advantages of our method are that: 1. No code modification of Lagrangian model is required, making it applicable to combination of global NIES TM and any Lagrangian model; 2. Once run, the Lagrangian output can be applied to any chemically neutral gas; 3. High-resolution results can be obtained over
Fully automatic adjoints: a robust and efficient mechanism for generating adjoint ocean models
Ham, D. A.; Farrell, P. E.; Funke, S. W.; Rognes, M. E.
2012-04-01
The problem of generating and maintaining adjoint models is sufficiently difficult that typically only the most advanced and well-resourced community ocean models achieve it. There are two current technologies which each suffer from their own limitations. Algorithmic differentiation, also called automatic differentiation, is employed by models such as the MITGCM [2] and the Alfred Wegener Institute model FESOM [3]. This technique is very difficult to apply to existing code, and requires a major initial investment to prepare the code for automatic adjoint generation. AD tools may also have difficulty with code employing modern software constructs such as derived data types. An alternative is to formulate the adjoint differential equation and to discretise this separately. This approach, known as the continuous adjoint and employed in ROMS [4], has the disadvantage that two different model code bases must be maintained and manually kept synchronised as the model develops. The discretisation of the continuous adjoint is not automatically consistent with that of the forward model, producing an additional source of error. The alternative presented here is to formulate the flow model in the high level language UFL (Unified Form Language) and to automatically generate the model using the software of the FEniCS project. In this approach it is the high level code specification which is differentiated, a task very similar to the formulation of the continuous adjoint [5]. However since the forward and adjoint models are generated automatically, the difficulty of maintaining them vanishes and the software engineering process is therefore robust. The scheduling and execution of the adjoint model, including the application of an appropriate checkpointing strategy is managed by libadjoint [1]. In contrast to the conventional algorithmic differentiation description of a model as a series of primitive mathematical operations, libadjoint employs a new abstraction of the simulation
Double-difference adjoint seismic tomography
Yuan, Yanhua O; Tromp, Jeroen
2016-01-01
We introduce a `double-difference' method for the inversion for seismic wavespeed structure based on adjoint tomography. Differences between seismic observations and model predictions at individual stations may arise from factors other than structural heterogeneity, such as errors in the assumed source-time function, inaccurate timings, and systematic uncertainties. To alleviate the corresponding nonuniqueness in the inverse problem, we construct differential measurements between stations, thereby reducing the influence of the source signature and systematic errors. We minimize the discrepancy between observations and simulations in terms of the differential measurements made on station pairs. We show how to implement the double-difference concept in adjoint tomography, both theoretically and in practice. We compare the sensitivities of absolute and differential measurements. The former provide absolute information on structure along the ray paths between stations and sources, whereas the latter explain relat...
Chiral transition of fundamental and adjoint quarks
Capdevilla, R.M. [Instituto de Física Teórica, UNESP – Universidade Estadual Paulista, Rua Dr. Bento T. Ferraz, 271, Bloco II, 01140-070 São Paulo, SP (Brazil); Doff, A., E-mail: agomes@utfpr.edu.br [Universidade Tecnológica Federal do Paraná – UTFPR – DAFIS, Av. Monteiro Lobato Km 04, 84016-210 Ponta Grossa, PR (Brazil); Natale, A.A., E-mail: natale@ift.unesp.br [Instituto de Física Teórica, UNESP – Universidade Estadual Paulista, Rua Dr. Bento T. Ferraz, 271, Bloco II, 01140-070 São Paulo, SP (Brazil); Centro de Ciências Naturais e Humanas, Universidade Federal do ABC, 09210-170 Santo André, SP (Brazil)
2014-01-20
The chiral symmetry breaking transition of quarks in the fundamental and adjoint representation is studied in a model where the gap equation contains two contributions, one containing a confining propagator and another corresponding to the exchange of one-dressed dynamically massive gluons. When quarks are in the fundamental representation the confinement effect dominates the chiral symmetry breaking while the gluon exchange is suppressed due to the dynamical gluon mass effect in the propagator and in the coupling constant. In this case the chiral and deconfinement transition temperatures are approximately the same. For quarks in the adjoint representation, due to the larger Casimir eigenvalue, the gluon exchange is operative and the chiral transition happens at a larger temperature than the deconfinement one.
Chiral transition of fundamental and adjoint quarks
The chiral symmetry breaking transition of quarks in the fundamental and adjoint representation is studied in a model where the gap equation contains two contributions, one containing a confining propagator and another corresponding to the exchange of one-dressed dynamically massive gluons. When quarks are in the fundamental representation the confinement effect dominates the chiral symmetry breaking while the gluon exchange is suppressed due to the dynamical gluon mass effect in the propagator and in the coupling constant. In this case the chiral and deconfinement transition temperatures are approximately the same. For quarks in the adjoint representation, due to the larger Casimir eigenvalue, the gluon exchange is operative and the chiral transition happens at a larger temperature than the deconfinement one
Effective freeness of adjoint line bundles
Heier, Gordon
2001-01-01
In this note we establish a new Fujita-type effective bound for the base point freeness of adjoint line bundles on a compact complex projective manifold of complex dimension $n$. The bound we obtain (approximately) differs from the linear bound conjectured by Fujita only by a factor of the cube root of $n$. As an application, a new effective statement for pluricanonical embeddings is derived.
Seismic imaging: From classical to adjoint tomography
Liu, Q.; Gu, Y. J.
2012-09-01
Seismic tomography has been a vital tool in probing the Earth's internal structure and enhancing our knowledge of dynamical processes in the Earth's crust and mantle. While various tomographic techniques differ in data types utilized (e.g., body vs. surface waves), data sensitivity (ray vs. finite-frequency approximations), and choices of model parameterization and regularization, most global mantle tomographic models agree well at long wavelengths, owing to the presence and typical dimensions of cold subducted oceanic lithospheres and hot, ascending mantle plumes (e.g., in central Pacific and Africa). Structures at relatively small length scales remain controversial, though, as will be discussed in this paper, they are becoming increasingly resolvable with the fast expanding global and regional seismic networks and improved forward modeling and inversion techniques. This review paper aims to provide an overview of classical tomography methods, key debates pertaining to the resolution of mantle tomographic models, as well as to highlight recent theoretical and computational advances in forward-modeling methods that spearheaded the developments in accurate computation of sensitivity kernels and adjoint tomography. The first part of the paper is devoted to traditional traveltime and waveform tomography. While these approaches established a firm foundation for global and regional seismic tomography, data coverage and the use of approximate sensitivity kernels remained as key limiting factors in the resolution of the targeted structures. In comparison to classical tomography, adjoint tomography takes advantage of full 3D numerical simulations in forward modeling and, in many ways, revolutionizes the seismic imaging of heterogeneous structures with strong velocity contrasts. For this reason, this review provides details of the implementation, resolution and potential challenges of adjoint tomography. Further discussions of techniques that are presently popular in
Chiral transition of fundamental and adjoint quarks
Capdevilla, R. M.; Doff, A.(Universidade Tecnológica Federal do Paraná – UTFPR – DAFIS, Av. Monteiro Lobato Km 04, 84016-210 Ponta Grossa, PR, Brazil); Natale, A. A.
2014-01-01
The chiral symmetry breaking transition of quarks in the fundamental and adjoint representation is studied in a model where the gap equation contains two contributions, one containing a confining propagator and another corresponding to the exchange of one-dressed dynamically massive gluons. When quarks are in the fundamental representation the confinement effect dominates the chiral symmetry breaking while the gluon exchange is suppressed due to the dynamical gluon mass effect in the propagat...
Fast Correlation Greeks by Adjoint Algorithmic Differentiation
Luca Capriotti; Mike Giles
2010-01-01
We show how Adjoint Algorithmic Differentiation (AAD) allows an extremely efficient calculation of correlation Risk of option prices computed with Monte Carlo simulations. A key point in the construction is the use of binning to simultaneously achieve computational efficiency and accurate confidence intervals. We illustrate the method for a copula-based Monte Carlo computation of claims written on a basket of underlying assets, and we test it numerically for Portfolio Default Options. For any...
Supersymmetric descendants of self-adjointly extended quantum mechanical Hamiltonians
Al-Hashimi, M. H.; Salman, M.; Shalaby, A.; Wiese, U.-J.
2013-10-01
We consider the descendants of self-adjointly extended Hamiltonians in supersymmetric quantum mechanics on a half-line, on an interval, and on a punctured line or interval. While there is a 4-parameter family of self-adjointly extended Hamiltonians on a punctured line, only a 3-parameter sub-family has supersymmetric descendants that are themselves self-adjoint. We also address the self-adjointness of an operator related to the supercharge, and point out that only a sub-class of its most general self-adjoint extensions is physical. Besides a general characterization of self-adjoint extensions and their supersymmetric descendants, we explicitly consider concrete examples, including a particle in a box with general boundary conditions, with and without an additional point interaction. We also discuss bulk-boundary resonances and their manifestation in the supersymmetric descendant.
Supersymmetric Descendants of Self-Adjointly Extended Quantum Mechanical Hamiltonians
Al-Hashimi, M H; Shalaby, A; Wiese, U -J
2013-01-01
We consider the descendants of self-adjointly extended Hamiltonians in supersymmetric quantum mechanics on a half-line, on an interval, and on a punctured line or interval. While there is a 4-parameter family of self-adjointly extended Hamiltonians on a punctured line, only a 3-parameter sub-family has supersymmetric descendants that are themselves self-adjoint. We also address the self-adjointness of an operator related to the supercharge, and point out that only a sub-class of its most general self-adjoint extensions is physical. Besides a general characterization of self-adjoint extensions and their supersymmetric descendants, we explicitly consider concrete examples, including a particle in a box with general boundary conditions, with and without an additional point interaction. We also discuss bulk-boundary resonances and their manifestation in the supersymmetric descendant.
D. A. Belikov
2015-07-01
Full Text Available We present the development of the Adjoint of the Global Eulerian–Lagrangian Coupled Atmospheric (A-GELCA model that consists of the National Institute for Environmental Studies (NIES model as an Eulerian three-dimensional transport model (TM, and FLEXPART (FLEXible PARTicle dispersion model as the Lagrangian plume diffusion model (LPDM. The tangent and adjoint components of the Eulerian model were constructed directly from the original NIES TM code using an automatic differentiation tool known as TAF (Transformation of Algorithms in Fortran; http://www.FastOpt.com, with additional manual pre- and post-processing aimed at improving the performance of the computing, including MPI (Message Passing Interface. As results, the adjoint of Eulerian model is discrete. Construction of the adjoint of the Lagrangian component did not require any code modification, as LPDMs are able to track a significant number of particles back in time and thereby calculate the sensitivity of observations to the neighboring emissions areas. Eulerian and Lagrangian adjoint components were coupled at the time boundary in the global domain.The results are verified using a series of test experiments. The forward simulation shown the coupled model is effective in reproducing the seasonal cycle and short-term variability of CO2 even in the case of multiple limiting factors, such as high uncertainty of fluxes and the low resolution of the Eulerian model. The adjoint model demonstrates the high accuracy compared to direct forward sensitivity calculations and fast performance. The developed adjoint of the coupled model combines the flux conservation and stability of an Eulerian discrete adjoint formulation with the flexibility, accuracy, and high resolution of a Lagrangian backward trajectory formulation.
Belikov, Dmitry A.; Maksyutov, Shamil; Yaremchuk, Alexey; Ganshin, Alexander; Kaminski, Thomas; Blessing, Simon; Sasakawa, Motoki; Gomez-Pelaez, Angel J.; Starchenko, Alexander
2016-02-01
We present the development of the Adjoint of the Global Eulerian-Lagrangian Coupled Atmospheric (A-GELCA) model that consists of the National Institute for Environmental Studies (NIES) model as an Eulerian three-dimensional transport model (TM), and FLEXPART (FLEXible PARTicle dispersion model) as the Lagrangian Particle Dispersion Model (LPDM). The forward tangent linear and adjoint components of the Eulerian model were constructed directly from the original NIES TM code using an automatic differentiation tool known as TAF (Transformation of Algorithms in Fortran; http://www.FastOpt.com, with additional manual pre- and post-processing aimed at improving transparency and clarity of the code and optimizing the performance of the computing, including MPI (Message Passing Interface). The Lagrangian component did not require any code modification, as LPDMs are self-adjoint and track a significant number of particles backward in time in order to calculate the sensitivity of the observations to the neighboring emission areas. The constructed Eulerian adjoint was coupled with the Lagrangian component at a time boundary in the global domain. The simulations presented in this work were performed using the A-GELCA model in forward and adjoint modes. The forward simulation shows that the coupled model improves reproduction of the seasonal cycle and short-term variability of CO2. Mean bias and standard deviation for five of the six Siberian sites considered decrease roughly by 1 ppm when using the coupled model. The adjoint of the Eulerian model was shown, through several numerical tests, to be very accurate (within machine epsilon with mismatch around to ±6 e-14) compared to direct forward sensitivity calculations. The developed adjoint of the coupled model combines the flux conservation and stability of an Eulerian discrete adjoint formulation with the flexibility, accuracy, and high resolution of a Lagrangian backward trajectory formulation. A-GELCA will be incorporated
Generalized uncertainty principle and self-adjoint operators
Balasubramanian, Venkat, E-mail: vbalasu8@uwo.ca [Department of Applied Mathematics, University of Western Ontario London, Ontario N6A 5B7 (Canada); Das, Saurya, E-mail: saurya.das@uleth.ca [Theoretical Physics Group, Department of Physics and Astronomy, University of Lethbridge, 4401 University Drive, Lethbridge, Alberta T1K 3M4 (Canada); Vagenas, Elias C., E-mail: elias.vagenas@ku.edu.kw [Theoretical Physics Group, Department of Physics, Kuwait University, P.O. Box 5969, Safat 13060 (Kuwait)
2015-09-15
In this work we explore the self-adjointness of the GUP-modified momentum and Hamiltonian operators over different domains. In particular, we utilize the theorem by von-Neumann for symmetric operators in order to determine whether the momentum and Hamiltonian operators are self-adjoint or not, or they have self-adjoint extensions over the given domain. In addition, a simple example of the Hamiltonian operator describing a particle in a box is given. The solutions of the boundary conditions that describe the self-adjoint extensions of the specific Hamiltonian operator are obtained.
On the essential self-adjointness of generalized Schroedinger operators
We give a necessary and sufficient condition for the generalized Schroedinger operator to be essentially self-adjoint in L2(Ω; rhodx), under general assumptions on rho and for arbitrary domains Ω in Rsup(n). In particular, if rho is strictly positive and locally Lipschitz continuous on Ω = Rsup(n), then A is essentially self-adjoint. We also give examples of non-essential self-adjointness and a complete discussion of the one-dimensional case. These results have applications to the problem of the essential self-adjointness of quantum Hamiltonians and to the uniqueness problem of Markov processes. (orig./WL)
The Schroedinger differential operator and intrinsic self adjointness
An attempt is made to establish a hermitization procedure for rendering any linear differential operator to be intrinsically self adjoint, independently of any prescribed representation. This is accomplished by introducing an associate differential operator, which is simply a linear combination of all the individual ordinary differential operators belonging to the adjoint of the given linear differential operator, that satisfies the criterion of intrinsic self adjointness. It turns out that the associate differential operator is capable of generating an infinite set of hermitized versions of any arbitrary linear differential operator. Both momentum and kinetic energy differential operators that belong to the Schroedinger wave equation are rendered self adjoint. (author)
Symmetries of linearized gravity from adjoint operators
Aksteiner, Steffen
2016-01-01
Using a covariant formulation it is shown that the Teukolsky equation and the Teukolsky-Starobinsky identities for spin-1 and linearized gravity on a vacuum type D background are self-adjoint. This fact is used to construct symmetry operators for each of the four cases. We find both irreducible second order symmetry operators for spin-1, a known fourth order, and a new sixth order symmetry operator for linearized gravity. The results are connected to Hertz and Debye potentials and to the separability of the Teukolsky equation.
Local Volatility Calibration Using An Adjoint Proxy
Gabriel TURINICI
2008-11-01
Full Text Available We document the calibration of the local volatility in a framework similar to Coleman, Li and Verma. The quality of a surface is assessed through a functional to be optimized; the specificity of the approach is to separate the optimization (performed with any suitable optimization algorithm from the computation of the functional where we use an adjoint (as in L. Jiang et. al. to obtain an approximation; moreover our main calibration variable is the implied volatility (the procedure can also accommodate the Greeks. The procedure performs well on benchmarks from the literature and on FOREX data.
Supersymmetric descendants of self-adjointly extended quantum mechanical Hamiltonians
We consider the descendants of self-adjointly extended Hamiltonians in supersymmetric quantum mechanics on a half-line, on an interval, and on a punctured line or interval. While there is a 4-parameter family of self-adjointly extended Hamiltonians on a punctured line, only a 3-parameter sub-family has supersymmetric descendants that are themselves self-adjoint. We also address the self-adjointness of an operator related to the supercharge, and point out that only a sub-class of its most general self-adjoint extensions is physical. Besides a general characterization of self-adjoint extensions and their supersymmetric descendants, we explicitly consider concrete examples, including a particle in a box with general boundary conditions, with and without an additional point interaction. We also discuss bulk-boundary resonances and their manifestation in the supersymmetric descendant. -- Highlights: •Self-adjoint extension theory and contact interactions. •Application of self-adjoint extensions to supersymmetry. •Contact interactions in finite volume with Robin boundary condition
Supersymmetric descendants of self-adjointly extended quantum mechanical Hamiltonians
Al-Hashimi, M.H., E-mail: hashimi@itp.unibe.ch [Albert Einstein Center for Fundamental Physics, Institute for Theoretical Physics, Bern University, Sidlerstrasse 5, CH-3012 Bern (Switzerland); Salman, M., E-mail: msalman@qu.edu.qa [Department of Mathematics, Statistics, and Physics, Qatar University, Al Tarfa, Doha 2713 (Qatar); Shalaby, A., E-mail: amshalab@qu.edu.qa [Department of Mathematics, Statistics, and Physics, Qatar University, Al Tarfa, Doha 2713 (Qatar); Physics Department, Faculty of Science, Mansoura University (Egypt); Wiese, U.-J., E-mail: wiese@itp.unibe.ch [Albert Einstein Center for Fundamental Physics, Institute for Theoretical Physics, Bern University, Sidlerstrasse 5, CH-3012 Bern (Switzerland); Center for Theoretical Physics, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA (United States)
2013-10-15
We consider the descendants of self-adjointly extended Hamiltonians in supersymmetric quantum mechanics on a half-line, on an interval, and on a punctured line or interval. While there is a 4-parameter family of self-adjointly extended Hamiltonians on a punctured line, only a 3-parameter sub-family has supersymmetric descendants that are themselves self-adjoint. We also address the self-adjointness of an operator related to the supercharge, and point out that only a sub-class of its most general self-adjoint extensions is physical. Besides a general characterization of self-adjoint extensions and their supersymmetric descendants, we explicitly consider concrete examples, including a particle in a box with general boundary conditions, with and without an additional point interaction. We also discuss bulk-boundary resonances and their manifestation in the supersymmetric descendant. -- Highlights: •Self-adjoint extension theory and contact interactions. •Application of self-adjoint extensions to supersymmetry. •Contact interactions in finite volume with Robin boundary condition.
Self-adjointness of the Gaffney Laplacian on Vector Bundles
We study the Gaffney Laplacian on a vector bundle equipped with a compatible metric and connection over a Riemannian manifold that is possibly geodesically incomplete. Under the hypothesis that the Cauchy boundary is polar, we demonstrate the self-adjointness of this Laplacian. Furthermore, we show that negligible boundary is a necessary and sufficient condition for the self-adjointness of this operator
A reduced adjoint approach to variational data assimilation
Altaf, Muhammad
2013-02-01
The adjoint method has been used very often for variational data assimilation. The computational cost to run the adjoint model often exceeds several original model runs and the method needs significant programming efforts to implement the adjoint model code. The work proposed here is variational data assimilation based on proper orthogonal decomposition (POD) which avoids the implementation of the adjoint of the tangent linear approximation of the original nonlinear model. An ensemble of the forward model simulations is used to determine the approximation of the covariance matrix and only the dominant eigenvectors of this matrix are used to define a model subspace. The adjoint of the tangent linear model is replaced by the reduced adjoint based on this reduced space. Thus the adjoint model is run in reduced space with negligible computational cost. Once the gradient is obtained in reduced space it is projected back in full space and the minimization process is carried in full space. In the paper the reduced adjoint approach to variational data assimilation is introduced. The characteristics and performance of the method are illustrated with a number of data assimilation experiments in a ground water subsurface contaminant model. © 2012 Elsevier B.V.
Self-adjointness of the Gaffney Laplacian on Vector Bundles
Bandara, Lashi, E-mail: lashi.bandara@chalmers.se [Chalmers University of Technology and University of Gothenburg, Mathematical Sciences (Sweden); Milatovic, Ognjen, E-mail: omilatov@unf.edu [University of North Florida, Department of Mathematics and Statistics (United States)
2015-12-15
We study the Gaffney Laplacian on a vector bundle equipped with a compatible metric and connection over a Riemannian manifold that is possibly geodesically incomplete. Under the hypothesis that the Cauchy boundary is polar, we demonstrate the self-adjointness of this Laplacian. Furthermore, we show that negligible boundary is a necessary and sufficient condition for the self-adjointness of this operator.
Radiation source reconstruction with known geometry and materials using the adjoint
We present a method to estimate an unknown isotropic source distribution, in space and energy, using detector measurements when the geometry and material composition are known. The estimated source distribution minimizes the difference between the measured and computed responses of detectors located at a selected number of points within the domain. In typical methods, a forward flux calculation is performed for each source guess in an iterative process. In contrast, we use the adjoint flux to compute the responses. Potential applications of the proposed method include determining the distribution of radio-contaminants following a nuclear event, monitoring the flow of radioactive fluids in pipes to determine hold-up locations, and retroactive reconstruction of radiation fields using workers' detectors' readings. After presenting the method, we describe a numerical test problem to demonstrate the preliminary viability of the method. As expected, using the adjoint flux reduces the number of transport solves to be proportional to the number of detector measurements, in contrast to methods using the forward flux that require a typically larger number proportional to the number of spatial mesh cells. (author)
Adjoint Monte Carlo simulation of fixed-energy secondary radiation
Fixed energy secondary generation for adjoint Monte Carlo methods constitutes certain difficulties because of zero probability of reaching fixed value from continuous distribution. This paper proposes a possible approach to adjoint Monte Carlo simulation with fixed energy secondary radiation which does not contain any simplifying restriction. This approach uses the introduced before generalized particle concept developed for description of mixed-type radiation transport and allows adjoint Monte Carlo simulation of such processes. It treats particle type as additional discrete coordinate and always considers only one particle even for the interactions with many particles outgoing from the collision. The adjoint fixed energy secondary radiation simulation is performed as local energy estimator through the intermediate state with fixed energy. The proposed algorithm is tested on the example of coupled gamma/electron/positron transport with generation of annihilation radiation. Forward and adjoint simulation according to generalized particle concept show statistically similar results. (orig.)
Shi, Lei; Wang, Z. J.
2015-08-01
Adjoint-based mesh adaptive methods are capable of distributing computational resources to areas which are important for predicting an engineering output. In this paper, we develop an adjoint-based h-adaptation approach based on the high-order correction procedure via reconstruction formulation (CPR) to minimize the output or functional error. A dual-consistent CPR formulation of hyperbolic conservation laws is developed and its dual consistency is analyzed. Super-convergent functional and error estimate for the output with the CPR method are obtained. Factors affecting the dual consistency, such as the solution point distribution, correction functions, boundary conditions and the discretization approach for the non-linear flux divergence term, are studied. The presented method is then used to perform simulations for the 2D Euler and Navier-Stokes equations with mesh adaptation driven by the adjoint-based error estimate. Several numerical examples demonstrate the ability of the presented method to dramatically reduce the computational cost comparing with uniform grid refinement.
Adjoint-based uncertainty quantification and sensitivity analysis for reactor depletion calculations
Stripling, Hayes Franklin
Depletion calculations for nuclear reactors model the dynamic coupling between the material composition and neutron flux and help predict reactor performance and safety characteristics. In order to be trusted as reliable predictive tools and inputs to licensing and operational decisions, the simulations must include an accurate and holistic quantification of errors and uncertainties in its outputs. Uncertainty quantification is a formidable challenge in large, realistic reactor models because of the large number of unknowns and myriad sources of uncertainty and error. We present a framework for performing efficient uncertainty quantification in depletion problems using an adjoint approach, with emphasis on high-fidelity calculations using advanced massively parallel computing architectures. This approach calls for a solution to two systems of equations: (a) the forward, engineering system that models the reactor, and (b) the adjoint system, which is mathematically related to but different from the forward system. We use the solutions of these systems to produce sensitivity and error estimates at a cost that does not grow rapidly with the number of uncertain inputs. We present the framework in a general fashion and apply it to both the source-driven and k-eigenvalue forms of the depletion equations. We describe the implementation and verification of solvers for the forward and ad- joint equations in the PDT code, and we test the algorithms on realistic reactor analysis problems. We demonstrate a new approach for reducing the memory and I/O demands on the host machine, which can be overwhelming for typical adjoint algorithms. Our conclusion is that adjoint depletion calculations using full transport solutions are not only computationally tractable, they are the most attractive option for performing uncertainty quantification on high-fidelity reactor analysis problems.
P. E. Haines
2014-01-01
Full Text Available A new methodology for the formulation of an adjoint to the transport component of the chemistry transport model TOMCAT is described and implemented in a new model RETRO-TOM. The Eulerian backtracking method is used, allowing the forward advection scheme (Prather's second-order moments, to be efficiently exploited in the backward adjoint calculations. Prather's scheme is shown to be time-symmetric suggesting the possibility of high accuracy. To attain this accuracy, however, it is necessary to make a careful treatment of the "density inconsistency" problem inherent to offline transport models. The results are verified using a series of test experiments. These demonstrate the high accuracy of RETRO-TOM when compared with direct forward sensitivity calculations, at least for problems in which flux-limiters in the advection scheme are not required. RETRO-TOM therefore combines the flexibility and stability of a "finite difference of adjoint" formulation with the accuracy of an "adjoint of finite difference" formulation.
Gauge Mediation Models with Adjoint Messengers
Gogoladze, Ilia; Shafi, Qaisar; Un, Cem Salih
2016-01-01
We present a class of models in the framework of gauge mediation supersymmetry breaking where the messenger fields transform in the adjoint representation of the Standard Model gauge symmetry. To avoid unacceptably light right-handed sleptons in the spectrum we introduce a non-zero U(1)_B-L D-term. This leads to an additional contribution to the soft supersymmetry breaking mass terms which makes the right-handed slepton masses compatible with the current experimental bounds. We show that in this framework the observed 125 GeV Higgs boson mass can be accommodated with the sleptons accessible at the LHC, while the squarks and gluinos lie in the multi-TeV range. We also discuss the issue of the fine-tuning and show that the desired relic dark matter abundance can also be accommodated.
Adjoint tomography of the southern California crust.
Tape, Carl; Liu, Qinya; Maggi, Alessia; Tromp, Jeroen
2009-08-21
Using an inversion strategy based on adjoint methods, we developed a three-dimensional seismological model of the southern California crust. The resulting model involved 16 tomographic iterations, which required 6800 wavefield simulations and a total of 0.8 million central processing unit hours. The new crustal model reveals strong heterogeneity, including local changes of +/-30% with respect to the initial three-dimensional model provided by the Southern California Earthquake Center. The model illuminates shallow features such as sedimentary basins and compositional contrasts across faults. It also reveals crustal features at depth that aid in the tectonic reconstruction of southern California, such as subduction-captured oceanic crustal fragments. The new model enables more realistic and accurate assessments of seismic hazard. PMID:19696349
GPU-accelerated adjoint algorithmic differentiation
Gremse, Felix; Höfter, Andreas; Razik, Lukas; Kiessling, Fabian; Naumann, Uwe
2016-03-01
Many scientific problems such as classifier training or medical image reconstruction can be expressed as minimization of differentiable real-valued cost functions and solved with iterative gradient-based methods. Adjoint algorithmic differentiation (AAD) enables automated computation of gradients of such cost functions implemented as computer programs. To backpropagate adjoint derivatives, excessive memory is potentially required to store the intermediate partial derivatives on a dedicated data structure, referred to as the "tape". Parallelization is difficult because threads need to synchronize their accesses during taping and backpropagation. This situation is aggravated for many-core architectures, such as Graphics Processing Units (GPUs), because of the large number of light-weight threads and the limited memory size in general as well as per thread. We show how these limitations can be mediated if the cost function is expressed using GPU-accelerated vector and matrix operations which are recognized as intrinsic functions by our AAD software. We compare this approach with naive and vectorized implementations for CPUs. We use four increasingly complex cost functions to evaluate the performance with respect to memory consumption and gradient computation times. Using vectorization, CPU and GPU memory consumption could be substantially reduced compared to the naive reference implementation, in some cases even by an order of complexity. The vectorization allowed usage of optimized parallel libraries during forward and reverse passes which resulted in high speedups for the vectorized CPU version compared to the naive reference implementation. The GPU version achieved an additional speedup of 7.5 ± 4.4, showing that the processing power of GPUs can be utilized for AAD using this concept. Furthermore, we show how this software can be systematically extended for more complex problems such as nonlinear absorption reconstruction for fluorescence-mediated tomography.
A PNJL Model for Adjoint Fermions with Periodic Boundary Conditions
Nishimura, Hiromichi; Ogilvie, Michael C.
2009-01-01
Recent work on QCD-like theories has shown that the addition of adjoint fermions obeying periodic boundary conditions to gauge theories on $R^{3}\\times S^{1}$ can lead to a restoration of center symmetry and confinement for sufficiently small circumference $L$ of $S^{1}$. At small $L$, perturbation theory may be used reliably to compute the effective potential for the Polyakov loop $P$ in the compact direction. Periodic adjoint fermions act in opposition to the gauge fields, which by themselv...
One-energy group, two-dimensional computer code was developed to calculate the response of a detector to a vibrating absorber in a reactor core. A concept of local/global components, based on the frequency dependent detector adjoint function, and a nodalization technique were utilized. The frequency dependent detector adjoint functions presented by complex equations were expanded into real and imaginary parts. In the nodalization technique, the flux is expanded into polynomials about the center point of each node. The phase angle and the magnitude of the one-energy group detector adjoint function were calculated for a detector located in the center of a 200x200 cm reactor using a two-dimensional nodalization technique, the computer code EXTERMINATOR, and the analytical solution. The purpose of this research was to investigate the applicability of a polynomial nodal model technique to the calculations of the real and the imaginary parts of the detector adjoint function for one-energy group two-dimensional polynomial nodal model technique. From the results as discussed earlier, it is concluded that the nodal model technique can be used to calculate the detector adjoint function and the phase angle. Using the computer code developed for nodal model technique, the magnitude of one energy group frequency dependent detector adjoint function and the phase angle were calculated for the detector located in the center of a 200x200 cm homogenous reactor. The real part of the detector adjoint function was compared with the results obtained from the EXTERMINATOR computer code as well as the analytical solution based on a double sine series expansion using the classical Green's Function solution. The values were found to be less than 1% greater at 20 cm away from the source region and about 3% greater closer to the source compared to the values obtained from the analytical solution and the EXTERMINATOR code. The currents at the node interface matched within 1% of the average
Utilisation de sources et d'adjoints dragon pour les calculs TRIPOLI
Camand, Corentin
usually non significant. The second method is to use of the adjoint neutron flux calculated by DRAGON as an importance function for Monte Carlo biaising in TRIPOLI. The objective is to improve the figure of merit of the detector response located far away of the neutron source. The neutron source initialisation of a TRIPOLI calculation required to develop the development of a module in DRAGON that generates a list of sources in the TRIPOLI syntaxe, including for each source, its intensity, its position and the energy domain it covers. We tested our method on a complete 17×17 PWR-UOX assembly and on a reduced 3×3 model. We first verified that the DRAGON and TRIPOLI models were consistent in order to ensure that TRIPOLI receives a coherent source distribution. Then we tested the use of DRAGON sources in TRIPOLI with neutron flux and the effective multiplying coefficient (keff). We observe slightly better standard deviations, of an order of 10 pcm, on keff for simulations using DRAGON sources distributions as compared to simulations with less precise initial sources. Flux convergence is also improved. However some incoherence were also observed in the results, some flux converging slower with DRAGON sources when fewer neutrons per batch are considered. In addition, a very large number of sources is too heavy to insert in TRIPOLI. It seems that our method is perfectible in order to improve implementation and convergence. Study of more complex geometries, with less regular sources distributions (for instance using MOX or irradiated fuel) may provide better performances using our method. For biaising TRIPOLI calculations using the DRAGON adjoint flux we created a module that produces importance maps readable by TRIPOLI. We tested our method on a source-detector shielding problem in one dimension. After checking the coherence of DRAGON and TRIPOLI models, we biaised TRIPOLI simulations using the DRAGON adjoint flux, and using INIPOND, the internal biaising option of TRIPOLI. We
Deterministic adjoint transport applications for He-3 neutron detector design
This work focuses on the determination of predicted neutron detector response accomplished using neutron importance derived from an adjoint discrete ordinates (SN) transport calculation. A hypothetical detector apparatus, intended to detect fast neutrons, was modeled using He-3 tubes with graphite moderation using the PENTRANTM 3-D multi-group discrete ordinates parallel transport code system. The detector geometry was modeled using z-axis symmetry and discretized into 30,280 3-D Cartesian cells. The material spatial mesh was generated using the PENMSHTM code in the PENTRAN system. The 47-group BUGLE-96 neutron cross section library was used for construction of macroscopic neutron cross sections. Results from an S8 angular quadrature using P3 anisotropy are presented. An adjoint transport source was established in the model using group dependent He-3 response cross sections. Each He-3 tube contained an adjoint source aliased to group He-3 absorption cross sections to permit assessment of detector performance. The spectrally dependent detector response from neutron capture in He-3 tubes from an arbitrary source can, therefore, be readily determined. This response comes from the complete integral of the actual source strength weighted by the adjoint function at the source location for any source distribution scenario. For selected neutron energies, an equivalent forward MCNP Monte Carlo model was used to demonstrate good agreement with the detector response determined from the adjoint calculation. The graphite used in this design has a large impact on detector performance due to the increasing sensitivity inherent in He-3 gas as neutrons thermalize. Computational adjoint results presented here predict a fast neutron detector design that yields efficiencies between 30 and 50% for neutron energies below 3 keV, and up to 30% efficiencies for neutron energies between 3 keV and 1 MeV. Overall, the methodology applied here highlights the elegant nature of an adjoint
Flux tubes at Finite Temperature
Bicudo, Pedro; Cardoso, Marco
2016-01-01
We show the flux tubes produced by static quark-antiquark, quark-quark and quark-gluon charges at finite temperature. The sources are placed in the lattice with fundamental and adjoint Polyakov loops. We compute the square densities of the chromomagnetic and chromoelectric fields above and below the phase transition. Our results are gauge invariant and produced in pure gauge SU(3). The codes are written in CUDA and the computations are performed with GPUs.
Adjoint transport methods for radiation-effects testing
Adjoint transport has been exploited for some time for neutral particle calculations. For charged particles, however, production adjoint capability was not available until Morel developed the ability to solve coupled-photon-electron transport problems with production discrete ordinates codes. This represents a significant advance for many problems of interest, such as predicting bremsstrahlung yield from flash X-ray machines, internal electromagnetic pulse (IEMP) for photons incident on printed circuit boards, shielding requirements for electron dosimetry, and dose enhancement from photon irradiation of printed circuit boards. The authors demonstrate here that adjoint photon-electron transport is at least an order of magnitude more efficient than forward transport for optimizing bremsstrahlung yield from flash X-ray machine converters. This problem is particularly interesting since adjoint transport provides a good approximation for a variable geometry in addition to a variable source, due to the highly forward-peaked nature of the electron scattering. Normally, neither forward nor adjoint transport is efficient for studying a variable-geometry problem
Adjoint-based Optimal Flow Control for Compressible DNS
Otero, J Javier; Sandberg, Richard D
2016-01-01
A novel adjoint-based framework oriented to optimal flow control in compressible direct numerical simulations is presented. Also, a new formulation of the adjoint characteristic boundary conditions is introduced, which enhances the stability of the adjoint simulations. The flow configuration chosen as a case study consists of a two dimensional open cavity flow with aspect ratio $L/H=3$ and Reynolds number $Re=5000$. This flow configuration is of particular interest, as the turbulent and chaotic nature of separated flows pushes the adjoint approach to its limit. The target of the flow actuation, defined as cost, is the reduction of the pressure fluctuations at the sensor location. To exploit the advantages of the adjoint method, a large number of control parameters is used. The control consists of an actuating sub-domain where a two-dimensional body force is applied at every point within the sub-volume. This results in a total of $2.256 \\cdot 10^6$ control parameters. The final actuation achieved a successful ...
Mesh-free adjoint methods for nonlinear filters
Daum, Fred
2005-09-01
We apply a new industrial strength numerical approximation, called the "mesh-free adjoint method", to solve the nonlinear filtering problem. This algorithm exploits the smoothness of the problem, unlike particle filters, and hence we expect that mesh-free adjoints are superior to particle filters for many practical applications. The nonlinear filter problem is equivalent to solving the Fokker-Planck equation in real time. The key idea is to use a good adaptive non-uniform quantization of state space to approximate the solution of the Fokker-Planck equation. In particular, the adjoint method computes the location of the nodes in state space to minimize errors in the final answer. This use of an adjoint is analogous to optimal control algorithms, but it is more interesting. The adjoint method is also analogous to importance sampling in particle filters, but it is better for four reasons: (1) it exploits the smoothness of the problem; (2) it explicitly minimizes the errors in the relevant functional; (3) it explicitly models the dynamics in state space; and (4) it can be used to compute a corrected value for the desired functional using the residuals. We will attempt to make this paper accessible to normal engineers who do not have PDEs for breakfast.
Andrews, Arlyn; Kawa, Randy; Zhu, Zhengxin; Burris, John; Abshire, Jim
2004-01-01
A detailed mechanistic understanding of the sources and sinks of CO2 will be required to reliably predict future CO2 levels and climate. A commonly used technique for deriving information about CO2 exchange with surface reservoirs is to solve an 'inverse problem', where CO2 observations are used with an atmospheric transport model to find the optimal distribution of sources and sinks. Synthesis inversion methods are powerful tools for addressing this question, but the results are disturbingly sensitive to the details of the calculation. Studies done using different atmospheric transport models and combinations of surface station data have produced substantially different distributions of surface fluxes. Adjoint methods are now being developed that will more effectively incorporate diverse datasets in estimates of surface fluxes of CO2. In an adjoint framework, it will be possible to combine CO2 concentration data from longterm surface and aircraft monitoring stations with data from intensive field campaigns and with proposed future satellite observations. We have recently developed an adjoint for the GSFC 3-D Parameterized Chemistry and Transport Model (PCTM). Here, we will present results from a PCTM Adjoint study comparing the sampling footprints of tall tower, aircraft and potential future lidar observations of CO2. The vertical resolution and extent of the profiles and the observation frequency will be considered for several sites in North America.
Adjoint-consistent formulations of slip models for coupled electroosmotic flow systems
Garg, Vikram V
2014-09-27
Background Models based on the Helmholtz `slip\\' approximation are often used for the simulation of electroosmotic flows. The objectives of this paper are to construct adjoint-consistent formulations of such models, and to develop adjoint-based numerical tools for adaptive mesh refinement and parameter sensitivity analysis. Methods We show that the direct formulation of the `slip\\' model is adjoint inconsistent, and leads to an ill-posed adjoint problem. We propose a modified formulation of the coupled `slip\\' model, which is shown to be well-posed, and therefore automatically adjoint-consistent. Results Numerical examples are presented to illustrate the computation and use of the adjoint solution in two-dimensional microfluidics problems. Conclusions An adjoint-consistent formulation for Helmholtz `slip\\' models of electroosmotic flows has been proposed. This formulation provides adjoint solutions that can be reliably used for mesh refinement and sensitivity analysis.
A Posteriori Analysis for Hydrodynamic Simulations Using Adjoint Methodologies
Woodward, C S; Estep, D; Sandelin, J; Wang, H
2009-02-26
This report contains results of analysis done during an FY08 feasibility study investigating the use of adjoint methodologies for a posteriori error estimation for hydrodynamics simulations. We developed an approach to adjoint analysis for these systems through use of modified equations and viscosity solutions. Targeting first the 1D Burgers equation, we include a verification of the adjoint operator for the modified equation for the Lax-Friedrichs scheme, then derivations of an a posteriori error analysis for a finite difference scheme and a discontinuous Galerkin scheme applied to this problem. We include some numerical results showing the use of the error estimate. Lastly, we develop a computable a posteriori error estimate for the MAC scheme applied to stationary Navier-Stokes.
SUPERSTABILITY OF ADJOINTABLE MAPPINGS ON HILBERT C*-MODULES
Mohammad Sal Moslehian
2009-02-01
Full Text Available We define the notion of $varphi$-perturbation of a densely definedadjointable mapping and prove that any such mapping $f$ betweenHilbert ${mathcal A}$-modules over a fixed $C^*$-algebra ${mathcalA}$ with densely defined corresponding mapping $g$ is ${mathcalA}$-linear and adjointable in the classical sense with adjoint $g$.If both $f$ and $g$ are everywhere defined then they are bounded.Our work concerns with the concept of {sc Hyers--Ulam--Rassias} stability originated from the {sc Th.~M.~Rassias}' stability theorem that appeared in his paper [{it On the stability of the linear mapping in Banach spaces}, Proc. Amer. Math. Soc., {f 72} (1978, 297--300]. We also indicate complementary results in the case where the {sc Hilbert} $C^*$-modules admit non-adjointable $C^*$-linear appings
Automatic differentiation, tangent linear models, and (pseudo) adjoints
Bischof, C.H.
1993-12-31
This paper provides a brief introduction to automatic differentiation and relates it to the tangent linear model and adjoint approaches commonly used in meteorology. After a brief review of the forward and reverse mode of automatic differentiation, the ADIFOR automatic differentiation tool is introduced, and initial results of a sensitivity-enhanced version of the MM5 PSU/NCAR mesoscale weather model are presented. We also present a novel approach to the computation of gradients that uses a reverse mode approach at the time loop level and a forward mode approach at every time step. The resulting ``pseudoadjoint`` shares the characteristic of an adjoint code that the ratio of gradient to function evaluation does not depend on the number of independent variables. In contrast to a true adjoint approach, however, the nonlinearity of the model plays no role in the complexity of the derivative code.
Reconstruction of ocean circulation from sparse data using the adjoint method: LGM and the present
Kurahashi-Nakamura, T.; Losch, M. J.; Paul, A.; Mulitza, S.; Schulz, M.
2010-12-01
Understanding the behavior of the Earth's climate system under different conditions in the past is the basis for more robust projections of future climate. It is thought that the ocean circulation plays a very important role in the climate system, because it can greatly affect climate by dynamic-thermodynamic (as a medium of heat transport) and biogeochemical processes (by affecting the global carbon cycle). In this context, studying the period of the Last Glacial Maximum (LGM) is particularly promising, as it represents a climate state that is very different from today. Furthermore the LGM, compared to other paleoperiods, is characterized by a relatively good paleo-data coverage. Unfortunately, the ocean circulation during the LGM is still uncertain, with a range of climate models estimating both a stronger and a weaker formation rate of North Atlantic Deep Water (NADW) as compared to the present rate. Here, we present a project aiming at reducing this uncertainty by combining proxy data with a numerical ocean model using variational techniques. Our approach, the so-called adjoint method, employs a quadratic cost function of model-data differences weighted by their prior error estimates. We seek an optimal state estimate at the global minimum of the cost function by varying the independent control variables such as initial conditions (e.g. temperature), boundary conditions (e.g. surface winds, heat flux), or internal parameters (e.g. vertical diffusivity). The adjoint or dual model computes the gradient of the cost function with respect to these control variables and thus provides the information required by gradient descent algorithms. The gradients themselves provide valuable information about the sensitivity of the system to perturbations in the control variables. We use the Massachusetts Institute of Technology ocean general circulation model (MITgcm) with a cubed-sphere grid system that avoids converging grid lines and pole singularities. This model code is
Improving the Fit of a Land-Surface Model to Data Using its Adjoint
Raoult, Nina; Jupp, Tim; Cox, Peter; Luke, Catherine
2016-04-01
Land-surface models (LSMs) are crucial components of the Earth System Models (ESMs) which are used to make coupled climate-carbon cycle projections for the 21st century. The Joint UK Land Environment Simulator (JULES) is the land-surface model used in the climate and weather forecast models of the UK Met Office. In this study, JULES is automatically differentiated using commercial software from FastOpt, resulting in an analytical gradient, or adjoint, of the model. Using this adjoint, the adJULES parameter estimation system has been developed, to search for locally optimum parameter sets by calibrating against observations. We present an introduction to the adJULES system and demonstrate its ability to improve the model-data fit using eddy covariance measurements of gross primary production (GPP) and latent heat (LE) fluxes. adJULES also has the ability to calibrate over multiple sites simultaneously. This feature is used to define new optimised parameter values for the 5 Plant Functional Types (PFTS) in JULES. The optimised PFT-specific parameters improve the performance of JULES over 90% of the FLUXNET sites used in the study. These reductions in error are shown and compared to reductions found due to site-specific optimisations. Finally, we show that calculation of the 2nd derivative of JULES allows us to produce posterior probability density functions of the parameters and how knowledge of parameter values is constrained by observations.
Normal and adjoint integral and integrodifferential neutron transport equations. Pt. 2
Using the simplifying hypotheses of the integrodifferential Boltzmann equations of neutron transport, given in JEN 334 report, several integral equations, and theirs adjoint ones, are obtained. Relations between the different normal and adjoint eigenfunctions are established and, in particular, proceeding from the integrodifferential Boltzmann equation it's found out the relation between the solutions of the adjoint equation of its integral one, and the solutions of the integral equation of its adjoint one (author)
Searching for Standard Model Adjoint Scalars with Diboson Resonance Signatures
Carpenter, Linda M
2015-01-01
We explore the phenomenology of scalar fields in the adjoint representation of SM gauge groups. We write a general set of dimension 5 effective operators in which SM adjoint scalars couple to pairs of standard model bosons. Using these effective operators, we explore new possible decay channels of a scalar color octet into a gluon and a Z boson/ gluon and a photon. We recast several analyses from Run I of the LHC to find constraints on an a scalar octet decaying into these channels, and we project the discovery potential of color octets in our gluon+photon channel for the 14 TeV run of LHC.
Hoteit, Ibrahim
2010-03-02
An eddy-permitting adjoint-based assimilation system has been implemented to estimate the state of the tropical Pacific Ocean. The system uses the Massachusetts Institute of Technology\\'s general circulation model and its adjoint. The adjoint method is used to adjust the model to observations by controlling the initial temperature and salinity; temperature, salinity, and horizontal velocities at the open boundaries; and surface fluxes of momentum, heat, and freshwater. The model is constrained with most of the available data sets in the tropical Pacific, including Tropical Atmosphere and Ocean, ARGO, expendable bathythermograph, and satellite SST and sea surface height data, and climatologies. Results of hindcast experiments in 2000 suggest that the iterated adjoint-based descent is able to significantly improve the model consistency with the multivariate data sets, providing a dynamically consistent realization of the tropical Pacific circulation that generally matches the observations to within specified errors. The estimated model state is evaluated both by comparisons with observations and by checking the controls, the momentum balances, and the representation of small-scale features that were not well sampled by the observations used in the assimilation. As part of these checks, the estimated controls are smoothed and applied in independent model runs to check that small changes in the controls do not greatly change the model hindcast. This is a simple ensemble-based uncertainty analysis. In addition, the original and smoothed controls are applied to a version of the model with doubled horizontal resolution resulting in a broadly similar “downscaled” hindcast, showing that the adjustments are not tuned to a single configuration (meaning resolution, topography, and parameter settings). The time-evolving model state and the adjusted controls should be useful for analysis or to supply the forcing, initial, and boundary conditions for runs of other models.
Large-volume results in SU(2) with adjoint fermions
Del Debbio, Luigi; Lucini, Biagio; Pica, Claudio; Patella, Agostino; Rago, Antonio; Roman, Sabin
2013-01-01
Taming finite-volume effects is a crucial ingredient in order to identify the existence of IR fixed points. We present the latest results from our numerical simulations of SU(2) gauge theory with 2 Dirac fermions in the adjoint representation on large volumes. We compare with previous results, and...
Approximate nonlinear self-adjointness and approximate conservation laws
In this paper, approximate nonlinear self-adjointness for perturbed PDEs is introduced and its properties are studied. Consequently, approximate conservation laws which cannot be obtained by the approximate Noether theorem are constructed by means of the method. As an application, a class of perturbed nonlinear wave equations is considered to illustrate the effectiveness. (paper)
Assimilating Remote Ammonia Observations with a Refined Aerosol Thermodynamics Adjoint"
Ammonia emissions parameters in North America can be refined in order to improve the evaluation of modeled concentrations against observations. Here, we seek to do so by developing and applying the GEOS-Chem adjoint nested over North America to conductassimilation of observations...
On self-adjointness of singular Floquet Hamiltonians
Duclos, Pierre; Jensen, Arne
2010-01-01
Schrödinger equations with time-dependent interactions are studied. We investigate how to define the Floquet Hamiltonian as a self-adjoint operator, when the interaction is singular in time or space. Using these results we establish the existence of a bounded propagator, by applying a result given...
Adjoint electron-photon transport Monte Carlo calculations with ITS
A general adjoint coupled electron-photon Monte Carlo code for solving the Boltzmann-Fokker-Planck equation has recently been created. It is a modified version of ITS 3.0, a coupled electronphoton Monte Carlo code that has world-wide distribution. The applicability of the new code to radiation-interaction problems of the type found in space environments is demonstrated
The Adjoint of the CMAQ Aqueous Chemistry Module
Baek, J.; Stanier, C.; Saide, P.; Carmichael, G.; Henze, D.; Turner, M.; Zhao, S.; Hakami, A.; Resler, Jaroslav; Sandu, A.; Russell, A. P.; Jeong, G.; Nenes, A.; Capps, S.; Percell, P.; Pinder, R.; Napelenok, S.; Bash, J.; Chai, T.; Byun, D
Chapel Hill : CMAS, 2012. s. 91-92. [Annual CMAS Conference /11./. Chapel Hill, 15.10.2012-17.10.2012] Institutional support: RVO:67985807 Keywords : air pollution * adjoint * aqueous chemistry Subject RIV: DG - Athmosphere Sciences, Meteorology http://www.cmascenter.org/conference/2012/agenda.cfm
Non-self-adjoint hamiltonians defined by Riesz bases
Bagarello, F., E-mail: fabio.bagarello@unipa.it [Dipartimento di Energia, Ingegneria dell' Informazione e Modelli Matematici, Facoltà di Ingegneria, Università di Palermo, I-90128 Palermo, Italy and INFN, Università di Torino, Torino (Italy); Inoue, A., E-mail: a-inoue@fukuoka-u.ac.jp [Department of Applied Mathematics, Fukuoka University, Fukuoka 814-0180 (Japan); Trapani, C., E-mail: camillo.trapani@unipa.it [Dipartimento di Matematica e Informatica, Università di Palermo, I-90123 Palermo (Italy)
2014-03-15
We discuss some features of non-self-adjoint Hamiltonians with real discrete simple spectrum under the assumption that the eigenvectors form a Riesz basis of Hilbert space. Among other things, we give conditions under which these Hamiltonians can be factorized in terms of generalized lowering and raising operators.
High Order Adjoint Derivatives using ESDIRK Methods for Oil Reservoir Production Optimization
Capolei, Andrea; Stenby, Erling Halfdan; Jørgensen, John Bagterp
2012-01-01
continuous adjoints . The high order integration scheme allows larger time steps and therefore faster solution times. We compare gradient computation by the continuous adjoint method to the discrete adjoint method and the finite-difference method. The methods are implemented for a two phase flow reservoir...... simulator. Computational experiments demonstrate that the accuracy of the sensitivities obtained by the adjoint methods are comparable to the accuracy obtained by the finite difference method. The continuous adjoint method is able to use a different time grid than the forward integration. Therefore, it can...
The activation probe is a robust tool to measure flux of fusion products from a magnetically confined plasma. A carefully chosen solid sample is exposed to the flux, and the impinging ions transmute the material making it radioactive. Ultra-low level gamma-ray spectroscopy is used post mortem to measure the activity and, thus, the number of fusion products. This contribution presents the numerical analysis of the first measurement in the ASDEX Upgrade tokamak, which was also the first experiment to measure a single discharge. The ASCOT suite of codes was used to perform adjoint/reverse Monte Carlo calculations of the fusion products. The analysis facilitates, for the first time, a comparison of numerical and experimental values for absolutely calibrated flux. The results agree to within a factor of about two, which can be considered a quite good result considering the fact that all features of the plasma cannot be accounted in the simulations.Also an alternative to the present probe orientation was studied. The results suggest that a better optimized orientation could measure the flux from a significantly larger part of the plasma. A shorter version of this contribution is due to be published in PoS at: 1st EPS conference on Plasma Diagnostics
Adjoint Formulation for an Embedded-Boundary Cartesian Method
Nemec, Marian; Aftosmis, Michael J.; Murman, Scott M.; Pulliam, Thomas H.
2004-01-01
Many problems in aerodynamic design can be characterized by smooth and convex objective functions. This motivates the use of gradient-based algorithms, particularly for problems with a large number of design variables, to efficiently determine optimal shapes and configurations that maximize aerodynamic performance. Accurate and efficient computation of the gradient, however, remains a challenging task. In optimization problems where the number of design variables dominates the number of objectives and flow- dependent constraints, the cost of gradient computations can be significantly reduced by the use of the adjoint method. The problem of aerodynamic optimization using the adjoint method has been analyzed and validated for both structured and unstructured grids. The method has been applied to design problems governed by the potential, Euler, and Navier-Stokes equations and can be subdivided into the continuous and discrete formulations. Giles and Pierce provide a detailed review of both approaches. Most implementations rely on grid-perturbation or mapping procedures during the gradient computation that explicitly couple changes in the surface shape to the volume grid. The solution of the adjoint equation is usually accomplished using the same scheme that solves the governing flow equations. Examples of such code reuse include multistage Runge-Kutta schemes coupled with multigrid, approximate-factorization, line-implicit Gauss-Seidel, and also preconditioned GMRES. The development of the adjoint method for aerodynamic optimization problems on Cartesian grids has been limited. In contrast to implementations on structured and unstructured grids, Cartesian grid methods decouple the surface discretization from the volume grid. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin e t al. developed an adjoint formulation for the TRANAIR code
Äkäslompolo, Simppa; Tardini, Giovanni; Kurki-Suonio, Taina
2015-01-01
The activation probe is a robust tool to measure flux of fusion products from a magnetically confined plasma. A carefully chosen solid sample is exposed to the flux, and the impinging ions transmute the material makig it radioactive. Ultra-low level gamma-ray spectroscopy is used post mortem to measure the activity and, thus, the number of fusion products. This contribution presents the numerical analysis of the first measurement in the ASDEX Upgrade tokamak, which was also the first experiment to measure a single discharge. The ASCOT suite of codes was used to perform adjoint/reverse Monte-Carlo calculations of the fusion products. The analysis facilitated, for the first time, a comparison of numerical and experimental values for absolutely calibrated flux. The results agree to within 40%, which can be considered remarkable considering the fact that all features of the plasma cannot be accounted in the simulations. Also an alternative probe orientation was studied. The results suggest that a better optimized...
Accurate adjoint design sensitivities for nano metal optics.
Hansen, Paul; Hesselink, Lambertus
2015-09-01
We present a method for obtaining accurate numerical design sensitivities for metal-optical nanostructures. Adjoint design sensitivity analysis, long used in fluid mechanics and mechanical engineering for both optimization and structural analysis, is beginning to be used for nano-optics design, but it fails for sharp-cornered metal structures because the numerical error in electromagnetic simulations of metal structures is highest at sharp corners. These locations feature strong field enhancement and contribute strongly to design sensitivities. By using high-accuracy FEM calculations and rounding sharp features to a finite radius of curvature we obtain highly-accurate design sensitivities for 3D metal devices. To provide a bridge to the existing literature on adjoint methods in other fields, we derive the sensitivity equations for Maxwell's equations in the PDE framework widely used in fluid mechanics. PMID:26368483
Refined topological vertex, cylindric partitions and U(1) adjoint theory
We study the partition function of the compactified 5D U(1) gauge theory (in the Ω-background) with a single adjoint hypermultiplet, calculated using the refined topological vertex. We show that this partition function is an example a periodic Schur process and is a refinement of the generating function of cylindric plane partitions. The size of the cylinder is given by the mass of adjoint hypermultiplet and the parameters of the Ω-background. We also show that this partition function can be written as a trace of operators which are generalizations of vertex operators studied by Carlsson and Okounkov. In the last part of the paper we describe a way to obtain (q,t) identities using the refined topological vertex.
Adjoint Fokker-Planck equation and runaway electron dynamics
Liu, Chang; Brennan, Dylan P.; Bhattacharjee, Amitava [Princeton University, Princeton, New Jersey 08544 (United States); Boozer, Allen H. [Columbia University, New York, New York 10027 (United States)
2016-01-15
The adjoint Fokker-Planck equation method is applied to study the runaway probability function and the expected slowing-down time for highly relativistic runaway electrons, including the loss of energy due to synchrotron radiation. In direct correspondence to Monte Carlo simulation methods, the runaway probability function has a smooth transition across the runaway separatrix, which can be attributed to effect of the pitch angle scattering term in the kinetic equation. However, for the same numerical accuracy, the adjoint method is more efficient than the Monte Carlo method. The expected slowing-down time gives a novel method to estimate the runaway current decay time in experiments. A new result from this work is that the decay rate of high energy electrons is very slow when E is close to the critical electric field. This effect contributes further to a hysteresis previously found in the runaway electron population.
Adjoint-based sensitivity analysis for reactor accident codes
This paper summarizes a recently completed study that identified and investigated the difficulties and limitations of applying first-order adjoint sensitivity methods to reactor accident codes. The work extends earlier adjoint sensitivity formulations and applications to consider problem/model discontinuities in a general fashion, provide for response (R) formulations required by reactor safety applications, and provide a scheme for accurately handling extremely time-sensitive reactor accident responses. The scheme involves partitioning (dividing) the model into submodels (with spearate defining equations and initial conditions) at the location of discontinuity. Successful partitioning moves the problem dependence on the discontinuity location from the whole model system equations to the initial conditions of the second submodel
Examination of Observation Impacts derived from OSEs and Adjoint Models
Gelaro, Ronald
2008-01-01
With the adjoint of a data assimilation system, the impact of any or all assimilated observations on measures of forecast skill can be estimated accurately and efficiently. The approach allows aggregation of results in terms of individual data types, channels or locations, all computed simultaneously. In this study, adjoint-based estimates of observation impact are compared with results from standard observing system experiments (OSEs) in the NASA Goddard Earth Observing System Model, Version 5 (GEOS-5) GEOS-5 system. The two approaches are shown to provide unique, but complimentary, information. Used together, they reveal both redundancies and dependencies between observing system impacts as observations are added or removed. Understanding these dependencies poses a major challenge for optimizing the use of the current observational network and defining requirements for future observing systems.
Adjoint equation of ADS sub-critical reactor
Compared with the critical reactor, the distributions of source neutron and fission neutron are asymmetric inside ADS (accelerator driven sub-critical system) sub-critical reactor, as well as the importance function is different. The multigroup-diffusion approximation was used to simplify the steady-state transport equation into multigroup equation. Then an adjoint equation normalized by the power of reactor core and an importance function associated with the relative power were derived. The physical significance of neutron importance in the sub-critical reactor was also derived. Finally, two different expressions of multiplication factor for sub-critical reactor with external neutron source were derived based on steady-state adjoint equations. (authors)
A comparison of adjoint and data-centric verification techniques.
Wildey, Timothy Michael; Cyr, Eric Christopher; Shadid, John Nicolas; Pawlowski, Roger Patrick; Smith, Thomas Michael
2013-03-01
This document summarizes the results from a level 3 milestone study within the CASL VUQ effort. We compare the adjoint-based a posteriori error estimation approach with a recent variant of a data-centric verification technique. We provide a brief overview of each technique and then we discuss their relative advantages and disadvantages. We use Drekar::CFD to produce numerical results for steady-state Navier Stokes and SARANS approximations. 3
Optimization of a neutron detector design using adjoint transport simulation
Yi, C.; Manalo, K.; Huang, M.; Chin, M.; Edgar, C.; Applegate, S.; Sjoden, G. [Georgia Inst. of Technology, Gilhouse Boggs Bldg., 770 State St, Atlanta, GA 30332-0745 (United States)
2012-07-01
A synthetic aperture approach has been developed and investigated for Special Nuclear Materials (SNM) detection in vehicles passing a checkpoint at highway speeds. SNM is postulated to be stored in a moving vehicle and detector assemblies are placed on the road-side or in chambers embedded below the road surface. Neutron and gamma spectral awareness is important for the detector assembly design besides high efficiencies, so that different SNMs can be detected and identified with various possible shielding settings. The detector assembly design is composed of a CsI gamma-ray detector block and five neutron detector blocks, with peak efficiencies targeting different energy ranges determined by adjoint simulations. In this study, formulations are derived using adjoint transport simulations to estimate detector efficiencies. The formulations is applied to investigate several neutron detector designs for Block IV, which has its peak efficiency in the thermal range, and Block V, designed to maximize the total neutron counts over the entire energy spectrum. Other Blocks detect different neutron energies. All five neutron detector blocks and the gamma-ray block are assembled in both MCNP and deterministic simulation models, with detector responses calculated to validate the fully assembled design using a 30-group library. The simulation results show that the 30-group library, collapsed from an 80-group library using an adjoint-weighting approach with the YGROUP code, significantly reduced the computational cost while maintaining accuracy. (authors)
Consistent Adjoint Driven Importance Sampling using Space, Energy and Angle
Peplow, Douglas E. [ORNL; Mosher, Scott W [ORNL; Evans, Thomas M [ORNL
2012-08-01
For challenging radiation transport problems, hybrid methods combine the accuracy of Monte Carlo methods with the global information present in deterministic methods. One of the most successful hybrid methods is CADIS Consistent Adjoint Driven Importance Sampling. This method uses a deterministic adjoint solution to construct a biased source distribution and consistent weight windows to optimize a specific tally in a Monte Carlo calculation. The method has been implemented into transport codes using just the spatial and energy information from the deterministic adjoint and has been used in many applications to compute tallies with much higher figures-of-merit than analog calculations. CADIS also outperforms user-supplied importance values, which usually take long periods of user time to develop. This work extends CADIS to develop weight windows that are a function of the position, energy, and direction of the Monte Carlo particle. Two types of consistent source biasing are presented: one method that biases the source in space and energy while preserving the original directional distribution and one method that biases the source in space, energy, and direction. Seven simple example problems are presented which compare the use of the standard space/energy CADIS with the new space/energy/angle treatments.
Unsteady adjoint of pressure loss for a fundamental transonic turbine vane
Talnikar, Chaitanya; Laskowski, Gregory M
2015-01-01
High fidelity simulations, e.g., large eddy simulation are often needed for accurately predicting pressure losses due to wake mixing in turbomachinery applications. An unsteady adjoint of such high fidelity simulations is useful for design optimization in these aerodynamic applications. In this paper we present unsteady adjoint solutions using a large eddy simulation model for a vane from VKI using aerothermal objectives. The unsteady adjoint method is effective in capturing the gradient for a short time interval aerothermal objective, whereas the method provides diverging gradients for long time-averaged thermal objectives. As the boundary layer on the suction side near the trailing edge of the vane is turbulent, it poses a challenge for the adjoint solver. The chaotic dynamics cause the adjoint solution to diverge exponentially from the trailing edge region when solved backwards in time. This results in the corruption of the sensitivities obtained from the adjoint solutions. An energy analysis of the unstea...
Adjoint-Based Sensitivity Maps for the Nearshore
Orzech, Mark; Veeramony, Jay; Ngodock, Hans
2013-04-01
The wave model SWAN (Booij et al., 1999) solves the spectral action balance equation to produce nearshore wave forecasts and climatologies. It is widely used by the coastal modeling community and is part of a variety of coupled ocean-wave-atmosphere model systems. A variational data assimilation system (Orzech et al., 2013) has recently been developed for SWAN and is presently being transitioned to operational use by the U.S. Naval Oceanographic Office. This system is built around a numerical adjoint to the fully nonlinear, nonstationary SWAN code. When provided with measured or artificial "observed" spectral wave data at a location of interest on a given nearshore bathymetry, the adjoint can compute the degree to which spectral energy levels at other locations are correlated with - or "sensitive" to - variations in the observed spectrum. Adjoint output may be used to construct a sensitivity map for the entire domain, tracking correlations of spectral energy throughout the grid. When access is denied to the actual locations of interest, sensitivity maps can be used to determine optimal alternate locations for data collection by identifying regions of greatest sensitivity in the mapped domain. The present study investigates the properties of adjoint-generated sensitivity maps for nearshore wave spectra. The adjoint and forward SWAN models are first used in an idealized test case at Duck, NC, USA, to demonstrate the system's effectiveness at optimizing forecasts of shallow water wave spectra for an inaccessible surf-zone location. Then a series of simulations is conducted for a variety of different initializing conditions, to examine the effects of seasonal changes in wave climate, errors in bathymetry, and variations in size and shape of the inaccessible region of interest. Model skill is quantified using two methods: (1) a more traditional correlation of observed and modeled spectral statistics such as significant wave height, and (2) a recently developed RMS
Inflow and initial conditions for direct numerical simulation based on adjoint data assimilation
Gronskis, A; Heitz, D.; Mémin, E.
2011-01-01
A method for generating inﬂow conditions for direct numerical simulations (DNS) of spatially-developing ﬂows is presented. The proposed method is based on variational data assimilation and adjoint-based optimization. The estimation is conducted through an iterative process involving a forward integration of a given dynamical model followed by a backward integration of an adjoint system deﬁned by the adjoint of the discrete scheme associated to the dynamical system. The ap...
Self-Adjoint Extensions of the Pauli Equation in the Presence of a Magnetic Monopole
Karat, Edwin R.; Schulz, Michael B.
1996-01-01
We discuss the Hamiltonian for a nonrelativistic electron with spin in the presence of an abelian magnetic monopole and note that it is not self-adjoint in the lowest two angular momentum modes. We then use von Neumann's theory of self-adjoint extensions to construct a self-adjoint operator with the same functional form. In general, this operator will have eigenstates in which the lowest two angular momentum modes mix, thereby removing conservation of angular momentum. However, consistency wi...
Integration of the adjoint gamma quantum transport equation by the Monte Carlo method
Comparative description and analysis of the direct and adjoint algorithms of calculation of gamma-quantum transmission in shielding using the Monte Carlo method have been carried out. Adjoint estimations for a number of monoenergetic sources have been considered. A brief description of ''COMETA'' program for BESM-6 computer reazaling direct and adjoint algorithms is presented. The program is modular-constructed which allows to widen it the new module-units being joined. Results of solution by the adjoint branch of two analog problems as compared to the analytical data are presented. These results confirm high efficiency of ''COMETA'' program
Adjoint optimization of natural convection problems: differentially heated cavity
Saglietti, Clio; Schlatter, Philipp; Monokrousos, Antonios; Henningson, Dan S.
2016-06-01
Optimization of natural convection-driven flows may provide significant improvements to the performance of cooling devices, but a theoretical investigation of such flows has been rarely done. The present paper illustrates an efficient gradient-based optimization method for analyzing such systems. We consider numerically the natural convection-driven flow in a differentially heated cavity with three Prandtl numbers (Pr=0.15{-}7 ) at super-critical conditions. All results and implementations were done with the spectral element code Nek5000. The flow is analyzed using linear direct and adjoint computations about a nonlinear base flow, extracting in particular optimal initial conditions using power iteration and the solution of the full adjoint direct eigenproblem. The cost function for both temperature and velocity is based on the kinetic energy and the concept of entransy, which yields a quadratic functional. Results are presented as a function of Prandtl number, time horizons and weights between kinetic energy and entransy. In particular, it is shown that the maximum transient growth is achieved at time horizons on the order of 5 time units for all cases, whereas for larger time horizons the adjoint mode is recovered as optimal initial condition. For smaller time horizons, the influence of the weights leads either to a concentric temperature distribution or to an initial condition pattern that opposes the mean shear and grows according to the Orr mechanism. For specific cases, it could also been shown that the computation of optimal initial conditions leads to a degenerate problem, with a potential loss of symmetry. In these situations, it turns out that any initial condition lying in a specific span of the eigenfunctions will yield exactly the same transient amplification. As a consequence, the power iteration converges very slowly and fails to extract all possible optimal initial conditions. According to the authors' knowledge, this behavior is illustrated here
Adjoint Monte Carlo techniques and codes for organ dose calculations
Adjoint Monte Carlo simulations can be effectively used for the estimation of doses in small targets when the sources are extended in large volumes or surfaces. The main features of two computer codes for calculating doses at free points or in organs of an anthropomorphic phantom are described. In the first program (REBEL-3) natural gamma-emitting sources are contained in the walls of a dwelling room; in the second one (POKER-CAMP) the user can specify arbitrary gamma sources with different spatial distributions in the environment: in (or on the surface of) the ground and in the air. 3 figures
Bimetric Gravity From Adjoint Frame Field In Four Dimensions
Guo, Zhi-Qiang
2015-01-01
We provide a novel model of gravity by using adjoint frame fields in four dimensions. It has a natural interpretation as a gravitational theory of a complex metric field, which describes interactions between two real metrics. The classical solutions establish three appealing features. The spherical symmetric black hole solution has an additional hair, which includes the Schwarzschild solution as a special case. The de Sitter solution is realized without introducing a cosmological constant. The constant flat background breaks the Lorentz invariance spontaneously, although the Lorentz breaking effect can be localized to the second metric while the first metric still respects the Lorentz invariance.
Running coupling in SU(2) with two adjoint fermions
Rantaharju, Jarno; Rummukainen, Kari; Tuominen, Kimmo
2015-01-01
We study SU(2) gauge theory with two Dirac fermions in the adjoint representation of the gauge group on the lattice. Using clover improved Wilson fermion action with hypercubic truncated stout smearing we perform simulations at larger coupling than earlier. We measure the evolution of the coupling constant using the Schr\\"odinger functional method. Extrapolating our lattice results to the continuum, we confirm the existence of a fixed point in the interval 2.2 $\\lesssim $ $g^{*2}$ $\\lesssim $ 3. We also measure the anomalous dimension and find its value at the fixed point is 0.18 $\\lesssim $ $\\gamma^*$ $\\lesssim $ 0.23.
Self-adjoint dirac equation for a fuzzy potential representation
In the present paper, the dependent-, as well as the independent-variables that belong to a self-adjointed relativistic scattering equation, are properly transformed, for a fuzzy potential representation of Woods and Saxon type, with the purpose of being amenable to analytic solutions. This has been accomplished, in terms of hypergeometric functions that can be analytically continued before-and after-the boundary separating the interior from the exterior. Studies of these analytic solutions about the singular points are performed and the corresponding asymptotic behavior is investigated. The corresponding scattering matrix is extracted, where the resonance energy Eigen-values can be identified
Large-N reduction with adjoint Wilson fermions
Bringoltz, Barak; Sharpe, Stephen R
2012-01-01
We analyze the large-N behavior of SU(N) lattice gauge theories with adjoint fermions by studying volume-reduced models, as pioneered by Eguchi and Kawai. We perform simulations on a single-site lattice for Nf = 1 and Nf = 2 Wilson Dirac fermions with values of N up to 53. We show for both values of Nf that in the large-N limit there is a finite region, containing both light and heavy fermions, of unbroken center symmetry where the theory exhibits volume independence. Using large-N reduction we attempt to calculate physical quantities such as the string tension and meson masses.
Inverse Modeling of Emissions using the CMAQ Adjoint Model
Resler, Jaroslav; Eben, Kryštof; Juruš, Pavel; Krč, Pavel
Chapel Hill : CMAS, 2008, s. 1-5. [Annual CMAS Conference /7./. Chapel Hill (US), 06.10.2008-08.10.2008] R&D Projects: GA AV ČR 1ET400300414; GA MŽP SP/1A4/107/07 Institutional research plan: CEZ:AV0Z10300504 Keywords : 4DVar * data assimilation * inverse modelling * emission * CMAQ adjoint * tropospheric column of NO2 * satellite instruments * GOME2 * OMI Subject RIV: IN - Informatics, Computer Science http://www.cmascenter.org/conference/2008/agenda.cfm
On the Norm Convergence of the Self-Adjoint Trotter–Kato Product Formula with Error Bound
Takashi Ichinose; Hideo Tamura
2002-02-01
The norm convergence of the Trotter–Kato product formula with error bound is shown for the semigroup generated by that operator sum of two nonnegative self-adjoint operators and which is self-adjoint.
On rational R-matrices with adjoint SU(n) symmetry
Stronks, Laurens; Schuricht, Dirk
2016-01-01
Using the representation theory of Yangians we construct the rational R-matrix which takes values in the adjoint representation of SU(n). From this we derive an integrable SU(n) spin chain with lattice spins transforming under the adjoint representation. However, the resulting Hamiltonian is found to be non-Hermitian.
An efficient adjoint design sensitivity analysis method is developed for reduced atomic systems. A reduced atomic system and the adjoint system are constructed in a locally confined region, utilizing generalized Langevin equation (GLE) for periodic lattice structures. Due to the translational symmetry of lattice structures, the size of time history kernel function that accounts for the boundary effects of the reduced atomic systems could be reduced to a single atom’s degrees of freedom. For the problems of highly nonlinear design variables, the finite difference method is impractical for its inefficiency and inaccuracy. However, the adjoint method is very efficient regardless of the number of design variables since one additional time integration is required for the adjoint GLE. Through numerical examples, the derived adjoint sensitivity turns out to be accurate and efficient through the comparison with finite difference sensitivity
Singular Potentials in Quantum Mechanics and Ambiguity in the Self-Adjoint Hamiltonian
Tamás Fülöp
2007-11-01
Full Text Available For a class of singular potentials, including the Coulomb potential (in three and less dimensions and $V(x = g/x^2$ with the coefficient $g$ in a certain range ($x$ being a space coordinate in one or more dimensions, the corresponding Schrödinger operator is not automatically self-adjoint on its natural domain. Such operators admit more than one self-adjoint domain, and the spectrum and all physical consequences depend seriously on the self-adjoint version chosen. The article discusses how the self-adjoint domains can be identified in terms of a boundary condition for the asymptotic behaviour of the wave functions around the singularity, and what physical differences emerge for different self-adjoint versions of the Hamiltonian. The paper reviews and interprets known results, with the intention to provide a practical guide for all those interested in how to approach these ambiguous situations.
Solar wind reconstruction from magnetosheath data using an adjoint approach
Nabert, C.; Othmer, C. [Technische Univ. Braunschweig (Germany). Inst. fuer Geophysik und extraterrestrische Physik; Glassmeier, K.H. [Technische Univ. Braunschweig (Germany). Inst. fuer Geophysik und extraterrestrische Physik; Max Planck Institute for Solar System Research, Goettingen (Germany)
2015-07-01
We present a new method to reconstruct solar wind conditions from spacecraft data taken during magnetosheath passages, which can be used to support, e.g., magnetospheric models. The unknown parameters of the solar wind are used as boundary conditions of an MHD (magnetohydrodynamics) magnetosheath model. The boundary conditions are varied until the spacecraft data matches the model predictions. The matching process is performed using a gradient-based minimization of the misfit between data and model. To achieve this time-consuming procedure, we introduce the adjoint of the magnetosheath model, which allows efficient calculation of the gradients. An automatic differentiation tool is used to generate the adjoint source code of the model. The reconstruction method is applied to THEMIS (Time History of Events and Macroscale Interactions during Substorms) data to calculate the solar wind conditions during spacecraft magnetosheath transitions. The results are compared to actual solar wind data. This allows validation of our reconstruction method and indicates the limitations of the MHD magnetosheath model used.
Solar wind reconstruction from magnetosheath data using an adjoint approach
We present a new method to reconstruct solar wind conditions from spacecraft data taken during magnetosheath passages, which can be used to support, e.g., magnetospheric models. The unknown parameters of the solar wind are used as boundary conditions of an MHD (magnetohydrodynamics) magnetosheath model. The boundary conditions are varied until the spacecraft data matches the model predictions. The matching process is performed using a gradient-based minimization of the misfit between data and model. To achieve this time-consuming procedure, we introduce the adjoint of the magnetosheath model, which allows efficient calculation of the gradients. An automatic differentiation tool is used to generate the adjoint source code of the model. The reconstruction method is applied to THEMIS (Time History of Events and Macroscale Interactions during Substorms) data to calculate the solar wind conditions during spacecraft magnetosheath transitions. The results are compared to actual solar wind data. This allows validation of our reconstruction method and indicates the limitations of the MHD magnetosheath model used.
Adjoint $SU(5)$ GUT model with $T_{7}$ flavor symmetry
Arbeláez, Carolina; Kovalenko, Sergey; Schmidt, Iván
2015-01-01
We propose an adjoint $SU(5)$ GUT model with a $T_{7}$ family symmetry and an extra $Z_{2}\\otimes Z_{2}^{\\prime }\\otimes Z_{3}\\otimes Z_{4}\\otimes Z_{12}$ discrete group, that successfully describes the prevailing Standard Model (SM) fermion mass and mixing pattern. The observed hierarchy of the charged fermion masses and the quark mixing angles arises from the $Z_{3}\\otimes Z_{4}\\otimes Z_{12}$ symmetry breaking, which occurs near the GUT scale. The light active neutrino masses are generated by type I and type III seesaw mechanisms mediated by the fermionic $SU(5)$ singlet and the adjoint $\\mathbf{24}$-plet. The model predicts the effective Majorana neutrino mass parameter of neutrinoless double beta decay to be $m_{\\beta \\beta }=$ 4 and 50 meV for the normal and the inverted neutrino spectrum, respectively. We construct several benchmark scenarios, which lead to $SU(5)$ gauge coupling unification and are compatible with the known phenomenological constraints originating from the lightness of neutrinos, prot...
Limitations of Adjoint-Based Optimization for Separated Flows
Otero, J. Javier; Sharma, Ati; Sandberg, Richard
2015-11-01
Cabin noise is generated by the transmission of turbulent pressure fluctuations through a vibrating panel and can lead to fatigue. In the present study, we model this problem by using DNS to simulate the flow separating off a backward facing step and interacting with a plate downstream of the step. An adjoint formulation of the full compressible Navier-Stokes equations with varying viscosity is used to calculate the optimal control required to minimize the fluid-structure-acoustic interaction with the plate. To achieve noise reduction, a cost function in wavenumber space is chosen to minimize the excitation of the lower structural modes of the structure. To ensure the validity of time-averaged cost functions, it is essential that the time horizon is long enough to be a representative sample of the statistical behaviour of the flow field. The results from the current study show how this scenario is not always feasible for separated flows, because the chaotic behaviour of turbulence surpasses the ability of adjoint-based methods to compute time-dependent sensitivities of the flow.
On the Self-adjointness of the Product Operators of Two mth-Order Differential Operators on [0, +∞)
Jian Ye AN; Jiong SUN
2004-01-01
In the present paper, the self-adjointness of the product of two mth-order differential operators on [0, +∞) is studied. By means of the construction theory of self-adjoint operators and matrix computation, we obtain a sufficient and necessary condition to ensure that the product operator is self-adjoint, which extends the results in the second order case.
Adjoint optimization scheme for lower hybrid current rampup and profile control in Tokamak
The purpose of this work is to take into account and study the effect of the electric field profiles on the Lower Hybrid (LH) current drive efficiency during transient phases such as rampup. As a complement to the full ray-tracing / Fokker Planck studies, and for the purpose of optimization studies, we developed a simplified 1-D model based on the adjoint Karney-Fisch numerical results. This approach allows us to estimate the LH power deposition profile which would be required for ramping the current with prescribed rate, total current density profile (q-profile) and surface loop voltage. For rampup optimization studies, we can therefore scan the whole parameter space and eliminate a posteriori those scenarios which correspond to unrealistic deposition profiles. We thus obtain the time evolution of the LH power, minor radius of the plasma, volt-second consumption and total energy dissipated. Optimization can thus be performed with respect to any of those criteria. This scheme is illustrated by some numerical simulations performed with TORE-SUPRA and NET/ITER parameters. We conclude with a derivation of a simple and general scaling law for the flux consumption during the rampup phase
Conformal versus confining scenario in SU(2) with adjoint fermions
The masses of the lowest-lying states in the meson and in the gluonic sector of an SU(2) gauge theory with two Dirac flavors in the adjoint representation are measured on the lattice at a fixed value of the lattice coupling β=4/g02=2.25 for values of the bare fermion mass m0 that span a range between the quenched regime and the massless limit, and for various lattice volumes. Even for light constituent fermions the lightest glueballs are found to be lighter than the lightest mesons. Moreover, the string tension between two static fundamental sources strongly depends on the mass of the dynamical fermions and becomes of the order of the inverse squared lattice linear size before the chiral limit is reached. The implications of these findings for the phase of the theory in the massless limit are discussed and a strategy for discriminating between the (near-)conformal and the confining scenario is outlined.
Infrared regime of SU(2) with one adjoint Dirac flavor
Athenodorou, Andreas; Bennett, Ed; Bergner, Georg; Lucini, Biagio
2015-06-01
SU(2) gauge theory with one Dirac flavor in the adjoint representation is investigated on a lattice. Initial results for the gluonic and mesonic spectrum, static potential from Wilson and Polyakov loops, and the anomalous dimension of the fermionic condensate from the Dirac mode number are presented. The results found are not consistent with conventional confining behavior, pointing instead tentatively towards a theory lying within or very near the onset of the conformal window, with the anomalous dimension of the fermionic condensate in the range 0.9 ≲γ*≲0.95 . The implications of our work for building a viable theory of strongly interacting dynamics beyond the standard model are discussed.
Adjoint Fokker-Planck equation and runaway electron dynamics
Liu, Chang; Boozer, Allen H; Bhattacharjee, Amitava
2016-01-01
A new method to obtain the runaway probability and the expected slowing-down time for runaway electrons is developed, by solving the adjoint Fokker-Planck equation in momentum space. The runaway probability function has a smooth transition at the runaway separatrix, which can be attributed to the effect of the pitch angle scattering term in the kinetic equation. The expected slowing-down time gives a new way to estimate the runaway current decay time in experiments. The result shows that the decay rate of high energy electron is very slow when E is close to the critical electric field, which helps elucidate the hysteresis effect seen in the runaway electron population. Given the same numerical accuracy, the new method is more efficient than the Monte Carlo simulation.
An Adjoint-Based Adaptive Ensemble Kalman Filter
Song, Hajoon
2013-10-01
A new hybrid ensemble Kalman filter/four-dimensional variational data assimilation (EnKF/4D-VAR) approach is introduced to mitigate background covariance limitations in the EnKF. The work is based on the adaptive EnKF (AEnKF) method, which bears a strong resemblance to the hybrid EnKF/three-dimensional variational data assimilation (3D-VAR) method. In the AEnKF, the representativeness of the EnKF ensemble is regularly enhanced with new members generated after back projection of the EnKF analysis residuals to state space using a 3D-VAR [or optimal interpolation (OI)] scheme with a preselected background covariance matrix. The idea here is to reformulate the transformation of the residuals as a 4D-VAR problem, constraining the new member with model dynamics and the previous observations. This should provide more information for the estimation of the new member and reduce dependence of the AEnKF on the assumed stationary background covariance matrix. This is done by integrating the analysis residuals backward in time with the adjoint model. Numerical experiments are performed with the Lorenz-96 model under different scenarios to test the new approach and to evaluate its performance with respect to the EnKF and the hybrid EnKF/3D-VAR. The new method leads to the least root-mean-square estimation errors as long as the linear assumption guaranteeing the stability of the adjoint model holds. It is also found to be less sensitive to choices of the assimilation system inputs and parameters.
Use of adjoint methods in the probabilistic finite element approach to fracture mechanics
Liu, Wing Kam; Besterfield, Glen; Lawrence, Mark; Belytschko, Ted
1988-01-01
The adjoint method approach to probabilistic finite element methods (PFEM) is presented. When the number of objective functions is small compared to the number of random variables, the adjoint method is far superior to the direct method in evaluating the objective function derivatives with respect to the random variables. The PFEM is extended to probabilistic fracture mechanics (PFM) using an element which has the near crack-tip singular strain field embedded. Since only two objective functions (i.e., mode I and II stress intensity factors) are needed for PFM, the adjoint method is well suited.
Freire, Igor Leite
2010-01-01
In this work we consider the problem on group classification and conservation laws of the general first order evolution equations. We obtain the subclasses of these general equations which are quasi-self-adjoint and self-adjoint. By using the recent Ibragimov's Theorem on conservation laws, we establish the conservation laws of the equations admiting self-adjoint equations. We illustrate our results applying them to the inviscid Burgers' equation. In particular an infinite number of new symmetries of these equations are found and their corresponding conservation laws are established.
Quantum cosmology of scalar-tensor theories and self-adjointness
Almeida, C R; Fabris, J C; Moniz, P V
2016-01-01
In this paper, the problem of the self-adjointness for the case of a quantum minisuperspace Hamiltonian retrieved from a Brans-Dicke (BD) action is investigated. Our matter content is presented in terms of a perfect fluid, onto which the Schutz's formalism will be applied. We use the von Neumann theorem and the similarity with the Laplacian operator in one of the variables to determine the cases where the Hamiltonian is self-adjoint and if it admits self-adjoint extensions. For the latter, we study which extension is physically more suitable.
Adjoint Assimilation in Marine Ecosystem Models and an Example of Application
XU Qing; LIU Yuguang; L(U) Xianqing
2005-01-01
This paper aims at a review of the work carried out to date on the adjoint assimilation of data in marine ecosystem models since 1995. The structure and feature of the adjoint assimilation in marine ecosystem models are also introduced.To illustrate the application of the adjoint technique and its merits, a 4-variable ecosystem model coupled with a 3-D physical model is established for the Bohai Sea and the Yellow Sea. The chlorophyll concentration data derived from the SeaWiFS ocean colour data are assimilated in the model with the technique. Some results are briefly presented.
Revisit boundary conditions for the self-adjoint angular flux formulation
Wang, Yaqi [Idaho National Lab. (INL), Idaho Falls, ID (United States); Gleicher, Frederick N. [Idaho National Lab. (INL), Idaho Falls, ID (United States)
2015-03-01
We revisit the boundary conditions for SAAF. We derived the equivalent parity variational form ready for coding up. The more rigorous approach of evaluating odd parity should be solving the odd parity equation coupled with the even parity. We proposed a symmetric reflecting boundary condition although neither positive definiteness nor even-odd decoupling is achieved. A simple numerical test verifies the validity of these boundary conditions.
MS S4.03.002 - Adjoint-Based Design for Configuration Shaping
Nemec, Marian; Aftosmis, Michael J.
2009-01-01
This slide presentation discusses a method of inverse design for low sonic boom using adjoint-based gradient computations. It outlines a method for shaping a configuration in order to match a prescribed near-field signature.
Self-adjoint extensions of the Pauli equation in the presence of a magnetic monopole
Karat, E R; Karat, Edwin R; Schulz, Michael B
1996-01-01
We discuss the Hamiltonian for a nonrelativistic electron with spin in the presence of a magnetic monopole and note that it is not self-adjoint in the lowest two angular momentum modes. We then use von Neumann's theory of self-adjoint extensions to construct a self-adjoint operator with the same functional form. In general, this operator will have eigenstates in which the lowest two angular momentum modes mix, thereby removing conservation of angular momentum. Because the same effect occurs for a spinless particle with a sufficiently attractive inverse square potential, we also study this system. We use this simpler Hamiltonian to compare the eigenfunctions corresponding to a particular self-adjoint extension with the eigenfunctions satisfying a boundary condition consistent with probability conservation.
A SYSTEMATIC FORMULATION OF THE CONTINUOUS ADJOINT METHOD APPLIED TO VISCOUS AERODYNAMIC DESIGN
C. Castro*, C. Lozano**, F. Palacios*** and E. Zuazua****
2009-01-01
Full Text Available A continuous adjoint approach to aerodynamic design for viscous compressible flows on unstructuredgrids is developed, and three important problems raised in the continuous adjoint literature are solved. First, using tools of shape deformation of boundary integrals a generic adjoint formulation is developed withindependence of the kind of mesh used. Then, a systematic way of reducing the 2nd order derivative terms which arise is presented which avoids the need of using higher order numerical solvers to obtain accurateapproximations of the 2nd order derivatives. And finally, the class of admissible optimization functionals isclarified. Several remarks are made concerning the longstanding discrete vs. continuous adjoint dichotomy, with the emphasis not on the advantages or disadvantages of each method, but rather on the well-posedness of the approaches. The accuracy of the sensitivity derivatives is assessed by comparison with finite-difference computations, and the validity of the overall methodology is illustrated with design examples under demanding subsonic conditions.
Self-adjoint extensions of the Pauli equation in the presence of a magnetic monopole
We discuss the Hamiltonian for a nonrelativistic electron with spin in the presence of an abelian magnetic monopole and note that it is not self-adjoint in the lowest two angular momentum modes. We then use von Neumann close-quote s theory of self-adjoint extensions to construct a self-adjoint operator with the same functional form. In general, this operator will have eigenstates in which the lowest two angular momentum modes mix, thereby removing conservation of angular momentum. However, consistency with the solutions of the Dirac equation limits the possibilities such that conservation of angular momentum is restored. Because the same effect occurs for a spinless particle with a sufficiently attractive inverse square potential, we also study this system. We use this simpler Hamiltonian to compare the eigenfunctions corresponding to a particular self-adjoint extension with the eigenfunctions satisfying a boundary condition consistent with probability conservation. copyright 1997 Academic Press, Inc
Self-adjoint extensions of the Pauli equation in the presence of a magnetic monopole
Karat, E.; Schulz, M. [Center for Theoretical Physics, Laboratory for Nuclear Science, and Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139-4307 (United States)
1997-02-01
We discuss the Hamiltonian for a nonrelativistic electron with spin in the presence of an abelian magnetic monopole and note that it is not self-adjoint in the lowest two angular momentum modes. We then use von Neumann{close_quote}s theory of self-adjoint extensions to construct a self-adjoint operator with the same functional form. In general, this operator will have eigenstates in which the lowest two angular momentum modes mix, thereby removing conservation of angular momentum. However, consistency with the solutions of the Dirac equation limits the possibilities such that conservation of angular momentum is restored. Because the same effect occurs for a spinless particle with a sufficiently attractive inverse square potential, we also study this system. We use this simpler Hamiltonian to compare the eigenfunctions corresponding to a particular self-adjoint extension with the eigenfunctions satisfying a boundary condition consistent with probability conservation. {copyright} 1997 Academic Press, Inc.
Variational characterizations for eigenfunctions of analytic self-adjoint operator functions
Georgios Katsouleas; John Maroulas
2013-01-01
In this paper we consider Rellich's diagonalization theorem for analytic self-adjoint operator functions and investigate variational principles for their eigenfunctions and interlacing statements. As an application, we present a characterization for the eigenvalues of hyperbolic operator polynomials.
Hermand, Jean-Pierre; Berrada, Mohamed; Meyer, Matthias; Asch, Mark
2005-09-01
Recently, an analytic adjoint-based method of optimal nonlocal boundary control has been proposed for inversion of a waveguide acoustic field using the wide-angle parabolic equation [Meyer and Hermand, J. Acoust. Soc. Am. 117, 2937-2948 (2005)]. In this paper a numerical extension of this approach is presented that allows the direct inversion for the geoacoustic parameters which are embedded in a spectral integral representation of the nonlocal boundary condition. The adjoint model is generated numerically and the inversion is carried out jointly across multiple frequencies. The paper further discusses the application of the numerical adjoint PE method for ocean acoustic tomography. To show the effectiveness of the implemented numerical adjoint, preliminary inversion results of water sound-speed profile and bottom acoustic properties will be shown for the YELLOW SHARK '94 experimental conditions.
Yaremchuk, Max; Martin, Paul; Koch, Andrey; Beattie, Christopher
2016-01-01
Performance of the adjoint and adjoint-free 4-dimensional variational (4dVar) data assimilation techniques is compared in application to the hydrographic surveys and velocity observations collected in the Adriatic Sea in 2006. Assimilating the data into the Navy Coastal Ocean Model (NCOM) has shown that both methods deliver similar reduction of the cost function and demonstrate comparable forecast skill at approximately the same computational expense. The obtained optimal states were, however, significantly different in terms of distance from the background state: application of the adjoint method resulted in a 30-40% larger departure, mostly due to the excessive level of ageostrophic motions in the southern basin of the Sea that was not covered by observations.
Adjoint Algorithm for CAD-Based Shape Optimization Using a Cartesian Method
Nemec, Marian; Aftosmis, Michael J.
2004-01-01
Adjoint solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape optimization. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (geometric parameters that control the shape). More recently, emerging adjoint applications focus on the analysis problem, where the adjoint solution is used to drive mesh adaptation, as well as to provide estimates of functional error bounds and corrections. The attractive feature of this approach is that the mesh-adaptation procedure targets a specific functional, thereby localizing the mesh refinement and reducing computational cost. Our focus is on the development of adjoint-based optimization techniques for a Cartesian method with embedded boundaries.12 In contrast t o implementations on structured and unstructured grids, Cartesian methods decouple the surface discretization from the volume mesh. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin et developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the Euler equations. In both approaches, a boundary condition is introduced to approximate the effects of the evolving surface shape that results in accurate gradient computation. Central to automated shape optimization algorithms is the issue of geometry modeling and control. The need to optimize complex, "real-life" geometry provides a strong incentive for the use of parametric-CAD systems within the optimization procedure. In previous work, we presented
Weighted $L^p$ estimates for the area integral associated to self-adjoint operators
Gong, Ruming; Yan, Lixin
2011-01-01
This article is concerned with some weighted norm inequalities for the so-called horizontal (i.e. involving time derivatives) area integrals associated to a non-negative self-adjoint operator satisfying a pointwise Gaussian estimate for its heat kernel, as well as the corresponding vertical (i.e. involving space derivatives) area integrals associated to a non-negative self-adjoint operator satisfying in addition a pointwise upper bounds for the gradient of the heat kernel. As applications, we...
Self-adjoint Time Operator is the Rule for Discrete Semibounded Hamiltonians
Galapon, E A
2002-01-01
We prove explicitly that to every discrete, semibounded Hamiltonian with constant degeneracy and with finite sum of the squares of the reciprocal of its eigenvalues and whose eigenvectors span the entire Hilbert space there exists a characteristic self-adjoint time operator which is canonically conjugate to the Hamiltonian in a dense subspace of the Hilbert space. Moreover, we show that each characteristic time operator generates an uncountable class of self- adjoint operators canonically conjugate with the same Hamiltonian.
Jacobs, Christian T.; Piggott, Matthew D.; Kramer, Stephan C; Funke, Simon W.
2016-01-01
Extracting the optimal amount of power from an array of tidal turbines requires an intricate understanding of tidal dynamics and the effects of turbine placement on the local and regional scale flow. Numerical models have contributed significantly towards this understanding, and more recently, adjoint-based modelling has been employed to optimise the positioning of the turbines in an array in an automated way and improve on simple, regular man-made configurations. Adjoint-based optimisation o...
Efficient parameter estimation in 2D transport models based on an adjoint formalism
An adjoint based optimization procedure is elaborated to estimate transport coefficients for plasma edge models based on a limited set of known profiles at different locations. It is shown that a set of adjoint equations can accurately determine all sensitivities towards transport coefficients at once. A proof of principle is provided on a simple geometry. The methodology is subsequently applied to assess whether a simple edge model can be tuned toward full B2-EIRENE profiles for a JET-configuration. (paper)
Highlights: ► Successful noise control of a 2D-planejet with DNS resolution and a 3D-planejet with LES-resolution using adjoint method. ► Validation of gradient-information obtained with the continuous-adjoint approach, by comparing gradient with finite differences. ► Extension of control-interval with the receding horizon algorithm. - Abstract: A control optimization technique using the continuous adjoint of the compressible Navier–Stokes equations is implemented for aeroacoustic optimization of plane jet flows. The purpose of the adjoint equations is to provide sensitivity information, which is afterwards used in a gradient-based minimization of a prescribed cost functional, designed to describe the far-field sound pressure level (SPL). The objective of the present paper is to demonstrate the ability to reduce the sound in the near far-field of plane jets. Furthermore, as the continuous adjoint approach can become inaccurate, due to inconsistencies between the continuous and the discretized system, the accuracy of the continuous adjoint approach is investigated. The considered cases exhibit a nozzle exit Reynolds number of Rejet = ρujetD/μ = 2000 and a Mach number of Mjet = 0.9, performed using two-dimensional direct numerical simulation and three-dimensional large-eddy simulation, respectively. A comparison of the obtained gradient via adjoint and finite differences is presented and it is shown, that in order to obtain reliable gradient directions, the length of the optimization time needs to be restricted. Furthermore, a receding horizon optimization for the two-dimensional plane jet simulation is used to obtain a sound reduction over much longer time intervals. The influence of different formulations of the viscosity in the adjoint equations is finally investigated.
A self-adjoint arrival time operator inspired by measurement models
Highlights: • Construction of a self-adjoint arrival time operator inspired by measurements. • Agreement with the strong measurement formula in the low momentum regime. • Review of self-adjoint and non-self-adjoint arrival time operators. • Discussion of the momentum operator on the half-line. • Discussion of the intuitive reasons obstructing self-adjointness. - Abstract: We introduce an arrival time operator which is self-adjoint and, unlike previously proposed arrival time operators, has a close link to simple measurement models. Its spectrum leads to an arrival time distribution which is a variant of the Kijowski distribution (a re-ordering of the current) in the large momentum regime but is proportional to the kinetic energy density in the small momentum regime, in agreement with measurement models. A brief derivation of the latter distribution is given. We make some simple observations about the physical reasons for self-adjointness, or its absence, in both arrival time operators and the momentum operator on the half-line and we also compare our operator with the dwell time operator
Plumes, Hotspot & Slabs Imaged by Global Adjoint Tomography
Bozdag, E.; Lefebvre, M. P.; Lei, W.; Peter, D. B.; Smith, J. A.; Komatitsch, D.; Tromp, J.
2015-12-01
We present the "first generation" global adjoint tomography model based on 3D wave simulations, which is the result of 15 conjugate-gradient iterations with confined transverse isotropy to the upper mantle. Our starting model is the 3D mantle and crustal models S362ANI (Kustowski et al. 2008) and Crust2.0 (Bassin et al. 2000), respectively. We take into account the full nonlinearity of wave propagation in numerical simulations including attenuation (both in forward and adjoint simulations), topography/bathymetry, etc., using the GPU version of the SPECFEM3D_GLOBE package. We invert for crust and mantle together without crustal corrections to avoid any bias in mantle structure. We started with an initial selection of 253 global CMT events within the magnitude range 5.8 ≤ Mw ≤ 7.0 with numerical simulations having resolution down to 27 s combining 30-s body and 60-s surface waves. After the 12th iteration we increased the resolution to 17 s, including higher-frequency body waves as well as going down to 45 s in surface-wave measurements. We run 180-min seismograms and assimilate all minor- and major-arc body and surface waves. Our 15th iteration model update shows a tantalisingly enhanced image of the Tahiti plume as well as various other plumes and hotspots, such as Caroline, Galapagos, Yellowstone, Erebus, etc. Furthermore, we see clear improvements in slab resolution along the Hellenic and Japan Arcs, as well as subduction along the East of Scotia Plate, which does not exist in the initial model. Point-spread function tests (Fichtner & Trampert 2011) suggest that we are close to the resolution of continental-scale studies in our global inversions and able to confidently map features, for instance, at the scale of the Yellowstone hotspot. This is a clear consequence of our multi-scale smoothing strategy, in which we define our smoothing operator as a function of the approximate Hessian kernel and smooth our gradients less wherever we have good ray coverage
Conventional collapsing for group cross sections used in multigroup nuclear reactor calculations is usually performed using normal (real; direct) flux weighting. The application of more advanced collapsing procedures using in an appropriate manner real, adjoint and bilinear weighting was in the past restricted in general to fundamental mode problems. Although the principles have been published for more than ten years, there seems to exist little recent experience on the merits and possible difficulties of these improved procedures for multidimensional diffusion problems for practical purposes, e.g. in the nuclear design and analysis of large Liquid Metal Fast Breeder Reactors (LMFBRs). The present work indicates the nature of the problems which could possibly be encountered in applying these procedures by tracing them back to the known close correspondence between group collapsing and synthesis methods. It tries to explain certain somewhat unusual features of the collapsed group constants obtained by adjoint and bilinear weighting and describes the experience gained in representative 1-dim. and 2-dim. test cases. It could be shown for criticality and perturbation calculations that in general it is advantageous to apply these improved collapsing methods if the necessary precautions are taken. Compared to the conventional collapsing procedures these improved procedures are especially useful for multidimensional problems because their application is well suited for that purpose. In the present study it could be proven that they are favorable with respect to computer time and storage needed due to the fact that the necessary number of coarse groups can be kept fairly small without deteriorating too much the accuracy and reliability of the coarse group results compared to reference results of corresponding fine group calculations with uncollapsed group constants. (orig.)
Gopalakrishnan, Ganesh; Cornuelle, Bruce D.; Hoteit, Ibrahim; Rudnick, Daniel L.; Owens, W. Brechner
2013-07-01
An ocean state estimate has been developed for the Gulf of Mexico (GoM) using the MIT general circulation model and its adjoint. The estimate has been tested by forecasting loop current (LC) evolution and eddy shedding in the GoM. The adjoint (or four-dimensional variational) method was used to match the model evolution to observations by adjusting model temperature and salinity initial conditions, open boundary conditions, and atmospheric forcing fields. The model was fit to satellite-derived along-track sea surface height, separated into temporal mean and anomalies, and gridded sea surface temperature for 2 month periods. The optimized state at the end of the assimilation period was used to initialize the forecast for 2 months. Forecasts explore practical LC predictability and provide a cross-validation test of the state estimate by comparing it to independent future observations. The model forecast was tested for several LC eddy separation events, including Eddy Franklin in May 2010 during the deepwater horizon oil spill disaster in the GoM. The forecast used monthly climatological open boundary conditions, atmospheric forcing, and run-off fluxes. The model performance was evaluated by computing model-observation root-mean-square difference (rmsd) during both the hindcast and forecast periods. The rmsd metrics for the forecast generally outperformed persistence (keeping the initial state fixed) and reference (forecast initialized using assimilated Hybrid Coordinate Ocean Model 1/12° global analysis) model simulations during LC eddy separation events for a period of 1˜2 months.
Gopalakrishnan, Ganesh
2013-07-01
An ocean state estimate has been developed for the Gulf of Mexico (GoM) using the MIT general circulation model and its adjoint. The estimate has been tested by forecasting loop current (LC) evolution and eddy shedding in the GoM. The adjoint (or four-dimensional variational) method was used to match the model evolution to observations by adjusting model temperature and salinity initial conditions, open boundary conditions, and atmospheric forcing fields. The model was fit to satellite-derived along-track sea surface height, separated into temporal mean and anomalies, and gridded sea surface temperature for 2 month periods. The optimized state at the end of the assimilation period was used to initialize the forecast for 2 months. Forecasts explore practical LC predictability and provide a cross-validation test of the state estimate by comparing it to independent future observations. The model forecast was tested for several LC eddy separation events, including Eddy Franklin in May 2010 during the deepwater horizon oil spill disaster in the GoM. The forecast used monthly climatological open boundary conditions, atmospheric forcing, and run-off fluxes. The model performance was evaluated by computing model-observation root-mean-square difference (rmsd) during both the hindcast and forecast periods. The rmsd metrics for the forecast generally outperformed persistence (keeping the initial state fixed) and reference (forecast initialized using assimilated Hybrid Coordinate Ocean Model 1/12° global analysis) model simulations during LC eddy separation events for a period of 1̃2 months.
A practical discrete-adjoint method for high-fidelity compressible turbulence simulations
Vishnampet, Ramanathan [Department of Mechanical Science and Engineering, University of Illinois at Urbana–Champaign, Urbana, IL 61801 (United States); Bodony, Daniel J. [Department of Aerospace Engineering, University of Illinois at Urbana–Champaign, Urbana, IL 61801 (United States); Freund, Jonathan B., E-mail: jbfreund@illinois.edu [Department of Mechanical Science and Engineering, University of Illinois at Urbana–Champaign, Urbana, IL 61801 (United States); Department of Aerospace Engineering, University of Illinois at Urbana–Champaign, Urbana, IL 61801 (United States)
2015-03-15
Methods and computing hardware advances have enabled accurate predictions of complex compressible turbulence phenomena, such as the generation of jet noise that motivates the present effort. However, limited understanding of underlying physical mechanisms restricts the utility of such predictions since they do not, by themselves, indicate a route to design improvements. Gradient-based optimization using adjoints can circumvent the flow complexity to guide designs, though this is predicated on the availability of a sufficiently accurate solution of the forward and adjoint systems. These are challenging to obtain, since both the chaotic character of the turbulence and the typical use of discretizations near their resolution limits in order to efficiently represent its smaller scales will amplify any approximation errors made in the adjoint formulation. Formulating a practical exact adjoint that avoids such errors is especially challenging if it is to be compatible with state-of-the-art simulation methods used for the turbulent flow itself. Automatic differentiation (AD) can provide code to calculate a nominally exact adjoint, but existing general-purpose AD codes are inefficient to the point of being prohibitive for large-scale turbulence simulations. Here, we analyze the compressible flow equations as discretized using the same high-order workhorse methods used for many high-fidelity compressible turbulence simulations, and formulate a practical space–time discrete-adjoint method without changing the basic discretization. A key step is the definition of a particular discrete analog of the continuous norm that defines our cost functional; our selection leads directly to an efficient Runge–Kutta-like scheme, though it would be just first-order accurate if used outside the adjoint formulation for time integration, with finite-difference spatial operators for the adjoint system. Its computational cost only modestly exceeds that of the flow equations. We confirm that
Use of adjoint transport solutions for inverse problems
Inverse problems represent a third major class of transport problems (after criticality-type problems and shielding-type problems). For the purpose of discussion here, an inverse problem is defined to be one in which a medium of interest emits some particular type of radiation, either spontaneously or when exposed to some other type of radiation. From a knowledge of the emitted radiation characteristics, one wishes to infer something about the medium that emits the detected radiation. An important example of a practical inverse problem is the geophysical logging of geological formations to determine the presence of minerals or hydrocarbons in the vicinity of a borehole. Proposals have also been made to detect hidden explosives by these methods. For a complete understanding of the apparatus used to solve inverse problems, it is important to know the characteristic interrogation distance of the method. Phrased more precisely, What is the relative contribution to the detector signal due to radiation particles emitted by different elements of volume in the medium being interrogated? As shown in the paper, it is possible to address this question with a computational procedure that utilizes adjoint transport solutions
Adjoint sensitivity analysis of hydrodynamic stability in cyclonic flows
Guzman Inigo, Juan; Juniper, Matthew
2015-11-01
Cyclonic separators are used in a variety of industries to efficiently separate mixtures of fluid and solid phases by means of centrifugal forces and gravity. In certain circumstances, the vortex core of cyclonic flows is known to precess due to the instability of the flow, which leads to performance reductions. We aim to characterize the unsteadiness using linear stability analysis of the Reynolds Averaged Navier-Stokes (RANS) equations in a global framework. The system of equations, including the turbulence model, is linearised to obtain an eigenvalue problem. Unstable modes corresponding to the dynamics of the large structures of the turbulent flow are extracted. The analysis shows that the most unstable mode is a helical motion which develops around the axis of the flow. This result is in good agreement with LES and experimental analysis, suggesting the validity of the approach. Finally, an adjoint-based sensitivity analysis is performed to determine the regions of the flow that, when altered, have most influence on the frequency and growth-rate of the unstable eigenvalues.
Adjoint sources, disconnected loops and other fruit of lattice QCD
Foster, M S
1998-01-01
eta' meson mass in full QCD. We introduce related source pairs to minimise the variance of the disconnected loop operators that are employed. We are able to obtain estimates of the mass from a very modest number of gauge configurations and purely local operators. We did not observe any evidence of unquenching from the measurements obtained, though statistical noise dominated the signal in the region where definitive effects would be seen. We undertake a comprehensive study of the gluelump mass spectrum, exploring the spin structure of the state. We use five lattice spacings from which we extract continuum values for the state splittings, with high statistics employed at beta = 6.0. We conduct a low statistics study of a related and previously unexamined lattice state, which we term the adjoint meson. We find that this state is fractionally more massive than the gluelump on the lattices considered, with indications that the splitting is greater than the pion mass at beta = 6.0. We investigate a sum rule approa...
Intensity-modulated radiation therapy (IMRT) is a new technique for administering external beam radiation therapy. This technology modulates the intensity and shape of the treatment beam as a function of source position and patient anatomy. This process of conforming the source to the patient requires the optimization of the independent variables of the source field. In this study, adjoint Monte Carlo methods were used to compute the sensitivity field that corresponds to a prescribed dose distribution. Given these data, linear and nonlinear optimization models were constructed with a simplified geometry to compute an optimized set of beams to deliver a desired dose distribution. The dose delivered to voxel i by beam j (Dij) influence matrix may be obtained from solutions to the adjoint transport equation. These solutions provide the sensitivity of the prescribed dose at a single point in the patient to all possible points in the source field. For this investigation, the source field consisted of 36 possible positions along a circular gantry. Each position had 21 possible directions to aim at the patient. Beam weights could vary continuously, and beam energy spectra matched that of a hospital-based linear accelerator. The MCNP Monte Carlo code was used to transport adjoint particles from each patient voxel to the 36 possible source locations where they were binned by direction and energy. The patient voxels (1 cm3) were defined within the central slice of a block phantom (31 x 31 x 11 cm) of unit-density water. The adjoint source for each voxel was the flux-to-dose conversion factor for tissue. The bin structures for the tallies matched the direction and energy structure of the forward source. Figure 1 shows the dose volume histograms (DVHs) for the optimized dose distributions for a ring-shaped tumor surrounding a sensitive structure. The DVH reports the fraction of each tissue type that is raised to each dose level. The lower dose limit prescribed for TU was 2
Adjoint-based airfoil shape optimization in transonic flow
Gramanzini, Joe-Ray
The primary focus of this work is efficient aerodynamic shape optimization in transonic flow. Adjoint-based optimization techniques are employed on airfoil sections and evaluated in terms of computational accuracy as well as efficiency. This study examines two test cases proposed by the AIAA Aerodynamic Design Optimization Discussion Group. The first is a two-dimensional, transonic, inviscid, non-lifting optimization of a Modified-NACA 0012 airfoil. The second is a two-dimensional, transonic, viscous optimization problem using a RAE 2822 airfoil. The FUN3D CFD code of NASA Langley Research Center is used as the ow solver for the gradient-based optimization cases. Two shape parameterization techniques are employed to study their effect and the number of design variables on the final optimized shape: Multidisciplinary Aerodynamic-Structural Shape Optimization Using Deformation (MASSOUD) and the BandAids free-form deformation technique. For the two airfoil cases, angle of attack is treated as a global design variable. The thickness and camber distributions are the local design variables for MASSOUD, and selected airfoil surface grid points are the local design variables for BandAids. Using the MASSOUD technique, a drag reduction of 72.14% is achieved for the NACA 0012 case, reducing the total number of drag counts from 473.91 to 130.59. Employing the BandAids technique yields a 78.67% drag reduction, from 473.91 to 99.98. The RAE 2822 case exhibited a drag reduction from 217.79 to 132.79 counts, a 39.05% decrease using BandAids.
Assessing the Impact of Observations on Numerical Weather Forecasts Using the Adjoint Method
Gelaro, Ronald
2012-01-01
The adjoint of a data assimilation system provides a flexible and efficient tool for estimating observation impacts on short-range weather forecasts. The impacts of any or all observations can be estimated simultaneously based on a single execution of the adjoint system. The results can be easily aggregated according to data type, location, channel, etc., making this technique especially attractive for examining the impacts of new hyper-spectral satellite instruments and for conducting regular, even near-real time, monitoring of the entire observing system. This talk provides a general overview of the adjoint method, including the theoretical basis and practical implementation of the technique. Results are presented from the adjoint-based observation impact monitoring tool in NASA's GEOS-5 global atmospheric data assimilation and forecast system. When performed in conjunction with standard observing system experiments (OSEs), the adjoint results reveal both redundancies and dependencies between observing system impacts as observations are added or removed from the assimilation system. Understanding these dependencies may be important for optimizing the use of the current observational network and defining requirements for future observing systems
STUDY ON THE ADJOINT METHOD IN DATA ASSIMILATION AND THE RELATED PROBLEMS
吕咸青; 吴自库; 谷艺; 田纪伟
2004-01-01
It is not reasonable that one can only use the adjoint of model in data assimilation.The simulated numerical experiment shows that for the tidal model,the result of the adjoint of equation is almost the same as that of the adjoint of model:the averaged absolute difference of the amplitude between observations and simulation is less than 5.0 cm and that of the phase-lag is less than 5.0°.The results are both in good agreement with the observed M2 tide in the Bohai Sea and the Yellow Sea.For comparison,the traditional methods also have been used to simulate M2 tide in the Bohai Sea and the Yellow Sea.The initial guess values of the boundary conditions are given first,and then are adjusted to acquire the simulated results that are as close as possible to the observations.As the boundary conditions contain 72 values,which should be adjusted and how to adjust them can only be partially solved by adjusting them many times.The satisfied results are hard to acquire even gigantic efforts are done.Here,the automation of the treatment of the open boundary conditions is realized.The method is unique and superior to the traditional methods.It is emphasized that if the adjoint of equation is used,tedious and complicated mathematical deduction can be avoided.Therefore the adjoint of equation should attract much attention.
Light Adjoint Quarks in the Instanton-Dyon Liquid Model IV
Liu, Yizhuang; Zahed, Ismail
2016-01-01
We discuss the instanton-dyon liquid model with $N_f$ Majorana quark flavors in the adjoint representation of color $SU_c(2)$ at finite temperature. We briefly recall the index theorem on $S^1\\times R^3$ for twisted adjoint fermions in a BPS dyon background of arbitrary holonomy, and use the ADHM construction to explicit the adjoint anti-periodic zero modes. We use these results to derive the partition function of an interacting instanton-dyon ensemble with $N_f$ light and anti-periodic adjoint quarks. We develop the model in details by mapping the theory on a 3-dimensional quantum effective theory with adjoint quarks with manifest $SU(N_f)\\times Z_{4N_f}$ symmetry. Using a mean-field analysis at weak coupling and strong screening, we show that center symmetry requires the spontaneous breaking of chiral symmetry, which is shown to only take place for $N_f=1$. For a sufficiently dense liquid, we find that the ground state is center symmetric and breaks spontaneously flavor symmetry through $SU(N_f)\\times Z_{4N...
Highlights: • We present a Monte Carlo method for computing the adjoint-weighted kinetics parameters via the IFP algorithm. • Extensive verification tests are performed on simple models. • Several validation tests are performed on the measured values of effective delayed neutron fraction and Rossi alpha. - Abstract: The analysis of neutron kinetics relies on the knowledge of adjoint-weighted kinetics parameters, which are key to safety issues in the context of transient or accidental reactor behavior. The Iterated Fission Probability (IFP) method allows the adjoint-weighted mean generation time and delayed neutron fraction to be computed within a Monte Carlo power iteration calculation. In this work we describe the specific features of the implementation of the IFP algorithm in the reference Monte Carlo code TRIPOLI-4® developed at CEA. Several verification and validation tests are discussed, and the impact of nuclear data libraries, IFP cycle length and inter-cycle correlations are analyzed in detail
Automated derivation of the adjoint of high-level transient finite element programs
Farrell, Patrick E; Funke, Simon F; Rognes, Marie E
2012-01-01
In this paper we demonstrate the capability of automatically deriving the discrete adjoint and tangent linear models from a forward model written in the high-level FEniCS finite element computing environment. In contrast to developing a model directly in Fortran or C++, high-level systems allow the developer to express the variational problems to be solved in near-mathematical notation. As such, these systems have a key advantage: since the mathematical structure of the problem is preserved, they are more amenable to automated analysis and manipulation. Our approach to automated adjoint derivation relies on run-time annotation of the temporal structure of the model, and employs the same finite element form compiler to automatically generate the low-level code for the derived models. The approach requires only trivial changes to a large class of forward models, including complicated time-dependent nonlinear models. The adjoint model automatically employs optimal checkpointing schemes to mitigate storage requir...
Adjoint Based A Posteriori Analysis of Multiscale Mortar Discretizations with Multinumerics
Tavener, Simon
2013-01-01
In this paper we derive a posteriori error estimates for linear functionals of the solution to an elliptic problem discretized using a multiscale nonoverlapping domain decomposition method. The error estimates are based on the solution of an appropriately defined adjoint problem. We present a general framework that allows us to consider both primal and mixed formulations of the forward and adjoint problems within each subdomain. The primal subdomains are discretized using either an interior penalty discontinuous Galerkin method or a continuous Galerkin method with weakly imposed Dirichlet conditions. The mixed subdomains are discretized using Raviart- Thomas mixed finite elements. The a posteriori error estimate also accounts for the errors due to adjoint-inconsistent subdomain discretizations. The coupling between the subdomain discretizations is achieved via a mortar space. We show that the numerical discretization error can be broken down into subdomain and mortar components which may be used to drive adaptive refinement.Copyright © by SIAM.
An adjoint-based approach for finding invariant solutions of Navier-Stokes equations
Farazmand, Mohammad
2015-01-01
We consider the incompressible Navier--Stokes equations with periodic boundary conditions and time-independent forcing. For this type of flow, we derive adjoint equations whose trajectories converge asymptotically to the equilibrium and traveling wave solutions of the Navier--Stokes equations. Using the adjoint equations, arbitrary initial conditions evolve to the vicinity of a (relative) equilibrium at which point a few Newton-type iterations yield the desired (relative) equilibrium solution. We apply this adjoint-based method to a chaotic two-dimensional Kolmogorov flow. A convergence rate of 100% is observed, leading to the discovery of 21 new steady state and traveling wave solutions at Reynolds number Re=40. Some of the new invariant solutions have spatially localized structures that were previously believed to only exist on domains with large aspect ratios. We show that one of the newly found steady state solutions underpins the temporal intermittencies, i.e., high energy dissipation episodes of the flo...
Tracking influential haze source areas in North China using an adjoint model, GRAPES-CUACE
An, X. Q.; Zhai, S. X.; Jin, M.; Gong, S. L.; Wang, Y.
2015-08-01
Based upon the adjoint theory, the adjoint of the aerosol module in the atmospheric chemical modeling system GRAPES-CUACE (Global/Regional Assimilation and PrEdiction System coupled with the CMA Unified Atmospheric Chemistry Environment) was developed and tested for its correctness. Through statistic comparison, BC (black carbon aerosol) concentrations simulated by GRAPES-CUACE were generally consistent with observations from Nanjiao (one urban observation station) and Shangdianzi (one rural observation station) stations. To track the most influential emission-sources regions and the most influential time intervals for the high BC concentration during the simulation period, the adjoint model was adopted to simulate the sensitivity of average BC concentration over Beijing at the highest concentration time point (referred to as the Objective Function) with respect to BC emission amount over Beijing-Tianjin-Hebei region. Four types of regions were selected based on administrative division and sensitivity coefficient distribution. The adjoint model was used to quantify the effects of emission-sources reduction in different time intervals over different regions by one independent simulation. Effects of different emission reduction strategies based on adjoint sensitivity information show that the more influential regions (regions with relatively larger sensitivity coefficients) do not necessarily correspond to the administrative regions, and the influence effectiveness of sensitivity-oriented regions was greater than the administrative divisions. The influence of emissions on the objective function decreases sharply approximately for the pollutants emitted 17-18 h ago in this episode. Therefore, controlling critical emission regions during critical time intervals on the basis of adjoint sensitivity analysis is much more efficient than controlling administrative specified regions during an experiential time period.
Yamaleev, N. K.; Diskin, B.; Nielsen, E. J.
2009-01-01
.We study local-in-time adjoint-based methods for minimization of ow matching functionals subject to the 2-D unsteady compressible Euler equations. The key idea of the local-in-time method is to construct a very accurate approximation of the global-in-time adjoint equations and the corresponding sensitivity derivative by using only local information available on each time subinterval. In contrast to conventional time-dependent adjoint-based optimization methods which require backward-in-time integration of the adjoint equations over the entire time interval, the local-in-time method solves local adjoint equations sequentially over each time subinterval. Since each subinterval contains relatively few time steps, the storage cost of the local-in-time method is much lower than that of the global adjoint formulation, thus making the time-dependent optimization feasible for practical applications. The paper presents a detailed comparison of the local- and global-in-time adjoint-based methods for minimization of a tracking functional governed by the Euler equations describing the ow around a circular bump. Our numerical results show that the local-in-time method converges to the same optimal solution obtained with the global counterpart, while drastically reducing the memory cost as compared to the global-in-time adjoint formulation.
The Hamiltonian of a system of quantum particles minimally coupled to a quantum field is considered for arbitrary coupling constants. The Hamiltonian has a translation invariant part. By means of functional integral representations the existence of an invariant domain under the action of the heat semigroup generated by a self-adjoint extension of the translation invariant part is shown. With a non-perturbative approach it is proved that the Hamiltonian is essentially self-adjoint on a domain. A typical example is the Pauli-Fierz model with spin 1/2 in nonrelativistic quantum electrodynamics for arbitrary coupling constants. (orig.)
Bounds on variation of spectral subspaces under J-self-adjoint perturbations
Albeverio, S.; Motovilov, A. K.; Shkalikov, A. A.
2008-01-01
Let $A$ be a self-adjoint operator on a Hilbert space $\\fH$. Assume that the spectrum of $A$ consists of two disjoint components $\\sigma_0$ and $\\sigma_1$. Let $V$ be a bounded operator on $\\fH$, off-diagonal and $J$-self-adjoint with respect to the orthogonal decomposition $\\fH=\\fH_0\\oplus\\fH_1$ where $\\fH_0$ and $\\fH_1$ are the spectral subspaces of $A$ associated with the spectral sets $\\sigma_0$ and $\\sigma_1$, respectively. We find (optimal) conditions on $V$ guaranteeing that the pertur...
Adjoint sensitivity studies of loop current and eddy shedding in the Gulf of Mexico
Gopalakrishnan, Ganesh
2013-07-01
Adjoint model sensitivity analyses were applied for the loop current (LC) and its eddy shedding in the Gulf of Mexico (GoM) using the MIT general circulation model (MITgcm). The circulation in the GoM is mainly driven by the energetic LC and subsequent LC eddy separation. In order to understand which ocean regions and features control the evolution of the LC, including anticyclonic warm-core eddy shedding in the GoM, forward and adjoint sensitivities with respect to previous model state and atmospheric forcing were computed using the MITgcm and its adjoint. Since the validity of the adjoint model sensitivities depends on the capability of the forward model to simulate the real LC system and the eddy shedding processes, a 5 year (2004–2008) forward model simulation was performed for the GoM using realistic atmospheric forcing, initial, and boundary conditions. This forward model simulation was compared to satellite measurements of sea-surface height (SSH) and sea-surface temperature (SST), and observed transport variability. Despite realistic mean state, standard deviations, and LC eddy shedding period, the simulated LC extension shows less variability and more regularity than the observations. However, the model is suitable for studying the LC system and can be utilized for examining the ocean influences leading to a simple, and hopefully generic LC eddy separation in the GoM. The adjoint sensitivities of the LC show influences from the Yucatan Channel (YC) flow and Loop Current Frontal Eddy (LCFE) on both LC extension and eddy separation, as suggested by earlier work. Some of the processes that control LC extension after eddy separation differ from those controlling eddy shedding, but include YC through-flow. The sensitivity remains stable for more than 30 days and moves generally upstream, entering the Caribbean Sea. The sensitivities of the LC for SST generally remain closer to the surface and move at speeds consistent with advection by the high-speed core of
Adaptive mesh refinement and adjoint methods in geophysics simulations
Burstedde, Carsten
2013-04-01
required by human intervention and analysis. Specifying an objective functional that quantifies the misfit between the simulation outcome and known constraints and then minimizing it through numerical optimization can serve as an automated technique for parameter identification. As suggested by the similarity in formulation, the numerical algorithm is closely related to the one used for goal-oriented error estimation. One common point is that the so-called adjoint equation needs to be solved numerically. We will outline the derivation and implementation of these methods and discuss some of their pros and cons, supported by numerical results.
The features and the algorithm of the program to calculate adjoint neutron cross sections on the basis of the continuous energy neutron cross sections as well as energy and angular distributions are described. The calculated adjoint cross sections are intended for Monte Carlo investigation of the nonuniform adjoint Boltzmann equation. 16 refs
Highlights: ► We develop an abstract framework for computing the adjoint to the neutron/nuclide burnup equations posed as a system of differential algebraic equations. ► We validate use of the adjoint for computing both sensitivity to uncertain inputs and for estimating global time discretization error. ► Flexibility of the framework is leveraged to add heat transfer physics and compute its adjoint without a reformulation of the adjoint system. ► Such flexibility is crucial for high performance computing applications. -- Abstract: We develop a general framework for computing the adjoint variable to nuclear engineering problems governed by a set of differential–algebraic equations (DAEs). The nuclear engineering community has a rich history of developing and applying adjoints for sensitivity calculations; many such formulations, however, are specific to a certain set of equations, variables, or solution techniques. Any change or addition to the physics model would require a reformulation of the adjoint problem and substantial difficulties in its software implementation. In this work we propose an abstract framework that allows for the modification and expansion of the governing equations, leverages the existing theory of adjoint formulation for DAEs, and results in adjoint equations that can be used to efficiently compute sensitivities for parametric uncertainty quantification. Moreover, as we justify theoretically and demonstrate numerically, the same framework can be used to estimate global time discretization error. We first motivate the framework and show that the coupled Bateman and transport equations, which govern the time-dependent neutronic behavior of a nuclear reactor, may be formulated as a DAE system with a power constraint. We then use a variational approach to develop the parameter-dependent adjoint framework and apply existing theory to give formulations for sensitivity and global time discretization error estimates using the adjoint
Möbius invariant BFKL equation for the adjoint representation in N=4 SUSY
Fadin, V.S., E-mail: fadin@inp.nsk.su [Budker Institute of Nuclear Physics of SD RAS, 630090 Novosibirsk (Russian Federation); Novosibirsk State University, 630090 Novosibirsk (Russian Federation); Fiore, R., E-mail: roberto.fiore@cs.infn.it [Dipartimento di Fisica, Università della Calabria, and Istituto Nazionale di Fisica Nucleare, Gruppo collegato di Cosenza, Arcavacata di Rende, I-87036 Cosenza (Italy); Lipatov, L.N., E-mail: lipatov@thd.pnpi.spb.ru [Petersburg Nuclear Physics Institute and St. Petersburg State University, Gatchina, 188300 St. Petersburg (Russian Federation); Papa, A., E-mail: alessandro.papa@cs.infn.it [Dipartimento di Fisica, Università della Calabria, and Istituto Nazionale di Fisica Nucleare, Gruppo collegato di Cosenza, Arcavacata di Rende, I-87036 Cosenza (Italy)
2013-09-01
It is shown that in the next-to-leading approximation of N=4 SUSY the BFKL equation for two-gluon composite states in the adjoint representation of the gauge group can be reduced to a form which is invariant under Möbius transformation in the momentum space. The corresponding similarity transformation of its integral kernel is constructed in an explicit way.
On a class of non-self-adjoint periodic boundary value problems with discrete real spectrum
Boulton, Lyonell; Levitin, Michael; Marletta, Marco
2010-01-01
In [arXiv:0801.0172] we examined a family of periodic Sturm-Liouville problems with boundary and interior singularities which are highly non-self-adjoint but have only real eigenvalues. We now establish Schatten class properties of the associated resolvent operator.
Coupling of Monte Carlo adjoint leakages with three-dimensional discrete ordinates forward fluences
A computer code, DRC3, has been developed for coupling Monte Carlo adjoint leakages with three-dimensional discrete ordinates forward fluences in order to solve a special category of geometrically-complex deep penetration shielding problems. The code extends the capabilities of earlier methods that coupled Monte Carlo adjoint leakages with two-dimensional discrete ordinates forward fluences. The problems involve the calculation of fluences and responses in a perturbation to an otherwise simple two- or three-dimensional radiation field. In general, the perturbation complicates the geometry such that it cannot be modeled exactly using any of the discrete ordinates geometry options and thus a direct discrete ordinates solution is not possible. Also, the calculation of radiation transport from the source to the perturbation involves deep penetration. One approach to solving such problems is to perform the calculations in three steps: (1) a forward discrete ordinates calculation, (2) a localized adjoint Monte Carlo calculation, and (3) a coupling of forward fluences from the first calculation with adjoint leakages from the second calculation to obtain the response of interest (fluence, dose, etc.). A description of this approach is presented along with results from test problems used to verify the method. The test problems that were selected could also be solved directly by the discrete ordinates method. The good agreement between the DRC3 results and the direct-solution results verify the correctness of DRC3
Coupling of Monte Carlo adjoint leakages with three-dimensional discrete ordinates forward fluences
Slater, C.O.; Lillie, R.A.; Johnson, J.O.; Simpson, D.B.
1998-04-01
A computer code, DRC3, has been developed for coupling Monte Carlo adjoint leakages with three-dimensional discrete ordinates forward fluences in order to solve a special category of geometrically-complex deep penetration shielding problems. The code extends the capabilities of earlier methods that coupled Monte Carlo adjoint leakages with two-dimensional discrete ordinates forward fluences. The problems involve the calculation of fluences and responses in a perturbation to an otherwise simple two- or three-dimensional radiation field. In general, the perturbation complicates the geometry such that it cannot be modeled exactly using any of the discrete ordinates geometry options and thus a direct discrete ordinates solution is not possible. Also, the calculation of radiation transport from the source to the perturbation involves deep penetration. One approach to solving such problems is to perform the calculations in three steps: (1) a forward discrete ordinates calculation, (2) a localized adjoint Monte Carlo calculation, and (3) a coupling of forward fluences from the first calculation with adjoint leakages from the second calculation to obtain the response of interest (fluence, dose, etc.). A description of this approach is presented along with results from test problems used to verify the method. The test problems that were selected could also be solved directly by the discrete ordinates method. The good agreement between the DRC3 results and the direct-solution results verify the correctness of DRC3.
Feynman's Operational Calculi: Spectral Theory for Noncommuting Self-adjoint Operators
Jefferies, Brian [University of New South Wales, School of Mathematics (Australia)], E-mail: b.jefferies@unsw.edu.au; Johnson, Gerald W. [333 Avery Hall, University of Nebraska, Lincoln, Department of Mathematics (United States)], E-mail: gjohnson@math.unl.edu; Nielsen, Lance [Creighton University, Department of Mathematics (United States)], E-mail: lnielsen@creighton.edu
2007-02-15
The spectral theorem for commuting self-adjoint operators along with the associated functional (or operational) calculus is among the most useful and beautiful results of analysis. It is well known that forming a functional calculus for noncommuting self-adjoint operators is far more problematic. The central result of this paper establishes a rich functional calculus for any finite number of noncommuting (i.e. not necessarily commuting) bounded, self-adjoint operators A{sub 1},..., A{sub n} and associated continuous Borel probability measures {mu}{sub 1}, ?, {mu}{sub n} on [0,1]. Fix A{sub 1},..., A{sub n}. Then each choice of an n-tuple ({mu}{sub 1},...,{mu}{sub n}) of measures determines one of Feynman's operational calculi acting on a certain Banach algebra of analytic functions even when A{sub 1}, ..., A{sub n} are just bounded linear operators on a Banach space. The Hilbert space setting along with self-adjointness allows us to extend the operational calculi well beyond the analytic functions. Using results and ideas drawn largely from the proof of our main theorem, we also establish a family of Trotter product type formulas suitable for Feynman's operational calculi.
Popov, A. V., E-mail: Popov.Barnaul@mail.ru [Polzunov Altai State Technical University (Russian Federation)
2016-01-15
Metallic lithium is used to demonstrate the possibilities of applying non-self-adjoint operators for quantitative description of orbital excitations of electrons in crystals. It is shown that, the nonequilibrium distribution function can be calculated when solving the spectral problem; therefore, the kinetic properties of a material can also be described with the unified band theory.
On the Finite Volume Element Method for Self-Adjoint Parabolic Integrodifferential Equations
Mohamed Bahaj; Anas Rachid
2013-01-01
Finite volume element schemes for non-self-adjoint parabolic integrodifferential equations are derived and stated. For the spatially discrete scheme, optimal-order error estimates in , , and , norms for are obtained. In this paper, we also study the lumped mass modification. Based on the Crank-Nicolson method, a time discretization scheme is discussed and related error estimates are derived.
Adjoint sensitivity in PDE constrained least squares problems as a multiphysics problem
Lahaye, D.; Mulckhuyse, W.F.W.
2012-01-01
Purpose - The purpose of this paper is to provide a framework for the implementation of an adjoint sensitivity formulation for least-squares partial differential equations constrained optimization problems exploiting a multiphysics finite elements package. The estimation of the diffusion coefficient
Brane Configurations for Nonsupersymmetric Meta-Stable Vacua in SQCD with Adjoint Matter
Ahn, C
2006-01-01
We present the configurations of intersecting branes in type IIA string theory corresponding to the meta-stable supersymmetry breaking vacua(hep-th/0608063) in the four-dimensional N}=1 supersymmetric Yang-Mills theory coupled massive flavors with adjoint matter where the superpotential has three deformed terms.
Using adjoint-based optimization to study wing flexibility in flapping flight
Wei, Mingjun; Xu, Min; Dong, Haibo
2014-11-01
In the study of flapping-wing flight of birds and insects, it is important to understand the impact of wing flexibility/deformation on aerodynamic performance. However, the large control space from the complexity of wing deformation and kinematics makes usual parametric study very difficult or sometimes impossible. Since the adjoint-based approach for sensitivity study and optimization strategy is a process with its cost independent of the number of input parameters, it becomes an attractive approach in our study. Traditionally, adjoint equation and sensitivity are derived in a fluid domain with fixed solid boundaries. Moving boundary is only allowed when its motion is not part of control effort. Otherwise, the derivation becomes either problematic or too complex to be feasible. Using non-cylindrical calculus to deal with boundary deformation solves this problem in a very simple and still mathematically rigorous manner. Thus, it allows to apply adjoint-based optimization in the study of flapping wing flexibility. We applied the ``improved'' adjoint-based method to study the flexibility of both two-dimensional and three-dimensional flapping wings, where the flapping trajectory and deformation are described by either model functions or real data from the flight of dragonflies. Supported by AFOSR.
Theory of Loops and Strings with Matter in the Adjoint Representation
Maharana, Jnanadeva; Singh, Lambodhar P.
1994-01-01
We have presented canonical and path integral formulations of a theory of loops and closed strings with the matter field quanta transforming in the adjoint representation of the SU(N) gauge group. The physical processes arising out of the interactions of loops and closed strings are discussed.
Wells, K. C.; Millet, D. B.; Bousserez, N.; Henze, D. K.; Chaliyakunnel, S.; Griffis, T. J.; Dlugokencky, E. J.; Prinn, R. G.; O'Doherty, S.; Weiss, R. F.; Dutton, G. S.; Elkins, J. W.; Krummel, P. B.; Langenfelds, R. L.; Steele, P.
2015-12-01
Nitrous oxide (N2O) is a long-lived greenhouse gas with a global warming potential approximately 300 times that of CO2, and plays a key role in stratospheric ozone depletion. Human perturbation of the nitrogen cycle has led to a rise in atmospheric N2O, but large uncertainties exist in the spatial and temporal distribution of its emissions. Here we employ a 4D-Var inversion framework for N2O based on the GEOS-Chem chemical transport model and its adjoint to derive new constraints on the space-time distribution of global land and ocean N2O fluxes. Based on an ensemble of global surface measurements, we find that emissions are overestimated over Northern Hemisphere land areas and underestimated in the Southern Hemisphere. Assigning these biases to particular land or ocean regions is more difficult given the long lifetime of N2O. To quantitatively evaluate where the current N2O observing network provides local and regional emission constraints, we apply a new, efficient information content analysis technique involving radial basis functions. The technique yields optimal state vector dimensions for N2O source inversions, with model grid cells grouped in space and time according to the resolution that can actually be provided by the network of global observations. We then use these optimal state vectors in an analytical inversion to refine current top-down emission estimates.
Healy, R.W.; Russell, T.F.
1992-01-01
A finite-volume Eulerian-Lagrangian local adjoint method for solution of the advection-dispersion equation is developed and discussed. The method is mass conservative and can solve advection-dominated ground-water solute-transport problems accurately and efficiently. An integrated finite-difference approach is used in the method. A key component of the method is that the integral representing the mass-storage term is evaluated numerically at the current time level. Integration points, and the mass associated with these points, are then forward tracked up to the next time level. The number of integration points required to reach a specified level of accuracy is problem dependent and increases as the sharpness of the simulated solute front increases. Integration points are generally equally spaced within each grid cell. For problems involving variable coefficients it has been found to be advantageous to include additional integration points at strategic locations in each well. These locations are determined by backtracking. Forward tracking of boundary fluxes by the method alleviates problems that are encountered in the backtracking approaches of most characteristic methods. A test problem is used to illustrate that the new method offers substantial advantages over other numerical methods for a wide range of problems.
The adjoint sensitivity method of global electromagnetic induction for CHAMP magnetic data
Complete text of publication follows. Martinec and McCreadie (2004) developed a time-domain spectral-finite element approach for the forward modelling of electromagnetic induction vector data as measured by the CHAMP satellite. Here, we present a new method of computing the sensitivity of the CHAMP electromagnetic induction data on the Earth's mantle electrical conductivity, which we term the adjoint sensitivity method. The forward and adjoint initial boundary-value problems, both solved in the time domain, are identical, except for the specification of prescribed boundary conditions. The respective boundary-value data at the satellite's altitude are the X magnetic component measured by the CHAMP vector magnetometer along satellite tracks for the forward method and the difference between the measured and predicted Z magnetic component for the adjoint method. The squares of these differences summed up over all CHAMP tracks determine the misfit. The sensitivity of the CHAMP data, that is the partial derivatives of the misfit function with respect to mantle conductivity parameters, are then determined by the scalar product of the forward and adjoint solutions, multiplied by the gradient of the conductivity and integrated over all CHAMP tracks. Such exactly determined sensitivities are checked against numerical differentiation of the misfit, and very good agreement is obtained. The attractiveness of the adjoint method lies in the fact that the adjoint sensitivities are calculated for little cost, regardless of the number of conductivity parameters. However, since the adjoint solution proceeds backwards in time, the forward solution must be stored at each time step, leading to memory requirements that are linear with respect to the number of steps undertaken. Having determined the sensitivities, we apply the conjugate gradient method to infer 1-D and 2-D conductivity structures of the Earth based on the CHAMP residual time serie (after the subtraction of static field
WANG Caixia; Paola Malanotte-Rizzoli
2014-01-01
The linkage between physical and biological processes, particularly the effect of the circulation field on the distribution of phytoplankton, is studied by applying a two-dimensional model and an adjoint data assimilation approach to the Gulf of Maine-Georges Bank region. The model results, comparing well with observation data, reveal seasonal and geographic variations of phytoplankton concentration and verify that the seasonal cycles of phytoplankton are controlled by both biological sources and ad-vection processes which are functions of space and time and counterbalance each other. Although advective flux divergences have greater magnitudes on Georges Bank than in the coastal region of the western Gulf of Maine, advection control over phytoplankton concentration is more significant in the coastal region of the western Gulf of Maine. The model results also suggest that the two separated populations in the coastal regions of the western Gulf of Maine and on Georges Bank are self-sustaining.
Marotzke, Jochem; Giering, Ralf; Zhang, Kate Q.; Stammer, Detlef; Hill, Chris; Lee, Tong
1999-12-01
We first describe the principles and practical considerations behind the computer generation of the adjoint to the Massachusetts Institute of Technology ocean general circulation model (GCM) using R. Giering's software tool Tangent-Linear and Adjoint Model Compiler (TAMC). The TAMC's recipe for (FORTRAN-) line-by-line generation of adjoint code is explained by interpreting an adjoint model strictly as the operator that gives the sensitivity of the output of a model to its input. Then, the sensitivity of 1993 annual mean heat transport across 29°N in the Atlantic, to the hydrography on January 1, 1993, is calculated from a global solution of the GCM. The "kinematic sensitivity" to initial temperature variations is isolated, showing how the latter would influence heat transport if they did not affect the density and hence the flow. Over 1 year the heat transport at 29°N is influenced kinematically from regions up to 20° upstream in the western boundary current and up to 5° upstream in the interior. In contrast, the dynamical influences of initial temperature (and salinity) perturbations spread from as far as the rim of the Labrador Sea to the 29°N section along the western boundary. The sensitivities calculated with the adjoint compare excellently to those from a perturbation calculation with the dynamical model. Perturbations in initial interior salinity influence meridional overturning and heat transport when they have propagated to the western boundary and can thus influence the integrated east-west density difference. Our results support the notion that boundary monitoring of meridional mass and heat transports is feasible.
Inversion of CO and NOx emissions using the adjoint of the IMAGES model
J.-F. Müller
2005-01-01
Full Text Available We use ground-based observations of CO mixing ratios and vertical column abundances together with tropospheric NO2 columns from the GOME satellite instrument as constraints for improving the global annual emission estimates of CO and NOx for the year 1997. The agreement between concentrations calculated by the global 3-dimensional CTM IMAGES and the observations is optimized using the adjoint modelling technique, which allows to invert for CO and NOx fluxes simultaneously, taking their chemical interactions into account. Our analysis quantifies a total of 39 flux parameters, comprising anthropogenic and biomass burning sources over large continental regions, soil and lightning emissions of NOx, biogenic emissions of CO and non-methane hydrocarbons, as well as the deposition velocities of both CO and NOx. Comparison between observed, prior and optimized CO mixing ratios at NOAA/CMDL sites shows that the inversion performs well at the northern mid- and high latitudes, and that it is less efficient in the Southern Hemisphere, as expected due to the scarsity of measurements over this part of the globe. The inversion, moreover, brings the model much closer to the measured NO2 columns over all regions. Sensitivity tests show that anthropogenic sources exhibit weak sensitivity to changes of the a priori errors associated to the bottom-up inventory, whereas biomass burning sources are subject to a strong variability. Our best estimate for the 1997 global top-down CO source amounts to 2760 Tg CO. Anthropogenic emissions increase by 28%, in agreement with previous inverse modelling studies, suggesting that the present bottom-up inventories underestimate the anthropogenic CO emissions in the Northern Hemisphere. The magnitude of the optimized NOx global source decreases by 14% with respect to the prior, and amounts to 42.1 Tg N, out of which 22.8 Tg N are due to anthropogenic sources. The NOx emissions increase over Tropical regions, whereas they decrease
Inversion of CO and NOx emissions using the adjoint of the IMAGES model
T. Stavrakou
2004-12-01
Full Text Available We use ground-based observations of CO mixing ratios and vertical column abundances together with tropospheric NO2 columns from the GOME satellite instrument as constraints for improving the global annual emission estimates of CO and NOx for the year 1997. The agreement between concentrations calculated by the global 3-dimensional CTM IMAGES and the observations is optimized using the adjoint modelling technique, which allows to invert for CO and NOx fluxes simultaneously, taking their chemical interactions into account. Our analysis quantifies a total of 39 flux parameters, comprising anthropogenic and biomass burning sources over large continental regions, soil and lightning emissions of NOx, biogenic emissions of CO and non-methane hydrocarbons, as well as the deposition velocities of both CO and NOx. Comparison between observed, prior and optimized CO mixing ratios at NOAA/CMDL sites shows that the inversion performs well at the northern mid- and high latitudes, and that it is less efficient in the Southern Hemisphere, as expected due to the scarsity of measurements over this part of the globe. The inversion, moreover, brings the model much closer to the measured NO2 columns over all regions. Sensitivity tests show that anthropogenic sources exhibit weak sensitivity to changes of the a priori errors associated to the bottom-up inventory, whereas biomass burning sources are subject to a strong variability. Our best estimate for the 1997 global top-down CO source amounts to 2760 Tg CO. Anthropogenic emissions increase by 28%, in agreement with previous inverse modelling studies, suggesting that the present bottom-up inventories underestimate the anthropogenic CO emissions in the Northern Hemisphere. The magnitude of the optimized NOx global source decreases by 14% with respect to the prior, and amounts to 42.1 Tg N, out of which 22.8 Tg N are due to anthropogenic sources. The NOx emissions increase over Tropical regions, whereas they
Grossman, Bernard
1999-01-01
The technical details are summarized below: Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. . An efficient surface parameterization based
Tber, Moulay Hicham; Vidard, Arthur; Dauvergne, Benjamin
2007-01-01
The ocean general circulation model OPA is developed by the LODYC team at Paris VI university. OPA has recently undergone a major rewriting, migrating to FORTRAN95, and its adjoint code needs to be rebuilt. For earlier versions, the adjoint of OPA was written by hand at a high development cost. We use the Automatic Differentiation tool TAPENADE to build mechanicaly the tangent and adjoint codes of OPA. We validate the differentiated codes by comparison with divided differences, and also with an identical twin experiment. We apply state-of-the-art methods to improve the performance of the adjoint code. In particular we implement the Griewank and Walther's binomial checkpointing algorithm which gives us an optimal trade-off between time and memory consumption. We apply a specific strategy to differentiate the iterative linear solver that comes from the implicit time stepping scheme
Mints, R G
1995-01-01
We consider the flux jump instability of the Bean's critical state arising in the flux creep regime in type-II superconductors. We find the flux jump field, $B_j$, that determines the superconducting state stability criterion. We calculate the dependence of $B_j$ on the external magnetic field ramp rate, magnetization experiments the slope of the current-voltage curve in the flux creep regime determines the stability of the Bean's critical state, {\\it i.e.}, the value of $B_j$. We show that a flux jump can be preceded by the magneto-thermal oscillations and find the frequency of these oscillations as a function of $\\dot B_e$.
Self-adjoint elliptic operators with boundary conditions on not closed hypersurfaces
Mantile, Andrea; Posilicano, Andrea; Sini, Mourad
2016-07-01
The theory of self-adjoint extensions of symmetric operators is used to construct self-adjoint realizations of a second-order elliptic differential operator on Rn with linear boundary conditions on (a relatively open part of) a compact hypersurface. Our approach allows to obtain Kreĭn-like resolvent formulae where the reference operator coincides with the "free" operator with domain H2 (Rn); this provides an useful tool for the scattering problem from a hypersurface. Concrete examples of this construction are developed in connection with the standard boundary conditions, Dirichlet, Neumann, Robin, δ and δ‧-type, assigned either on a (n - 1) dimensional compact boundary Γ = ∂ Ω or on a relatively open part Σ ⊂ Γ. Schatten-von Neumann estimates for the difference of the powers of resolvents of the free and the perturbed operators are also proven; these give existence and completeness of the wave operators of the associated scattering systems.
Deconfinement and chiral symmetry restoration in an SU(3) gauge theory with adjoint fermions
We analyze the finite temperature phase diagram of QCD with fermions in the adjoint representation. The simulations performed with four dynamical Majorana fermions show that the deconfinement and chiral phase transitions occur at two distinct temperatures. While the deconfinement transition is first-order at Td we find evidence for a continuous chiral transition at a higher temperature Tc ≅ 8 Td. We observe a rapid change of bulk thermodynamic observables at Td which reflects the increase in the number of degrees of freedom. However, these show little variation at Tc, where the fermion condensate vanishes. We also analyze the potential between static fundamental and adjoint charges in all three phases and extract the corresponding screening masses above Td
Adjoint-based linear analysis in reduced-order thermo-acoustic models
Magri, Luca
2014-01-01
This paper presents the linear theory of adjoint equations as applied to thermo-acoustics. The purpose is to describe the mathematical foundations of adjoint equations for linear sensitivity analysis of thermo-acoustic systems, recently developed by Magri and Juniper (J. Fluid Mech. (2013), vol. 719, pp. 183--202). This method is applied pedagogically to a damped oscillator, for which analytical solutions are available, and then for an electrically heated Rijke tube with a mean-flow temperature discontinuity induced by the compact heat source. Passive devices that most affect the growth rate / frequency of the electrical Rijke-tube system are presented, including a discussion about the effect of modelling the mean-flow temperature discontinuity.
Adjoint-Based a Posteriori Error Estimation for Coupled Time-Dependent Systems
Asner, Liya
2012-01-01
We consider time-dependent parabolic problem s coupled across a common interface which we formulate using a Lagrange multiplier construction and solve by applying a monolithic solution technique. We derive an adjoint-based a posteriori error representation for a quantity of interest given by a linear functional of the solution. We establish the accuracy of our error representation formula through numerical experimentation and investigate the effect of error in the adjoint solution. Crucially, the error representation affords a distinction between temporal and spatial errors and can be used as a basis for a blockwise time-space refinement strategy. Numerical tests illustrate the efficacy of the refinement strategy by capturing the distinctive behavior of a localized traveling wave solution. The saddle point systems considered here are equivalent to those arising in the mortar finite element technique for parabolic problems. © 2012 Society for Industrial and Applied Mathematics.
Weyl theorems for the polluted set of self-adjoint operators in Galerkin approximations
Boulton, Lyonell; Lewin, Mathieu
2010-01-01
Let A be a self-adjoint operator on a separable Hilbert space H and let (L_n) be a sequence of finite dimensional subspaces of the domain of A, approximating H in the large n limit. Denote by A_n the compression of A to L_n. In general the spectrum of A is only a subset of the limit of the spectra of A_n and the latter might differ from the former in a non-trivial "polluted set". In this paper we show that this polluted set is determined by the existence of particular Weyl sequences of singular type. This characterization allows us to identify verifiable conditions on self-adjoint perturbations B, ensuring that the polluted set of B is identical to that of A. The results reported are illustrated by means of several canonical examples and they reveal the many subtleties involved in the systematic study of spectral pollution.
Quantitative photoacoustic tomography using forward and adjoint Monte Carlo models of radiance
Hochuli, Roman; Arridge, Simon; Cox, Ben
2016-01-01
Forward and adjoint Monte Carlo (MC) models of radiance are proposed for use in model-based quantitative photoacoustic tomography. A 2D radiance MC model using a harmonic angular basis is introduced and validated against analytic solutions for the radiance in heterogeneous media. A gradient-based optimisation scheme is then used to recover 2D absorption and scattering coefficients distributions from simulated photoacoustic measurements. It is shown that the functional gradients, which are a challenge to compute efficiently using MC models, can be calculated directly from the coefficients of the harmonic angular basis used in the forward and adjoint models. This work establishes a framework for transport-based quantitative photoacoustic tomography that can fully exploit emerging highly parallel computing architectures.
Active adjoint modeling method in microwave induced thermoacoustic tomography for breast tumor.
Zhu, Xiaozhang; Zhao, Zhiqin; Wang, Jinguo; Chen, Guoping; Liu, Qing Huo
2014-07-01
To improve the model-based inversion performance of microwave induced thermoacoustic tomography for breast tumor imaging, an active adjoint modeling (AAM) method is proposed. It aims to provide a more realistic breast acoustic model used for tumor inversion as the background by actively measuring and reconstructing the structural heterogeneity of human breast environment. It utilizes the reciprocity of acoustic sensors, and adapts the adjoint tomography method from seismic exploration. With the reconstructed acoustic model of breast environment, the performance of model-based inversion method such as time reversal mirror is improved significantly both in contrast and accuracy. To prove the advantage of AAM, a checkerboard pattern model and anatomical realistic breast models have been used in full wave numerical simulations. PMID:24956614
Adjoint Airfoil Optimization of Darrieus-Type Vertical Axis Wind Turbine
Fuchs, Roman; Nordborg, Henrik
2012-11-01
We present the feasibility of using an adjoint solver to optimize the torque of a Darrieus-type vertical axis wind turbine (VAWT). We start with a 2D cross section of a symmetrical airfoil and restrict us to low solidity ratios to minimize blade vortex interactions. The adjoint solver of the ANSYS FLUENT software package computes the sensitivities of airfoil surface forces based on a steady flow field. Hence, we find the torque of a full revolution using a weighted average of the sensitivities at different wind speeds and angles of attack. The weights are computed analytically, and the range of angles of attack is given by the tip speed ratio. Then the airfoil geometry is evolved, and the proposed methodology is evaluated by transient simulations.
Neutrino masses in SU(5) x U(1){sub F} with adjoint flavons
Nardi, Enrico [INFN, Laboratori Nazionali di Frascati, C.P. 13, Frascati (Italy); IFT-UAM/CSIC, Madrid (Spain); Universidad Autonoma de Madrid, Departamento de Fisica Teorica, Madrid (Spain); Restrepo, Diego; Velasquez, Mauricio [Universidad de Antioquia, Instituto de Fisica, Medellin (Colombia)
2012-03-15
We present a SU(5) x U(1){sub F} supersymmetric model for neutrino masses and mixings that implements the seesaw mechanism by means of the heavy SU(2) singlets and triplets states contained in three adjoints of SU(5). We discuss how Abelian U(1){sub F} symmetries can naturally yield non-hierarchical light neutrinos even when the heavy states are strongly hierarchical, and how it can also ensure that R-parity arises as an exact accidental symmetry. By assigning two flavons that break U(1){sub F} to the adjoint representation of SU(5) and assuming universality for all the fundamental couplings, the coefficients of the effective Yukawa and Majorana mass operators become calculable in terms of group theoretical quantities. There is a single free parameter in the model, however, at leading order the structure of the light neutrinos mass matrix is determined in a parameter independent way. (orig.)
Method of tallying adjoint fluence and calculating kinetics parameters in Monte Carlo codes
A method of using iterated fission probability to estimate the adjoint fluence during particles simulation, and using it as the weighting function to calculate kinetics parameters βeff and A in Monte Carlo codes, was introduced in this paper. Implements of this method in continuous energy Monte Carlo code MCNP and multi-group Monte Carlo code MCMG are both elaborated. Verification results show that, with regardless additional computing cost, using this method, the adjoint fluence accounted by MCMG matches well with the result computed by ANISN, and the kinetics parameters calculated by MCNP agree very well with benchmarks. This method is proved to be reliable, and the function of calculating kinetics parameters in Monte Carlo codes is carried out effectively, which could be the basement for Monte Carlo codes' utility in the analysis of nuclear reactors' transient behavior. (authors)
Topology Optimization of Turbulent Fluid Flow with a Sensitive Porosity Adjoint Method (SPAM)
Philippi, B
2015-01-01
A sensitive porosity adjoint method (SPAM) for optimizing the topology of fluid machines has been proposed. A sensitivity function with respect to the porosity has been developed. In the first step of the optimization process, porous media are introduced into the flow regime according to the sensitivity function. Then the optimized porous media are transformed to solid walls. The turbulent flow in porous media is accounted for by a modified eddy-viscosity based turbulence model. Its influence on the adjoint equations is nevertheless neglected, which refers to the so called frozen turbulence assumption. A test case of application in terms of the turbulent rough wall channel flow shows that a considerable reduction of the objective function can be obtained by this method. The transformation from porous media to solid walls may have important effect on the optimization results.
Solving Large-Scale Inverse Magnetostatic Problems using the Adjoint Method
Bruckner, Florian; Wautischer, Gregor; Huber, Christian; Vogler, Christoph; Hinze, Michael; Suess, Dieter
2016-01-01
An efficient algorithm for the reconstruction of the magnetization state within magnetic components is presented. The occurring inverse magnetostatic problem is solved by means of an adjoint approach, based on the Fredkin-Koehler method for the solution of the forward problem. Due to the use of hybrid FEM-BEM coupling combined with matrix compression techniques the resulting algorithm is well suited for large-scale problems. Furthermore the reconstruction of the magnetization state within a permanent magnet is demonstrated.
Adjoint gradient-based approach for aerodynamic optimization of transport aircraft
Ilic, Caslav
2013-01-01
Aerodynamic design of transport aircraft has been steadily improved over past several decades, to the point where today highly-detailed shape control is needed to achieve further improvements. Aircraft manufacturers are therefore increasingly looking into formal optimization methods, driving high-fidelity CFD analysis of finely-parametrized candidate designs. We present an adjoint gradient-based approach for maximizing the aerodynamic performance index relevant to cruise-climb mission segment...
Spectrum of SU(2) gauge theory with two fermions in the adjoint representation
Hietanen, Ari; Rantaharju, Jarno; Rummukainen, Kari; Tuominen, Kimmo
2008-01-01
We present preliminary results of lattice simulations of SU(2) gauge theory with two Wilson fermions in the adjoint representation. This theory has recently attracted considerable attention because it might possess an infrared fixed point (or an almost-fixed-point), and hence be a candidate for a walking technicolor theory. In this work we study the particle spectrum of the theory, and compare it with more familiar spectrum of the theory with SU(2) gauge fields and two flavors of fundamental ...
Serre's reduction of linear systems of partial differential equations with holonomic adjoints
Cluzeau, Thomas; Quadrat, Alban
2010-01-01
Given a linear functional system (e.g., ordinary/partial di erential system, di erential time-delay system, di erence system), Serre's reduction aims at nding an equivalent linear functional system which contains fewer equations and fewer unknowns. The purpose of this paper is to study Serre's reduction of underdetermined linear systems of partial di erential equations with either polynomial, formal power series or analytic coe cients and with holonomic adjoints in the sense of algebraic anal...
Coupling Unification and Dark Matter in a Standard Model Extension with Adjoint Majorana Fermions
Aizawa, Tasuku; Ibe, Masahiro; Kaneta, Kunio
2014-01-01
We revisit an extension of the Standard Model with Majorana fermions in the adjoint representations. There, a precise coupling unification and the good candidate for dark matter (the $SU(2)_L$ triplet fermion) are achieved simultaneously. In particular, we show that the $SU(3)_c$ octet fermion which is required for successful unification can be a good non-thermal source of the triplet fermion dark matter. We also show that the scenario predicts a rather short lifetime of the proton compared w...
Adjoint BFKL at finite coupling: a short-cut from the collinear limit
Basso, Benjamin; Sever, Amit
2014-01-01
In the high energy Regge limit, the six gluons scattering amplitude is controlled by the adjoint BFKL eigenvalue and impact factor. In this paper we determine these two building blocks at any value of the 't Hooft coupling in planar $\\cal{N}$=4 SYM theory. This is achieved by means of analytic continuations from the collinear limit, where similar all loops expressions were recently established. We check our predictions against all available data at weak and strong coupling.
Scattering and self-adjoint extensions of the Aharonov-Bohm Hamiltonian
De Oliveira, Cesar R [Departamento de Matematica-UFSCar, Sao Carlos, Sao Paulo 13560-970 (Brazil); Pereira, Marciano, E-mail: marciano@uepg.b [Departamento de Matematica e EstatIstica-UEPG, Ponta Grossa, Parana 84030-900 (Brazil)
2010-09-03
We consider the Hamiltonian operator associated with planar sections of infinitely long cylindrical solenoids and with a homogeneous magnetic field in their interior. First, in the Sobolev space H{sup 2}, we characterize all generalized boundary conditions on the solenoid border compatible with quantum mechanics, i.e. the boundary conditions, so that the corresponding Hamiltonian operators are self-adjoint. Then we study and compare the scattering of the most usual boundary conditions, that is, Dirichlet, Neumann and Robin.
Edwards, S.; Reuther, J.; Chattot, J. J.
The objective of this paper is to present a control theory approach for the design of airfoils in the presence of viscous compressible flows. A coupled system of the integral boundary layer and the Euler equations is solved to provide rapid flow simulations. An adjoint approach consistent with the complete coupled state equations is employed to obtain the sensitivities needed to drive a numerical optimization algorithm. Design to a target pressure distribution is demonstrated on an RAE 2822 airfoil at transonic speeds.
Some results on the dynamics and transition probabilities for non self-adjoint hamiltonians
We discuss systematically several possible inequivalent ways to describe the dynamics and the transition probabilities of a quantum system when its hamiltonian is not self-adjoint. In order to simplify the treatment, we mainly restrict our analysis to finite dimensional Hilbert spaces. In particular, we propose some experiments which could discriminate between the various possibilities considered in the paper. An example taken from the literature is discussed in detail
The Adjoint Monte Carlo - a viable option for efficient radiotherapy treatment planning
In cancer therapy using collimated beams of photons, the radiation oncologist must determine a set of beams that delivers the required dose to each point in the tumor and minimizes the risk of damage to the healthy tissue and vital organs. Currently, the oncologist determines these beams iteratively, by using a sequence of dose calculations using approximate numerical methods. In this paper, a more accurate and potentially faster approach, based on the Adjoint Monte Carlo method, is presented (authors)
U.S. Geological Survey, Department of the Interior — Methane (CH4) flux is the net rate of methane exchange between an ecosystem and the atmosphere. Data of this variable were generated by the USGS LandCarbon project...
Source attribution of particulate matter pollution over North China with the adjoint method
We quantify the source contributions to surface PM2.5 (fine particulate matter) pollution over North China from January 2013 to 2015 using the GEOS-Chem chemical transport model and its adjoint with improved model horizontal resolution (1/4° × 5/16°) and aqueous-phase chemistry for sulfate production. The adjoint method attributes the PM2.5 pollution to emissions from different source sectors and chemical species at the model resolution. Wintertime surface PM2.5 over Beijing is contributed by emissions of organic carbon (27% of the total source contribution), anthropogenic fine dust (27%), and SO2 (14%), which are mainly from residential and industrial sources, followed by NH3 (13%) primarily from agricultural activities. About half of the Beijing pollution originates from sources outside of the city municipality. Adjoint analyses for other cities in North China all show significant regional pollution transport, supporting a joint regional control policy for effectively mitigating the PM2.5 air pollution. (letter)
First-arrival traveltime tomography for anisotropic media using the adjoint-state method
Waheed, Umair bin
2016-05-27
Traveltime tomography using transmission data has been widely used for static corrections and for obtaining near-surface models for seismic depth imaging. More recently, it is also being used to build initial models for full-waveform inversion. The classic traveltime tomography approach based on ray tracing has difficulties in handling large data sets arising from current seismic acquisition surveys. Some of these difficulties can be addressed using the adjoint-state method, due to its low memory requirement and numerical efficiency. By coupling the gradient computation to nonlinear optimization, it avoids the need for explicit computation of the Fréchet derivative matrix. Furthermore, its cost is equivalent to twice the solution of the forward-modeling problem, irrespective of the size of the input data. The presence of anisotropy in the subsurface has been well established during the past few decades. The improved seismic images obtained by incorporating anisotropy into the seismic processing workflow justify the effort. However, previous literature on the adjoint-state method has only addressed the isotropic approximation of the subsurface. We have extended the adjoint-state technique for first-arrival traveltime tomography to vertical transversely isotropic (VTI) media. Because δ is weakly resolvable from surface seismic alone, we have developed the mathematical framework and procedure to invert for vNMO and η. Our numerical tests on the VTI SEAM model demonstrate the ability of the algorithm to invert for near-surface model parameters and reveal the accuracy achievable by the algorithm.
An adjoint-based approach for finding invariant solutions of Navier–Stokes equations
Farazmand, M.
2016-05-01
We consider the incompressible Navier--Stokes equations with periodic boundary conditions and time-independent forcing. For this type of flow, we derive adjoint equations whose trajectories converge asymptotically to the equilibrium and traveling wave solutions of the Navier--Stokes equations. Using the adjoint equations, arbitrary initial conditions evolve to the vicinity of a (relative) equilibrium at which point a few Newton-type iterations yield the desired (relative) equilibrium solution. We apply this adjoint-based method to a chaotic two-dimensional Kolmogorov flow. A convergence rate of 100% is observed, leading to the discovery of 21 new steady state and traveling wave solutions at Reynolds number Re=40. Some of the new invariant solutions have spatially localized structures that were previously believed to only exist on domains with large aspect ratios. We show that one of the newly found steady state solutions underpins the temporal intermittencies, i.e., high energy dissipation episodes of the flow. More precisely, it is shown that each intermittent episode of a generic turbulent trajectory corresponds to its close passage to this equilibrium solution.
Adjoint-based deviational Monte Carlo methods for phonon transport calculations
Péraud, Jean-Philippe M.; Hadjiconstantinou, Nicolas G.
2015-06-01
In the field of linear transport, adjoint formulations exploit linearity to derive powerful reciprocity relations between a variety of quantities of interest. In this paper, we develop an adjoint formulation of the linearized Boltzmann transport equation for phonon transport. We use this formulation for accelerating deviational Monte Carlo simulations of complex, multiscale problems. Benefits include significant computational savings via direct variance reduction, or by enabling formulations which allow more efficient use of computational resources, such as formulations which provide high resolution in a particular phase-space dimension (e.g., spectral). We show that the proposed adjoint-based methods are particularly well suited to problems involving a wide range of length scales (e.g., nanometers to hundreds of microns) and lead to computational methods that can calculate quantities of interest with a cost that is independent of the system characteristic length scale, thus removing the traditional stiffness of kinetic descriptions. Applications to problems of current interest, such as simulation of transient thermoreflectance experiments or spectrally resolved calculation of the effective thermal conductivity of nanostructured materials, are presented and discussed in detail.
Source attribution of particulate matter pollution over North China with the adjoint method
Zhang, Lin; Liu, Licheng; Zhao, Yuanhong; Gong, Sunling; Zhang, Xiaoye; Henze, Daven K.; Capps, Shannon L.; Fu, Tzung-May; Zhang, Qiang; Wang, Yuxuan
2015-08-01
We quantify the source contributions to surface PM2.5 (fine particulate matter) pollution over North China from January 2013 to 2015 using the GEOS-Chem chemical transport model and its adjoint with improved model horizontal resolution (1/4° × 5/16°) and aqueous-phase chemistry for sulfate production. The adjoint method attributes the PM2.5 pollution to emissions from different source sectors and chemical species at the model resolution. Wintertime surface PM2.5 over Beijing is contributed by emissions of organic carbon (27% of the total source contribution), anthropogenic fine dust (27%), and SO2 (14%), which are mainly from residential and industrial sources, followed by NH3 (13%) primarily from agricultural activities. About half of the Beijing pollution originates from sources outside of the city municipality. Adjoint analyses for other cities in North China all show significant regional pollution transport, supporting a joint regional control policy for effectively mitigating the PM2.5 air pollution.
Neural Network Training by Integration of Adjoint Systems of Equations Forward in Time
Toomarian, Nikzad (Inventor); Barhen, Jacob (Inventor)
1999-01-01
A method and apparatus for supervised neural learning of time dependent trajectories exploits the concepts of adjoint operators to enable computation of the gradient of an objective functional with respect to the various parameters of the network architecture in a highly efficient manner. Specifically. it combines the advantage of dramatic reductions in computational complexity inherent in adjoint methods with the ability to solve two adjoint systems of equations together forward in time. Not only is a large amount of computation and storage saved. but the handling of real-time applications becomes also possible. The invention has been applied it to two examples of representative complexity which have recently been analyzed in the open literature and demonstrated that a circular trajectory can be learned in approximately 200 iterations compared to the 12000 reported in the literature. A figure eight trajectory was achieved in under 500 iterations compared to 20000 previously required. Tbc trajectories computed using our new method are much closer to the target trajectories than was reported in previous studies.
Jacobs, Christian T; Kramer, Stephan C; Funke, Simon W
2016-01-01
Extracting the optimal amount of power from an array of tidal turbines requires an intricate understanding of tidal dynamics and the effects of turbine placement on the local and regional scale flow. Numerical models have contributed significantly towards this understanding, and more recently, adjoint-based modelling has been employed to optimise the positioning of the turbines in an array in an automated way and improve on simple, regular man-made configurations. Adjoint-based optimisation of high-resolution and ideally 3D transient models is generally a very computationally expensive problem. As a result, existing work on the adjoint optimisation of tidal turbine placement has been mostly limited to steady-state simulations in which very high, non-physical values of the background viscosity are required to ensure that a steady-state solution exists. However, such compromises may affect the reliability of the modelled turbines, their wakes and interactions, and thus bring into question the validity of the co...
Flux mapping system for AHWR critical facility
A software for flux mapping system (FMS) for AHWR critical facility has been developed. The system consists of 25 LEU based pulse detectors and associated software. The objective of the FMS is to obtain the flux profiles over the central 5 x 5 lattice locations. For development of flux mapping system it is required to compute the higher harmonics of the diffusion equation. These harmonics (also called λ-modes) are the eigen functions of multi-group diffusion equation. The fundamental mode is found by power iteration method. Apart from the fundamental mode, other higher modes are evaluated by subtraction technique. In the present paper, fundamental eigenvalue and eigenfunction are evaluated by finite difference method. The bi-orthogonality relations between the direct and adjoint eigenvectors are used for this purpose. The reference core for AHWR has been simulated by computer code FINSQR. Two group lattice cell data have been generated using transport theory code WIMSD and its associated 69-group nuclear data library. The reactor core along with the surrounding radial and axial reflector was represented using 3243 mesh points. In total 5 λ-modes and corresponding eigenvalues have been estimated. The computer code FMS has been developed specifically for AHWR/PHWR critical facility. Flux construction at 5 x 5 lattice locations of the core has been achieved by using observed fluxes at 25 detector locations and linear combinations of pre calculated eigen functions. Combining coefficients have been computed by least square method. To validate the code, we have used computer code based on Monte Carlo method for estimation of thermal fluxes at a few discrete locations. The fluxes were estimated using ENDF/B-VI point nuclear data library. (author)
Reentry-Vehicle Shape Optimization Using a Cartesian Adjoint Method and CAD Geometry
Nemec, Marian; Aftosmis, Michael J.
2006-01-01
A DJOINT solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (e.g., geometric parameters that control the shape). Classic aerodynamic applications of gradient-based optimization include the design of cruise configurations for transonic and supersonic flow, as well as the design of high-lift systems. are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric computer-aided design (CAD). In previous work on Cartesian adjoint solvers, Melvin et al. developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the two-dimensional Euler equations using a ghost-cell method to enforce the wall boundary conditions. In Refs. 18 and 19, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm were the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The accuracy of the gradient computation was verified using several three-dimensional test cases, which included design
Asorey, M.; Ibort, A.; Marmo, G.
2015-06-01
The theory of self-adjoint extensions of first- and second-order elliptic differential operators on manifolds with boundary is studied via its most representative instances: Dirac and Laplace operators. The theory is developed by exploiting the geometrical structures attached to them and, by using an adapted Cayley transform on each case, the space {M} of such extensions is shown to have a canonical group composition law structure. The obtained results are compared with von Neumann's theorem characterizing the self-adjoint extensions of densely defined symmetric operators on Hilbert spaces. The 1D case is thoroughly investigated. The geometry of the submanifold of elliptic self-adjoint extensions {M}ellip is studied and it is shown that it is a Lagrangian submanifold of the universal Grassmannian Gr. The topology of {M}ellip is also explored and it is shown that there is a canonical cycle whose dual is the Maslov class of the manifold. Such cycle, called the Cayley surface, plays a relevant role in the study of the phenomena of topology change. Self-adjoint extensions of Laplace operators are discussed in the path integral formalism, identifying a class of them for which both treatments leads to the same results. A theory of dissipative quantum systems is proposed based on this theory and a unitarization theorem for such class of dissipative systems is proved. The theory of self-adjoint extensions with symmetry of Dirac operators is also discussed and a reduction theorem for the self-adjoint elliptic Grassmannian is obtained. Finally, an interpretation of spontaneous symmetry breaking is offered from the point of view of the theory of self-adjoint extensions.
Adjoint sensitivity of global cloud droplet number to aerosol and dynamical parameters
V. A. Karydis
2012-10-01
Full Text Available We present the development of the adjoint of a comprehensive cloud droplet formation parameterization for use in aerosol-cloud-climate interaction studies. The adjoint efficiently and accurately calculates the sensitivity of cloud droplet number concentration (CDNC to all parameterization inputs (e.g., updraft velocity, water uptake coefficient, aerosol number and hygroscopicity with a single execution. The adjoint is then integrated within three dimensional (3-D aerosol modeling frameworks to quantify the sensitivity of CDNC formation globally to each parameter. Sensitivities are computed for year-long executions of the NASA Global Modeling Initiative (GMI Chemical Transport Model (CTM, using wind fields computed with the Goddard Institute for Space Studies (GISS Global Circulation Model (GCM II', and the GEOS-Chem CTM, driven by meteorological input from the Goddard Earth Observing System (GEOS of the NASA Global Modeling and Assimilation Office (GMAO. We find that over polluted (pristine areas, CDNC is more sensitive to updraft velocity and uptake coefficient (aerosol number and hygroscopicity. Over the oceans of the Northern Hemisphere, addition of anthropogenic or biomass burning aerosol is predicted to increase CDNC in contrast to coarse-mode sea salt which tends to decrease CDNC. Over the Southern Oceans, CDNC is most sensitive to sea salt, which is the main aerosol component of the region. Globally, CDNC is predicted to be less sensitive to changes in the hygroscopicity of the aerosols than in their concentration with the exception of dust where CDNC is very sensitive to particle hydrophilicity over arid areas. Regionally, the sensitivities differ considerably between the two frameworks and quantitatively reveal why the models differ considerably in their indirect forcing estimates.
Adjoint Sensitivity Computations for an Embedded-Boundary Cartesian Mesh Method and CAD Geometry
Nemec, Marian; Aftosmis,Michael J.
2006-01-01
Cartesian-mesh methods are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric Computer-Aided Design (CAD) tools. Our goal is to combine the automation capabilities of Cartesian methods with an eficient computation of design sensitivities. We address this issue using the adjoint method, where the computational cost of the design sensitivities, or objective function gradients, is esseutially indepeudent of the number of design variables. In previous work, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm included the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Central to this development is the computation of volume-mesh sensitivities to obtain a reliable approximation of the objective finction gradient. Motivated by the success of mesh-perturbation schemes commonly used in body-fitted unstructured formulations, we propose an approach based on a local linearization of a mesh-perturbation scheme similar to the spring analogy. This approach circumvents most of the difficulties that arise due to non-smooth changes in the cut-cell layer as the boundary shape evolves and provides a consistent approximation tot he exact gradient of the discretized abjective function. A detailed gradient accurace study is presented to verify our approach
Self-adjointness and the Casimir effect with confined quantized spinor matter
Sitenko, Yurii A
2015-01-01
A generalization of the MIT bag boundary condition for spinor matter is proposed basing on the requirement that the Dirac hamiltonian operator be self-adjoint. An influence of a background magnetic field on the vacuum of charged spinor matter confined between two parallel material plates is studied. Employing the most general set of boundary conditions at the plates in the case of the uniform magnetic field directed orthogonally to the plates, we find the pressure from the vacuum onto the plates. In physically plausible situations, the Casimir effect is shown to be repulsive, independently of a choice of boundary conditions and of a distance between the plates.
Infrared fixed point in SU(2) gauge theory with adjoint fermions
We apply Schroedinger-functional techniques to the SU(2) lattice gauge theory with Nf=2 flavors of fermions in the adjoint representation. Our use of hypercubic smearing enables us to work at stronger couplings than did previous studies, before encountering a critical point and a bulk phase boundary. Measurement of the running coupling constant gives evidence of an infrared fixed point g* where 1/g*2=0.20(4)(3). At the fixed point, we find a mass anomalous dimension γm(g*)=0.31(6).
Infrared fixed point in SU(2) gauge theory with adjoint fermions
DeGrand, Thomas; Shamir, Yigal; Svetitsky, Benjamin
2011-01-01
We apply Schrodinger-functional techniques to the SU(2) lattice gauge theory with N_f=2 flavors of fermions in the adjoint representation. Our use of hypercubic smearing enables us to work at stronger couplings than did previous studies, before encountering a critical point and a bulk phase boundary. Measurement of the running coupling constant gives evidence of an infrared fixed point g* where 1/g*^2 = 0.20(4)(3). At the fixed point, we find a mass anomalous dimension gamma_m(g*) = 0.31(6).
Lu, A.H.
1991-09-01
The adjoint method is applied to groundwater flow-mass transport coupled equations in variably saturated media. The sensitivity coefficients derived by this method can be calculated by a single execution for each performance measure regardless of the number of parameters in question. The method provides an efficient and effective way to rank the importance of the parameters, so that data collection can be guided in support of site characterization programs. The developed code will facilitate the sensitivity/uncertainty analysis in both model prediction and model calibration/validation. 13 refs., 1 tab.
The adjoint method is applied to groundwater flow-mass transport coupled equations in variably saturated media. The sensitivity coefficients derived by this method can be calculated by a single execution for each performance measure regardless of the number of parameters in question. The method provides an efficient and effective way to rank the importance of the parameters, so that data collection can be guided in support of site characterization programs. The developed code will facilitate the sensitivity/uncertainty analysis in both model prediction and model calibration/validation. 13 refs., 1 tab
Self-adjointness and the Casimir effect with confined quantized spinor matter
Sitenko, Yurii A.
2016-01-01
A generalization of the MIT bag boundary condition for spinor matter is proposed basing on the requirement that the Dirac hamiltonian operator be self-adjoint. An influence of a background magnetic field on the vacuum of charged spinor matter confined between two parallel material plates is studied. Employing the most general set of boundary conditions at the plates in the case of the uniform magnetic field directed orthogonally to the plates, we find the pressure from the vacuum onto the plates. In physically plausible situations, the Casimir effect is shown to be repulsive, independently of a choice of boundary conditions and of a distance between the plates.
The Adjoint Method Formulation for an Inverse Problem in the Generalized Black-Scholes Model
PIERRE NGNEPIEBA
2006-08-01
Full Text Available A general framework is developed to treat optimal control problems for a generalized Black-Scholes model, which is used for option pricing. The volatility function is retrieved from a set of market observations. The optimal volatility function is found by minimizing the cost functional measuring the discrepancy between the model solution (pricing and the observed market price, via the unconstrained minimization algorithm of the quasi-Newton limited memory type. The gradient is computed via the adjoint method. The effectiveness of the method is demonstrated on an European call option.
Self-adjoint Extensions of Schrödinger Operators with ?-magnetic Fields on Riemannian Manifolds
T. Mine
2010-01-01
Full Text Available We consider the magnetic Schr¨odinger operator on a Riemannian manifold M. We assume the magnetic field is given by the sum of a regular field and the Dirac δ measures supported on a discrete set Γ in M. We give a complete characterization of the self-adjoint extensions of the minimal operator, in terms of the boundary conditions. The result is an extension of the former results by Dabrowski-Šťoviček and Exner-Šťoviček-Vytřas.
On Maximal Abelian Self-adjoint Subalgebras of Factors of Type Ⅱ1
Li Guang WANG
2005-01-01
In this note, we show that if (N) is a proper subfactor of a factor (M) of type Ⅱ1 with finite Jones index, then there is a maximal abelian self-adjoint subalgebra (masa) (A) of (N) that is not a masa in (M). Popa showed that there is a proper subfactor (R)O of the hyperfinite type Ⅱ1 factor (R) such that each masa in (R)O is also a masa in (R). We shall give a detailed proof of Popa's result.