SEARCH · Engineering Papers
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.
This poster discusses the modeling, performance and optimization of a solvent absorption system for CO2 capture with towers utilizing intensified packing. The packing is an alternative to traditional structured packing by incorporating cooling channel for simultaneous mass and heat transfer.
The user's manual for a computer program that performs system optimization of gasdynamic lasers is provided. Detailed input/output formats are CDC 7600/6600 computers using a dialect of FORTRAN. Sample input/output data are provided to verify correct program operation along with a program listing.
The thermal-control requirements for design-optimized aeromaneuvering performance are determined for space-based applications and low-earth orbit sorties involving large, multiple plane-inclination changes. The leading-edge heating analysis is the most advanced developed for hypersonic-rarefied flow over lifting surfaces at incidence. The effects of leading-edge bluntness, low-density viscous phenomena, and finite-rate flow-field chemistry and surface catalysis are accounted for. The predicted aerothermodynamic heating characteristics are correlated with thermal-control and flight-performance capabilities. The mission payload capability for delivery, retrieval, and combined operations is determined for round-trip sorties extending to polar orbits. Recommendations are given for future design refinements. The results help to identify technology issues required to develop prototype operational systems.
Solving the Alternating Current Optimal Power Flow (ACOPF) problem is key to economically efficient and reliable power networks with a good solution potentially saving utilities tens of billions of dollars annually (according to FERC). The Grid Optimization (GO) Competition set up by ARPA-E saw several promising solutions in Challenge 1. While team GOT-TJU-OPF placed top 10 in Division 3 and 4, this was not a satisfactory performance, and the team has identified specific areas to improve and will be adding more members to round out the necessary skills and expertise needed to be more competitive. The team set out for redemption during GO Challenge 2 with an improved High-Performance Solution for Security Constrained Optimal Power Flow. The (renamed) BSI-GOT-OPF Team ended up finishing top 2 overall in the competition. The algorithms developed have potential impacts for the electric power markets that are enormous. Optimal Power Flow technology can be an enabling technology to achieve energy-efficient power grids while enhancing renewable energy penetration, among others.
The quest for reduced levelized cost of energy has driven significant growth in wind turbine size; however, larger rotors face significant technical and logistical challenges. The largest published rotor design is 25 MW, and here we consider an even larger 50 MW design with blade length over 250 m. This paper shows that a 50 MW design is indeed possible from a detailed engineering perspective and presents a series of aero-structural blade designs, and critical assessment of technology pathways and challenges for extreme-scale rotors. The 50 MW rotor design begins with Monte Carlo simulations focused on optimizing carbon spar cap and root design. A baseline design resulted in a 250-m blade with mass of 502 tonnes. Subsequently, an aero-structural design and optimization were performed to reduce the blade mass/cost with more than 25% mass reduction and 30% cost reduction by determining optimal blade chord and airfoil thickness for best aero-structural performance.
Abstract While some engineering fields have benefited from systematic design optimization studies, wave energy converters have yet to successfully incorporate such analyses into practical engineering workflows. The current iterative approach to wave energy converter design leads to sub-optimal solutions. This short paper presents an open-source MATLAB toolbox for performing design optimization studies on wave energy converters where power take-off behavior and realistic constraints can be easily included. This tool incorporates an adaptable control co-design approach, in that a constrained optimal controller is used to simulate device dynamics and populate an arbitrary objective function of the user’s choosing. A brief explanation of the tool’s structure and underlying theory is presented. To demonstrate the capabilities of the tool, verify its functionality, and begin to explore some basic wave energy converter design relationships, three conceptual case studies are presented. In particular, the importance of considering (and constraining) the magnitudes of device motion and forces in design optimization is shown.
Here, using a thick NMC811 (LiNi 0.8 Mn 0.1 Co 0.1 O 2 ) electrode as an example, we present a macro- to nanoscale 2D and 3D imaging analysis approach coupled with 4D (space + time) computational modeling to probe its degradation mechanism in a lithium-ion battery cell. Particle cracking increases and contact loss between particles and carbon-binder domain are observed to correlate with the cell degradation. This study unravels that the reaction heterogeneity within the thick cathode caused by the unbalanced electron conduction is the main cause of the battery degradation over cycling. The increased heterogeneity in the system will entail more cathode regions where the degree of active material utilization is uneven, leading to higher probabilities of particle cracking. These findings shed light on the crucial role of the electronic and ionic transportation networks in the performance deterioration of the thick cathode. They also provide guidance for cathode architecture optimization and performance improvement.
Layered nickel-rich materials are promising next-generation cathode materials for lithium ion batteries due to their high capacity and low cost. However, the poor thermal stability and longtime cycling performance hinders the commercial applications of high nickel materials. Doping with heteroatoms has been an effective approach for improving electrochemical performance of cathode materials. Controlling doping concentration and geometrical distribution is desired for optimal electrochemical performance, but it is challenging in traditional co-precipitation methods. In this work, controlled dysprosium (Dy) doping to NCM811 was studied in an aerosol synthesis method by controlling the precursor concentrations and heating parameters. The obtained materials were characterized by SEM, XRD, and XPS, and their electrochemical properties and thermal stability were evaluated. By controlling the doping concentration (1.5%), Dy-doped NCM811 was improved simultaneously in long-term cycling and high-rate performance. Here, the thermal-chemical stability of the Dy-doped cathode materials was examined in a microflow reactor with a mass spectrometer. The results showed that Dy-doping shifted the O 2 onset temperature to a higher temperature and reduced O 2 release by 80%, thus dramatically increasing the thermal-chemical stability and improving the fire safety of cathode materials. Since high temperature aerosol synthesis is a low-cost and scalable method, the findings in this work have broad implications for commercial synthesis of novel materials with controlled doping modification to achieve high electrochemical performance and safety in lithium ion batteries.
Here, dielectric microstructures have generated much interest in recent years as a means of accelerating charged particles when powered by solid state lasers. The acceleration gradient (or particle energy gain per unit length) is an important figure of merit. To design structures with high acceleration gradients, we explore the adjoint variable method, a highly efficient technique used to compute the sensitivity of an objective with respect to a large number of parameters. With this formalism, the sensitivity of the acceleration gradient of a dielectric structure with respect to its entire spatial permittivity distribution is calculated by the use of only two full-field electromagnetic simulations, the original and ‘adjoint’. The adjoint simulation corresponds physically to the reciprocal situation of a point charge moving through the accelerator gap and radiating. Using this formalism, we perform numerical optimizations aimed at maximizing acceleration gradients, which generate fabricable structures of greatly improved performance in comparison to previously examined geometries.
This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.
Numerical calculations of enhancement factors offered by dynamic nuclear polarization in solids under magic angle spinning (DNP-MAS) were performed to determine the optimal EPR parameters for a dinitroxide polarizing agent. We found that the DNP performance of a biradical is more tolerant to the relative orientation of the two nitroxide moieties than previously thought. In general, any condition in which the gyy tensor components of both radicals are perpendicular to one another is expected to have near-optimal DNP performance. These results highlight the important role of the exchange coupling, which can lessen the sensitivity of DNP performance to the inter-radical distance, but also lead to lower enhancements when the number of atoms in the linker becomes less than three. Finally, the calculations showed that the electron T1e value should be near 500μs to yield optimal performance. Importantly, the newest polarizing agents already feature all of the qualities of the optimal polarizing agent, leaving little room for further improvement. Further research into DNP polarizing agents should then target non-nitroxide radicals, as well as improvements in sample formulations to advance high-temperature DNP and limit quenching and reactivity.
This work aims to combine an aerodynamic model based on the unsteady vortex lattice method with an acoustic model provided by Farassat's formulation 1A to perform gradient-based optimizations of a proprotor with aerodynamic and acoustic constraints. The resulting combination of tools is applied to the problem of designing a single proprotor operating at a cruise condition, with and without an acoustic constraint and wing placed downstream of the proprotor rotation plane. Results are compared to a baseline design studied previously, and to similar optimizations performed with a simpler blade element momentum theory aerodynamic model. Each optimization case achieved feasibility and made significant improvements in the objective function, but the optimality criterion was not satisfied. Overall the designs the optimizer found roughly comported with our previous experience with similar problems, with some discrepancies that are discussed.
This work aims to combine an aerodynamic model based on the unsteady vortex lattice method with an acoustic model provided by Farassat's formulation 1A to perform gradient-based optimizations of a proprotor with aerodynamic and acoustic constraints. The resulting combination of tools is applied to the problem of designing a single proprotor operating at a cruise condition, with and without an acoustic constraint and wing placed downstream of the proprotor rotation plane. Results are compared to a baseline design studied previously, and to similar optimizations performed with a simpler blade element momentum theory aerodynamic model. Each optimization case achieved feasibility and made significant improvements in the objective function, but the optimality criterion was not satisfied. Overall the designs the optimizer found roughly comported with our previous experience with similar problems, with some discrepancies that are discussed.
In this work we perform quantum optimal control simulations, based on the Time-Dependent Hartree (TDH) approximation, for systems of three to five dipole-dipole coupled OCS rotors. A control electric field is used to steer all of the individual rotors, arranged in chains and regular polygons in a plane, toward either identical or unique objectives. The goal is to explore the utility of the TDH approximation to model the field-induced dynamics of multiple interacting rotors in the weak dipole-dipole coupling regime. A stochastic hill climbing approach is employed to seek an optimal control field that achieves the desired objectives at a specified target time. We first show that multiple rotors in chain and polygon geometries can be identically oriented in the same direction; these cases do not significantly depend on the presence of the dipole-dipole interaction. Additionally, in particular geometrical arrangements, we demonstrate that individual rotors can be uniquely manipulated toward different objectives with the same field. Specifically, it is shown that for a three rotor chain, the two end rotors can be identically oriented in a specific direction while keeping the middle rotor in its ground state, and for an equilateral triangle, two rotors can be identically oriented in a specific direction while the third rotor is oriented in the opposite direction. These multirotor unique objective cases exploit the shape of the field in coordination with dipole-dipole coupling between the rotors. Comparisons to numerically exact calculations, utilizing the TDH-determined fields, are given for all optimal control studies involving systems of three rotors.
Rational design of interfaces with proper physico-chemical properties is necessary to optimize the performance of electrochemical devices, which requires fundamental understanding of the interfaces. In this project, we leveraged quantum mechanical and molecular dynamics simulations as well as machine learning (ML) technique to correlate the relationships between structure and properties including ion transport and electronic structures of cathode host materials in lithium-sulfur (Li-S) and carbon anode in sodium (Na) batteries for improved electrochemical performance.
The cloud microphysics scheme, CASIM, and the radiation scheme, SOCRATES, are two computationally intensive parts within the Met Office's Unified Model (UM). This study enables CASIM and SOCRATES to use accelerated multi-core systems for optimal computational performance of the UM. Using profiling to guide our efforts, we refactored the code for optimal threading and kernel arrangement and implemented OpenACC directives manually or through the CLAW source-to-source translator. Initial porting results achieved 10.02x and 9.25x speedup in CASIM and SOCRATES respectively on 1 GPU compared with 1 CPU core. A granular performance analysis of the strategy and bottlenecks are discussed. These improvements will enable UM to run on heterogeneous computers and a path forward for further improvements is provided.
Lignin-based carbon fiber represents a renewable and low-cost alternative to petroleum-based carbon fiber. However, the poor mechanical performance of current lignin carbon fiber hinders its application. We hypothesized that the lower optimal mechanical performance is caused by the inherent heterogeneity of lignin. It is still unknown how the molecular weights (MWs) and lignin uniformity will impact the performance of lignin carbon fiber. We thereby addressed this hypothesis by fractionating lignin into fractions with different MWs and polydispersity indices (PDIs). An enzyme–mediator-based method and a dialysis method were developed to derive lignin fractions with increased MW and decreased PDI. Lignin fractions were electro-spun into fibers after blending with polyacrylonitrile (PAN) at 1 : 1 (w/w) ratio. The fractionation in general improved the spinnability of lignin to allow us to obtain finer lignin-based carbon fibers. The elastic modulus of lignin carbon fibers, as measured by nanoindentation, was increased as the lignin MW increased and as PDI decreased. The scatter plot and linear regression revealed very good correlation between the elastic modulus and PDI, as well as certain correlation between the elastic modulus and MW. XRD and Raman spectroscopy revealed that the crystallite size and the content of the pre-graphitic turbostratic carbon increased with higher lignin MW and lower PDI, revealing the mechanism of the improvement in carbon fiber mechanical performance. Finally, this study elucidated the impacts of lignin MW and uniformity on the mechanical properties of carbon fiber, and could thus guide the development of lignin processing technologies for quality lignin-based carbon fiber.
Recent advances in nanotechnology have enabled researchers to manipulate small collections of quantum-mechanical objects with unprecedented accuracy. In semiconductor quantum-dot qubits, this manipulation requires controlling the dot orbital energies, the tunnel couplings, and the electron occupations. These properties all depend on the voltages placed on the metallic electrodes that define the device, the positions of which are fixed once the device is fabricated. While there has been much success with small numbers of dots, as the number of dots grows, it will be increasingly useful to control these systems with as few electrode voltage changes as possible. Here, we introduce a protocol, which we call the "compressed optimization of device architectures" (CODA), in order both to efficiently identify sparse sets of voltage changes that control quantum systems and to introduce a metric that can be used to compare device designs. As an example of the former, we apply this method to simulated devices with up to 100 quantum dots and show that CODA automatically tunes devices more efficiently than other common nonlinear optimizers. To demonstrate the latter, we determine the optimal lateral scale for a triple quantum dot, yielding a simulated device that can be tuned with small voltage changes on a limited number of electrodes.