Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geometry optimization”

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,099 records · Page 61

Linear Approximation to Optimal Control Allocation for Rocket Nozzles with Elliptical Constraints

In this paper we present a straightforward technique for assessing and realizing the maximum control moment effectiveness for a launch vehicle with multiple constrained rocket nozzles, where elliptical deflection limits in gimbal axes are expressed as an ensemble of independent quadratic constraints. A direct method of determining an approximating ellipsoid that inscribes the set of attainable angular accelerations is derived. In the case of a parameterized linear generalized inverse, the geometry of the attainable set is computationally expensive to obtain but can be approximated to a high degree of accuracy with the proposed method. A linear inverse can then be optimized to maximize the volume of the true attainable set by maximizing the volume of the approximating ellipsoid. The use of a linear inverse does not preclude the use of linear methods for stability analysis and control design, preferred in practice for assessing the stability characteristics of the inertial and servoelastic coupling appearing in large boosters. The present techniques are demonstrated via application to the control allocation scheme for a concept heavy-lift launch vehicle.

Orr, Jeb S.↗

Fast ion confinement in quasi-axisymmetric stellarator equilibria

This report presents an initial analysis of the fast ion confinement and losses within quasi-axisymmetric stellarator equilibria in consideration by Thea Energy. The equilibria have not yet been explicitly optimized for fast particle confinement and require validation. Modeling with the ASCOT5 code is used to directly examine the fast ion transport. The particle tracking simulations are purely (neo)classical in nature and simply contain the supplied equilibrium and collisions (pitch-angle, energy slowing, and velocity diffusion) from supplied thermal profiles. Uniform marker deposition is used to probe the general confinement properties of the equilibria while a realistic beam-born population is provided from the BEAMS3D code and an alpha particle population is calculated from a fusion source integrator. A first wall is included and defines the loss boundary. Analysis for NBI ions within Thea Energy’s conceptual Eos neutron source are presented along with alpha particles in an enlarged DT-plasma. The fast ions are assessed in regards to their confinement time, pitch, energy, and spatial coordinates. For each population, the impact of collisions and orbit drifts are discussed. It is found that NBI born ions in Eos are strongly confined until slowing-down, owing largely to the tangential injection geometry, while 22% of DT-born alpha energy is lost in the scaled device, indicating that any fusion pilot plant design optimization should include metrics for fast ion confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation of Liquid Lithium Divertor Geometry Using SOLPS-ITER

Here, plasma-facing component (PFC) geometries are evaluated for a Fusion Nuclear Science Facility (FNSF) design to find solutions that are compatible with flowing liquid lithium components while satisfying requirements from the perspective of plasma and neutral particle transport. Flowing liquid metal (LM) divertor systems offer important advantages due to continual regeneration of the surface material, but may require modifications from standard optimized tokamak divertor designs. Scrape Off Layer Plasma Simulator for the International Thermonuclear Experimental Reactor (SOLPS-ITER) simulations are used to compare standard vertical target, open, and balanced baffled geometries. It is found that the open geometry offers little control of neutral particles and has narrow operating windows where the divertor fluxes can be lowered to acceptable levels, while maintaining the required upstream density. The addition of neutral baffles allows greater control over the divertor conditions with less impact on the upstream density. The divertor neutral pressure in the baffled geometry is slightly higher than the open divertor, but not as large with a vertical target. A first-pass coupling of the SOLPS solution with sheath and implantation/erosion models to compute the lithium emission shows that the lithium is largely confined to the near-surface region due to the small ionization mean free path and large main ion flow toward the surface. The lithium source was varied by three orders of magnitude. At the highest rate, when the sourced lithium approaches the level of the main ion recycling flux, a high core lithium concentration (~10%) is observed with only moderate power dissipation to affect the divertor plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spline-Based Smoothing of Airfoil Curvatures

