SEARCH · Search NASA

Results for “convex 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

628 records · Page 35

All sources NASA NTRS DOE OSTI

DOE OSTI2018

A centralized-based method to determine the local voltage control strategies of distributed generator operation in active distribution networks

The increasing penetration of distributed generators (DGs) exacerbates the risk of voltage violations in active distribution networks (ADNs). The var capacity provided by DG inverters is a potential solution for voltage regulation. The conventional centralized control strategies limited by computation and communication burdens are difficult to meet the requirement of rapid voltage control. However, the local control strategies based on real-time measurements have a fast response to the frequent fluctuations of DG outputs. The performance of these local controllers depends on the tuning of the control parameters. This paper proposes a centralized-based method to determine the local voltage control strategies for DGs. A centralized parameter tuning model of control curves is built, in which and curves of DG inverters are mathematically formulated based on piecewise linearization. The original mixed-integer nonlinear programming (MINLP) model is converted into an effectively solved mixed-integer second-order cone programming (MISOCP) model using convex relaxation. The potential benefits of DG inverters are explored to regulate both reactive and curtailed active power, based on local voltage measurements. Case studies on a modified PG&E 69-node distribution system are carried out to verify the effectiveness of the proposed method. Results show that power losses of ADNs are significantly reduced and voltage profiles are also improved.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2016

Polynuclear Aromatic Hydrocarbons with Curved Surfaces: Buckyballs

The discovery of a new allotropic form of elemental carbon – the fullerenes – and subsequently other novel forms of elemental carbon with pyramidalized surfaces, most notably single-walled and multi-walled carbon nanotubes, introduced a novel structural motif to the polycyclic aromatic hydrocarbons (PAHs) with nonplanar surfaces. Our research program supported by BES DOE grant DE-FG02-04ER15514 has dealt with the synthesis, structural studies, and chemistry of the novel curved-surface PAHs with carbon frameworks structurally related to fullerenes. They are referred to as “buckybowls”. We prepared several new buckybowls and, even more importantly, developed the efficient, gram-scale synthetic methodologies for the preparation of small buckybowls, most notably corannulene (C20H10) and its derivatives. In addition, the employment of the corannulene-based synthons previously developed in our laboratory led to a number of highly nonplanar molecular architectures with two or more corannulene subunits with a potential for the applications as novel materials in separation sciences, nanoelectronics, photovoltaics and catalysis. In collaboration with Professor Angelici (Iowa State) we prepared and characterized several transition metal complexes of corannulene, providing the first structural characterization of η6 metal complexes of buckybowls by a single crystal X-ray diffraction. In addition to the definitive structural characterization of the complexes we demonstrated that the (η6-C6Me6)Ru2+ unit in some relatively stable complexes activate the corannulene ligand to react with proper nucleophiles suggesting that such complexex may be used in catalysis. (Section C). We have explored the efficiency of the dispersion-based interactions of curved-surface conjugated carbon networks by high-level computational models. We showed that the curvature of such networks does not reduce the van der Waals attractions as compared to the planar systems of similar size. We than concentrated on the design, synthesis and testing of the previously unknown molecular receptors with corannulene pincers which are capable of forming the strong ball-and-socket inclusion complexes with fullerenes in solution and in the solid state. Buckycatcher I (a molecular receptor with two corannulene pincers preorganized on a tetrabensocyclooctatetraene tether) inclusion complexes with fullerenes provided the first experimental evidence for the importance of concave-convex π – π interactions in the supramolecular chemistry of fullerene carbon cages with buckybowls. C60@Buckycatcher I has become the prototypical supramolecular system formed by the relatively weak (on the atom-to-atom bases) and has been extensively used by the computational community to test the quality of the theoretical models. In addition, Buckycatcher I shows an exceptional ability to adopt other guest molecules and several its inclusion complexes and/or solvates have been structurally characterized. These unusual structures show the potential of the buckycatchers to be tested as novel organic materials. “Intelligent design” of molecular receptors with corannulene receptors using computational approach allowing for a priori prediction of their binding potentials led to the preparation of other buckycatchers. Both bi- and tridentate corannulene pincers were tested and two of the bidentate receptors (Buckycatchers II and III) exhibited the outstanding affinity toward fullerenes, significantly exciding Buckycatcher I performance. Careful tailoring of the tethers resulted in the optimization of binding energies of the receptors with guest carbon cages and in the reduction of the entropy/solvation penalties. Practical preparation of dicorannulenopentacene opens a new avenue for the synthesis of a pool of bis-corannulene receptors possessing polar groups on their dibenzobarrelene tether. The polar “anchors” can be used to attach these efficient fullerene receptors to solid supports in order to modify their surfaces with the potential applications in separation sciences, catalysis and photovoltaic materials. additionally, the presence of such groups will improve the solubility of the receptors in polar solvents, expanding the scope of their supramolecular chemistry in solutions. Rigorous studies of thermodynamics of host-guest association of molecular clips and tweezers in organic solvents with C60/C70 inclusion complexes with Buckycatcher I as a model system have provided a full set of thermodynamical data including ΔG, ΔH, and –TΔS, for the formation of the inclusion complexes over a range of temperatures and in several solvents. Both 1H NMR and Isothermal Calorimetry (ITC) titrations provide virtually identical association constants, confirming the reliability of the results. WE have shown that thermodynamics of the host-guest complex formation is more complicated than anticipated and proposed by the existing computational models. These data are of premium importance for the testing of the computational solvation models and for a deeper understanding of the supramolecular interactions in organic solvents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

