Search NASA⌕ Search

SEARCH · Search NASA

Results for “Optimal sizing method”

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 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

High-Temperature MEMS Based Venus Seismometer

Seismology is the method of choice for studying a planet’s interior and assessing current tectonic activity. With a single strategically placed seismometer, or preferably a network of seismometers, a variety of key aspects of a planet can be evaluated. The level and nature of seismicity observed provides a gauge of current geologic, especially tectonic, activity. Where tectonic activity is and is not occurring illuminates what regions and geologic features on the planet are currently active, thus, constraining global geodynamics. Larger distal earthquakes can be used to determine major compositional and structural boundaries within a planet’s interior, providing major constraints on planet formation and interior evolution. Seismometers have been deployed successfully on the Moon and Mars. Only small subsets of the scientific knowledge that can be gained from seismology can be obtained from alternative scientific approaches, such as global gravity solutions or repeat-pass interferometry. By virtue of its similar size to Earth, Venus is likely to be highly seismically active (Fig. 1). Even at Earth-like levels of seismology, returning meaningful data would require observation over a period that is at least 2–3 orders of magnitude longer than the 1–2 h lifetime of previous Soviet landers. Clever insulation could extend the lifetime of standard electronics in a modern lander to perhaps 24 h. While active cooling of electronics for a Venusian lander is plausible, this approach would require considerable technological advancement and will likely come at a high monetary and mass cost. The high density of the Venusian atmosphere at the surface should mean that there is good coupling of seismic energy into the atmosphere, such that seismology using infrasound from a balloon platform is currently being explored, and observations of the upper levels of the atmosphere from orbit has also been proposed [3]. Besides losing all shear wave information with these two methods, calibrating and interpreting seismic information that has been distorted by, and convolved with, atmospheric phenomena will be challenging without existing seismology data collected through surface seismometers (optimally, simultaneous surface seismological observations would be made).

silicon carbide↗
DOE OSTI2022

An automated approach to the alignment of compound refractive lenses

Compound refractive lenses (CRLs) are established X-ray focusing optics, and are used to focus the beam or image the sample in many beamlines at X-ray facilities. While CRLs are quite established, the stack of single lens elements affords a very small numerical aperture because of the thick lens profile, making them far more difficult to align than classical optical lenses that obey the thin-lens approximation. This means that the alignment must be very precise and is highly sensitive to changes to the incident beam, often requiring regular readjustments. Some groups circumvent the full realignment procedure by using engineering controls ( e.g. mounting optics) that sacrifice some of the beam's focusing precision, i.e. spot size, or resolution. While these choices minimize setup time, there are clear disadvantages. This work presents a new automated approach to align CRLs using a simple alignment apparatus that is easy to adapt and install at different types of X-ray experiments or facilities. This approach builds on recent CRL modeling efforts, using an approach based on the Stochastic Nelder–Mead (SNM) simplex method. This method is outlined and its efficacy is demonstrated with numerical simulation that is tested in real experiments conducted at the Advanced Photon Source to confirm its performance with a synchrotron beam. This work provides an opportunity to automate key instrumentation at X-ray facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

F-Hash: Feature-Based Hash Design for Time-Varying Volume Visualization via Multi-Resolution Tesseract Encoding

Interactive time-varying volume visualization is challenging due to its complex spatiotemporal features and sheer size of the dataset. Recent works transform the original discrete time-varying volumetric data into continuous Implicit Neural Representations (INR) to address the issues of compression, rendering, and super-resolution in both spatial and temporal domains. However, training the INR takes a long time to converge, especially when handling large-scale time-varying volumetric datasets. In this work, we proposed F-Hash, a novel feature-based multi-resolution Tesseract encoding architecture to greatly enhance the convergence speed compared with existing input encoding methods for modeling time-varying volumetric data. The proposed design incorporates multi-level collision-free hash functions that map dynamic 4D multi-resolution embedding grids without bucket waste, achieving high encoding capacity with compact encoding parameters. Our encoding method is agnostic to time-varying feature detection methods, making it a unified encoding solution for feature tracking and evolution visualization. Experiments show the F-Hash achieves state-of-the-art convergence speed in training various time-varying volumetric datasets for diverse features. We also proposed an adaptive ray marching algorithm to optimize the sample streaming for faster rendering of the time-varying neural representation.

deep learning↗
NASA NTRS2011

