Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shape 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 541 records · Page 30

Slot-coupled CW standing wave accelerating cavity

A slot-coupled CW standing wave multi-cell accelerating cavity. To achieve high efficiency graded beta acceleration, each cell in the multi-cell cavity may include different cell lengths. Alternatively, to achieve high efficiency with acceleration for particles with beta equal to 1, each cell in the multi-cell cavity may include the same cell design. Coupling between the cells is achieved with a plurality of axially aligned kidney-shaped slots on the wall between cells. The slot-coupling method makes the design very compact. The shape of the cell, including the slots and the cone, are optimized to maximize the power efficiency and minimize the peak power density on the surface. The slots are non-resonant, thereby enabling shorter slots and less power loss.

43 PARTICLE ACCELERATORS↗

Compliant heat exchangers, heat pipes and methods for making same

A method for fabricating heat exchangers using additive manufacturing technologies. Additive manufacturing enables the manufacture of heat exchangers with complex geometries and/or with internal and external integral surface features. Additive manufacture also facilitates the manufacture of heat exchangers with regional variations, such as changes in size, shape and surface features. In one embodiment, the present invention provides a heat exchanger with a helicoidal shape that provides axial elastic compliance. In one embodiment, the internal channel of the heat exchanger varies along its length. The internal channel may have a cross-sectional area that increases progressively from one end to the other. In one embodiment, the external shape of the tubular structure may be non-circular to optimize heat transfer with an external heat transfer fluid. In one embodiment, the present invention provides a heat pipe with an internal wicking structure formed as an integral part of the additive manufacturing process.

Lara-Curzio, Edgar↗

Effect of Powder Shape and Size on Rheological, Thermal, and Segregation Properties of Low-Pressure Powder Injection Molding Feedstocks

Optimization of molding parameters at different filling stages requires that the characteristics of low-viscosity feedstocks must be properly known. In this study, the influence of the particle size (3, 7, and 12 μm) and powder morphology (gas- and water-atomized) was quantified using different feedstock formulations based on a 17-4 PH stainless steel powder, all combined with the same binder system. The specific heat capacity, viscosity, and segregation levels were assessed using a differential scanning calorimeter, a rotational rheometer, and a thermogravimetric analyzer, respectively. It was then shown that the feedstock prepared with the gas-atomized powder exhibits a higher specific heat capacity value when compared to a water-atomized powder (i.e., Cp varying from 0.6 to 0.3 J/g K, respectively). It was also shown that the feedstock viscosity profiles and the intensity of segregation both depend significantly on the powder shape and size used in the feedstock formulation. For short processing times (e.g., < 1 min spent in molten state), a feedstock formulated with coarse gas-atomized powder can be considered to be the best candidate for an injection process because its viscosity and segregation potential are relatively low (with a viscosity varying from 0.5 to 2 Pa s, and a volume fraction of powder close to 60 vol.%). Following this, the moldability index was used to predict that coarser water-atomized powder (i.e., 12 versus 7 μm) will not have a significant effect on the molding properties, and can be seen as a potential means of decreasing shrinkage during sintering without impacting the injection properties. For very long processing times, the magnitude of segregation remains insignificant for feedstock formulated with fine powder (i.e., a volume fraction of powder remaining close to 60 vol.%). However, this important gain (i.e., no segregation) comes with a significant increase in viscosity varying from 2 to 10 Pa s, resulting in a decrease in the molding properties.

36 MATERIALS SCIENCE↗

Pushing the limits of pulse shape discrimination in a large liquid xenon detector

Abstract The LUX-ZEPLIN (LZ) experiment is a direct-detection dark matter experiment, optimized to search for weakly interacting massive particles (WIMPs) through WIMP-nucleon interactions. The main challenge in dark matter detection is differentiating between WIMP signals and background events. In LZ, the ratio of ionization to scintillation signals (charge-to-light) is the primary method for rejecting electronic recoil (ER) background. Pulse shape discrimination (PSD) offers a method for additional ER backgrounds rejection in liquid xenon detectors. In this paper, the discrimination power of PSD with the LZ experiment is discussed. To precisely characterize the scintillation pulse shape, an analysis framework is developed to reconstruct the detection time of individual photons. Using LZ calibration data, the photon-timing tail fraction discriminator is optimized and achieves ER leakage as low as $$15\%$$ 15 % . For specific background processes such as $$^{124}$$ 124 Xe double electron capture, the leakage is reduced further to about $$5\%$$ 5 % . PSD is combined with charge-to-light to form two-factor discrimination (TFD). The optimized TFD performance is compared with the performance of the charge-to-light method, with the corresponding false positive rate reduced by up to a factor of two for large scintillation pulses. Finally, PSD and TFD are applied to data from LZ’s WS2024 run and their performance is summarized.

