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 361 records · Page 20

Introducing Tropical Geometric Approaches to Delay Tolerant Networking Optimization

Delay Tolerant Networking (DTN) is the standard approach to the networking of space systems with the goal of supporting the Solar System Internet (SSI). Current space networks have a small scale and often depend on rigorously scheduled (pre-determined) contact opportunities; this manual approach inhibits scalability. The goal of this paper is to recast these scheduling problems in order to apply the optimization machinery of tropical geometry. Contact opportunities in space are dependent on such factors as orbital mechanics and asset availability, which induce time-varying connectivity; indeed, end-to-end connectivity might never occur. Routing optimization within this structure is classically difficult and typically utilizes Dijkstra's algorithm as applied to contact graphs. Alternatively, we follow the successes of tropical geometry in train schedule optimization, job assignments, and even traditional networking, by extending this approach to this more general (i.e. disconnected) problem space. These successes imply tropical geometry provides a useful framework in the context of DTNs, starting with applications to queuing theory and long-haul links. Recently, tropical geometry has been applied to parametric path optimization on graphs with variable edge weights. In this work, we extend these advances to account for the problem of routing in a space network, and find that tropical geometry is well-suited to the challenges offered by this new setting, including contact schedules featuring probabilities. Our approach leverages the combinatorial nature of the problem to give feasible shortest path trees in the presence of variable channel conditions and latency, evolving topologies, and uncertainty inherent in space routing. We discuss our tropical approach to DTN for two Python implementations, a Verilog Tropical ALU implementation, tropical frameworks for other parametric graph problems, and solution stability. Lastly, a program for future work is included to illuminate the path ahead.

Delay Tolerant Networking↗

Reentry-Vehicle Shape Optimization Using a Cartesian Adjoint Method and CAD Geometry

A DJOINT solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (e.g., geometric parameters that control the shape). Classic aerodynamic applications of gradient-based optimization include the design of cruise configurations for transonic and supersonic flow, as well as the design of high-lift systems. are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric computer-aided design (CAD). In previous work on Cartesian adjoint solvers, Melvin et al. developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the two-dimensional Euler equations using a ghost-cell method to enforce the wall boundary conditions. In Refs. 18 and 19, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm were the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The accuracy of the gradient computation was verified using several three-dimensional test cases, which included design variables such as the free stream parameters and the planform shape of an isolated wing. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Factors under consideration include the computation of mesh sensitivities that provide a reliable approximation of the objective function gradient, as well as the computation of surface shape sensitivities based on a direct-CAD interface. We present detailed gradient verification studies and then focus on a shape optimization problem for an Apollo-like reentry vehicle. The goal of the optimization is to enhance the lift-to-drag ratio of the capsule by modifying the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design.

Nemec, Marian↗

Melt Pool and Heat Treatment Optimization for the Fabrication of High-Strength and High-Toughness Additively Manufactured 4340 Steel

Additively manufactured (AM) components offer superior design flexibility compared to their conventionally manufactured counterparts, and optimizing processing parameters is key to achieving high-quality depositions with desirable and predictable mechanical properties. This study was focused on 4340 steel fabricated using laser powder bed fusion (LPBF), and 42 laser power and scan speed combinations have been systematically investigated to determine an optimized melt pool geometry that would ensure fully-dense parts. The AM material was compared with a wrought 4340 equivalent and studied in two customized heat treated conditions, optimized for strength and toughness, respectively. The microstructures of the as-fabricated and heat treated AM and wrought materials were characterized to assess differences introduced by the layer-by-layer fabrication process and subsequent heat treatment. Tensile properties of both materials were also evaluated and demonstrate that the AM materials offer equal or superior properties compared to the wrought equivalents. Differences in fracture surface morphologies indicate the distinct failure mechanisms associated with the materials’ characteristic microstructures, and the role of inclusions in the failures was studied to elucidate these differences. Complementary to the experimental investigations, the dataset was leveraged to make recommendations for future design of experiments to optimize AM build parameters in other material systems. A statistical Monte Carlo analysis was used to predict the interpolation error produced using reduced datasets and to enable informed processing parameters selection. These findings are discussed to make recommendations for the use of AM materials for high-integrity structural applications.

36 MATERIALS SCIENCE↗

Hexagonal Geometries in MPACT