On the Solution of ℓ 0 -Constrained Sparse Inverse Covariance Estimation Problems

The sparse inverse covariance matrix is used to model conditional dependencies between variables in a graphical model to fit a multivariate Gaussian distribution. Estimating the matrix from data are well known to be computationally expensive for large-scale problems. Sparsity is employed to handle noise in the data and to promote interpretability of a learning model. Although the use of a convex ℓ 1 regularizer to encourage sparsity is common practice, the combinatorial ℓ 0 penalty often has more favorable statistical properties. In this paper, we directly constrain sparsity by specifying a maximally allowable number of nonzeros, in other words, by imposing an ℓ 0 constraint. Here, we introduce an efficient approximate Newton algorithm using warm starts for solving the nonconvex ℓ 0 -constrained inverse covariance learning problem. Numerical experiments on standard data sets show that the performance of the proposed algorithm is competitive with state-of-the-art methods.

$\ell_0$-Constrained↗

DOE OSTI2026

Regional surrogates for predictive control of digital twins

Digital twins of complex systems must involve a model that is fast, generalizable, and usable for real-time control. For example, high-fidelity nonlinear multiphysics simulations can capture laser-material interactions, but are too slow for optimization or model predictive control (MPC). Reduced-order models, used to accelerate such computation, frequently fail to generalize to unseen inputs or control states. We show theoretically that this failure is intrinsic, i.e., that a learned model is non-unique outside the sampled subspace when its low-rank structure arises from limited excitation and clustered eigenvalues, rather than from a user-imposed truncation alone. Motivated by this result, we propose a control-ready regional surrogate-construction framework for both autonomous and nonautonomous dynamics; it employs Koopman lifting to represent nonlinearities, while preserving spatial locality. We illustrate our approach by constructing a control-ready surrogate for the digital twin of a thermal component of additive-manufacturing process. Our surrogate, localized in space through a von Neumann stencil, is learned from noisy high-fidelity simulations that emulate thermal-camera images collected during the manufacturing. It is linear in thermo-physically augmented states so that MPC reduces to a convex quadratic program. The surrogate requires no online correction, generalizes to unseen scan paths and power profiles of the laser, and is more than three orders of magnitude faster than a finite-difference solver. Furthermore, when the MPC sequence computed on the digital twin is applied to this solver, closed-loop temperature regulation is recovered, showing that the surrogate preserves control-relevant input-output behavior.

Data-driven model↗

DOE OSTI2017

High-Throughput Combinatorial Development of High-Entropy Alloys For Light-Weight Structural Applications

