Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shape constraints”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2024

Evolving Particles in the 2022 Hunga-Tonga Hunga-Ha’apai Volcano Eruption Plume

The Multi-angle Imaging SpectroRadiometer (MISR) aboard NASA’s Terra satellite observed the Hunga Tonga—Hunga Ha’apai (HTHH) 15 January eruption plume on eight occasions between 15 and 23 January 2022. From the MISR multi-angle, multi-spectral imagery we retrieve aerosol plume height geometrically, along with plume-level motion vectors, and derive radiometrically constraints on particle effective size, shape, and light-absorption properties. Parts of two downwind aerosol layers were observed in different places and times, one concentrated in the upper troposphere (11-18 km ASL), and a mid-stratosphere layer~23 –30+ km ASL. After the initial day (1/15), the retrievals identified only spherical, non-light-absorbing particles, typical of volcanic sulfate/water particles. The near-tropopause plume particles show constant, medium-small (several tenths of a micron) effective size over four days. The mid-stratosphere particles were consistently smaller, but retrieved effective particle size increased between 1/17 and 1/23, though they might have decreased slightly on 1/22. As a vast amount of water was also injected into the stratosphere by this eruption, models predicted relatively rapid sulfate particle growth from the modest amounts of SO2 gas injected by the eruption to high altitudes along with the water (Zhu et al, 2022). MISR observations up to ten days after the eruption are consistent with these model predictions. The possible decrease in stratospheric particle size after initial growth was likely caused by evaporation, as the plume mixed with drier, ambient air. Particles in the lower-elevation plume observed on 1/15 were larger than all the downwind aerosols and contained significant non-spherical (likely ash) particles.

aerosol plume height↗
NASA NTRS2013

Tracker Toolkit

This software can track multiple moving objects within a video stream simultaneously, use visual features to aid in the tracking, and initiate tracks based on object detection in a subregion. A simple programmatic interface allows plugging into larger image chain modeling suites. It extracts unique visual features for aid in tracking and later analysis, and includes sub-functionality for extracting visual features about an object identified within an image frame. Tracker Toolkit utilizes a feature extraction algorithm to tag each object with metadata features about its size, shape, color, and movement. Its functionality is independent of the scale of objects within a scene. The only assumption made on the tracked objects is that they move. There are no constraints on size within the scene, shape, or type of movement. The Tracker Toolkit is also capable of following an arbitrary number of objects in the same scene, identifying and propagating the track of each object from frame to frame. Target objects may be specified for tracking beforehand, or may be dynamically discovered within a tripwire region. Initialization of the Tracker Toolkit algorithm includes two steps: Initializing the data structures for tracked target objects, including targets preselected for tracking; and initializing the tripwire region. If no tripwire region is desired, this step is skipped. The tripwire region is an area within the frames that is always checked for new objects, and all new objects discovered within the region will be tracked until lost (by leaving the frame, stopping, or blending in to the background).

Lewis, Steven J.↗
NASA NTRS1984

X-ray emission from active galactic nuclei

The advantages of X-ray astronomy for studying active galactic nuclei (AGN) are discussed. If the X-ray emission from AGN arises from a region close to the central energy source X-ray observations may provide the best constraints on the central engine. The shape of the X-ray continuum gives information about the mechanism for photon generation; X-ray time variability data can constrain the size and mass of the continuum source; and X-ray occultation data gives constraints on the relative sizes of the continuum source and the intervening absorbing material. Since a fair fraction of the total energy of an AGN is emitted at X-ray wavelengths, direct measurement of the amount and spectral form of this radiation is important for modeling of the optically emitting clouds.

Mushotzky, R.↗
NASA NTRS2009

Modeling and Simulation of Variable Mass, Flexible Structures

