Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LOADS AND STRESSES, STRUCTURAL”

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 289 records · Page 16

Data-Driven Constitutive Model for the Inelastic Response of Metals: Application to 316H Steel

Here, predictions of the mechanical response of structural elements are conditioned by the accuracy of constitutive models used at the engineering length-scale. In this regard, a prospect of mechanistic crystal-plasticity-based constitutive models is that they could be used for extrapolation beyond regimes in which they are calibrated. However, their use for assessing the performance of a component is computationally onerous. To address this limitation, a new approach is proposed whereby a surrogate constitutive model (SM) of the inelastic response of 316H steel is derived from a mechanistic crystal plasticity-based polycrystal model tracking the evolution of dislocation densities on all slip systems. The latter is used to generate a database of the expected plastic response and dislocation content evolution associated with several instances of creep loading. From the database, a SM is developed. It relies on the use of orthogonal polynomial regression to describe the evolution of the dislocation content. The SM is then validated against predictions of the dead load creep response given by the polycrystal model across a range of temperatures and stresses. When the SM is used to predict the response of 316H during complex non monotonic loading, extrapolating to new loading conditions, it is found that predictions compare particularly well against those from the physics-based polycrystal model.

36 MATERIALS SCIENCE↗

In situ indentation and high cycle tapping deformation responses in a nanolaminate crystalline/amorphous metal composite

The incorporation of nanostructured and amorphous metals into modern applications is reliant on the understanding of deformation and failure modes in constrained conditions. To study this, a 105 nm crystalline Cu/160 nm amorphous Cu 45 Zr 55 (at.%) multilayer structure was fabricated with the two crystalline layers sputter deposited between the top-middle-bottom amorphous layers and prepared to electron transparency. The multilayer was then in situ indented either under a single load to a depth of ~ 100 nm (max load of ~ 100 μN) or held at 20 μN and then repeatedly indented with an additional 5 μN up to 20,000 cycles in a transmission electron microscope to compare the deformation responses in the nanolaminate. For the single indentation test, the multilayer showed serrated load-displacement behavior upon initial indentation inductive of shear banding. Additionally, at an indentation depth of ~ 32 nm, the multilayer exhibited perfect plastic behavior and no strain hardening. Both indented and fatigue-indented films revealed diffraction contrast changes with deformation. Subsequent Automated Crystal Orientation Mapping (ACOM) measurements confirmed and quantified global texture changes in the crystalline layers with specifically identified grains revealing rotation. Using a finite element model, the in-plane displacement vectors under the indent mapped conditions where ACOM determined grain rotation was observed, indicating the stress flow induced grain rotation. Furthermore, the single indented Cu layers also exhibited evidence of deformation induced grain growth, which was not evident in the fatigue-indented Cu based multilayer. Finally, the single indented multilayer retained a significant plastic crater in the upper most amorphous layer that directly contacted the indenter; a negligible crater impression in the same region was observed in the fatigued tested multilayer. These differences are explained by the different loading methods, applied load, and deformation mechanisms experienced in the multilayers.

36 MATERIALS SCIENCE↗

Draft Rules for Alloy 617 Creep-Fatigue Design Using an EPP+SMT Approach

The report describes an improved new creep-fatigue design method for structural components in high temperature nuclear service. The new method uses an integrated elastic-perfectly plastic (EPP) analysis and Simplified Model Test (SMT) approach that greatly simplifies the design process by avoiding the separate evaluation of creep and fatigue damages and eliminating the requirement of stress classification. The report includes draft design rules, presented in a format compatible with an ASME nuclear Code Case, as well as a commentary on the rules and the validation of the rules using pressurized SMT (p-SMT) test data. Areas of potential improvement have been identified for further evaluation. The report also develops EPP+SMT design charts for Alloy 617 at temperatures between 800°C and 950°C by extrapolating high strain range, short hold test data to low strain range and long hold times, typical for structural components in high temperature nuclear service. The design charts are currently limited to 10,000 cycles to avoid excessive extrapolation outside the test database. More robust alternate extrapolation procedures will be considered in future work. The report includes several verification problems comparing the new EPP+SMT design method with current ASME creep-fatigue design methods. The comparison demonstrates that the EPP+SMT method significantly reduces the over-conservatism in current methods. Finally, the report includes a fully-documented sample problem with multiple load cases detailing the application of proposed EPP+SMT design rules.

36 MATERIALS SCIENCE↗

Rapid approach for structural design of the tower and monopile for a series of 25 MW offshore turbines

