Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shape from gradient”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Adjoint approach to calculating shape gradients for three-dimensional magnetic confinement equilibria

The shape gradient quantifies the change in some figure of merit resulting from differential perturbations to a shape. Shape gradients can be applied to gradient-based optimization, sensitivity analysis and tolerance calculation. An efficient technique for computing the shape gradient for toroidal three-dimensional magnetohydrodynamic (MHD) equilibria is introduced. The method is based on the self-adjoint property of the equations for driven perturbations of MHD equilibria and is similar to the Onsager symmetry of transport coefficients. Two versions of the shape gradient are considered. One describes the change in a figure of merit due to an arbitrary displacement of the outer flux surface; the other describes the change in the figure of merit due to the displacement of a coil. The method is implemented for several example figures of merit and compared with direct calculation of the shape gradient. In these examples the adjoint method reduces the number of equilibrium computations by factors of $O(N)$ , where $N$ is the number of parameters used to describe the outer flux surface or coil shapes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Computing the shape gradient of stellarator coil complexity with respect to the plasma boundary

Coil complexity is a critical consideration in stellarator design. The traditional two-step optimization approach, in which the plasma boundary is optimized for physics properties and the coils are subsequently optimized to be consistent with this boundary, can result in plasma shapes which cannot be produced with sufficiently simple coils. To address this challenge, we propose a method to incorporate considerations of coil complexity in the optimization of the plasma boundary. Coil complexity metrics are computed from the current potential solution obtained with the REGCOIL code (Landreman, Nucl. Fusion , vol. 57, 2017, 046003). While such metrics have previously been included in derivative-free fixed-boundary optimization (Drevlak et al. , Nucl. Fusion , vol. 59, 2018, 016010), we compute the local sensitivity of these metrics with respect to perturbations of the plasma boundary using the shape gradient (Landreman & Paul, Nucl. Fusion , vol. 58, 2018, 076023). We extend REGCOIL to compute derivatives of these metrics with respect to parameters describing the plasma boundary. In keeping with previous research on winding surface optimization (Paul et al. , Nucl. Fusion , vol. 58, 2018, 076015), the shape derivatives are computed with a discrete adjoint method. In contrast with the previous work, derivatives are computed with respect to the plasma surface parameters rather than the winding surface parameters. To further reduce the resolution required to compute the shape gradient, we present a more efficient representation of the plasma surface which uses a single Fourier series to describe the radial distance from a coordinate axis and a spectrally condensed poloidal angle. This representation is advantageous over the standard cylindrical representation used in the VMEC code (Hirshman & Whitson, Phys. Fluids , vol. 26, 1983, pp. 3553–3568), as it provides a uniquely defined poloidal angle, eliminating a null space in the optimization of the plasma surface. In comparison with previous spectral condensation methods (Hirshman & Breslau, Phys. Plasmas , vol. 5, 1998, p. 2664), the modified poloidal angle is obtained algebraically rather than through the solution of a nonlinear optimization problem. The resulting shape gradient highlights features of the plasma boundary that are consistent with simple coils and can be used to couple coil and fixed-boundary optimization.

Physics↗
NASA NTRS1995

Effect of modulator sorption on gradient shape in ion-exchange chromatography

Mobile phase additives, or modulators, are used in gradient elution chromatography to facilitate separation and reduce separation time. The modulators are usually assumed to be linearly adsorbed or unadsorbed. Here, the consequences of nonlinear modulator adsorption are examined for ion-exchange gradient elution through a series of simulations. Even when the buffer salt is identical to the modulator salt, gradient deformation is observed; the extent of deformation increases as the volume of the feed is increased. When the modulator salt is different from the buffer salt, unusual effects are observed, and the chromatograms are quite different from those predicted by classical gradient elution theory. In particular, local increases in the buffer concentration are found between feed bands, and serve to improve the separation. These effects become more pronounced as the feed volume increases, and could therefore prove valuable in preparative applications.

NASA Discipline Life Support Systems↗
NASA NTRS1997

Efficient Gradient-Based Shape Optimization Methodology Using Inviscid/Viscous CFD