The MPACT code is a high-fidelity light-water reactor analysis code using whole-core pin-resolved neutron transport calculations on modern parallel-computing hardware. MPACT uses the 2D/1D method to solve 3D neutron transport problems by decomposing the problem into a stack of 2D slices, each of which is solved independently using the method of characteristics (MOC). The slices are then coupled axially using the P3 nodal expansion method (NEM-P3) for the 1D axial calculations. MPACT also employs the coarse mesh finite difference (CMFD) method to accelerate calculations. This manuscript details work supporting advanced reactor designs using hexagonal pins and hexagonal assemblies such as the VVER-1000. If performed correctly, MOC is geometry agnostic. However, MPACT previously had optimizations in place for Cartesian geometries, specifically in the modularization and current calculations. Sections 2 and 3 detail the changes made to MPACT to support MOC and CMFD calculations on hexagonal geometries. Section 4 reports results demonstrating solution consistency for problems run with and without CMFD acceleration, results demonstrating solution consistency when run in serial and parallel, and pincell results using the Monte Carlo code, McCard’s benchmark results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Coupled Target-Beam-Moderator Optimization for the Second Target Station

This report describes the results for a coupled target-beam-moderator optimization analysis for the Second Target Station (STS) at ORNL's Spallation Neutron Source. This study is a continuation of the optimization analysis for the moderators in the preliminary design of STS performed in 2022. In the 2022 analysis the dimensions of the moderators are parameterized, while the target and the proton beam profile are kept constant. In this analysis the target height and the proton beam profile are added as parameters. This allows to study the coupled effects of changing target, moderator and beam dimensions. Similar to the 2022 analysis, this work is performed with an automated optimization workflow that uses the optimization toolbox DAKOTA, parameterized geometries in CREO and SpaceClaim, the unstructured mesh generation in Attila4MC, and the particle transport code MCNP6.2©. This workflow enables an efficient optimization using high-fidelity geometries. The main conclusions of this analysis are the following: • Coupled beam-target-moderator optimization provides a few additional percent performance gain over stand-alone moderator optimization. • The moderator performance is not very sensitive to the target height (between ≈60 and ≈80 mm) as long as the beam profile is chosen adequately. • The moderator performance is sensitive to the choice of beam spatial standard deviations, even when the footprint is kept constant. • The optimal moderator radius is the same for a beam footprint of 30 cm 2 , 62.5 cm 2 , and 90 cm 2 . Also the slope of the super-gaussian beam profile does not significantly impact the optimal moderator radius. • The optimal parameters and sensitivities are very similar to the 2022 optimization analysis. These results only indicate a a difference in the optimal radius of the cylindrical moderator, however, this has been corrected in the final design moderator optimization. The main purpose of this report is to document the simulations, results and lessons learned. The most impactful results are summarized in. We also note that the target geometry used in this work is not the final design.

43 PARTICLE ACCELERATORS↗

Innovations in optimization and control of accelerators using methods of differential geometry and genetic algorithms (Final Report)

Online tuning of particle accelerators is necessary in order to achieve optimal machine performance. However, it is also a major challenge due to the large parameter space which must be searched and the fact that many of the desir- able objectives compete with one another, and so will not reach their optimal values simultaneously. In order to mitigate these issues, we have explored using dimension-reduction techniques to reduce the size of the parameter space which must be searched and multi-objective genetic algorithms to obtain the sets of tuning parameters which provide pareto-optimal values for the objectives. These methods have enabled us to obtain improved values of the vertical emittance at the Cornell Electron Storage Ring (CESR), and to do so with greater control of orbit errors. Algorithmic tuning is a multidisciplinary endeavour, requiring expertise in beam dynamics, diagnostics, control systems and computer science, and thus a key practical problem is to formulate a common language in which experts with different specialties can communicate. We developed a solution to this problem in the form of a generic accelerator software interface that allows for rapid prototyping of optimization and control algorithms. Our interface is built on the Experimental Physics and Industrial Control Systems (EPICS) and makes possible testing control code in simulation before deployment on real accelerators, as well as deployment of third-party optimization code.

43 PARTICLE ACCELERATORS↗

On the Use of Parmetric-CAD Systems and Cartesian Methods for Aerodynamic Design

