Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical integration”

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 451 records · Page 25

Sorbmatech{sup R}Cs, an Innovative Highly Selective Cs Sorbent for Wastewater Decontamination - 20065

The treatment of radioactive wastewater from the operating nuclear industry (energy production, reprocessing of spent fuel, decommissioning of shut-down facilities, etc.) or in case of incidental or accidental situation (such as Fukushima Daichii NPP accident) is an important issue requiring constant improvement of the technologies and techniques used to prevent radioactive releases into the environment and/or to reduce the production of secondary waste. CEA, Orano and CTI are collaborating since 2014, with the support of the French government's program 'investments for the future', to develop Sorbmatech{sup R}Cs, an innovative and patented selective sorbent (Grandjean and Barre 2014) to deal with the removal of radioactive {sup 137}Cs, one of the most hazardous isotope to be managed in nuclear waste, highly efficient either in agitated reactor or in column process. Sorbmatech{sup R}Cs sorbent consists of mesoporous silica particles loaded with potassium-copper hexacyanoferrate nanoparticles (KCu-HCF). This peculiar structure enables much faster sorption kinetics compared to denser crystalline structures, such as Silicotitanate or Zeolite sorbents, whose drawback consists in hindering the diffusion process through the material, and whose larger sorbent particles slow down the ion exchange process. The equilibrium of Sorbmatech{sup R}Cs is reached in few minutes, even at high flow rates in a fixed-bed column process (up to 10 m.h{sup -1} Darcy velocity) with low drop pressure. Moreover, Sorbmatech{sup R}Cs can be used over a wide pH range from 2 to 10 and has a very high selectivity towards {sup 137}Cs even in brine solution (saline, mineral or sea water). This new sorbent also has a K{sub d} higher than 10{sup 5} mL.g{sup -1} depending on the salinity of the wastewater, enabling a highly efficient removal of Cs trace in radioactive wastewater compared to Cs selective sorbent such as Zeolites and Silicotitanates. Its radiological resistance enables the treatment of high, medium and low activity wastewater. The CEA-Orano-CTI/ALSYS partnership leads to the development of a simple and low-cost synthetic route, which demonstrated large-scale manufacturing of the product on a batch of 100 kg. Two products have been developed and are ready for short run market introduction so far with grain size well adapted either for column processes or agitated reactors.. Finally, the industrial implementation has been studied and largely benefits from the properties of the sorbent. Indeed, its high sorption kinetics enable high flow rates resulting in a decrease of operation times, unit sizes, and volumes of secondary waste produced. Moreover, Sorbmatech{sup R}Cs is designed to be combined into a full and integrated solution, ranging from preliminary optimization of operating conditions through numerical dimensioning model, industrial implementation and treatment of the wastewater (through cartridge or preconditioned columns systems developed thanks to Orano's extensive expertise), to long-term stabilization of generated waste (under progress). (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Thermal Management System of an Outer-Rotor-Motor-Based Traction Drive With Integrated Power Electronics in its Central Cavity

Increasing power densities of electric vehicle traction drive systems necessitates combining the electric motor and the power electronics into one unit. A compact, integrated traction drive unit with fewer components also drives production costs down, enabling wider adoption of electric vehicles. However, the integration of power electronics in the electric machine is associated with challenges of designing an effective thermal management solution for the combined traction drive system. This paper focuses on the thermal management approach selected for Oak Ridge National Laboratory’s outer-rotor-motor-based integrated traction drive and evaluates its potential performance. The outer-rotor-motor configuration provided an opportunity for integration of the six-phase inverter in the available space in the central cavity of the internal stator. A cylindrical inverter enclosure with integrated coolant (water-ethylene glycol) channels in its walls was designed to enable heat removal from the power electronics. Numerical thermal-fluid modeling and initial channel/fin optimization results for the cylindrical heat sink are presented here. As permanent magnets are integrated in a high-speed (20 000 RPM) outer rotor, forced air convection provides cooling for the magnets. The magnets were segmented axially to mitigate the eddy current losses. Heat generated in the stator windings and laminations is removed by water-ethylene glycol coolant circulating in interpolar T-shape ceramic heat exchangers inserted between windings. General design concepts and numerical thermal-fluid simulations illustrating the electric motor thermal management solution are also presented.

Kekelia, Bidzina↗

Robust Medium-Voltage Distribution System State Estimation using Multi-Source Data

Due to the lack of sufficient online measurements for distribution system observability, pseudo-measurements from short-term load or distributed renewable energy resources (DERs) forecasting are used. However, the accuracy of them is low and thus significantly limits the performance of distribution system state estimation (DSSE). In this paper, a robust DSSE that integrates multi-source measurement data is proposed. Specifically, the historical low-voltage (LV) side smart meters are used to forecast load and DERs injections via the support vector machine (SVM) with optimally tuned parameters. By contrast, the online smart meters at LV side are utilized to derive equivalent power injections at the MV/LV transformers, yielding more accurate pseudo-measurements compared to the forecasted injections. Furthermore, to deal with bad data caused by communication loss, instrumental errors and cyber attacks, robust DSSE that relies on generalized maximum-likelihood (GM)-estimation criterion is developed. The projection statistics are developed to adjust the weights of each measurement, leading to better balance between pseudo- and real-time measurements. Numerical results conducted on modified IEEE 33-bus system with DG integration demonstrate the effectiveness and robustness of the proposed method.

distribution system state estimation↗

Fully implicit crystal plasticity models representing orientations with modified Rodrigues parameters

Here, this work describes a crystal plasticity formulation combining several mathematical, numerical, and implementation choices to produce a highly efficient model. Specifically, the key choices in the implementation are (1) representing orientations with modified Rodrigues parameters, (2) implementing a fully coupled implicit time integration for the elastic stretch, the crystal orientations, and the model internal variables, (3) implementing the model in the NEML2 constitutive modeling framework, based on PyTorch, to vectorize the calculations and port the computation to GPUs and other hardware accelerators, and (4) an exact implementation of the consistent tangent matrix, even for arbitrary coupling to other field variables beyond the displacements, like temperature, neutron fluence, etc. The first two features of the model are, to our knowledge, novel. The paper considers each of these choices individually as well as the final model as a whole. This includes a full description of modified Rodrigues parameters, their advantages over other representations of orientations, the mathematical formulae and tools required to implement a model with modified Rodrigues parameters, and a detailed description of the geometry of the space of modified Rodrigues parameters (in an appendix). It also includes a description of a fully implicit time integration scheme for the orientations and the advantages in representing orientations with modified Rodrigues parameters in implementing such a model. The work then assess, via numerical examples, the advantages of fully coupled implicit time integration versus more common decoupled and explicit time integration schemes. These studies demonstrate the computational advantages of fully coupled integration versus other time integration algorithms, though the performance of the competing models depends on the complexity of the underlying single crystal model. The study concludes by demonstrating that the choice of time integration method affects the sharpness of the predicted texture, with explicit methods for integrating the orientations overestimating texture sharpness and implicit methods underestimating texture sharpness.

Crystal plasticity↗

Large-Eddy Simulation Study of Flow and Combustion Dynamics in a Full-Scale Hydrogen–Air Rotating Detonation Combustor-Stator Integrated System

In the present work, a first-of-its-kind three-dimensional (3D) large-eddy simulation (LES) study is conducted to numerically investigate the combustion dynamics as well as aero-thermal phenomena in a full-scale nonpremixed hydrogen–air rotating detonation engine (RDE) (with a diverging-shaped lower-end wall), when integrated with nozzle guide vanes (NGV) acting as the turbine stator. The wall-modeled LES framework incorporates hydrogen–air detailed chemical kinetics and adaptive mesh refinement (AMR). A comparative analysis is carried out for two operating conditions with different fuel/air mass flow rates but global equivalence ratio of unity, and considering RDE configurations without and with stator. The LES model is validated against available experimental data for the low mass flux condition with respect to detonation wave speed/height, wave dynamics, and axial static pressure distribution. Numerical results indicate significant deflagrative combustion occurring in the fill region near the inner wall due to formation of recirculation zones in the injection near-field driven by the backward facing step. The leading detonation wave is found to be trailed by an azimuthal reflected-shock combustion (ARSC) wave, consistent with experimental observations, which consumes unburned vitiated reactants that leak through the main detonation wave. The main detonation wave characteristics, such as detonation wave speed/height and combustion efficiency, do not change appreciably with the presence of NGV. A novel combustion diagnostic technique based on chemical explosive mode analysis (CEMA) is employed to quantify the fraction of heat release occurring in the detonative mode versus deflagrative mode for the simulated conditions. The exit flow is found to be nearly fully subsonic and supersonic for the low and high mass flux conditions, respectively. Further analysis of the exit flow profiles shows that the presence of NGV renders the flow more axial and significantly impacts the exit Mach number and total pressure, while the total temperature shows negligible change. In addition, the low mass flux operating point, despite exhibiting more deflagrative losses within the combustor, yields overall lower pressure drop from plenum to exhaust, which is mainly attributed to lower pressure drop across the injectors. Lastly, the rotating detonation engine-nozzle guide vanes (RDE-NGV) configuration exhibits higher total pressure loss compared to rotating detonation engine (RDE) without stator across both the mass flux conditions. In conclusion, this study extends the state-of-the-art in numerical modeling of pressure gain combustion (PGC) systems by demonstrating high-fidelity 3D reacting LES of full-scale RDE-NGV systems relevant to RDE-turbine integration for stationary power generation.

Combustion↗

Computation of the Biot–Savart line integral with higher-order convergence using straight segments

One common approach to computing the magnetic field produced by a filamentary current-carrying coil is to approximate the coil as a series of straight segments. The Biot–Savart field from each straight segment is analytically known. However, if the endpoints of the straight segments are chosen to lie on the coil, then the accuracy of the Biot–Savart computation is generally only the second order in the number of endpoints. In this work, we propose a simple modification: shift each end point of the coil in the outward normal direction by an amount proportional to the local curvature. With this modification, the Biot–Savart accuracy increases to the fourth order and the numerical error is dramatically reduced for a given number of discretization points.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improving Solution Accuracy and Convergence for Stochastic Physics Parameterizations with Colored Noise

Stochastic parameterizations are used in numerical weather prediction and climate modeling to help capture the uncertainty in the simulations and improve their statistical properties. Convergence issues can arise when time integration methods originally developed for deterministic differential equations are applied naively to stochastic problems. In previous studies, it has been demonstrated that a correction term, known in stochastic analysis as the Itô correction, can help improve solution accuracy for various deterministic numerical schemes and ensure convergence to the physically relevant solution without substantial computational overhead. The usual formulation of the Itô correction is valid only when the stochasticity is represented by white noise. In this study, a generalized formulation of the Itô correction is derived for noises of any color. The formulation is applied to a test problem described by an advection–diffusion equation forced with a spectrum of fast processes. We present numerical results for cases with both constant and spatially varying advection velocities to show that, for the same time step sizes, the introduction of the generalized Itô correction helps to substantially reduce time integration error and significantly improve the convergence rate of the numerical solutions when the forcing term in the governing equation is rough (fast varying); alternatively, for the same target accuracy, the generalized Itô correction allows for the use of significantly longer time steps and, hence, helps to reduce the computational cost of the numerical simulation.

54 ENVIRONMENTAL SCIENCES↗

Performance Investigation of High-temperature Latent Heat Storage for Integration with Advanced Nuclear Reactors

Integrating thermal energy storage (TES) with advanced nuclear reactors enhances the flexible use of nuclear energy, facilitating its expansion beyond conventional electricity generation. High-temperature TES, capable of storing heat above 400 degree C, emerges as a vital carbon-free energy solution for decarbonizing industrial sector, especially when combined with low-emission energy sources like advanced nuclear reactors. Idaho National Laboratory (INL) has recently developed a novel design for high-temperature latent heat storage system, called Heat pipe Integrated Thermal Battery (HITB), and efforts are underway to experimentally demonstrate the concept. HITB employs liquid-metal heat pipes to establish thermal linkage between TES and nuclear systems without direct fluid exchange, minimizing the potential risk of the integrated systems such as cross-contamination. HITB is designed to achieve high charging and discharging efficiency as well as high energy storage density by employing metal alloys as Phase Change Material (PCM). Initial proof-of-concept experiments are being undertaken using an aluminum alloy (Al59%-Mg35%-Zn6%). This paper discusses current progress of the HITB project at INL, seeking to develop high-temperature TES for versatile integration with advanced nuclear reactors, and shares insights from the performance evaluations via numerical modeling and analysis.

25 ENERGY STORAGE↗

On numerical errors to the fields surrounding a relativistically moving particle in PIC codes

The particle-in-cell (PIC) method is widely used to model the self-consistent interaction between discrete particles and electromagnetic fields. It has been successfully applied to problems across plasma physics including plasma based acceleration, inertial confinement fusion, magnetically confined fusion, space physics, astrophysics, high energy density plasmas. In many cases the physics involves how relativistic particles (those with high relativistic γ factors) are generated and interact with plasmas. However, when relativistic particles stream across the grid, both in “vacuum” and in plasma, many numerical issues may arise which can lead to unphysical results. We present a detailed analysis of how discretized Maxwell solvers used in PIC codes can lead to numerical errors to the fields that surround particles that move at relativistic speeds across the grid. Expressions for the axial electric field as integrals in k space are presented that reveal two types of errors. The first arises from errors to the numerator of the integrand and leads to unphysical fields that are antisymmetric about the particle. Furthermore, the second arises from errors to the denominator of the integrand and lead to Cherenkov like radiation in “vacuum”. These fields are not anti-symmetric, extend behind the particle, and cause the particle to accelerate or decelerate depending on the solver and parameters. The unphysical fields are studied in detail for two representative solvers - the Yee solver and the FFT based solver. Although the Cherenkov fields are absent, the space charge fields are still present in the fundamental Brillouin zone for the FFT based solvers. In addition, the Cherenkov fields are present in higher order zones for the FFT based solvers. Comparison between the analytical solutions and PIC simulation results are presented. A solution for eliminating these unphysical fields by modifying the k operator in the axial direction is also presented. Using a customized finite difference solver, this solution was successfully implemented into OSIRIS [1]. Results from the customized solver are also presented. Additionally, this solution will be useful for a beam of particles that all move in one direction with a small angular divergence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Tunneling electroresistance effects in epitaxial complex oxides on silicon

Complex oxide-based ferroelectric tunnel junctions (FTJs) show excellent nonvolatile memory characteristics promising for emerging technology. However, integration of these epitaxially grown FTJs electrically with a silicon substrate remains challenging due to their incompatible lattice structures and poor electronic interfaces resulting from the direct synthesis techniques. Here, we present an epitaxial SrRuO 3 /PbZr 0.2 Ti 0.8 O 3 /SrRuO 3 FTJ integrated electrically with a doped silicon substrate after a layer transfer process. The tunnel currents of the FTJ on silicon show a large tunneling electroresistance (~1×10 5 %) effect, which is explained by a numerical FTJ model incorporating pinned dipoles at the interfaces. This proof of concept of the integration of functional oxide heterostructures with silicon opens a pathway to beyond-CMOS computing devices using unconventional materials.

36 MATERIALS SCIENCE↗

Computing neutrino cross sections from Euclidian responses

Energy integrated neutrino cross sections are integrals of nuclear responses weighted with kinematic prefactors. We decompose the prefactors into a limited set of functions of energy transfer and show the relevant integrals are the moments of the responses, and integrals weighted with $1/(a+ω)^n$ with $n\leq 2$. These can be directly obtained from the Euclidean response, avoiding the need for inversion of the Laplace transform. As a proof of concept we study the procedure with toy-model responses for the quasielastic peak. We show that the different contributions can be straightforwardly organized in terms of relative importance, and how flux-averaged cross sections can be obtained. Using a realistic model for the response and numerical uncertainty we show that it is feasible to obtain the required integrals from the Euclidean response, with large uncertainties only for the third moment. Due to kinematic restrictions, the integrals contain contributions from the unphysical region for neutrino scattering, coming from high-momentum nucleons. We show that (in the absence of two-body currents) robust corrections for this contamination are obtained from the single-nucleon momentum distribution. These results present an opportunity to compute certain neutrino cross sections with ab-initio methods with controlled uncertainties.

Nikolakopoulos, A. [Washington U., Seattle]↗

Computing neutrino cross sections from Euclidian responses

Energy integrated neutrino cross sections are integrals of nuclear responses weighted with kinematic prefactors. We decompose the prefactors into a limited set of functions of energy transfer and show the relevant integrals are the moments of the responses, and integrals weighted with $1/(a+ω)^n$ with $n\leq 2$. These can be directly obtained from the Euclidean response, avoiding the need for inversion of the Laplace transform. As a proof of concept we study the procedure with toy-model responses for the quasielastic peak. We show that the different contributions can be straightforwardly organized in terms of relative importance, and how flux-averaged cross sections can be obtained. Using a realistic model for the response and numerical uncertainty we show that it is feasible to obtain the required integrals from the Euclidean response, with large uncertainties only for the third moment. Due to kinematic restrictions, the integrals contain contributions from the unphysical region for neutrino scattering, coming from high-momentum nucleons. We show that (in the absence of two-body currents) robust corrections for this contamination are obtained from the single-nucleon momentum distribution. These results present an opportunity to compute certain neutrino cross sections with ab-initio methods with controlled uncertainties.

Nikolakopoulos, A. [Washington U., Seattle]↗

Sparse Cholesky factorization for solving nonlinear PDEs via Gaussian processes

In recent years, there has been widespread adoption of machine learning-based approaches to automate the solving of partial differential equations (PDEs). Among these approaches, Gaussian processes (GPs) and kernel methods have garnered considerable interest due to their flexibility, robust theoretical guarantees, and close ties to traditional methods. They can transform the solving of general nonlinear PDEs into solving quadratic optimization problems with nonlinear, PDE-induced constraints. However, the complexity bottleneck lies in computing with dense kernel matrices obtained from pointwise evaluations of the covariance kernel, and its partial derivatives, a result of the PDE constraint and for which fast algorithms are scarce. The primary goal of this paper is to provide a near-linear complexity algorithm for working with such kernel matrices. We present a sparse Cholesky factorization algorithm for these matrices based on the near-sparsity of the Cholesky factor under a novel ordering of pointwise and derivative measurements. The near-sparsity is rigorously justified by directly connecting the factor to GP regression and exponential decay of basis functions in numerical homogenization. We then employ the Vecchia approximation of GPs, which is optimal in the Kullback-Leibler divergence, to compute the approximate factor. This enables us to compute ϵ-approximate inverse Cholesky factors of the kernel matrices with complexity O(N log d (N/ϵ)) in space and O(N log 2d (N/ϵ)) in time. We integrate sparse Cholesky factorizations into optimization algorithms to obtain fast solvers of the nonlinear PDE. We numerically illustrate our algorithm’s near-linear space/time complexity for a broad class of nonlinear PDEs such as the nonlinear elliptic, Burgers, and Monge-Ampère equations. In summary, we provide a fast, scalable, and accurate method for solving general PDEs with GPs and kernel methods.

97 MATHEMATICS AND COMPUTING↗

A universal and facile approach for building multifunctional conjugated polymers for human-integrated electronics

Polymer semiconductors have shown distinct promise for the development of human- integrated electronics, owing to their solution processability and mechanical softness. However, numerous functional properties required for this application domain face synthetic challenges to be imparted onto conjugated polymers and thus combined with efficient charge- transport property. Here, we develop a “click-to-polymer” (CLIP) synthesis strategy for conjugated polymers, which uses a click reaction for the facile and versatile attachment of diverse types of functional units to a pre-synthesized conjugated-polymer precursor. With four types of functional groups, we show that functionalized polymers from this CLIP method can still retain good charge-carrier mobility. Here, we take two realized polymers to showcase the photo-patternable property and biochemical sensing function, both of which advance the state of the art of realizing these two types of functions on conjugated polymers. We expect the expanded use of this synthesis approach can largely enrich the functional properties from conjugated polymers.

36 MATERIALS SCIENCE↗

Managing Subsurface Pressure Buildup and Interference in Commercial-Scale CO 2 Storage Project with Proximal Injection Wells

Large-scale decarbonization using carbon capture and storage (CCS) is likely to involve many commercial-scale CO 2 storage projects located in close proximity to each other. This close proximity raises concerns over pressure interference among the storage projects. Pressure interference between injection and storage efforts can reduce the practicable CO 2 storage resource and force wells to inject CO 2 at a lower rate to avoid the fracture pressure thresholds per United States Environmental Protection Agency (EPA) Class VI well regulations to preserve injection and confining zone integrity and potentially mitigate against inducing seismic activity. These analyses employ numerical full-physics reservoir modeling to evaluate how pressure buildup fronts and CO 2 plumes evolve under commercial-scale injection volumes of CO 2 in which multiple storage sites located in close proximity occur in tandem. The simulation models mimic injection at pseudo basin-scale and assume homogeneous saline formation(s) as storage targets with a pair of upper and lower homogeneous seal layer/s. These analyses specifically investigate the efficacy of two basin-wide reservoir pressure management strategies in addressing the technical challenges associated with pressure buildup and CO 2 plume commingling. The strategies explored include: 1) enlarging the area of injection well spacing (WS) and 2) storing CO 2 in a stacked sequence (SSS) of saline formations compared to a single formation. The storage and confining zones properties assumed were held common across the scenarios, unless specified otherwise. Analyses results show that after injecting 4 million tons per year for 30 years using 4 separate wells (each injecting 1 million metric tons per year), the radius of CO 2 plume extends to a mere 3 km or less from injection wells. Meanwhile, the radius of pressure buildup ranges on the order of tens to a few hundreds of kilometers, depending on the magnitude of pressure buildup threshold that one would use to define the front. CO 2 plume commingling from different injection wells appears to occur 50 years post-injection, especially under scenarios with narrowly spaced (i.e., < 5 km apart) injection well locations. Findings from sensitivity cases on the well spacing suggest that storage formations modeled would require different well spacing to avoid fracture pressure thresholds. For instance, modeled storage formations with high fracture gradients (i.e., 0.8 psi/ft) would need less than 5–km well spacing, whereas those with lower fracture gradients (i.e., 0.7 psi/ft) would need approximately 20–km well spacing, based on assumed modeling parameters. Under stacked injection, the pressure challenges (described above) still exist but are more alleviated due to distributing the same injection volume across more available reservoir volume. These analyses demonstrate that stacked-sequence storage can effectively address the challenges, while still providing the same target CO 2 storage volumes and allowing a large number of storage projects to be deployed in the same basin by better utilizing the available storage resource across different reservoir depths. Among cases modeled, the resulting pressure buildup front is most suppressed when each storage project distributes injection volumes over several wells, each of which injects a portion of the total CO 2 across the stacked sequence. This strategy results in the smallest CO 2 aerial footprint amongst scenarios evaluated but also shows the largest reduction in the pressure buildup at the top of perforation at the injection wells (upwards of approximately 42 percent compared to the commercial-scale single-formation storage), the result of which is crucial to maintain caprock integrity. The findings presented by this research draw attention to the importance of greater coordination among storage operators and regulatory stakeholders to foster the upscaling and deployment of CCS. These analyses provide insights into required decision-making when considering multi-project deployment in a shared basin. Because these analyses evaluate a very specific geologic situation, they bear further investigations across other geologic situations.