The formerly developed preconditioned-biconjugate-gradient (PBCG) solvers for the analysis and the sensitivity equations had resulted in very large error reductions per iteration; quadratic convergence was achieved whenever the solution entered the domain of attraction to the root. Its memory requirement was also lower as compared to a direct inversion solver. However, this memory requirement was high enough to preclude the realistic, high grid-density design of a practical 3D geometry. This limitation served as the impetus to the first-year activity (March 9, 1995 to March 8, 1996). Therefore, the major activity for this period was the development of the low-memory methodology for the discrete-sensitivity-based shape optimization. This was accomplished by solving all the resulting sets of equations using an alternating-direction-implicit (ADI) approach. The results indicated that shape optimization problems which required large numbers of grid points could be resolved with a gradient-based approach. Therefore, to better utilize the computational resources, it was recommended that a number of coarse grid cases, using the PBCG method, should initially be conducted to better define the optimization problem and the design space, and obtain an improved initial shape. Subsequently, a fine grid shape optimization, which necessitates using the ADI method, should be conducted to accurately obtain the final optimized shape. The other activity during this period was the interaction with the members of the Aerodynamic and Aeroacoustic Methods Branch of Langley Research Center during one stage of their investigation to develop an adjoint-variable sensitivity method using the viscous flow equations. This method had algorithmic similarities to the variational sensitivity methods and the control-theory approach. However, unlike the prior studies, it was considered for the three-dimensional, viscous flow equations. The major accomplishment in the second period of this project (March 9, 1996 to March 8, 1997) was the extension of the shape optimization methodology for the Thin-Layer Navier-Stokes equations. Both the Euler-based and the TLNS-based analyses compared with the analyses obtained using the CFL3D code. The sensitivities, again from both levels of the flow equations, also compared very well with the finite-differenced sensitivities. A fairly large set of shape optimization cases were conducted to study a number of issues previously not well understood. The testbed for these cases was the shaping of an arrow wing in Mach 2.4 flow. All the final shapes, obtained either from a coarse-grid-based or a fine-grid-based optimization, using either a Euler-based or a TLNS-based analysis, were all re-analyzed using a fine-grid, TLNS solution for their function evaluations. This allowed for a more fair comparison of their relative merits. From the aerodynamic performance standpoint, the fine-grid TLNS-based optimization produced the best shape, and the fine-grid Euler-based optimization produced the lowest cruise efficiency.

Baysal, Oktay↗
DOE OSTI2019

Study of non-isothermal liquid evaporation in synthetic micro-pore structures with hybrid lattice Boltzmann model

Non-isothermal liquid evaporation in micro-pore structures is studied experimentally and numerically using the lattice Boltzmann method. A hybrid thermal entropic multiple-relaxation-time multiphase lattice Boltzmann model (T-EMRT-MP LBM) is implemented and validated with experiments of droplet evaporation on a heated hydrophobic substrate. Then liquid evaporation is investigated in two specific pore structures, i.e. spiral-shaped and gradient-shaped micro-pillar cavities, referred to as SMS and GMS, respectively. In SMS, the liquid receding front follows the spiral pattern; while in GMS, the receding front moves layer by layer from the pillar rows with large pitch to the rows with small one. Both simulations agree well with experiments. Moreover, evaporative cooling effects in liquid and vapour are observed and explained with simulation results. Quantitatively, in both SMS and GMS, the change of liquid mass with time coincides with experimental measurements. The evaporation rate generally decreases slightly with time mainly because of the reduction of liquid–vapour interface. Isolated liquid films in SMS increase the evaporation rate temporarily resulting in local peaks in evaporation rate. Reynolds and capillary numbers show that the liquid internal flow is laminar and that the capillary forces are dominant resulting in menisci pinned to the pillars. Similar Péclet number is found in simulations and experiments, indicating a diffusive type of heat, liquid and vapour transport. Our numerical and experimental studies indicate a method for controlling liquid evaporation paths in micro-pore structures and maintaining high evaporation rate by specific geometry designs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Investigating profile stiffness and critical gradients in shaped TCV discharges using local gyrokinetic simulations of turbulent transport