PING Gamma Ray and Neutron Measurements of a Meter-Sized Carbonaceous Asteroid Analog

Determining the elemental composition of carbonaceous (spectral type C) asteroids is still one of the basic problems when studying these objects. The only main source of elemental composition information for asteroids is from their optical, NIR and IR properties, which include their spectral reflectance characteristics, albedo, polarization, and the comparison of optical spectroscopy with meteorite groups corresponding to asteroids of every spectral type. Unfortunately, these sources reflect observations from widely contrasting spatial scales that presently yield a void in the continuum of microscopic and macroscopic evidence, a lack of in situ measurement confirmation, and require deeper sensing techniques to discern the nature of these asteroids. The Probing In situ with Neutrons and Gamma rays (PING) instrument is ideally suited to address this problem because it can be used to determine the bulk elemental composition, H and C content, the average atomic weight and density of the surface and subsurface layers of C-type asteroids, and can provide measurements used to determine the difference between and distinguish between different types of asteroids. We are currently developing the PING instrument that combines gamma ray and neutron detectors with a 14 Me V pulsed neutron generator to determine the in-situ bulk elemental abundances and geochemistry of C-type asteroids with a spatial resolution of 1 m down to depths of tens of cm to 1 m. One aspect of the current work includes experimentally testing and optimizing PING on a known meter-sized Columbia River basalt C-type asteroid analog sample that has a similar composition and the same neutron response as that of a C-type asteroid. An important part of this effort focuses on utilizing timing measurements to isolate gamma rays produced by neutron inelastic scattering, neutron capture and delayed activation processes. Separating the gamma ray spectra by nuclear processes results in higher precision and sensitivity elemental composition measurements. Using gated data acquisition techniques allows for the unambiguous identification of gamma ray lines from different isotopes and nuclear processes, especially in situations when limited detector resolution results in overlapping gamma ray lines that cannot be individually resolved. In this paper, we will present the PING basalt layering experimental data, taken at the test facility at NASA Goddard Space Flight Center using the time tagged event-by-event data analysis technique, compared to our MCNPX computer simulation results for the C-type asteroid and basalt layering simulant models. Comparison of these data will show the advantages, validity, and measurement sensitivity of PING's nuclear interrogation methods to obtain more precise and sensitive in situ bulk elemental composition and density measurements of the subsurface of asteroids.

Bodnarik, J.↗
NASA NTRS2004

Modular, Reconfigurable, and Rapid Response Space Systems: The Remote Sensing Advanced Technology Microsatellite

Modular, Reconfigurable, and Rapid-response (MR(sup 2)) space systems represent a paradigm shift in the way space assets of all sizes are designed, manufactured, integrated, tested, and flown. This paper will describe the MR(sup 2) paradigm in detail, and will include guidelines for its implementation. The Remote Sensing Advanced Technology microsatellite (RSAT) is a proposed flight system test-bed used for developing and implementing principles and best practices for MR(sup 2) spacecraft, and their supporting infrastructure. The initial goal of this test-bed application is to produce a lightweight (approx. 100 kg), production-minded, cost-effective, and scalable remote sensing micro-satellite capable of high performance and broad applicability. Such applications range from future distributed space systems, to sensor-webs, and rapid-response satellite systems. Architectures will be explored that strike a balance between modularity and integration while preserving the MR(sup 2) paradigm. Modularity versus integration has always been a point of contention when approaching a design: whereas one-of-a-kind missions may require close integration resulting in performance optimization, multiple and flexible application spacecraft benefit &om modularity, resulting in maximum flexibility. The process of building spacecraft rapidly (< 7 days), requires a concerted and methodical look at system integration and test processes and pitfalls. Although the concept of modularity is not new and was first developed in the 1970s by NASA's Goddard Space Flight Center (Multi-Mission Modular Spacecraft), it was never modernized and was eventually abandoned. Such concepts as the Rapid Spacecraft Development Office (RSDO) became the preferred method for acquiring satellites. Notwithstanding, over the past 30 years technology has advanced considerably, and the time is ripe to reconsider modularity in its own right, as enabler of R(sup 2), and as a key element of transformational systems. The MR2 architecture provides a competitive advantage over the old modular approach in its rapid response to market needs that are difficult to predict both from the perspectives of evolving technology, as well as mission and application requirements.

Esper, Jaime↗
NASA NTRS1986

Sensitivity analysis and approximation methods for general eigenvalue problems