The goal of further reducing the Levelized Cost of Energy (LCOE) has driven the investigation of large-scale wind turbines. This work presents a simple, rapid and detailed approach for the structural design of the tower and monopile without a controller, but with frequency and high fidelity structural verification. The approach uses an optimization to reduce the mass of the structures while meeting strength, buckling and geometric constraints by using analytical equations. A verification of frequency constraints is performed with BModes, and ANSYS Mechanical APDL is used for high fidelity verification of stress and buckling. The approach is applied to study the design space of three 25 MW offshore wind turbines with different rotor diameters and cone angles, and to evaluate the nacelle center of mass fore-aft location effect. Results obtained show that the tower and monopile are more susceptible to changes in the rotor thrust than the overturning moment even for designs with high pre-cone angle and large distance of the nacelle center of mass from the tower axis. But it is possible to obtain structurally feasible tower and monopile designs for the three 25 MW turbines studied while not exceeding diameter and wall thickness limits. However, mass penalties can be decreased by 0.8-14%, to further reduce the cost of energy, by increasing the diameter limit which may require manufacturing technology development. The approach applied and studies serve to understand the design space of the tower and monopile for a 25 MW turbine, and provide baseline designs that can be used in the development of a controller and evaluation of a full suite of design load cases.

17 WIND ENERGY↗

Status Report on Design of In-situ Thermomechanical Testing at LANSCE

Nuclear fuel encounters severe thermomechanical environments in which its mechanical response is determined by its microstructure, temperature and stress level histories. Simulating the response of such microstructures is crucial for predicting both performance and transient fuel mechanical responses and experimental verification of such predictions is therefore of great interest. While most of the deformation in a nuclear fuel rod occurs in the cladding, deformation of the fuel itself is still of interest with deformation mechanisms at operating temperature and above including creep, swelling, cracking as well as pellet-clad interaction. Characterization of these properties and understanding of the underlying deformation phenomena at operating or excursion temperatures is therefore of great importance for development and ultimately licensing of improved and novel nuclear fuel forms. Diffraction techniques offer unique insight on the atomistic (e.g. crystal structure) and microstructure (e.g. phase transformations, texture, defects) length scales and have a long history of providing unique data to inform relevant deformation models that enable the required predictive capabilities. For example, dislocations lead to diffraction peak broadening that can be characterized to estimate the dislocation density and study the role of dislocations on the deformation while measuring lattice strains allows to studie load sharing in two phase materials. In this report the requirements for a sample environment for high temperature deformation of nuclear fuels are defined. The HIPPO neutron time-of-flight diffractometer at LANSCE will host this sample environment and is also described. This instrument covers diffraction angles from 140° to 40° and is also equipped with an event-mode neutron imaging detector system, enabling energy-resolved neutron imaging in parallel with the diffraction that could measure sample temperature from Doppler broadening of neutron absorption resonances or measure pore densities from changes in the attenuation. Designs of devices to characterize thermomechanical properties of nuclear fuel without diffraction are also considered to guide the design. While this report is focused on applications for nuclear fuels, the device can also characterize cladding, moderator or structural materials and therefore contribute to other fields of research and development for advanced reactors. The temperatures planned to be reached are above 2000℃, thus enabling characterization of LWR reactor fuels under accident scenarios but also reaching temperatures of fuels developed for nuclear thermal propulsion and providing opportunities to characterize those. In conjunction with the energy-resolved neutron imaging detector, this setup would allow to measure neutron cross-sections at high temperatures, filling a gap towards development of reactors operating at high temperatures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling of Stress and Temperature Effects on Creep of Reduced Activation Ferritic-Martensitic Steel Alloy F82H (Tertiary Creep Modeling of RAFM Steel)

A Bayesian optimization procedure is presented for calibrating a multi-mechanism micromechanical model for creep to experimental data of F82H steel. Reduced activation ferritic martensitic (RAFM) steels based on are the most promising candidates for some fusion reactor structures. Although there are indications that RAFM steel could be viable for fusion applications at temperatures up to 600 °C, the maximum operating temperature will be determined by the creep properties of the structural material and the breeder material compatibility with the structural material. Due to the relative paucity of available creep data on F82H steel compared to other alloys such as Grade 91 steel, micromechanical models are sought for simulating creep based on relevant deformation mechanisms. As a point of departure, this work recalibrates a model form that was previously proposed for Grade 91 steel to match creep curves for F82H steel. Due to the large number of parameters (9) and cost of the nonlinear simulations, an automated approach for tuning the parameters is pursued using a recently developed Bayesian optimization for functional output (BOFO) framework [1]. Incorporating extensions such as batch sequencing and weighted experimental load cases into BOFO, a reasonably small error between experimental and simulated creep curves at two load levels is achieved in a reasonable number of iterations. Validation with an additional creep curve provides confidence in the fitted parameters obtained from the automated calibration procedure to describe the creep behavior of F82H steel at 600 °C. The model is further extended using a temperature dependent scaling law approach to simulate creep response between 550 °C and 650 °C. The efficacy of this extension is compared with the previously used scaling law approach for Grade 91 steel.