The advent of the new Ares I launch vehicle has highlighted the need for advanced dynamic analysis tools for variable mass, flexible structures. This system is composed of interconnected flexible stages or components undergoing rapid mass depletion through the consumption of solid or liquid propellant. In addition to large rigid body configuration changes, the system simultaneously experiences elastic deformations. In most applications, the elastic deformations are compatible with linear strain-displacement relationships and are typically modeled using the assumed modes technique. The deformation of the system is approximated through the linear combination of the products of spatial shape functions and generalized time coordinates. Spatial shape functions are traditionally composed of normal mode shapes of the system or even constraint modes and static deformations derived from finite element models of the system. Equations of motion for systems undergoing coupled large rigid body motion and elastic deformation have previously been derived through a number of techniques [1]. However, in these derivations, the mode shapes or spatial shape functions of the system components were considered constant. But with the Ares I vehicle, the structural characteristics of the system are changing with the mass of the system. Previous approaches to solving this problem involve periodic updates to the spatial shape functions or interpolation between shape functions based on system mass or elapsed mission time. These solutions often introduce misleading or even unstable numerical transients into the system. Plus, interpolation on a shape function is not intuitive. This paper presents an approach in which the shape functions are held constant and operate on the changing mass and stiffness matrices of the vehicle components. Each vehicle stage or component finite element model is broken into dry structure and propellant models. A library of propellant models is used to describe the distribution of mass in the fuel tank or Solid Rocket Booster (SRB) case for various propellant levels. Based on the mass consumed by the liquid engine or SRB, the appropriate propellant model is coupled with the dry structure model for the stage. Then using vehicle configuration data, the integrated vehicle model is assembled and operated on by the constant system shape functions. The system mode shapes and frequencies can then be computed from the resulting generalized mass and stiffness matrices for that mass configuration. The rigid body mass properties of the vehicle are derived from the integrated vehicle model. The coupling terms between the vehicle rigid body motion and elastic deformation are also updated from the constant system shape functions and the integrated vehicle model. This approach was first used to analyze variable mass spinning beams and then prototyped into a generic dynamics simulation engine. The resulting code was tested against Crew Launch Vehicle (CLV-)class problems worked in the TREETOPS simulation package and by Wilson [2]. The Ares I System Integration Laboratory (SIL) is currently being developed at the Marshall Space Flight Center (MSFC) to test vehicle avionics hardware and software in a hardware-in-the-loop (HWIL) environment and certify that the integrated system is prepared for flight. The Ares I SIL utilizes the Ares Real-Time Environment for Modeling, Integration, and Simulation (ARTEMIS) tool to simulate the launch vehicle and stimulate avionics hardware. Due to the presence of vehicle control system filters and the thrust oscillation suppression system, which are tuned to the structural characteristics of the vehicle, ARTEMIS must incorporate accurate structural models of the Ares I launch vehicle. The ARTEMIS core dynamics simulation models the highly coupled nature of the vehicle flexible body dynamics, propellant slosh, and vehicle nozzle inertia effects combined with mass and flexible body properties that vary significant with time during the flight. All forces that act on the vehicle during flight must be simulated, including deflected engine thrust force, spatially distributed aerodynamic forces, gravity, and reaction control jet thrust forces. These forces are used to excite an integrated flexible vehicle, slosh, and nozzle dynamics model for the vehicle stack that simulates large rigid body translations and rotations along with small elastic deformations. Highly effective matrix math operations on a distributed, threaded high-performance simulation node allow ARTEMIS to retain up to 30 modes of flex for real-time simulation. Stage elements that separate from the stack during flight are propagated as independent rigid six degrees of freedom (6DOF) bodies. This paper will present the formulation of the resulting equations of motion, solutions to example problems, and describe the resulting dynamics simulation engine within ARTEMIS.

Tobbe, Patrick A.↗
DOE OSTI2021

Two Dimensional Topology Optimization of Heat Exchangers with the Density and Level-Set Methods

We design heat exchangers using two topology optimization approaches: the density, i.e. volume fraction and level set methods. Our goal is to maximize the heat exchange between two fluids in separate channels while constraining the pressure drop across each channel. The heat exchanger is modeled with a coupled thermal-flow formulation. The flow is governed by an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation with high Peclet number. We solve one set of Stokes-Brinkman equations per fluid. Each Brinkman term in the flow equation serves to model the other phase as a solid, thereby preventing mixing. We first represent the solid and fluid phases using a volume fraction variable and apply a SIMP-like penalization in the Brinkman term to drive the optimization to a discrete design. The cost and constraint function derivatives are automatically calculated with the library pyadjoint and the optimization is performed by the Method of Moving Asymptotes. In a second optimization formulation, we use the level set approach to define the interface that separates the two fluids. Pyadjoint calculates the shape derivatives of the cost and constraint functions and the Hamilton-Jacobi advects the interface, allowing for topological changes. We present results in two dimensions and discuss the advantages and disadvantages of each approach.

