Consistent Adjoint Driven Importance Sampling using Space, Energy and Angle
Peplow, Douglas E.; Mosher, Scott W; Evans, Thomas M
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.
Characterization of a Source Importance Function in an Accelerator-Driven System
Kim, Yonghee; Park, Won Seok; Park, Chang Kue
2003-07-15
An importance function of the external spallation neutrons in an accelerator-driven system (ADS) has been introduced and characterized to address the source multiplication in a subcritical blanket. For a model ADS problem with a central external source, the source importance function is evaluated with a neutron transport code system. For a homogeneous core, essential characteristics of the importance are identified from the viewpoint of spatial distributions and energy dependency, etc. The importance function is evaluated for two different beam tube diameters, and its dependency on the buffer thickness is also addressed. In order to assess the impact of the power distribution on the importance function, a heterogeneous core is considered, and its importance function is evaluated. The analyses show that the peak importance occurs in the inner fuel blanket zone, not in the central source region, and the neutron importance in a high-energy regime, above 7 to 20 MeV, is high and increases with the energy. Also, the effects of a neutron absorber on the source importance are studied, and it is found that the source importance could be drastically reduced by surrounding the source with a strong neutron absorber such as B{sub 4}C. In addition, the source importance function is compared with the conventional {lambda}-mode adjoint flux, which is used as an importance function of fission neutrons in critical reactors. The comparison reveals that the inhomogeneous source importance function could be quite similar to the homogeneous {lambda}-mode adjoint flux in both spatial and spectral distributions for a wide range of subcriticality.
Introduction to Adjoint Models
NASA Technical Reports Server (NTRS)
Errico, Ronald M.
2015-01-01
In this lecture, some fundamentals of adjoint models will be described. This includes a basic derivation of tangent linear and corresponding adjoint models from a parent nonlinear model, the interpretation of adjoint-derived sensitivity fields, a description of methods of automatic differentiation, and the use of adjoint models to solve various optimization problems, including singular vectors. Concluding remarks will attempt to correct common misconceptions about adjoint models and their utilization.
FAST TRACK COMMUNICATION Quasi self-adjoint nonlinear wave equations
NASA Astrophysics Data System (ADS)
Ibragimov, N. H.; Torrisi, M.; Tracinà, R.
2010-11-01
In this paper we generalize the classification of self-adjoint second-order linear partial differential equation to a family of nonlinear wave equations with two independent variables. We find a class of quasi self-adjoint nonlinear equations which includes the self-adjoint linear equations as a particular case. The property of a differential equation to be quasi self-adjoint is important, e.g. for constructing conservation laws associated with symmetries of the differential equation.
Importance-driven feature enhancement in volume visualization.
Viola, Ivan; Kanitsar, Armin; Gröller, M Eduard
2005-01-01
This paper presents importance-driven feature enhancement as a technique for the automatic generation of cut-away and ghosted views out of volumetric data. The presented focus+context approach removes or suppresses less important parts of a scene to reveal more important underlying information. However, less important parts are fully visible in those regions, where important visual information is not lost, i.e., more relevant features are not occluded. Features within the volumetric data are first classified according to a new dimension, denoted as object importance. This property determines which structures should be readily discernible and which structures are less important. Next, for each feature, various representations (levels of sparseness) from a dense to a sparse depiction are defined. Levels of sparseness define a spectrum of optical properties or rendering styles. The resulting image is generated by ray-casting and combining the intersected features proportional to their importance (importance compositing). The paper includes an extended discussion on several possible schemes for levels of sparseness specification. Furthermore, different approaches to importance compositing are treated. PMID:16138551
Diagnositcs With Adjoint Modelling
NASA Astrophysics Data System (ADS)
Blessing, S.; Fraedrich, K.; Kirk, E.; Lunkeit, F.
The potential usefulness of an adjoint primitive equations global atmospheric circu- lation model for climate diagnostics is demonstrated in a feasibility study. A daily NAO-type index is calculated as one-point correlation of the 300 hPa streamfunction anomaly. By application of the adjoint model we diagnose its temperature forcing on short timescales in terms of spatial temperature sensitivity patterns at different time lags, which, in a first order approximation, induce growth of the index. The dynamical relevance of these sensitivity patterns is confirmed by lag-correlating the index time series and the projection time series of the model temperature on these sensitivity patterns.
Adjoint affine fusion and tadpoles
NASA Astrophysics Data System (ADS)
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-Based Uncertainty Quantification with MCNP
Seifried, Jeffrey E.
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.
Adjoint methods for external beam inverse treatment planning
NASA Astrophysics Data System (ADS)
Kowalok, Michael E.
Forward and adjoint radiation transport methods may both be used to determine the dosimetric relationship between source parameters and voxel elements of a phantom. Forward methods consider one specific tuple of source parameters and calculate the response in all voxels of interest. This response is often cast as the dose delivered per unit source-weight. Adjoint transport methods, conversely, consider one particular voxel and calculate the response of that voxel in relation to all possible source parameters. In this regard, adjoint methods provide an "adjoint function" in addition to a dose value. Although the dose is for a single voxel only, the adjoint function illustrates the source parameters, (e.g. beam positions and directions) that are most important to delivering the dose to that voxel. In this regard, adjoint methods of analysis lend themselves in a natural way to optimization problems and perturbation studies. This work investigates the utility of adjoint analytic methods for treatment planning and for Monte Carlo dose calculations. Various methods for implementing this approach are discussed, along with their strengths and weaknesses. The complementary nature of adjoint and forward techniques is illustrated and exploited. Also, several features of the Monte Carlo codes MCNP and MCNPX are reviewed for treatment planning applications.
MCNP: Multigroup/adjoint capabilities
Wagner, J.C.; Redmond, E.L. II; Palmtag, S.P.; Hendricks, J.S.
1994-04-01
This report discusses various aspects related to the use and validity of the general purpose Monte Carlo code MCNP for multigroup/adjoint calculations. The increased desire to perform comparisons between Monte Carlo and deterministic codes, along with the ever-present desire to increase the efficiency of large MCNP calculations has produced a greater user demand for the multigroup/adjoint capabilities. To more fully utilize these capabilities, we review the applications of the Monte Carlo multigroup/adjoint method, describe how to generate multigroup cross sections for MCNP with the auxiliary CRSRD code, describe how to use the multigroup/adjoint capability in MCNP, and provide examples and results indicating the effectiveness and validity of the MCNP multigroup/adjoint treatment. This information should assist users in taking advantage of the MCNP multigroup/adjoint capabilities.
Adjoint Function: Physical Basis of Variational & Perturbation Theory in Transport
Energy Science and Technology Software Center (ESTSC)
2009-07-27
Version 00 Dr. J.D. Lewins has now released the following legacy book for free distribution: Importance: The Adjoint Function: The Physical Basis of Variational and Perturbation Theory in Transport and Diffusion Problems, North-Holland Publishing Company - Amsterdam, 582 pages, 1966 Introduction: Continuous Systems and the Variational Principle 1. The Fundamental Variational Principle 2. The Importance Function 3. Adjoint Equations 4. Variational Methods 5. Perturbation and Iterative Methods 6. Non-Linear Theory
NASA Astrophysics Data System (ADS)
Bozdag, Ebru; Lefebvre, Matthieu; Lei, Wenjie; Peter, Daniel; Smith, James; Komatitsch, Dimitri; Tromp, Jeroen
2015-04-01
We will present our initial results of global adjoint tomography based on 3D seismic wave simulations which is one of the most challenging examples in seismology in terms of intense computational requirements and vast amount of high-quality seismic data that can potentially be assimilated in inversions. Using a spectral-element method, we incorporate full 3D wave propagation in seismic tomography by running synthetic seismograms and adjoint simulations to compute exact sensitivity kernels in realistic 3D background models. We run our global simulations on the Oak Ridge National Laboratory's Cray XK7 "Titan" system taking advantage of the GPU version of the SPECFEM3D_GLOBE package. We have started iterations with initially selected 253 earthquakes within the magnitude range of 5.5 < Mw < 7.0 and numerical simulations having resolution down to ~27 s to invert for a transversely isotropic crust and mantle model using a non-linear conjugate gradient algorithm. The measurements are currently based on frequency-dependent traveltime misfits. We use both minor- and major-arc body and surface waves by running 200 min simulations where inversions are performed with more than 2.6 million measurements. Our initial results after 12 iterations already indicate several prominent features such as enhanced slab (e.g., Hellenic, Japan, Bismarck, Sandwich), plume/hotspot (e.g., the Pacific superplume, Caroline, Yellowstone, Hawaii) images, etc. To improve the resolution and ray coverage, particularly in the lower mantle, our aim is to increase the resolution of numerical simulations first going down to ~17 s and then to ~9 s to incorporate high-frequency body waves in inversions. While keeping track of the progress and illumination of features in our models with a limited data set, we work towards to assimilate all available data in inversions from all seismic networks and earthquakes in the global CMT catalogue.
NASA Astrophysics Data System (ADS)
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.
Towards Global Adjoint Tomography
NASA Astrophysics Data System (ADS)
Bozdag, E.; Zhu, H.; Peter, D. B.; Tromp, J.
2012-12-01
Seismic tomography is at a stage where we can harness entire seismograms using the opportunities offered by advances in numerical wave propagation solvers and high-performance computing. Adjoint methods provide an efficient way for incorporating full nonlinearity of wave propagation and 3D Fréchet kernels in iterative seismic inversions which have so far given promising results at continental and regional scales. Our goal is to take adjoint tomography forward to image the entire planet. Using an iterative conjugate gradient scheme, we initially set the aim to obtain a global crustal and mantle model with confined transverse isotropy in the upper mantle. We have started with around 255 global CMT events having moment magnitudes between 5.8 and 7, and used GSN stations as well as some local networks such as USArray, European stations etc. Prior to the structure inversion, we reinvert global CMT solutions by computing Green functions in our 3D reference model to take into account effects of crustal variations on source parameters. Using the advantages of numerical simulations, our strategy is to invert crustal and mantle structure together to avoid any bias introduced into upper-mantle images due to "crustal corrections", which are commonly used in classical tomography. 3D simulations dramatically increase the usable amount of data so that, with the current earthquake-station setup, we perform each iteration with more than two million measurements. Multi-resolution smoothing based on ray density is applied to the gradient to better deal with the imperfect source-station distribution on the globe and extract more information underneath regions with dense ray coverage and vice versa. Similar to frequency domain approach, we reduce nonlinearities by starting from long periods and gradually increase the frequency content of data after successive model updates. To simplify the problem, we primarily focus on the elastic structure and therefore our measurements are based on
Mesh-free adjoint methods for nonlinear filters
NASA Astrophysics Data System (ADS)
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.
Learning a trajectory using adjoint functions and teacher forcing
NASA Technical Reports Server (NTRS)
Toomarian, Nikzad B.; Barhen, Jacob
1992-01-01
A new methodology for faster supervised temporal learning in nonlinear neural networks is presented which builds upon the concept of adjoint operators to allow fast computation of the gradients of an error functional with respect to all parameters of the neural architecture, and exploits the concept of teacher forcing to incorporate information on the desired output into the activation dynamics. The importance of the initial or final time conditions for the adjoint equations is discussed. A new algorithm is presented in which the adjoint equations are solved simultaneously (i.e., forward in time) with the activation dynamics of the neural network. We also indicate how teacher forcing can be modulated in time as learning proceeds. The results obtained show that the learning time is reduced by one to two orders of magnitude with respect to previously published results, while trajectory tracking is significantly improved. The proposed methodology makes hardware implementation of temporal learning attractive for real-time applications.
The relative importance of the wind-driven and tidal circulations in Malacca Strait
NASA Astrophysics Data System (ADS)
Chen, Haoliang; Malanotte-Rizzoli, Paola; Koh, Tieh-Yong; Song, Guiting
2014-10-01
The Malacca Strait is traditionally treated as a typical tidally-driven channel with the wind-driven and other components considered negligible. However, the strait is frequently affected by intense tropical weather events distorting the background monsoon winds. The variable winds can create large wind-stress curl at the surface level. To answer the question of how significant the wind-driven circulation is to the total circulation, numerical simulations are carried out by isolating or superimposing the different driving mechanisms. Comparison of the time series at selected points reveals that the winds significantly affect the tidal currents in different ways in the northern and southern strait. In the northern wide strait, the tidal current is enhanced while in the southern narrow channel it is weakened. Experiments with uniform water depth confirm that the weakening is mainly due to the interaction among tidal current, wind-driven current and bathymetry in the southern strait. Spectral analysis of the currents in the whole MS quantifies that the wind-driven current energy is more significant in the northern channel than in the southern one. Furthermore, winds with high intensity and large wind-stress curl can produce an eddy as large as the northern channel width which significantly distorts the tidal circulation especially during the neap tide. Vorticity analysis shows that the eddy in the northern Malacca Strait is purely wind-driven. Our study highlights that the wind stress, which has been ignored in previous studies in this region, is an important driver of the circulation in the Malacca Strait even when tidal forcing is strong.
Self-adjointness of deformed unbounded operators
Much, Albert
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.
Numerical Computation of Sensitivities and the Adjoint Approach
NASA Technical Reports Server (NTRS)
Lewis, Robert Michael
1997-01-01
We discuss the numerical computation of sensitivities via the adjoint approach in optimization problems governed by differential equations. We focus on the adjoint problem in its weak form. We show how one can avoid some of the problems with the adjoint approach, such as deriving suitable boundary conditions for the adjoint equation. We discuss the convergence of numerical approximations of the costate computed via the weak form of the adjoint problem and show the significance for the discrete adjoint problem.
ANISORROPIA: the adjoint of the aerosol thermodynamic model ISORROPIA
NASA Astrophysics Data System (ADS)
Capps, S. L.; Henze, D. K.; Hakami, A.; Russell, A. G.; Nenes, A.
2011-08-01
We present the development of ANISORROPIA, the discrete adjoint of the ISORROPIA thermodynamic equilibrium model that treats the Na+-SO42--HSO4--NH4+-NO3--Cl--H2O aerosol system, and we demonstrate its sensitivity analysis capabilities. ANISORROPIA calculates sensitivities of an inorganic species in aerosol or gas phase with respect to the total concentrations of each species present with only a two-fold increase in computational time over the forward model execution. Due to the highly nonlinear and discontinuous solution surface of ISORROPIA, evaluation of the adjoint required a new, complex-variable version of the the model, which determines first-order sensitivities with machine precision and avoids cancellation errors arising from finite difference calculations. The adjoint is verified over an atmospherically relevant range of concentrations, temperature, and relative humidity. We apply ANISORROPIA to recent field campaign results from Atlanta, GA, USA, and Mexico City, Mexico, to characterize the inorganic aerosol sensitivities of these distinct urban air masses. The variability in the relationship between PM2.5 mass and precursor concentrations shown has important implications for air quality and climate. ANISORROPIA enables efficient elucidation of aerosol concentration dependence on aerosol precursor emissions in the context of atmospheric chemical transport model adjoints.
Adjoints and Low-rank Covariance Representation
NASA Technical Reports Server (NTRS)
Tippett, Michael K.; Cohn, Stephen E.
2000-01-01
Quantitative measures of the uncertainty of Earth System estimates can be as important as the estimates themselves. Second moments of estimation errors are described by the covariance matrix, whose direct calculation is impractical when the number of degrees of freedom of the system state is large. Ensemble and reduced-state approaches to prediction and data assimilation replace full estimation error covariance matrices by low-rank approximations. The appropriateness of such approximations depends on the spectrum of the full error covariance matrix, whose calculation is also often impractical. Here we examine the situation where the error covariance is a linear transformation of a forcing error covariance. We use operator norms and adjoints to relate the appropriateness of low-rank representations to the conditioning of this transformation. The analysis is used to investigate low-rank representations of the steady-state response to random forcing of an idealized discrete-time dynamical system.
Adjoint Error Estimation for Linear Advection
Connors, J M; Banks, J W; Hittinger, J A; Woodward, C S
2011-03-30
An a posteriori error formula is described when a statistical measurement of the solution to a hyperbolic conservation law in 1D is estimated by finite volume approximations. This is accomplished using adjoint error estimation. In contrast to previously studied methods, the adjoint problem is divorced from the finite volume method used to approximate the forward solution variables. An exact error formula and computable error estimate are derived based on an abstractly defined approximation of the adjoint solution. This framework allows the error to be computed to an arbitrary accuracy given a sufficiently well resolved approximation of the adjoint solution. The accuracy of the computable error estimate provably satisfies an a priori error bound for sufficiently smooth solutions of the forward and adjoint problems. The theory does not currently account for discontinuities. Computational examples are provided that show support of the theory for smooth solutions. The application to problems with discontinuities is also investigated computationally.
Hekmat, Mohamad Hamed; Mirzaei, Masoud
2015-01-01
In the present research, we tried to improve the performance of the lattice Boltzmann (LB) -based adjoint approach by utilizing the mesoscopic inherent of the LB method. In this regard, two macroscopic discrete adjoint (MADA) and microscopic discrete adjoint (MIDA) approaches are used to answer the following two challenging questions. Is it possible to extend the concept of the macroscopic and microscopic variables of the flow field to the corresponding adjoint ones? Further, similar to the conservative laws in the LB method, is it possible to find the comparable conservation equations in the adjoint approach? If so, then a definite framework, similar to that used in the flow solution by the LB method, can be employed in the flow sensitivity analysis by the MIDA approach. This achievement can decrease the implementation cost and coding efforts of the MIDA method in complicated sensitivity analysis problems. First, the MADA and MIDA equations are extracted based on the LB method using the duality viewpoint. Meanwhile, using an elementary case, inverse design of a two-dimensional unsteady Poiseuille flow in a periodic channel with constant body forces, the procedure of analytical evaluation of the adjoint variables is described. The numerical results show that similar correlations between the distribution functions can be seen between the corresponding adjoint ones. Besides, the results are promising, emphasizing the flow field adjoint variables can be evaluated via the adjoint distribution functions. Finally, the adjoint conservative laws are introduced. PMID:25679735
Southern California Adjoint Source Inversions
NASA Astrophysics Data System (ADS)
Tromp, J.; Kim, Y.
2007-12-01
Southern California Centroid-Moment Tensor (CMT) solutions with 9 components (6 moment tensor elements, latitude, longitude, and depth) are sought to minimize a misfit function computed from waveform differences. The gradient of a misfit function is obtained based upon two numerical simulations for each earthquake: one forward calculation for the southern California model, and an adjoint calculation that uses time-reversed signals at the receivers. Conjugate gradient and square-root variable metric methods are used to iteratively improve the earthquake source model while reducing the misfit function. The square-root variable metric algorithm has the advantage of providing a direct approximation to the posterior covariance operator. We test the inversion procedure by perturbing each component of the CMT solution, and see how the algorithm converges. Finally, we demonstrate full inversion capabilities using data for real Southern California earthquakes.
Adjoint calculations for multiple scattering of Compton and Rayleigh effects
NASA Astrophysics Data System (ADS)
Fernández, J. E.; Sumini, M.
1992-08-01
As is well known, the experimental determination of the Compton profile requires a particular geometry with a scattering angle close to π. That situation involves a narrow multiple-scattering spectrum that overlaps the Compton peak, making it difficult to analyze the different contributions to the profile. We show how the solution of the adjoint problem can help in devising more useful experimental configurations, giving, through its classical "importance" meaning, a formally clear picture of the whole problem.
Adjoint optimization of natural convection problems: differentially heated cavity
NASA Astrophysics Data System (ADS)
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 Algorithm for CAD-Based Shape Optimization Using a Cartesian Method
NASA Technical Reports Server (NTRS)
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
Receptivity in parallel flows: An adjoint approach
NASA Technical Reports Server (NTRS)
Hill, D. Christopher
1993-01-01
Linear receptivity studies in parallel flows are aimed at understanding how external forcing couples to the natural unstable motions which a flow can support. The vibrating ribbon problem models the original Schubauer and Skramstad boundary layer experiment and represents the classic boundary layer receptivity problem. The process by which disturbances are initiated in convectively-unstable jets and shear layers has also received attention. Gaster was the first to handle the boundary layer analysis with the recognition that spatial modes, rather than temporal modes, were relevant when studying convectively-unstable flows that are driven by a time-harmonic source. The amplitude of the least stable spatial mode, far downstream of the source, is related to the source strength by a coupling coefficient. The determination of this coefficient is at the heart of this type of linear receptivity study. The first objective of the present study was to determine whether the various wave number derivative factors, appearing in the coupling coefficients for linear receptivity problems, could be reexpressed in a simpler form involving adjoint eigensolutions. Secondly, it was hoped that the general nature of this simplification could be shown; indeed, a rather elegant characterization of the receptivity properties of spatial instabilities does emerge. The analysis is quite distinct from the usual Fourier-inversion procedures, although a detailed knowledge of the spectrum of the Orr-Sommerfeld equation is still required. Since the cylinder wake analysis proved very useful in addressing control considerations, the final objective was to provide a foundation upon which boundary layer control theory may be developed.
ANISORROPIA: the adjoint of the aerosol thermodynamic model ISORROPIA
NASA Astrophysics Data System (ADS)
Capps, S. L.; Henze, D. K.; Hakami, A.; Russell, A. G.; Nenes, A.
2012-01-01
We present the development of ANISORROPIA, the discrete adjoint of the ISORROPIA thermodynamic equilibrium model that treats the Na+-SO42-- HSO4--NH4+ -NO3--Cl--H2O aerosol system, and we demonstrate its sensitivity analysis capabilities. ANISORROPIA calculates sensitivities of an inorganic species in aerosol or gas phase with respect to the total concentrations of each species present with less than a two-fold increase in computational time over the concentration calculations. Due to the highly nonlinear and discontinuous solution surface of ISORROPIA, evaluation of the adjoint required a new, complex-variable version of the model, which determines first-order sensitivities with machine precision and avoids cancellation errors arising from finite difference calculations. The adjoint is verified over an atmospherically relevant range of concentrations, temperature, and relative humidity. We apply ANISORROPIA to recent field campaign results from Atlanta, GA, USA, and Mexico City, Mexico, to characterize the inorganic aerosol sensitivities of these distinct urban air masses. The variability in the relationship between fine mode inorganic aerosol mass and precursor concentrations shown has important implications for air quality and climate.
Double-difference adjoint seismic tomography
NASA Astrophysics Data System (ADS)
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
Worley, B.A.; Pin, F.G.; Horwedel, J.E.; Oblow, E.M.
1989-05-01
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
NASA Technical Reports Server (NTRS)
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.
ERIC Educational Resources Information Center
King, Jean A.; Rohmer-Hirt, Johnna A.
2011-01-01
From the 1980s to the present, educational accountability in the United States has grown dramatically. Such accountability in U.S. school districts, although driven primarily by external demands, has internal manifestations as well. The chapter traces the historical development of internal evaluation in American school districts, then highlights…
Adjoint sensitivity studies of loop current and eddy shedding in the Gulf of Mexico
NASA Astrophysics Data System (ADS)
Gopalakrishnan, Ganesh; Cornuelle, Bruce D.; Hoteit, Ibrahim
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
Martien, Philip T; Harley, Robert A; Cacuci, Dan G
2006-04-15
Photochemical air pollution forms when emissions of nitrogen oxides (NO(x)) and volatile organic compounds (VOC) react in the atmosphere in the presence of sunlight. The goal of applying three-dimensional photochemical air quality models is usually to conduct sensitivity analysis: for example, to predict changes in an ozone response due to changes in NO(x) and VOC emissions or other model data. Forward sensitivity analysis methods are best suited to investigating sensitivities of many model responses to changes in a few inputs or parameters. Here we develop a continuous adjoint model and demonstrate an adjoint sensitivity analysis procedure that is well-suited to the complementary case of determining sensitivity of a small number of model responses to many parameters. Sensitivities generated using the adjoint method agree with those generated using other methods. Compared to the forward method, the adjoint method had large disk storage requirements but was more efficient in terms of computer processor time for receptor-based investigations focused on a single response at a specified site and time. The adjoint method also generates sensitivity apportionment fields, which reveal when and where model data are important to the target response. PMID:16683606
Assessing the Impact of Observations on Numerical Weather Forecasts Using the Adjoint Method
NASA Technical Reports Server (NTRS)
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
Important requirements for RF generators for Accelerator-Driven Transmutation Technologies (ADTT)
Lynch, M.T.; Tallerico, P.J.; Lawrence, G.P.
1994-09-01
All Accelerator-Driven Transmutation applications require very large amounts of RF Power. For example, one version of a Plutonium burning system requires an 800-MeV, 80-mA, proton accelerator running at 100% duty factor. This accelerator requires approximately 110-MW of continuous RF power if one assumes only 10% reserve power for control of the accelerator fields. In fact, to minimize beam spill, the RF controls may need as much as 15 to 20% of reserve power. In addition, unlike an electron accelerator in which the beam is relativistic, a failed RF station can disturb the synchronism of the beam, possibly shutting down the entire accelerator. These issues and more lead to a set of requirements for the RF generators which are stringent, and in some cases, conflicting. In this paper, we will describe the issues and requirements, and outline a plan for RF generator development to meet the needs of the Accelerator-Driven Transmutation Technologies. The key issues which will be discussed include: operating efficiency, operating linearity, effect on the input power grid, bandwidth, gain, reliability, operating voltage, and operating current.
Fully automatic adjoints: a robust and efficient mechanism for generating adjoint ocean models
NASA Astrophysics Data System (ADS)
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
Uematsu, Mikio; Kurosawa, Masahiko
2005-01-01
A generalised and convenient skyshine dose analysis method has been developed based on forward-adjoint folding technique. In the method, the air penetration data were prepared by performing an adjoint DOT3.5 calculation with cylindrical air-over-ground geometry having an adjoint point source (importance of unit flux to dose rate at detection point) in the centre. The accuracy of the present method was certified by comparing with DOT3.5 forward calculation. The adjoint flux data can be used as generalised radiation skyshine data for all sorts of nuclear facilities. Moreover, the present method supplies plenty of energy-angular dependent contribution flux data, which will be useful for detailed shielding design of facilities. PMID:16604693
On the Importance of the Equation of State for the Neutrino-driven Supernova Explosion Mechanism
NASA Astrophysics Data System (ADS)
Suwa, Yudai; Takiwaki, Tomoya; Kotake, Kei; Fischer, Tobias; Liebendörfer, Matthias; Sato, Katsuhiko
2013-02-01
By implementing the widely used equations of state (EOS) from Lattimer & Swesty (LS) and H. Shen et al. (SHEN) in core-collapse supernova simulations, we explore possible impacts of these EOS on the post-bounce dynamics prior to the onset of neutrino-driven explosions. Our spherically symmetric (1D) and axially symmetric (2D) models are based on neutrino radiation hydrodynamics including spectral transport, which is solved by the isotropic diffusion source approximation. We confirm that in 1D simulations neutrino-driven explosions cannot be obtained for any of the employed EOS. Impacts of the EOS on the post-bounce hydrodynamics are more clearly visible in 2D simulations. In 2D models of a 15 M ⊙ progenitor using the LS EOS, the stalled bounce shock expands to increasingly larger radii, which is not the case when using the SHEN EOS. Keeping in mind that the omission of the energy drain by heavy-lepton neutrinos in the present scheme could facilitate explosions, we find that 2D models of an 11.2 M ⊙ progenitor produce neutrino-driven explosions for all the EOS under investigation. Models using the LS EOS are slightly more energetic compared with those with the SHEN EOS. The more efficient neutrino heating in the LS models coincides with a higher electron antineutrino luminosity and a larger mass that is enclosed within the gain region. The models based on the LS EOS also show a more vigorous and aspherical downflow of accreting matter to the surface of the protoneutron star (PNS). The accretion pattern is essential for the production and strength of outgoing pressure waves, which can push in turn the shock to larger radii and provide more favorable conditions for the explosion. Based on our models, we investigate several diagnostic indicators of the explosion that have been suggested in the literature, e.g., the amplitude of the standing accretion shock instability mode, the mass-weighted average entropy in the gain region, the PNS radius, the antesonic
Bhogal, Moninder S; Lanyon-Hogg, Thomas; Johnston, Katherine A; Warriner, Stuart L; Baker, Alison
2016-01-29
Peroxisomes are vital metabolic organelles found in almost all eukaryotic organisms, and they rely exclusively on import of their matrix protein content from the cytosol. In vitro import of proteins into isolated peroxisomal fractions has provided a wealth of knowledge on the import process. However, the common method of protease protection garnered no information on the import of an N-terminally truncated PEX5 (PEX5C) receptor construct or peroxisomal malate dehydrogenase 1 (pMDH1) cargo protein into sunflower peroxisomes because of high degrees of protease susceptibility or resistance, respectively. Here we present a means for analysis of in vitro import through a covalent biotin label transfer and employ this method to the import of PEX5C. Label transfer demonstrates that the PEX5C construct is monomeric under the conditions of the import assay. This technique was capable of identifying the PEX5-PEX14 interaction as the first interaction of the import process through competition experiments. Labeling of the peroxisomal protein import machinery by PEX5C demonstrated that this interaction was independent of added cargo protein, and, strikingly, the interaction between PEX5C and the import machinery was shown to be ATP-dependent. These important mechanistic insights highlight the power of label transfer in studying interactions, rather than proteins, of interest and demonstrate that this technique should be applied to future studies of peroxisomal in vitro import. PMID:26567336
Advances in Global Adjoint Tomography -- Massive Data Assimilation
NASA Astrophysics Data System (ADS)
Ruan, Y.; Lei, W.; Bozdag, E.; Lefebvre, M. P.; Smith, J. A.; Krischer, L.; Tromp, J.
2015-12-01
the conventional multi-taper method to obtain better frequency-dependent measurements of surface-wave phase and amplitude anomalies, and therefore more accurate adjoint sources, which are particularly important for anelastic tomography. We present a summary of these data culling and processing procedures for global adjoint tomography.
Wing planform optimization via an adjoint method
NASA Astrophysics Data System (ADS)
Leoviriyakit, Kasidit
This dissertation focuses on the problem of wing planform optimization for transonic aircraft based on flow simulation using Computational Fluid Dynamics (CFD) combined with an adjoint-gradient based numerical optimization procedure. The adjoint method, traditionally used for wing section design has been extended to cover planform variations and to compute the sensitivities of the structural weight of both the wing section and planform variations. The two relevant disciplines accounted for are the aerodynamics and structural weight. A simplified structural weight model is used for the optimization. Results of a variety of long range transports indicate that significant improvement in both aerodynamics and structures can be achieved simultaneously. The proof-of-concept optimal results indicate large improvements for both drag and structural weight. The work is an "enabling step" towards a realistic automated wing designed by a computer.
Development of CO2 inversion system based on the adjoint of the global coupled transport model
NASA Astrophysics Data System (ADS)
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
NASA Astrophysics Data System (ADS)
Bora, Tanujjal; Zoepfl, David; Dutta, Joydeep
2016-05-01
Herein we explore the role of localized plasmonic heat generated by resonantly excited gold (Au) NPs on visible light driven photocatalysis process. Au NPs are deposited on the surface of vertically aligned zinc oxide nanorods (ZnO NRs). The localized heat generated by Au NPs under 532 nm continuous laser excitation (SPR excitation) was experimentally probed using Raman spectroscopy by following the phonon modes of ZnO. Under the resonant excitation the temperature at the surface of the Au-ZnO NRs reaches up to about 300 °C, resulting in almost 6 times higher apparent quantum yield (AQY) for photocatalytic degradation of methylene blue (MB) compared to the bare ZnO NRs. Under solar light irradiation the Au-ZnO NRs demonstrated visible light photocatalytic activity twice that of what was achieved with bare ZnO NRs, while significantly reduced the activation energy required for the photocatalytic reactions allowing the reactions to occur at a faster rate.
Bora, Tanujjal; Zoepfl, David; Dutta, Joydeep
2016-01-01
Herein we explore the role of localized plasmonic heat generated by resonantly excited gold (Au) NPs on visible light driven photocatalysis process. Au NPs are deposited on the surface of vertically aligned zinc oxide nanorods (ZnO NRs). The localized heat generated by Au NPs under 532 nm continuous laser excitation (SPR excitation) was experimentally probed using Raman spectroscopy by following the phonon modes of ZnO. Under the resonant excitation the temperature at the surface of the Au-ZnO NRs reaches up to about 300 °C, resulting in almost 6 times higher apparent quantum yield (AQY) for photocatalytic degradation of methylene blue (MB) compared to the bare ZnO NRs. Under solar light irradiation the Au-ZnO NRs demonstrated visible light photocatalytic activity twice that of what was achieved with bare ZnO NRs, while significantly reduced the activation energy required for the photocatalytic reactions allowing the reactions to occur at a faster rate. PMID:27242172
Bora, Tanujjal; Zoepfl, David; Dutta, Joydeep
2016-01-01
Herein we explore the role of localized plasmonic heat generated by resonantly excited gold (Au) NPs on visible light driven photocatalysis process. Au NPs are deposited on the surface of vertically aligned zinc oxide nanorods (ZnO NRs). The localized heat generated by Au NPs under 532 nm continuous laser excitation (SPR excitation) was experimentally probed using Raman spectroscopy by following the phonon modes of ZnO. Under the resonant excitation the temperature at the surface of the Au-ZnO NRs reaches up to about 300 °C, resulting in almost 6 times higher apparent quantum yield (AQY) for photocatalytic degradation of methylene blue (MB) compared to the bare ZnO NRs. Under solar light irradiation the Au-ZnO NRs demonstrated visible light photocatalytic activity twice that of what was achieved with bare ZnO NRs, while significantly reduced the activation energy required for the photocatalytic reactions allowing the reactions to occur at a faster rate. PMID:27242172
Dual of QCD with one adjoint fermion
Mojaza, Matin; Nardecchia, Marco; Pica, Claudio; Sannino, Francesco
2011-03-15
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 anomalous dimension of the Dirac fermion mass operator to be less than one in the conformal window.
Supersymmetric descendants of self-adjointly extended quantum mechanical Hamiltonians
NASA Astrophysics Data System (ADS)
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.
Argument-Driven Inquiry to Promote the Understanding of Important Concepts & Practices in Biology
ERIC Educational Resources Information Center
Sampson, Victor; Gleim, Leeanne
2009-01-01
Inquiry is an integral part of the teaching and learning of science. However, many science teachers are unsure of how to promote and support inquiry in the classroom or how to design lessons that engage students in inquiry in a way that improves students' understanding of important concepts and practices in biology. In this article, the authors…
Generalized uncertainty principle and self-adjoint operators
Balasubramanian, Venkat; Das, Saurya; Vagenas, Elias C.
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.
Laser-driven isomerization of HCN → HNC: the importance of rotational excitation.
Sun, Zhaopeng; Zheng, Yujun
2015-03-26
We report a time-dependent quantum wave packet theory, which is employed to interpret the isomerization dynamics of HCN molecules induced by an intense picosecond infrared laser field. Considering the molecular rotational degrees of the freedom, the wave functions are expanded in terms of molecular rotational bases. Our full-dimensional quantum model includes the full Coriolis coupling in the molecular kinetic energy Hamiltonian and dipole approximation in interaction terms. The numerical results show that the field-induced molecule rotational excitation plays an important role in the isomerization dynamical process. Some phenomena appear such as two-step two-photon absorption and highly oscillatory structure in rotational state distributions. The centrifugal sudden (CS) approximation calculation is also carried out and compared in this work; it is shown that the Coriolis couplings may lead to a significant decrease in the isomerization rate but highly enhanced molecular rotational excitation. PMID:25746130
Importance of a flexible hinge near the motor domain in kinesin-driven motility.
Grummt, M; Woehlke, G; Henningsen, U; Fuchs, S; Schleicher, M; Schliwa, M
1998-01-01
Conventional kinesin is a molecular motor consisting of an N-terminal catalytic motor domain, an extended stalk and a small globular C-terminus. Whereas the structure and function of the catalytic motor domain has been investigated, little is known about the function of domains outside the globular head. A short coiled-coil region adjacent to the motor domain, termed the neck, is known to be important for dimerization and may be required for kinesin processivity. We now provide evidence that a helix-disrupting hinge region (hinge 1) that separates the neck from the first extended coiled-coil of the stalk plays an essential role in basic motor activity. A fast fungal kinesin from Syncephalastrum racemosum was used for these studies. Deletion, substitution by a coiled-coil and truncation of the hinge 1 region all reduce motor speed and uncouple ATP turnover from gliding velocity. Insertion of hinge 1 regions from two conventional kinesins, Nkin and DmKHC, fully restores motor activity, whereas insertion of putative flexible linkers of other proteins does not, suggesting that hinge 1 regions of conventional kinesins can functionally replace each other. We suggest that this region is essential for kinesin movement in its promotion of chemo-mechanical coupling of the two heads and therefore the functional motor domain should be redefined to include not only the catalytic head but also the adjacent neck and hinge 1 domains. PMID:9755154
Using adjoint-based optimization to study wing flexibility in flapping flight
NASA Astrophysics Data System (ADS)
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.
Adjoint sensitivity analysis of an ultrawideband antenna
Stephanson, M B; White, D A
2011-07-28
The frequency domain finite element method using H(curl)-conforming finite elements is a robust technique for full-wave analysis of antennas. As computers become more powerful, it is becoming feasible to not only predict antenna performance, but also to compute sensitivity of antenna performance with respect to multiple parameters. This sensitivity information can then be used for optimization of the design or specification of manufacturing tolerances. In this paper we review the Adjoint Method for sensitivity calculation, and apply it to the problem of optimizing a Ultrawideband antenna.
The importance of the transmission flux in evaluating the preheat effect in x-ray driven ablation
Li, Liling; Zhang, Lu; Zheng, Jianhua; Qing, Bo; Zhang, Jiyan; Kuang, Longyu; Li, Hang; Jiang, Shaoen
2015-02-15
In x-ray driven ablation, the preheat effect is caused by the high energy x-rays that pass through the ablator. Thus, the transmission flux can be used to characterize preheat effect in a certain degree. With the radiation temperature being 200 eV, the transmission flux and preheat temperature of pure polymer (CH) have been studied by using the one-dimensional multi-group radiation hydrodynamic code MULTI-1D. By studying the spectrum of the transmitted x-rays, it is found that the energy of the transmitted x-rays is in the range of 2–5 keV for pure CH ablator. This is of importance for selecting a dopant for CH ablator. We also calculated both the preheat temperature of CH near the surface of thick target (47.66 μm) and the transmission flux of a thinner target (38.66 μm). It is found that the more transmission flux leads to the higher preheat temperature. Preheat effect of graded Si-doped CH targets with different doped concentrations has also been studied. The results are consistent with this phenomenon. By analyzing the relationship between the transmission flux and the preheat temperature, we have presented a novel method to evaluate preheat effect in x-ray driven ablation.
Improved Adjoint-Operator Learning For A Neural Network
NASA Technical Reports Server (NTRS)
Toomarian, Nikzad; Barhen, Jacob
1995-01-01
Improved method of adjoint-operator learning reduces amount of computation and associated computational memory needed to make electronic neural network learn temporally varying pattern (e.g., to recognize moving object in image) in real time. Method extension of method described in "Adjoint-Operator Learning for a Neural Network" (NPO-18352).
Supersymmetric descendants of self-adjointly extended quantum mechanical Hamiltonians
Al-Hashimi, M.H.; Salman, M.; Shalaby, A.; Wiese, U.-J.
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.
The compressible adjoint equations in geodynamics: equations and numerical assessment
NASA Astrophysics Data System (ADS)
Ghelichkhan, Siavash; Bunge, Hans-Peter
2016-04-01
The adjoint method is a powerful means to obtain gradient information in a mantle convection model relative to past flow structure. While the adjoint equations in geodynamics have been derived for the conservation equations of mantle flow in their incompressible form, the applicability of this approximation to Earth is limited, because density increases by almost a factor of two from the surface to the Core Mantle Boundary. Here we introduce the compressible adjoint equations for the conservation equations in the anelastic-liquid approximation. Our derivation applies an operator formulation in Hilbert spaces, to connect to recent work in seismology (Fichtner et al (2006)) and geodynamics (Horbach et al (2014)), where the approach was used to derive the adjoint equations for the wave equation and incompressible mantle flow. We present numerical tests of the newly derived equations based on twin experiments, focusing on three simulations. A first, termed Compressible, assumes the compressible forward and adjoint equations, and represents the consistent means of including compressibility effects. A second, termed Mixed, applies the compressible forward equation, but ignores compressibility effects in the adjoint equations, where the incompressible equations are used instead. A third simulation, termed Incompressible, neglects compressibility effects entirely in the forward and adjoint equations relative to the reference twin. The compressible and mixed formulations successfully restore earlier mantle flow structure, while the incompressible formulation yields noticeable artifacts. Our results suggest the use of a compressible formulation, when applying the adjoint method to seismically derived mantle heterogeneity structure.
Adjoint-based uncertainty quantification and sensitivity analysis for reactor depletion calculations
NASA Astrophysics Data System (ADS)
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.
Fisch, N.J.
1985-03-01
A plasma in contact with an external source of power, especially a source that interacts specifically with high-velocity electrons, exhibits transport properties, such as conductivity, different from those of an isolated plasma near thermal equilibrium. This is true even when the bulk of the particles in the driven plasma are near thermal equilibrium. To describe the driven plasma we derive an adjoint equation to the inhomogeneous, linearized, dynamic Boltzmann equation. The Green's functions for a variety of plasma responses can then be generated. It is possible to modify the Chapman-Enskog expansion in order to incorporate the response functions derived here.
Aerodynamic design optimization by using a continuous adjoint method
NASA Astrophysics Data System (ADS)
Luo, JiaQi; Xiong, JunTao; Liu, Feng
2014-07-01
This paper presents the fundamentals of a continuous adjoint method and the applications of this method to the aerodynamic design optimization of both external and internal flows. General formulation of the continuous adjoint equations and the corresponding boundary conditions are derived. With the adjoint method, the complete gradient information needed in the design optimization can be obtained by solving the governing flow equations and the corresponding adjoint equations only once for each cost function, regardless of the number of design parameters. An inverse design of airfoil is firstly performed to study the accuracy of the adjoint gradient and the effectiveness of the adjoint method as an inverse design method. Then the method is used to perform a series of single and multiple point design optimization problems involving the drag reduction of airfoil, wing, and wing-body configuration, and the aerodynamic performance improvement of turbine and compressor blade rows. The results demonstrate that the continuous adjoint method can efficiently and significantly improve the aerodynamic performance of the design in a shape optimization problem.
Adjoint active surfaces for localization and imaging.
Cook, Daniel A; Mueller, Martin Fritz; Fedele, Francesco; Yezzi, Anthony J
2015-01-01
This paper addresses the problem of localizing and segmenting regions embedded within a surrounding medium by characterizing their boundaries, as opposed to imaging the entirety of the volume. Active surfaces are used to directly reconstruct the shape of the region of interest. We describe the procedure for finding the optimal surface, which is computed iteratively via gradient descent that exploits the sensitivity of an error minimization functional to changes of the active surface. In doing so, we introduce the adjoint model to compute the sensitivity, and in this respect, the method shares common ground with several other disciplines, such as optimal control. Finally, we illustrate the proposed active surface technique in the framework of wave propagation governed by the scalar Helmholtz equation. Potential applications include electromagnetics, acoustics, geophysics, nondestructive testing, and medical imaging. PMID:25438311
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
NASA Astrophysics Data System (ADS)
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.
GPU-Accelerated Adjoint Algorithmic Differentiation
Gremse, Felix; Höfter, Andreas; Razik, Lukas; Kiessling, Fabian; Naumann, Uwe
2015-01-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
Self-adjointness and conservation laws of difference equations
NASA Astrophysics Data System (ADS)
Peng, Linyu
2015-06-01
A general theorem on conservation laws for arbitrary difference equations is proved. The theorem is based on an introduction of an adjoint system related with a given difference system, and it does not require the existence of a difference Lagrangian. It is proved that the system, combined by the original system and its adjoint system, is governed by a variational principle, which inherits all symmetries of the original system. Noether's theorem can then be applied. With some special techniques, e.g. self-adjointness properties, this allows us to obtain conservation laws for difference equations, which are not necessary governed by Lagrangian formalisms.
Dupree, S. A.
1980-06-01
The use of adjoint techniques to determine the interaction of externally incident collimated beams of particles with cylindrical targets is a convenient means of examining a class of problems important in radiation transport studies. The theory relevant to such applications is derived, and a simple example involving a fissioning target is discussed. Results from both discrete ordinates and Monte Carlo transport-code calculations are presented, and comparisons are made with results obtained from forward calculations. The accuracy of the discrete ordinates adjoint results depends on the order of angular quadrature used in the calculation. Reasonable accuracy by using EQN quadratures can be expected from order S/sub 16/ or higher.
Sensitivity of Lumped Constraints Using the Adjoint Method
NASA Technical Reports Server (NTRS)
Akgun, Mehmet A.; Haftka, Raphael T.; Wu, K. Chauncey; Walsh, Joanne L.
1999-01-01
Adjoint sensitivity calculation of stress, buckling and displacement constraints may be much less expensive than direct sensitivity calculation when the number of load cases is large. Adjoint stress and displacement sensitivities are available in the literature. Expressions for local buckling sensitivity of isotropic plate elements are derived in this study. Computational efficiency of the adjoint method is sensitive to the number of constraints and, therefore, the method benefits from constraint lumping. A continuum version of the Kreisselmeier-Steinhauser (KS) function is chosen to lump constraints. The adjoint and direct methods are compared for three examples: a truss structure, a simple HSCT wing model, and a large HSCT model. These sensitivity derivatives are then used in optimization.
Surface wave sensitivity: mode summation versus adjoint SEM
NASA Astrophysics Data System (ADS)
Zhou, Ying; Liu, Qinya; Tromp, Jeroen
2011-12-01
We compare finite-frequency phase and amplitude sensitivity kernels calculated based on frequency-domain surface wave mode summation and a time-domain adjoint method. The adjoint calculations involve a forward wavefield generated by an earthquake and an adjoint wavefield generated at a seismic receiver. We determine adjoint sources corresponding to frequency-dependent phase and amplitude measurements made using a multitaper technique, which may be applied to any single-taper measurement, including box car windowing. We calculate phase and amplitude sensitivity kernels using an adjoint method based on wave propagation simulations using a spectral element method (SEM). Sensitivity kernels calculated using the adjoint SEM are in good agreement with kernels calculated based on mode summation. In general, the adjoint SEM is more computationally expensive than mode summation in global studies. The advantage of the adjoint SEM lies in the calculation of sensitivity kernels in 3-D earth models. We compare surface wave sensitivity kernels computed in 1-D and 3-D reference earth models and show that (1) lateral wave speed heterogeneities may affect the geometry and amplitude of surface wave sensitivity; (2) sensitivity kernels of long-period surface waves calculated in 1-D model PREM and 3-D models S20RTS+CRUST2.0 and FFSW1+CRUST2.0 do not show significant differences, indicating that the use of a 1-D reference model is adequate in global inversions of long-period surface waves (periods of 50 s and longer); and (3) the differences become significant for short-period Love waves when mode coupling is sensitive to large differences in reference crustal structure. Finally, we show that sensitivity kernels in anelastic earth models may be calculated in purely elastic earth models provided physical dispersion is properly accounted for.
Universal Racah matrices and adjoint knot polynomials: Arborescent knots
NASA Astrophysics Data System (ADS)
Mironov, A.; Morozov, A.
2016-04-01
By now it is well established that the quantum dimensions of descendants of the adjoint representation can be described in a universal form, independent of a particular family of simple Lie algebras. The Rosso-Jones formula then implies a universal description of the adjoint knot polynomials for torus knots, which in particular unifies the HOMFLY (SUN) and Kauffman (SON) polynomials. For E8 the adjoint representation is also fundamental. We suggest to extend the universality from the dimensions to the Racah matrices and this immediately produces a unified description of the adjoint knot polynomials for all arborescent (double-fat) knots, including twist, 2-bridge and pretzel. Technically we develop together the universality and the "eigenvalue conjecture", which expresses the Racah and mixing matrices through the eigenvalues of the quantum R-matrix, and for dealing with the adjoint polynomials one has to extend it to the previously unknown 6 × 6 case. The adjoint polynomials do not distinguish between mutants and therefore are not very efficient in knot theory, however, universal polynomials in higher representations can probably be better in this respect.
Unsteady adjoint of a gas turbine inlet guide vane
NASA Astrophysics Data System (ADS)
Talnikar, Chaitanya; Wang, Qiqi
2015-11-01
Unsteady fluid flow simulations like large eddy simulation have been shown to be crucial in accurately predicting heat transfer in turbomachinery applications like transonic flow over an inlet guide vane. To compute sensitivities of aerothermal objectives for a vane with respect to design parameters an unsteady adjoint is required. In this talk we present unsteady adjoint solutions for a vane from VKI using pressure loss and heat transfer over the vane surface as the objectives. The boundary layer on the suction side near the trailing edge of the vane is turbulent and this poses a challenge for an adjoint solver. The chaotic dynamics cause the adjoint solution to diverge exponentially to infinity from that region when simulated backwards in time. The prospect of adding artificial viscosity to the adjoint equations to dampen the adjoint fields is investigated. Results for the vane from simulations performed on the Titan supercomputer will be shown and the effect of the additional viscosity on the accuracy of the sensitivities will be discussed.
Baryogenesis via leptogenesis in adjoint SU(5)
Blanchet, Steve; Fileviez Perez, Pavel E-mail: fileviez@physics.wisc.edu
2008-08-15
The possibility of explaining the baryon asymmetry in the Universe through the leptogenesis mechanism in the context of adjoint SU(5) is investigated. In this model neutrino masses are generated through the type I and type III seesaw mechanisms, and the field responsible for the type III seesaw, called {rho}{sub 3}, generates the B-L asymmetry needed to satisfy the observed value of the baryon asymmetry in the Universe. We find that the CP asymmetry originates only from the vertex correction, since the self-energy contribution is not present. When neutrino masses have a normal hierarchy, successful leptogenesis is possible for 10{sup 11} GeV{approx}
Big Data Challenges in Global Seismic 'Adjoint Tomography' (Invited)
NASA Astrophysics Data System (ADS)
Tromp, J.; Bozdag, E.; Krischer, L.; Lefebvre, M.; Lei, W.; Smith, J.
2013-12-01
The challenge of imaging Earth's interior on a global scale is closely linked to the challenge of handling large data sets. The related iterative workflow involves five distinct phases, namely, 1) data gathering and culling, 2) synthetic seismogram calculations, 3) pre-processing (time-series analysis and time-window selection), 4) data assimilation and adjoint calculations, 5) post-processing (pre-conditioning, regularization, model update). In order to implement this workflow on modern high-performance computing systems, a new seismic data format is being developed. The Adaptable Seismic Data Format (ASDF) is designed to replace currently used data formats with a more flexible format that allows for fast parallel I/O. The metadata is divided into abstract categories, such as "source" and "receiver", along with provenance information for complete reproducibility. The structure of ASDF is designed keeping in mind three distinct applications: earthquake seismology, seismic interferometry, and exploration seismology. Existing time-series analysis tool kits, such as SAC and ObsPy, can be easily interfaced with ASDF so that seismologists can use robust, previously developed software packages. ASDF accommodates an automated, efficient workflow for global adjoint tomography. Manually managing the large number of simulations associated with the workflow can rapidly become a burden, especially with increasing numbers of earthquakes and stations. Therefore, it is of importance to investigate the possibility of automating the entire workflow. Scientific Workflow Management Software (SWfMS) allows users to execute workflows almost routinely. SWfMS provides additional advantages. In particular, it is possible to group independent simulations in a single job to fit the available computational resources. They also give a basic level of fault resilience as the workflow can be resumed at the correct state preceding a failure. Some of the best candidates for our particular workflow
Determining scaling laws from geodynamic simulations using adjoint gradients.
NASA Astrophysics Data System (ADS)
Reuber, Georg; Kaus, Boris; Popov, Anton
2016-04-01
Whereas significant progress has been made in modelling of lithospheric and crustal scale processes in recent years, it often remains a challenge to understand which of the many model parameters is of key importance for a particular simulation. Determining this is usually done by manually changing the model input parameters and performing new simulations. For a few cases, such as for crustal-scale folding instabilities (with viscous rheologies, e.g. [1]) or for Rayleigh-Taylor instabilities, one can use existing scaling laws to obtain such insights. Yet, for a more general case, it is not straightforward to do this (apart from running many simulations). Here, we test a different approach which computes gradients of the model parameters using adjoint based methods, which has the advantage that we can test the influence of an independent number of parameters on the system by computing and analysing the covariance matrix and the gradient of the parameter space. This method might give us the chance to get insights on which parameters affect for example subduction processes and how strong the system depends on their influence. [1] Fernandez, N., & Kaus, B. J. (2014). Fold interaction and wavelength selection in 3D models of multilayer detachment folding. Tectonophysics, 632, 199-217.
Multigrid methods for bifurcation problems: The self adjoint case
NASA Technical Reports Server (NTRS)
Taasan, Shlomo
1987-01-01
This paper deals with multigrid methods for computational problems that arise in the theory of bifurcation and is restricted to the self adjoint case. The basic problem is to solve for arcs of solutions, a task that is done successfully with an arc length continuation method. Other important issues are, for example, detecting and locating singular points as part of the continuation process, switching branches at bifurcation points, etc. Multigrid methods have been applied to continuation problems. These methods work well at regular points and at limit points, while they may encounter difficulties in the vicinity of bifurcation points. A new continuation method that is very efficient also near bifurcation points is presented here. The other issues mentioned above are also treated very efficiently with appropriate multigrid algorithms. For example, it is shown that limit points and bifurcation points can be solved for directly by a multigrid algorithm. Moreover, the algorithms presented here solve the corresponding problems in just a few work units (about 10 or less), where a work unit is the work involved in one local relaxation on the finest grid.
Adjoint tomography of the Middle East
NASA Astrophysics Data System (ADS)
Peter, D. B.; Savage, B.; Rodgers, A. J.; Tromp, J.
2010-12-01
Improvements in nuclear explosion monitoring require refined seismic models of the target region. In our study, we focus on the Middle East, spanning a region from Turkey to the west and West India to the east. This area represents a complex geologic and tectonic setting with sparse seismic data coverage. This has lead to diverging interpretations of crustal and underlying upper-mantle structure by different research groups, complicating seismic monitoring of the Middle East at regional distances. We evaluated an initial 3D seismic model of this region by computing full waveforms for several regional earthquakes by a spectral-element method. We measure traveltime and multitaper phase shifts between observed broadband data and synthetic seismograms for distinct seismic phases within selected time windows using a recently developed automated measurement algorithm. Based on the remaining misfits, we setup an iterative inversion procedure for a fully numerical 3D seismic tomography approach. In order to improve the initial 3D seismic model, the sensitivity to seismic structure of the traveltime and multitaper phase measurements for all available seismic network recordings is computed. As this represents a computationally very intensive task, we take advantage of a fully numerical adjoint approach by using the efficient software package SPECFEM3D_GLOBE on a dedicated cluster. We show examples of such sensitivity kernels for different seismic events and use them in a steepest descent approach to update the 3D seismic model, starting at longer periods between 60 s and up to 200 s and moving towards shorter periods of 11 s. We highlight various improvements in the initial seismic structure during the iterations in order to better fit regional seismic waveforms in the Middle East.
Adjoint tomography of the Middle East
NASA Astrophysics Data System (ADS)
Peter, D. B.; Savage, B.; Rodgers, A.; Morency, C.; Tromp, J.
2011-12-01
Improvements in nuclear explosion monitoring require refined seismic models of the target region. In our study, we focus on the Middle East, spanning a region from Turkey to the west and West India to the east. This area represents a complex geologic and tectonic setting with sparse seismic data coverage. This has lead to diverging interpretations of crustal and underlying upper-mantle structure by different research groups, complicating seismic monitoring of the Middle East at regional distances. We evaluated an initial 3D seismic model of this region by computing full waveforms for several regional earthquakes based on a spectral-element method. We measure traveltime and multitaper phase differences between observed broadband data and synthetic seismograms for distinct seismic phases within selected time windows using a recently developed automated measurement algorithm. Based on the remaining misfits, we setup an iterative inversion procedure for a fully numerical 3D seismic tomography approach. In order to improve the initial 3D seismic model, sensitivity to seismic structures of traveltime and multitaper phase measurements for all available seismic network recordings is computed. As this represents a computationally very intensive task, we take advantage of a fully numerical adjoint approach by using the efficient software package SPECFEM3D_GLOBE on a dedicated cluster. We show examples of such sensitivity kernels for different seismic events. All these `event kernels' are then summed, smoothed and further used in a preconditioned conjugate-gradient approach. Thus we iteratively update the 3D seismic model, starting at longer periods between 60~s and up to 150~s and moving towards shorter periods of 11~s. We highlight various improvements in the initial seismic structure during the iterations in order to better fit regional seismic waveforms in the Middle East.
Adjoint simulation of stream depletion due to aquifer pumping.
Neupauer, Roseanna M; Griebling, Scott A
2012-01-01
If an aquifer is hydraulically connected to an adjacent stream, a pumping well operating in the aquifer will draw some water from aquifer storage and some water from the stream, causing stream depletion. Several analytical, semi-analytical, and numerical approaches have been developed to estimate stream depletion due to pumping. These approaches are effective if the well location is known. If a new well is to be installed, it may be desirable to install the well at a location where stream depletion is minimal. If several possible locations are considered for the location of a new well, stream depletion would have to be estimated for all possible well locations, which can be computationally inefficient. The adjoint approach for estimating stream depletion is a more efficient alternative because with one simulation of the adjoint model, stream depletion can be estimated for pumping at a well at any location. We derive the adjoint equations for a coupled system with a confined aquifer, an overlying unconfined aquifer, and a river that is hydraulically connected to the unconfined aquifer. We assume that the stage in the river is known, and is independent of the stream depletion, consistent with the assumptions of the MODFLOW river package. We describe how the adjoint equations can be solved using MODFLOW. In an illustrative example, we show that for this scenario, the adjoint approach is as accurate as standard forward numerical simulation methods, and requires substantially less computational effort. PMID:22182421
Adjoint-Based Methodology for Time-Dependent Optimization
NASA Technical Reports Server (NTRS)
Yamaleev, N. K.; Diskin, B.; Nielsen, E. J.
2008-01-01
This paper presents a discrete adjoint method for a broad class of time-dependent optimization problems. The time-dependent adjoint equations are derived in terms of the discrete residual of an arbitrary finite volume scheme which approximates unsteady conservation law equations. Although only the 2-D unsteady Euler equations are considered in the present analysis, this time-dependent adjoint method is applicable to the 3-D unsteady Reynolds-averaged Navier-Stokes equations with minor modifications. The discrete adjoint operators involving the derivatives of the discrete residual and the cost functional with respect to the flow variables are computed using a complex-variable approach, which provides discrete consistency and drastically reduces the implementation and debugging cycle. The implementation of the time-dependent adjoint method is validated by comparing the sensitivity derivative with that obtained by forward mode differentiation. Our numerical results show that O(10) optimization iterations of the steepest descent method are needed to reduce the objective functional by 3-6 orders of magnitude for test problems considered.
Investigating Sensitivity to Saharan Dust in Tropical Cyclone Formation Using Nasa's Adjoint Model
NASA Technical Reports Server (NTRS)
Holdaway, Daniel
2015-01-01
As tropical cyclones develop from easterly waves coming of the coast of Africa they interact with dust from the Sahara desert. There is a long standing debate over whether this dust inhibits or advances the developing storm and how much influence it has. Dust can surround the storm and absorb incoming solar radiation, cooling the air below. As a result an energy source for the system is potentially diminished, inhibiting growth of the storm. Alternatively dust may interact with clouds through micro-physical processes, for example by causing more moisture to condense, potentially increasing the strength. As a result of climate change, concentrations and amount of dust in the atmosphere will likely change. It it is important to properly understand its effect on tropical storm formation. The adjoint of an atmospheric general circulation model provides a very powerful tool for investigating sensitivity to initial conditions. The National Aeronautics and Space Administration (NASA) has recently developed an adjoint version of the Goddard Earth Observing System version 5 (GEOS-5) dynamical core, convection scheme, cloud model and radiation schemes. This is extended so that the interaction between dust and radiation is also accounted for in the adjoint model. This provides a framework for examining the sensitivity to dust in the initial conditions. Specifically the set up allows for an investigation into the extent to which dust affects cyclone strength through absorption of radiation. In this work we investigate the validity of using an adjoint model for examining sensitivity to dust in hurricane formation. We present sensitivity results for a number of systems that developed during the Atlantic hurricane season of 2006. During this period there was a significant outbreak of Saharan dust and it is has been argued that this outbreak was responsible for the relatively calm season. This period was also covered by an extensive observation campaign. It is shown that the
Adjoint sensitivity analysis of plasmonic structures using the FDTD method.
Zhang, Yu; Ahmed, Osman S; Bakr, Mohamed H
2014-05-15
We present an adjoint variable method for estimating the sensitivities of arbitrary responses with respect to the parameters of dispersive discontinuities in nanoplasmonic devices. Our theory is formulated in terms of the electric field components at the vicinity of perturbed discontinuities. The adjoint sensitivities are computed using at most one extra finite-difference time-domain (FDTD) simulation regardless of the number of parameters. Our approach is illustrated through the sensitivity analysis of an add-drop coupler consisting of a square ring resonator between two parallel waveguides. The computed adjoint sensitivities of the scattering parameters are compared with those obtained using the accurate but computationally expensive central finite difference approach. PMID:24978258
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.
Sonic Boom Mitigation Through Aircraft Design and Adjoint Methodology
NASA Technical Reports Server (NTRS)
Rallabhandi, Siriam K.; Diskin, Boris; Nielsen, Eric J.
2012-01-01
This paper presents a novel approach to design of the supersonic aircraft outer mold line (OML) by optimizing the A-weighted loudness of sonic boom signature predicted on the ground. The optimization process uses the sensitivity information obtained by coupling the discrete adjoint formulations for the augmented Burgers Equation and Computational Fluid Dynamics (CFD) equations. This coupled formulation links the loudness of the ground boom signature to the aircraft geometry thus allowing efficient shape optimization for the purpose of minimizing the impact of loudness. The accuracy of the adjoint-based sensitivities is verified against sensitivities obtained using an independent complex-variable approach. The adjoint based optimization methodology is applied to a configuration previously optimized using alternative state of the art optimization methods and produces additional loudness reduction. The results of the optimizations are reported and discussed.
Higginbottom, Nick; Knigge, Christian; Matthews, James H.; Proga, Daniel; Long, Knox S.; Sim, Stuart A.
2014-07-01
Accretion disk winds are thought to produce many of the characteristic features seen in the spectra of active galactic nuclei (AGNs) and quasi-stellar objects (QSOs). These outflows also represent a natural form of feedback between the central supermassive black hole and its host galaxy. The mechanism for driving this mass loss remains unknown, although radiation pressure mediated by spectral lines is a leading candidate. Here, we calculate the ionization state of, and emergent spectra for, the hydrodynamic simulation of a line-driven disk wind previously presented by Proga and Kallman. To achieve this, we carry out a comprehensive Monte Carlo simulation of the radiative transfer through, and energy exchange within, the predicted outflow. We find that the wind is much more ionized than originally estimated. This is in part because it is much more difficult to shield any wind regions effectively when the outflow itself is allowed to reprocess and redirect ionizing photons. As a result, the calculated spectrum that would be observed from this particular outflow solution would not contain the ultraviolet spectral lines that are observed in many AGN/QSOs. Furthermore, the wind is so highly ionized that line driving would not actually be efficient. This does not necessarily mean that line-driven winds are not viable. However, our work does illustrate that in order to arrive at a self-consistent model of line-driven disk winds in AGN/QSO, it will be critical to include a more detailed treatment of radiative transfer and ionization in the next generation of hydrodynamic simulations.
Reentry-Vehicle Shape Optimization Using a Cartesian Adjoint Method and CAD Geometry
NASA Technical Reports Server (NTRS)
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
Approximations of Strongly Continuous Families of Unbounded Self-Adjoint Operators
NASA Astrophysics Data System (ADS)
Ben-Artzi, Jonathan; Holding, Thomas
2016-07-01
The problem of approximating the discrete spectra of families of self-adjoint operators that are merely strongly continuous is addressed. It is well-known that the spectrum need not vary continuously (as a set) under strong perturbations. However, it is shown that under an additional compactness assumption the spectrum does vary continuously, and a family of symmetric finite-dimensional approximations is constructed. An important feature of these approximations is that they are valid for the entire family uniformly. An application of this result to the study of plasma instabilities is illustrated.
Approximations of Strongly Continuous Families of Unbounded Self-Adjoint Operators
NASA Astrophysics Data System (ADS)
Ben-Artzi, Jonathan; Holding, Thomas
2016-05-01
The problem of approximating the discrete spectra of families of self-adjoint operators that are merely strongly continuous is addressed. It is well-known that the spectrum need not vary continuously (as a set) under strong perturbations. However, it is shown that under an additional compactness assumption the spectrum does vary continuously, and a family of symmetric finite-dimensional approximations is constructed. An important feature of these approximations is that they are valid for the entire family uniformly. An application of this result to the study of plasma instabilities is illustrated.
Ginsparg-Wilson relation on a fuzzy 2-sphere for adjoint matter
Aoki, Hajime
2010-10-15
We formulate a Ginsparg-Wilson relation on a fuzzy 2-sphere for matter in the adjoint representation of the gauge group. Because of the Ginsparg-Wilson relation, an index theorem is satisfied. Our formulation is applicable to topologically nontrivial configurations as monopoles. It gives a solid basis for obtaining chiral fermions, which are an important ingredient of the standard model, from matrix model formulations of the superstring theory, such as the IIB matrix model, by considering topological configurations in the extra dimensions. We finally discuss whether this mechanism really works.
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...
Adjoint operator approach to shape design for internal incompressible flows
NASA Technical Reports Server (NTRS)
Cabuk, H.; Sung, C.-H.; Modi, V.
1991-01-01
The problem of determining the profile of a channel or duct that provides the maximum static pressure rise is solved. Incompressible, laminar flow governed by the steady state Navier-Stokes equations is assumed. Recent advances in computational resources and algorithms have made it possible to solve the direct problem of determining such a flow through a body of known geometry. It is possible to obtain a set of adjoint equations, the solution to which permits the calculation of the direction and relative magnitude of change in the diffuser profile that leads to a higher pressure rise. The solution to the adjoint problem can be shown to represent an artificially constructed flow. This interpretation provides a means to construct numerical solutions to the adjoint equations that do not compromise the fully viscous nature of the problem. The algorithmic and computational aspects of solving the adjoint equations are addressed. The form of these set of equations is similar but not identical to the Navier-Stokes equations. In particular some issues related to boundary conditions and stability are discussed.
Adjoint-Operator Learning For A Neural Network
NASA Technical Reports Server (NTRS)
Barhen, Jacob; Toomarian, Nikzad
1993-01-01
Electronic neural networks made to synthesize initially unknown mathematical models of time-dependent phenomena or to learn temporally evolving patterns by use of algorithms based on adjoint operators. Algorithms less complicated, involve less computation and solve learning equations forward in time possibly simultaneously with equations of evolution of neural network, thereby both increasing computational efficiency and making real-time applications possible.
Adjoint Based Data Assimilation for an Ionospheric Model
NASA Astrophysics Data System (ADS)
Rosen, I. G.; Hajj, G. A.; Hajj, G. A.; Pi, X.; Pi, X.; Wang, C.; Wilson, B. D.
2001-05-01
The success of ionospheric modeling depends primarily on accurate knowledge of the forces (drivers) which enter into the collisional plasma hydrodynamic equations for the ionosphere and control the ionization as well as other dynamical and chemical processes. These include solar EUV and UV radiation, magnetospheric electric fields, particle precipitation, dynamo electric fields, thermospheric winds, neutral densities, and temperature. The determination of these model parameters from observational data is known as data assimilation. The data assimilation problem is formulated as a problem of minimizing a nonlinear functional, J (typically least squares) under a system of constraints consisting primarily of the underlying model equations. The performance index, J, can, in principle, be minimized using standard techniques such as the Newton's steepest decent method. There are however major technical challenges in practice. Since J is highly nonlinear and each evaluation of the functional requires the integration of the ionospheric model equations, computing the gradient vector of J with respect to the unknown parameters is time consuming. This problem is solved by use of the adjoint method. The ionospheric model used in this effort is for mid- and low-latitudes and consists of solving the continuity and momentum partial differential equations in four dimensional (three spatial dimensions and time) to compute the O+ density in the ionosphere and plasmasphere. We have developed codes for solving the forward model on a fixed grid. This makes it relatively straight forward to apply the adjoint method for computing gradients when doing nonlinear least squares based data assimilation. Because of the significant cost (in computational effort and CPU time) involved in performing a forward integration of the underlying 3-D model at any reasonable grid resolution, the use of the adjoint method for computing the gradients is indispensable. The adjoint method provides an elegant
A user's manual for MASH 1. 0: A Monte Carlo Adjoint Shielding Code System
Johnson, J.O.
1992-03-01
The Monte Carlo Adjoint Shielding Code System, MASH, calculates neutron and gamma-ray environments and radiation protection factors for armored military vehicles, structures, trenches, and other shielding configurations by coupling a forward discrete ordinates air-over-ground transport calculation with an adjoint Monte Carlo treatment of the shielding geometry. Efficiency and optimum use of computer time are emphasized. The code system include the GRTUNCL and DORT codes for air-over-ground transport calculations, the MORSE code with the GIFT5 combinatorial geometry package for adjoint shielding calculations, and several peripheral codes that perform the required data preparations, transformations, and coupling functions. MASH is the successor to the Vehicle Code System (VCS) initially developed at Oak Ridge National Laboratory (ORNL). The discrete ordinates calculation determines the fluence on a coupling surface surrounding the shielding geometry due to an external neutron/gamma-ray source. The Monte Carlo calculation determines the effectiveness of the fluence at that surface in causing a response in a detector within the shielding geometry, i.e., the dose importance'' of the coupling surface fluence. A coupling code folds the fluence together with the dose importance, giving the desired dose response. The coupling code can determine the dose response a a function of the shielding geometry orientation relative to the source, distance from the source, and energy response of the detector. This user's manual includes a short description of each code, the input required to execute the code along with some helpful input data notes, and a representative sample problem (input data and selected output edits) for each code.
NASA Astrophysics Data System (ADS)
Kovalskyy, V.; Henebry, G. M.
2011-05-01
Evapotranspiration (ET) flux constitutes a major component of both the water and energy balances at the land surface. Among the many factors that control evapotranspiration, phenology poses a major source of uncertainty in attempts to predict ET. Contemporary approaches to ET modeling and monitoring frequently summarize the complexity of the seasonal development of vegetation cover into static phenological trajectories (or climatologies) that lack sensitivity to changing environmental conditions. The Event Driven Phenology Model (EDPM) offers an alternative, interactive approach to representing phenology. This study presents the results of an experiment designed to illustrate the differences in ET arising from various techniques used to mimic phenology in models of land surface processes. The experiment compares and contrasts two realizations of static phenologies derived from long-term satellite observations of the Normalized Difference Vegetation Index (NDVI) against canopy trajectories produced by the interactive EDPM trained on flux tower observations. The assessment was carried out through validation of predicted ET against records collected by flux tower instruments. The VegET model (Senay, 2008) was used as a framework to estimate daily actual evapotranspiration and supplied with seasonal canopy trajectories produced by the EDPM and traditional techniques. The interactive approach presented the following advantages over phenology modeled with static climatologies: (a) lower prediction bias in crops; (b) smaller root mean square error in daily ET - 0.5 mm per day on average; (c) stable level of errors throughout the season similar among different land cover types and locations; and (d) better estimation of season duration and total seasonal ET.
NASA Astrophysics Data System (ADS)
Kovalskyy, V.; Henebry, G. M.
2012-01-01
Evapotranspiration (ET) flux constitutes a major component of both the water and energy balances at the land surface. Among the many factors that control evapotranspiration, phenology poses a major source of uncertainty in attempts to predict ET. Contemporary approaches to ET modeling and monitoring frequently summarize the complexity of the seasonal development of vegetation cover into static phenological trajectories (or climatologies) that lack sensitivity to changing environmental conditions. The Event Driven Phenology Model (EDPM) offers an alternative, interactive approach to representing phenology. This study presents the results of an experiment designed to illustrate the differences in ET arising from various techniques used to mimic phenology in models of land surface processes. The experiment compares and contrasts two realizations of static phenologies derived from long-term satellite observations of the Normalized Difference Vegetation Index (NDVI) against canopy trajectories produced by the interactive EDPM trained on flux tower observations. The assessment was carried out through validation of predicted ET against records collected by flux tower instruments. The VegET model (Senay, 2008) was used as a framework to estimate daily actual evapotranspiration and supplied with seasonal canopy trajectories produced by the EDPM and traditional techniques. The interactive approach presented the following advantages over phenology modeled with static climatologies: (a) lower prediction bias in crops; (b) smaller root mean square error in daily ET - 0.5 mm per day on average; (c) stable level of errors throughout the season similar among different land cover types and locations; and (d) better estimation of season duration and total seasonal ET.
NASA Astrophysics Data System (ADS)
Penenko, Vladimir; Tsvetova, Elena; Penenko, Alexey
2015-04-01
The proposed method is considered on an example of hydrothermodynamics and atmospheric chemistry models [1,2]. In the development of the existing methods for constructing numerical schemes possessing the properties of total approximation for operators of multiscale process models, we have developed a new variational technique, which uses the concept of adjoint integrating factors. The technique is as follows. First, a basic functional of the variational principle (the integral identity that unites the model equations, initial and boundary conditions) is transformed using Lagrange's identity and the second Green's formula. As a result, the action of the operators of main problem in the space of state functions is transferred to the adjoint operators defined in the space of sufficiently smooth adjoint functions. By the choice of adjoint functions the order of the derivatives becomes lower by one than those in the original equations. We obtain a set of new balance relationships that take into account the sources and boundary conditions. Next, we introduce the decomposition of the model domain into a set of finite volumes. For multi-dimensional non-stationary problems, this technique is applied in the framework of the variational principle and schemes of decomposition and splitting on the set of physical processes for each coordinate directions successively at each time step. For each direction within the finite volume, the analytical solutions of one-dimensional homogeneous adjoint equations are constructed. In this case, the solutions of adjoint equations serve as integrating factors. The results are the hybrid discrete-analytical schemes. They have the properties of stability, approximation and unconditional monotony for convection-diffusion operators. These schemes are discrete in time and analytic in the spatial variables. They are exact in case of piecewise-constant coefficients within the finite volume and along the coordinate lines of the grid area in each
Adjoint Formulation for an Embedded-Boundary Cartesian Method
NASA Technical Reports Server (NTRS)
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
Reconstruction of ocean circulation from sparse data using the adjoint method: LGM and the present
NASA Astrophysics Data System (ADS)
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
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
Examination of Observation Impacts derived from OSEs and Adjoint Models
NASA Technical Reports Server (NTRS)
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.
Three-Dimensional Turbulent RANS Adjoint-Based Error Correction
NASA Technical Reports Server (NTRS)
Park, Michael A.
2003-01-01
Engineering problems commonly require functional outputs of computational fluid dynamics (CFD) simulations with specified accuracy. These simulations are performed with limited computational resources. Computable error estimates offer the possibility of quantifying accuracy on a given mesh and predicting a fine grid functional on a coarser mesh. Such an estimate can be computed by solving the flow equations and the associated adjoint problem for the functional of interest. An adjoint-based error correction procedure is demonstrated for transonic inviscid and subsonic laminar and turbulent flow. A mesh adaptation procedure is formulated to target uncertainty in the corrected functional and terminate when error remaining in the calculation is less than a user-specified error tolerance. This adaptation scheme is shown to yield anisotropic meshes with corrected functionals that are more accurate for a given number of grid points then isotropic adapted and uniformly refined grids.
On improving storm surge forecasting using an adjoint optimal technique
NASA Astrophysics Data System (ADS)
Li, Yineng; Peng, Shiqiu; Yan, Jing; Xie, Lian
2013-12-01
A three-dimensional ocean model and its adjoint model are used to simultaneously optimize the initial conditions (IC) and the wind stress drag coefficient (Cd) for improving storm surge forecasting. To demonstrate the effect of this proposed method, a number of identical twin experiments (ITEs) with a prescription of different error sources and two real data assimilation experiments are performed. Results from both the idealized and real data assimilation experiments show that adjusting IC and Cd simultaneously can achieve much more improvements in storm surge forecasting than adjusting IC or Cd only. A diagnosis on the dynamical balance indicates that adjusting IC only may introduce unrealistic oscillations out of the assimilation window, which can be suppressed by the adjustment of the wind stress when simultaneously adjusting IC and Cd. Therefore, it is recommended to simultaneously adjust IC and Cd to improve storm surge forecasting using an adjoint technique.
Seismic Window Selection and Misfit Measurements for Global Adjoint Tomography
NASA Astrophysics Data System (ADS)
Lei, W.; Bozdag, E.; Lefebvre, M.; Podhorszki, N.; Smith, J. A.; Tromp, J.
2013-12-01
Global Adjoint Tomography requires fast parallel processing of large datasets. After obtaing the preprocessed observed and synthetic seismograms, we use the open source software packages FLEXWIN (Maggi et al. 2007) to select time windows and MEASURE_ADJ to make measurements. These measurements define adjoint sources for data assimilation. Previous versions of these tools work on a pair of SAC files---observed and synthetic seismic data for the same component and station, and loop over all seismic records associated with one earthquake. Given the large number of stations and earthquakes, the frequent read and write operations create severe I/O bottlenecks on modern computing platforms. We present new versions of these tools utilizing a new seismic data format, namely the Adaptive Seismic Data Format(ASDF). This new format shows superior scalability for applications on high-performance computers and accommodates various types of data, including earthquake, industry and seismic interferometry datasets. ASDF also provides user-friendly APIs, which can be easily integrated into the adjoint tomography workflow and combined with other data processing tools. In addition to solving the I/O bottleneck, we are making several improvements to these tools. For example, FLEXWIN is tuned to select windows for different types of earthquakes. To capture their distinct features, we categorize earthquakes by their depths and frequency bands. Moreover, instead of only picking phases between the first P arrival and the surface-wave arrivals, our aim is to select and assimilate many other later prominent phases in adjoint tomography. For example, in the body-wave band (17 s - 60 s), we include SKS, sSKS and their multiple, while in the surface-wave band (60 s - 120 s) we incorporate major-arc surface waves.
Forward and adjoint sensitivity computation of chaotic dynamical systems
Wang, Qiqi
2013-02-15
This paper describes a forward algorithm and an adjoint algorithm for computing sensitivity derivatives in chaotic dynamical systems, such as the Lorenz attractor. The algorithms compute the derivative of long time averaged “statistical” quantities to infinitesimal perturbations of the system parameters. The algorithms are demonstrated on the Lorenz attractor. We show that sensitivity derivatives of statistical quantities can be accurately estimated using a single, short trajectory (over a time interval of 20) on the Lorenz attractor.
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
Monopole condensation in two-flavor adjoint QCD
Cossu, Guido; D'Elia, Massimo; Di Giacomo, Adriano; Lacagnina, Giuseppe; Pica, Claudio
2008-04-01
In QCD with adjoint fermions, the deconfining transition takes place at a lower temperature than the chiral transition. We study the two transitions by use of the Polyakov loop, the monopole order parameter, and the chiral condensate. The deconfining transition is first order, the chiral is a crossover. The order parameters for confinement are not affected by the chiral transition. We conclude that the degrees of freedom relevant to confinement are different from those describing chiral symmetry.
Spectral monodromy of non-self-adjoint operators
NASA Astrophysics Data System (ADS)
Phan, Quang Sang
2014-01-01
In the present paper, we build a combinatorial invariant, called the "spectral monodromy" from the spectrum of a single (non-self-adjoint) h-pseudodifferential operator with two degrees of freedom in the semi-classical limit. Our inspiration comes from the quantum monodromy defined for the joint spectrum of an integrable system of n commuting self-adjoint h-pseudodifferential operators, given by S. Vu Ngoc ["Quantum monodromy in integrable systems," Commun. Math. Phys. 203(2), 465-479 (1999)]. The first simple case that we treat in this work is a normal operator. In this case, the discrete spectrum can be identified with the joint spectrum of an integrable quantum system. The second more complex case we propose is a small perturbation of a self-adjoint operator with a classical integrability property. We show that the discrete spectrum (in a small band around the real axis) also has a combinatorial monodromy. The main difficulty in this case is that we do not know the description of the spectrum everywhere, but only in a Cantor type set. In addition, we also show that the corresponding monodromy can be identified with the classical monodromy, defined by J. Duistermaat ["On global action-angle coordinates," Commun. Pure Appl. Math. 33(6), 687-706 (1980)].
Adjoint-based sensitivity analysis for reactor-safety applications
Parks, C.V.
1985-01-01
The application and usefulness of an adjoint-based methodology for performing sensitivity analysis on reactor safety computer codes is investigated. The adjoint-based methodology, referred to as differential sensitivity theory (DST), provides first-order derivatives of the calculated quantities of interest (responses) with respect to the input parameters. The basic theoretical development of DST is presented along with the needed general extensions for consideration of model discontinuities and a variety of useful response definitions. A simple analytic problem is used to highlight the general DST procedures. Finally, DST procedures presented in this work are applied to two highly nonlinear reactor accident analysis codes: (1) FASTGAS, a relatively small code for analysis of loss-of-decay-heat-removal accident in a gas-cooled fast reactor, and (2) an existing code called VENUS-II which is typically employed for analyzing the core disassembly phase of a hypothetical fast reactor accident. The two codes are different both in terms of complexity and in terms of the facets of DST which can be illustrated. Sensitivity results from the adjoint codes ADJGAS and VENUS-ADJ are verified with direct recalculations using perturbed input parameters. The effectiveness of the DST results for parameter ranking, prediction of response changes, and uncertainty analysis are illustrated. The conclusion drawn from this study is that DST is a viable, cost-effective methodology for accurate sensitivity analysis.
Self-adjoint time operators and invariant subspaces
NASA Astrophysics Data System (ADS)
Gómez, Fernando
2008-02-01
The question of existence of self-adjoint time operators for unitary evolutions in classical and quantum mechanics is revisited on the basis of Halmos-Helson theory of invariant subspaces, Sz.-Nagy-Foiaş dilation theory and Misra-Prigogine-Courbage theory of irreversibility. It is shown that the existence of self-adjoint time operators is equivalent to the intertwining property of the evolution plus the existence of simply invariant subspaces or rigid operator-valued functions for its Sz.-Nagy-Foiaş functional model. Similar equivalent conditions are given in terms of intrinsic randomness in the context of statistical mechanics. The rest of the contents are mainly a unifying review of the subject scattered throughout an unconnected literature. A well-known extensive set of equivalent conditions is derived from the above results; such conditions are written in terms of Schrrdinger couples, the Weyl commutation relation, incoming and outgoing subspaces, innovation processes, Lax-Phillips scattering, translation and spectral representations, and spectral properties. Also the natural procedure dealing with symmetric time operators in standard quantum mechanics involving their self-adjoint extensions is illustrated by considering the quantum Aharonov-Bohm time-of-arrival operator.
Unsteady Adjoint Approach for Design Optimization of Flapping Airfoils
NASA Technical Reports Server (NTRS)
Lee, Byung Joon; Liou, Meng-Sing
2012-01-01
This paper describes the work for optimizing the propulsive efficiency of flapping airfoils, i.e., improving the thrust under constraining aerodynamic work during the flapping flights by changing their shape and trajectory of motion with the unsteady discrete adjoint approach. For unsteady problems, it is essential to properly resolving time scales of motion under consideration and it must be compatible with the objective sought after. We include both the instantaneous and time-averaged (periodic) formulations in this study. For the design optimization with shape parameters or motion parameters, the time-averaged objective function is found to be more useful, while the instantaneous one is more suitable for flow control. The instantaneous objective function is operationally straightforward. On the other hand, the time-averaged objective function requires additional steps in the adjoint approach; the unsteady discrete adjoint equations for a periodic flow must be reformulated and the corresponding system of equations solved iteratively. We compare the design results from shape and trajectory optimizations and investigate the physical relevance of design variables to the flapping motion at on- and off-design conditions.
Spectral monodromy of non-self-adjoint operators
Phan, Quang Sang
2014-01-15
In the present paper, we build a combinatorial invariant, called the “spectral monodromy” from the spectrum of a single (non-self-adjoint) h-pseudodifferential operator with two degrees of freedom in the semi-classical limit. Our inspiration comes from the quantum monodromy defined for the joint spectrum of an integrable system of n commuting self-adjoint h-pseudodifferential operators, given by S. Vu Ngoc [“Quantum monodromy in integrable systems,” Commun. Math. Phys. 203(2), 465–479 (1999)]. The first simple case that we treat in this work is a normal operator. In this case, the discrete spectrum can be identified with the joint spectrum of an integrable quantum system. The second more complex case we propose is a small perturbation of a self-adjoint operator with a classical integrability property. We show that the discrete spectrum (in a small band around the real axis) also has a combinatorial monodromy. The main difficulty in this case is that we do not know the description of the spectrum everywhere, but only in a Cantor type set. In addition, we also show that the corresponding monodromy can be identified with the classical monodromy, defined by J. Duistermaat [“On global action-angle coordinates,” Commun. Pure Appl. Math. 33(6), 687–706 (1980)].
Improving the Fit of a Land-Surface Model to Data Using its Adjoint
NASA Astrophysics Data System (ADS)
Raoult, N.; Jupp, T. E.; Cox, P. M.; Luke, C.
2015-12-01
Land-surface models (LSMs) are of growing importance in the world of climate prediction. They are crucial components of larger Earth system models that are aimed at understanding the effects of land surface processes on the global carbon cycle. The Joint UK Land Environment Simulator (JULES) is the land-surface model used by the UK Met Office. It has been 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. adJULES presents an opportunity to confront JULES with many different observations, and make improvements to the model parameterisation. In the newest version of adJULES, multiple sites can be used in the calibration, to giving a generic set of parameters that can be generalised over plant functional types. We present an introduction to the adJULES system and its applications to data from a variety of flux tower sites. 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.
Solid oxide fuel cell simulation and design optimization with numerical adjoint techniques
NASA Astrophysics Data System (ADS)
Elliott, Louie C.
This dissertation reports on the application of numerical optimization techniques as applied to fuel cell simulation and design. Due to the "multi-physics" inherent in a fuel cell, which results in a highly coupled and non-linear behavior, an experimental program to analyze and improve the performance of fuel cells is extremely difficult. This program applies new optimization techniques with computational methods from the field of aerospace engineering to the fuel cell design problem. After an overview of fuel cell history, importance, and classification, a mathematical model of solid oxide fuel cells (SOFC) is presented. The governing equations are discretized and solved with computational fluid dynamics (CFD) techniques including unstructured meshes, non-linear solution methods, numerical derivatives with complex variables, and sensitivity analysis with adjoint methods. Following the validation of the fuel cell model in 2-D and 3-D, the results of the sensitivity analysis are presented. The sensitivity derivative for a cost function with respect to a design variable is found with three increasingly sophisticated techniques: finite difference, direct differentiation, and adjoint. A design cycle is performed using a simple optimization method to improve the value of the implemented cost function. The results from this program could improve fuel cell performance and lessen the world's dependence on fossil fuels.
Paulot, Fabien; Jacob, Daniel J; Henze, Daven K
2013-04-01
Anthropogenic enrichment of reactive nitrogen (Nr) deposition is an ecological concern. We use the adjoint of a global 3-D chemical transport model (GEOS-Chem) to identify the sources and processes that control Nr deposition to an ensemble of biodiversity hotspots worldwide and two U.S. national parks (Cuyahoga and Rocky Mountain). We find that anthropogenic sources dominate deposition at all continental sites and are mainly regional (less than 1000 km) in origin. In Hawaii, Nr supply is controlled by oceanic emissions of ammonia (50%) and anthropogenic sources (50%), with important contributions from Asia and North America. Nr deposition is also sensitive in complicated ways to emissions of SO2, which affect Nr gas-aerosol partitioning, and of volatile organic compounds (VOCs), which affect oxidant concentrations and produce organic nitrate reservoirs. For example, VOC emissions generally inhibit deposition of locally emitted NOx but significantly increase Nr deposition downwind. However, in polluted boreal regions, anthropogenic VOC emissions can promote Nr deposition in winter. Uncertainties in chemical rate constants for OH + NO2 and NO2 hydrolysis also complicate the determination of source-receptor relationships for polluted sites in winter. Application of our adjoint sensitivities to the representative concentration pathways (RCPs) scenarios for 2010-2050 indicates that future decreases in Nr deposition due to NOx emission controls will be offset by concurrent increases in ammonia emissions from agriculture. PMID:23458244
NASA Astrophysics Data System (ADS)
Yu, Jia; Ji, Lucheng; Li, Weiwei; Yi, Weilin
2016-06-01
Adjoint method is an important tool for design refinement of multistage compressors. However, the radial static pressure distribution deviates during the optimization procedure and deteriorates the overall performance, producing final designs that are not well suited for realistic engineering applications. In previous development work on multistage turbomachinery blade optimization using adjoint method and thin shear-layer N-S equations, the entropy production is selected as the objective function with given mass flow rate and total pressure ratio as imposed constraints. The radial static pressure distribution at the interfaces between rows is introduced as a new constraint in the present paper. The approach is applied to the redesign of a five-stage axial compressor, and the results obtained with and without the constraint on the radial static pressure distribution at the interfaces between rows are discussed in detail. The results show that the redesign without the radial static pressure distribution constraint (RSPDC) gives an optimal solution that shows deviations on radial static pressure distribution, especially at rotor exit tip region. On the other hand, the redesign with the RSPDC successfully keeps the radial static pressure distribution at the interfaces between rows and make sure that the optimization results are applicable in a practical engineering design.
NASA Astrophysics Data System (ADS)
Zheng, J. P.; Lee, C.-T. A.; Lu, J. G.; Zhao, J. H.; Wu, Y. B.; Xia, B.; Li, X. Y.; Zhang, J. F.; Liu, Y. S.
2015-01-01
Continents are underlain by thick, cold thermal boundary layers. Thermal contraction should render these boundary layers negatively buoyant and unstable; this is why old, cold oceanic lithospheres subduct. However, the ancient lithospheric roots of many continents appear to have existed for billions of years. In the common view, this preservation is due to the fact that the thermal boundary layers are compositionally distinct from the ambient mantle in that they are highly melt-depleted and dehydrated; the former provides positive buoyancy and the latter provides strength. Here, we show using mantle xenoliths that the Precambrian South China Block originally was underlain by highly depleted mantle, but has been refertilized via silicate melts generated from the asthenosphere. It is now more fertile than the ambient convecting mantle and is intrinsically denser by more than 1.5%. Achieving sufficient melt generation for refertilization is only possible if the lithosphere is thin enough to provide "headspace" for decompression melting. Thus, continental boundary layers thinner than the maximum depth of melting should experience refertilization, whereas thicker continents would altogether suppress melting and hence the potential for refertilization. We propose that refertilization, once initiated, will destabilize the base of the continent; this in turn will increase the amount of "headspace" and promote further refertilization, resulting in a positive feedback that could culminate in lithospheric destruction. By contrast, continents that are thick enough may not experience significant refertilization. This suggests that initial lithospheric thickness, as well as lithospheric composition, may be important for defining the fate of continents.
A new mathematical adjoint for the modified SAAF_{-SN} equations
Schunert, Sebastian; Wang, Yaqi; Martineau, Richard; DeHart, Mark D.
2015-01-01
We present a new adjoint FEM weak form, which can be directly used for evaluating the mathematical adjoint, suitable for perturbation calculations, of the self-adjoint angular flux SN equations (SAAF_{-SN}) without construction and transposition of the underlying coefficient matrix. Stabilization schemes incorporated in the described SAAF_{-SN} method make the mathematical adjoint distinct from the physical adjoint, i.e. the solution of the continuous adjoint equation with SAAF_{-SN} . This weak form is implemented into RattleSnake, the MOOSE (Multiphysics Object-Oriented Simulation Environment) based transport solver. Numerical results verify the correctness of the implementation and show its utility both for fixed source and eigenvalue problems.
Nonlinear self-adjointness and conservation laws of Klein-Gordon-Fock equation with central symmetry
NASA Astrophysics Data System (ADS)
Abdulwahhab, Muhammad Alim
2015-05-01
The concept of nonlinear self-adjointness, introduced by Ibragimov, has significantly extends approaches to constructing conservation laws associated with symmetries since it incorporates the strict self-adjointness, the quasi self-adjointness as well as the usual linear self-adjointness. Using this concept, the nonlinear self-adjointness condition for the Klein-Gordon-Fock equation was established and subsequently used to construct simplified but infinitely many nontrivial and independent conserved vectors. The Noether's theorem was further applied to the Klein-Gordon-Fock equation to explore more distinct first integrals, result shows that conservation laws constructed through this approach are exactly the same as those obtained under strict self-adjointness of Ibragimov's method.
Adjoint-Based Sensitivity Maps for the Nearshore
NASA Astrophysics Data System (ADS)
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
NASA Astrophysics Data System (ADS)
Rock, B. N.; Hale, S. R.; Graham, K. J.; Hayden, L.; Barber, L.; Perry, C.; Schloss, J.; Sullivan, E.; Yuan, J.; Abebe, E.; Mitchell, L.; Abrams, E.; Gagnon, M.
2008-12-01
Watershed Watch (NSF 0525433) engages early undergraduate students from two-year and four-year colleges in student-driven full inquiry-based instruction in the biogeosciences. Program goals for Watershed Watch are to test if inquiry-rich student-driven projects sufficiently engage undeclared students (or noncommittal STEM majors) to declare a STEM major (or remain with their STEM major). A significant component of this program is an intensive two-week Summer course, in which undeclared freshmen research various aspects of a local watershed. Students develop their own research questions and study design, collect and analyze data, and produce a scientific or an oral poster presentation. The course objectives, curriculum and schedule are presented as a model for dissemination for other institutions and programs seeking to develop inquiry-rich courses designed to attract students into biogeoscience disciplines. Data from self-reported student feedback indicated the most important factors explaining high-levels of student motivation and research excellence in the course are 1) working with committed, energetic, and enthusiastic faculty mentors; and 2) faculty mentors demonstrating high degrees of teamwork and coordination.
Adjoint-based sensitivity analysis for reactor safety applications
Parks, C.V.
1986-08-01
The application and usefulness of an adjoint-based methodology for performing sensitivity analysis on reactor safety computer codes is investigated. The adjoint-based methodology, referred to as differential sensitivity theory (DST), provides first-order derivatives of the calculated quantities of interest (responses) with respect to the input parameters. The basic theoretical development of DST is presented along with the needed general extensions for consideration of model discontinuities and a variety of useful response definitions. A simple analytic problem is used to highlight the general DST procedures. finally, DST procedures presented in this work are applied to two highly nonlinear reactor accident analysis codes: (1) FASTGAS, a relatively small code for analysis of a loss-of-decay-heat-removal accident in a gas-cooled fast reactor, and (2) an existing code called VENUS-II which has been employed for analyzing the core disassembly phase of a hypothetical fast reactor accident. The two codes are different both in terms of complexity and in terms of the facets of DST which can be illustrated. Sensitivity results from the adjoint codes ADJGAS and VENUS-ADJ are verified with direct recalcualtions using perturbed input parameters. The effectiveness of the DST results for parameter ranking, prediction of response changes, and uncertainty analysis are illustrated. The conclusion drawn from this study is that DST is a viable, cost-effective methodology for accurate sensitivity analysis. In addition, a useful sensitivity tool for use in the fast reactor safety area has been developed in VENUS-ADJ. Future work needs to concentrate on combining the accurate first-order derivatives/results from DST with existing methods (based solely on direct recalculations) for higher-order response surfaces.
Adjoint-based optimization for understanding and suppressing jet noise
NASA Astrophysics Data System (ADS)
Freund, Jonathan B.
2011-08-01
Advanced simulation tools, particularly large-eddy simulation techniques, are becoming capable of making quality predictions of jet noise for realistic nozzle geometries and at engineering relevant flow conditions. Increasing computer resources will be a key factor in improving these predictions still further. Quality prediction, however, is only a necessary condition for the use of such simulations in design optimization. Predictions do not themselves lead to quieter designs. They must be interpreted or harnessed in some way that leads to design improvements. As yet, such simulations have not yielded any simplifying principals that offer general design guidance. The turbulence mechanisms leading to jet noise remain poorly described in their complexity. In this light, we have implemented and demonstrated an aeroacoustic adjoint-based optimization technique that automatically calculates gradients that point the direction in which to adjust controls in order to improve designs. This is done with only a single flow solutions and a solution of an adjoint system, which is solved at computational cost comparable to that for the flow. Optimization requires iterations, but having the gradient information provided via the adjoint accelerates convergence in a manner that is insensitive to the number of parameters to be optimized. This paper, which follows from a presentation at the 2010 IUTAM Symposium on Computational Aero-Acoustics for Aircraft Noise Prediction, reviews recent and ongoing efforts by the author and co-workers. It provides a new formulation of the basic approach and demonstrates the approach on a series of model flows, culminating with a preliminary result for a turbulent jet.
Adjoint Techniques for Topology Optimization of Structures Under Damage Conditions
NASA Technical Reports Server (NTRS)
Akgun, Mehmet A.; Haftka, Raphael T.
2000-01-01
The objective of this cooperative agreement was to seek computationally efficient ways to optimize aerospace structures subject to damage tolerance criteria. Optimization was to involve sizing as well as topology optimization. The work was done in collaboration with Steve Scotti, Chauncey Wu and Joanne Walsh at the NASA Langley Research Center. Computation of constraint sensitivity is normally the most time-consuming step of an optimization procedure. The cooperative work first focused on this issue and implemented the adjoint method of sensitivity computation (Haftka and Gurdal, 1992) in an optimization code (runstream) written in Engineering Analysis Language (EAL). The method was implemented both for bar and plate elements including buckling sensitivity for the latter. Lumping of constraints was investigated as a means to reduce the computational cost. Adjoint sensitivity computation was developed and implemented for lumped stress and buckling constraints. Cost of the direct method and the adjoint method was compared for various structures with and without lumping. The results were reported in two papers (Akgun et al., 1998a and 1999). It is desirable to optimize topology of an aerospace structure subject to a large number of damage scenarios so that a damage tolerant structure is obtained. Including damage scenarios in the design procedure is critical in order to avoid large mass penalties at later stages (Haftka et al., 1983). A common method for topology optimization is that of compliance minimization (Bendsoe, 1995) which has not been used for damage tolerant design. In the present work, topology optimization is treated as a conventional problem aiming to minimize the weight subject to stress constraints. Multiple damage configurations (scenarios) are considered. Each configuration has its own structural stiffness matrix and, normally, requires factoring of the matrix and solution of the system of equations. Damage that is expected to be tolerated is local
A self-adjoint decomposition of the radial momentum operator
NASA Astrophysics Data System (ADS)
Liu, Q. H.; Xiao, S. F.
2015-12-01
With acceptance of the Dirac's observation that the canonical quantization entails using Cartesian coordinates, we examine the operator erPr rather than Pr itself and demonstrate that there is a decomposition of erPr into a difference of two self-adjoint but noncommutative operators, in which one is the total momentum and another is the transverse one. This study renders the operator Pr indirectly measurable and physically meaningful, offering an explanation of why the mean value of Pr over a quantum mechanical state makes sense and supporting Dirac's claim that Pr "is real and is the true momentum conjugate to r".
Examining Tropical Cyclone - Kelvin Wave Interactions using Adjoint Diagnostics
NASA Astrophysics Data System (ADS)
Reynolds, C. A.; Doyle, J. D.; Hong, X.
2015-12-01
Adjoint-based tools can provide valuable insight into the mechanisms that influence the evolution and predictability of atmospheric phenomena, as they allow for the efficient and rigorous computation of forecast sensitivity to changes in the initial state. We apply adjoint-based tools from the non-hydrostatic Coupled Atmosphere/Ocean Mesoscale Prediction System (COAMPS) to explore the initial-state sensitivity and interactions between a tropical cyclone and atmospheric equatorial waves associated with the Madden Julian Oscillation (MJO) in the Indian Ocean during the DYNAMO field campaign. The development of Tropical Cyclone 5 (TC05) coincided with the passage of an equatorial Kelvin wave and westerly wind burst associated with an MJO that developed in the Indian Ocean in late November 2011, but it was unclear if and how one affected the other. COAMPS 24-h and 36-h adjoint sensitivities are analyzed for both TC05 and the equatorial waves to understand how the evolution of each system is sensitive to the other. The sensitivity of equatorial westerlies in the western Indian Ocean on 23 November shares characteristics with the classic Gill (1980) Rossby and Kelvin wave response to symmetric heating about the equator, including symmetric cyclonic circulations to the north and south of the westerlies, and enhanced heating in the area of convergence between the equatorial westerlies and easterlies. In addition, there is sensitivity in the Bay of Bengal associated with the cyclonic circulation that eventually develops into TC05. At the same time, the developing TC05 system shows strongest sensitivity to local wind and heating perturbations, but sensitivity to the equatorial westerlies is also clear. On 24 November, when the Kelvin wave is immediately south of the developing tropical cyclone, both phenomena are sensitive to each other. On 25 November TC05 no longer shows sensitivity to the Kelvin wave, while the Kelvin Wave still exhibits some weak sensitivity to TC05. In
NASA Astrophysics Data System (ADS)
Xie, Zhinan; Komatitsch, Dimitri; Martin, Roland; Matzen, René
2014-09-01
In recent years, the application of time-domain adjoint methods to improve large, complex underground tomographic models at the regional scale has led to new challenges for the numerical simulation of forward or adjoint elastic wave propagation problems. An important challenge is to design an efficient infinite-domain truncation method suitable for accurately truncating an infinite domain governed by the second-order elastic wave equation written in displacement and computed based on a finite-element (FE) method. In this paper, we make several steps towards this goal. First, we make the 2-D convolution formulation of the complex-frequency-shifted unsplit-field perfectly matched layer (CFS-UPML) derived in previous work more flexible by providing a new treatment to analytically remove singular parameters in the formulation. We also extend this new formulation to 3-D. Furthermore, we derive the auxiliary differential equation (ADE) form of CFS-UPML, which allows for extension to higher order time schemes and is easier to implement. Secondly, we rigorously derive the CFS-UPML formulation for time-domain adjoint elastic wave problems, which to our knowledge has never been done before. Thirdly, in the case of classical low-order FE methods, we show numerically that we achieve long-time stability for both forward and adjoint problems both for the convolution and the ADE formulations. In the case of higher order Legendre spectral-element methods, we show that weak numerical instabilities can appear in both formulations, in particular if very small mesh elements are present inside the absorbing layer, but we explain how these instabilities can be delayed as much as needed by using a stretching factor to reach numerical stability in practice for applications. Fourthly, in the case of adjoint problems with perfectly matched absorbing layers we introduce a computationally efficient boundary storage strategy by saving information along the interface between the CFS-UPML and
Comparison of the Monte Carlo adjoint-weighted and differential operator perturbation methods
Kiedrowski, Brian C; Brown, Forrest B
2010-01-01
Two perturbation theory methodologies are implemented for k-eigenvalue calculations in the continuous-energy Monte Carlo code, MCNP6. A comparison of the accuracy of these techniques, the differential operator and adjoint-weighted methods, is performed numerically and analytically. Typically, the adjoint-weighted method shows better performance over a larger range; however, there are exceptions.
A User's Manual for MASH V1.5 - A Monte Carlo Adjoint Shielding Code System
C. O. Slater; J. M. Barnes; J. O. Johnson; J.D. Drischler
1998-10-01
The Monte Carlo ~djoint ~ielding Code System, MASH, calculates neutron and gamma- ray environments and radiation protection factors for armored military vehicles, structures, trenches, and other shielding configurations by coupling a forward discrete ordinates air- over-ground transport calculation with an adjoint Monte Carlo treatment of the shielding geometry. Efficiency and optimum use of computer time are emphasized. The code system includes the GRTUNCL and DORT codes for air-over-ground transport calculations, the MORSE code with the GIFT5 combinatorial geometry package for adjoint shielding calculations, and several peripheral codes that perform the required data preparations, transformations, and coupling functions. The current version, MASH v 1.5, is the successor to the original MASH v 1.0 code system initially developed at Oak Ridge National Laboratory (ORNL). The discrete ordinates calculation determines the fluence on a coupling surface surrounding the shielding geometry due to an external neutron/gamma-ray source. The Monte Carlo calculation determines the effectiveness of the fluence at that surface in causing a response in a detector within the shielding geometry, i.e., the "dose importance" of the coupling surface fluence. A coupling code folds the fluence together with the dose importance, giving the desired dose response. The coupling code can determine the dose response as a function of the shielding geometry orientation relative to the source, distance from the source, and energy response of the detector. This user's manual includes a short description of each code, the input required to execute the code along with some helpful input data notes, and a representative sample problem.
Kim, Min-Geun; Jang, Hong-Lae; Cho, Seonho
2013-05-01
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.
Probability density adjoint for sensitivity analysis of the Mean of Chaos
Blonigan, Patrick J. Wang, Qiqi
2014-08-01
Sensitivity analysis, especially adjoint based sensitivity analysis, is a powerful tool for engineering design which allows for the efficient computation of sensitivities with respect to many parameters. However, these methods break down when used to compute sensitivities of long-time averaged quantities in chaotic dynamical systems. This paper presents a new method for sensitivity analysis of ergodic chaotic dynamical systems, the density adjoint method. The method involves solving the governing equations for the system's invariant measure and its adjoint on the system's attractor manifold rather than in phase-space. This new approach is derived for and demonstrated on one-dimensional chaotic maps and the three-dimensional Lorenz system. It is found that the density adjoint computes very finely detailed adjoint distributions and accurate sensitivities, but suffers from large computational costs.
Limitations of Adjoint-Based Optimization for Separated Flows
NASA Astrophysics Data System (ADS)
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.
Adelman, James S; Carter, Amanda W; Hopkins, William A; Hawley, Dana M
2013-10-23
Although ambient temperature has diverse effects on disease dynamics, few studies have examined how temperature alters pathogen transmission by changing host physiology or behaviour. Here, we test whether reducing ambient temperature alters host foraging, pathology and the potential for fomite transmission of the bacterial pathogen Mycoplasma gallisepticum (MG), which causes seasonal outbreaks of severe conjunctivitis in house finches (Haemorhous mexicanus). We housed finches at temperatures within or below the thermoneutral zone to manipulate food intake by altering energetic requirements of thermoregulation. We predicted that pathogen deposition on bird feeders would increase with temperature-driven increases in food intake and with conjunctival pathology. As expected, housing birds below the thermoneutral zone increased food consumption. Despite this difference, pathogen deposition on feeders did not vary across temperature treatments. However, pathogen deposition increased with conjunctival pathology, independently of temperature and pathogen load, suggesting that MG could enhance its transmission by increasing virulence. Our results suggest that in this system, host physiological responses are more important for transmission potential than temperature-dependent alterations in feeding. Understanding such behavioural and physiological contributions to disease transmission is critical to linking individual responses to climate with population-level disease dynamics. PMID:23966599
A weighted adjoint-source for weight-window generation by means of a linear tally combination
Sood, Avneet; Booth, Thomas E; Solomon, Clell J
2009-01-01
A new importance estimation technique has been developed that allows weight-window optimization for a linear combination of tallies. This technique has been implemented in a local version of MCNP and effectively weights the adjoint source term for each tally in the combination. Optimizing weight window parameters for the linear tally combination allows the user to optimize weight windows for multiple regions at once. In this work, we present our results of solutions to an analytic three-tally-region test problem and a flux calculation on a 100,000 voxel oil-well logging tool problem.
Seismic Imaging of VTI, HTI and TTI based on Adjoint Methods
NASA Astrophysics Data System (ADS)
Rusmanugroho, H.; Tromp, J.
2014-12-01
Recent studies show that isotropic seismic imaging based on adjoint method reduces low-frequency artifact caused by diving waves, which commonly occur in two-wave wave-equation migration, such as Reverse Time Migration (RTM). Here, we derive new expressions of sensitivity kernels for Vertical Transverse Isotropy (VTI) using the Thomsen parameters (ɛ, δ, γ) plus the P-, and S-wave speeds (α, β) as well as via the Chen & Tromp (GJI 2005) parameters (A, C, N, L, F). For Horizontal Transverse Isotropy (HTI), these parameters depend on an azimuthal angle φ, where the tilt angle θ is equivalent to 90°, and for Tilted Transverse Isotropy (TTI), these parameters depend on both the azimuth and tilt angles. We calculate sensitivity kernels for each of these two approaches. Individual kernels ("images") are numerically constructed based on the interaction between the regular and adjoint wavefields in smoothed models which are in practice estimated through Full-Waveform Inversion (FWI). The final image is obtained as a result of summing all shots, which are well distributed to sample the target model properly. The impedance kernel, which is a sum of sensitivity kernels of density and the Thomsen or Chen & Tromp parameters, looks crisp and promising for seismic imaging. The other kernels suffer from low-frequency artifacts, similar to traditional seismic imaging conditions. However, all sensitivity kernels are important for estimating the gradient of the misfit function, which, in combination with a standard gradient-based inversion algorithm, is used to minimize the objective function in FWI.
NASA Astrophysics Data System (ADS)
Kubina, Filip; Moczo, Peter; Kristek, Jozef; Michlik, Filip
2016-04-01
Adjoint tomography has proven an irreplaceable useful tool in exploring Earth's structure in the regional and global scales. It has not been widely applied for improving models of local surface sedimentary structures (LSSS) in numerical predictions of earthquake ground motion (EGM). Anomalous earthquake motions and corresponding damage in earthquakes are often due to site effects in local surface sedimentary basins. Because majority of world population is located atop surface sedimentary basins, it is important to predict EGM at these sites during future earthquakes. A major lesson learned from dedicated international tests focused on numerical prediction of EGM in LSSS is that it is hard to reach better agreement between data and synthetics without an improved structural model. If earthquake records are available for sites atop a LSSS it is natural to consider them for improving the structural model. Computationally efficient adjoint tomography might be a proper tool. A seismic wavefield in LSSS is relatively very complex due to diffractions, conversions, interference and often also resonant phenomena. In shallow basins, the first arrivals are not suitable for inversion due to almost vertical incidence and thus insufficient vertical resolution. Later wavefield consists mostly of local surface waves often without separated wave groups. Consequently, computed kernels are complicated and not suitable for inversion without pre-processing. The spatial complexity of a kernel can be dramatic in a typical situation with relatively low number of sources (local earthquakes) and surface receivers. This complexity can be simplified by directionally-dependent smoothing and spatially-dependent normalization that condition reasonable convergence. A multiscale approach seems necessary given the usual difference between the available and true models. Interestingly, only a successive inversion of μ and λ elastic moduli, and different scale sequences lead to good results.
Dixit, Prakash N; Telleria, Roberto
2015-04-01
Inter-annual and seasonal variability in climatic parameters, most importantly rainfall, have potential to cause climate-induced risk in long-term crop production. Short-term field studies do not capture the full nature of such risk and the extent to which modifications to crop, soil and water management recommendations may be made to mitigate the extent of such risk. Crop modeling studies driven by long-term daily weather data can predict the impact of climate-induced risk on crop growth and yield however, the availability of long-term daily weather data can present serious constraints to the use of crop models. To tackle this constraint, two weather generators namely, LARS-WG and MarkSim, were evaluated in order to assess their capabilities of reproducing frequency distributions, means, variances, dry spell and wet chains of observed daily precipitation, maximum and minimum temperature, and solar radiation for the eight locations across cropping areas of Northern Syria and Lebanon. Further, the application of generated long-term daily weather data, with both weather generators, in simulating barley growth and yield was also evaluated. We found that overall LARS-WG performed better than MarkSim in generating daily weather parameters and in 50 years continuous simulation of barley growth and yield. Our findings suggest that LARS-WG does not necessarily require long-term e.g., >30 years observed weather data for calibration as generated results proved to be satisfactory with >10 years of observed data except in area with higher altitude. Evaluating these weather generators and the ability of generated weather data to perform long-term simulation of crop growth and yield is an important first step to assess the impact of future climate on yields, and to identify promising technologies to make agricultural systems more resilient in the given region. PMID:25590537
A user`s manual for MASH 1.0: A Monte Carlo Adjoint Shielding Code System
Johnson, J.O.
1992-03-01
The Monte Carlo Adjoint Shielding Code System, MASH, calculates neutron and gamma-ray environments and radiation protection factors for armored military vehicles, structures, trenches, and other shielding configurations by coupling a forward discrete ordinates air-over-ground transport calculation with an adjoint Monte Carlo treatment of the shielding geometry. Efficiency and optimum use of computer time are emphasized. The code system include the GRTUNCL and DORT codes for air-over-ground transport calculations, the MORSE code with the GIFT5 combinatorial geometry package for adjoint shielding calculations, and several peripheral codes that perform the required data preparations, transformations, and coupling functions. MASH is the successor to the Vehicle Code System (VCS) initially developed at Oak Ridge National Laboratory (ORNL). The discrete ordinates calculation determines the fluence on a coupling surface surrounding the shielding geometry due to an external neutron/gamma-ray source. The Monte Carlo calculation determines the effectiveness of the fluence at that surface in causing a response in a detector within the shielding geometry, i.e., the ``dose importance`` of the coupling surface fluence. A coupling code folds the fluence together with the dose importance, giving the desired dose response. The coupling code can determine the dose response a a function of the shielding geometry orientation relative to the source, distance from the source, and energy response of the detector. This user`s manual includes a short description of each code, the input required to execute the code along with some helpful input data notes, and a representative sample problem (input data and selected output edits) for each code.
Mass anomalous dimension in SU(2) with two adjoint fermions
Bursa, Francis; Del Debbio, Luigi; Keegan, Liam; Pica, Claudio; Pickup, Thomas
2010-01-01
We study SU(2) lattice gauge theory with two flavors of Dirac fermions in the adjoint representation. We measure the running of the coupling in the Schroedinger functional scheme and find it is consistent with existing results. We discuss how systematic errors affect the evidence for an infrared fixed point (IRFP). We present the first measurement of the running of the mass in the Schroedinger functional scheme. The anomalous dimension of the chiral condensate, which is relevant for phenomenological applications, can be easily extracted from the running of the mass, under the assumption that the theory has an IRFP. At the current level of accuracy, we can estimate 0.05<{gamma}<0.56 at the IRFP.
Infrared regime of SU(2) with one adjoint Dirac flavor
NASA Astrophysics Data System (ADS)
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.
Optimizing spectral wave estimates with adjoint-based sensitivity maps
NASA Astrophysics Data System (ADS)
Orzech, Mark; Veeramony, Jay; Flampouris, Stylianos
2014-04-01
A discrete numerical adjoint has recently been developed for the stochastic wave model SWAN. In the present study, this adjoint code is used to construct spectral sensitivity maps for two nearshore domains. The maps display the correlations of spectral energy levels throughout the domain with the observed energy levels at a selected location or region of interest (LOI/ROI), providing a full spectrum of values at all locations in the domain. We investigate the effectiveness of sensitivity maps based on significant wave height ( H s ) in determining alternate offshore instrument deployment sites when a chosen nearshore location or region is inaccessible. Wave and bathymetry datasets are employed from one shallower, small-scale domain (Duck, NC) and one deeper, larger-scale domain (San Diego, CA). The effects of seasonal changes in wave climate, errors in bathymetry, and multiple assimilation points on sensitivity map shapes and model performance are investigated. Model accuracy is evaluated by comparing spectral statistics as well as with an RMS skill score, which estimates a mean model-data error across all spectral bins. Results indicate that data assimilation from identified high-sensitivity alternate locations consistently improves model performance at nearshore LOIs, while assimilation from low-sensitivity locations results in lesser or no improvement. Use of sub-sampled or alongshore-averaged bathymetry has a domain-specific effect on model performance when assimilating from a high-sensitivity alternate location. When multiple alternate assimilation locations are used from areas of lower sensitivity, model performance may be worse than with a single, high-sensitivity assimilation point.
NASA Technical Reports Server (NTRS)
Arian, Eyal; Salas, Manuel D.
1997-01-01
We derive the adjoint equations for problems in aerodynamic optimization which are improperly considered as "inadmissible." For example, a cost functional which depends on the density, rather than on the pressure, is considered "inadmissible" for an optimization problem governed by the Euler equations. We show that for such problems additional terms should be included in the Lagrangian functional when deriving the adjoint equations. These terms are obtained from the restriction of the interior PDE to the control surface. Demonstrations of the explicit derivation of the adjoint equations for "inadmissible" cost functionals are given for the potential, Euler, and Navier-Stokes equations.
Use of adjoint methods in the probabilistic finite element approach to fracture mechanics
NASA Technical Reports Server (NTRS)
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.
NASA Technical Reports Server (NTRS)
Nguyen, Nhan T.; Hornby, Gregory; Ishihara, Abe
2013-01-01
This paper describes two methods of trajectory optimization to obtain an optimal trajectory of minimum-fuel- to-climb for an aircraft. The first method is based on the adjoint method, and the second method is based on a direct trajectory optimization method using a Chebyshev polynomial approximation and cubic spine approximation. The approximate optimal trajectory will be compared with the adjoint-based optimal trajectory which is considered as the true optimal solution of the trajectory optimization problem. The adjoint-based optimization problem leads to a singular optimal control solution which results in a bang-singular-bang optimal control.
Time domain adjoint sensitivity analysis of electromagnetic problems with nonlinear media.
Bakr, Mohamed H; Ahmed, Osman S; El Sherif, Mohamed H; Nomura, Tsuyoshi
2014-05-01
In this paper, we propose a theory for wideband adjoint sensitivity analysis of problems with nonlinear media. We show that the sensitivities of the desired response with respect to all shape and material parameters are obtained through one extra adjoint simulation. Unlike linear problems, the system matrices of this adjoint simulation are time varying. Their values are determined during the original simulation. The proposed theory exploits the time-domain transmission line modeling (TLM) and provides an efficient AVM approach for sensitivity analysis of general time domain objective functions. The theory has been illustrated through a number of examples. PMID:24921783
Adjoint-Based Methods for Estimating CO2 Sources and Sinks from Atmospheric Concentration Data
NASA Technical Reports Server (NTRS)
Andrews, Arlyn E.
2003-01-01
Work to develop adjoint-based methods for estimating CO2 sources and sinks from atmospheric concentration data was initiated in preparation for last year's summer institute on Carbon Data Assimilation (CDAS) at the National Center for Atmospheric Research in Boulder, CO. The workshop exercises used the GSFC Parameterized Chemistry and Transport Model and its adjoint. Since the workshop, a number of simulations have been run to evaluate the performance of the model adjoint. Results from these simulations will be presented, along with an outline of challenges associated with incorporating a variety of disparate data sources, from sparse, but highly precise, surface in situ observations to less accurate, global future satellite observations.
On basic conditions to generate multi-adjoint concept lattices via Galois connections
NASA Astrophysics Data System (ADS)
Díaz-Moreno, J. C.; Medina, J.; Ojeda-Aciego, M.
2014-02-01
This paper introduces sufficient and necessary conditions with respect to the fuzzy operators considered in a multi-adjoint frame under which the standard combinations of multi-adjoint sufficiency, possibility, and necessity operators form (antitone or isotone) Galois connections. The underlying idea is to study the minimal algebraic requirements so that the concept-forming operators (defined using the same syntactical form than the extension and intension operators of multi-adjoint concept lattices) form a Galois connection. As a consequence, given a relational database, we have much more possibilities to construct concept lattices associated with it, so that we can choose the specific version which better suits the situation.
Adjoint-Based Algorithms for Adaptation and Design Optimizations on Unstructured Grids
NASA Technical Reports Server (NTRS)
Nielsen, Eric J.
2006-01-01
Schemes based on discrete adjoint algorithms present several exciting opportunities for significantly advancing the current state of the art in computational fluid dynamics. Such methods provide an extremely efficient means for obtaining discretely consistent sensitivity information for hundreds of design variables, opening the door to rigorous, automated design optimization of complex aerospace configuration using the Navier-Stokes equation. Moreover, the discrete adjoint formulation provides a mathematically rigorous foundation for mesh adaptation and systematic reduction of spatial discretization error. Error estimates are also an inherent by-product of an adjoint-based approach, valuable information that is virtually non-existent in today's large-scale CFD simulations. An overview of the adjoint-based algorithm work at NASA Langley Research Center is presented, with examples demonstrating the potential impact on complex computational problems related to design optimization as well as mesh adaptation.
Discrete Adjoint-Based Design Optimization of Unsteady Turbulent Flows on Dynamic Unstructured Grids
NASA Technical Reports Server (NTRS)
Nielsen, Eric J.; Diskin, Boris; Yamaleev, Nail K.
2009-01-01
An adjoint-based methodology for design optimization of unsteady turbulent flows on dynamic unstructured grids is described. The implementation relies on an existing unsteady three-dimensional unstructured grid solver capable of dynamic mesh simulations and discrete adjoint capabilities previously developed for steady flows. The discrete equations for the primal and adjoint systems are presented for the backward-difference family of time-integration schemes on both static and dynamic grids. The consistency of sensitivity derivatives is established via comparisons with complex-variable computations. The current work is believed to be the first verified implementation of an adjoint-based optimization methodology for the true time-dependent formulation of the Navier-Stokes equations in a practical computational code. Large-scale shape optimizations are demonstrated for turbulent flows over a tiltrotor geometry and a simulated aeroelastic motion of a fighter jet.
Application of Adjoint Methodology in Various Aspects of Sonic Boom Design
NASA Technical Reports Server (NTRS)
Rallabhandi, Sriram K.
2014-01-01
One of the advances in computational design has been the development of adjoint methods allowing efficient calculation of sensitivities in gradient-based shape optimization. This paper discusses two new applications of adjoint methodology that have been developed to aid in sonic boom mitigation exercises. In the first, equivalent area targets are generated using adjoint sensitivities of selected boom metrics. These targets may then be used to drive the vehicle shape during optimization. The second application is the computation of adjoint sensitivities of boom metrics on the ground with respect to parameters such as flight conditions, propagation sampling rate, and selected inputs to the propagation algorithms. These sensitivities enable the designer to make more informed selections of flight conditions at which the chosen cost functionals are less sensitive.
MS S4.03.002 - Adjoint-Based Design for Configuration Shaping
NASA Technical Reports Server (NTRS)
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.
NASA Astrophysics Data System (ADS)
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.
NASA Astrophysics Data System (ADS)
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.
The Θ-KMS adjoint and time reversed quantum Markov semigroups
NASA Astrophysics Data System (ADS)
Bolaños-Servin, Jorge R.; Quezada, Roberto
2015-08-01
We introduce the notion of Θ-KMS adjoint of a quantum Markov semigroup, which is identified with the time reversed semigroup. The break of Θ-KMS symmetry, or Θ-standard quantum detailed balance in the sense of Fagnola-Umanità,11 is measured by means of the von Neumann relative entropy of states associated with the semigroup and its Θ-KMS adjoint.
Preliminary Results from the Application of Automated Adjoint Code Generation to CFL3D
NASA Technical Reports Server (NTRS)
Carle, Alan; Fagan, Mike; Green, Lawrence L.
1998-01-01
This report describes preliminary results obtained using an automated adjoint code generator for Fortran to augment a widely-used computational fluid dynamics flow solver to compute derivatives. These preliminary results with this augmented code suggest that, even in its infancy, the automated adjoint code generator can accurately and efficiently deliver derivatives for use in transonic Euler-based aerodynamic shape optimization problems with hundreds to thousands of independent design variables.
Global adjoint tomography: Perspectives, initial results and future directions
NASA Astrophysics Data System (ADS)
Bozdag, Ebru; Zhu, Hejun; Peter, Daniel; Tromp, Jeroen
2013-04-01
Adjoint methods provide an efficient way for incorporating the full nonlinearity of wave propagation and 3D Fréchet kernels in iterative seismic inversions. Our goal is to take adjoint tomography forward to image the entire planet using the opportunities offered by advances in numerical wave propagation solvers and high-performance computing. Using an iterative pre-conditioned conjugate gradient scheme, we initially set the aim to obtain a global crustal and mantle model with confined transverse isotropy in the upper mantle. Our strategy is to invert crustal and mantle structure together to avoid any bias introduced into upper-mantle images due to "crustal corrections", which are commonly used in classical tomography. We have started with 255 global CMT events (5.8 ≤ Mw ≤ 7) and used GSN stations as well as some local networks such as USArray, European stations, etc. We have demonstrated the feasibility of global scale inversions by performing two iterations based on numerical simulations accurate down to ~27 s. To simplify the problem, we primarily focus on elastic structure, and therefore our measurements are based on multitaper traveltime differences between observed and synthetic seismograms. We compute 3D sensitivity kernels for the selected events combining long-period surface waves (initially T > 60 s), where it is easier to handle nonlinearities due to the crust, with shorter-period body waves (initially T > 27 s), which are more sensitive to deeper parts of the mantle. 3D simulations dramatically increase the usable amount of data so that, with the current earthquake-station setup, we perform each iteration with more than two million measurements. Our initial results are promising to improve images from the upper mantle all the way down to the core-mantle boundary. Recent improvements in our 3D solvers (e.g., a GPU version) and access to high-performance computational centers (e.g., ORNL's Cray XK7 "Titan" system) now enable us to perform iterations
A Generalized Adjoint Approach for Quantifying Reflector Assembly Discontinuity Factor Uncertainties
Yankov, Artem; Collins, Benjamin; Jessee, Matthew Anderson; Downar, Thomas
2012-01-01
Sensitivity-based uncertainty analysis of assembly discontinuity factors (ADFs) can be readily performed using adjoint methods for infinite lattice models. However, there is currently no adjoint-based methodology to obtain uncertainties for ADFs along an interface between a fuel and reflector region. To accommodate leakage effects in a reflector region, a 1D approximation is usually made in order to obtain the homogeneous interface flux required to calculate the ADF. Within this 1D framework an adjoint-based method is proposed that is capable of efficiently calculating ADF uncertainties. In the proposed method the sandwich rule is utilized to relate the covariance of the input parameters of 1D diffusion theory in the reflector region to the covariance of the interface ADFs. The input parameters covariance matrix can be readily obtained using sampling-based codes such as XSUSA or adjoint-based codes such as TSUNAMI. The sensitivity matrix is constructed using a fixed-source adjoint approach for inputs characterizing the reflector region. An analytic approach is then used to determine the sensitivity of the ADFs to fuel parameters using the neutron balance equation. A stochastic approach is used to validate the proposed adjoint-based method.
Adjoint Data Assimilative Model Study of the Gulf of Maine Coastal Circulation
NASA Astrophysics Data System (ADS)
He, R.; McGillicuddy, D. J.; Lynch, D. R.
2004-12-01
Data assimilation (DA) in the coastal ocean can be divided into category of either sequential estimation or variational adjoint. Sequential estimation techniques blend models with observations directly, using a variety of algorithms with which the relative weights of data and model are calculated. Variational adjoint techniques infer model control variables (e.g. parameters, forcing functions, boundary conditions, etc.) that minimize the misfit between observations and predictions. The advantage of the latter techniques over the former is that the resulting model solutions obey model dynamics. In this study, the Gulf of Maine coastal circulation and the material property transport are investigated with the Dartmouth variational adjoint DA modeling system, which assimilates in-situ data via inversion for the unknown sea level elevation at open boundaries. In-situ observations include ADCP currents and coastal sea levels. The adjoint DA model skill is evaluated by the inter-comparisons between modeled and observed drifter trajectories. Excellent model skill is found, demonstrating the utility and effectiveness of the adjoint DA modeling system in bridging in-situ observations with coastal ocean model simulations. Implications of the adjoint DA strategy on the emergent coastal ocean observing systems are discussed.
Plumes, Hotspot & Slabs Imaged by Global Adjoint Tomography
NASA Astrophysics Data System (ADS)
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
A practical discrete-adjoint method for high-fidelity compressible turbulence simulations
NASA Astrophysics Data System (ADS)
Vishnampet, Ramanathan; Bodony, Daniel J.; Freund, Jonathan B.
2015-03-01
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 its
A practical discrete-adjoint method for high-fidelity compressible turbulence simulations
Vishnampet, Ramanathan; Bodony, Daniel J.; Freund, Jonathan B.
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
Adjoint sensitivity analysis of hydrodynamic stability in cyclonic flows
NASA Astrophysics Data System (ADS)
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.
Periodic differential equations with self-adjoint monodromy operator
NASA Astrophysics Data System (ADS)
Yudovich, V. I.
2001-04-01
A linear differential equation \\dot u=A(t)u with p-periodic (generally speaking, unbounded) operator coefficient in a Euclidean or a Hilbert space \\mathbb H is considered. It is proved under natural constraints that the monodromy operator U_p is self-adjoint and strictly positive if A^*(-t)=A(t) for all t\\in\\mathbb R.It is shown that Hamiltonian systems in the class under consideration are usually unstable and, if they are stable, then the operator U_p reduces to the identity and all solutions are p-periodic.For higher frequencies averaged equations are derived. Remarkably, high-frequency modulation may double the number of critical values.General results are applied to rotational flows with cylindrical components of the velocity a_r=a_z=0, a_\\theta=\\lambda c(t)r^\\beta, \\beta<-1, c(t) is an even p-periodic function, and also to several problems of free gravitational convection of fluids in periodic fields.
Conformal versus confining scenario in SU(2) with adjoint fermions
Del Debbio, L.; Pica, C.; Lucini, B.; Patella, A.; Rago, A.
2009-10-01
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 {beta}=4/g{sub 0}{sup 2}=2.25 for values of the bare fermion mass m{sub 0} 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.
Adjoint tomography of crust and upper-mantle structure beneath Continental China
NASA Astrophysics Data System (ADS)
Chen, M.; Niu, F.; Liu, Q.; Tromp, J.
2013-12-01
Four years of regional earthquake recordings from 1,869 seismic stations are used for high-resolution and high-fidelity seismic imaging of the crust and upper-mantle structure beneath Continental China. This unprecedented high-density dataset is comprised of seismograms recorded by the China Earthquake Administration Array (CEArray), NorthEast China Extended SeiSmic Array (NECESSArray), INDEPTH-IV Array, F-net and other global and regional seismic networks, and involves 1,326,384 frequency-dependent phase measurements. Adjoint tomography is applied to this unprecedented dataset, aiming to resolve detailed 3D maps of compressional and shear wavespeeds, and radial anisotropy. Contrary to traditional ray-theory based tomography, adjoint tomography takes into account full 3D wave propagation effects and off-ray-path sensitivity. In our implementation, it utilizes a spectral-element method for precise wave propagation simulations. The tomographic method starts with a 3D initial model that combines smooth radially anisotropic mantle model S362ANI and 3D crustal model Crust2.0. Traveltime and amplitude misfits are minimized iteratively based on a conjugate gradient method, harnessing 3D finite-frequency kernels computed for each updated 3D model. After 17 iterations, our inversion reveals strong correlations of 3D wavespeed heterogeneities in the crust and upper mantle with surface tectonic units, such as the Himalaya Block, the Tibetan Plateau, the Tarim Basin, the Ordos Block, and the South China Block. Narrow slab features emerge from the smooth initial model above the transition zone beneath the Japan, Ryukyu, Philippine, Izu-Bonin, Mariana and Andaman arcs. 3D wavespeed variations appear comparable to or much sharper than in high-frequency P-and S-wave models from previous studies. Moreover our results include new information, such as 3D variations of radial anisotropy and the Vp/Vs ratio, which are expected to shed new light to the composition, thermal state, flow
NASA Astrophysics Data System (ADS)
Komatitsch, Dimitri; Xie, Zhinan; Bozdağ, Ebru; Sales de Andrade, Elliott; Peter, Daniel; Liu, Qinya; Tromp, Jeroen
2016-06-01
We introduce a technique to compute exact anelastic sensitivity kernels in the time domain using parsimonious disk storage. The method is based on a reordering of the time loop of time-domain forward/adjoint wave propagation solvers combined with the use of a memory buffer. It avoids instabilities that occur when time-reversing dissipative wave propagation simulations. The total number of required time steps is unchanged compared to usual acoustic or elastic approaches. The cost is reduced by a factor of 4/3 compared to the case in which anelasticity is partially accounted for by accommodating the effects of physical dispersion. We validate our technique by performing a test in which we compare the Kα sensitivity kernel to the exact kernel obtained by saving the entire forward calculation. This benchmark confirms that our approach is also exact. We illustrate the importance of including full attenuation in the calculation of sensitivity kernels by showing significant differences with physical-dispersion-only kernels.
Focus point gauge mediation with incomplete adjoint messengers and gauge coupling unification
NASA Astrophysics Data System (ADS)
Bhattacharyya, Gautam; Yanagida, Tsutomu T.; Yokozaki, Norimi
2015-10-01
As the mass limits on supersymmetric particles are gradually pushed to higher values due to their continuing non-observation at the CERN LHC, looking for focus point regions in the supersymmetric parameter space, which shows considerably reduced fine-tuning, is increasingly more important than ever. We explore this in the context of gauge mediated supersymmetry breaking with messengers transforming in the adjoint representation of the gauge group, namely, octet of color SU(3) and triplet of weak SU(2). A distinctive feature of this scenario is that the focus point is achieved by fixing a single combination of parameters in the messenger sector, which is invariant under the renormalization group evolution. Because of this invariance, the focus point behavior is well under control once the relevant parameters are fixed by a more fundamental theory. The observed Higgs boson mass is explained with a relatively mild fine-tuning Δ = 60- 150. Interestingly, even in the presence of incomplete messenger multiplets of the SU(5) GUT group, the gauge couplings still unify perfectly, but at a scale which is one or two orders of magnitude above the conventional GUT scale. Because of this larger unification scale, the colored Higgs multiplets become too heavy to trigger proton decay at a rate larger than the experimentally allowed limit.
Adjoint-based optimization for the understanding of the aerodynamics of a flapping plate
NASA Astrophysics Data System (ADS)
Wei, Mingjun; Xu, Min
2015-11-01
An adjoint-based optimization is applied on a rigid flapping plate and a flexible flapping plate for drag reduction and for propulsive efficiency. Non-cylindrical calculus is introduced to handle the moving boundary. The rigid plate has a combined plunging and pitching motion with incoming flow, the control parameter is the phase delay which is considered first as a constant then as an arbitrary time-varying function. The optimal controls with different cost functions provide different strategies to reach maximum drag reduction or propulsive efficiency. The flexible plate has plunging, pitching, and deformation which is defined by the first two natural modes. With the same optimization goals, the control is instead the amplitude and phase delay of the pitching, the first eigen mode, and the second eigen mode. Similar analyses are taken to understand the conditions for drag reduction and propulsive efficiency when flexibility is involved. It is also shown that the flexibility plays a more important role at lower Reynolds number. Supported by AFOSR.
NASA Astrophysics Data System (ADS)
Shadid, J. N.; Smith, T. M.; Cyr, E. C.; Wildey, T. M.; Pawlowski, R. P.
2016-09-01
A critical aspect of applying modern computational solution methods to complex multiphysics systems of relevance to nuclear reactor modeling, is the assessment of the predictive capability of specific proposed mathematical models. In this respect the understanding of numerical error, the sensitivity of the solution to parameters associated with input data, boundary condition uncertainty, and mathematical models is critical. Additionally, the ability to evaluate and or approximate the model efficiently, to allow development of a reasonable level of statistical diagnostics of the mathematical model and the physical system, is of central importance. In this study we report on initial efforts to apply integrated adjoint-based computational analysis and automatic differentiation tools to begin to address these issues. The study is carried out in the context of a Reynolds averaged Navier-Stokes approximation to turbulent fluid flow and heat transfer using a particular spatial discretization based on implicit fully-coupled stabilized FE methods. Initial results are presented that show the promise of these computational techniques in the context of nuclear reactor relevant prototype thermal-hydraulics problems.
NASA Astrophysics Data System (ADS)
Ito, Shin-Ichi; Nagao, Hiromichi; Yamanaka, Akinori; Tsukada, Yuhki; Koyama, Toshiyuki; Inoue, Junya
Phase field (PF) method, which phenomenologically describes dynamics of microstructure evolutions during solidification and phase transformation, has progressed in the fields of hydromechanics and materials engineering. How to determine, based on observation data, an initial state and model parameters involved in a PF model is one of important issues since previous estimation methods require too much computational cost. We propose data assimilation (DA), which enables us to estimate the parameters and states by integrating the PF model and observation data on the basis of the Bayesian statistics. The adjoint method implemented on DA not only finds an optimum solution by maximizing a posterior distribution but also evaluates the uncertainty in the estimations by utilizing the second order information of the posterior distribution. We carried out an estimation test using synthetic data generated by the two-dimensional Kobayashi's PF model. The proposed method is confirmed to reproduce the true initial state and model parameters we assume in advance, and simultaneously estimate their uncertainties due to quality and quantity of the data. This result indicates that the proposed method is capable of suggesting the experimental design to achieve the required accuracy.
Adjoint-based airfoil shape optimization in transonic flow
NASA Astrophysics Data System (ADS)
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.
Generalized adjoint consistent treatment of wall boundary conditions for compressible flows
NASA Astrophysics Data System (ADS)
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.
Self-adjoint Operators as Functions I. Lattices, Galois Connections, and the Spectral Order
NASA Astrophysics Data System (ADS)
Döring, Andreas; Dewitt, Barry
2014-06-01
Observables of a quantum system, described by self-adjoint operators in a von Neumann algebra or affiliated with it in the unbounded case, form a conditionally complete lattice when equipped with the spectral order. Using this order-theoretic structure, we develop a new perspective on quantum observables. In this first paper (of two), we show that self-adjoint operators affiliated with a von Neumann algebra can equivalently be described as certain real-valued functions on the projection lattice of the algebra, which we call q-observable functions. Bounded self-adjoint operators correspond to q-observable functions with compact image on non-zero projections. These functions, originally defined in a similar form by de Groote (Observables II: quantum observables, 2005), are most naturally seen as adjoints (in the categorical sense) of spectral families. We show how they relate to the daseinisation mapping from the topos approach to quantum theory (Döring and Isham , New Structures for Physics, Springer, Heidelberg, 2011). Moreover, the q-observable functions form a conditionally complete lattice which is shown to be order-isomorphic to the lattice of self-adjoint operators with respect to the spectral order. In a subsequent paper (Döring and Dewitt, 2012, preprint), we will give an interpretation of q-observable functions in terms of quantum probability theory, and using results from the topos approach to quantum theory, we will provide a joint sample space for all quantum observables.
NASA Astrophysics Data System (ADS)
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 (<10 pcm) obtained in these experiments, a reference calculation need to be achieved. This calculation may be accomplished using the continuous-energy Monte Carlo code 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
Using an Adjoint Approach to Eliminate Mesh Sensitivities in Computational Design
NASA Technical Reports Server (NTRS)
Nielsen, Eric J.; Park, Michael A.
2006-01-01
An algorithm for efficiently incorporating the effects of mesh sensitivities in a computational design framework is introduced. The method is based on an adjoint approach and eliminates the need for explicit linearizations of the mesh movement scheme with respect to the geometric parameterization variables, an expense that has hindered practical large-scale design optimization using discrete adjoint methods. The effects of the mesh sensitivities can be accounted for through the solution of an adjoint problem equivalent in cost to a single mesh movement computation, followed by an explicit matrix-vector product scaling with the number of design variables and the resolution of the parameterized surface grid. The accuracy of the implementation is established and dramatic computational savings obtained using the new approach are demonstrated using several test cases. Sample design optimizations are also shown.
Instantons and the 5D U(1) gauge theory with extra adjoint
NASA Astrophysics Data System (ADS)
Poghossian, Rubik; Samsonyan, Marine
2009-07-01
In this paper, we compute the partition function of 5D supersymmetric U(1) gauge theory with extra adjoint matter in general Ω background. It is well known that such partition functions encode very rich topological information. We show in particular that unlike the case with no extra matter, the partition function with extra adjoint at some special values of the parameters directly reproduces the generating function for the Poincare polynomial of the moduli space of instantons. We compare our results with those recently obtained by Iqbal et al (Refined topological vertex, cylindric partitions and the U(1) adjoint theory, arXiv:0803.2260), who used the so-called refined topological vertex method.
Using an Adjoint Approach to Eliminate Mesh Sensitivities in Computational Design
NASA Technical Reports Server (NTRS)
Nielsen, Eric J.; Park, Michael A.
2005-01-01
An algorithm for efficiently incorporating the effects of mesh sensitivities in a computational design framework is introduced. The method is based on an adjoint approach and eliminates the need for explicit linearizations of the mesh movement scheme with respect to the geometric parameterization variables, an expense that has hindered practical large-scale design optimization using discrete adjoint methods. The effects of the mesh sensitivities can be accounted for through the solution of an adjoint problem equivalent in cost to a single mesh movement computation, followed by an explicit matrix-vector product scaling with the number of design variables and the resolution of the parameterized surface grid. The accuracy of the implementation is established and dramatic computational savings obtained using the new approach are demonstrated using several test cases. Sample design optimizations are also shown.
Tracking influential haze source areas in North China using an adjoint model, GRAPES-CUACE
NASA Astrophysics Data System (ADS)
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.
NASA Technical Reports Server (NTRS)
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.
Stratospheric Water Vapor and the Asian Monsoon: An Adjoint Model Investigation
NASA Technical Reports Server (NTRS)
Olsen, Mark A.; Andrews, Arlyn E.
2003-01-01
A new adjoint model of the Goddard Parameterized Chemistry and Transport Model is used to investigate the role that the Asian monsoon plays in transporting water to the stratosphere. The adjoint model provides a unique perspective compared to non-diffusive and non-mixing Lagrangian trajectory analysis. The quantity of water vapor transported from the monsoon and the pathways into the stratosphere are examined. The emphasis is on the amount of water originating from the monsoon that contributes to the tropical tape recorder signal. The cross-tropopause flux of water from the monsoon to the midlatitude lower stratosphere will also be discussed.
NEMOTAM: tangent and adjoint models for the ocean modelling platform NEMO
NASA Astrophysics Data System (ADS)
Vidard, A.; Bouttier, P.-A.; Vigilant, F.
2014-10-01
The tangent linear and adjoint model (TAM) are efficient tools to analyse and to control dynamical systems such as NEMO. They can be involved in a large range of applications such as sensitivity analysis, parameter estimation or the computation of characteristics vectors. TAM is also required by the 4-D-VAR algorithm which is one of the major method in Data Assimilation. This paper describes the development and the validation of the Tangent linear and Adjoint Model for the NEMO ocean modelling platform (NEMOTAM). The diagnostic tools that are available alongside NEMOTAM are detailed and discussed and several applications are also presented.
Global Adjoint Tomography: Combining Big Data with HPC Simulations
NASA Astrophysics Data System (ADS)
Bozdag, E.; Lefebvre, M. P.; Lei, W.; Peter, D. B.; Smith, J. A.; Komatitsch, D.; Tromp, J.
2014-12-01
The steady increase in data quality and the number of global seismographic stations have substantially grown the amount of data available for construction of Earth models. Meanwhile, developments in the theory of wave propagation, numerical methods and HPC systems have enabled unprecedented simulations of seismic wave propagation in realistic 3D Earth models which lead the extraction of more information from data, ultimately culminating in the use of entire three-component seismograms.Our aim is to take adjoint tomography further to image the entire planet which is one of the extreme cases in seismology due to its intense computational requirements and vast amount of high-quality seismic data that can potentially be assimilated in inversions. We have started low resolution (T > 27 s, soon will be > 17 s) global inversions with 253 earthquakes for a transversely isotropic crust and mantle model on Oak Ridge National Laboratory's Cray XK7 "Titan" system. Recent improvements in our 3D solvers, such as the GPU version of the SPECFEM3D_GLOBE package, will allow us perform higher-resolution (T > 9 s) and longer-duration (~180 m) simulations to take the advantage of high-frequency body waves and major-arc surface waves to improve imbalanced ray coverage as a result of uneven distribution of sources and receivers on the globe. Our initial results after 10 iterations already indicate several prominent features reported in high-resolution continental studies, such as major slabs (Hellenic, Japan, Bismarck, Sandwich, etc.) and enhancement in plume structures (the Pacific superplume, the Hawaii hot spot, etc.). Our ultimate goal is to assimilate seismic data from more than 6,000 earthquakes within the magnitude range 5.5 ≤ Mw ≤ 7.0. To take full advantage of this data set on ORNL's computational resources, we need a solid framework for managing big data sets during pre-processing (e.g., data requests and quality checks), gradient calculations, and post-processing (e
Adaptive mesh refinement and adjoint methods in geophysics simulations
NASA Astrophysics Data System (ADS)
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.
NASA Astrophysics Data System (ADS)
Wang, Brian; Goldstein, Moshe; Xu, X. George; Sahoo, Narayan
2005-03-01
Recently, the theoretical framework of the adjoint Monte Carlo (AMC) method has been developed using a simplified patient geometry. In this study, we extended our previous work by applying the AMC framework to a 3D anatomical model called VIP-Man constructed from the Visible Human images. First, the adjoint fluxes for the prostate (PTV) and rectum and bladder (organs at risk (OARs)) were calculated on a spherical surface of 1 m radius, centred at the centre of gravity of PTV. An importance ratio, defined as the PTV dose divided by the weighted OAR doses, was calculated for each of the available beamlets to select the beam angles. Finally, the detailed doses in PTV and OAR were calculated using a forward Monte Carlo simulation to include the electron transport. The dose information was then used to generate dose volume histograms (DVHs). The Pinnacle treatment planning system was also used to generate DVHs for the 3D plans with beam angles obtained from the AMC (3D-AMC) and a standard six-field conformal radiation therapy plan (3D-CRT). Results show that the DVHs for prostate from 3D-AMC and the standard 3D-CRT are very similar, showing that both methods can deliver prescribed dose to the PTV. A substantial improvement in the DVHs for bladder and rectum was found for the 3D-AMC method in comparison to those obtained from 3D-CRT. However, the 3D-AMC plan is less conformal than the 3D-CRT plan because only bladder, rectum and PTV are considered for calculating the importance ratios. Nevertheless, this study clearly demonstrated the feasibility of the AMC in selecting the beam directions as a part of a treatment planning based on the anatomical information in a 3D and realistic patient anatomy.
Estimation of ex-core detector responses by adjoint Monte Carlo
Hoogenboom, J. E.
2006-07-01
Ex-core detector responses can be efficiently calculated by combining an adjoint Monte Carlo calculation with the converged source distribution of a forward Monte Carlo calculation. As the fission source distribution from a Monte Carlo calculation is given only as a collection of discrete space positions, the coupling requires a point flux estimator for each collision in the adjoint calculation. To avoid the infinite variance problems of the point flux estimator, a next-event finite-variance point flux estimator has been applied, witch is an energy dependent form for heterogeneous media of a finite-variance estimator known from the literature. To test the effects of this combined adjoint-forward calculation a simple geometry of a homogeneous core with a reflector was adopted with a small detector in the reflector. To demonstrate the potential of the method the continuous-energy adjoint Monte Carlo technique with anisotropic scattering was implemented with energy dependent absorption and fission cross sections and constant scattering cross section. A gain in efficiency over a completely forward calculation of the detector response was obtained, which is strongly dependent on the specific system and especially the size and position of the ex-core detector and the energy range considered. Further improvements are possible. The method works without problems for small detectors, even for a point detector and a small or even zero energy range. (authors)
Eulerian-Lagrangian localized adjoint methods for reactive transport in groundwater
Ewing, R.E.; Wang, Hong
1996-12-31
In this paper, we present Eulerian-Lagrangian localized adjoint methods (ELLAM) to solve convection-diffusion-reaction equations governing contaminant transport in groundwater flowing through an adsorbing porous medium. These ELLAM schemes can treat various combinations of boundary conditions and conserve mass. Numerical results are presented to demonstrate the strong potential of ELLAM schemes.
Sensitivity analysis of a model of CO2 exchange in tundra ecosystems by the adjoint method
NASA Technical Reports Server (NTRS)
Waelbroek, C.; Louis, J.-F.
1995-01-01
A model of net primary production (NPP), decomposition, and nitrogen cycling in tundra ecosystems has been developed. The adjoint technique is used to study the sensitivity of the computed annual net CO2 flux to perturbation in initial conditions, climatic inputs, and model's main parameters describing current seasonal CO2 exchange in wet sedge tundra at Barrow, Alaska. The results show that net CO2 flux is most sensitive to parameters characterizing litter chemical composition and more sensitive to decomposition parameters than to NPP parameters. This underlines the fact that in nutrient-limited ecosystems, decomposition drives net CO2 exchange by controlling mineralization of main nutrients. The results also indicate that the short-term (1 year) response of wet sedge tundra to CO2-induced warming is a significant increase in CO2 emission, creating a positive feedback to atmosphreic CO2 accumulation. However, a cloudiness increase during the same year can severely alter this response and lead to either a slight decrease or a strong increase in emitted CO2, depending on its exact timing. These results demonstrate that the adjoint method is well suited to study systems encountering regime changes, as a single run of the adjoint model provides sensitivities of the net CO2 flux to perturbations in all parameters and variables at any time of the year. Moreover, it is shown that large errors due to the presence of thresholds can be avoided by first delimiting the range of applicability of the adjoint results.
Analysis of Correlated Coupling of Monte Carlo Forward and Adjoint Histories
Ueki, Taro; Hoogenboom, J.E.; Kloosterman, J. L.
2001-02-15
In Monte Carlo correlated coupling, forward and adjoint particle histories are initiated in exactly opposite directions at an arbitrarily placed surface between a physical source and a physical detector. It is shown that this coupling calculation can become more efficient than standard forward calculations. In many cases, the basic form of correlated coupling is less efficient than standard forward calculations. This inherent inefficiency can be overcome by applying a black absorber perturbation to either the forward or the adjoint problem and by processing the product of batch averages as one statistical entity. The usage of the black absorber is based on the invariance of the response flow integral with a material perturbation in either the physical detector side volume in the forward problem or the physical source side volume in the adjoint problem. The batch-average product processing makes use of a quadratic increase of the nonzero coupled-score probability. All the developments have been done in such a way that improved efficiency schemes available in widely distributed Monte Carlo codes can be applied to both the forward and adjoint simulations. Also, the physical meaning of the black absorber perturbation is interpreted based on surface crossing and is numerically validated. In addition, the immediate reflection at the intermediate surface with a controlled direction change is investigated within the invariance framework. This approach can be advantageous for a void streaming problem.
Canonical quantization of lattice Higgs-Maxwell-Chern-Simons fields: Krein Self-adjointness
Bowman, Daniel A.; Challifour, John L.
2006-10-15
It is shown how techniques from constructive quantum field theory may be applied to indefinite metric gauge theories in Hilbert space for the case of a Higgs-Maxwell-Chern-Simons theory on a lattice. The Hamiltonian operator is shown to be Krein essentially self-adjoint by means of unbounded but Krein unitary transformations relating the Hamiltonian to an essentially maximal accretive operator.
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.
NASA Astrophysics Data System (ADS)
Ren, L.; Liu, Q.
2012-12-01
We present multiple moment-tensor solution of the December 26, 2004 Sumatra earthquake based upon adjoint methods. An objective function Φ that measures the goodness of waveform fit between data and synthetics is minimized. Synthetics are calculated by spectral-element simulations (SPECFEM3D_GLOBE) in a 3D global earth model S362ANI to reduce the effect of heterogeneous structures. The Fréchet derivatives of Φ in the form δΦ = ∫T ∫VI(ɛ †ij)(x,T-t) δ(m_dot)ij(x,t)d3xdt, where δmij is the perturbation of moment density function and I(ɛ†ij)(x,T-t) denotes the time-integrated adjoint strain tensor, are calculated based upon adjoint methods implemented in SPECFEM3D_GLOBE. Our initial source model is obtained by monitoring the time-integrated adjoint strain tensors in the vicinity of the presumed source region. Source model parameters are iteratively updated by a preconditioned conjugate-gradient method to iteratively utilizing the calculated Φ and δΦ values. Our final inversion results show both similarities to and differences from previous source inversion results based on 1D background models.
NASA Astrophysics Data System (ADS)
Ren, L.; Liu, Q.; Hjörleifsdóttir, V.
2010-12-01
We present multiple moment-tensor solution of the Dec 26, 2004 Sumatra earthquake based upon the adjoint methods. An objective function Φ(m), where m is the multiple source model, measures the goodness of waveform fit between data and synthetics. The Fréchet derivatives of Φ in the form δΦ = ∫∫I(ɛ†)(x,T-t)δmij_dot(x,t)dVdt, where δmij is the source model perturbation and I(ɛ†)(x,T-t) denotes the time-integrated adjoint strain tensor, are calculated based upon adjoint methods and spectral-element simulations (SPECFEM3D_GLOBE) in a 3D global earth model S362ANI. Our initial source model is obtained independently by monitoring the time-integrated adjoint strain tensors around the presumed source region. We then utilize the Φ and δΦ calculations in a conjugate-gradient method to iteratively invert for the source model. Our final inversion results show both similarities with and differences to previous source inversion results based on 1D earth models.
On the proper treatment of grid sensitivities in continuous adjoint methods for shape optimization
NASA Astrophysics Data System (ADS)
Kavvadias, I. S.; Papoutsis-Kiachagias, E. M.; Giannakoglou, K. C.
2015-11-01
The continuous adjoint method for shape optimization problems, in flows governed by the Navier-Stokes equations, can be formulated in two different ways, each of which leads to a different expression for the sensitivity derivatives of the objective function with respect to the control variables. The first formulation leads to an expression including only boundary integrals; it, thus, has low computational cost but, when used with coarse grids, its accuracy becomes questionable. The second formulation comprises a sum of boundary and field integrals; due to the field integrals, it has noticeably higher computational cost, obtaining though higher accuracy. In this paper, the equivalence of the two formulations is revisited from the mathematical and, particularly, the numerical point of view. Internal and external aerodynamics cases, in which the objective function is either the total pressure losses or the force exerted on a solid body, are examined and differences in the computed gradients are discussed. After identifying the reason behind these discrepancies, the adjoint formulation is enhanced by the adjoint to a (hypothetical) grid displacement model and the new approach is proved to reproduce the accuracy of the second adjoint formulation while maintaining the low cost of the first one.
NASA Astrophysics Data System (ADS)
Peter, Daniel; Videau, Brice; Pouget, Kevin; Komatitsch, Dimitri
2015-04-01
Improving the resolution of tomographic images is crucial to answer important questions on the nature of Earth's subsurface structure and internal processes. Seismic tomography is the most prominent approach where seismic signals from ground-motion records are used to infer physical properties of internal structures such as compressional- and shear-wave speeds, anisotropy and attenuation. Recent advances in regional- and global-scale seismic inversions move towards full-waveform inversions which require accurate simulations of seismic wave propagation in complex 3D media, providing access to the full 3D seismic wavefields. However, these numerical simulations are computationally very expensive and need high-performance computing (HPC) facilities for further improving the current state of knowledge. During recent years, many-core architectures such as graphics processing units (GPUs) have been added to available large HPC systems. Such GPU-accelerated computing together with advances in multi-core central processing units (CPUs) can greatly accelerate scientific applications. There are mainly two possible choices of language support for GPU cards, the CUDA programming environment and OpenCL language standard. CUDA software development targets NVIDIA graphic cards while OpenCL was adopted mainly by AMD graphic cards. In order to employ such hardware accelerators for seismic wave propagation simulations, we incorporated a code generation tool BOAST into an existing spectral-element code package SPECFEM3D_GLOBE. This allows us to use meta-programming of computational kernels and generate optimized source code for both CUDA and OpenCL languages, running simulations on either CUDA or OpenCL hardware accelerators. We show here applications of forward and adjoint seismic wave propagation on CUDA/OpenCL GPUs, validating results and comparing performances for different simulations and hardware usages.
NASA Astrophysics Data System (ADS)
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.
Aerosol Health Impact Source Attribution Studies with the CMAQ Adjoint Air Quality Model
NASA Astrophysics Data System (ADS)
Turner, M. D.
Fine particulate matter (PM2.5) is an air pollutant consisting of a mixture of solid and liquid particles suspended in the atmosphere. Knowledge of the sources and distributions of PM2.5 is important for many reasons, two of which are that PM2.5 has an adverse effect on human health and also an effect on climate change. Recent studies have suggested that health benefits resulting from a unit decrease in black carbon (BC) are four to nine times larger than benefits resulting from an equivalent change in PM2.5 mass. The goal of this thesis is to quantify the role of emissions from different sectors and different locations in governing the total health impacts, risk, and maximum individual risk of exposure to BC both nationally and regionally in the US. We develop and use the CMAQ adjoint model to quantify the role of emissions from all modeled sectors, times, and locations on premature deaths attributed to exposure to BC. From a national analysis, we find that damages resulting from anthropogenic emissions of BC are strongly correlated with population and premature death. However, we find little correlation between damages and emission magnitude, suggesting that controls on the largest emissions may not be the most efficient means of reducing damages resulting from BC emissions. Rather, the best proxy for locations with damaging BC emissions is locations where premature deaths occur. Onroad diesel and nonroad vehicle emissions are the largest contributors to premature deaths attributed to exposure to BC, while onroad gasoline emissions cause the highest deaths per amount emitted. Additionally, emissions in fall and winter contribute to more premature deaths (and more per amount emitted) than emissions in spring and summer. From a regional analysis, we find that emissions from outside each of six urban areas account for 7% to 27% of the premature deaths attributed to exposure to BC within the region. Within the region encompassing New York City and Philadelphia
Non-self-adjoint Hamiltonians with complex eigenvalues
NASA Astrophysics Data System (ADS)
Bagarello, F.
2016-05-01
Motivated by what one observes dealing with PT-symmetric quantum mechanics, we discuss what happens if a physical system is driven by a diagonalizable Hamiltonian with not all real eigenvalues. In particular, we consider the functional structure related to systems living in finite-dimensional Hilbert spaces, and we show that certain intertwining relations can be deduced also in this case if we introduce suitable antilinear operators. We also analyze a simple model, computing the transition probabilities in the broken and in the unbroken regime.
Importance of a Specific Amino Acid Pairing for Murine MLL Leukemias Driven by MLLT1/3 or AFF1/4
Lokken, Alyson A.; Achille, Nicholas J.; Chang, Ming-jin; Lin, Jeffrey J.; Kuntimaddi, Aravinda; Leach, Benjamin I.; Malik, Bhavna; Nesbit, Jacqueline B.; Zhang, Shubin; Bushweller, John H.; Zeleznik-Le, Nancy J.; Hemenway, Charles S.
2014-01-01
Acute leukemias caused by translocations of the MLL gene at chromosome 11 band q23 (11q23) are characterized by a unique gene expression profile. More recently, data from several laboratories indicate that the most commonly encountered MLL fusion proteins, MLLT1, MLLT3, and AFF1 are found within a molecular complex that facilitates the elongation phase of mRNA transcription. Mutational analyses suggest that interaction between the MLLT1/3 proteins and AFF family proteins are required for experimental transformation of hematopoietic progenitor cells (HPCs). Here, we define a specific pairing of two amino acids that creates a salt bridge between MLLT1/3 and AFF proteins that is critically important for MLL-mediated transformation of HPCs. Our findings, coupled with the newly defined structure of MLLT3 in complex with AFF1, should facilitate the development of small molecules that block this amino acid interaction and interfere with the activity of the most common MLL oncoproteins. PMID:25282333
Chang, Fiona T M; Chan, F Lyn; R McGhie, James D; Udugama, Maheshi; Mayne, Lynne; Collas, Philippe; Mann, Jeffrey R; Wong, Lee H
2015-03-11
Human ALT cancers show high mutation rates in ATRX and DAXX. Although it is well known that the absence of ATRX/DAXX disrupts H3.3 deposition at heterochromatin, its impact on H3.3 deposition and post-translational modification in the global genome remains unclear. Here, we explore the dynamics of phosphorylated H3.3 serine 31 (H3.3S31ph) in human ALT cancer cells. While H3.3S31ph is found only at pericentric satellite DNA repeats during mitosis in most somatic human cells, a high level of H3.3S31ph is detected on the entire chromosome in ALT cells, attributable to an elevated CHK1 activity in these cells. Drug inhibition of CHK1 activity during mitosis and expression of mutant H3.3S31A in these ALT cells result in a decrease in H3.3S31ph levels accompanied with increased levels of phosphorylated H2AX serine 139 on chromosome arms and at the telomeres. Furthermore, the inhibition of CHK1 activity in these cells also reduces cell viability. Our findings suggest a novel role of CHK1 as an H3.3S31 kinase, and that CHK1-mediated H3.3S31ph plays an important role in the maintenance of chromatin integrity and cell survival in ALT cancer cells. PMID:25690891
Demonstration of Automatically-Generated Adjoint Code for Use in Aerodynamic Shape Optimization
NASA Technical Reports Server (NTRS)
Green, Lawrence; Carle, Alan; Fagan, Mike
1999-01-01
Gradient-based optimization requires accurate derivatives of the objective function and constraints. These gradients may have previously been obtained by manual differentiation of analysis codes, symbolic manipulators, finite-difference approximations, or existing automatic differentiation (AD) tools such as ADIFOR (Automatic Differentiation in FORTRAN). Each of these methods has certain deficiencies, particularly when applied to complex, coupled analyses with many design variables. Recently, a new AD tool called ADJIFOR (Automatic Adjoint Generation in FORTRAN), based upon ADIFOR, was developed and demonstrated. Whereas ADIFOR implements forward-mode (direct) differentiation throughout an analysis program to obtain exact derivatives via the chain rule of calculus, ADJIFOR implements the reverse-mode counterpart of the chain rule to obtain exact adjoint form derivatives from FORTRAN code. Automatically-generated adjoint versions of the widely-used CFL3D computational fluid dynamics (CFD) code and an algebraic wing grid generation code were obtained with just a few hours processing time using the ADJIFOR tool. The codes were verified for accuracy and were shown to compute the exact gradient of the wing lift-to-drag ratio, with respect to any number of shape parameters, in about the time required for 7 to 20 function evaluations. The codes have now been executed on various computers with typical memory and disk space for problems with up to 129 x 65 x 33 grid points, and for hundreds to thousands of independent variables. These adjoint codes are now used in a gradient-based aerodynamic shape optimization problem for a swept, tapered wing. For each design iteration, the optimization package constructs an approximate, linear optimization problem, based upon the current objective function, constraints, and gradient values. The optimizer subroutines are called within a design loop employing the approximate linear problem until an optimum shape is found, the design loop
NASA Astrophysics Data System (ADS)
Kavvadias, I. S.; Papoutsis-Kiachagias, E. M.; Dimitrakopoulos, G.; Giannakoglou, K. C.
2015-11-01
In this article, the gradient of aerodynamic objective functions with respect to design variables, in problems governed by the incompressible Navier-Stokes equations coupled with the k-ω SST turbulence model, is computed using the continuous adjoint method, for the first time. Shape optimization problems for minimizing drag, in external aerodynamics (flows around isolated airfoils), or viscous losses in internal aerodynamics (duct flows) are considered. Sensitivity derivatives computed with the proposed adjoint method are compared to those computed with finite differences or a continuous adjoint variant based on the frequently used assumption of frozen turbulence; the latter proves the need for differentiating the turbulence model. Geometries produced by optimization runs performed with sensitivities computed by the proposed method and the 'frozen turbulence' assumption are also compared to quantify the gain from formulating and solving the adjoint to the turbulence model equations.
Dietze, K.; Fort, E.; Rahlfs, S.; Rimpault, G.; Salvatores, M.
1994-12-31
The Rossendorf RRRJSEG configurations characterized by, energy-independent or continuously rising adjoint spectra have been recalculated using the full European scheme JEF2/ECCO/ERANOS. C/E-values are given for structural materials and fission product nuclides using the results of sample reactivity measurements at the central position of these configurations. Due to the specially designed adjoint spectra, capture or scattering cross- sections can be checked separately. Recommendations for data corrections are given based on perturbation theory calculations.
NASA Astrophysics Data System (ADS)
Gaševic, Dragan; Djuric, Dragan; Devedžic, Vladan
A relevant initiative from the software engineering community called Model Driven Engineering (MDE) is being developed in parallel with the Semantic Web (Mellor et al. 2003a). The MDE approach to software development suggests that one should first develop a model of the system under study, which is then transformed into the real thing (i.e., an executable software entity). The most important research initiative in this area is the Model Driven Architecture (MDA), which is Model Driven Architecture being developed under the umbrella of the Object Management Group (OMG). This chapter describes the basic concepts of this software engineering effort.
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
NASA Astrophysics Data System (ADS)
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.
Efficient Construction of Discrete Adjoint Operators on Unstructured Grids Using Complex Variables
NASA Technical Reports Server (NTRS)
Nielsen, Eric J.; Kleb, William L.
2005-01-01
A methodology is developed and implemented to mitigate the lengthy software development cycle typically associated with constructing a discrete adjoint solver for aerodynamic simulations. The approach is based on a complex-variable formulation that enables straightforward differentiation of complicated real-valued functions. An automated scripting process is used to create the complex-variable form of the set of discrete equations. An efficient method for assembling the residual and cost function linearizations is developed. The accuracy of the implementation is verified through comparisons with a discrete direct method as well as a previously developed handcoded discrete adjoint approach. Comparisons are also shown for a large-scale configuration to establish the computational efficiency of the present scheme. To ultimately demonstrate the power of the approach, the implementation is extended to high temperature gas flows in chemical nonequilibrium. Finally, several fruitful research and development avenues enabled by the current work are suggested.
NASA Technical Reports Server (NTRS)
Nielsen, Eric J.; Kleb, William L.
2005-01-01
A methodology is developed and implemented to mitigate the lengthy software development cycle typically associated with constructing a discrete adjoint solver for aerodynamic simulations. The approach is based on a complex-variable formulation that enables straightforward differentiation of complicated real-valued functions. An automated scripting process is used to create the complex-variable form of the set of discrete equations. An efficient method for assembling the residual and cost function linearizations is developed. The accuracy of the implementation is verified through comparisons with a discrete direct method as well as a previously developed handcoded discrete adjoint approach. Comparisons are also shown for a large-scale configuration to establish the computational efficiency of the present scheme. To ultimately demonstrate the power of the approach, the implementation is extended to high temperature gas flows in chemical nonequilibrium. Finally, several fruitful research and development avenues enabled by the current work are suggested.
Adjoint-based optimal control of an airfoil in gusting flows
NASA Astrophysics Data System (ADS)
Choi, Jeesoon; Colonius, Tim; California Institute of Technology Team
2015-11-01
In this study, we apply optimal control to an airfoil in gusting flow to investigate the possibility of extracting energy. The gradients of an objective function are obtained via the adjoint method and used to minimize the cost. The immersed boundary projection method is used for our forward solver, and the relevant adjoint equations are derived by the discrete-then-differentiate approach. Translational gusts are generated by a body force in the computational domain upstream to the body, and the method finds the optimal angles of the airfoil that exploits the greatest amount of energy. The influence of a vortex traversing an airfoil is also investigated and optimized to reduce the fluctuating lift.
An Exact Dual Adjoint Solution Method for Turbulent Flows on Unstructured Grids
NASA Technical Reports Server (NTRS)
Nielsen, Eric J.; Lu, James; Park, Michael A.; Darmofal, David L.
2003-01-01
An algorithm for solving the discrete adjoint system based on an unstructured-grid discretization of the Navier-Stokes equations is presented. The method is constructed such that an adjoint solution exactly dual to a direct differentiation approach is recovered at each time step, yielding a convergence rate which is asymptotically equivalent to that of the primal system. The new approach is implemented within a three-dimensional unstructured-grid framework and results are presented for inviscid, laminar, and turbulent flows. Improvements to the baseline solution algorithm, such as line-implicit relaxation and a tight coupling of the turbulence model, are also presented. By storing nearest-neighbor terms in the residual computation, the dual scheme is computationally efficient, while requiring twice the memory of the flow solution. The scheme is expected to have a broad impact on computational problems related to design optimization as well as error estimation and grid adaptation efforts.
Application of Adjoint Methodology to Supersonic Aircraft Design Using Reversed Equivalent Areas
NASA Technical Reports Server (NTRS)
Rallabhandi, Sriram K.
2013-01-01
This paper presents an approach to shape an aircraft to equivalent area based objectives using the discrete adjoint approach. Equivalent areas can be obtained either using reversed augmented Burgers equation or direct conversion of off-body pressures into equivalent area. Formal coupling with CFD allows computation of sensitivities of equivalent area objectives with respect to aircraft shape parameters. The exactness of the adjoint sensitivities is verified against derivatives obtained using the complex step approach. This methodology has the benefit of using designer-friendly equivalent areas in the shape design of low-boom aircraft. Shape optimization results with equivalent area cost functionals are discussed and further refined using ground loudness based objectives.
NASA Astrophysics Data System (ADS)
Zhao, Xiao-Feng; Huang, Si-Xun; Du, Hua-Dong
2011-02-01
This paper puts forward possibilities of refractive index profile retrieval using field measurements at an array of radio receivers in terms of variational adjoint approach. The derivation of the adjoint model begins with the parabolic wave equation for a smooth, perfectly conducting surface and horizontal polarization conditions. To deal with the ill-posed difficulties of the inversion, the regularization ideas are introduced into the establishment of the cost function. Based on steepest descent iterations, the optimal value of refractivity could be retrieved quickly at each point over height. Numerical experiments demonstrate that the method works well for low-distance signals, while it is not accurate enough for long-distance propagations. Through curve fitting processing, however, giving a good initial refractivity profile could generally improve the inversions.
Self-adjoint elliptic operators with boundary conditions on not closed hypersurfaces
NASA Astrophysics Data System (ADS)
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.
Realization of Center Symmetry in Two Adjoint Flavor Large-N Yang-Mills
Catterall, Simon; Galvez, Richard; Unsal, Mithat; /SLAC /Stanford U., Phys. Dept.
2010-08-26
We report on the results of numerical simulations of SU(N) lattice Yang Mills with two flavors of (light) Wilson fermion in the adjoint representation. We analytically and numerically address the question of center symmetry realization on lattices with {Lambda} sites in each direction in the large-N limit. We show, by a weak coupling calculation that, for massless fermions, center symmetry realization is independent of {Lambda}, and is unbroken. Then, we extend our result by conducting simulations at non zero mass and finite gauge coupling. Our results indicate that center symmetry is intact for a range of fermion mass in the vicinity of the critical line on lattices of volume 2{sup 4}. This observation makes it possible to compute infinite volume physical observables using small volume simulations in the limit N {yields} {infinity}, with possible applications to the determination of the conformal window in gauge theories with adjoint fermions.
Adjoint equations and analysis of complex systems: Application to virus infection modelling
NASA Astrophysics Data System (ADS)
Marchuk, G. I.; Shutyaev, V.; Bocharov, G.
2005-12-01
Recent development of applied mathematics is characterized by ever increasing attempts to apply the modelling and computational approaches across various areas of the life sciences. The need for a rigorous analysis of the complex system dynamics in immunology has been recognized since more than three decades ago. The aim of the present paper is to draw attention to the method of adjoint equations. The methodology enables to obtain information about physical processes and examine the sensitivity of complex dynamical systems. This provides a basis for a better understanding of the causal relationships between the immune system's performance and its parameters and helps to improve the experimental design in the solution of applied problems. We show how the adjoint equations can be used to explain the changes in hepatitis B virus infection dynamics between individual patients.
NASA Astrophysics Data System (ADS)
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.
Seismic wave-speed structure beneath the metropolitan area of Japan based on adjoint tomography
NASA Astrophysics Data System (ADS)
Miyoshi, T.; Obayashi, M.; Tono, Y.; Tsuboi, S.
2015-12-01
We have obtained a three-dimensional (3D) model of seismic wave-speed structure beneath the metropolitan area of Japan. We applied the spectral-element method (e.g. Komatitsch and Tromp 1999) and adjoint method (Liu and Tromp 2006) to the broadband seismograms in order to infer the 3D model. We used the travel-time tomography result (Matsubara and Obara 2011) as an initial 3D model and used broadband waveforms recorded at the NIED F-net stations. We selected 147 earthquakes with magnitude of larger than 4.5 from the F-net earthquake catalog and used their bandpass filtered seismograms between 5 and 20 second with a high S/N ratio. The 3D model used for the forward and adjoint simulations is represented as a region of approximately 500 by 450 km in horizontal and 120 km in depth. Minimum period of theoretical waveforms was 4.35 second. For the adjoint inversion, we picked up the windows of the body waves from the observed and theoretical seismograms. We used SPECFEM3D_Cartesian code (e.g. Peter et al. 2011) for the forward and adjoint simulations, and their simulations were implemented by K-computer in RIKEN. Each iteration required about 0.1 million CPU hours at least. The model parameters of Vp and Vs were updated by using the steepest descent method. We obtained the fourth iterative model (M04), which reproduced observed waveforms better than the initial model. The shear wave-speed of M04 was significantly smaller than the initial model at any depth. The model of compressional wave-speed was not improved by inversion because of small alpha kernel values. Acknowledgements: This research was partly supported by MEXT Strategic Program for Innovative Research. We thank to the NIED for providing seismological data.
Towards Adjoint Finite Source Inversion: Application to the 2011 M9 Tohoku Earthquake
NASA Astrophysics Data System (ADS)
Somala, S.; Galvez, P.; Inbal, A.; Ampuero, J. P.; Lapusta, N.
2011-12-01
The recent 2011 M9 Tohoku, Japan, earthquake was recorded by thousands of sensors at near-fault distance, including broad band, strong motion and continuous GPS sensors. This event provides a unique opportunity to image the earthquake rupture process with high resolution. In order to enable the exploitation of the immense dataset available, orders of magnitude larger than in previous earthquakes, we are developing a scalable source inversion procedure based on time-reversal adjoint inversion. We adopt the linear least squares formulation of the source inversion problem, whose basic unknown is the spatio-temporal distribution of slip rate. We formulate an iterative conjugate gradient procedure to minimize the L2 norm of ground velocity residuals between data and synthetics. Each iteration involves one time-reversal (adjoint) and one forward simulation. Exploiting the time-reversal symmetry and the reciprocity principle of elastodynamics, the adjoint is computed by a wave propagation simulation in which time-reversed seismogram residuals are imposed as point forces at the stations simulated. The resulting fault tractions on a locked fault are the adjoint fields, related to the gradient of the misfit function with respect to the model. The simulations are performed with a recent extension of the SPECFEM3D spectral element code to dynamic and kinematic finite sources on unstructured meshes (Galvez et al, session S24 of this meeting). The non-planar geometry of the megathrust fault is accounted for in the spectral element mesh (generated with CUBIT). The subsurface structure is incorporated, on a coarse scale, using regional 3D velocity models, e.g. from the Japan Seismic Hazard Information Station (J-SHIS) website. We will report on the results of our initial efforts, focused on exploiting the continuous 1 Hz GPS signals recorded in Japan to understand the low frequency aspects of the rupture process of the 2011 Tohoku earthquake.
Some results on the dynamics and transition probabilities for non self-adjoint hamiltonians
Bagarello, F.
2015-05-15
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.
Nonlinear self adjointness, conservation laws and exact solutions of ill-posed Boussinesq equation
NASA Astrophysics Data System (ADS)
Yaşar, Emrullah; San, Sait; Özkan, Yeşim Sağlam
2016-01-01
In this work, we consider the ill-posed Boussinesq equation which arises in shallow water waves and non-linear lattices. We prove that the ill-posed Boussinesq equation is nonlinearly self-adjoint. Using this property and Lie point symmetries, we construct conservation laws for the underlying equation. In addition, the generalized solitonary, periodic and compact-like solutions are constructed by the exp-function method.
Inequivalence of unitarity and self-adjointness: An example in quantum cosmology
Lemos, N.A. )
1990-02-15
An example of a quantum cosmological model is presented whose dynamics is unitary although the time-dependent Hamiltonian operator fails to be self-adjoint (because it is not defined) for a particular value of {ital t}. The model is shown to be singular, and this disproves a conjecture put forward by Gotay and Demaret to the effect that unitary quantum dynamics in a slow-time'' gauge is always nonsingular.
NASA Technical Reports Server (NTRS)
Andrews, A.
2002-01-01
A detailed mechanistic understanding of the sources and sinks of CO2 will be required to reliably predict future COS 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 long-term surface monitoring stations with data from intensive field campaigns and with proposed future satellite observations. A major advantage of the adjoint approach is that meteorological and surface data, as well as data for other atmospheric constituents and pollutants can be efficiently included in addition to observations of CO2 mixing ratios. This presentation will provide an overview of potentially useful datasets for carbon cycle research in general with an emphasis on planning for the North American Carbon Project. Areas of overlap with ongoing and proposed work on air quality/air pollution issues will be highlighted.
Source attribution of particulate matter pollution over North China with the adjoint method
NASA Astrophysics Data System (ADS)
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.
Towards efficient backward-in-time adjoint computations using data compression techniques
Cyr, E. C.; Shadid, J. N.; Wildey, T.
2014-12-16
In the context of a posteriori error estimation for nonlinear time-dependent partial differential equations, the state-of-the-practice is to use adjoint approaches which require the solution of a backward-in-time problem defined by a linearization of the forward problem. One of the major obstacles in the practical application of these approaches, we found, is the need to store, or recompute, the forward solution to define the adjoint problem and to evaluate the error representation. Our study considers the use of data compression techniques to approximate forward solutions employed in the backward-in-time integration. The development derives an error representation that accounts for themore » difference between the standard-approach and the compressed approximation of the forward solution. This representation is algorithmically similar to the standard representation and only requires the computation of the quantity of interest for the forward solution and the data-compressed reconstructed solution (i.e. scalar quantities that can be evaluated as the forward problem is integrated). This approach is then compared with existing techniques, such as checkpointing and time-averaged adjoints. Lastly, we provide numerical results indicating the potential efficiency of our approach on a transient diffusion–reaction equation and on the Navier–Stokes equations. These results demonstrate memory compression ratios up to 450×450× while maintaining reasonable accuracy in the error-estimates.« less
Towards efficient backward-in-time adjoint computations using data compression techniques
Cyr, E. C.; Shadid, J. N.; Wildey, T.
2014-12-16
In the context of a posteriori error estimation for nonlinear time-dependent partial differential equations, the state-of-the-practice is to use adjoint approaches which require the solution of a backward-in-time problem defined by a linearization of the forward problem. One of the major obstacles in the practical application of these approaches, we found, is the need to store, or recompute, the forward solution to define the adjoint problem and to evaluate the error representation. Our study considers the use of data compression techniques to approximate forward solutions employed in the backward-in-time integration. The development derives an error representation that accounts for the difference between the standard-approach and the compressed approximation of the forward solution. This representation is algorithmically similar to the standard representation and only requires the computation of the quantity of interest for the forward solution and the data-compressed reconstructed solution (i.e. scalar quantities that can be evaluated as the forward problem is integrated). This approach is then compared with existing techniques, such as checkpointing and time-averaged adjoints. Lastly, we provide numerical results indicating the potential efficiency of our approach on a transient diffusion–reaction equation and on the Navier–Stokes equations. These results demonstrate memory compression ratios up to 450×450× while maintaining reasonable accuracy in the error-estimates.
Adjoint-based deviational Monte Carlo methods for phonon transport calculations
NASA Astrophysics Data System (ADS)
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.
Self-adjointness of the Fourier expansion of quantized interaction field Lagrangians
Paneitz, S. M.; Segal, I. E.
1983-01-01
Regularity properties significantly stronger than were previously known are developed for four-dimensional non-linear conformally invariant quantized fields. The Fourier coefficients of the interaction Lagrangian in the interaction representation—i.e., evaluated after substitution of the associated quantized free field—is a densely defined operator on the associated free field Hilbert space K. These Fourier coefficients are with respect to a natural basis in the universal cosmos ˜M, to which such fields canonically and maximally extend from Minkowski space-time M0, which is covariantly a submanifold of ˜M. However, conformally invariant free fields over M0 and ˜M are canonically identifiable. The kth Fourier coefficient of the interaction Lagrangian has domain inclusive of all vectors in K to which arbitrary powers of the free hamiltonian in ˜M are applicable. Its adjoint in the rigorous Hilbert space sense is a-k in the case of a hermitian Lagrangian. In particular (k = 0) the leading term in the perturbative expansion of the S-matrix for a conformally invariant quantized field in M0 is a self-adjoint operator. Thus, e.g., if ϕ(x) denotes the free massless neutral scalar field in M0, then ∫M0:ϕ(x)4:d4x is a self-adjoint operator. No coupling constant renormalization is involved here. PMID:16593346
Neural network training by integration of adjoint systems of equations forward in time
NASA Technical Reports Server (NTRS)
Toomarian, Nikzad (Inventor); Barhen, Jacob (Inventor)
1992-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. The trajectories computed using our new method are much closer to the target trajectories than was reported in previous studies.
Neural Network Training by Integration of Adjoint Systems of Equations Forward in Time
NASA Technical Reports Server (NTRS)
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.
Strong Adjoint Sensitivities in Tropical Eddy-Permitting Variational Data Assimilation
NASA Astrophysics Data System (ADS)
Cornuelle, B.; Hoteit, I.; Koehl, A.; Stammer, D.
2007-05-01
A variational data assimilation system has been implemented for the tropical Pacific Ocean for an eddy-permitting regional implementation of the MIT general circulation model (MITgcm). The model uses realistic topography with parameterizations for the surface boundary layer (KPP) and open boundaries at the south and north, as well as in the Indonesian throughflow. The adjoint method is used to adjust the model to observations in the tropical Pacific region using control parameters which include initial temperature and salinity, temperature, salinity and horizontal velocities at the open boundaries, and twice-daily surface fluxes of momentum, heat and freshwater. The model is constrained with most of the available datasets in the tropical Pacific, including climatologies, TAO, ARGO, XBT, and satellite SST and SSH data. The forward model runs exhibit strongly growing flow instabilities in the regions of high kinetic energy and low planetary potential vorticity gradient. The growth of these perturbations is limited by nonlinearities once they reach finite size, meaning that the high linear growth rates do not apply for long time periods. This poses a technical problem for adjoint-based assimilation, which depends on the linearized sensitivities to adjust the controls. Relative to the forward model runs, increased viscosity and diffusivity terms are used in the adjoint model runs to avoid large sensitivities related to the flow instabilities present in the high-resolution model. This talk will discuss some of the technical aspects and show results for 1 year assimilation period.
A modified Monte Carlo 'local importance function transform' method
Keady, K. P.; Larsen, E. W.
2013-07-01
The Local Importance Function Transform (LIFT) method uses an approximation of the contribution transport problem to bias a forward Monte-Carlo (MC) source-detector simulation [1-3]. Local (cell-based) biasing parameters are calculated from an inexpensive deterministic adjoint solution and used to modify the physics of the forward transport simulation. In this research, we have developed a new expression for the LIFT biasing parameter, which depends on a cell-average adjoint current to scalar flux (J{sup *}/{phi}{sup *}) ratio. This biasing parameter differs significantly from the original expression, which uses adjoint cell-edge scalar fluxes to construct a finite difference estimate of the flux derivative; the resulting biasing parameters exhibit spikes in magnitude at material discontinuities, causing the original LIFT method to lose efficiency in problems with high spatial heterogeneity. The new J{sup *}/{phi}{sup *} expression, while more expensive to obtain, generates biasing parameters that vary smoothly across the spatial domain. The result is an improvement in simulation efficiency. A representative test problem has been developed and analyzed to demonstrate the advantage of the updated biasing parameter expression with regards to solution figure of merit (FOM). For reference, the two variants of the LIFT method are compared to a similar variance reduction method developed by Depinay [4, 5], as well as MC with deterministic adjoint weight windows (WW). (authors)
NASA Astrophysics Data System (ADS)
Cacuci, Dan G.
2015-03-01
This work presents an illustrative application of the second-order adjoint sensitivity analysis methodology (2nd-ASAM) to a paradigm neutron diffusion problem, which is sufficiently simple to admit an exact solution, thereby making transparent the underlying mathematical derivations. The general theory underlying 2nd-ASAM indicates that, for a physical system comprising Nα parameters, the computation of all of the first- and second-order response sensitivities requires (per response) at most (2Nα + 1) "large-scale" computations using the first-level and, respectively, second-level adjoint sensitivity systems (1st-LASS and 2nd-LASS). Very importantly, however, the illustrative application presented in this work shows that the actual number of adjoint computations needed for computing all of the first- and second-order response sensitivities may be significantly less than (2Nα + 1) per response. For this illustrative problem, four "large-scale" adjoint computations sufficed for the complete and exact computations of all 4 first- and 10 distinct second-order derivatives. Furthermore, the construction and solution of the 2nd-LASS requires very little additional effort beyond the construction of the adjoint sensitivity system needed for computing the first-order sensitivities. Very significantly, only the sources on the right-sides of the diffusion (differential) operator needed to be modified; the left-side of the differential equations (and hence the "solver" in large-scale practical applications) remained unchanged. All of the first-order relative response sensitivities to the model parameters have significantly large values, of order unity. Also importantly, most of the second-order relative sensitivities are just as large, and some even up to twice as large as the first-order sensitivities. In the illustrative example presented in this work, the second-order sensitivities contribute little to the response variances and covariances. However, they have the
NASA Astrophysics Data System (ADS)
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.
NASA Astrophysics Data System (ADS)
Guerrette, J. J.; Henze, D. K.
2015-02-01
Here we present the online meteorology and chemistry adjoint and tangent linear model, WRFPLUS-Chem, which incorporates modules to treat boundary layer mixing, emission, aging, dry deposition, and advection of black carbon aerosol. We also develop land surface and surface layer adjoints to account for coupling between radiation and vertical mixing. Model performance is verified against finite difference derivative approximations. A second order checkpointing scheme is created to reduce computational costs and enable simulations longer than six hours. The adjoint is coupled to WRFDA-Chem, in order to conduct a sensitivity study of anthropogenic and biomass burning sources throughout California during the 2008 Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) field campaign. A cost function weighting scheme was devised to increase adjoint sensitivity robustness in future inverse modeling studies. Results of the sensitivity study show that, for this domain and time period, anthropogenic emissions are over predicted, while wildfire emissions are under predicted. We consider the diurnal variation in emission sensitivities to determine at what time sources should be scaled up or down. Also, adjoint sensitivities for two choices of land surface model indicate that emission inversion results would be sensitive to forward model configuration. The tools described here are the first step in conducting four-dimensional variational data assimilation in a coupled meteorology-chemistry model, which will potentially provide new constraints on aerosol precursor emissions and their distributions. Such analyses will be invaluable to assessments of particulate matter health and climate impacts.
NASA Technical Reports Server (NTRS)
Reuther, James; Jameson, Antony; Alonso, Juan Jose; Rimlinger, Mark J.; Saunders, David
1997-01-01
An aerodynamic shape optimization method that treats the design of complex aircraft configurations subject to high fidelity computational fluid dynamics (CFD), geometric constraints and multiple design points is described. The design process will be greatly accelerated through the use of both control theory and distributed memory computer architectures. Control theory is employed to derive the adjoint differential equations whose solution allows for the evaluation of design gradient information at a fraction of the computational cost required by previous design methods. The resulting problem is implemented on parallel distributed memory architectures using a domain decomposition approach, an optimized communication schedule, and the MPI (Message Passing Interface) standard for portability and efficiency. The final result achieves very rapid aerodynamic design based on a higher order CFD method. In order to facilitate the integration of these high fidelity CFD approaches into future multi-disciplinary optimization (NW) applications, new methods must be developed which are capable of simultaneously addressing complex geometries, multiple objective functions, and geometric design constraints. In our earlier studies, we coupled the adjoint based design formulations with unconstrained optimization algorithms and showed that the approach was effective for the aerodynamic design of airfoils, wings, wing-bodies, and complex aircraft configurations. In many of the results presented in these earlier works, geometric constraints were satisfied either by a projection into feasible space or by posing the design space parameterization such that it automatically satisfied constraints. Furthermore, with the exception of reference 9 where the second author initially explored the use of multipoint design in conjunction with adjoint formulations, our earlier works have focused on single point design efforts. Here we demonstrate that the same methodology may be extended to treat
Mapping Emissions that Contribute to Air Pollution Using Adjoint Sensitivity Analysis
NASA Astrophysics Data System (ADS)
Bastien, L. A. J.; Mcdonald, B. C.; Brown, N. J.; Harley, R.
2014-12-01
The adjoint of the Community Multiscale Air Quality model (CMAQ) is used to map emissions that contribute to air pollution at receptors of interest. Adjoint tools provide an efficient way to calculate the sensitivity of a model response to a large number of model inputs, a task that would require thousands of simulations using a more traditional forward sensitivity approach. Initial applications of this technique, demonstrated here, are to benzene and directly-emitted diesel particulate matter, for which atmospheric reactions are neglected. Emissions of these pollutants are strongly influenced by light-duty gasoline vehicles and heavy-duty diesel trucks, respectively. We study air quality responses in three receptor areas where populations have been identified as especially susceptible to, and adversely affected by air pollution. Population-weighted air basin-wide responses for each pollutant are also evaluated for the entire San Francisco Bay area. High-resolution (1 km horizontal grid) emission inventories have been developed for on-road motor vehicle emission sources, based on observed traffic count data. Emission estimates represent diurnal, day of week, and seasonal variations of on-road vehicle activity, with separate descriptions for gasoline and diesel sources. Emissions that contribute to air pollution at each receptor have been mapped in space and time using the adjoint method. Effects on air quality of both relative (multiplicative) and absolute (additive) perturbations to underlying emission inventories are analyzed. The contributions of local versus upwind sources to air quality in each receptor area are quantified, and weekday/weekend and seasonal variations in the influence of emissions from upwind areas are investigated. The contribution of local sources to the total air pollution burden within the receptor areas increases from about 40% in the summer to about 50% in the winter due to increased atmospheric stagnation. The effectiveness of control
Imaging the slab beneath central Chile using the Spectral Elements Method and adjoint techniques
NASA Astrophysics Data System (ADS)
Mercerat, E. D.; Nolet, G.; Marot, M.; Deshayes, P.; Monfret, T.
2010-12-01
This work focuses on imaging the subducting slab beneath Central Chile using novel inversion techniques based on the adjoint method and accurate wave propagation simulations using the Spectral Elements Method. The study area comprises the flat slab portion of the Nazca plate between 29 S and 34 S subducting beneath South America. We will use a database of regional seismicity consisting of both crustal and deep slab earthquakes with magnitude 3 < Mw < 6 recorded by different temporary and permanent seismological networks. Our main goal is to determine both the kinematics and the geometry of the subducting slab in order to help the geodynamical interpretation of such particular active margin. The Spectral Elements Method (SPECFEM3D code) is used to generate the synthetic seismograms and it will be applied for the iterative minimization based on adjoint techniques. The numerical mesh is 600 km x 600 km in horizontal coordinates and 220 km depth. As a first step, we are faced to well-known issues concerning mesh generation (resolution, quality, absorbing boundary conditions). In particular, we must evaluate the influence of free surface topography, as well as the MOHO and other geological interfaces in the synthetic seismograms. The initial velocity model from a previous travel-time tomography study, is linearly interpolated to the Gauss-Lobatto-Legendre grid. The comparison between the first forward simulations (up to 4 seconds minimum period) validate the initial velocity model of the study area, although many features not reproduced by the initial model have already been identified. Next step will concentrate in the comparison between finite-frequency kernels calculated by travel-time methods with ones based on adjoint methods, in order to highlight advantages and disadvantages in terms of resolution, accuracy, but also computational cost.
Adjoint Sensitivity Computations for an Embedded-Boundary Cartesian Mesh Method and CAD Geometry
NASA Technical Reports Server (NTRS)
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
NASA Astrophysics Data System (ADS)
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.
Diffusion Acceleration Schemes for Self-Adjoint Angular Flux Formulation with a Void Treatment
Yaqi Wang; Hongbin Zhang; Richard C. Martineau
2014-02-01
A Galerkin weak form for the monoenergetic neutron transport equation with a continuous finite element method and discrete ordinate method is developed based on self-adjoint angular flux formulation. This weak form is modified for treating void regions. A consistent diffusion scheme is developed with projection. Correction terms of the diffusion scheme are derived to reproduce the transport scalar flux. A source iteration that decouples the solution of all directions with both linear and nonlinear diffusion accelerations is developed and demonstrated. One-dimensional Fourier analysis is conducted to demonstrate the stability of the linear and nonlinear diffusion accelerations. Numerical results of these schemes are presented.
Pricing of American style options with an adjoint process correction method
NASA Astrophysics Data System (ADS)
Jaekel, Uwe
2005-07-01
Pricing of American options is a more complicated problem than pricing of European options. In this work a formula is derived that allows the computation of the early exercise premium, i.e. the price difference between these two option types in terms of an adjoint process evolving in the reversed time direction of the original process determining the evolution of the European price. We show how this equation can be utilised to improve option price estimates from numerical schemes like finite difference or Monte Carlo methods.
Adjoint transport calculations for sensitivity analysis of the Hiroshima air-over-ground environment
Broadhead, B.L.; Cacuci, D.G.; Pace, J.V. III
1984-01-01
A major effort within the US Dose Reassessment Program is aimed at recalculating the transport of initial nuclear radiation in an air-over-ground environment. This paper is the first report of results from adjoint calculations in the Hiroshima air-over-ground environment. The calculations use a Hiroshima/Nagasaki multi-element ground, ENDF/B-V nuclear data, one-dimensional ANISN flux weighting for neutron and gamma cross sections, a source obtained by two-dimensional hydrodynamic and three-dimensional transport calculations, and best-estimate atmospheric conditions from Japanese sources. 7 references, 2 figures.
A forward operator and its adjoint for GPS slant total delays
NASA Astrophysics Data System (ADS)
Zus, Florian; Dick, Galina; Heise, Stefan; Wickert, Jens
2015-05-01
In a recent study we developed a fast and accurate algorithm to compute Global Positioning System (GPS) Slant Total Delay (STDs) utilizing numerical weather model data. Having developed a forward operator we construct in this study the tangent linear (adjoint) operator by application of the chain rule of differential calculus in forward (reverse) mode. Armed with these operators we show in a simulation study the potential benefit of GPS STDs in inverse modeling. We conclude that the developed operators are tailored for three (four)-dimensional variational data assimilation and/or travel time tomography.
Infrared fixed point in SU(2) gauge theory with adjoint fermions
DeGrand, Thomas; Shamir, Yigal; Svetitsky, Benjamin
2011-04-01
We apply Schroedinger-functional techniques to the SU(2) lattice gauge theory with N{sub 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{sub *} where 1/g{sub *}{sup 2}=0.20(4)(3). At the fixed point, we find a mass anomalous dimension {gamma}{sub m}(g{sub *})=0.31(6).
On the Neutron Kinetics and Control of Accelerator-Driven Systems
Cacuci, D.G.
2004-09-15
This work addresses fundamental aspects of the time- and space-dependent behavior of an Accelerator-Driven Subcritical Core System (ADS) and presents a paradigm ADS neutron kinetics model that is solved exactly. Thus, this paradigm model can serve for benchmarking two- and/or three-dimensional computational tools. Furthermore, this work also proposes a global optimal control theory framework for the operation and control of an ADS. This framework encompasses conceptually the time- and space-dependent behavior of the ADS coupled neutron kinetics/thermal-hydraulic balance equations and aims at the optimal control of ADS operational objectives, which would include minimization of local flux disturbances, load and source following, etc. Importantly, this new conceptual framework makes no use of a 'fictitious ADS steady state' and yields the correct and complete (i.e., including sources) adjoint equations, without leaving any room for ambiguities. Thus, this new conceptual framework provides a natural basis for developing new computational methods and corresponding verification experiments specifically tailored for the control and operation of ADS.
Hilbert-Schmidt Inner Product for an Adjoint Representation of the Quantum Algebra U⌣Q(SU2)
NASA Astrophysics Data System (ADS)
Fakhri, Hossein; Nouraddini, Mojtaba
2015-10-01
The Jordan-Schwinger realization of quantum algebra U⌣q(su2) is used to construct the irreducible submodule Tl of the adjoint representation in two different bases. The two bases are known as types of irreducible tensor operators of rank l which are related to each other by the involution map. The bases of the submodules are equipped with q-analogues of the Hilbert-Schmidt inner product and it is also shown that the adjoint representation corresponding to one of those submodules is a *-representation.
Nievaart, V. A.; Legrady, D.; Moss, R. L.; Kloosterman, J. L.; Hagen, T. H. J. J. van der; Dam, H. van
2007-04-15
This paper deals with the application of the adjoint transport theory in order to optimize Monte Carlo based radiotherapy treatment planning. The technique is applied to Boron Neutron Capture Therapy where most often mixed beams of neutrons and gammas are involved. In normal forward Monte Carlo simulations the particles start at a source and lose energy as they travel towards the region of interest, i.e., the designated point of detection. Conversely, with adjoint Monte Carlo simulations, the so-called adjoint particles start at the region of interest and gain energy as they travel towards the source where they are detected. In this respect, the particles travel backwards and the real source and real detector become the adjoint detector and adjoint source, respectively. At the adjoint detector, an adjoint function is obtained with which numerically the same result, e.g., dose or flux in the tumor, can be derived as with forward Monte Carlo. In many cases, the adjoint method is more efficient and by that is much quicker when, for example, the response in the tumor or organ at risk for many locations and orientations of the treatment beam around the patient is required. However, a problem occurs when the treatment beam is mono-directional as the probability of detecting adjoint Monte Carlo particles traversing the beam exit (detector plane in adjoint mode) in the negative direction of the incident beam is zero. This problem is addressed here and solved first with the use of next event estimators and second with the application of a Legendre expansion technique of the angular adjoint function. In the first approach, adjoint particles are tracked deterministically through a tube to a (adjoint) point detector far away from the geometric model. The adjoint particles will traverse the disk shaped entrance of this tube (the beam exit in the actual geometry) perpendicularly. This method is slow whenever many events are involved that are not contributing to the point
Sensitivity analysis of numerically-simulated convective storms using direct and adjoint methods
Park, S.K.; Droegemeier, K.K.; Bischof, C.; Knauff, T.
1994-06-01
The goal of this project is to evaluate the sensitivity of numerically modeled convective storms to control parameters such as the initial conditions, boundary conditions, environment, and various physical and computational parameters. In other words, the authors seek the gradient of the solution vector with respect to specified parameters. One can use two approaches to accomplish this task. In the first or so-called brute force method, one uses a fully nonlinear model to generate a control forecast starting from a specified initial state. Then, a number of other forecasts are made in which chosen parameters (e.g., initial conditions) are systematically varied. The obvious drawback is that a large number of full model predictions are needed to examine the effects of only a single parameter. The authors describe herein an alternative, essentially automated method (ADIFOR, or Automatic DIfferentiation of FORtran) for obtaining the solution gradient that bypasses the adjoint altogether yet provides even more information about the gradient. (ADIFOR, like the adjoint technique, is constrained by the linearity assumption.) Applied to a 1-D moist cloud model, the authors assess the utility of ADIFOR relative to the brute force approach and evaluate the validity of the tangent linear approximation in the context of deep convection.
Self-Adjoint Angular Flux Equation for Coupled Electron-Photon Transport
Liscum-Powell, J.L.; Lorence, L.J. Jr.; Morel, J.E.; Prinja, A.K.
1999-07-08
Recently, Morel and McGhee described an alternate second-order form of the transport equation called the self adjoint angular flux (SAAF) equation that has the angular flux as its unknown. The SAAF formulation has all the advantages of the traditional even- and odd-parity self-adjoint equations, with the added advantages that it yields the full angular flux when it is numerically solved, it is significantly easier to implement reflective and reflective-like boundary conditions, and in the appropriate form it can be solved in void regions. The SAAF equation has the disadvantage that the angular domain is the full unit sphere and, like the even- and odd- parity form, S{sub n} source iteration cannot be implemented using the standard sweeping algorithm. Also, problems arise in pure scattering media. Morel and McGhee demonstrated the efficacy of the SAAF formulation for neutral particle transport. Here we apply the SAAF formulation to coupled electron-photon transport problems using multigroup cross-sections from the CEPXS code and S{sub n} discretization.
Improving the Fit of a Land-Surface Model to Data Using its Adjoint
NASA Astrophysics Data System (ADS)
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.
Ocean acoustic tomography from different receiver geometries using the adjoint method.
Zhao, Xiaofeng; Wang, Dongxiao
2015-12-01
In this paper, an ocean acoustic tomography inversion using the adjoint method in a shallow water environment is presented. The propagation model used is an implicit Crank-Nicolson finite difference parabolic equation solver with a non-local boundary condition. Unlike previous matched-field processing works using the complex pressure fields as the observations, here, the observed signals are the transmission losses. Based on the code tests of the tangent linear model, the adjoint model, and the gradient, the optimization problem is solved by a gradient-based minimization algorithm. The inversions are performed in numerical simulations for two geometries: one in which hydrophones are sparsely distributed in the horizontal direction, and another in which the hydrophones are distributed vertically. The spacing in both cases is well beyond the half-wavelength threshold at which beamforming could be used. To deal with the ill-posedness of the inverse problem, a linear differential regularization operator of the sound-speed profile is used to smooth the inversion results. The L-curve criterion is adopted to select the regularization parameter, and the optimal value can be easily determined at the elbow of the logarithms of the residual norm of the measured-predicted fields and the norm of the penalty function. PMID:26723329