Akerib, D. S. [SLAC National Accelerator Laborator↗

Controlling circuitry underlies the growth optimization of Saccharomyces cerevisiae

Microbial growth emerges from coordinated synthesis of various cellular components from limited resources. In Saccharomyces cerevisiae, cyclic AMP (cAMP)-mediated signaling is shown to orchestrate cellular metabolism; however, it remains unclear quantitatively how the controlling circuit drives resource partition and subsequently shapes biomass growth. Here we combined experiment with mathematical modeling to dissect the signaling-mediated growth optimization of S. cerevisiae. We showed that, through cAMP-mediated control, the organism achieves maximal or nearly maximal steady-state growth during the utilization of multiple tested substrates as well as under perturbations impairing glucose uptake. However, the optimal cAMP concentration varies across cases, suggesting that different modes of resource allocation are adopted for varied conditions. Under settings with nutrient alterations, S. cerevisiae tunes its cAMP level to dynamically reprogram itself to realize rapid adaptation. Moreover, to achieve growth maximization, cells employ additional regulatory systems such as the GCN2-mediated amino acid control. This study establishes a systematic understanding of global resource allocation in S. cerevisiae, providing insights into quantitative yeast physiology as well as metabolic strain engineering for biotechnological applications.

59 BASIC BIOLOGICAL SCIENCES↗

Data for Controlling Circuitry Underlies the Growth Optimization of Saccharomyces cerevisiae

Microbial growth emerges from coordinated synthesis of various cellular components from limited resources. In Saccharomyces cerevisiae , cyclic AMP (cAMP)-mediated signaling is shown to orchestrate cellular metabolism; however, it remains unclear quantitatively how the controlling circuit drives resource partition and subsequently shapes biomass growth. Here we combined experiment with mathematical modeling to dissect the signaling-mediated growth optimization of S. cerevisiae . We showed that, through cAMP-mediated control, the organism achieves maximal or nearly maximal steady-state growth during the utilization of multiple tested substrates as well as under perturbations impairing glucose uptake. However, the optimal cAMP concentration varies across cases, suggesting that different modes of resource allocation are adopted for varied conditions. Under settings with nutrient alterations, S. cerevisiae tunes its cAMP level to dynamically reprogram itself to realize rapid adaptation. Moreover, to achieve growth maximization, cells employ additional regulatory systems such as the GCN2-mediated amino acid control. This study establishes a systematic understanding of global resource allocation in S. cerevisiae , providing insights into quantitative yeast physiology as well as metabolic strain engineering for biotechnological applications.

Conversion↗

Fast automatic 3D liver segmentation based on a three-level AdaBoost-guided active shape model

Purpose: A robust, automatic, and rapid method for liver delineation is urgently needed for the diagnosis and treatment of liver disorders. Until now, the high variability in liver shape, local image artifacts, and the presence of tumors have complicated the development of automatic 3D liver segmentation. In this study, an automatic three-level AdaBoost-guided active shape model (ASM) is proposed for the segmentation of the liver based on enhanced computed tomography images in a robust and fast manner, with an emphasis on the detection of tumors. Methods: The AdaBoost voxel classifier and AdaBoost profile classifier were used to automatically guide three-level active shape modeling. In the first level of model initialization, fast automatic liver segmentation by an AdaBoost voxel classifier method is proposed. A shape model is then initialized by registration with the resulting rough segmentation. In the second level of active shape model fitting, a prior model based on the two-class AdaBoost profile classifier is proposed to identify the optimal surface. In the third level, a deformable simplex mesh with profile probability and curvature constraint as the external force is used to refine the shape fitting result. In total, three registration methods—3D similarity registration, probability atlas B-spline, and their proposed deformable closest point registration—are used to establish shape correspondence. Results: The proposed method was evaluated using three public challenge datasets: 3Dircadb1, SLIVER07, and Visceral Anatomy3. The results showed that our approach performs with promising efficiency, with an average of 35 s, and accuracy, with an average Dice similarity coefficient (DSC) of 0.94 ± 0.02, 0.96 ± 0.01, and 0.94 ± 0.02 for the 3Dircadb1, SLIVER07, and Anatomy3 training datasets, respectively. The DSC of the SLIVER07 testing and Anatomy3 unseen testing datasets were 0.964 and 0.933, respectively. Conclusions: The proposed automatic approach achieves robust, accurate, and fast liver segmentation for 3D CTce datasets. The AdaBoost voxel classifier can detect liver area quickly without errors and provides sufficient liver shape information for model initialization. The AdaBoost profile classifier achieves sufficient accuracy and greatly decreases segmentation time. These results show that the proposed segmentation method achieves a level of accuracy comparable to that of state-of-the-art automatic methods based on ASM.

60 APPLIED LIFE SCIENCES↗

Formulation of Conductive Inks for DIW Printed Battery Electrodes

Ink formulations for DIW printing of Li ion battery electrodes were optimized for high conductivity, surface area and printability. Our focus is particularly on optimizing Li-ion battery electrodes for low temperature applications by 3-D printing the cathodes into high surface area shapes. The sluggish kinetics of ion transport at the solid-electrolyte phase boundary can be a major problem in batteries at low temperatures, so optimizing for ion transport can help combat these issues. Two families of DIW ink were investigated; aqueous conductive polymer suspension with graphene nano platelets (GNP) ink, and “oil in water” graphene oxide (GO) micro emulsion ink which carbonizes to an open cell graphene aerogel. By iterative formulation changes, the conductivity of PEDOT:PSS based ink was increased over tenfold, and stiffness of GO emulsion ink was increased to provide better printability without sacrificing porosity.

36 MATERIALS SCIENCE↗

Integrating novel stellarator single-stage optimization algorithms to design the Columbia stellarator experiment

Abstract The Columbia Stellarator eXperiment (CSX), currently being designed at Columbia University, aims to test theoretical predictions related to QA plasma behavior, and to pioneer the construction of an optimized stellarator using three-dimensional, non-insulated high-temperature superconducting (NI-HTS) coils. The magnetic configuration is generated by a combination of two circular planar poloidal field (PF) coils and two 3D-shaped interlinked (IL) coils, with the possibility to add windowpane coils to enhance shaping and experimental flexibility. The PF coils and vacuum vessel are repurposed from the former Columbia Non-Neutral Torus experiment, while the IL coils will be custom-wound in-house using NI-HTS tapes. To obtain a plasma shape that meets the physics objectives with a limited number of coils, novel single-stage optimization techniques are employed, optimizing both the plasma and coils concurrently, in particular targeting a tight aspect ratio QA plasma and minimized strain on the HTS tape. Despite the increased complexity due to the expanded degrees of freedom, these methods successfully identify optimized plasma geometries that can be realized by coils meeting engineering specifications. This paper discusses the derivation of the constraints and objectives specific to CSX, and describe how two recently developed single-stage optimization methodologies are applied to the design of CSX. A set of selected configurations for CSX is then described in detail.

Baillod, A. (ORCID:0000000303529180)↗

TH-AB-209-01: Making Benchtop X-Ray Fluorescence Computed Tomography (XFCT) Practical for in Vivo Imaging by Integration of a Dedicated High-Performance X-Ray Source in Conjunction with Micro-CT Functionality

Purpose: To make benchtop x-ray fluorescence computed tomography (XFCT) practical for routine preclinical imaging tasks with gold nanoparticles (GNPs) by deploying, integrating, and characterizing a dedicated high-performance x-ray source and addition of simultaneous micro-CT functionality. Methods: Considerable research effort is currently under way to develop a polychromatic benchtop cone-beam XFCT system capable of imaging GNPs by stimulation and detection of gold K-shell x-ray fluorescence (XRF) photons. Recently, an ad hoc high-power x-ray source was incorporated and used to image the biodistribution of GNPs within a mouse, postmortem. In the current work, a dedicated x-ray source system featuring a liquid-cooled tungsten-target x-ray tube (max 160 kVp, ∼3 kW power) was deployed. The source was operated at 125 kVp, 24 mA. The tube’s compact dimensions allowed greater flexibility for optimizing both the irradiation and detection geometries. Incident x-rays were shaped by a conical collimator and filtered by 2 mm of tin. A compact “OEM” cadmium-telluride x-ray detector was implemented for detecting XRF/scatter spectra. Additionally, a flat panel detector was installed to allow simultaneous transmission CT imaging. The performance of the system was characterized by determining the detection limit (10-second acquisition time) for inserts filled with water/GNPs at various concentrations (0 and 0.010–1.0 wt%) and embedded in a small-animal-sized phantom. The phantom was loaded with 0.5, 0.3, and 0.1 wt% inserts and imaged using XFCT and simultaneous micro-CT. Results: An unprecedented detection limit of 0.030 wt% was experimentally demonstrated, with a 33% reduction in acquisition time. The reconstructed XFCT image accurately localized the imaging inserts. Micro-CT imaging did not provide enough contrast to distinguish imaging inserts from the phantom under the current conditions. Conclusion: The system is immediately capable of in vivo preclinical XFCT imaging with GNPs. Micro-CT imaging will require optimization of irradiation parameters to improve contrast. Supported by NIH/NCI grant R01CA155446; This investigation was supported by NIH/NCI grant R01CA155446.

60 APPLIED LIFE SCIENCES↗

Knowledge Oriented Graph Unified Transformer (KOGUT) v0.1

KOGUT — Knowledge Oriented Graph Unified Transformer KOGUT implements the Relational Graph Transformer (RelGT) architecture for knowledge graph link prediction in biological domains, with a primary focus on microbial growth media prediction. While the original RelGT (arXiv:2505.10960) targets relational tables, time series, and multi-table databases, KOGUT adapts this architecture for heterogeneous biological knowledge graphs, providing first-in-class AI predictive models for microbial cultivation. Key Adaptations Beyond Original RelGT: - Knowledge Graph Focus: Applied to biological KGs with semantic node types (taxa, chemicals, media, phenotypes, environments) versus generic relational database tables, trained on the KG-Microbe knowledge graph (1.3M entities, 2.9M edges, 24 relation types). - Multimodal Node Encoding: Integrates node labels, categories, descriptions, and synonyms from KG metadata through learned embedding layers—adapting relational column features to graph node attributes with textual semantics. - Extended K-Hop Subgraph Strategy: Optimized neighborhood sampling (3-hop default, configurable up to 200 nodes) tuned for sparse biological networks, building on the original local-global attention framework with biological relation preservation. - Biolink Predicate Preservation: Type-specific transformations for 24 biological edge semantics (occurs_in, consumes, produces, has_phenotype, subclass_of) beyond standard relational foreign keys, enabling multi-relation link prediction. - Inductive Learning Support: Enables zero-shot predictions for novel taxa through feature-based embeddings (temperature, oxygen requirements, gram stain, cell shape), extending the original transductive relational benchmark scope to uncultured microorganisms. CheapSOTA Performance Optimizations (This Distribution): - VQ-EMA Centroid Attention: Vector quantization with exponential moving average for improved global context modeling (+5-10% MRR improvement). - HDF5 Precomputed Data Loading: One-time preprocessing of k-hop subgraphs to eliminate redundant graph traversals (2-5× training speedup). - Distributed Data Parallel Training: Multi-GPU support for scaling to larger knowledge graphs (tested on 4× NVIDIA A100 GPUs at NERSC Perlmutter). - Mixed Precision Training: Automatic mixed precision (AMP) for memory efficiency and faster training. Advantages Over Standard Knowledge Graph Embedding Models: Combines RelGT's proven multi-element tokenization (features, type, hop, structure) with graph-native biological representations, enabling interpretable link prediction across heterogeneous entities that standard embedding models (TransE, RotatE, ComplEx) and table-based transformers cannot directly model. Achieves near-perfect performance on microbial growth media prediction (MRR: 0.9966, Precision@1: 0.9932, Hit@10: 1.0000) while maintaining explainability through attention-based reasoning over biological pathways. Training Data: - KG-Microbe merged knowledge graph: 1,379,337 nodes, 2,960,472 edges - 24 biological relation types including taxonomic hierarchies, metabolic interactions, phenotype associations, and environmental relationships - Primary prediction task: Growth media suitability for microbial taxa (biolink:occurs_in, 50K edges) - Multi-relation capability: Predicts links for any of the 24 relation types, including chemical consumption/production, phenotype associations, and taxonomic classification Citation: Original RelGT Architecture: Dwivedi et al., "Relational Graph Transformer", arXiv:2505.10960, 2025 KOGUT Implementation: Knowledge Oriented Graph Unified Transformer for Microbial Growth Media Prediction Developed at Lawrence Berkeley National Laboratory (LBNL) Trained on NERSC Perlmutter supercomputer

Joachimiak, Marcin [Lawrence Berkeley National Lab↗

CRNT4SBML: a Python package for the detection of bistability in biochemical reaction networks

Motivation: Signaling pathways capable of switching between two states are ubiquitous within living organisms. They provide the cells with the means to produce reversible or irreversible decisions. Switchlike behavior of biological systems is realized through biochemical reaction networks capable of having two or more distinct steady states which are dependent on initial conditions. Investigation of whether a certain signaling pathway can confer bistability involves a substantial amount of hypothesis testing. The cost of direct experimental testing can be prohibitive. Therefore, constraining the hypothesis space is highly bene?cial. One such methodology is based on Chemical Reaction Network Theory which uses computational techniques to rule out pathways that are not capable of bistability regardless of kinetic constant values and molecule concentrations. Although useful, these methods are complicated from both pureandcomputationalmathematicsperspectives.Thus,theiradoptionisverylimitedamongstbiologists. Results: We brought Chemical Reaction Network Theory approaches closer to experimental biologists by automating all the necessary steps in CRNT4SMBL. The input is based on SBML format which is the community standard for biological pathway communication. The tool parses SBML and derives C-graph representations of the biological pathway with mass action kinetics. Next steps involve an ef?cient search for potential saddle-node bifurcation points using an optimization technique. This type of bifurcation is important as it has the potential of acting as a switching point between two steady states. Finally, if any bifurcation points are present, numerical continuation analysis extends the equilibria branches for generating the diagram. Presence of an S-shaped bifurcation diagram indicates that the pathway acts as a bistable switch for the given optimization parameters. Availability: CRNT4SBML is available via the Python Package Index. The documentation can be found at https://crnt4sbml.readthedocs.io. CRNT4SBML is licensed under the Apache Software License 2.0. Contact: vladislav.petyuk@pnnl.gov

Reyes, Brandon C.↗

RL-HEMS: Reinforcement Learning Based Home Energy Management System for HVAC Energy Optimization

Heating, ventilation, and air conditioning (HVAC) is one of the major energy consumers in the residential sector. It is important to be able to monitor and control the energy consumed to provide utility services such as load shaping while satisfying the comfort and economic constraints of the homeowner. The objective of this work is to create the optimal schedule for HVAC operation to reduce the cost while satisfying the home-owner and equipment’s constraints using a model-free Reinforcement Learning (RL)-based optimization. The specific goal is to find the right balance between reducing energy cost, consumption, and customer comfort level. Our research effort addresses this optimization problem using multiple components: the development of initial learning test-bed and implementation of RL techniques on a real home. This will enable the rapid evaluation of the RL techniques and provide an early baseline to train before implementation on site. The RL algorithm is designed to learn the energy use patterns and generate the optimized schedule for HVAC within an acceptable time-interval to satisfy the homeowner’s comfort and minimize the energy usage. Our preliminary results are promising and we have observed a 17% reduction in the total cost and a 15% reduction in the power utilization using our RL-based HVAC model—RL-HEMS.

Kotevska, Olivera↗

Acoustic transfer of protein crystals from agarose pedestals to micromeshes for high-throughput screening

An acoustic high-throughput screening method is described for harvesting protein crystals and combining the protein crystals with chemicals such as a fragment library. Acoustic droplet ejection (ADE) is an emerging technology with broad applications in serial crystallography such as growing, improving and manipulating protein crystals. One application of this technology is to gently transfer crystals onto MiTeGen micromeshes with minimal solvent. Once mounted on a micromesh, each crystal can be combined with different chemicals such as crystal-improving additives or a fragment library. Acoustic crystal mounting is fast (2.33 transfers s{sup −1}) and all transfers occur in a sealed environment that is in vapor equilibrium with the mother liquor. Here, a system is presented to retain crystals near the ejection point and away from the inaccessible dead volume at the bottom of the well by placing the crystals on a concave agarose pedestal (CAP) with the same chemical composition as the crystal mother liquor. The bowl-shaped CAP is impenetrable to crystals. Consequently, gravity will gently move the crystals into the optimal location for acoustic ejection. It is demonstrated that an agarose pedestal of this type is compatible with most commercially available crystallization conditions and that protein crystals are readily transferred from the agarose pedestal onto micromeshes with no loss in diffraction quality. It is also shown that crystals can be grown directly on CAPs, which avoids the need to transfer the crystals from the hanging drop to a CAP. This technology has been used to combine thermolysin and lysozyme crystals with an assortment of anomalously scattering heavy atoms. The results point towards a fast nanolitre method for crystal mounting and high-throughput screening.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Thermodynamic analysis of a MEMS based differential scanning calorimeter model

Thin film based differential scanning calorimeter (DSC) was developed and utilized for the study of phase-transition kinetics in protein samples. In this paper, we proposed a MEMS based DSC model. The distance between the sample and reference cell was optimized based on the thermal analysis result. The influence of the response time to the shape of the DSC thermograms and the transition temperature was also studied systematically. Meanwhile, Lumry-Eyring model was used to investigate the relationship between the scanning rate and the transition temperature. The modeling result was verified by the lysozyme denaturation tests. This study provided a pathway for the optimization of MEMS based DSC array and the correction of the DSC thermograms.

47 OTHER INSTRUMENTATION↗

Survey on stochastic distribution systems: A full probability density function control theory with potential applications

Complex systems seen either in general engineering practice or economics are subjected to ever increased uncertainties that are mostly represented as random variables or parameters, and the characteristics of random variables are represented by their probability density functions (PDFs). Controlling their PDFs means to shape their stochastic distributions and in general it would provide a full treatment for system analysis and operational control and optimization. This leads to the development of stochastic distribution control (SDC) systems theory in the past decades, where the original aim of the controller design is to realize a shape control of the distributions of certain random variables in their PDFs sense for some engineering processes. Indeed, once the PDFs of these random variables or parameters are used to describe their distribution characters, the control task is to obtain control signals so that the output PDFs of stochastic systems are made to follow their target PDFs. The subject of SDC was initially originated for non-Gaussian stochastic control systems design but has found a wide spectrum of applications in general systems in terms of data-driven modeling, analysis, signal processing (filtering), data mining via multivariable statistics, decision-making (optimization) for systems subjected to uncertainties and even in economics. In this context, SDC constitutes an effective primer tool for complex system analysis, control and operational optimizations. In this review paper, a detailed survey of the developments on the research of SDC systems will be made together with their wide spectrum applications and future perspectives.

42 ENGINEERING↗

Design considerations of highly efficient D-shaped plasmonic biosensor

Design considerations of photonic crystal fiber (PCF) surface plasmon biosensor with high sensitivity to monitor glucose concentration is reported and studied. Based on a well-known large mode area (LMA) single mode PCF, two different configurations have been studied to investigate the impact of the etching process on the biosensor sensitivity. Furthermore, the possibility of infiltrating one hole by a plasmonic material to increases the biosensor sensitivity is studied. A full vectorial finite element method (FVFEM) is used to carry out the analysis of the reported biosensor. In addition, the geometrical parameters of the suggested biosensors have been optimized to obtain the highest sensitivity. The suggested biosensor based on a D-shaped PCF with plasmonic rod achieves sensitivity as high as 13,600 nm/RIU with corresponding resolution of 7.35×10{sup −6} RIU.The analysis also reveals that the proposed biosensor has a linear performance which is needed practically. Therefore, the reported biosensor has advantages in terms of fabrication feasibility and high linear sensitivity.

36 MATERIALS SCIENCE↗

Stochastic Programming Approach for Resource Selection under Demand Uncertainty

Cost-efficient selection and scheduling of a subset of geographically distributed resources to meet the demands of a scientific workflow is a challenging problem. The problem is exacerbated by uncertainties in demand and availability of resources. In this paper, we present a stochastic optimization based framework for robust decision making in the selection of distributed resources over a planning horizon under demand uncertainty. We present a novel two-stage stochastic programming model for resource selection, and implement an L-shaped decomposition algorithm to solve this model. A Sample Average Approximation algorithm is integrated to enable stochastic optimization to solve problems with a large number of scenarios. Using the metric of stochastic solution, we demonstrate up to 30% cost reduction relative to solutions without explicit consideration of demand uncertainty for a 24-month problem. We also demonstrate up to 54% cost reduction relative to a previously developed solution for a 36-month problem. We further argue that the composition of resources selected is superior to solutions computed without explicit consideration of uncertainties. Given the importance of resource selection and scheduling of complex scientific workflows, especially in the context of commercial cloud computing, we believe that our novel stochastic programming framework will benefit many researchers as well as users of distributed computing resources.

scientific workflows, task scheduling, optimizatio↗