42 ENGINEERING↗
NASA NTRS2002

Thickness Constraints on the Icy Shells of the Galilean Satellites from a Comparison of Crater Shapes

A thin outer ice shell on Jupiter's large moon Europa would imply easy exchange between the surface and any organic or biotic material in its putative subsurface ocean. The thickness of the outer ice shell is poorly constrained, however, with model-dependent estimates ranging from a few kilometers of depths of impact craters on Europa, Ganymede and Callisto that reveal two anomalous transitions in crater shape with diameter. The first transition is probably related to temperature-dependent ductility of the crust at shallow depths (7-8 km on Europa). The second transition is attributed to the influence of subsurface oceans on all three satellites, which constrains Europa's icy shell to be at least 19 km thick. The icy lithospheres of Ganymede and Callisto are equally ice-rich, but Europa's icy shell has a thermal structure about 0.25-0.5 times the thickness of Ganymede's or Callisto's shells, depending on epoch. The appearances of the craters on Europa are inconsistent with thin-ice-shell models and indicate that exchange of oceanic and surface material could be difficult.

Schenk, Paul M.↗
NASA NTRS1991

On the theory of the synthesis of offset dual-shaped reflectors - Case examples

In an earlier paper by Galindo-Israel et al. (1987), the geometrical optics (GO) principles, constraints, and requirements of the dual- and single-offset-shaped reflector synthesis problem were collected and developed into a set of nonlinear first-order PDEs. An extension of the methods by which solutions to the PDEs can be obtained is presented, together with several case examples. These examples are independently analyzed by GO and physical optics diffraction methods. The starting point for the integration over each reflector can be taken on the outer rim, at the center, or at an intermediate point-the intermediate starting point being the more general case. The utility of the speed of this synthesis method is demonstrated. For example, the program utilized completes a synthesis in less than 45 s on a 386/20 personal computer and within a fraction of a second on a mainframe Cray. This makes practical the incorporation of the synthesis into a search algorithm that can optimize one or more parameters of the reflector system. As an example, the optimization of the mapping equations for low cross polarization is discussed.

Galindo-Israel, Victor↗
NASA NTRS1999

Demonstration of Automatically-Generated Adjoint Code for Use in Aerodynamic Shape Optimization

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 limit is reached, or no further design improvement is possible due to active design variable bounds and/or constraints. The resulting shape parameters are then used by the grid generation code to define a new wing surface and computational grid. The lift-to-drag ratio and its gradient are computed for the new design by the automatically-generated adjoint codes. Several optimization iterations may be required to find an optimum wing shape. Results from two sample cases will be discussed. The reader should note that this work primarily represents a demonstration of use of automatically- generated adjoint code within an aerodynamic shape optimization. As such, little significance is placed upon the actual optimization results, relative to the method for obtaining the results.

Green, Lawrence↗
NASA NTRS2012

Computational Optimization of a Natural Laminar Flow Experimental Wing Glove

Computational optimization of a natural laminar flow experimental wing glove that is mounted on a business jet is presented and discussed. The process of designing a laminar flow wing glove starts with creating a two-dimensional optimized airfoil and then lofting it into a three-dimensional wing glove section. The airfoil design process does not consider the three dimensional flow effects such as cross flow due wing sweep as well as engine and body interference. Therefore, once an initial glove geometry is created from the airfoil, the three dimensional wing glove has to be optimized to ensure that the desired extent of laminar flow is maintained over the entire glove. TRANAIR, a non-linear full potential solver with a coupled boundary layer code was used as the main tool in the design and optimization process of the three-dimensional glove shape. The optimization process uses the Class-Shape-Transformation method to perturb the geometry with geometric constraints that allow for a 2-in clearance from the main wing. The three-dimensional glove shape was optimized with the objective of having a spanwise uniform pressure distribution that matches the optimized two-dimensional pressure distribution as closely as possible. Results show that with the appropriate inputs, the optimizer is able to match the two dimensional pressure distributions practically across the entire span of the wing glove. This allows for the experiment to have a much higher probability of having a large extent of natural laminar flow in flight.