The primary limitation of today’s lightweight structural alloys is that specific yield strengths (SYS) higher than 200MPa x cc/g (typical value for titanium alloys) are extremely difficult to achieve. This holds true especially at a cost lower than 5dollars/kg (typical value for magnesium alloys). Recently, high-entropy alloys (HEA) have shown promising SYS, yet the large composition space of HEA makes screening compositions complex and time-consuming. Over the course of this 2-year project we started from 150 billion compositions and reduced the number of potential low-density (<5g/cc), low-cost (<5dollars/kg) high-entropy alloy (LDHEA) candidates that are single-phase, disordered, solid-solution (SPSS) to a few thousand compositions. This was accomplished by means of machine learning to guide design for SPSS LDHEA based on a combination of recursive partitioning, an extensive, experimental HEA database compiled from 24 literature sources, and 91 calculated parameters serving as phenomenological selection rules. Machine learning shows an accuracy of 82% in identifying which compositions of a separate, smaller, experimental HEA database are SPSS HEA. Calculation of Phase Diagrams (CALPHAD) shows an accuracy of 71-77% for the alloys supported by the CALPHAD database, where 30% of the compiled HEA database is not supported by CALPHAD. In addition to machine learning, and CALPHAD, a third tool was developed to aid design of SPSS LDHEA. Phase diagrams were calculated by constructing the Gibbs-free energy convex hull based on easily accessible enthalpy and entropy terms. Surprisingly, accuracy was 78%. Pursuing these LDHEA candidates by high-throughput experimental methods resulted in SPSS LDHEA composed of transition metals (e.g. Cr, Mn, Fe, Ni, Cu) alloyed with Al, yet the high concentration of Al, necessary to bring the mass density below 5.0g/cc, makes these materials hard and brittle, body-centered-cubic (BCC) alloys. A related, yet multi-phase BCC alloy, based on Al-Cr-Fe-Ni, shows compressive strain >10% and specific compressive yield strength of 229 MPa x cc/g, yet does not show ductility in tensile tests due to cleavage. When replacing Cr in Al-Cr-Fe-based 4- and 5-element LDHEA with Mn, hardness drops 2x. Combined with compression test results, including those on the ternaries Al-Cr-Fe and Al-Mn-Fe suggest that Al-Mn-Fe-based LDHEA are still worth pursuing. These initial results only represent one compressive stress-strain curve per composition without any property optimization. As such, reproducibility needs to be followed by optimization to show their full potential. When including Li, Mg, and Zn, single-phase Li-Mg-Al-Ti-Zn LDHEA has been found with a specific ultimate compressive strength of 289MPa x cc/g. Al-Ti-Mn-Zn showed a specific ultimate compressive strength of 73MPa x cc/g. These initial results after hot isostatic pressing (HIP) of the ball-milled powders represent the lower end of what is possible, since no secondary processing (e.g. extrusion) has been performed to optimize strength and ductility. Compositions for multi-phase (e.g. dual-phase) LDHEA were identified largely by automated searches through CALPHAD databases, while screening for large face-centered-cubic (FCC) volume fractions, followed by experimental verification. This resulted in several new alloys. Li-Mg-Al-Mn-Fe and Mg-Mn-Fe-Co ball-milled powders upon HIP show specific ultimate compressive strengths of 198MPa x cc/g and 45MPa x cc/g, respectively. Several malleable quarternary Al-Zn-based alloys have been found upon arc/induction melting, yet with limited specific compressive yield strength (<75 MPa x cc/g). These initial results are all without any optimization for strength and/or ductility. High-throughput experimentation allowed us to triple the existing experimental HEA database as published in the past 10 years in less than 2 years which happened at a rate 10x higher than previous methods. Furthermore, we showed that high-throughput thin-film combinatorial methods can be used to get insight in isothermal phase diagram slices. Although it is straightforward to map hardness as a function of composition for sputtered, thin-film, compositional gradients by nano-indentation and compare the results to micro-indentation on bulk samples, the simultaneous impact of composition, roughness, film density, and microstructure on hardness requires monitoring all these properties as a function of location on the compositional gradient, including dissecting the impact of these 4 factors on the hardness map. These additional efforts impact throughput significantly. This work shows that a lot of progress has been made over the years in predicting phase formation that aids the discovery of new alloys, yet that a lot of work needs to be done to predict phases more accurately for LDHEA, whether done by CALPHAD or by other means. More importantly, more work needs to be done to predict mechanical properties of novel alloys, like yield strength, and ductility. Furthermore, this work shows that there is a need for the generation of an empirical alloy database covering strategic points in a multi-dimensional composition space to allow for faster and more accurate predictive interpolations to identify the oasis in the dessert more quickly. Finally, this work suggests that it is worth pursuing a ductile alloy with a SYS > 300 MPa x cc/g in a mass density range of 6-7 g/cc, since the chances for a single-phase or majority-phase FCC increase significantly. Today’s lightweight steels are in this density range.