The experimental observation made on the TCV tokamak of a significant confinement improvement in plasmas with negative triangularity ( δ < 0) compared to those with standard positive triangularity has been interpreted in terms of different degrees of profile stiffness (Sauter et al. Phys. Plasmas 21 055906, 2014) and/or different critical gradients. Employing the Eulerian gyrokinetic code GENE (Jenko et al. Phys. Plasmas 7 1904, 2000), profile stiffness and critical gradients are studied under TCV relevant conditions. For the considered experimental discharges, trapped electron modes (TEMs) and electron temperature gradient (ETG) modes are the dominant microinstabilities, with the latter providing a significant contribution to the non-linear electron heat fluxes near the plasma edge. Two series of simulations with different levels of realism are performed, addressing the question of profile stiffness at various radial locations. Retaining finite collisionality, impurities and electromagnetic effects, as well as the physical electron-to-ion mass ratio are all necessary in order to approach the experimental flux measurements. However, flux-tube simulations are unable to fully reproduce the TCV results, pointing towards the need to carry out radially nonlocal (global) simulations, i.e. retaining finite machine size effects, in a future study. Some conclusions about the effect of triangularity can nevertheless be drawn based on the flux-tube results. In particular, the importance of considering the sensitivity to both temperature and density gradient is shown. The flux tube results show an increase of the critical gradients towards the edge, further enhanced when δ < 0, and they also appear to indicate a reduction of profile stiffness towards plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1997

Methods for Sensitivity Analysis in Gradient-Based Shape Optimization: A Review

CFD is now being widely used in the design of mechanical and aerospace systems and their components. Improving the accuracy and efficiency of this process can reduce its cost and increase its effectiveness. Other motivating factors that push and pull CFD into design optimization are; (1) the need for design methodologies (single- and multi-disciplinary) to reduce cost and time in developing new fluid and aero systems for fast response to market changes; (2) new design requirements outside existing database; (3) new CFD algorithms and growth in computing power; (4) need for increased product quality and reliability; (5) lack of test facilities; (6) limitations and cost of testing techniques; (7) limitations of existing design practices; and (8) prohibitive cost of and conflicting designs derived from isolated component analyses and design, which are repetitive in nature and omit mutual interference.

Baysal, Oktay↗
DOE OSTI2020

Modular gradient-free shaped charge

A shaped charge produces a constant velocity jet. The shaped charge is comprised of individual modules which can be assembled to produce a constant velocity jet of arbitrary length. The resulting jet speed is approximately twice the detonation velocity and independent of the liner material. The modular design also allows different liner materials to be used sequentially in the same jet.

42 ENGINEERING↗
NASA NTRS1997

Gradient-Based Aerodynamic Shape Optimization Using ADI Method for Large-Scale Problems

A gradient-based shape optimization methodology, that is intended for practical three-dimensional aerodynamic applications, has been developed. It is based on the quasi-analytical sensitivities. The flow analysis is rendered by a fully implicit, finite volume formulation of the Euler equations.The aerodynamic sensitivity equation is solved using the alternating-direction-implicit (ADI) algorithm for memory efficiency. A flexible wing geometry model, that is based on surface parameterization and platform schedules, is utilized. The present methodology and its components have been tested via several comparisons. Initially, the flow analysis for for a wing is compared with those obtained using an unfactored, preconditioned conjugate gradient approach (PCG), and an extensively validated CFD code. Then, the sensitivities computed with the present method have been compared with those obtained using the finite-difference and the PCG approaches. Effects of grid refinement and convergence tolerance on the analysis and shape optimization have been explored. Finally the new procedure has been demonstrated in the design of a cranked arrow wing at Mach 2.4. Despite the expected increase in the computational time, the results indicate that shape optimization, which require large numbers of grid points can be resolved with a gradient-based approach.

Pandya, Mohagna J.↗
DOE OSTI2014

Comparison of two-dimensional integration methods for shape reconstruction from gradient data

As a requisite and key step in some gradient-based measurement techniques, the reconstruction of the shape, more generally the scalar potential, from the measured gradient data has been studied for many years. Here, in this work, three types of two-dimensional integration methods are compared under various conditions. The merits and drawbacks of each integration method are consequently revealed to provide suggestions in selection of a proper integration method for a particular application.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Molecular concentration field design using closed-form steady-state solutions