Hartshom, Fletcher↗
DOE OSTI2024

Development of a Method for Shape Optimization for a Gas Turbine Fuel Injector Design Using Metal-Additive Manufacturing

Adjoint shape optimization has enabled physics-based optimal designs for aerodynamic surfaces. Additive manufacturing (AM) makes it possible to manufacture complex shapes. However, there has been a gap between optimal and manufacturable surfaces due to the inherent limitations of commercial computational fluid dynamics (CFD) codes to implement geometric constraints during adjoint computation. In such cases, the design sensitivities are exported and used to perform constrained shape modifications using parametric information stored in computer aided design (CAD) files to satisfy manufacturability constraints. However, modifying the design using adjoint methods in CFD solvers and performing constrained shape modification in CAD can lead to inconsistencies due to different shape parameterization schemes. This paper describes a method to enable the simultaneous optimization of the fluid domain and impose AM manufacturability constraints, resolving one of the key issues of geometry definition for isogeometric analysis. Similar to a grid convergence study, the proposed method verifies the consistencies between shape parameterization techniques present within commercial CAD and CFD software during mesh movement as a part of the adjoint shape optimization routine. By identifying the appropriate parameters essential to a shape optimization study, the error metric between the different parameterization techniques converges to demonstrate sufficient consistencies for justifiable exchange of data between CAD and CFD. For the identified shape optimization parameters, the error metric to measure the deviation between the two parameterization schemes lies within the AM laser-powder bed fusion (L-PBF) process tolerance. Additionally, comparison for subsequent objective function calculations between iterations of the optimization loop showed acceptable differences within 1% variation between the modified geometries obtained using the two parameterization schemes. This method provides justification for the use of multiphysics guided adjoint design sensitivities computed in CFD software to perform shape modifications in CAD to incorporate AM manufacturability constraints during the shape optimization loop such that optimal designs are also additively manufacturable.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

More of the bulk from extremal area variations

It was shown recently in (Bao N et al 2019 Class. Quantum Grav. 36 185002), building on work of Alexakis, Balehowksy, and Nachman (Alexakis S et al 2017 arXiv:1711.09379), that the geometry of (some portion of) a manifold with boundary is uniquely fixed by the areas of a foliation of two-dimensional disk-shaped surfaces anchored to the boundary. In the context of AdS/CFT, this implies that (a portion of) a four-dimensional bulk geometry can be fixed uniquely from the entanglement entropies of disk-shaped boundary regions, subject to several constraints. Here, we loosen some of these constraints, in particular allowing for the bulk foliation of extremal surfaces to be local and removing the constraint of disk topology; these generalizations ensure uniqueness of more of the deep bulk geometry by allowing for e.g. surfaces anchored on disconnected asymptotic boundaries, or HRT surfaces past a phase transition. We also explore in more depth the generality of the local foliation requirement, showing that even in a highly dynamical geometry like AdS-Vaidya it is satisfied.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1979

Minimization theory of induced drag subject to constraint conditions

Exact analytical solutions in terms of induced drag influence coefficients can be attained which define the spanwise loading with minimized induced drag, subject to specified constraint conditions, for any nonplanar wing shape or number of lift plus wing bending moment about a given wing span station. Example applications of the theory are made to a biplane, a wing in ground effect, a cruciform wing, a V-wing, a planar-wing winglet, and linked wingtips in formation flying. For minimal induced drag, the spanwise loading, relative to elliptic, is outboard for the biplane and is inboard for the wing in ground effect and for the planar-wing winglet. A spinoff of the triplane solution provides mathematically exact equations for downwash and sidewash about a planar vorticity sheet having an arbitrary loading distribution.

Deyoung, J.↗
NASA NTRS1981

Radar investigation of asteroids