36 MATERIALS SCIENCE↗

Dynamic Bottlebrush Polymer Networks: Self-Healing in Super-Soft Materials

Here, we introduce a design strategy to expand the range of accessible mechanical properties in covalent adaptable networks (CANs) using bottlebrush polymer building blocks. Well-defined bottlebrush polymers with rubbery poly(4-methylcaprolactone) side-chains were crosslinked in formulations that include a bislactone and strong Lewis acid (tin ethylhexanoate). The resulting materials exhibit tunable stress relaxation rates at elevated temperatures (160–180 °C) due to dynamic ester crosslinks that undergo transesterification with residual hydroxy groups. Varying the crosslinker loading or bottlebrush backbone degree of polymerization yields predictable low-frequency shear moduli ca. 10–100 kPa, well below values typical of linear polymer CANs (1 MPa). These extensible networks can be stretched to strains as large as 350% before failure and undergo efficient self-healing to recover >85% of their original toughness upon repeated fracture and melt processing. In summary, molecular architecture creates new opportunities to tailor the mechanical properties of CANs in ways that are otherwise difficult to achieve.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Air Blast Mesh Sensitivity and Pressure Mapping Study

Nose cone structural and aerodynamic characteristics are essential for intelligent design of aircraft, spacecraft, and ballistic systems. Finite element analysis can be used to help understand the structural integrity and flight characteristics of different nose cones. A mesh sensitivity study was undertaken for a particular nose cone geometry that was used in tests at LANL facilities in order to confirm the integrity of the meshed geometry. A simple cone that best matched closed-form theoretical solutions was modeled, and received good correlation to the theory. Complexity was then added back to the nose cone. Parameters applied to the simple cone were then implemented in the nose cone geometry giving assurance of accuracy after the geometry was changed. Nose cone results averaged 6.3% error for radial displacement when compared with the theoretical. Hoop stress averaged 6.0% error and meridional stress averaged 5.7% error at the finest mesh level. Meshes showed signs of convergence when compared to all three theoretical solutions. Finally, pressure time-history data from LANL computational fluid dynamics simulations was applied to the surface of the final nose cone geometry. The pressure data was interpolated from pressure gauge locations onto nearby meshed elements, which allowed for FEA software to run simulations on the cone with the pressure data as a loading condition. The pressure mapping resulted in the ability to understand the nose cone’s rigid body motion that in turn can inform design of future nose cones.

42 ENGINEERING↗

Variationally consistent Maxwell stress in flexoelectric structures under finite deformation and immersed in free space

Maxwell stress refers to the mechanical stress exerted on a dielectric material due to the presence of electric fields. It plays a significant role in the interaction between a dielectric material and the surrounding free space under finite deformation. Previous research on finite deformation of flexoelectricity mainly adopted a modified form of Maxwell stress, potentially not able to correctly capture some physical phenomena, such as the compression of a dielectric droplet in an electric field. In this work, we propose a consistent and complete variational principle for flexoelectricity, in which the Maxwell stress emerges naturally from the derivation, without introducing additional assumptions. An Isogeometric analysis-based numerical framework is developed accordingly and verified by both linear and nonlinear benchmark cases compared with experimental results. The present framework successfully captures and quantifies the behaviors of conductive liquids and soft dielectric solids subjected to an external electric field. Finally, a novel scenario is investigated in which a flexoelectric beam immersed in free space is analyzed, showing the interesting distribution of Maxwell stress-induced tractions at opposing boundaries. The test demonstrates that a higher dielectric constant can effectively enhance the material's stiffness in response to the external electric loading.

36 MATERIALS SCIENCE↗

Evaluation of Graph Analytics Frameworks Using the GAP Benchmark Suite