42 ENGINEERING↗

Enhancement of phase change material hysteresis model: A case study of modeling building envelope in EnergyPlus

Nowadays, buildings are expected to offer demand side services to the power grid to enhance the electrical load flexibility, which leads to the concepts of grid-interactive efficient buildings (GEBs). Phase change material (PCM)-based thermal energy storage has seen increasing attention in recent years for peak load shifting of grid-interactive efficient buildings (GEBs). Numerical models are critical tools for design and evaluation of PCM-integrated systems. Most industrial-grade PCMs are reported to melt/freeze over a temperature range instead of at a unique temperature. Such thermal hysteresis effect significantly affects the reliability of simulation results because not only the heat transfer process depends on melting and freezing temperatures, the PCM thermal properties change significantly during the phase change process as well. This study is aimed to develop a model for the PCMs used in the building envelope with the capability to accurately simulate hysteretic behaviors. Further, this model is based on a two-phase assumption and is implemented in a whole building energy performance simulation program (i.e., EnergyPlus). A comparison between numerical results and experimental data shows that during a complete phase transition, the two-phase model could achieve a good agreement with the experimental data. During a partial phase transition, the two-phase model could lead to significant improvements compared to other alternative PCM models, including the existing PCM model in EnergyPlus. Last, whole building simulations were performed to study this model's performance regarding heating/cooling loads and zone mean air temperature of a given building. The results show that the difference in hourly heating/cooling loads introduced by the models was less than 1% in design conditions, while significant changes were observed in both hourly heating/cooling loads and zone mean air temperature when the PCM envelope underwent partial phase transition processes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Conservative High-Order Time Integration for Lagrangian Hydrodynamics

