Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variability study”

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 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Intra-annual variation in microclimatic conditions in relation to vegetation type and structure in two tropical dry forests undergoing secondary succession

Microclimate acts as a strong filter on species performance in restored and regenerating forests, particularly in seasonally dry tropical forests (SDTF). Yet few studies have measured microclimate patterns across succession in SDTF. Furthermore, although dynamic vegetation models simulate microclimate, evaluation of these simulated variables with field observations has been relatively uncommon. Here, we investigated the seasonal patterns of soil temperature and soil water in naturally regenerated and planted successional vegetation in SDTF in Costa Rica and Puerto Rico, using complementary approaches of intensive field observations and simulation modeling with the Ecosystem Demography model. We found that plots representing later successional stages were wetter on average, but only during the dry season. During the wet season, mean soil water did not differ across vegetation types, but open, early successional vegetation experienced more frequent extreme wet and dry conditions than older forest and plantations. Soil temperature tended to decline with forest structure, and later successional vegetation also experienced less extreme daily temperature fluctuations. Basal area and leaf area index were the best predictors of differences in soil water and temperature across plots. Model simulations were consistent with observations of wet season soil temperature and soil water, but the model failed to reproduce dry season soil moisture dynamics, suggesting that further work is needed to reduce model biases in microclimate variables. Collectively, our results imply that common assumptions about how microclimates influence successional processes in SDTF should be revisited.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Metal additively manufactured wavy fin cold-plate architecture for improved thermal-hydraulic performance

Rapid growth in artificial intelligence and data center workloads demands high-performance liquid cooling to manage increasing chip power. This study presents two metal-additive-manufactured cold plates with sinusoidal fins, constant-amplitude wavy fins and linearly variable-amplitude wavy fins and compares them against metal-additive-manufactured straight fins using experiments conducted at 1 kW heat dissipation as well as high-fidelity 3D conjugate computational fluid dynamic simulations. The cold plates were printed in AlSi10Mg material and underwent design using a Python-automated workflow prior to manufacture and testing. The experiments show that wavy fins reduce the normalized thermal resistance by 35 to 45 % at water flow rates from 1 to 4 LPM. At a fixed 20 kPa pressure drop, the variable-waviness design lowered peak surface temperature by 9 °C and thermal resistance by 51 %, while edge-channel maldistribution in the constant wavy fin design limited gains. A thermal resistance breakdown revealed that 55–63 % of the total thermal resistance in wavy designs comes from base heat conduction, 27–33 % from fin heat conduction, and 9–13 % from fin heat convection, indicating the need to address conduction bottlenecks. Parametric sweeps identify a 3 mm fin pitch as optimal, and that horizontal inlet/outlet manifolds further reduce pressure drop by 30–60 % and thermal resistance by 9–16 % relative to vertical inlet-outlet manifolds. The results yield comprehensive guidelines for fin geometry, manifold alignment, material selection and additive-manufacturing constraints to realize high-performance liquid-cooled cold plates for power-dense electronics.

3d printing↗
DOE OSTI2023

Toggling Stereochemical Activity through Interstitial Positioning of Cations between 2D V 2 O 5 Double Layers