The analysis of connected data is an increasingly important application in high-performance computing. Such analyses can reveal fraudulent patterns in financial transactions, optimize telecommunications networks, predict information flow in social networks, etc. However, the landscape of graph analytics is highly diverse. Graph algorithms stress processor architectures differently, and no one graph can represent all topologies. Consequently, no single approach or framework is expected to be optimal for all graph analytics problems. To help make sense of this diverse landscape, we evaluated four approaches to graph analytics: GraphBLAS, Galois, BGL17, GraphIt; and compare them against hand-tuned implementations that take advantage of hardware features on our test platform. Graph- BLAS formulates graph analytics as sparse linear algebra. Galois provides syntactic constructs for data parallelism over irregular data structures. BGL17 is a generic C++ template library for implementing graph algorithms. GraphIt provides a domain- specific language to describe and optimize graph algorithms. We use the GAP Benchmark Suite to establish baseline performance and guide the side-by-side evaluation of each framework. GAP consists of 30 tests: six graph analytics algorithms (breadth- first search, single-source shortest path, PageRank, betweenness centrality, connected components, and triangle counting) run on five graphs, each with different topological characteristics (e.g., high diameter, skewed degree distribution, high average degree). High-performance reference implementations are included for each benchmark algorithm. Because a graph can be loaded into memory a number of ways (e.g., flat file on disk, compressed sparse format, data frames, retrieved from SQL or NoSQL databases), our evaluation focused on computational performance rather than I/O. Our results show the relative strengths of each framework.

Graph algorithms, Benchmarking, shared-memory prog↗

Rheology of concentrated fiber suspensions with a load-dependent friction coefficient

Here we investigate the effects of fiber aspect ratio, roughness, flexibility, and volume fraction on the rheology of concentrated suspensions in a steady shear flow using direct numerical simulations. We model the fibers as inextensible continuous flexible slender bodies with the Euler-Bernoulli beam equation governing their dynamics suspended in an incompressible Newtonian fluid. The fiber dynamics and fluid flow coupling is achieved using the immersed boundary method. In addition, the fiber surface roughness might lead to interfiber contacts, resulting in normal and tangential forces between the fibers, which follow Coulomb’s law of friction. The surface roughness is modeled as hemispherical protrusions on the fiber surfaces. Their deformation results in a normal load-dependent friction coefficient. Our simulations accurately predict the experimentally observed shear thinning in fiber suspensions. Furthermore, we find that the suspension viscosity η increases with increasing the volume fraction, roughness, fiber rigidity, and aspect ratio. The increase in η is the macroscopic manifestation of a similar increase in the microscopic contact contribution to the total stress with these parameters. In addition, we observe positive and negative first N 1 and second N 2 normal stress differences, respectively, with |N 2 | < |N 1 |, in agreement with previous experiments. Last, we propose a modified Maron-Pierce law to quantify the reduction in the jamming volume fraction by increasing the fiber aspect ratio and roughness. Our results and analysis establish the use of fiber surface tribology to tune the suspension flow behavior.

42 ENGINEERING↗

A method for assessing economic, environmental, and reliability tradeoffs of interregional transmission connecting ERCOT (the Texas grid) to the eastern and western grids

Reliable development of the power grid is an evolving concern for humanity due to extreme weather that frequently threatens power sector infrastructure. The state of Texas is a uniquely structured testbed for grid planners to study when looking for solutions to development, innovation, and overcoming such challenges. Because of its size and islanded structure, Texas is small enough to model, but big enough to matter. Texas is a global leader in energy production, energy consumption, and maintains an unusually diverse fuel mix. In addition, the state has experienced winter freezes, heat waves, wind storms, droughts and floods that have threatened power sector infrastructure or caused recent blackouts and calls for demand side conservation. One of the most devastating of these events was the North American winter storm, dubbed “Winter Storm Uri” by the Weather Channel, that froze the region in February 2021 and led to an extended power outage event that put the majority of Texan residents in darkness for days. While preparing to avoid such outage events in the future, various tools have been proposed to improve grid reliability, including energy efficiency, demand response, and distributed energy resources. An additional option would be to develop interregional transmission that connects the Texas grid to other national grids. To assess the merits of this idea, we developed a novel, universally-applicable and internationally-relevant framework to study how the Texas grid would evolve alongside access to various interregional ties. This method allows us to stress the synthetic grid structure and analyze how it would respond to the shock of a simulated winter storm event. Our method leverages open-source modeling tools, such as PowerGenome, pyGRETA, and GenX to synthesize unique zonal grid data, construct a consolidated network of model regions, and simulate different developmental pathways of capacity expansion and operational dispatch. We demonstrate our method with an analysis connecting the Electric Reliability Council of Texas (ERCOT), the grid that serves most of Texas, the Western Electricity Coordinating Council (WECC), the grid that serves the western half of the contiguous U.S., and the Eastern Interconnect, the grid that serves the eastern half of the contiguous U.S. Our results indicate that the cost-optimal capacity of interregional transmission connecting the ERCOT grid to other grids lies between 9–13 GW assuming baseline conditions. Building this amount of connecting capacity in one or multiple directions lowers the costs and emissions of development and operation by up to $16 billion and 257 million metric tonnes (MMT) respectively. Additionally, our results show that the interregional connections between ERCOT and other national grids reduce the amount of total load shed required through mild winter storm events. However, our results also show that there is a threshold of very extreme winter storm conditions, spanning multiple service areas, above which the connections exacerbate resource adequacy problems. Therefore, the results indicate that the connections need to be carefully planned alongside the rest of the grid infrastructure to avoid over-reliance on specific resources or technology options.

