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 55 records · Page 3

Optimization of a High Pressure Turbine Blade and Sector-Based Annular Rig Design for Supercritical CO2 Power Cycle Representative Testing

Abstract As part of the ongoing research into the design of hardware for zero emission cycles, a first-stage high-pressure turbine (HPT) blade is optimized for a 300 MWe supercritical CO2 (sCO2) power cycle using the surrogate-assisted genetic algorithm optimizer in Numeca FINE/Design three-dimensional with objectives of increasing efficiency and decreasing heat load to the blade. Supercritical CO2 property tables are constructed from NIST REFPROP data for the condensable gas simulation in FINE/Turbo. A detailed mesh sensitivity study is performed for a baseline design to identify the proper-grid refinement and efficiently allocate resources for the optimization. Seventy design variables are selected for the initial population generation. Self-organizing maps are then used to focus the design variables on the most important ones affecting the objective functions. The optimization results in approximately 3000 three-dimensional Reynolds Averaged Navier Stokes simulations of different blade shapes with increases in efficiency of up to 0.85% and decreases in heat load of 14%. Families of blade shapes are identified for experimental testing in an annular rig at the Purdue Experimental Turbine Aerothermal Laboratory. A design to adapt the annular cascade for testing optimized geometries is introduced, which features eccentric radius sectors allowing for scaled-up geometries of sCO2 optimized blade profiles to be tested at design cycle representative conditions at high Reynolds numbers in dry air. Analysis into the effects of Reynolds number, working fluid, and geometric relations are presented to prove the efficacy of the test method.

Engineering↗

First survey of centimeter-scale AC-LGAD strip sensors with a 120 GeV proton beam

We will present the first beam test results with centimeter-scale AC-LGAD strip sensors, using the Fermilab Test Beam Facility, and a study of the performance of AC-LGAD sensors as a function of their thickness. Sensors of this type are envisioned for applications that require large-area precision 4D tracking coverage with economical channel counts, including timing layers for the Electron Ion Collider (EIC), and space-based particle experiments. Long strip sensors with sparse readout offer better cost and performance for applications where channel count or electrical power density is a constraint. Thanks to the excellent signal to noise ratio in AC-LGADs, sparse readout can be exploited without significant degradation of spatial or time resolution, which is demonstrated in our studies. A survey of sensor designs is presented, with the aim of optimizing the electrode geometry for spatial resolution and timing performance. We will present our studies of the sensor geometry optimization to maintain the desirable sensor performance characteristics with increasingly larger electrodes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Augmenting Molecular Graphs with Geometries via Machine Learning Interatomic Potentials

Accurate molecular property predictions require 3D geometries, which are typically obtained using expensive methods such as density functional theory (DFT). Here, we attempt to obtain molecular geometries by relying solely on machine learning interatomic potential (MLIP) models. To this end, we first curate a large-scale molecular relaxation dataset comprising 3.5 million molecules and 300 million snapshots. Then MLIP pre-trained models are trained with supervised learning to predict energy and forces given 3D molecular structures. Once trained, we show that the pre-trained models can be used in different ways to obtain geometries either explicitly or implicitly. First, it can be used to obtain approximate low-energy 3D geometries via geometry optimization. While these geometries do not consistently reach DFT-level chemical accuracy or convergence, they can still improve downstream performance compared to non-relaxed structures. To mitigate potential biases and enhance downstream predictions, we introduce geometry fine-tuning based on the relaxed 3D geometries. Second, the pre-trained models can be directly fine-tuned for property prediction when ground truth 3D geometries are available. Our results demonstrate that MLIP pre-trained models trained on relaxation data can learn transferable molecular representations to improve downstream molecular property prediction and can provide practically valuable but approximate molecular geometries that benefit property predictions. Our code is publicly available at: https://github.com/divelab/AIRS/.