The 5/6s 2 lone-pair electrons of p-block cations in their lower oxidation states are a versatile electronic and geometric structure motif that can underpin lattice anharmonicity and often engender electronic and structural instabilities that underpin the function of active elements in nonlinear optics, thermochromics, thermoelectrics, neuromorphic computing, and photocatalysis. In contrast to periodic solids where lone-pair-bearing cations are part of the structural framework, installing lone-pair-bearing cations in the interstitial sites of intercalation hosts provides a means of a systematically modulating electronic structure through the choice of the group and the period of the inserted cation while preserving the overall framework connectivity. The extent of stereochemical activity and the energy positioning of lone-pair-derived mid-gap states depend on the cation identity, stoichiometry, and strength of anion hybridization. V 2 O 5 polymorphs are versatile insertion hosts that can accommodate a broad range of s-, p-, and d-block cations. However, the insertion of lone-pair-bearing cations remains largely underexplored. In this article, we examine the implications of varying the 6s 2 cations situated in interlayer sites between condensed [V 4 O 10 ] n double layers. Systematic modulations of lattice distortions, electronic structure, and magnetic ordering are observed with increasing strength of stereochemical activity from group 12 to group 14 cations. We compare and contrast p-block-layered M x V 2 O 5 (M = Hg, Tl, and Pb) compounds and map the significance of local off-centering arising from the stereochemical activity of lone-pair cations to the emergence of filled antibonding lone-pair 6s 2 –O 2p-hybridized mid-gap states mediated by second-order Jahn–Teller distortions. Crystallographic studies of cation coordination environments and the resulting modulation of V–V interactions have been used in conjunction with variable-energy hard X-ray photoelectron spectroscopy measurements, first-principles electronic structure calculations, and crystal orbital Hamilton population analyses to decipher the origins of stereochemical activity. Magnetic susceptibility measurements reveal antiferromagnetic signatures for all the three compounds. However, the differences in V–V interactions significantly affect the energy balance of the superexchange interactions, resulting in an ordering temperature of 160 and 260 K for Hg 0.5 V 2 O 5 and δ-Tl 0.5 V 2 O 5 , respectively, as compared to 7 K for δ-Pb 0.5 V 2 O 5 . In δ-Pb 0.5 V 2 O 5 , the strong stereochemical activity of electron lone pairs and the resulting electrostatic repulsions enforce superlattice ordering, which strongly modifies the electronic localization patterns along the [V 4 O 10 ] slabs, resulting in disrupted magnetic ordering and an anomalously low ordering temperature. The results demonstrate a versatile strategy for toggling the stereochemical activity of electron lone pairs to modify the electronic structure near the Fermi level and to mediate superexchange interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Operando X-Ray Tomography Imaging of Solid-State Electrolyte Response to Li Evolution under Realistic Operating Conditions

Solid-state Li-ion conductors are a next-generation battery technology that are particularly promising for electric vehicles, offering the capacitive benefits of Li metal anodes with nonflammable electrolytes. Microstructural evolution in these solid-state batteries, especially the Li anode, requires deeper understanding of the conditions under which certain undesired behaviors are more likely to occur. This study utilizes operando X-ray computed tomography to visualize the behavior of lithium in response to device operation, considering variables that are relevant to a realistic battery: stack pressure, microscale defects such as pores, and temperature. This work demonstrates that experimental and operational conditions, especially temperature, affect the fundamental driving forces of realistic solid electrolyte systems and provides visual insight into possible mechanisms of Li migration. In particular, the effect of these variables on Li propagation from the anode through pre-existing defects and Li nucleation within the imperfectly electron-insulating electrolyte is highlighted.

battery↗
DOE OSTI2025

Diurnal Trends and Meteorological Factors Influencing the Variability of Fluorescent Bioaerosol in Mt. Crested Butte, Colorado During SAIL

Abstract Despite the crucial role of bioaerosol particles (BAP) in our climate system, local ecosystems, and human health, our grasp on their atmospheric interactions is hampered by a lack of high‐resolution and long‐term data, which is essential for understanding their abundance and variability in response to meteorological conditions. To discern these relationships in a high‐altitude mountainous terrain is even less well‐studied. Therefore, we deployed a Wideband Integrated Bioaerosol Sensor (WIBS‐Neo) for three months during the first biologically active season of the Surface Atmosphere Integrated Field Laboratory (SAIL) campaign in summer 2022 on Mt. Crested Butte in Colorado (elevation 3,140 m ASL). Here, we report real‐time fluorescent BAP (FBAP) data collected from 15 June to 13 September 2022, sampled from within the ARM Aerosol Observing System (AOS). To our knowledge, these are the first high altitude (>3,000 m ASL) continuous measurements of FBAP made in North America. During this deployment, we observed on average 21% and as many as 48% (hourly maximum) of particles within the detection size range ( to ) of the WIBS as FBAP. Our analysis presents the diurnal cycles for seven distinct types of FBAP, showing unique patterns and, in some cases, correlations with temperature and solar radiation cycles. Abundance and composition of FBAP varied with relative humidity and precipitation. Precipitation events appeared to both cause emission and removal of FBAP, whereas dust and smoke events had no significant effect, highlighting the critical role of meteorology on FBAP at high altitudes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Intensification of future subsurface marine heatwaves in an eddy-resolving model