24 POWER TRANSMISSION AND DISTRIBUTION↗

MoorDyn V2: New Capabilities in Mooring System Components and Load Cases

MoorDyn, an open-source mooring dynamics model, is being expanded with capabilities for additional mooring system features and load cases. As floating wind turbine technology matures, mooring systems are becoming more sophisticated and more complex scenarios need to be considered in the design process. Mooring systems may have synthetic line materials, ballast/buoyancy bodies along the lines, or interconnections between platforms. Failure modes may involve multiple cascading line failures that depend on mooring system dynamics. Features recently added to MoorDyn aim to address these emerging needs. MoorDyn’s linear elasticity model has been supplemented to support user-defined stress-strain curves, which can be adjusted to represent synthetic mooring materials. Rigid six-degree-of-freedom bodies in the mooring system can now be modeled using two new model objects. “Rod” objects provide an option for rigid cylindrical bodies. They use the existing Morison equation-based hydrodynamics model and can be connected to mooring lines at either end. “Body” objects provide a generic six-degree-of-freedom rigid-body representation based on a lumped-parameter model of translational and rotational properties. Rod objects can be added to Body objects and mooring lines can be attached at any location, allowing a wide variety of submerged structures to be integrated into the mooring system. Lastly, a means of dynamically simulating mooring line failures has been implemented. These new features, currently in the C++ version of MoorDyn, are described and then demonstrated on a two-turbine shared-mooring array. A qualitative view of the results suggests the new features are functioning as expected.