36 MATERIALS SCIENCE↗

NASA NTRSDate not supplied

Volumetric Assessment of UPRITE Exercises From Marker-Based Motion Capture

BACKGROUND Lack of volumetric data on full-body movement of exercises presents a challenge to ensuring the fit of crew member’s full range of motion on the International Space Station (ISS). The Upright Proprioception Retention via In-flight Training and Evaluation (UPRITE) is a sensorimotor countermeasure device designed for maintaining crew members’ proprioception in a microgravity environment. A footplate—attached to a static base—rotates in two degrees of freedom (pitch and roll) up to a 20 deg angle. An initial volumetric assessment assuming an upright standing posture produced a cone-like shape with a narrow bottom and wide top. Such general volumetric assessments risk creating an overly conservative volume estimate, taking up more space than is necessary on the already limited interior space of the ISS, and neglecting necessary volume due to oversimplifying assumptions. Rather, higher-fidelity volumetric assessments offer more comprehensive insights in an environment where every area counts. The main objective of this work is to provide the spatial parameters of exercises on the UPRITE such that it is placed on the ISS according to its volumetric demands or that usage is adjusted to fit the available space. METHODS In 2023, a data collection was performed originally to inform loads and dynamics of system use and was recently leveraged for volumetric assessment. Three human subjects representing different body types (~63-76 inches in stature) performed a variety of board manipulations using UPRITE with body weight offload. The test collected the 3D positional data of a modified full-body Plug-in Gait marker set [1] via a 16-camera OptiTrack MoCap system. After processing – filling marker gaps and trimming data – in OptiTrack Motive, the recorded marker location data, which included device markers, was exported to a readable trajectory file. To accurately represent the full volume defining landmarks, additional markers were digitally added to an unscaled Modified Full Body Model [2]. The model was then scaled according to its subject parameters upon which an inverse kinematics analysis was performed. A custom plugin yielded model marker location data files. Volumetric analyses were performed on the recorded trajectory and model trajectory files using a custom Python-built tool that extracted the marker location data and plotted it in a 3D space. Concerned with only the maximum volume of the motion, a 3D convex hull analysis was applied to the plot, extracting the vertices or external points of the eventual 3D CAD output, dubbed aptly as a “volume shell”. This overall approach was based on guidance in a NASA-STD-3001 Technical Brief [3]. RESULTS AND DISCUSSION Batch volumetric assessment on the exercises for each subject was performed, producing high-fidelity volume shells in minimal time. Preliminary results highlighted the value in higher-fidelity volumes based on collected data when possible. For example, revolving a single posture in the cone assessment would not have sufficiently captured a single leg stance; rather, it would need to involve swinging the leg both forward and back. Additional observations and the maximal dimensions of the volumes, including those based on scaled data for ISS anthropometric requirements, will be presented at the Human Research Program Investigator’s Workshop. CONCLUSIONS While this work’s primary objective was for the UPRITE-to-ISS integration, the tool built to conduct this analysis has wide applications for future exercise systems as an informational tool for optimal device placement. The tool and its findings also have implications for exercise device design and spacecraft interior considerations on Gateway, the Lunar Pressurized Rover, and beyond. REFERENCES [1] Bell, C. A., et al. (2023) Recent Improvements and Verification of a Full Body Model in OpenSim. NASA Human Research Program Investigator’s Workshop. https://ntrs.nasa.gov/citations/20230001080 [2] Lostroscio, K., et al (2023) The Digital Astronaut Simulation. AHFE International Conference on Human Factors in Design, Engineering, and Computing for All. [3] Exercise Overview. (2023) NASA-STD-3001 Technical Brief. https://www.nasa.gov/wp-content/uploads/2023/12/ochmo-tb-031-exercise-overview.pdf?emrc=9d454c?emrc=9d454c

L D Quinto↗

DOE OSTI2020

Proppant embedment in coal and shale: Impacts of stress hardening and sorption