A shift in depth range enables marine organisms to adapt to marine heatwaves (MHWs). Subsurface MHWs could limit this pathway, yet their response to climate warming remains unclear. Here, using an eddy-resolving Earth system model forced under a high emission scenario, we project a robust global increase in subsurface MHWs driven by rising subsurface mean temperatures and enhanced temperature variability. Historically, maximum MHW intensity occurs around 100m depth, which shifts to the faster-warming surface under greenhouse warming. However, removing the long-term warming trend yields an increase in subsurface MHW intensity and annual days greater than that at the surface, especially in large marine ecosystem regions, primarily due to increased variability. Additionally, days of the surface and subsurface concurrent event increase ten times more than those of individual events. Our study highlights a heightened threat to marine organisms under global warming, as the increased subsurface heatwaves reduce their refuge options.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Variability and origins of grain boundary electric potential detected by electron holography and atom-probe tomography

A number of grain boundary phenomena in ionic materials, in particular, anomalous (either depressed or enhanced) charge transport, have been attributed to space charge effects. Developing effective strategies to manipulate transport behaviour requires deep knowledge of the origins of the interfacial charge, as well as its variability within a polycrystalline sample with millions of unique grain boundaries. Electron holography is a powerful technique uniquely suited for studying the electric potential profile at individual grain boundaries, whereas atom-probe tomography provides access to the chemical identify of essentially every atom at individual grain boundaries. Using these two techniques, we show in this paper that the space charge potential at grain boundaries in lightly doped, high-purity ceria can vary by almost an order of magnitude. We further find that trace impurities (<25 ppm), rather than inherent thermodynamic factors, may be the ultimate source of grain boundary charge. These insights suggest chemical tunability of grain boundary transport properties. A number of grain boundary phenomena in ionic materials, such as anomalous charge transport, have been attributed to space charge effects. Space charge potential at grain boundaries in lightly doped, high-purity ceria is now shown to vary by almost an order of magnitude.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Physics-based compact modeling of electro-thermal memristors: Negative differential resistance, local activity, and non-local dynamical bifurcations

Leon Chua's Local Activity theory quantitatively relates the compact model of an isolated nonlinear circuit element, such as a memristor, to its potential for desired dynamical behaviors when externally coupled to passive elements in a circuit. However, the theory's use has often been limited to potentially unphysical toy models and analyses of small-signal linear circuits containing pseudo-elements (resistors, capacitors, and inductors), which provide little insight into required physical, material, and device properties. Furthermore, the Local Activity concept relies on a local analysis and must be complemented by examining dynamical behavior far away from the steady-states of a circuit. In this work, we review and study a class of generic and extended one-dimensional electro-thermal memristors (i.e., temperature is the sole state variable), re-framing the analysis in terms of physically motivated definitions and visualizations to derive intuitive compact models and simulate their dynamical behavior in terms of experimentally measurable properties, such as electrical and thermal conductance and capacitance and their derivatives with respect to voltage and temperature. Within this unified framework, we connect steady-state phenomena, such as negative differential resistance, and dynamical behaviors, such as instability, oscillations, and bifurcations, through a set of dimensionless nonlinearity parameters. In particular, we reveal that the reactance associated with electro-thermal memristors is the result of a phase shift between oscillating current and voltage induced by the dynamical delay and coupling between the electrical and thermal variables. We thus, demonstrate both the utility and limitations of local analyses to understand non-local dynamical behavior. Critically for future experimentation, the analyses show that external coupling of a memristor to impedances within modern sourcing and measurement instruments can dominate the response of the total circuit, making it impossible to characterize the response of an uncoupled circuit element for which a compact model is desired. However, these effects can be minimized by proper understanding of the Local Activity theory to design and utilize purpose-built instruments.

Brown, Timothy D. (ORCID:0000000347068355)↗
DOE OSTI2022

Mechanical loading prediction through accelerometry data during walking and running

ABSTRACT Currently, there is no way to assess mechanical loading variables such as peak ground reaction forces (pGRF) and peak loading rate (pLR) in clinical settings. The purpose of this study was to develop accelerometry‐based equations to predict both pGRF and pLR during walking and running. One hundred and thirty one subjects (79 females; 76.9 ± 19.6 kg) walked and ran at different speeds (2–14 km·h −1 ) on a force plate–instrumented treadmill while wearing accelerometers at their ankle, lower back and hip. Regression equations were developed to predict pGRF and pLR from accelerometry data. Leave‐one‐out cross‐validation was used to calculate prediction accuracy and Bland–Altman plots. Our pGRF prediction equation was compared with a reference equation previously published. Body mass and peak acceleration were included for pGRF prediction and body mass and peak acceleration rate for pLR prediction. All pGRF equation coefficients of determination were above 0.96, and a good agreement between actual and predicted pGRF was observed, with a mean absolute percent error (MAPE) below 7.3%. Accuracy indices from our equations were better than previously developed equations. All pLR prediction equations presented a lower accuracy compared to those developed to predict pGRF. Walking and running pGRF can be predicted with high accuracy by accelerometry‐based equations, representing an easy way to determine mechanical loading in free‐living conditions. The pLR prediction equations yielded a somewhat lower prediction accuracy compared with the pGRF equations.

