Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “optimization methods”

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 847 records · Page 47

Probabilistic Unsteady Aerodynamic Analysis

Probabilistic CFD design is needed because we are asked to do more with less. To cost effectively accomplish the design task, we need to formally quantify the effect of uncertainties (variables) in the design. Probabilistic design is one effective method to formally quantify the effect of uncertainties. Our objective is to establish a revolutionary new early design process, by developing non-deterministic physics-based probabilistic design tools, which will include all the life cycle processes. This work was concerned with the usefulness of parametric optimization method coupled with a Navier-Stokes analysis code for the aero-thermodynamic design of turbomachinery combustor liner. The interconnection between the CFD code and NESSUS codes will facilitate the coupling between the thermal profiles and structural design. We have developed new concepts for reducing the computational cost of unsteady, three-dimensional, compressible aerodynamic analyses for multistage turbomachinery flows. The flow was modeled by the three-dimensional Favre-Reynolds-averaged Navier-Stokes equations using the k-E turbulence closure, which was integrated using an implicit third-order upwind solver. The methodology developed in this work is expected to lead to the design optimization of turbomachinery blades.

Gorla, Rama S. R.↗

On making things the best - Aeronautical uses of optimization /Wright Bros. lecture/

The paper's purpose is to summarize and evaluate the results of an investigation into the degree to which formal optimization methods have contributed practically to the design and operation of atmospheric flight vehicles. The nature of this technology is reviewed and illustrated with simple structural examples. A series of published successful applications is described, from the fields of aerodynamics, structures, guidance and control, optimal trajectories and vehicle configuration optimization. The corresponding improvements over conventional analysis are assessed. Speculations are offered as to why these tools have made such little headway toward acceptance by designers. The growing need for their use in the future is explained; they hold out an unparalleled opportunity for improved efficiencies.

Ashley, H.↗

Input-response space-filling designs

Traditional space-filling designs are a convenient way to explore throughout an input space of flexible dimension and have design points close to any region where future predictions might be of interest. In some applications, there may be a model connecting the input factors to the response(s), which provides an opportunity to consider the spacing not only in the input space but also in the response space. In this paper, we present an approach for leveraging current understanding of the relationship between inputs and responses to generate designs that allow the experimenter to flexibly balance the spacing in these two regions to find an appropriate design for the experimental goals. Applications where good spacing of the observed response values include calibration problems where the goal is to demonstrate the adequacy of the model across the range of the responses, sensitivity studies where the outputs from a submodel may be used as inputs for subsequent models, and inverse problems where the outputs of a process will be used in the inverse prediction for the unknown inputs. Here, we use the multi-objective optimization method of Pareto fronts to generate multiple non-dominated designs with different emphases on the input and response space-filling criteria from which the experimenter can choose. The methods are illustrated through several examples and a chemical engineering case study.

42 ENGINEERING↗

Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations

Additive Runge-Kutta (ARK) methods are investigated for application to the spatially discretized one- dimensional convection-diffusion-reaction (CDR) equations. Accuracy, stability, conservation, and dense-output are first considered for the general case when N different Runge-Kutta methods are grouped into a single composite method. Then, implicit-explicit, (N = 2), additive Runge-Kutta (ARK(sub 2)) methods from third- to fifth-order are presented that allow for integration of stiff terms by an L-stable, stiffly-accurate explicit, singly diagonally implicit Runge-Kutta (ESDIRK) method while the nonstiff terms are integrated with a traditional explicit Runge-Kutta method (ERK). Coupling error terms of the partitioned method are of equal order to those of the elemental methods. Derived ARK(sub 2) methods have vanishing stability functions for very large values of the stiff scaled eigenvalue, z['] yields -infinity, and retain high stability efficiency in the absence of stiffness, z['] yield 0. Extrapolation-type stage- value predictors are provided based on dense-output formulae. Optimized methods minimize both leading order ARK(sub 2) error terms and Butcher coefficient magnitudes as well as maximize conservation properties. Numerical tests of the new schemes on a CDR problem show negligible stiffness leakage and near classical order convergence rates. However, tests on three simple singular-perturbation problems reveal generally predictable order reduction. Error control is best managed with a PID-controller. While results for the fifth-order method are disappointing, both the new third- and fourth-order methods are at least as efficient as existing ARK(sub 2) methods.

