Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “adjoint problem”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Adaptive mesh refinement for steady flows in Nek5000

Adaptive mesh refinement is performed in the framework of the spectral element method augmented by approaches to error estimation and control. The h-refinement technique is used for adapting the mesh, where selected grid elements are split by a quadtree (2D) or octree (3D) structure. Continuity between parent-child elements is enforced by high-order interpolation of the solution across the common faces. Parallel mesh partitioning and grid management respectively, are taken care of by the external libraries ParMETIS and p4est. Two methods are considered for estimating and controlling the error of the solution. The first error estimate is local and based on the spectral properties of the solution on each element. This method gives a local measure of the L-2-norm of the solution over the entire computational domain. The second error estimate uses the dual-weighted residuals method - it is based on and takes into account both the local properties of the solution and the global dependence of the error in the solution via an adjoint problem. The objective of this second approach is to optimize the computation of a given functional of physical interest. The simulations are performed by using the code Nek5000 and three steady-state test cases are studied: a two-dimensional lid-driven cavity at Re = 7, 500, a two-dimensional flow past a cylinder at Re = 40, and a three-dimensional lid-driven cavity at Re = 2,000 with a moving lid tilted by an angle of 30 degrees. Finally, the efficiency of both error estimators is compared in terms of refinement patterns and accuracy on the functional of interest. In the case of the adjoint error estimators, the trend on the error of the functional is shown to be correctly represented up to a multiplicative constant.

42 ENGINEERING↗
DOE OSTI2025

OpenSn: A massively parallel, open-source simulation environment for discrete ordinates radiation transport

OpenSn is an open-source, massively parallel deterministic radiation transport code for solving the discrete-ordinates ( S N ) form of the Boltzmann transport equation on unstructured, arbitrary polyhedral meshes. It supports high-fidelity simulations involving steady-state, eigenvalue, and adjoint problems for neutral particles (e.g., neutrons, photons, multi-particles), using the multigroup approximation in energy. OpenSn combines angular discretization via discrete ordinates with a discontinuous Galerkin finite element method (DGFEM) in space, enabling accurate resolution of transport physics on arbitrary polyhedral cells, included locally refined spatial grids. It includes multiple angular quadrature types, including locally refined angular quadratures. Written in modern C++ with a Python API, OpenSn runs efficiently on platforms ranging from laptops to supercomputers. The transport sweep algorithm is implemented using a task-based, directed-acyclic-graph (DAG) approach for each angle and supports asynchronous parallelism across thousands of MPI ranks. Group-set aggregation improves compute intensity, and synthetic acceleration techniques (e.g., diffusion synthetic acceleration, second-moment method) enhance solver convergence. OpenSn has been verified on reactor physics problems and demonstrated excellent weak and strong scaling performance on more than 32,768 processes, making it a versatile and robust platform for large-scale transport simulations in complex geometries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Longitudinal sound waves in a collisionless, quasineutral plasma

The time evolution of slow sound waves in a homogeneous, collisionless and quasineutral plasma, in particular their Landau damping, is investigated using the kinetic-magnetohydrodynamics formulation of Ramos ( J. Plasma Phys. vol. 81, 2015 p. 905810325; vol. 82, 2016 p. 905820607). In this approach, the electric field is eliminated from a closed, hybrid fluid-kinetic system that ensures automatically the fulfilment of the charge neutrality condition. Considering the time dependence of a spatial-Fourier-mode linear perturbation with wavevector parallel to the equilibrium magnetic field, this can be cast as a second-order self-adjoint problem with a continuum spectrum of real and positive squared frequencies. Therefore, a conventional resolution of the identity with a continuum basis of singular normal modes is guaranteed, which simplifies significantly a Van Kampen-like treatment of the Landau damping. The explicit form of such singular normal modes is obtained, along with their orthogonality relations. These are used to derive the damped time evolution of the fluid moments of solutions of initial-value problems, for the most general kinds of initial conditions. The non-zero parallel electric field is not used explicitly in this analysis, but it is calculated from any given solution after the later has been obtained.