Optimization of dynamic systems involving complex non-hermitian matrices is often computationally expensive. Major contributors to the computational expense are the sensitivity analysis and reanalysis of a modified design. The present work seeks to alleviate this computational burden by identifying efficient sensitivity analysis and approximate reanalysis methods. For the algebraic eigenvalue problem involving non-hermitian matrices, algorithms for sensitivity analysis and approximate reanalysis are classified, compared and evaluated for efficiency and accuracy. Proper eigenvector normalization is discussed. An improved method for calculating derivatives of eigenvectors is proposed based on a more rational normalization condition and taking advantage of matrix sparsity. Important numerical aspects of this method are also discussed. To alleviate the problem of reanalysis, various approximation methods for eigenvalues are proposed and evaluated. Linear and quadratic approximations are based directly on the Taylor series. Several approximation methods are developed based on the generalized Rayleigh quotient for the eigenvalue problem. Approximation methods based on trace theorem give high accuracy without needing any derivatives. Operation counts for the computation of the approximations are given. General recommendations are made for the selection of appropriate approximation technique as a function of the matrix size, number of design variables, number of eigenvalues of interest and the number of design points at which approximation is sought.

Murthy, D. V.↗
NASA NTRS2001

Rock Driller

The next series of planetary exploration missions require a method of extracting rock and soil core samples. Therefore a prototype ultrasonic core driller (UTCD) was developed to meet the constraints of Small Bodies Exploration and Mars Sample Return Missions. The constraints in the design are size, weight, power, and axial loading. The ultrasonic transducer requires a relatively low axial load, which is one of the reasons this technology was chosen. The ultrasonic generator breadboard section can be contained within the 5x5x3 limits and weighs less than two pounds. Based on results attained the objectives for the first phase were achieved. A number of transducer probes were made and tested. One version only drills, and the other will actually provide a small core from a rock. Because of a more efficient transducer/probe, it will run at very low power (less than 5 Watts) and still drill/core. The prototype generator was built to allow for variation of all the performance-effecting elements of the transducer/probe/end effector, i.e., pulse, duty cycle, frequency, etc. The heart of the circuitry is what will be converted to a surface mounted board for the next phase, after all the parameters have been optimized and the microprocessor feedback can be installed.

Peterson, Thomas M.↗
DOE OSTI2024

Optimization and Multimachine Learning Algorithms to Predict Nanometal Surface Area Transfer Parameters for Gold and Silver Nanoparticles

Interactions between gold metallic nanoparticles and molecular dyes have been well described by the nanometal surface energy transfer (NSET) mechanism. However, the expansion and testing of this model for nanoparticles of different metal composition is needed to develop a greater variety of nanosensors for medical and commercial applications. In this study, the NSET formula was slightly modified in the size-dependent dampening constant and skin depth terms to allow for modeling of different metals as well as testing the quenching effects created by variously sized gold, silver, copper, and platinum nanoparticles. Overall, the metal nanoparticles followed more closely the NSET prediction than for Förster resonance energy transfer, though scattering effects began to occur at 20 nm in the nanoparticle diameter. To further improve the NSET theoretical equation, an attempt was made to set a best-fit line of the NSET theoretical equation curve onto the Au and Ag data points. An exhaustive grid search optimizer was applied in the ranges for two variables, 0.1≤C≤2.0 and 0≤α≤4, representing the metal dampening constant and the orientation of donor to the metal surface, respectively. Three different grid searches, starting from coarse (entire range) to finer (narrower range), resulted in more than one million total calculations with values C=2.0 and α=0.0736. The results improved the calculation, but further analysis needed to be conducted in order to find any additional missing physics. With that motivation, two artificial intelligence/machine learning (AI/ML) algorithms, multilayer perception and least absolute shrinkage and selection operator regression, gave a correlation coefficient, R2, greater than 0.97, indicating that the small dataset was not overfitting and was method-independent. This analysis indicates that an investigation is warranted to focus on deeper physics informed machine learning for the NSET equations.

Demers, Steven M. E. (ORCID:0000000192213246)↗
DOE OSTI2023

Exploring Parameter Redundancy in the Unitary Coupled-Cluster Ansätze for Hybrid Variational Quantum Computing