Kennedy, Christopher A.↗

Explosively driven Richtmyer–Meshkov instability jet suppression and enhancement via coupling machine learning and additive manufacturing

The ability to control the behavior of fluid instabilities at material interfaces, such as the shock-driven Richtmyer–Meshkov instability, is a grand technological challenge with a broad number of applications ranging from inertial confinement fusion experiments to explosively driven shaped charges. In this work, we use a linear-geometry shaped charge as a means of studying methods for controlling material jetting that results from the Richtmyer–Meshkov instability. A shaped charge produces a high-velocity jet by focusing the energy from the detonation of high explosives. The interaction of the resulting detonation wave with a hollowed cavity lined with a thin metal layer produces the unstable jetting effect. By modifying the characteristics of the detonation wave prior to striking the lined cavity, the kinetic energy of the jet can be enhanced or reduced. Modifying the geometry of the liner material can also be used to alter jetting properties. We apply optimization methods to investigate several design parameterizations for both enhancing or suppressing the shaped-charge jet. This is accomplished using 2D and 3D hydrodynamic simulations to investigate the design space that we consider. We also apply new additive manufacturing methods for producing the shaped-charge assemblies, which allow for the experimental testing of complicated design geometries obtained through computational optimization. We present a direct comparison of our optimized designs with experimental results carried out at the High Explosives Application Facility at Lawrence Livermore National Laboratory.

36 MATERIALS SCIENCE↗

Accelerated Discovery of Large Electrostrains in BaTiO 3 ‐Based Piezoelectrics Using Active Learning

Abstract A key challenge in guiding experiments toward materials with desired properties is to effectively navigate the vast search space comprising the chemistry and structure of allowed compounds. Here, it is shown how the use of machine learning coupled to optimization methods can accelerate the discovery of new Pb‐free BaTiO 3 (BTO‐) based piezoelectrics with large electrostrains. By experimentally comparing several design strategies, it is shown that the approach balancing the trade‐off between exploration (using uncertainties) and exploitation (using only model predictions) gives the optimal criterion leading to the synthesis of the piezoelectric (Ba 0.84 Ca 0.16 )(Ti 0.90 Zr 0.07 Sn 0.03 )O 3 with the largest electrostrain of 0.23% in the BTO family. Using Landau theory and insights from density functional theory, it is uncovered that the observed large electrostrain is due to the presence of Sn, which allows for the ease of switching of tetragonal domains under an electric field.

Yuan, Ruihao↗

A promising intersection of excited‐state‐specific methods from quantum chemistry and quantum Monte Carlo

Abstract We present a discussion of recent progress in excited‐state‐specific quantum chemistry and quantum Monte Carlo alongside a demonstration of how a combination of methods from these two fields can offer reliably accurate excited state predictions across singly excited, doubly excited, and charge transfer states. Both of these fields have seen important advances supporting excited state simulation in recent years, including the introduction of more effective excited‐state‐specific optimization methods, improved handling of complicated wave function forms, and ways of explicitly balancing the quality of wave functions for ground and excited states. To emphasize the promise that exists at this intersection, we provide demonstrations using a combination of excited‐state‐specific complete active space self‐consistent field theory, selected configuration interaction, and state‐specific variance minimization. These demonstrations show that combining excited‐state‐specific quantum chemistry and variational Monte Carlo can be more reliably accurate than either equation of motion coupled cluster theory or multi‐reference perturbation theory, and that it can offer new clarity in cases where existing high‐level methods do not agree. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Software > Quantum Chemistry

Otis, Leon↗

Convolutional Neural Network Surrogate Models for the Mechanical Properties of Periodic Structures