Automated, high-fidelity tools for aerodynamic design face critical issues in attempting to optimize real-life geometry arid in permitting radical design changes. Success in these areas promises not only significantly shorter design- cycle times, but also superior and unconventional designs. To address these issues, we investigate the use of a parmetric-CAD system in conjunction with an embedded-boundary Cartesian method. Our goal is to combine the modeling capabilities of feature-based CAD with the robustness and flexibility of component-based Cartesian volume-mesh generation for complex geometry problems. We present the development of an automated optimization frame-work with a focus on the deployment of such a CAD-based design approach in a heterogeneous parallel computing environment.

Nemec, Marian↗

Emission from dust in visual reflection nebulae at infrared and submillimeter wavelengths

Far-IR and submillimeter images of five bright visual reflection nebulae, IC 446, NGC2247, NGC 2245, NGC 7023, and CED 201 are presented and used with composite IR spectra to derive parameters such as the fraction of nebulae emission attributed to molecule-sized grains, the range of nebulae grain albedos, gas densities, and gas cloud geometries. The results show that 30-45 percent of the nebulae emission lies at wavelengths of less than 30 microns. The variation in IR luminosity may be related to variations in nebulae gas density and less than optimal gas cloud geometries rather than to anomalous grain albedos. Relative extinction efficincies of roughly 1000-5000 are inferred for stellar photospheres with effective wavelengths of roughly 0.25-0.5 micron. The resultant mass-extinction coefficient is roughly 10-50 sq cm/g.

Casey, Sean C.↗

CAD-Based Aerodynamic Design of Complex Configurations using a Cartesian Method

A modular framework for aerodynamic optimization of complex geometries is developed. By working directly with a parametric CAD system, complex-geometry models are modified nnd tessellated in an automatic fashion. The use of a component-based Cartesian method significantly reduces the demands on the CAD system, and also provides for robust and efficient flowfield analysis. The optimization is controlled using either a genetic or quasi-Newton algorithm. Parallel efficiency of the framework is maintained even when subject to limited CAD resources by dynamically re-allocating the processors of the flow solver. Overall, the resulting framework can explore designs incorporating large shape modifications and changes in topology.

Nemec, Marian↗

Rational Design of Nanoplasmonic Array Geometries for Biosensing

Background: Molecular diagnostics provide early and accurate diagnosis, which is essential for the prevention and treatment of infectious as well as chronic diseases. These tests are designed to detect disease-specific bioanalytes such as nucleic acid (DNA or RNA) or protein (antigens, antibodies) biomarkers. In the context of infectious disease diagnosis, nucleic acid-based detection methods are known to provide more specific and sensitive results. Here, the presence of a unique sequence belonging to the pathogenic genomic material is targeted to identify species, organism, genera and/or antimicrobial resistant gene markers. The majority of the common nucleic acid based diagnostic techniques require amplification (polymerase chain reaction, isothermal amplification etc.) of the pathogenic genetic material prior to detection impacting diagnostic speed, complexity, and cost thereby limiting ease of use. Thus, the development of simplified nucleic acid-based diagnostics that can be even used in resource-poor settings may hugely benefit patients across the globe. Nanopath is a molecular diagnostics company utilizing a solid-state nanosensor to enable sequence-specific detection of target nucleic acids without the need of amplification. These nanostructures enable ultra-sensitive biomarker detection using geometric, feature-dependent properties highly dependent on the local dielectric environment, allowing them to be sensitive to low concentration binding events. This paper describes an application of this approach to provide highly relevant clinical information within a single doctor’s office visit. Intro: The Nanopath team is in collaboration with NASA (National Aeronautics and Space Administration) and NIST (National Institute of Standards and Technology) to push the bounds of the fundamental physics associated with their biosensing platform. The ability of metals to support electromagnetic surface waves gives rise to surface plasmons when optically illuminated. This property, and its strong sensitivity to changes in the local refractive index, allows for the use of metal nanoparticles as ultra-sensitive transducers. In prior work by members of this team, ensembles of randomly oriented nanoparticles (i.e., colloidal nanorods dispersed on chip) were employed for sequence-specific nucleic acid sensing (1-3). While these particle sensors have the advantage of rapid fabrication, they suffer from low sensitivity and quality factor due to the random particle dispersity. In contrast, in this study we employ ordered array nanoparticle ensembles which can be used to improve sensor sensitivity and figure-of-merit. Study Methods Overview: In this talk, we detail the results of sensing experiments and computational simulations to outline a rational design of the structure of these plasmonic nanoparticle arrays for biomolecular sensing. Through simulation and experiment, we iteratively tailor nanostructure dimension to provide high quality signal and large resonance shifts upon modeled nucleic acid binding. In particular, full-wave electromagnetic simulations were conducted using Lumerical photonic simulation software in which periodic boundary conditions were applied in the x- and y- dimensions for each of the nanoplasmonic sensor geometries. To simulate the resonance response to changes in the bulk solution in contact with the sensor surface, the refractive index of the surrounding media was changed appropriately. Nucleic acid hybridization events were modeled using either using spherical structures approximating the relevant radius of genomic material as estimated by polymer models, or as conformal layers with the known refractive indices for nucleic acids. On the basis of initial simulations, nanosensors were fabricated using traditional electron-beam lithography protocols at NIST. To evaluate consensus between simulations and experiments, bulk sensing experiments were carried out in which the resonance peaks were obtained by submerging the sensors in refractive index standards. Key nanosensor characteristics including resonance peak locations, resonance peak shifts as a function of refractive index, and figure of merit (FOM) of extinction curves were examined between the experimental and simulation results prior to proceeding with simulations on additional geometries and more complex solution conditions, and further device fabrication. This iterative process is repeated toward a rational design of nanoplasmonic array geometries for biosensing optimizing response for targeted disease detection. In summary, this study puts forth a methodology for rational design and characterization of regularly spaced nanoparticle arrays for optics-based biosensing. The results of this study will allow for more informed design of nanostructure geometries towards sequence-specific nucleic acid detection. These improved designs have the potential to improve clinical sensitivity and limit-of-detection across disease indication.