One of the commonly used chemical-inspired approaches in variational quantum computing is the unitary coupled-cluster (UCC) ansatze. Despite being a systematic way of approaching the exact limit, the number of parameters in the standard UCC ansatze exhibits unfavorable scaling with respect to the system size, hindering its practical use on near-term quantum devices. Efforts have been taken to propose some variants of UCC ansatze with better scaling. In this paper we explore the parameter redundancy in the preparation of unitary coupled-cluster singles and doubles (UCCSD) ansatze employing spin-adapted formulation, small amplitude filtration, and entropy-based orbital selection approaches. Numerical results of using our approach on some small molecules have exhibited a significant cost reduction in the number of parameters to be optimized and in the time to convergence compared with conventional UCCSD-VQE simulations. Further, we also discuss the potential application of some machine learning techniques in further exploring the parameter redundancy, providing a possible direction for future studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Technical and Imaging Outcomes from the UK Registry of Prostate Artery Embolization (UK-ROPE) Study: Focusing on Predictors of Clinical Success

Introduction: The UK Registry of Prostate Artery Embolization (UK-ROPE) was a prospective, multicentre study comparing PAE against surgical therapies for symptomatic benign prostatic hyperplasia (BPH). A wealth of data was collected supplementary to the main study outcomes which provide a snapshot of UK PAE practice. We aimed to interpret these data in the hope of providing insight into factors which affect clinical outcome and radiation dose. Methods: 216 patients (mean age 66, mean IPSS 21.3) undergoing PAE at 20 British centres from July 2014 to January 2016 were prospectively followed up to 12 months with retrospective analysis of the data. Technical outcome was evaluated based on procedural and fluoroscopy times, skin dose and dose area product (DAP). Clinical outcome was evaluated through collection of Qmax, IPSS reduction and prostate volume reduction. Multiple analysis of variance (MANOVA) was used to assess the significance of various patients and procedural factors on clinical outcome and patient dose. Results: Significant predictors of technical outcome which affected patient skin dose included severity of CTA-detected atheroma (p < 0.001), the practitioner (p < 0.001) and use of protective coil embolization (p = 0.019). Predictors of clinical outcome included initial prostate size (dichotomized into groups > 80 ml and = <80 ml, d = 1, p = 0.0138), embolic agent (spherical particles < 300 nm performed best, p = 0.01) and number of arteries embolized (IPSS reduction of 32.9% in unilateral PAE versus 54.4% for bilateral PAE, p = 0.026). Conclusion: We have identified several important factors which are associated with improved clinical outcome and increased patient dose which we hope will facilitate optimal patient selection and encourage improved embolization technique.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRSDate not supplied

A Cost and Benefit Analysis of Orbital Debris Remediation, Mitigation,Tracking, and Characterization

Orbital debris may collide with crewed and robotic spacecraft, placing them at risk. The wide range of debris, from 9,000-kilogram rocket bodies to millions of millimeter-size debris, has led to a similarly wide range of proposed actions for addressing the risks posed by debris. However, the costs and benefits of these actions have historically been unknown. This is a challenge for decision makers who are choosing which actions to support through technology development or policy changes. NASA’s Office of Technology, Policy, and Strategy is addressing these technical and economic uncertainties by building a capability to (1) complete rigorous calculations of the net present value of each action, (2) identify an optimal portfolio of actions to reduce risk, and (3) quantitatively analyze policies related to space sustainability. This report describes our progress toward that capability and to solicit feedback from the space and economic communities. Our previous work(Colvin, Karcz, and Wusk. 2023), referred to here as Phase 1, assessed the costs and benefits of performing debris remediation on operationally relevant timescales. The current paper summarizes the work of Locke and Colvin (2024), which contains major updates to the risk model used in Phase 1 and expands the breadth of actions considered to include mitigating the creation of debris, improving the ability to track debris, and more methods for cleaning up existing debris. We demonstrate that our approach of measuring risks in dollars allows for the effectiveness of seemingly incommensurate actions to be compared and generates insights that other approaches to measuring risk have missed.

Orbital Debris↗
NASA NTRS1993

Analysis of multigrid methods on massively parallel computers: Architectural implications