The work in this paper describes neural network surrogate models for calculating the effective mechanical properties of a periodic composites. The models achieve good accuracy even when only provided with training data sampling a small portion of the design space. As an example, the surrogate models are applied to solving the inverse design problem of finding structures with optimal mechanical properties. The surrogate models are sufficiently accurate to recover optimal solutions in general agreement with established topology optimization methods. However, improvements will be required to develop robust, efficient neural network-based surrogate models and several directions for future research are highlighted here.

42 ENGINEERING↗

Hormone purification by isoelectric focusing in space

The objective of the program was the definition and development of optimal methods for electrophoretic separations in microgravity. The approach is based on a triad consisting of ground based experiments, mathematical modeling and experiments in microgravity. Zone electrophoresis is a rate process, where separation is achieved in uniform buffers on the basis of differences in electrophoretic mobilities. Optimization and modeling of continuous flow electrophoresis mainly concern the hydrodynamics of the flow process, including gravity dependent fluid convection due to density gradients and gravity independent electroosmosis. Optimization of focusing requires a more complex model describing the molecular transport processes involved in electrophoresis of interacting systems. Three different focusing instruments were designed, embodying novel principles of fluid stabilization. Fluid stability was achieved by: (1) flow streamlining by means of membrane elements in combination with rapid fluid recycling; (2) apparatus rotation in combination with said membrane elements; and (3) shear stress induced by rapid recycling through a narrow gap channel.

Bier, M.↗

Genomic clocks and evolutionary timescales

For decades, molecular clocks have helped to illuminate the evolutionary timescale of life, but now genomic data pose a challenge for time estimation methods. It is unclear how to integrate data from many genes, each potentially evolving under a different model of substitution and at a different rate. Current methods can be grouped by the way the data are handled (genes considered separately or combined into a 'supergene') and the way gene-specific rate models are applied (global versus local clock). There are advantages and disadvantages to each of these approaches, and the optimal method has not yet emerged. Fortunately, time estimates inferred using many genes or proteins have greater precision and appear to be robust to different approaches.

Evolution, Molecular↗

A global approach of the representativity concept: Application on a high-conversion light water reactor MOX lattice case

The development of new types of reactor and the increase in the safety specifications and requirements induce an enhancement in both nuclear data knowledge and a better understanding of the neutronic properties of the new systems. This enhancement is made possible using ad hoc critical mock-up experiments. The main difficulty is to design these experiments in order to obtain the most valuable information. Its quantification is usually made by using representativity and transposition concepts. These theories enable to extract some information about a quantity of interest (an integral parameter) on a configuration, but generally a posteriori. This paper presents a more global approach of this theory, with the idea of optimizing the representativity of a new experiment, and its transposition a priori, based on a multiparametric approach. Using a quadratic sum, we show the possibility to define a global representativity which permits to take into account several quantities of interest at the same time. The maximization of this factor gives information about all quantities of interest. An optimization method of this value in relation to technological parameters (over-clad diameter, atom concentration) is illustrated on a high-conversion light water reactor MOX lattice case. This example tackles the problematic of plutonium experiment for the plutonium aging and a solution through the optimization of both the over-clad and the plutonium content. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

SU-E-T-101: An Arc Based Treatment Technique for Proton PBS Beam Utilizing Overlaid Bragg Peak Spots