Fu, Cong [Texas A & M Univ., College Station, TX (↗

The G2(B3LYP/MP2/CC) Approach: A Modification of the G2(MP2) Approach and a Comparison with B3LYP Results

The quadratic configuration interaction calculation in the G2(MP2) approach is replaced by a coupled-cluster singles and doubles calculation including a perturbational estimate of the triples excitations. In addition, the SCF and MP2 geometry optimizations and SCF frequency calculation in the G2(MP2) approach are replaced by a B3LYP geometry optimization and frequency calculation in the proposed G2(B3LYP/MP2/CC) approach. This simplification does not affect the average absolute deviation from experiment, but decreases the maximum error compared with the G2(MP2) approach. The G2(B3LYP/MP2/CC) results are compared with those obtained using the B3LYP approach, and the G2(B3LYP/MP2/CC) model is found to be more reliable, even if the B3LYP calculations are performed using a large basis set.

Bauschlicher, Charles W., Jr.↗

Quadrupolar NMR crystallography guided crystal structure prediction (QNMRX-CSP) of zwitterionic organic HCl salts

In this work, we benchmark quadrupolar NMR crystallography guided crystal structure prediction (QNMRX-CSP) for determining the crystal structures of two zwitterionic organic HCl salts, L-ornithine HCl ( Orn ) and L-histidine HCl·H 2 O ( Hist ). These salts present an interesting challenge for QNMRX-CSP, as gas-phase geometry optimizations used to generate starting structures for the organic zwitterionic fragments fail to capture their correct solid-state geometries. To overcome this limitation, geometry optimizations using the COSMO water-solvation model are employed to generate initial structural models. Using this approach, QNMRX-CSP yields structural models of the two zwitterionic organic HCl salts that closely match experimentally determined crystal structures. In addition, the application of QNMRX-CSP to Hist represents a further step toward the de novo structural determination of solvated organic HCl salts, as Hist is the first benchmark system of this type to include a water molecule as a component of its crystal structure. This work is significant for its potential application to the structural determination of active pharmaceutical ingredients, which often feature complex organic components and solvated solid forms.

Fleischer, Carl H. [Florida State Univ., Tallahass↗

Comparison of Machine Learning Approaches for Prediction of the Equivalent Alkane Carbon Number for Microemulsions Based on Molecular Properties

The chemical properties of oils are vital in the design of microemulsion systems. The hydrophilic–lipophilic difference equation used to predict microemulsions’ phase behavior expresses the oils’ physiochemical properties as the equivalent alkane carbon number (EACN). The experimental determination of EACN requires knowledge of the temperature dependence of the microemulsion system and the effects of different surfactant concentrations. Thus, the experimental determination is time-intensive and tedious, requiring days to months for proper separations. Furthermore, the experiments require high purity of chemicals because microemulsions are sensitive to impurities. Our work focuses on the quick and reliable predictions of the EACN with machine learning (ML) models. Due to the immaturity of ML chemical predictions, we compare three graph neural networks (GNNs) and a gradient-boosted tree algorithm, known as XGBoost. The GNNs use the molecular structures represented as simplified molecular-input line-entry system (SMILES) codes for the initial input, which allows us to assess whether geometry optimization is necessary for reliable results. The XGBoost model also begins with the SMILES representations of the molecules but uses molecular descriptors instead of geometry optimizations. As a result, the best model tested (crystal graph convolutional neural network with Merck molecular force field-94) has an error of 1.15 EACN units of the true EACN for unknown data with the errors skewed toward zero and an R² score of 0.9

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Controlling Ligand Excimer Formation with Dipole Changes in Emissive Rare-Earth/Phosphonic Acid Complexes

The interactions between substituted arylvinyl phosphonic acid (AVPA) ligands within a Eu-AVPA complex are shown to influence the outcomes of excited state evolution after photoexcitation. Compared with unfunctionalized AVPAs, pairs of ligands functionalized with CF 3 in the para position preassociate in the ground state of complexes with Eu 3+ according to calculated geometry optimizations. The CF 3 -substituted AVPA complexes show evidence of red-shifted optical absorption and undergo more efficient excimer formation, as revealed by transient absorption spectroscopy. We rationalize this behavior through simulations of excited-state geometry optimizations that reveal evolution toward interligand phenyl-phenyl planarity for specific excited states. Emission from complexed Eu 3+ after energy transfer from the ligand is found to be weaker with CF 3 substitution, which we hypothesize is due to intracomplex, interligand aggregates with excimer-promoting geometries. These observations point to the need to consider ground-state geometries as well as dynamic excited-state processes to understand the flow of energy in rare earth coordination complexes.

09 BIOMASS FUELS↗

Controlling Molecular Structure and Spin with Multiconfigurational Quantum Chemistry (Final Technical Report)

For many first-row transition metal complexes, structure-property relationships can be obtained from high-level molecular geometry optimizations and subsequent electronic structure studies. The project funded under this award utilized newly implemented fully internally contracted (FIC) nuclear gradients for extended multi-state (XMS) complete active space second-order multireference perturbation theory (CASPT2) to explore geometry changes in first-row transition metal coordination complexes. At the start of the project, only one full geometry optimization using FIC-CASPT2 analytical gradients had been reported (J. Chem. Theory Comput., 2016, 12 (8), 3781), and many open-questions regarding the performance and achievable accuracy in applying such computations to larger complexes persisted. Key deliverables in the report include the implementation of the numerical Hessian and subsequent vibrational analysis in the BAGEL program package. This includes both full Hessian vibrational analysis (FHVA) and a partial Hessian vibrational analysis (PHVA). The latter of which allows us to reduce the significant cost of computing the vibrational frequencies in a large molecule by focusing on the modes of highest interest. The gradient code also proved useful in the development of an approach to improve the Hubbard-U correction by removing bias in predicting spin-splitting in Fe(II) complexes via plane-wave density functional theory (DFT). Finally, we showed for three families of complexes (spin-crossover (SCO) complexes, metallocorroles complexes, and those with metal-metal bonds) that CASPT2 can result in good molecular geometries and established best practices in undertaking this work. We are now using this approach in collaborative projects with experimental groups, which would not have been possible at the start of this project. A common theme also arose through this work that static and dynamic correlation must be recovered in a balanced way to yield quantitative results. By systematically varying how both effects were included, we were able to make key chemical insights. This work supported the training of two postdoctoral scholars, one graduate student, and one undergraduate student in applying high-level wavefunction based methods to challenging systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Topology of modified helical gears and Tooth Contact Analysis (TCA) program

The contents of this report covers: (1) development of optimal geometries for crowned helical gears; (2) a method for their generation; (3) tooth contact analysis (TCA) computer programs for the analysis of meshing and bearing contact of the crowned helical gears; and (4) modelling and simulation of gear shaft deflection. The developed method for synthesis was used to determine the optimal geometry for a crowned helical pinion surface and was directed to localize the bearing contact and guarantee favorable shape and a low level of transmission errors. Two new methods for generation of the crowned helical pinion surface are proposed. One is based on the application of a tool with a surface of revolution that slightly deviates from a regular cone surface. The tool can be used as a grinding wheel or as a shaver. The other is based on a crowning pinion tooth surface with predesigned transmission errors. The pinion tooth surface can be generated by a computer-controlled automatic grinding machine. The TCA program simulates the meshing and bearing contact of the misaligned gears. The transmission errors are also determined. The gear shaft deformation was modelled and investigated. It was found that the deflection of gear shafts has the same effect as gear misalignment.

Litvin, Faydor L.↗

blastforge

BlastForge is a python module meant to aid reinforcement learning research for geometry optimization research projects. The code is built to use PyTorch as a backend and will contain several model architectures and reinforcement learning training loops as well as helper python functions to evaluate a model’s performance during and after training. BlastForge is meant to be a small, focused python project to study moderate-complexity geometry optimization problems

Hickmann, Kyle [Los Alamos National Laboratory]↗

QMMM 2023: A program for combined quantum mechanical and molecular mechanical modeling and simulations

Combined quantum mechanical and molecular mechanical (QM/MM) methods play an important role in multiscale modeling and simulations. QMMM 2023 is a general-purpose program for single-point calculations, geometry optimizations, transition-state optimizations, and molecular dynamics (MD) at the QM/MM level. It calls a QM package and an MM package to perform the required single-level calculations and combines them into a QM/MM energy by a variety of schemes. QMMM 2023 supports GAMESS-US, Gaussian, and ORCA as QM packages and Tinker as the MM package. Four types of treatments are available for embedding the QM subsystem in the MM environment: mechanical embedding with gas-phase calculations of the QM region, electronic embedding that allows polarization of the QM region by the MM environment, polarizable embedding for mutual polarization of the QM and MM regions, and flexible embedding for both mutual polarization and partial charge transfer between the QM and MM regions. Boundaries between QM and MM regions that pass through covalent bonds can be treated by several methods, including the redistributed charge (RC) scheme, redistributed charge and dipole (RCD) scheme, balanced-RC scheme, balanced-RCD scheme, screened charge scheme that takes account of charge penetration effects, and smeared charge scheme that delocalizes the MM charges near the QM–MM boundary. Geometry optimization can be done using the optimizer implemented in QMMM 2023 or the Berny optimizer in Gaussian through external calls to Gaussian. Molecular dynamics simulations can be performed at the pure-MM level, pure-QM level, fixed-partitioning QM/MM level, and adaptive-partitioning QM/MM level. As a result, the adaptive-partitioning treatments permit on-the-fly relocation of the QM–MM boundary by dynamically reclassifying atoms or groups into the QM or MM subsystems.

97 MATHEMATICS AND COMPUTING↗

Review of TEAMER Awards for WEC-Sim Support: Preprint

Testing Expertise and Access for Marine Energy Research (TEAMER) is a U.S. Department of Energy Water Power Technologies Office sponsored program, overseen by the Pacific Ocean Energy Trust, which aims to advance the state of marine energy technologies. The program connects technology developers with experts at U.S. facilities, including numerical modeling and analysis facilities, like WEC-Sim. The WEC-Sim facility is supported by the WEC-Sim development team at Sandia National Laboratories and the National Renewable Energy Laboratory. WEC-Sim (Wave Energy Converter SIMulator) is an open-source software for simulating wave energy converters. WEC-Sim can model the multi-body dynamics of devices comprised of bodies, joints, power take-off systems, and mooring systems. Since TEAMER's first round of support in 2020, there have been eighteen TEAMER awards focused on numerical model development in WEC-Sim. TEAMER awards for WEC-Sim support have modeled a wide range of wave energy converter archetypes, including point absorbers, attenuators, oscillating water columns, and many other novel architectures. A wide variety of studies have been conducted, leading to important insights for TEAMER partners and software improvements for WECSim. This article highlights several successful WEC-Sim TEAMER awards. The awards described herein include TEAMER recipients Ocean Motion Technologies, AquaHarmonics, iProTech, East Carolina University, Virginia Tech, Maiden Wave Energy, and the University of Massachusetts Dartmouth. The awards of these seven partners contain a wide range of investigations and cover the creation of baseline hydrodynamic models, PTO modeling, geometry optimization in both boundary element methods and WECSim, and model tuning and validation.

industry support↗

Phase-Space Geometry and Optimal State Preparation in Quantum Metrology with Collective Spins

We revisit well-known protocols in quantum metrology using collective spins and propose a unifying picture for optimal state preparation based on a semiclassical description in phase space. We show how this framework allows for quantitative predictions of the timescales required to prepare various metrologically useful states, and that these predictions remain accurate even for moderate system sizes, surprisingly far from the classical limit. Furthermore, this framework allows us to build a geometric picture that relates optimal (exponentially fast) entangled probe preparation to the existence of separatrices connecting saddle points in phase space. We illustrate our results with the paradigmatic examples of the two-axis countertwisting and twisting-and-turning Hamiltonians, where we provide analytical expressions for all the relevant optimal timescales. Finally, we propose a generalization of these models to include p-body collective interaction (or p-order twisting), beyond the usual case of p = 2. Using our geometric framework, we prove a no-go theorem for the local optimality of these models for p > 2.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

FIREFLY: heat load and particle exhaust approximations for rapid evaluation of divertor designs

The divertor in a magnetic confinement fusion reactor is an essential component for power dissipation and particle removal. The FIREFLY package for rapid evaluation of divertor designs is presented as an extension of the FLARE code for field line reconstruction from a flux tube mesh. First, divertor loads are approximated with a simplified heat transport model. Neutralized particles are then sampled from the resulting load distribution, and the EIRENE code is used to track molecules and atoms in a plasma background while accounting for dissociation, charge exchange and ionization. Particles are removed on pumping surfaces in order to estimate the exhaust efficiency for a given divertor geometry. Optimization of the divertor geometry for more efficient particle exhaust is explored by using W7-X as an example, and the sensitivity to model parameters for the plasma background in the proxy calculations is evaluated.

mesh generation, magnetic field lines, scrape-off ↗

Optimizing view/illumination geometry for terrestrial features using Space Shuttle and aerial polarimetry

This paper describes to relationship of polarimetric observations from orbital and aerial platforms and the determination optimum sun-target-sensor geometry. Polarimetric observations were evaluated for feature discrimination. The Space Shuttle experiment was performed using two boresighted Hasselblad 70 mm cameras with identical settings with linear polarizing filters aligned orthogonally about the optic axis. The aerial experiment was performed using a single 35 mm Nikon FE2 and rotating the linear polarizing filter 90 deg to acquire both minimum and maximum photographs. Characteristic curves were created by covertype and waveband for both aerial and Space Shuttle imagery. Though significant differences existed between the two datasets, the observed polarimetric signatures were unique and separable.

Israel, Steven A.↗

Optimized Radiator Geometries for Hot Lunar Thermal Environments

The optimum radiator configuration in hot lunar thermal environments is one in which the radiator is parallel to the ground and has no view to the hot lunar surface. However, typical spacecraft configurations have limited real estate available for top-mounted radiators, resulting in a desire to use the spacecraft's vertically oriented sides. Vertically oriented, flat panel radiators will have a large view factor to the lunar surface, and thus will be subjected to significant incident lunar infrared heat. Consequently, radiator fluid temperatures will need to exceed approximately 325 K (assuming standard spacecraft radiator optical properties) in order to provide positive heat rejection at lunar noon. Such temperatures are too high for crewed spacecraft applications in which a heat pump is to be avoided. A recent study of vertically oriented radiator configurations subjected to lunar noon thermal environments led to the discovery of a novel radiator concept that yielded positive heat rejection at lower fluid temperatures. This radiator configuration, called the Intense Thermal Infrared Reflector (ITIR), has exhibited superior performance to all previously analyzed concepts in terms of heat rejection in the lunar noon thermal environment. A key benefit of ITIR is the absence of louvers or other moving parts and its simple geometry (no parabolic shapes). ITIR consists of a specularly reflective shielding surface and a diffuse radiating surface joined to form a horizontally oriented V-shape (shielding surface on top). The point of intersection of these surfaces is defined by two angles, those which define the tilt of each surface with respect to the local horizontal. The optimum set of these angles is determined on a case-by-case basis. The idea assumes minimal conductive heat transfer between shielding and radiating surfaces, and a practical design would likely stack sets of these surfaces on top of one another to reduce radiator thickness.

Ochoa, Dustin↗

Importance of Higher Fidelity Model Geometries during Optimization of Critical Experiments

PARADIGM, PARallel Approach of Differential and InteGral Measurements, is a cross-collaborative effort at Los Alamos National Laboratory between nuclear data theorists, differential and integral experimenters, as well as machine learning statisticians to tackle uncertainties in the intermediate region of 239 Pu. In essence, the idea behind PARADIGM is to remove the linear conceptualization of the nuclear data pipeline, shown in Figure 1, and replace it with a far more parallelized approach. The novel approach leverages machine learning to guide which differential measurements and integral experiments will result in the largest decrease in uncertain ties for a nuclide reaction pair in a given energy range. The concept builds off earlier work, EUCLID, which focused on the fast region of 239 Pu. The practical benefit of having evaluation, differential measurement, and integral experiment personnel in collaboration with machine learning is to represent the entire nuclear data in one snapshot. This enable large reduction in the time to deliver improved nuclear data, which using the PARADIGM approach could be done in 3 years. A general outline of PARADIGM and specific topics are available in other papers. The discussion here will pertain directly to the integral experiment design. More specifically, the process of taking a rough design and transforming it into a finalized neutronic model will be discussed.

97 MATHEMATICS AND COMPUTING↗

ELISA: A Tool for Optimization of Rotor Hover Performance at Low Reynolds Number in the Mars Atmosphere

The Evolutionary aLgorithm for Iterative Studies of Aeromechanics (ELISA) was developed in support of the Rotorcraft Optimization for the Advancement of Mars eXploration (ROAMX) project. ELISA was developed to enable aerodynamic rotor hover optimization for low Reynolds number flows in the Mars atmosphere. The first objective of the algorithm allows for unconventional airfoil parameterization and multi-objective airfoil geometry optimization using OVERFLOW. The Pareto-optimal airfoil sets are converted to a set of Pareto-optimal airfoil decks, providing the lowest drag air foil geometry for each angle of attack, removing the need to arbitrarily select the airfoils to be used in the rotor optimization. The second objective allows for rotor geometry optimization with simultaneous maximization of blade loading and minimization of rotor power using the comprehensive analysis CAMRADII. The result is a Pareto-optimal rotor set, providing the lowest power rotor for each attainable blade loading, and one of the first tools for hover-optimized rotors for high-subsonic low Reynolds number conditions. The airfoil thickness can be modified after the airfoil optimization is complete, allowing for a post-airfoil-optimization adjustment of blade thickness to facilitate conforming to external structural analyses requirements. The relevance of the code is demonstrated with case studies for the ROAMX rotor optimization for Ingenuity-sized single rotors in the Mars atmosphere, a performance study optimizing the chord and twist of Ingenuity’s coaxial rotor resulting in the Sample Recovery Helicopters candidate rotor, and high-subsonic low Reynolds number airfoil optimization providing novel insights for higher-efficiency low Reynolds number airfoil geometries and flow physics.

ELISA↗