During methane production in CBM reservoirs, the influence of proppant embedment and permeability damage cannot be neglected – especially where the wall-rock is soft. Effective stresses are elevated during methane recovery, increasing both normal loading stress and confinement and simultaneously overprinting sorption-induced volumetric strains. Experiments and analytic modeling are conducted to define key mechanisms controlling these competitive effects. We independently measure overall sample compaction (external LVDT) and local strain (strain gauge) in the matrix to deconvolve proppant embedment in a propped fracture for different conditions of confining stress. The results show symptomatic behaviors of elastic (shale) and elastoplastic (coal) responses of embedment. Different from shale, the evolution of embedment is convex upwards with increased stress where indented depth increases more rapidly as loading stress increases under constant confinement. In addition, a stress-hardening effect is found to play a pivotal role in determining the characteristics of indentation, which are examined in terms of evolution profiles, deformation regimes, embedment slopes, curvatures, yield points and irreversible indentations. Based on the experimental observations a semianalytical model predicts indentation and the evolution of propped permeability under recreated in-situ stress conditions. A simplified case study is conducted to further illustrate the evolution of aperture and permeability of a propped fracture in CBM reservoirs. The modeling results suggest that proppant embedment is significantly overestimated if the variable stress-hardening (VSH) effect is neglected, especially when effective stress is large. Moreover, a decrease in indentation depth possibly occurs during late stage methane production, resulting in a reversal/recovery in fracture closure. This is because desorption-induced shrinkage becomes the predominant effect, causing an increase in aperture and a reduction in the indented volume of proppant. The resulting recovery in permeability implies that the propped coal fracture has the potential to optimally facilitate methane production as a pathway, even at high closure stresses generated by methane drainage.

01 COAL, LIGNITE, AND PEAT↗

NASA NTRS2016

Average Cross-Sectional Area of DebriSat Fragments Using Volumetrically Constructed 3D Representations

Debris fragments from the hypervelocity impact testing of DebriSat are being collected and characterized for use in updating existing satellite breakup models. One of the key parameters utilized in these models is the ballistic coefficient of the fragment which is directly related to its area‐to‐mass ratio. However, since the attitude of fragments varies during their orbital lifetime, it is customary to use the average cross‐sectional area in the calculation of the area‐to‐mass ratio. The average cross‐sectional area is defined as the average of the projected surface areas perpendicular to the direction of motion and has been shown to be equal to one‐fourth of the total surface area of a convex object. Unfortunately, numerous fragments obtained from the DebriSat experiment show significant concavity (i.e., shadowing) and thus we have explored alternate methods for computing the average cross‐sectional area of the fragments. An imaging system based on the volumetric reconstruction of a 3D object from multiple 2D photographs of the object was developed for use in determining the size characteristic (i.e., characteristics length) of the DebriSat fragments. For each fragment, the imaging system generates N number of images from varied azimuth and elevation angles and processes them using a space‐carving algorithm to construct a 3D point cloud of the fragment. This paper describes two approaches for calculating the average cross‐sectional area of debris fragments based on the 3D imager. Approach A utilizes the constructed 3D object to generate equally distributed cross‐sectional area projections and then averages them to determine the average cross‐sectional area. Approach B utilizes a weighted average of the area of the 2D photographs to directly compute the average cross‐sectional area. A comparison of the accuracy and computational needs of each approach is described as well as preliminary results of an analysis to determine the "optimal" number of images needed for the 3D imager to accurately measure the average cross sectional area of objects with known dimensions.

Scruggs, T.↗

NASA NTRS1993

Techniques for estimating Space Station aerodynamic characteristics

A method was devised and calculations were performed to determine the effects of reflected molecules on the aerodynamic force and moment coefficients for a body in free molecule flow. A procedure was developed for determining the velocity and temperature distributions of molecules reflected from a surface of arbitrary momentum and energy accommodation. A system of equations, based on momentum and energy balances for the surface, incident, and reflected molecules, was solved by a numerical optimization technique. The minimization of a 'cost' function, developed from the set of equations, resulted in the determination of the defining properties of the flow reflected from the arbitrary surface. The properties used to define both the incident and reflected flows were: average temperature of the molecules in the flow, angle of the flow with respect to a vector normal to the surface, and the molecular speed ratio. The properties of the reflected flow were used to calculate the contribution of multiply reflected molecules to the force and moments on a test body in the flow. The test configuration consisted of two flat plates joined along one edge at a right angle to each other. When force and moment coefficients of this 90 deg concave wedge were compared to results that did not include multiple reflections, it was found that multiple reflections could nearly double lift and drag coefficients, with nearly a 50 percent increase in pitching moment for cases with specular or nearly specular accommodation. The cases of diffuse or nearly diffuse accommodation often had minor reductions in axial and normal forces when multiple reflections were included. There were several cases of intermediate accommodation where the addition of multiple reflection effects more than tripled the lift coefficient over the convex technique.