Veras, Lucas↗
DOE OSTI2024

An open-source framework for balancing computational speed and fidelity in production cost models

Studies of bulk power system operations need to incorporate uncertainty and sensitivity analyses, especially around exposure to weather and climate variability and extremes, but this remains a computational modeling challenge. Commercial production cost models (PCMs) have shorter runtimes, but also important limitations (opacity, license restrictions) that do not fully support stochastic simulation. Open-source PCMs represent a potential solution. They allow for multiple, simultaneous runs in high-performance computing environments and offer flexibility in model parameterization. Yet, developers must balance computational speed (i.e. runtime) with model fidelity (i.e. accuracy). In this paper, we present Grid Operations (GO), a framework for instantiating open-source, scale-adaptive PCMs. GO allows users to search across parameter spaces to identify model versions that appropriately balance computational speed and fidelity based on experimental needs and resource limits. Results provide generalizable insights on how to navigate the fidelity and computational speed tradeoff through parameter selection. We show that models with coarser network topologies can accurately mimic market operations, sometimes better than higher-resolution models. It is thus possible to conduct large simulation experiments that characterize operational risks related to climate and weather extremes while maintaining sufficient model accuracy.

42 ENGINEERING↗
DOE OSTI2023

The use of synchrotron X-ray fluorescent imaging to study distribution and content of elements in chemically fixed single cells: a case study using mouse pancreatic beta-cells

Synchrotron X-ray fluorescence microscopy (SXRF) presents a valuable opportunity to study the metallome of single cells because it simultaneously provides high-resolution subcellular distribution and quantitative cellular content of multiple elements. Different sample preparation techniques have been used to preserve cells for observations with SXRF, with a goal to maintain fidelity of the cellular metallome. In this case study, mouse pancreatic beta-cells have been preserved with optimized chemical fixation. We show that cell-to-cell variability is normal in the metallome of beta-cells due to heterogeneity and should be considered when interpreting SXRF data. Additionally, we determined the impact of several immunofluorescence (IF) protocols on metal distribution and quantification in chemically fixed beta-cells and found that the metallome of beta-cells was not well preserved for quantitative analysis. However, zinc and iron qualitative analysis could be performed after IF with certain limitations. To help minimize metal loss using samples that require IF, we describe a novel IF protocol that can be used with chemically fixed cells after the completion of SXRF.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Spin-piston problem for a ferromagnetic thin film: Shock waves and solitons

In this report the unsteady, nonlinear magnetization dynamics induced by spin injection in an easy-plane ferromagnetic channel subject to an external magnetic field are studied analytically. Leveraging a dispersive hydrodynamic description, the Landau-Lifshitz equation is recast in terms of hydrodynamic-like variables for the magnetization's perpendicular component (spin density) and azimuthal phase gradient (fluid velocity). Spin injection acts as a moving piston that generates nonlinear, dynamical spin textures in the ferromagnetic channel with downstream quiescent spin density set by the external field. In contrast to the classical problem of a piston accelerating a compressible gas, here, variable spin injection and field lead to a rich variety of nonlinear wave phenomena from oscillatory spin shocks to solitons and rarefaction (expansion) waves. A full classification of solutions is provided using nonlinear wave modulation theory by identifying two key aspects of the fluid-like dynamics: subsonic/supersonic conditions and convex/nonconvex hydrodynamic flux. Familiar waveforms from the classical piston problem such as rarefaction waves and shocks manifest in their spin-based counterparts as smooth and highly oscillatory transitions, respectively, between two distinct magnetic states. The spin shock is an example of a dispersive shock wave, which arises in many physical systems. New features without a gas dynamics counterpart include composite wave complexes with “contact” spin shocks and rarefactions. Magnetic supersonic conditions lead to two pronounced piston edge behaviors including a stationary soliton and an oscillatory wave train. These coherent wave structures have physical implications for the generation of high frequency spin waves from pulsed injection and persistent, stable stationary and/or propagating solitons in the presence of magnetic damping. The analytical results are favorably compared with numerical simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Compensation of heat load deformations using adaptive optics for the ALS upgrade: a wave optics study