Radar investigations were conducted of selected minor planets, including: (1) observations during 1981-82 of 10 potential targets (2 Pallas, 8 Flora, 12 Victoria, 15 Eunomia, 19 Fortuna, 22 Kalliope, 132 Aethra, 219 Thusnelda, 433 Eros, and 2100 Ra-Shalom); and (2) continued analyses of observational data obtained during 1980-81 for 10 other asteroids (4 Vesta, 7 Iris, 16 Psyche, 75 Eurydike, 97 Klotho, 216 Kleopatra, 1685 Toro, 1862 Apollo, 1865 Cerberus, and 1915 Quetzalcoatl). Scientific objectives include estimation of echo strength, polarization, spectral shape, spectral bandwidth, and Doppler shift. These measurements: (1) yield estimates of target size, shape, and spin vector; (2) place constraints on topography, morphology, and composition of the planetary surface; (3) yield refined estimates of target orbital parameters; (4) reveal the presence of asteroidal satellites.

Ostro, S. J.↗
NASA NTRS1990

Surface characterizations of color threshold

The paper evaluates how well three different parametric shapes, ellipsoids, rectangles, and parallelograms, serve as models of three-dimentional detection contours. The constraints of the procedures for deriving the best-fitting shapes on inferences about the theoretical visual detection mechanisms are described. Results of two statistical tests show that only the parallelogram fits the data with more precision than the variance in repeated threshold measurements, and thus provides a slightly better fit than the other two shapes. Nevertheless it does not serve as a better guide than the ellipsoidal model for interpolating from the measurements to thresholds in novel color directions.

Poirson, Allen B.↗
NASA NTRS1990

Optimization of linear, controlled structures

An approach to the integrated design of linear controlled structures, which is being investigated as part of NASA's controlled structures interaction (CSI) methodology development program, is presented. The approach is integrated in the sense that the structure and its controller are simultaneously designed. The design methodology uses constrained nonlinear optimization procedures based on analytically obtained gradients of the structural responses. Design of the controller is based on the so-called Q-parameterization theory, which parameterizes all closed-loop input/output maps achievable with stabiizing linear controllers. Very general objective and constraint functions are possible, and structural shape can be included in the design variables. This method has been partially implemented and demonstrated; early findings are reported.

Hunziker, K. Scott↗
NASA NTRS1993

Ionization states of metallic absorption-line systems in continua of quasars

Ionization states of metallic absorption-line systems in continua of quasars (QSO's) are studied, assuming that the metallic lines arise in gaseous halos of high-redshift galaxies in photoionization equilibrium under the background UV radiation, and constraints on the intensity and spectral shape of the UV radiation are obtained. Then a structure of absorbers suitable for all of the metallic absorption line systems are discussed.

Denda, Kiyomi↗
NASA NTRS1992

Motion planning of mobile multi-limb robotic systems subject to force and friction constraints

A method is presented to generate motions for a class of multilimb robotic systems enabling them to apply large static forces over large ranges of motion without saturating actuator effort limits, system-environment friction constraints, kinematic joint limits, or geometric workspace obstacles. The approach, termed the force-workspace (FW) approach, maps these constraints into the system C-space to form constraint obstacles using a recursive subdivision process. To generate motions along which actuator efforts can be specified without violating system constraints, paths are planned that avoid these constraint obstacles. The method permits the shape of the paths to be controlled using any configuration-dependent performance criterion. The FW approach was applied to a proposed three-limb planar climbing robot whose task is to climb upwards between two vertical walls by pushing outwards to generate frictional support. Motions were planned automatically within the system FW, enabling it to lift itself upwards using two limbs at a time, and a gait was planned to enable it to switch limbs and climb continuously.

Madhani, Akhil↗
NASA NTRS1993

Real-time collision avoidance in teleoperated whole-sensitive robot arm manipulators

A hybrid robot teleoperation system is presented which makes use of the methodology of motion planning for whole-sensitive robots to assist the operator in generating collision-free motion in a master-slave robot arm manipulator system. The system combines operator commands with data from the sensitive skin to guarantee safe motion for the entire body of the robot arm. The arm avoids obstacles automatically and in real time and moves in a collision-free manner although no prior knowledge of the objects in the environment is available to the motion planning system and no constraints are imposed on the obstacle shapes. The operator is thus relieved of the task of providing safety of the robot arm and surrounding objects.

Lumelsky, Vladimir J.↗