Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “optimal start”

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 199 records · Page 11

An optimization program based on the method of feasible directions: Theory and users guide

The theory and user instructions for an optimization code based on the method of feasible directions are presented. The code was written for wide distribution and ease of attachment to other simulation software. Although the theory of the method of feasible direction was developed in the 1960's, many considerations are involved in its actual implementation as a computer code. Included in the code are a number of features to improve robustness in optimization. The search direction is obtained by solving a quadratic program using an interior method based on Karmarkar's algorithm. The theory is discussed, focusing on the important and often overlooked role played by the various parameters guiding the iterations within the program. Also discussed is a robust approach for handling infeasible starting points. The code was validated by solving a variety of structural optimization test problems that have known solutions obtained by other optimization codes. It has been observed that this code is accurate and robust: it has solved a variety of problems from different starting points. However, the code is inefficient in that it takes considerable CPU time as compared with certain other available codes. Further work is required to improve its efficiency while retaining its robustness.

Ashok D. Belegundu↗

A Mixed Integer Efficient Global Optimization Algorithm with Multiple Infill Strategy - Applied to a Wing Topology Optimization Problem

With the advancement in high performance computing and numerical optimization techniques,engineering design optimization problems are becoming more complex, larger scale,higher fidelity, and computationally more demanding, requiring longer run times than ever before. There exists methodologies and techniques that can address some of these challenges but very few can address all, and most are limited in the extent that these concerns can be addressed. With the goal of addressing such challenging engineering problems, we developed anew optimization framework, named AMIEGO, that combines concepts from surrogate-based optimization approaches, gradient-based numerical methods, Partial Least Squares, evolutionary algorithms, and Branch-and-Bound, providing newer capabilities that were not previouslyperceived. However, the original version of this framework, in the process of adaptive samplingto explore and exploit the design space, finds only a single sample point per iteration. The efforthere builds upon this previously developed optimization framework to include multiple infillsampling capability that combines the concept of generalized expected improvement function,unsupervised learning, and multi-objective evolutionary technique. To demonstrate, AMIEGOwith the multiple infill capability (called AMIEGO-MIMOS) solves a series of increasingly difficultengineering design optimization problems. The results reveal the performance of the newapproach is problem dependent. When applied to a ten-bar truss problem, the newly proposedmultiple infill strategy consistently leads to a better design solutions when compared to theexisting CPTV method (implemented with the context of the AMIEGO framework). On theother hand, when applied to a mixed-integer high fidelity wing topology optimization problem- MIMOS, despite showing a steeper convergence at the start, eventually leads to an inferiorsolution as compared to CPTV approach. These results also reveal that a small number ofstarting points, in general, are sufficient to lead to a good overall solution.

Mixed-integer optimization↗

Conservation of Helium while Maintaining High System Purity

Recent helium shortages and helium price increases have lead to an increased emphasis being placed on conserving helium. The need to conserve helium must be balanced with need to maintain the high levels of purity necessary to prevent operational problems caused by contamination. Helium losses and contamination control are especially important for test stands that have cryogenic distribution systems operating continuously with frequent changeover of cryogenic temperature components that are being tested. This paper describes a mathematical model to estimate the quantity of helium lost and the purity of the helium after the pump and backfill procedure is complete. The process to determine the optimal time during pump down to cut off pumping and start backfilling is described. There is a tradeoff between trying to achieve the lowest possible pressure during pumping and the quantity of air leaking into the volume while pumping is occurring. An additional benefit of careful selection of pump and backfill parameters in conjunction with real-time pressure monitoring can reduce the labor and time required to complete a successful pump and backfill procedure. This paper is intended to be a tool for engineers to review their pump and backfill procedures and measured data to optimize helium losses, system purity, and labor required.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Conservation of Helium while Maintaining High System Purity