Thomas, Richard E.↗

DOE OSTI2017

On extreme points of the diffusion polytope

Here, we consider a class of diffusion problems defined on simple graphs in which the populations at any two vertices may be averaged if they are connected by an edge. The diffusion polytope is the convex hull of the set of population vectors attainable using finite sequences of these operations. A number of physical problems have linear programming solutions taking the diffusion polytope as the feasible region, e.g. the free energy that can be removed from plasma using waves, so there is a need to describe and enumerate its extreme points. We also review known results for the case of the complete graph Kn, and study a variety of problems for the path graph Pn and the cyclic graph Cn. Finall, we describe the different kinds of extreme points that arise, and identify the diffusion polytope in a number of simple cases. In the case of increasing initial populations on Pn the diffusion polytope is topologically an n-dimensional hypercube.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Computational optimal transport for molecular spectra: The fully continuous case

Computational optimal transport is used to analyze the difference between pairs of continuous molecular spectra. It is demonstrated that transport distances which are derived from this approach may be a more appropriate measure of the difference between two continuous spectra than more familiar measures of distance under many common circumstances. Associated with the transport distances is the transport map which provides a detailed analysis of the difference between two molecular spectra and is a key component of our study of quantitative differences between two continuous spectra. The use of optimal transport for comparing molecular spectra is developed in detail here with a set of model spectra, so that the discussion is self-contained. The difference between the transport distance and more common definitions of distance is elucidated for some well-chosen examples and it is shown where transport distances may be very useful alternatives to standard definitions of distance. The transport distance between a theoretical and experimental electronic absorption spectrum for SO 2 is studied and it is shown how the theoretical spectrum can be modified to fit the experimental spectrum better adjusting the theoretical band origin and the resolution of the theoretical spectrum. In conclusion, this analysis includes the calculation of transport maps between the theoretical and experimental spectra suggesting future applications of the methodology.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

Novel artificial neural network model for instantaneous power losses and operational efficiency mapping of MW-scale vanadium redox flow battery for improved technoeconomic analysis

A novel data-driven, machine-learning-based method for modeling the instantaneous power losses of a distribution-sited 2 MW/8MWh vanadium redox flow battery (VRFB), a grid-scale electrochemical storage technology, is introduced and compared against benchmark empirical modeling approaches, including symmetric and asymmetric models, as well as a recent convex hull modeling approach. The novel loss modeling method introduces several advantages over the benchmark models and over simplistic efficiency estimates, the most significant of which is that the model can accurately reflect the stepwise and non-linear parasitic losses associated with the duty cycles of mechanical auxiliary systems like pump motor drives and blower fans. Residuals of the models are compared; the proposed data driven model features significantly improved accuracy over the benchmark models. The model's coefficient of determination is also improved relative to that of the benchmark models. Furthermore, a novel method for visualization of operational efficiency of the grid-scale storage technology is introduced. To demonstrate the benefits of the novel data-driven method for modeling the VRFB, the benchmark models and the proposed models are embedded into an Open DSS distribution network model to study two applications of the grid-scale electrical storage system: load leveling for grid support and energy arbitrage. This article demonstrates that the accuracy of the instantaneous power loss model significantly impacts the understanding of the state of charge of the VRFB. In turn, the accuracy of the efficiency modeling of the VRFB impacts the understanding of the potential economic value and technical benefits to the distribution network operators. In conclusion, the presented power loss modeling approach is, therefore, highly relevant for utility-stakeholders, battery asset owners, system engineers, system designers, and financial planners interested in evaluating or optimizing the operation of grid-scale VRFBs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2014

SU-E-T-379: Concave Approximations of Target Volume Dose Metrics for Intensity- Modulated Radiotherapy Treatment Planning