Physics↗
DOE OSTI2018

Kinetic-magnetohydrodynamic description of longitudinal sound waves in a collisionless, quasineutral plasma

The time evolution of slow sound waves in a homogeneous, collisionless and quasineutral plasma, in particular their Landau damping, is investigated using the kinetic magnetohydrodynamics formulation of J.J. Ramos (J . Plasma Phys. vol. 81, 2015 p. 905810325; vol. 82, 2016 p. 905820607). In this approach, the electric field is eliminated from a closed, hybrid fluid-kinetic system that ensures automatically the fulfillment of the charge neutrality condition. Considering the time dependence of a spatial-Fourier-mode linear perturbation with wavevector parallel to the equilibrium magnetic field, this can be cast as a second-order self-adjoint problem with a continuum spectrum of real and positive squared frequencies. Therefore, a conventional resolution of the identity with a continuum basis of singular normal modes is guaranteed, which simplifies significantly a Van Kampen-like treatment of the Landau damping. The explicit form of such singular normal modes is obtained, along with their orthogonality relations. These are used to derive the damped time evolution of the fluid moments of solutions of initial-value problems, for the most general kinds of initial conditions. The non-zero parallel electric field is not used explicitly in this analysis, but it is calculated from any given solution after the later has been obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1984

The solution of non-linear hyperbolic equation systems by the finite element method

A finite-element method for the solution of nonlinear hyperbolic systems of equations, such as those encountered in non-self-adjoint problems of transient phenomena in convection-diffusion or in the mixed representation of wave problems, is developed and demonstrated. The problem is rewritten in moving coordinates and reinterpolated to the original mesh by a Taylor expansion prior to a standard Galerkin spatial discretization, and it is shown that this procedure is equivalent to the time-discretization approach of Donea (1984). Numerical results for sample problems are presented graphically, including such shallow-water problems as the breaking of a dam, the shoaling of a wave, and the outflow of a river; compressible flows such as the isothermal flow in a nozzle and the Riemann shock-tube problem; and the two-dimensional scalar-advection, nonlinear-shallow-water, and Euler equations.

Loehner, R.↗
NASA NTRS1989

Local error estimates for discontinuous solutions of nonlinear hyperbolic equations

Let u(x,t) be the possibly discontinuous entropy solution of a nonlinear scalar conservation law with smooth initial data. Suppose u sub epsilon(x,t) is the solution of an approximate viscosity regularization, where epsilon greater than 0 is the small viscosity amplitude. It is shown that by post-processing the small viscosity approximation u sub epsilon, pointwise values of u and its derivatives can be recovered with an error as close to epsilon as desired. The analysis relies on the adjoint problem of the forward error equation, which in this case amounts to a backward linear transport with discontinuous coefficients. The novelty of this approach is to use a (generalized) E-condition of the forward problem in order to deduce a W(exp 1,infinity) energy estimate for the discontinuous backward transport equation; this, in turn, leads one to an epsilon-uniform estimate on moments of the error u(sub epsilon) - u. This approach does not follow the characteristics and, therefore, applies mutatis mutandis to other approximate solutions such as E-difference schemes.

Tadmor, Eitan↗
NASA NTRS2002

Adjoint-Based, Three-Dimensional Error Prediction and Grid Adaptation

Engineering computational fluid dynamics (CFD) analysis and design applications focus on output functions (e.g., lift, drag). Errors in these output functions are generally unknown and conservatively accurate solutions may be computed. Computable error estimates can offer the possibility to minimize computational work for a prescribed error tolerance. Such an estimate can be computed by solving the flow equations and the linear adjoint problem for the functional of interest. The computational mesh can be modified to minimize the uncertainty of a computed error estimate. This robust mesh-adaptation procedure automatically terminates when the simulation is within a user specified error tolerance. This procedure for estimating and adapting to error in a functional is demonstrated for three-dimensional Euler problems. An adaptive mesh procedure that links to a Computer Aided Design (CAD) surface representation is demonstrated for wing, wing-body, and extruded high lift airfoil configurations. The error estimation and adaptation procedure yielded corrected functions that are as accurate as functions calculated on uniformly refined grids with ten times as many grid points.