Recent helium shortages and helium price increases have lead to an increased emphasis being placed on conserving helium. The need to conserve helium must be balanced with need to maintain the high levels of purity necessary to prevent operational problems caused by contamination. Helium losses and contamination control are especially important for test stands that have cryogenic distribution systems operating continuously with frequent changeover of cryogenic temperature components that are being tested. This paper describes a mathematical model to estimate the quantity of helium lost and the purity of the helium after the pump and backfill procedure is complete. The process to determine the optimal time during pump down to cut off pumping and start backfilling is described. There is a tradeoff between trying to achieve the lowest possible pressure during pumping and the quantity of air leaking into the volume while pumping is occurring. An additional benefit of careful sel ection of pump and backfill parameters in conjunction with real-time pressure monitoring can reduce the labor and time required to complete a successful pump and backfill procedure. This paper is intended to be a tool for engineers to review their pump and backfill procedures and measured data to optimize helium losses, system purity, and labor required.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Experiments and modeling to evaluate global reaction kinetics of three-way catalyst light off for net-zero carbon fuels and selected fuel blends

Internal combustion engine exhaust composition can have a significant impact on the light-off characteristics of CO, NO, and organics in the three-way-catalyst during cold-start. The U.S. Department of Energy’s Co-Optimization of Fuels and Engines initiative was aimed at developing new high-performance fuels, which raises the question of how current three-way-catalysts will perform if the exhaust composition is significantly changed. This work is aimed at quantifying the performance of the catalyst under a variety of different fuel components and mixtures. To this end, a reduced-order model is developed wherein the set of organics examined are placed into molecular groups based on their reactivities and to model each organic within that group with a common kinetic expression. In conclusion, this strategy allows for a first approximation to the light-off characteristics of all organics that fall within that group and to make predictions for fuel blends.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Autonomous alloy composition optimization using molecular dynamics guided by a large language model

Here, we present an autonomous materials discovery framework that couples a large language model (LLM) with molecular dynamics (MD) simulations to optimize Fe–Cr–Mn alloy compositions for tensile strength. Starting from six distinct compositions, the LLM operated as an intelligent agent, iteratively proposing changes based on prior simulation results and constraints. Over 50 iterations per case, the LLM adaptively explored the composition space, identifying high-strength regions, not easily accessible by conventional methods. The highest strength, 18.7 GPa, was achieved with Fe 71 Cr 25 Mn 4 composition, identified from a Fe 75 Cr 20 Mn 5 starting point. The LLM autonomously adjusted its strategy in real time, demonstrating closed-loop decision-making using commodity hardware. This approach showcases the potential of LLMs as scientific co-pilots, capable of accelerating materials discovery and generalizable to other domains like biology and drug design.

Autonomy↗

How Many Trip Requests Could We Support? An Activity-Travel Based Vehicle Scheduling Approach

In a world of ever-changing travel behavior and ever-increasing modal options, is vital to have integrated models that could capture the interactions between supply and demand layers of travel. Addressing this need, we propose three different versions of network representation and mathematical models for the activity-based vehicle routing problem to connect activity-travel graphs of passengers (demand layer) to spatio-temporal networks of vehicles (supply layer). Versions I and II are arc-based, while version III is path-based. In version I, we introduce the concept of activity-travel graphs for passengers. For vehicles, we construct space–time networks and add a new dimension, called “under-service state”, to track the execution status of trip requests at any location and time. In version II, we reduce the complexity of the network structure by eliminating the state dimension and some other modifications in the structure of the passengers’ and vehicles’ network. Although both versions can capture various behavioral constraints of the activity-based vehicle routing problem (e.g., mandatory and optimal activities, duration of activities, chain of activities, preferred starting and ending times of activities), due to the high level of complexity of the network structure, both versions can only solve small-sized problems. To tackle the computational complexity, we propose a path-based network representation in version III, and to make a balance between the disutility of passengers and vehicles, we present a tolled user equilibrium problem. Mathematical models are coded in C and GAMS and implemented on real-world Phoenix regional transportation network with more than 39 million trip requests, which demonstrate the effectiveness of the proposed solution for the original and restricted master problems.

33 ADVANCED PROPULSION SYSTEMS↗

Direct computation of magnetic surfaces in Boozer coordinates and coil optimization for quasisymmetry