We study the potential performance of multigrid algorithms running on massively parallel computers with the intent of discovering whether presently envisioned machines will provide an efficient platform for such algorithms. We consider the domain parallel version of the standard V cycle algorithm on model problems, discretized using finite difference techniques in two and three dimensions on block structured grids of size 10(exp 6) and 10(exp 9), respectively. Our models of parallel computation were developed to reflect the computing characteristics of the current generation of massively parallel multicomputers. These models are based on an interconnection network of 256 to 16,384 message passing, 'workstation size' processors executing in an SPMD mode. The first model accomplishes interprocessor communications through a multistage permutation network. The communication cost is a logarithmic function which is similar to the costs in a variety of different topologies. The second model allows single stage communication costs only. Both models were designed with information provided by machine developers and utilize implementation derived parameters. With the medium grain parallelism of the current generation and the high fixed cost of an interprocessor communication, our analysis suggests an efficient implementation requires the machine to support the efficient transmission of long messages, (up to 1000 words) or the high initiation cost of a communication must be significantly reduced through an alternative optimization technique. Furthermore, with variable length message capability, our analysis suggests the low diameter multistage networks provide little or no advantage over a simple single stage communications network.

Matheson, Lesley R.↗
NASA NTRS2015

Evolving Design Criteria for Very Large Aperture Space-Based Telescopes and Their Influence on the Need for Integrated Tools in the Optimization Process

NASA's Advanced Mirror Technology Development (AMTD) program has been developing the means to design and build the future generations of space based telescopes. With the nearing completion of the James Webb Space Telescope (JWST), the astrophysics community is already starting to define the requirements for follow on observatories. The restrictions of available launch vehicles and the possibilities of planned future vehicles have fueled the competition between monolithic primaries (with better optical quality) and segmented primaries (with larger apertures, but with diffraction, costs and figure control issues). Regardless of the current shroud sizes and lift capacities, these competing architectures share the need for rapid design tools. As part of the AMTD program a number of tools have been developed and tested to speed up the design process. Starting with the Arnold Mirror Modeler (which creates Finite Element Models (FEM) for structural analysis) and now also feeds these models into thermal stability analyses. They share common file formats and interchangeable results. During the development of the program, numerous trade studies were created for 4 meter and 8 meter monolithic primaries, complete with support systems. Evaluation of these results has led to a better understanding of how the specification drives the results. This paper will show some of the early trade studies for typical specification requirements such as lowest mirror bending frequency and suspension system lowest frequency. The results use representative allowable stress values for each mirror substrate material and construction method and generic material properties. These studies lead to some interesting relationships between feasible designs and the realities of actually trying to build these mirrors. Much of the traditional specifications were developed for much smaller systems, where the mass and volume of the primary where a small portion of the overall satellite. JWST shows us that as the aperture grows, the primary takes up the majority of the mass and volume and the established rules need to be adjusted. For example, a small change in lowest frequency requirement can change the cost by millions of dollars.

Space-based telescopes↗
NASA NTRS2015

Evolving Design Criteria for Very Large Aperture Space Based Telescopes and Their Influence on the Need for Integrated Tools in the Optimization Process

NASA's Advanced Mirror Technology Development (AMTD) program has been developing the means to design and build the future generations of space based telescopes. With the nearing completion of the James Webb Space Telescope (JWST), the astrophysics community is already starting to define the requirements for follow-on observatories. The restrictions of available launch vehicles and the possibilities of planned future vehicles have fueled the competition between monolithic primaries (with better optical quality) and segmented primaries (with larger apertures, but with diffraction, costs and figure control issues). Regardless of the current shroud sizes and lift capacities, these competing architectures share the need for rapid design tools. As part of the AMTD program a number of tools have been developed and tested to speed up the design process. Starting with the Arnold Mirror Modeler (which creates Finite Element Models (FEM) for structural analysis) and now also feeds these models into thermal stability analyses. They share common file formats and interchangeable results. During the development of the program, numerous trade studies were created for 4-meter and 8-meter monolithic primaries, complete with support systems. Evaluation of these results has led to a better understanding of how the specification drives the results. This paper will show some of the early trade studies for typical specification requirements such as lowest mirror bending frequency and suspension system lowest frequency. The results use representative allowable stress values for each mirror substrate material and construction method and generic material properties. These studies lead to some interesting relationships between feasible designs and the realities of actually trying to build these mirrors. Much of the traditional specifications were developed for much smaller systems, where the mass and volume of the primary where a small portion of the overall satellite. JWST shows us that as the aperture grows, the primary takes up the majority of the mass and volume and the established rules need to be adjusted. For example, a small change in lowest frequency requirement can change the cost by millions of dollars. The paper uses numerous trade studies created during the software development phase of the Arnold Mirror Modeler to illustrate the influences of system specifications on the design space. The future telescopes will require better performance, stability and documented feasibility to meet the hurdles of today's budget and schedules realities. AMTD is developing the tools, but the basic system planning mentality also has to adopt to the requirements of these very large and complex physical structures.