Purpose: When Bragg peaks from proton PBS beams are overlaid from a large number of directions, it creates a localized dose distribution with rapid falloff. By positioning these overlaid Bragg peak (OBP) spots strategically throughout the targeted volume, a uniform and conformal dose distribution can be achieved. When PBS spots from all the directions are delivered in an arc fashion, it becomes the Volume Modulated Proton Therapy (VMPT) technique which corresponds to the VMAT technique in photon therapy. Methods: PBS beam from an IBA machine and a 20cm radius cylindrical phantom were simulated with GATE/GEANT4. To position the OBP spot at the center of the phantom, single energy PBS beam from 10/36/360 equally-spaced directions around the phantom was used. To position the OBP spot at an off-center position, PBS beams with various energies were used. To form a uniform dose distribution in 1D/2D, OBP spots are positioned in arrays with spot intervals and weight optimized. Results: Dose profile for single energy OBP spot improved when the number of beam angles increased from 10 to 36, but remained same from 36 to 360 angles. Multi-energy OBP spot had a slightly non-symmetrical profile due to the different energy proton beams used. By optimizing each OBP spot’s weight, 11 OBP spots with 10mm interval size created a uniform dose distribution covering a length of 10cm. Similarly, 11×11 OBP spots in an array of 10mm interval size created a uniform dose distribution covering an area of 10×10cm{sup 2}. Conclusion: We demonstrated the technique of using OBP spots to deliver uniform dose in 1D/2D. This technique can be easily extended through the optimization of position and weight of OBP spots in 3D and becomes the arc based VMPT technique. Further studies are needed to explore the optimization methods and full potential of the VMPT technique.

60 APPLIED LIFE SCIENCES↗

Two Point Exponential Approximation Method for structural optimization of problems with frequency constraints

The point exponential approximation method was introduced by Fadel et al. (Fadel, 1990), and tested on structural optimization problems with stress and displacement constraints. The reports in earlier papers were promising, and the method, which consists of correcting Taylor series approximations using previous design history, is tested in this paper on optimization problems with frequency constraints. The aim of the research is to verify the robustness and speed of convergence of the two point exponential approximation method when highly non-linear constraints are used.

Fadel, G. M.↗

Hydrothermal solubility of Dy hydroxide as a function of pH and stability of Dy hydroxyl aqueous complexes from 25 to 250 °C

The rare earth elements (REE) have important applications in green energy technologies. The formation of mineral deposits in geologic systems commonly involves hydrothermal fluids which can mobilize the REE. However, the REE speciation is not well known as a function of pH. The thermodynamic properties of REE hydroxyl complexes used in geochemical models are based on the Helgeson-Kirkham-Flowers (HKF) equation of state parameters which were derived by extrapolation of low temperature experimental and estimated data. In this study, Dy hydroxide solubility experiments are combined with available literature data to improve these models from 25 to 250 °C and optimize the thermodynamic properties of Dy 3+ and Dy hydroxyl complexes using GEMSFITS. Batch-type solubility experiments were conducted from 150 to 250 °C and at saturated water vapor pressure in perchloric acid solutions with initial pH values of 2 to 5 in 0.5 pH unit increments. The measured solubility of Dy hydroxide is retrograde with temperature and decreases with pH. The logarithm of total dissolved Dy molality ranges from –2.3 to –5.3 at 150 °C (pH 4.7–5.5), from –2.4 to –5.6 at 200 °C (pH 3.9–5.1), and from –3.7 to –6.9 at 250 °C (pH of 3.4 and 5.0). The optimized standard partial molal Gibbs energies of formation (Δ f G° T ) derived for Dy 3+ and DyOH 2+ display a close to linear relationship with temperature, fitting with previous optimizations based on DyPO 4 solubility data in the literature. A comparison of the optimized ΔfG°T values for aqueous Dy species with predictions from available HKF parameters indicates significant differences ranging from +11 to –26 kJ/mol between 25 and 250 °C. The experimental fits are used to derive the Dy hydroxide solubility products (K s0 ) and formation constants for the hydrolysis of Dy (β n with n = 1 to 3; Dy 3+ + nOH – = DyOH n 3-n ) as a function of temperature. The optimization method presented yields accurate thermodynamic properties for the Dy 3+ aqua ions and the DyOH 2+ species at the acidic to mildly acidic pH studied whereas more experimental work is needed at near-neutral and alkaline conditions to better constrain the other hydroxyl complexes. Furthermore, the optimized thermodynamic data have a significant impact on geochemical modeling of the mobility and solubility of REE minerals in acidic hydrothermal fluids.

58 GEOSCIENCES↗

Robust Decentralized Secondary Frequency Control in Power Systems: Merits and Trade-Offs