We propose a new method to compute magnetic surfaces that are parametrized in Boozer coordinates for vacuum magnetic fields. We also propose a measure for quasisymmetry on the computed surfaces and use it to design coils that generate a magnetic field that is quasisymmetric on those surfaces. The rotational transform of the field and complexity measures for the coils are also controlled in the design problem. Using an adjoint approach, we are able to obtain analytic derivatives for this optimization problem, yielding an efficient gradient-based algorithm. Starting from an initial coil set that presents nested magnetic surfaces for a large fraction of the volume, our method converges rapidly to coil systems generating fields with excellent quasisymmetry and low particle losses. In particular for low complexity coils, we are able to significantly improve the performance compared with coils obtained from the standard two-stage approach, e.g. reduce losses of fusion-produced alpha particles born at half-radius from $17.7\,\%$ to $6.6\,\%$ . We also demonstrate 16-coil configurations with alpha loss ${<}1\,\%$ and neoclassical transport magnitude $\epsilon _{\text {eff}}^{3/2}$ less than approximately $5\times 10^{-9}$ .

Physics↗

De novo design of high-affinity binders of bioactive helical peptides

Many peptide hormones form an α-helix on binding their receptors, and sensitive methods for their detection could contribute to better clinical management of disease. De novo protein design can now generate binders with high affinity and specificity to structured proteins. However, the design of interactions between proteins and short peptides with helical propensity is an unmet challenge. Here we describe parametric generation and deep learning-based methods for designing proteins to address this challenge. We show that by extending RFdiffusion to enable binder design to flexible targets, and to refining input structure models by successive noising and denoising (partial diffusion), picomolar-affinity binders can be generated to helical peptide targets by either refining designs generated with other methods, or completely de novo starting from random noise distributions without any subsequent experimental optimization. The RFdiffusion designs enable the enrichment and subsequent detection of parathyroid hormone and glucagon by mass spectrometry, and the construction of bioluminescence-based protein biosensors. The ability to design binders to conformationally variable targets, and to optimize by partial diffusion both natural and designed proteins, should be broadly useful.

59 BASIC BIOLOGICAL SCIENCES↗

Design and Characterization of Highly Diffusive Turbine Vanes Suitable for Transonic Rotating Detonation Combustors

In rotating detonation engines the turbine inlet conditions may be transonic with unprecedented unsteady fluctuations. To ensure an acceptable engine performance, the turbine passages must be suited to these conditions. This article focuses on designing and characterizing highly diffusive turbine vanes to operate at any inlet Mach number up to Mach 1. First, the effect of pressure loss on the starting limit is presented. Afterward, a multi-objective optimization with steady RANS simulations, including the endwall and 3D vane design is performed. Compared to previous research, significant reductions in pressure loss and stator-induced rotor forcing are obtained, with an extended operating range and preserving high flow turning. Finally, the influence of the inlet boundary layer thickness on the vane performance is evaluated, inducing remarkable increases in pressure loss and downstream pressure distortion. Employing an optimization with a thicker inlet boundary layer, specific endwall design recommendations are found, providing a notable improvement in both objective functions.

Grasa, Sergio↗

Optimization of cascade blade mistuning. II - Global optimum and numerical optimization

The values of the mistuning which yield the most stable eigenvectors are analytically determined, using the simplified equations of motion which were developed in Part I of this work. It is shown that random mistunings, if large enough, may lead to the maximal stability, whereas the alternate mistunings cannot. The problem of obtaining maximum stability for minimal mistuning is formulated, based on numerical optimization techniques. Several local minima are obtained using different starting mistuning vectors. The starting vectors which lead to the global minimum are identified. It is analytically shown that all minima appear in multiplicities which are equal to the number of compressor blades. The effect of mistuning on the flutter speed is studied using both an optimum mistuning vector and an alternate mistuning vector. Effects of mistunings in elastic axis locations are shown to have a negligible effect on the eigenvalues. Finally, it is shown that any general two-dimensional bending-torsion system can be reduced to an equivalent uncoupled torsional system.

Nissim, E.↗

Validation of a New Procedure for Impedance Eduction in Flow