sensor↗

Wing design code using three-dimensional Euler equations and optimization

This paper describes a new wing design code which is based on the Euler equations and a constrained numerical optimization technique. The geometry modification is based on a set of fundamental modes define on the unit interval. A design example involving a high speed civil transport wing is presented to demonstrate the usefulness of the design code. It is shown that the use of an Euler solver in the direct numerical optimization procedures is affordable on the current generation of supercomputers.

Chang, I-Chung↗

Optimization of Wing-Body Configurations by the Euler Equations

This paper describes a new wing-body design procedure which is based on the Euler equations and a constrained numerical optimization technique. The geometry modification is based on a set of fundamental modes defined on the unit interval. A design example involving a generic wing-body model is presented to demonstrate the usefulness of the design program. It is shown that the use of an Euler solver coupled with a direct numerical optimization procedure is affordable on the current generation of supercomputers.

Chang, I.-Chung↗

Displacement Energies for Nanorods

The geometries have been optimized and harmonic frequencies computed for three small hydrocarbon rods. The two lowest unique frequencies correspond to a torsion and a bend of the rods. The energy to bend or twist the rods has been computed by adding a fraction of the normal coordinate corresponding to these two modes to the equilibrium geometry. The density functional theory (using B3LYP) and self-consistent field (SCF) results are very similar for all properties considered.

Bauschlicher, Charles W., Jr.↗

Enhancing Header Shape Through Computational Fluid Dynamics for Improved Performance

Shape optimization in power plant design is crucial for maximizing efficiency and minimizing energy losses It impacts performance, cost effectiveness, and environmental sustainability Our project focuses on optimizing header pipe geometry using a method that considers temperature, flow, and pressure distributions, along with structural analysis This approach ensures structural integrity while minimizing material costs

20 FOSSIL-FUELED POWER PLANTS↗

RDE Nozzle Computational Design Methodology Development and Application

The Open National Combustion Code (OpenNCC) is used to simulate the aerospike exhaust nozzle region of a rotating detonation engine (RDE).The main objectives of the current study are twofold. The first goal is to validate the proposed computational methodology using the experimental data. The second goal is to demonstrate how the validated prediction tool can be used to optimize the nozzle geometry. To achieve a significant speed-up of computational time, the computational domain in divided into two parts: a combustion region and a throat-nozzle section. In the combustion region, a validated quasi-two-dimensional in-house code is utilized to generate an unsteady RDE flow field solution just upstream of the combustor throat. Subsequently, the unsteady flow data is fed into the three-dimensional throat-nozzle section as the inflow boundary condition. This facilitates the design optimization process since the unsteady inflow can be reused, and a relatively coarse mesh (i.e., larger time-step) can be used to analyze the flow fields around the nozzle. Five nozzle designs were studied and the predicted performance (i.e., thrust)compared. The optimized nozzle was found to produce 3.2% more overall thrust than a baseline nozzle design. Without the nozzle, there is a large low-pressure region at the throat exit, which significantly reduces the overall performance. This methodology is shown to be a promising approach to explore a wide variety of nozzle geometries in a relatively short amount of time.