Constrained fitting for airfoil curvature smoothing (CFACS) is a splinebased method of interpolating airfoil surface coordinates (and, concomitantly, airfoil thicknesses) between specified discrete design points so as to obtain smoothing of surface-curvature profiles in addition to basic smoothing of surfaces. CFACS was developed in recognition of the fact that the performance of a transonic airfoil is directly related to both the curvature profile and the smoothness of the airfoil surface. Older methods of interpolation of airfoil surfaces involve various compromises between smoothing of surfaces and exact fitting of surfaces to specified discrete design points. While some of the older methods take curvature profiles into account, they nevertheless sometimes yield unfavorable results, including curvature oscillations near end points and substantial deviations from desired leading-edge shapes. In CFACS as in most of the older methods, one seeks a compromise between smoothing and exact fitting. Unlike in the older methods, the airfoil surface is modified as little as possible from its original specified form and, instead, is smoothed in such a way that the curvature profile becomes a smooth fit of the curvature profile of the original airfoil specification. CFACS involves a combination of rigorous mathematical modeling and knowledge-based heuristics. Rigorous mathematical formulation provides assurance of removal of undesirable curvature oscillations with minimum modification of the airfoil geometry. Knowledge-based heuristics bridge the gap between theory and designers best practices. In CFACS, one of the measures of the deviation of an airfoil surface from smoothness is the sum of squares of the jumps in the third derivatives of a cubicspline interpolation of the airfoil data. This measure is incorporated into a formulation for minimizing an overall deviation- from-smoothness measure of the airfoil data within a specified fitting error tolerance. CFACS has been extensively tested on a number of supercritical airfoil data sets generated by inverse design and optimization computer programs. All of the smoothing results show that CFACS is able to generate unbiased smooth fits of curvature profiles, trading small modifications of geometry for increasing curvature smoothness by eliminating curvature oscillations and bumps (see figure).

Li, W.↗

Nonlinear zero-sum differential game analysis by singular perturbation methods

A class of nonlinear, zero-sum differential games, exhibiting time-scale separation properties, can be analyzed by singular-perturbation techniques. The merits of such an analysis, leading to an approximate game solution, as well as the 'well-posedness' of the formulation, are discussed. This approach is shown to be attractive for investigating pursuit-evasion problems; the original multidimensional differential game is decomposed to a 'simple pursuit' (free-stream) game and two independent (boundary-layer) optimal-control problems. Using multiple time-scale boundary-layer models results in a pair of uniformly valid zero-order composite feedback strategies. The dependence of suboptimal strategies on relative geometry and own-state measurements is demonstrated by a three dimensional, constant-speed example. For game analysis with realistic vehicle dynamics, the technique of forced singular perturbations and a variable modeling approach is proposed. Accuracy of the analysis is evaluated by comparison with the numerical solution of a time-optimal, variable-speed 'game of two cars' in the horizontal plane.

Sinar, J.↗

The role of service areas in the optimization of FSS orbital and frequency assignments

A relationship is derived, on a single-entry interference basis, for the minimum allowable spacing between two satellites as a function of electrical parameters and service-area geometries. For circular beams, universal curves relate the topocentric satellite spacing angle to the service-area separation angle measured at the satellite. The corresponding geocentric spacing depends only weakly on the mean longitude of the two satellites, and this is true also for alliptical antenna beams. As a consequence, if frequency channels are preassigned, the orbital assignment synthesis of a satellite system can be formulated as a mixed-integer programming (MIP) problem or approximated by a linear programming (LP) problem, with the interference protection requirements enforced by constraints while some linear function is optimized. Possible objective-function choices are discussed and explicit formulations are presented for the choice of the sum of the absolute deviations of the orbital locations from some prescribed ideal location set. A test problem is posed consisting of six service areas, each served by one satellite, all using elliptical antenna beams and the same frequency channels. Numerical results are given for the three ideal location prescriptions for both the MIP and LP formulations. The resulting scenarios also satisfy reasonable aggregate interference protection requirements.

Levis, C. A.↗

Aerodynamic Rotor Design for a 25 MW Offshore Downwind Turbine

Continuously increasing offshore wind turbine scales require rotor designs that maximize power and performance. Downwind rotors offer advantages in lower mass due to reduced potential for tower strike, and is especially true at large scales, e.g., for a 25 MW turbine. In this study, three 25 MW downwind rotors, each with different prescribed lift coefficient distributions were designed (chord, geometry, and twist) and compared to maximize power production at unprecedented scales and Reynolds numbers, including a new approach to optimize rotor tilt and coning based on aeroelastic effects. To achieve this objective the design process was focused on achieving high power coefficients, while maximizing swept area and minimizing blade mass. Maximizing swept area was achieved by prescribing pre-cone and shaft tilt angles to ensure the aeroelastic orientation when the blades point upwards was nearly vertical at nearly rated conditions. Maximizing the power coefficient was achieved by prescribing axial induction factor and lift coefficient distributions which were then used as inputs for an inverse rotor design tool. The resulting rotors were then simulated to compare performance and subsequently optimized for minimum rotor mass. To achieve these goals, a high Reynolds number design space was developed using computational predictions as well as new empirical correlations for flatback airfoil drag and maximum lift. Within this design space, three rotors of small, medium and large chords were considered for clean airfoil conditions (effects of premature transition were also considered but did not significantly modify the design space). The results indicated that the medium chord design provided the best performance, producing the highest power in Region 2 from simulations while resulting in the lowest rotor mass, both of which support minimum LCOE. The methodology developed herein can be used for the design of other extreme-scale (upwind and downwind) turbines.