A new impedance eduction procedure is validated by comparing the educed impedance spectrum to that of an older but well-tested eduction procedure. The older procedure requires the installation of a microphone array in the liner test section but the new procedure removes this requirement. A 12.7-mm stainless steel plate and a conventional liner consisting of a perforated plate bonded to a honeycomb core are tested. Test data is acquired from a grazing flow, impedance tube facility for a range of source frequencies and mean flow Mach numbers for which only plane waves are cut on. For the stainless steel plate, the educed admittance spectrum using the new procedure shows an improvement over that of the old procedure. This improvement shows up primarily in the educed conductance spectrum. Both eduction procedures show discrepancies in educed admittance in the mid-frequency range. Indications are that this discrepancy is triggered by an inconsistency between the measured eduction data (that contains boundary layer effects) and the two eduction models (for which the boundary layer is neglected). For the conventional liner, both eduction procedures are in very good agreement with each other. Small discrepancies occur for one or two frequencies in the mid-frequency range and for frequencies beyond the cut on frequency of higher-order duct modes. This discrepancy in the midfrequency range occurs because an automated optimizer is used to educe the impedance and the objective function used by the optimizer is extremely flat and therefore sensitive to initial starting values. The discrepancies at frequencies beyond the cut on frequency of higher order duct modes are due to the assumption of only plane waves in the impedance eduction model, although higher order modes are propagating in the impedance tube facility.

Watson, W. R.↗

Integrated Lens Antennas for Multi-Pixel Receivers

Future astrophysics and planetary experiments are expected to require large focal plane arrays with thousands of detectors. Feedhorns have excellent performance, but their mass, size, fabrication challenges, and expense become prohibitive for very large focal plane arrays. Most planar antenna designs produce broad beam patterns, and therefore require additional elements for efficient coupling to the telescope optics, such as substrate lenses or micromachined horns. An antenna array with integrated silicon microlenses that can be fabricated photolithographically effectively addresses these issues. This approach eliminates manual assembly of arrays of lenses and reduces assembly errors and tolerances. Moreover, an antenna array without metallic horns will reduce mass of any planetary instrument significantly. The design has a monolithic array of lens-coupled, leaky-wave antennas operating in the millimeter- and submillimeter-wave frequencies. Electromagnetic simulations show that the electromagnetic fields in such lens-coupled antennas are mostly confined in approximately 12 15 . This means that one needs to design a small-angle sector lens that is much easier to fabricate using standard lithographic techniques, instead of a full hyper-hemispherical lens. Moreover, this small-angle sector lens can be easily integrated with the antennas in an array for multi-pixel imager and receiver implementation. The leaky antenna is designed using double-slot irises and fed with TE10 waveguide mode. The lens implementation starts with a silicon substrate. Photoresist with appropriate thickness (optimized for the lens size) is spun on the substrate and then reflowed to get the desired lens structure. An antenna array integrated with individual lenses for higher directivity and excellent beam profile will go a long way in realizing multi-pixel arrays and imagers. This technology will enable a new generation of compact, low-mass, and highly efficient antenna arrays for use in multi-pixel receivers and imagers for future planetary and astronomical instruments. These antenna arrays can also be used in radars and imagers for contraband detection at stand-off distances. This will be enabling technology for future balloon-borne, smaller explorer class mission (SMEX), and other missions, and for a wide range of proposed planetary sounders and radars for planetary bodies.

Lee, Choonsup↗

Integrated data-driven and experimental approaches to accelerate lead optimization targeting SARS-CoV- 2 main protease

Identification of potential therapeutic candidates can be expedited by integrating computational modeling with domain aware machine learning (ML) approaches followed by experimental validation. Generative deep learning models have been recently developed that can generate thousands of new candidates, but their physiochemical properties are typically not optimized. Using our deep learning models and a scaffold as a starting point, we generated tens of thousands of compounds for SARS-CoV-2 M pro that preserve the core scaffold. Here we utilized and implemented several computational tools such as structural alert and toxicity analysis, high throughput virtual screening, ML-based 3D quantitative structure–activity relationships, multi-parameter optimization, and graph neural networks on libraries of generated candidates to predict biological activity and binding affinity a priori. From these collective computational results, eight promising candidates were identified and tested experimentally using Native Mass Spectrometry (MS) and FRET-based functional assays. Two compounds, with quinazoline-2-thiol and acetylpiperidine core moiety showed IC 50 values in the low micromolar range: 2.95±0.0017 µM and 3.41±0.0015 µM, respectively. The molecular dynamics simulations further highlight that binding of these compounds results in allosteric modulations in the chain B and the interface domains of the M pro . The key fragments from these top hits can be used as input for closed loop lead optimization in the integrated pipeline.

60 APPLIED LIFE SCIENCES↗

Optimization Through Multi-Fidelity Modeling