Park, Michael A.↗
NASA NTRS2005

Using an Adjoint Approach to Eliminate Mesh Sensitivities in Computational Design

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.

Nielsen, Eric J.↗
NASA NTRS2006

Using an Adjoint Approach to Eliminate Mesh Sensitivities in Computational Design

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.

Nielsen, Eric J.↗
NASA NTRS1972

Adjoint variational methods in nonconservative stability problems.

A general nonself-adjoint eigenvalue problem is examined and it is shown that the commonly employed approximate methods, such as the Galerkin procedure, the method of weighted residuals and the least square technique lack variational descriptions. When used in their previously known forms they do not yield stationary eigenvalues and eigenfunctions. With the help of an adjoint system, however, several analogous variational descriptions may be developed and it is shown in the present study that by properly restating the method of least squares, stationary eigenvalues may be obtained. Several properties of the adjoint eigenvalue problem, known only for a restricted group, are shown to exist for the more general class selected for study.

Prasad, S. N.↗
NASA NTRS1997

Numerical Computation of Sensitivities and the Adjoint Approach

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.

Lewis, Robert Michael↗
DOE OSTI2017

Revisiting Boundary Perturbation Theory for Inhomogeneous Transport Problems

Adjoint-based first-order perturbation theory is applied again to boundary perturbation problems. Rahnema developed a perturbation estimate that gives an accurate first-order approximation of a flux or reaction rate within a radioactive system when the boundary is perturbed. When the response of interest is the flux or leakage current on the boundary, the Roussopoulos perturbation estimate has long been used. The Rahnema and Roussopoulos estimates differ in one term. Our paper shows that the Rahnema and Roussopoulos estimates can be derived consistently, using different responses, from a single variational functional (due to Gheorghiu and Rahnema), resolving any apparent contradiction. In analytic test problems, Rahnema’s estimate and the Roussopoulos estimate produce exact first derivatives of the response of interest when appropriately applied. We also present a realistic, nonanalytic test problem.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1984

A solution of Mx(double dot) + Cx(dot) + Kx = 0 applicable to the design of active dampers

A solution is presented for the equations of motion for the damped linear oscillator, Mx(double dot) + Cx(dot) + Kx = 0. The algorithm solves a transformed set of equations in terms of the modal variables of the undamped system and, at the same time, solves the adjoint equation of the transformed problem. The adjoint solution is normalized to give the inverse of the solution matrix of the transformed problem. The normalized inverse is useful in design for direct computation of sensitivity derivatives of damping ratios with respect to damping rates. The algorithm is programmed to reduce storage requirements by a factor of three-fourths compared to standard complex eigenvalue subroutines. A numerical example is included.

Thurston, G. A.↗
NASA NTRS1992

A theoretical approach for analyzing the restabilization of wakes

A linear stability analysis provides the basis for an analysis of the effects of prescribed changes upon the global stability characteristics of a bluff body wake. The analysis is applied in particular to the problem of placing a small cylinder in the wake of a larger one in order to suppress the vortex street. A prediction of the alteration in the critical Reynolds number of the system as a function of the location of the second cylinder relative to the first is presented and is shown to compare favorably with experimental results. The effect upon shedding frequency is also considered. The approach relies on the solution of an adjoint eigenvalue problem and an inhomogeneous adjoint boundary value problem. These adjoint solutions are shown to be a key component in the analysis of a variety of wake control strategies.

Hill, D. C.↗
NASA NTRS2003

Three-Dimensional Turbulent RANS Adjoint-Based Error Correction

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.

Park, Michael A.↗