downwind rotors↗

Advances in high-pressure and high-temperature heat exchangers: Innovations via additive manufacturing and their applications

This paper provides a comprehensive overview of technological developments and applications for high-pressure and high-temperature (HPHT) heat exchangers produced through additive manufacturing (AM) techniques. As the demand for high-efficiency energy systems has been rising, the role of heat exchangers for thermal management and transport is becoming increasingly crucial. In this context, AM has emerged as a promising manufacturing solution to address the challenges posed by higher operating temperatures and/or pressures. This paper delves into three primary areas: the exploration of HPHT regimes supported by conventional heat exchangers, the impact of AM processes on thermal transport devices, and the discussion focused on HPHT heat exchangers produced via AM. The study examines the traditional manufacturing processes, materials, structures and geometries, and limitations of the commonly used heat exchangers for HPHT applications. Furthermore, it highlights the potential of AM in optimizing HPHT heat exchanger design, by reviewing studies that enable uniform flow distribution, minimization of fouling, and thermo-hydraulic performance augmentation without compromising structural stability. While addressing the limitations of cost, size, and other constraints, this study offers insights into the possibilities of AM in revolutionizing HPHT heat exchanger technology. The goal has been to fill the gap in literature by summarizing key studies, findings, and limitations, presenting a detailed examination, and proposing new areas for research. In conclusion, this review contributes to a deeper understanding of AM’s role in advancing HPHT heat exchangers for future energy systems.

36 MATERIALS SCIENCE↗

Compositionally Tuning Electron Transfer from Photoexcited Core/Shell Quantum Dots via Cation Exchange

It is critical to find methods to control the thermodynamic driving force for photoexcited charge transfer from quantum dots (QDs) and explore how this affects charge transfer rates, since the efficiency of QD-based photovoltaic and photocatalysis technologies depends on both this rate and the associated energetic losses. In this work, we introduce a single-pot shell growth and Cu-catalyzed cation exchange method to synthesize Cd x Zn 1-x Se/Cd y Zn 1-y S QDs with tunable driving forces for electron transfer. Functionalizing them with two molecular electron acceptors—naphthalenediimide (NDI) and anthraquinone (AQ)—allowed us to probe nearly 1 eV of driving forces. For AQ, at lower driving forces, we find that higher Zn content results in a 130-fold increase of electron transfer rate constants. However, at higher driving forces electron transfer dynamics are unaltered. Here, the data are understood using an Auger-assisted electron transfer model and analyzed with computational work to determine approximate binding geometries of these electron acceptors. Our work provides a method to tune QD reducing power and produces useful metrics for optimizing QD charge transfer systems that maximize rates of electron transfer while minimizing energetic losses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advanced Compact Generation Module with Fish Safe Runner Technology

Hydropower is an important contributor of stable, load-leveling renewable energy to our national grid. Increasingly, the development of new hydropower facilities or the retrofit of existing ones hinges not only on minimizing costs but also environmental impact. In this report, the analysis and testing of a modular and scalable low-head hydropower generation design using Natel Energy’s Restoration Hydro Turbine is described. This design leverages many established industry approaches for compactness and efficiency while simultaneously allowing for safe downstream fish passage through the turbines themselves. This unique approach reduces overall hydropower facility costs and enables a simpler inclusive method of project design and operation. To assess this design, mechanical and fluid computational analyses were used to study and optimize key parameters. Passage tests of important migratory species (salmonids, American eel) were conducted through representative turbines. The unique propeller geometry of the fish-safe Restoration Hydro Turbine was subjected to detailed design and testing using advanced manufacturing composite techniques. Comprehensive module cost models were developed and assessed alongside hydraulic efficiency. The results of this project show promising and economical applications for downstream passage offish through Restoration Hydro Turbine modules.

13 HYDRO ENERGY↗

Printed Metal Molds for Mainstream Automotive Production

The main focus of this work was to evaluate additive manufacturing of metal molds, with significantly reduced cost and/or procurement lead time, for production of large composite components for mainstream (> 100k units/yr) automobiles. A two-component large tool with conformal heating channels was designed, optimized, printed, machined, and tested in production. The resulting analysis indicates that large scale AM is suitable for manufacturing complex geometry large scale metal molds.

36 MATERIALS SCIENCE↗

Advancing In-Situ Plasma Processing for SRF Cavities at Fermilab