In tis study, we develop novel time integration methods for the compressible Euler equations in the Lagrangian frame that are of arbitrary high order and exactly preserve the mass, momentum, and total energy of the system. The equations are considered in nonconservative form, that is, common for staggered grid hydrodynamics (SGH) methods; namely, the evolved quantities are mass, momentum, and internal energy. A general family of time integration schemes is formulated, and practical pairs for orders three and four are derived. Numerical results on standard hydrodynamics benchmarks confirm the high-order convergence on smooth problems and the exact numerical preservation of all physically conserved quantities.

97 MATHEMATICS AND COMPUTING↗

Harnessing the predicted maize pan-interactome for putative gene function prediction and prioritization of candidate genes for important traits

Abstract The recent assembly and annotation of the 26 maize nested association mapping population founder inbreds have enabled large-scale pan-genomic comparative studies. These studies have expanded our understanding of agronomically important traits by integrating pan-transcriptomic data with trait-specific gene candidates from previous association mapping results. In contrast to the availability of pan-transcriptomic data, obtaining reliable protein–protein interaction (PPI) data has remained a challenge due to its high cost and complexity. We generated predicted PPI networks for each of the 26 genomes using the established STRING database. The individual genome-interactomes were then integrated to generate core- and pan-interactomes. We deployed the PPI clustering algorithm ClusterONE to identify numerous PPI clusters that were functionally annotated using gene ontology (GO) functional enrichment, demonstrating a diverse range of enriched GO terms across different clusters. Additional cluster annotations were generated by integrating gene coexpression data and gene description annotations, providing additional useful information. We show that the functionally annotated PPI clusters establish a useful framework for protein function prediction and prioritization of candidate genes of interest. Our study not only provides a comprehensive resource of predicted PPI networks for 26 maize genomes but also offers annotated interactome clusters for predicting protein functions and prioritizing gene candidates. The source code for the Python implementation of the analysis workflow and a standalone web application for accessing the analysis results are available at https://github.com/eporetsky/PanPPI.

Genetics & Heredity↗