CFD↗

Optimization of three-dimensional metamaterials for terahertz energy harvesting

This work focuses on finite element modeling (FEM) of a three-dimensional metamaterial used as an absorber for terahertz energy harvesting. The metamaterial consists of patterned pillars of an SU-8 dielectric photoresist coupled to a copper metal overlayer. Here, our study shows that the electromagnetic performance of the metamaterial is dependent on the following characteristic design parameters of the SU-8 dielectric: pillar height, bottom side length, and spacing between adjacent pillars. Using FEM, the metamaterial geometry is successfully optimized and the surface plasmon can be tuned to a peak frequency of 1.2 THz and a maximum terahertz absorption amplitude of 30%.

36 MATERIALS SCIENCE↗

Aerodynamic configuration design using response surface methodology analysis

An investigation has been conducted to determine a set of optimal design parameters for a single-stage-to-orbit reentry vehicle. Several configuration geometry parameters which had a large impact on the entry vehicle flying characteristics were selected as design variables: the fuselage fineness ratio, the nose to body length ratio, the nose camber value, the wing planform area scale factor, and the wing location. The optimal geometry parameter values were chosen using a response surface methodology (RSM) technique which allowed for a minimum dry weight configuration design that met a set of aerodynamic performance constraints on the landing speed, and on the subsonic, supersonic, and hypersonic trim and stability levels. The RSM technique utilized, specifically the central composite design method, is presented, along with the general vehicle conceptual design process. Results are presented for an optimized configuration along with several design trade cases.

Engelund, Walter C.↗

Analysis of Conduction Cooling Strategies for Wire Arc Additive Manufacturing

Metal additive manufacturing (AM) processing consists of numerous parameters which take time to optimize for various geometries. One aspect of the metal AM process that continues to be explored is the control of thermal energy accumulation during component manufacturing due to the melting and solidification of the feedstock. Excessive energy accumulation causes thermal failure of the component while minimal energy accumulation causes lack of fusion with the build plate or previous layer. The ability to simulate the thermal response of an AM component can increase research efficiency by reducing the time to optimize thermal energy accumulation. This paper presents an effective implementation of finite element analysis to determine the thermal response of a wire arc additive manufactured component with various build plate sizes and cooling methods including, integral build plate cooling, oversized build plates with passive cooling, and non-integral build plate cooling. The use of integral build plate cooling channels was shown to decrease the interpass temperature at the conclusion of the build process by 55% and build plate temperature by 96% compared to the conventionally deposited sample with 20 second dwell time. The use of a tall build plate with passive cooling was shown to reduce the interpass temperature by 32% as compared to the conventionally deposited sample with 20 second dwell time. Each cooling strategy evaluated decreased the interpass temperature within a range of 20–55% which enables higher deposition rates and decreased dwell times during depositions. The cooling strategies are designed to be implemented in a hybrid or retrofit AM platform to mitigate concerns of the thermal input from the additive process having detrimental effects on the precision of the machining process. This paper shows that accurate simulations of all strategies can be used to accurately predict the thermal response of the various strategies discussed. These cooling strategies will allow for increased deposition rates with comparable interpass temperature and decreased dwell time, increasing deposition efficiency. This model and these simulations are verified by experimental results. It is concluded that passive strategies, such as the over-sized tall build plate, can be used when liquid coolant in the AM environment could negatively affect the deposition process. Active cooling strategies, such as the integral build plate cooling could be used if low thermal conductivity materials are deposited or higher material deposition rates are desired. This paper discusses the use of active and passive cooling used during AM and shows how a simulation model can be used to make design choices for cooling strategies. The model also enables verification of select critical process parameters such as dwell times for a desired interpass temperature.

Heinrich, Lauren↗