This talk presents Fermilab’s experience with plasma processing applied to both elliptical and low-beta superconducting RF cavities. The development and implementation of the technique are compared across different geometries, with particular emphasis on constraints imposed by cryomodule assembly configurations. The talk also reports on systematic studies aimed at optimizing the processing recipe, including the exploration of various inert gas mixtures and oxygen concentrations. These results aim at providing guidance for tailoring plasma processing conditions to different cavity types and advancing plasma processing applicability in operational cryomodules.

Giaccone, Bianca [Fermilab] (ORCID:000000027275846↗

Advancing In-Situ Plasma Processing for SRF Cavities at Fermilab

This talk presents Fermilab’s experience with plasma processing applied to both elliptical and low-beta superconducting RF cavities. The development and implementation of the technique are compared across different geometries, with particular emphasis on constraints imposed by cryomodule assembly configurations. The talk also reports on systematic studies aimed at optimizing the processing recipe, including the exploration of various inert gas mixtures and oxygen concentrations. These results aim at providing guidance for tailoring plasma processing conditions to different cavity types and advancing plasma processing applicability in operational cryomodules.

Giaccone, Bianca [Fermilab] (ORCID:000000027275846↗

FASTER 3: A generalized-geometry Monte Carlo computer program for the transport of neutrons and gamma rays. Volume 2: Users manual

A description of the FASTER-III program for Monte Carlo Carlo calculation of photon and neutron transport in complex geometries is presented. Major revisions include the capability of calculating minimum weight shield configurations for primary and secondary radiation and optimal importance sampling parameters. The program description includes a users manual describing the preparation of input data cards, the printout from a sample problem including the data card images, definitions of Fortran variables, the program logic, and the control cards required to run on the IBM 7094, IBM 360, UNIVAC 1108 and CDC 6600 computers.

Jordan, T. M.↗

A function space approach to state and model error estimation for elliptic systems

An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal estimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed.

Rodriguez, G.↗

A function space approach to state and model error estimation for elliptic systems

An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal esimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed. While this paper focuses on theoretical developments, applications arising in the area of large structure static shape determination are contained in a closely related paper (Rodriguez and Scheid, 1982).

Rodriguez, G.↗

Numerical Modeling of a Vortex Stabilized Arcjet

Arcjet thrusters are being actively considered for use in Earth orbit maneuvering applications. Experimental studies are currently the chief means of determining an optimal thruster configuration. Earlier numerical studies have failed to include all of the effects found in typical arcjets including complex geometries, viscosity, and swirling flow. Arcjet geometries are large area ratio converging nozzles with centerbodies in the subsonic portion of the nozzle. The nozzle walls serve as the anode while the centerbody functions as the cathode. Viscous effects are important because the Reynolds number, based on the throat radius, is typically less than 1,000. Experimental studies have shown that a swirl or circumferential velocity component stabilizes a constricted arc. This dissertation describes the equations governing flow through a constricted arcjet thruster. An assumption that the flowfield is in local thermodynamic equilibrium leads to a single fluid plasma temperature model. An order of magnitude analysis reveals the governing fluid mechanics equations are uncoupled from the electromagnetic field equations. A numerical method is developed to solve the governing fluid mechanics equations, the Thin Layer Navier-Stokes equations. A coordinate transformation is employed in deriving the governing equations to simplify the application of boundary conditions in complex geometries. An axisymmetric formulation is employed to include the swirl velocity component as well as the axial and radial velocity components. The numerical method is an implicit finite-volume technique and allows for large time steps to reach a converged steady-state solution. The inviscid fluxes are flux-split, and Gauss-Seidel line relaxation is used to accelerate convergence. Converging-diverging nozzles with exit-to-throat area ratios up to 100:1 and annular nozzles were examined. Quantities examined included Mach number and static wall pressure distributions, and oblique shock structures. As the level of swirl and viscosity in the flowfield increased the mass flow rate and thrust decreased. The technique was used to predict the flow through a typical arcjet thruster geometry. Results indicate swirl and viscosity play an important role in the complex geometry of an arcjet.

Pawlas, Gary E.↗

The isolation limits of stochastic vibration

The vibration isolation problem is formulated as a 1D kinematic problem. The geometry of the stochastic wall trajectories arising from the stroke constraint is defined in terms of their significant extrema. An optimal control solution for the minimum acceleration return path determines a lower bound on platform mean square acceleration. This bound is expressed in terms of the probability density function on the significant maxima and the conditional fourth moment of the first passage time inverse. The first of these is found analytically while the second is found using a Monte Carlo simulation. The rms acceleration lower bound as a function of available space is then determined through numerical quadrature.

Knopse, C. R.↗