In this paper, a realistic wave optics simulation method has been developed to study how wavefront distortions originating from heat load deformations can be corrected using adaptive X-ray optics. Several planned soft X-ray and tender X-ray insertion-device beamlines in the Advanced Light Source upgrade rely on a common design principle. A flat, first mirror intercepts the white beam; vertical focusing is provided by a variable-line-space monochromator; and horizontal focusing comes from a single, pre-figured, adaptive mirror. A variety of scenarios to cope with thermal distortion in the first mirror are studied by finite-element analysis. The degradation of the intensity distribution at the focal plane is analyzed and the adaptive optics that correct it is modeled. The range of correctable wavefront errors across the operating range of the beamlines is reported in terms of mirror curvature and spatial frequencies. The software developed is a one-dimensional wavefront propagation package made available in the OASYS suite, an adaptable, customizable and efficient beamline modeling platform.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

From Reproducible Edge–Cloud Experimentation to Real-World Practice: The E2Clab Experience

Reproducibility is already difficult in distributed systems; on the computing continuum, it becomes substantially harder. Applications that span sensing devices, edge and fog resources, and cloud platforms must be evaluated across heterogeneous hardware, variable network conditions, cross-layer orchestration decisions, and long-running workflow lifecycles. We use E2Clab as a case study to examine these challenges and their implications for experimental methodology. We explain why reproducible experimentation is harder on the continuum, then revisit E2Clab as an initial response based on explicit modeling of infrastructure, workflow lifecycle, and artifacts. Lastly, we discuss how its evolution toward more realistic application settings can be understood through the lens of Translational Computer Science. We argue that reproducible continuum experimentation requires methods that are rigorous enough for research while remaining adaptable to real-world practice.

42 ENGINEERING↗
DOE OSTI2023

A Multi-Fidelity Gaussian Process Regression Method for Probabilistic Wind Farm Power Curve Estimation

Accurate estimation of the power curve for wind turbines or wind farms is crucial to ensure their efficient operation and management. However, conventional methods for power curve estimation rely either on expensive and infrequent measurements or on low-quality numerical simulations. Moreover, the majority of previous studies on power curve estimation for wind turbines or wind farms focused on deterministic estimation, which provides a point estimate of the relationship between wind speed and power generation. Nevertheless, the deterministic approach fails to consider the inherent uncertainty associated with wind energy production resulting from varying turbine characteristics. This can lead to inaccurate power generation estimation and suboptimal decisions regarding energy management. In this paper, a kernel density estimation (KDE) based Multi-Fidelity Gaussian Process Regression (MFGPR) model is proposed to fuse theoretical power curve data and the ground true measurements to create a mapping of wind speed and wind power. By conducting a case study on an actual wind farm in China, the efficacy of the proposed MFGPR model was demonstrated in characterizing the variability of wind power. The probabilistic MFGPR model was also able to generate confidence intervals that encompassed the measured power, thereby improving the accuracy and confidence in wind power estimation or wind resource assessment. Overall, the proposed MFGPR model offers a reliable approach to integrate high-fidelity ground measurements and theoretical power curve data, resulting in precise wind resource assessment and power estimation.

Gaussian process regression↗
DOE OSTI2024

Sea Ice and Cloud Processes Mediating Compensation between Atmospheric and Oceanic Meridional Heat Transports across the CMIP6 Preindustrial Control Experiment