49 EE - Wind and Water Power Program - Wind (EE-4W↗

MoorDyn V2: New Capabilities in Mooring System Components and Load Cases: Preprint

MoorDyn, an open-source mooring dynamics model, is being expanded with capabilities for additional mooring system features and load cases. As floating wind turbine technology matures, mooring systems are becoming more sophisticated and more complex scenarios need to be considered in the design process. Mooring systems may have synthetic line materials, ballast/buoyancy bodies along the lines, or interconnections between platforms. Failure modes may involve multiple cascading line failures that depend on mooring system dynamics. Features recently added to MoorDyn aim to address these emerging needs. MoorDyn’s linear elasticity model has been supplemented to support user-defined stress-strain curves, which can be adjusted to represent synthetic mooring materials. Rigid six-degree-of-freedom bodies in the mooring system can now be modeled using two new model objects. “Rod” objects provide an option for rigid cylindrical bodies. They use the existing Morison equation-based hydrodynamics model and can be connected to mooring lines at either end. “Body” objects provide a generic six-degree-of-freedom rigid-body representation based on a lumped-parameter model of translational and rotational properties. Rod objects can be added to Body objects and mooring lines can be attached at any location, allowing a wide variety of submerged structures to be integrated into the mooring system. Lastly, a means of dynamically simulating mooring line failures has been implemented. These new features, currently in the C++ version of MoorDyn, are described and then demonstrated on a two-turbine shared-mooring array. A qualitative view of the results suggests the new features are functioning as expected.

49 EE - Wind and Water Power Program - Wind (EE-4W↗

Analysis of thermal and mechanical properties with inventory level of the molten salt storage tank in central receiver concentrating solar power plants

Molten salt thermal energy storage (TES) tanks ensure steady power output of concentrating solar power (CSP) plants; however, recent tank failures have highlighted the need for further analysis. Current studies primarily focus on analyzing the molten salt flow, heat transfer, and thermal efficiency. Additionally, research on the latest tank structures is limited and lacks newest experimental validation. This study measures temperature and molten salt inventory levels in the high-temperature tank at a 50 MW central receiver CSP plant, connected to the power grid in 2019. A multi-physics model was developed to evaluate thermal and mechanical properties of TES tanks by combining computational fluid dynamics and finite element modeling using real plant data. Heat loss, temperature, displacement, and stress distribution of the tank at different inventory levels were investigated. Results show that ambient air velocity near the tank roof reaches 2.14 m/s, much higher than 0.2 m/s near the wall. The temperatures of inventory fluid and tank are close, varying slightly at different levels due to thermal conduction and radiation. Because the heat loss strongly depends on temperature, the total tank loss remains nearly constant across inventory levels. Larger temperature gradients and thermal stresses are primarily localized along the tank floor edge and the air-salt interface. Notably, the maximum thermal stress at the tank edge is three times higher than that at the interface. The magnitude of total stress changes by less than 5 MPa with and without thermal load, indicating that high temperatures exert only a minor impact on tank stress. In contrast, thermal load significantly affects tank deformation, particularly at the roof edge, where values exceed 150 mm. Despite the large variation in molten salt levels, tank wall temperatures and displacements present a minor change, suggesting a weak correlation with inventory levels. In conclusion, the findings obtained in this study provide important insights on the TES tank that could be used to optimize tank design and operation strategies.

14 SOLAR ENERGY↗

A conformal heat-drying direct ink writing 3D printing for high-performance lithium-ion batteries

High areal capacity electrodes hold great potential for high-energy density lithium-ion batteries (LIBs), but their poor electrochemical kinetics limit their power density. Here, in this study, high areal capacity 3D-structured LiNi 0.8 Mn 0.1 Co 0.1 O 2 cathodes (4.3 mAh cm -2 ) are prepared via 3D printing with a manner of direct ink writing. The electrodes had an enlarged electrode–electrolyte contact area, shortened diffusion pathway, and reduced intercalation-induced stress, thereby delivering enhanced rate capability and cyclability in LIBs, which is 143.6 mAh g -1 at 3C and a 60.2 % capacity retention over 800 cycles at 1C. Moreover, at electrode level, the 3D-NMC exhibits an energy and power densities of 313.1 Wh kg -1 and 657.9 W kg -1 , respectively. Furthermore, the theoretical calculation suggests that reducing the gap width will be highly beneficial to the energy and power densities. This work establishes a milestone in understanding the cycling effect on the electrode local structure, including the void area and the LiNi 0.8 Mn 0.1 Co 0.1 O 2 region, which confirms the effectiveness of 3D printing for electrode preparation.

25 ENERGY STORAGE↗

COMPILE: a GWAS computational pipeline for gene discovery in complex genomes

Abstract Background Genome-Wide Association Studies (GWAS) are used to identify genes and alleles that contribute to quantitative traits in large and genetically diverse populations. However, traits with complex genetic architectures create an enormous computational load for discovery of candidate genes with acceptable statistical certainty. We developed a streamlined computational pipeline for GWAS (COMPILE) to accelerate identification and annotation of candidate maize genes associated with a quantitative trait, and then matches maize genes to their closest rice and Arabidopsis homologs by sequence similarity. Results COMPILE executed GWAS using a Mixed Linear Model that incorporated, without compression, recent advancements in population structure control, then linked significant Quantitative Trait Loci (QTL) to candidate genes and RNA regulatory elements contained in any genome. COMPILE was validated using published data to identify QTL associated with the traits of α-tocopherol biosynthesis and flowering time, and identified published candidate genes as well as additional genes and non-coding RNAs. We then applied COMPILE to 274 genotypes of the maize Goodman Association Panel to identify candidate loci contributing to resistance of maize stems to penetration by larvae of the European Corn Borer ( Ostrinia nubilalis ). Candidate genes included those that encode a gene of unknown function, WRKY and MYB-like transcriptional factors, receptor-kinase signaling, riboflavin synthesis, nucleotide-sugar interconversion, and prolyl hydroxylation. Expression of the gene of unknown function has been associated with pathogen stress in maize and in rice homologs closest in sequence identity. Conclusions The relative speed of data analysis using COMPILE allowed comparison of population size and compression. Limitations in population size and diversity are major constraints for a trait and are not overcome by increasing marker density. COMPILE is customizable and is readily adaptable for application to species with robust genomic and proteome databases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