Frequency restoration in power systems is conven- tionally performed by broadcasting a centralized signal to local controllers. As a result of the energy transition, technological advances, and the scientific interest in distributed control and optimization methods, a plethora of distributed frequency control strategies have been proposed recently that rely on communication amongst local controllers. In this paper we propose a fully decentralized leaky integral controller for frequency restoration that is derived from a classic lag element. We study steady-state, asymptotic optimality, nominal stability, input-to-state stability, noise rejection, transient performance, and robustness properties of this controller in closed loop with a nonlinear and multivariable power system model. Here, we demonstrate that the leaky integral controller can strike an acceptable trade-off between performance and robustness as well as between asymptotic disturbance rejection and transient convergence rate by tuning its DC gain and time constant. We compare our findings to conventional decentralized integral control and distributed- averaging-based integral control in theory and simulations.

24 POWER TRANSMISSION AND DISTRIBUTION↗

New frequency domain methods for system identification

This paper presents two new techniques for frequency domain identification of linear system parameters. The first technique uses the instrumental variables approach. The frequency domain formulation is shown to give a considerable insight into the selection of efficient and convergent instrumental variables. The new maximum likelihood formulation affords simplier numerical solution and provides a way to select parameter starting values in the gradient based optimization method.

Gupta, N. K.↗

Rapid and semi-analytical design and simulation of a toroidal magnet made with YBCO and MgB 2 superconductors

Recent advances in second-generation (YBCO) high-temperature superconducting wire could potentially enable the design of super high performance energy storage devices that combine the high energy density of chemical storage with the high power of superconducting magnetic storage. However, the high aspect ratio and the considerable filament size of these wires require the concomitant development of dedicated optimization methods that account for the critical current density in type-II superconductors. In this study, we report on the novel application and results of a CPU-efficient semianalytical computer code based on the Radia 3-D magnetostatics software package. Our algorithm is used to simulate and optimize the energy density of a superconducting magnetic energy storage device model, based on design constraints, such as overall size and number of coils. The rapid performance of the code is pivoted on analytical calculations of the magnetic field based on an efficient implementation of the Biot-Savart law for a large variety of 3-D “base” geometries in the Radia package. The significantly reduced CPU time and simple data input in conjunction with the consideration of realistic input variables, such as material-specific, temperature, and magnetic-field-dependent critical current densities, have enabled the Radia -based algorithm to outperform finite-element approaches in CPU time at the same accuracy levels. Comparative simulations of MgB 2 and YBCO-based devices are performed at 4.2 K, in order to ascertain the realistic efficiency of the design configurations.

2G↗

SU-E-T-618: Dosimetric Comparison of Manual and Beam Angle Optimization of Gantry Angles in IMRT for Cervical Cancer

Purpose: To evaluate the dosimetric characteristics of intensity-modulated radiotherapy (IMRT) treatment plan with beam angle optimization. Methods: Ten post-operation patients with cervical cancer were included in this analysis. Two IMRT plans using seven beams were designed in each patient. A standard coplanar equi-space beam angles were used in the first plan (plan 1), whereas the selection of beam angle was optimized by beam angle optimization algorithm in Varian Eclipse treatment planning system for the same number of beams in the second plan (plan 2). Two plans were designed for each patient with the same dose-volume constraints and prescription dose. All plans were normalized to the mean dose to PTV. The dose distribution in the target, the dose to the organs at risk and total MU were compared. Results: For conformity and homogeneity in PTV, no statistically differences were observed in the two plans. For the mean dose in bladder, plan 2 were significantly lower than plan 1(p<0.05). No statistically significant differences were observed between two plans for the mean doses in rectum, left and right femur heads. Compared with plan1, the average monitor units reduced 16% in plan 2. Conclusion: The IMRT plan based on beam angle optimization for cervical cancer could reduce the dose delivered to bladder and also reduce MU. Therefore there were some dosimetric advantages in the IMRT plan with beam angle optimization for cervical cancer.

60 APPLIED LIFE SCIENCES↗