Purpose: The widely used treatment plan metric Dx (mimimum dose to the hottest x% by volume of the target volume) is simple to interpret and use, but is computationally poorly behaved (non-convex), this impedes its use in computationally efficient intensity-modulated radiotherapy (IMRT) treatment planning algorithms. We therefore searched for surrogate metrics that are concave, computationally efficient, and accurately correlated to Dx values in IMRT treatment plans. Methods: To find concave surrogates of D95—and more generally, Dx values with variable x values—we tested equations containing one or two generalized equivalent uniform dose (gEUD) functions. Fits were obtained by varying gEUD ‘a’ parameter values, as well as the linear equation coefficients. Fitting was performed using a dataset of dose-volume histograms from 498 de-identified head and neck IMRT treatment plans. Fit characteristics were tested using a crossvalidation process. Reported root-mean-square error values were averaged over the cross-validation shuffles. Results: As expected, the two-gEUD formula provided a superior fit, compared to the single-gEUD formula. The best approximation uses two gEUD terms: 16.25 x gEUD[a=0.45] – 15.30 x gEUD[a=1.75] – 0.69. The average root-mean-square error on repeated (70/30) cross validation was 0.94 Gy. In addition, a formula was found that reasonably approximates Dx for x between 80% and 96%. Conclusion: A simple concave function using two gEUD terms was found that correlates well with PTV D95s for these head and neck treatment plans. More generally, a formula was found that represents well the Dx for x values from 80% to 96%, thus providing a computationally efficient formula for use in treatment planning optimization. The formula may need to be adjusted for other institutions with different treatment planning protocols. We conclude that the strategy of replacing Dx values with gEUD-based formulas is promising.

60 APPLIED LIFE SCIENCES↗

NASA NTRS2003

Long Focal Length Large Mirror Fabrication System

The goals of this ambitious program are (1) to develop systems to make large superpolished optical mirrors, (2) to develop low scatter polishing techniques using centrifugal elutriation, (3) to develop a means of measuring scatter at any point on the mirror, (4) to polish a Hindle sphere to measure the optical figure of a one meter diameter convex mandrel, and (5) to fabricate low scatter, large adaptive optic graphite filled, cyanate ester replica transfer mirrors using these mandrels. Deliverables are a 30 cm diameter superpolished composite AO mirror. We fabricated a 1/3rd meter superpolished zerodur flat mandrel and with the support of our major subcontractor, Composite Mirror Applications Inc (CMA) we have demonstrated a 30 cm lightweight cyanate ester mirror with an rms microroughness between 0.6 and 0.8 nm and 8 faceplate influence function of 5 cm. The influence function was chosen to be comparable to the atmospheric correlation coefficient r(sub 0) which is about 5 cm at sea level. There was no print-thru of the graphite fibers in the cyanate ester surface (the bane of many previous efforts to use cyanate ester mirrors). Our subcontractor has devised a means for developing a 30-50 nm thick layer of graphite free pure ester resin on the surface of the mirrors. This graphite fiber filled material has a thermal expansion coefficient in the 10(exp -8) centimeter per Kelvin range (the same range of expansion coefficient as Zerodur and ULE glasses) and does not take up water and swell, so it is a nearly ideal mirror material in these areas. Unfortunately for these 0.8mm thick faceplates, the number of plies is not enough to result in isometric coverage. Isolated figure irregularities can appear, making it necessary to go to thicker faceplates. The influence function will then only approximate the length of r(sub 0), at higher altitudes or longer wavelengths. The influence function goes as the cube of the thickness, so we are now making a faceplate optimized for an infrared wavelength, possibly that used for the AO system of the Keck telescope, instead of 0.63 microns. We have polished a 55 cm diameter mandrel to better than 1/20th wave optical figure in the visible using centrifugal elutriation. CMA has just told us that it needs to retool to get optimum mirror faceplate quality in this size, so implementing the 55 cm AO mirror may be delayed somewhat. We expect to complete our 1/3 rd meter AO mirror on time using novel piezoelectric actuators with a throw of one micrometer per volt, as compared to 0.005 micrometers per volt for conventional piezoelectric actuators. We will then demonstrate its AO performance interferometrically.

Bennett, H. E.↗

DOE OSTI2026

Bio-inspired alula-based winglet design for enhanced heat transfer in high temperature fin-and-tube heat exchangers