We present a novel method for optimizing parameter selection for simulations with an evaluation budget. We start with an existing method for building a multi-fidelity model out of many low-fidelity simulations and few high-fidelity simulations. We propose a novel method to simplify parameter selection without sacrificing performance. We verify these results and compare with existing literature. Next, we propose a novel algorithm which uses this difference model to suggest new points in the parameter design space to simulate. We add each point we simulate to the model to improve its quality for the next iteration. The algorithm trades off reducing the uncertainty of the existing model with optimization of the objective. The first is more useful when a large fraction of the computation budget remains. The second is more useful when a small fraction of the computation budget remains. Our method converges to the optimum by using a high-fidelity evaluation for just 16 of the 427 points. Our method is general enough to work if there is no low-fidelity model. Furthermore, it is agnostic to the underlying physics of the problem. Therefore, both the low-fidelity and high-fidelity models can be generated by any arbitrary function, including simulations and physical experiments.

97 MATHEMATICS AND COMPUTING↗

Computation of optimal feedback strategies for interception in a horizontal plane

The problem of minimum-time interception of a target moving in a horizontal plane is studied. The target may fly a trajectory that is known at the start of the interception or try to evade the intercepting aircraft. The interceptor's optimal motion is described independently of the target in terms of an extremal trajectory map. The latter is used to develop guidelines for suboptimal approximations to the extremals. A method of constructing the feedback solution by drawing isochrones (constant minimum-time loci) is developed. Sections of the feedback solution for the interception of a straight-flying target are presented, and the construction of an isochrone for the interception of an actively evading target is demonstrated.

Rajan, N.↗

Computational Role of Tunneling in a Programmable Quantum Annealer

Quantum tunneling is a phenomenon in which a quantum state tunnels through energy barriers above the energy of the state itself. Tunneling has been hypothesized as an advantageous physical resource for optimization. Here we present the first experimental evidence of a computational role of multiqubit quantum tunneling in the evolution of a programmable quantum annealer. We developed a theoretical model based on a NIBA Quantum Master Equation to describe the multi-qubit dissipative cotunneling effects under the complex noise characteristics of such quantum devices.We start by considering a computational primitive, the simplest non-convex optimization problem consisting of just one global and one local minimum. The quantum evolutions enable tunneling to the global minimum while the corresponding classical paths are trapped in a false minimum. In our study the non-convex potentials are realized by frustrated networks of qubit clusters with strong intra-cluster coupling. We show that the collective effect of the quantum environment is suppressed in the critical phase during the evolution where quantum tunneling decides the right path to solution. In a later stage dissipation facilitates the multiqubit cotunneling leading to the solution state. The predictions of the model accurately describe the experimental data from the D-WaveII quantum annealer at NASA Ames. In our computational primitive the temperature dependence of the probability of success in the quantum model is opposite to that of the classical paths with thermal hopping. Specially, we provide an analysis of an optimization problem with sixteen qubits,demonstrating eight qubit cotunneling that increases success probabilities. Furthermore, we report results for larger problems with up to 200 qubits that contain the primitive as subproblems.

hard problems↗

Downward quantum learning from element 118: Automated generation of Fermi–Löwdin orbitals for all atoms

A new algorithm based on a rigorous theorem and quantum data computationally mined from element 118 guarantees automated construction of initial Fermi–Löwdin-Orbital (FLO) starting points for all elements in the Periodic Table. It defines a means for constructing a small library of scalable FLOs for universal use in molecular and solid-state calculations. The method can be systematically improved for greater efficiency and for applications to excited states such as x-ray excitations and optically silent excitations. FLOs were introduced to recast the Perdew–Zunger self-interaction correction (PZSIC) into an explicit unitarily invariant form. The FLOs are generated from a set of N quasi-classical electron positions, referred to as Fermi-Orbital descriptors (FODs), and a set of N-orthonormal single-electron orbitals. FOD positions, when optimized, minimize the PZSIC total energy. However, creating sets of starting FODs that lead to a positive definite Fermi orbital overlap matrix has proven to be challenging for systems composed of open-shell atoms and ions. The proof herein guarantees the existence of a FLOSIC solution and further guarantees that if a solution for N electrons is found, it can be used to generate a minimum of N – 1 and a maximum of 2 N – 2 initial starting points for systems composed of a smaller number of electrons. As a result, applications to heavy and super-heavy atoms are presented. All starting solutions reported here were obtained from a solution for element 118, Oganesson.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