Abstract Bjerknes compensation (BJC) refers to the anticorrelation observed between atmospheric and oceanic heat transport (AHT/OHT) variability, particularly on decadal to longer time scales that may be important to the predictability of the climate system. This study investigates the spread in BJC across fully coupled simulations of phase 6 of the Coupled Model Intercomparison Project (CMIP6) and critical processes (particularly related to sea ice and clouds) that may contribute to that spread. BJC on decadal to longer time scales is confirmed across all the simulations evaluated, and it is strongest in the Northern Hemisphere (NH) between 60° and 70°N. At these latitudes, BJC appears to be primarily driven by the exchange of turbulent fluxes (sensible and latent) in the Greenland, Iceland, and Barents Seas. Metrics to break down how sea ice and clouds uniquely modify the radiative balance of the polar atmosphere during anomalous OHT events are presented. These metrics quantify the impacts of sea ice and clouds on surface and top of atmosphere (latent, sensible, longwave, and shortwave radiative) energy fluxes. Cloud responses tend to counter the clear sky impacts over the Marginal Ice Zone (MIZ). It is further shown that the degree of BJC present in a simulation at high latitudes is heavily influenced by the sensitivity of the sea ice to OHT, which is most influential over the MIZ. These results are qualitatively robust across models and explain the intermodel spread in NH BJC in the preindustrial control experiment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Effects of Horizontal Resolution, Domain Size, Boundary Conditions, and Surface Heterogeneity on Coarse LES of a Convective Boundary Layer

Atmospheric properties in a convective boundary layer vary over a wide range of spatial scales and are commonly studied using large-eddy simulations (LES) in various configurations. We examine how the boundary layer depth and distribution of variability across scales are affected by LES grid spacing, domain size, inhomogeneity of surface properties, and external forcing. Two different setups of the Weather Research and Forecasting (WRF) model are analyzed. A semi-idealized configuration uses a periodic domain, flat surface, prescribed homogeneous surface heat fluxes, and horizontally uniform profiles of large-scale advective tendencies. A nested LES setup employs a larger domain and realistic initial and boundary conditions, including an interactive land surface model with representative topography and vegetation and soil types. Subdomains of identical size are analyzed for all simulations. Characteristic structure sizes are quantified using the variability scales L 50 and L 95 , defined such that features smaller than that contain 50% and 95% of the total variance, respectively. Progressive increase in L 50 from vertical velocity to temperature and moisture structures is systematically reproduced in all simulation configurations. This dependence of L 50 on the considered variable complicates the development of scale-aware parameterizations for models with grid spacing in the “terra incognita”. In simulations using a larger domain with heterogeneous surface properties, the development of internal mesoscale patterns significantly affects variance distributions inside analyzed subdomains. Sizes of boundary layer structures also strongly depend on the LES grid spacing and, in case of heterogeneous surface and topography, on location of the subdomain inside a larger computational domain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Profiling the BLAST bioinformatics application for load balancing on high-performance computing clusters

Abstract Background The Basic Local Alignment Search Tool (BLAST) is a suite of commonly used algorithms for identifying matches between biological sequences. The user supplies a database file and query file of sequences for BLAST to find identical sequences between the two. The typical millions of database and query sequences make BLAST computationally challenging but also well suited for parallelization on high-performance computing clusters. The efficacy of parallelization depends on the data partitioning, where the optimal data partitioning relies on an accurate performance model. In previous studies, a BLAST job was sped up by 27 times by partitioning the database and query among thousands of processor nodes. However, the optimality of the partitioning method was not studied. Unlike BLAST performance models proposed in the literature that usually have problem size and hardware configuration as the only variables, the execution time of a BLAST job is a function of database size, query size, and hardware capability. In this work, the nucleotide BLAST application BLASTN was profiled using three methods: shell-level profiling with the Unix “time” command, code-level profiling with the built-in “profiler” module, and system-level profiling with the Unix “gprof” program. The runtimes were measured for six node types, using six different database files and 15 query files, on a heterogeneous HPC cluster with 500+ nodes. The empirical measurement data were fitted with quadratic functions to develop performance models that were used to guide the data parallelization for BLASTN jobs. Results Profiling results showed that BLASTN contains more than 34,500 different functions, but a single function, RunMTBySplitDB, takes 99.12% of the total runtime. Among its 53 child functions, five core functions were identified to make up 92.12% of the overall BLASTN runtime. Based on the performance models, static load balancing algorithms can be applied to the BLASTN input data to minimize the runtime of the longest job on an HPC cluster. Four test cases being run on homogeneous and heterogeneous clusters were tested. Experiment results showed that the runtime can be reduced by 81% on a homogeneous cluster and by 20% on a heterogeneous cluster by re-distributing the workload. Discussion Optimal data partitioning can improve BLASTN’s overall runtime 5.4-fold in comparison with dividing the database and query into the same number of fragments. The proposed methodology can be used in the other applications in the BLAST+ suite or any other application as long as source code is available.

59 BASIC BIOLOGICAL SCIENCES↗