Fin-and-tube heat exchangers (FTHEs) are widely used for high-temperature flue-gas heat recovery, but their performance is often limited by wake regions and non-uniform fin-surface temperatures. This study proposes and numerically evaluates four bio-inspired longitudinal vortex generator (VG) configurations in a high-temperature FTHE with flue-gas inlet temperature ∼1230 K: double-delta, curved double-delta, alula, and a new curved-alula geometry. The reference fin is not hydraulically plain; it already incorporates leading-edge separation columns and convex protrusions, so the alula-type winglets are assessed as downstream add-ons acting on a strongly disturbed flow. In a second step, perforations (one, two and three circular holes) are introduced into the curved-alula VGs to further tailor the flow field. Three-dimensional simulations with the Shear Stress Transpor (SST) $k - ω$ model, temperature-dependent flue-gas properties and conjugate conduction are carried out for gas-side Reynolds numbers $Re_g ≈ 8.0$ x $10^2 - 3.6$ x $10^3$ (mass flow rates 0.5 – 2.5 g/s), and the designs are compared in terms of surface heat flux, Nusselt number, friction factor and hydrothermal performance factor (HTPF). For this already-promoted fin, the additional downstream winglets provide moderate, incremental hydrothermal gains. At the highest Reynolds number, the best non-perforated design (curved-alula) increases surface heat flux from 1630.9 to 1794.7 kW/m² (∼ 10 % gain) and the Nusselt number from 227.6 to 242.6 (∼ 7 % gain), while the friction factor rises from 0.26 to about 0.30, yielding HTPF values close to unity (∼ 0.9 – 1.0). Introducing circular perforations into the curved-alula winglets acts mainly as a wake-bleeding refinement: the three-hole configuration provides a heat flux of 1824.7 kW/m² and a pressure drop of 127.9 Pa, with HTPF in the range ∼ 1.03 – 1.14 and a small (∼ 1 – 3 %) improvement over the solid curved-alula design. Flow-field analysis shows that the perforated curved-alula VGs shrink tube-wake regions, thin the thermal boundary layer and homogenize the fin-surface temperature (outlet-gas temperature ∼ 510 – 520 K and fin-surface temperature ∼ 420 – 421 K for the three-hole case). An optimal flue-gas mass flow rate of ∼ 1 g/s ($Re_g ≈ 1.5$ x $10^3$) is identified, beyond which additional heat-transfer gains are offset by rapidly increasing pressure losses. Overall, the results highlight that initial fin geometry and VG placement are as important as VG shape: alula-based winglets are expected to yield larger relative gains on simpler flat-fin layouts or when positioned closer to the fin leading edge and tube

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRSDate not supplied

Rollable Corrugated Tubular Booms

A promising candidate for a new class of deployable composite boom is the COrrugated ROllable TUbular Boom (COROTUB), which is to be employed on future large space structures by the National Aeronautics and Space Administration (NASA). The COROTUB’s two corrugated thin shells form a closed section, which yields high bending and torsional stiffness, allowing for high dimensional stability. The corrugation geometry that dictates the boom cross-section shape was completely defined by closed-form analytical equations given a set of key design parameters. Parametric studies of two geometric variables were used to evaluate which parameters most influence the cross-section’s area moment of inertia and torsional constant. The number of corrugations, flattened height of the boom, and concave and convex corrugation radii were varied to study the complete boom design space. Deeper corrugations increased the local buckling (crimpling) loads of the thin-shell boom but reduced the area moment of inertia and therefore the global buckling load of the beam-column member, and thus a fine balance was sought. Finite element analyses of pure bending and compression load cases for short, 1.22-m-long fully-deployed booms ofeither the same size (flattened height) or similar inertia properties, were first carried out to efficiently study this design trade-off in detail.A second computational study evaluated the boom transition shape and the structural response of a 6-m-long partially-deployed boom, where the effects of partially bracing the root of the boom were also assessed. It was found that bracing length could greatly affect the strength and stiffness of this boom.Bracing type (as in the discrete number of clamps at the root) had no effect on stiffness, but particularly impacted strength for loading cases and boundary conditions that promoted the onset of local web buckling at the root. The corrugated designs were also found to improve, by as much as 20%, the boom bending and axial strength and boom transition length from flat to deployed when compared to a state-of-the-art boom alternative of the same size, for short bracing distances, such as those that would be required by small satellite applications. In addition, they showed a comparable structural response for longer bracing supports.

Booms↗