Arnold, William R., Sr.↗
DOE OSTI2024

Tuning Surface Adhesion Using Grayscale Electron-beam Lithography

Surface texturing of manufactured products tailors their properties, such as friction, adhesion, biocompatibility, or fluid interactions. However, advancements in this area are largely the result of trial-and-effort testing and generally lack a science-guided framework for determining the surface topography that will optimize performance. The present investigation explores grayscale electron-beam lithography as a means to create multiscale surface patterns to control surface performance. Here, we created and characterized a set of surface textures on a silicon wafer; the textures were superpositions of sine waves of varying wavelengths and amplitudes. First, the multiscale topography of the patterned surface was characterized, using profilometry and atomic force microscopy, to understand its fidelity to the designed-in pattern. The results of this analysis demonstrated how grayscale lithography accurately controlled the lateral size of features but was less precise on the vertical height of the surface, and also introduced inherent roughness below the scale of patterning. Second, a micromechanical tester was used to characterize the adhesion of the surfaces with large-scale polished silicon spheres. The results showed that adhesion could be tailored, with significant contribution from all of the designed-in length scales of topography. The strength of adhesion did not correlate with conventional roughness parameters but could be accurately modeled using simple numerical integration. Taken together, this investigation demonstrates the promise and challenges of grayscale e-beam lithography with multiscale patterns as a method for the tailoring of surface performance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

CSPlib: A performance portable parallel software toolkit for analyzing complex kinetic mechanisms

Computational singular perturbation (CSP) is a method to analyze dynamical systems. It targets the decoupling of fast and slow dynamics using an alternate linear expansion of the right-hand side of the governing equations based on eigenanalysis of the associated Jacobian matrix. This representation facilitates diagnostic analysis, detection and control of stiffness, and the development of simplified models. For this work, we have implemented CSP in a C++ open-source library CSPlib using the Kokkos parallel programming model to address portability across diverse heterogeneous computing platforms, i.e., multi/many-core CPUs and GPUs. We describe the CSPlib implementation and present its computational performance across different computing platforms using several test problems. Specifically, we test the CSPlib performance for a constant pressure ignition reactor model on different architectures, including IBM Power 9, Intel Xeon Skylake, and NVIDIA V100 GPU. The size of the chemical kinetic mechanism is varied in these tests. As expected, the Jacobian matrix evaluation, the eigensolution of the Jacobian matrix, and matrix inversion are the most expensive computational tasks. When considering the higher throughput characteristic of GPUs, GPUs performs better for small matrices with higher occupancy rate. CPUs gain more advantages from the higher performance of well-tuned and optimized linear algebra libraries such as OpenBLAS.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Reverse-Mode Algorithmic Differentiation of an OpenMP-Parallel Compressible Flow Solver

y Reverse-mode algorithmic differentiation (AD) is an established method for obtaining adjoint derivatives of computer simulation applications. In computational fluid dynamics (CFD), adjoint derivatives of a cost function output such as drag or lift with respect to design parameters such as surface coordinates or geometry control points are a key ingredient for shape optimization, uncertainty quantification and flow control. The computational cost of CFD applications and their derivatives makes it essential to use high-performance computing hardware efficiently, including multi- and many-core architectures. Nevertheless, OpenMP is not supported in most AD tools, and previously shown methods achieve poor scalability of the derivative code. We present the AD of an OpenMP-parallelized finite volume compressible flow solver for unstructured meshes. Our approach enables us to reuse the parallelization of the original code in the computation of adjoint derivatives. The method works by identifying code segments that can be differentiated in reverse-mode without changing their memory access pattern. The OpenMP parallelization is integrated into the derivative code during the build process in a way that is robust to modifications of the original code and independent of the OpenMP support of the differentiation tool. We show the scalability of our adjoint CFD solver on test cases ranging from thousands to millions of finite volume mesh cells on CPUs with up to 16 threads as well as on an Intel XeonPhi card with 236 threads. We demonstrate that our approach is more practical to implement for production-sized CFD codes and produces more efficient adjoint derivative code than previously shown AD methods.

Huckelheim, Jan↗