Control over spatial concentration fields represents a fundamental challenge in designing synthetic biological systems and programmable soft materials. While nature creates morphogen gradients that orchestrate complex developmental processes, synthetic approaches have largely relied on empirical optimization and computationally intensive simulations. Here, we present an analytical framework for steady-state concentration fields generated by finite-sized localized sources in diffusion–degradation systems and derive closed-form solutions for one-, two-, and three-dimensional geometries. By expressing these solutions in dimensionless form, we show that gradient steepness and spatial structure are organized by the Thiele modulus, which captures the competition between diffusion and degradation length scales. The analysis reveals distinct design regimes: in degradation-dominated systems, gradient shape is governed by exponential decay and becomes dimension-independent, whereas in diffusion-dominated systems, gradient magnitude and extent follow dimension-dependent power-law scaling. Building on these results, we introduce a quantitative design strategy that uses threshold-based criteria to program concentration ranges by tuning physically accessible parameters, most directly the production rate, while holding transport and degradation properties fixed. Comparisons with numerical solutions and reported experimental systems demonstrate consistency with the predicted scaling behavior. Together, this work provides a generalizable and physically transparent framework for designing steady-state concentration fields in synthetic biological and soft matter systems, enabling predictive control of gradient-mediated organization without reliance on extensive numerical optimization.

Kim, Dong Woo [Johns Hopkins University, Baltimore↗
DOE OSTI2021

An adjoint method for determining the sensitivity of island size to magnetic field variations

An adjoint method to calculate the gradient of island width in stellarators is presented and applied to a set of magnetic field configurations. The underlying method for calculation of the island width is that of Cary & Hanson ( Phys. Fluids B, vol. 3, issue 4, 1991, pp. 1006–1014) (with a minor modification), and requires that the residue of the island centre be small. Therefore, the gradient of the residue is calculated in addition. Both the island width and the gradient calculations are verified using an analytical magnetic field configuration introduced by Reiman & Greenside ( Comput. Phys. Commun. , vol. 43, issue 1, 1986, pp. 157–167). The method is also applied to the calculation of the shape gradient of the width of a magnetic island in a National Compact Stellarator Experiment (NCSX) vacuum configuration with respect to positions on a coil. A gradient-based optimization is applied to a magnetic field configuration studied by Hanson & Cary ( Phys. Fluids , vol. 27, issue 4, 1984, pp. 767–769) to minimize stochasticity by adding perturbations to a pair of helical coils. Although only vacuum magnetic fields and an analytical magnetic field model are considered in this work, the adjoint calculation of the island width gradient could also be applied to a magnetohydrodynamic (MHD) equilibrium if the derivative of the magnetic field, with respect to the equilibrium parameters, is known. Using the island width gradient calculation presented here, more general gradient-based optimization methods can be applied to design stellarators with small magnetic islands. Moreover, the sensitivity of the island size may itself be optimized to ensure that coil tolerances, with respect to island size, are kept as high as possible.

Physics↗
DOE OSTI2021

Dataset from: "Adjoint methods for quasisymmetry of vacuum fields on a surface"

Adjoint methods can speed up stellarator optimisation by providing gradient information more efficiently compared to finite-difference evaluations. Adjoint methods are herein applied to vacuum magnetic fields, with objective functions targeting quasisymmetry and a rotational transform value on a surface. To measure quasisymmetry, a novel way of evaluating approximate flux coordinates on a single flux surface without the assumption of a neighbourhood of flux surfaces is proposed. The shape gradients obtained from the adjoint formalism are evaluated numerically and verified against finite-difference evaluations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Adjoint methods for quasi-symmetry of vacuum fields on a surface

Adjoint methods can speed up stellarator optimisation by providing gradient information more efficiently compared with finite-difference evaluations. Adjoint methods are herein applied to vacuum magnetic fields, with objective functions targeting quasi-symmetry and a rotational transform value on a surface. Here to measure quasi-symmetry, a novel way of evaluating approximate flux coordinates on a single flux surface without the assumption of a neighbourhood of flux surfaces is proposed. The shape gradients obtained from the adjoint formalism are evaluated numerically and verified against finite-difference evaluations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