Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulator 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,297 records · Page 72

Correlations for the specific heat capacity of ( U x Pu 1 - x ) 1 - y Gd y O 2 - z derived from molecular dynamics

We report UO 2 is the primary conventional fuel used in most nuclear reactors with Gd 2 O 3 commonly added as a burnable absorber to produce a more level power distribution in the reactor core at the beginning of operation. It can also be mixed with other actinide oxides to produce mixed oxide (MOx) fuel. In this study, molecular dynamics simulations were used to predict the specific heat capacity of Gd-doped PuO 2 , UO 2 and (U, Pu)O 2 MOx accommodating Gd 3+ substituted at cation sites via two charge compensation mechanisms - oxygen vacancy formation and the oxidation of U 4+ to U 5+ . The specific heat capacity values for PuO 2 and UO 2 are in good agreement with other studies showing a distinct peak at high temperatures - above 1800 K. As Gd 3+ is added, the peak height reduces for each composition considered. An analytical fit was applied to the data where Gd 3+ was fully charge compensated by either oxygen vacancies or U 5+ . The expression was then validated by predicting the specific heat capacity for three compositions of (Ux Pu 1-x ) 1-y Gd y O 2-z containing both oxygen vacancies and U 5+ , and compared to molecular dynamics data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Simulating Fracture Sealing by Granular LCM Particles in Geothermal Drilling

Lost circulation occurs when the returned fluid is less than what is pumped into the well due to loss of fluid to pores or fractures. A lost-circulation event is a common occurrence in a geothermal well. Typical geothermal reservoirs are often under-pressured and have larger fracture apertures. A severe lost-circulation event is costly and may lead to stuck pipe, well instability, and well abandonment. One typical treatment is adding lost-circulation materials (LCMs) to seal fractures. Conventional LCMs fail to properly seal fractures because their mechanical limit is exceeded at elevated temperatures. In this paper, parametric studies in numerical simulations are conducted to better understand different thermal effects on the sealing mechanisms of LCMs. The computational fluid dynamics (CFDs) and the discrete element method (DEM) are coupled to accurately capture the true physics of sealing by granular materials. Due to computational limits, the traditional Eulerian–Eulerian approach treats solid particles as a group of continuum matter. With the advance of modern computational power, particle bridging is achievable with DEM to track individual particles by modeling their interactive forces between each other. Particle–fluid interactions can be modeled by coupling CFD algorithms. Fracture sealing capability is investigated by studying the effect of four individual properties including fluid viscosity, particle size, friction coefficient, and Young’s modulus. It is found that thermally degraded properties lead to inefficient fracture sealing.

15 GEOTHERMAL ENERGY↗

Scattering Matrix Determination in Crystalline Materials from 4D Scanning Transmission Electron Microscopy at a Single Defocus Value

Recent work has revived interest in the scattering matrix formulation of electron scattering in transmission electron microscopy as a stepping stone toward atomic-resolution structure determination in the presence of multiple scattering. Here, we discuss ways of visualizing the scattering matrix that make its properties clear. Through a simulation-based case study incorporating shot noise, we shown how regularizing on this continuity enables the scattering matrix to be reconstructed from 4D scanning transmission electron microscopy (STEM) measurements from a single defocus value. Intriguingly, for crystalline samples, this process also yields the sample thickness to nanometer accuracy with no a priori knowledge about the sample structure. The reconstruction quality is gauged by using the reconstructed scattering matrix to simulate STEM images at defocus values different from that of the data from which it was reconstructed.

36 MATERIALS SCIENCE↗

Dynamic Ride-Matching for Large-Scale Transportation Systems

Efficient dynamic ride-matching (DRM) in large-scale transportation systems is a key driver in transport simulations to yield answers to challenging problems. Although the DRM problem is simple to solve, it quickly becomes a computationally challenging problem in large-scale transportation system simulations. Therefore, this study thoroughly examines the DRM problem dynamics and proposes an optimization-based solution framework to solve the problem efficiently. To benefit from parallel computing and reduce computational times, the problem’s network is divided into clusters utilizing a commonly used unsupervised machine learning algorithm along with a linear programming model. Then, these sub-problems are solved using another linear program to finalize the ride-matching. At the clustering level, the framework allows users adjusting cluster sizes to balance the trade-off between the computational time savings and the solution quality deviation. A case study in the Chicago Metropolitan Area, U.S., illustrates that the framework can reduce the average computational time by 58% at the cost of increasing the average pick up time by 26% compared with a system optimum, that is, non-clustered, approach. Another case study in a relatively small city, Bloomington, Illinois, U.S., shows that the framework provides quite similar results to the system-optimum approach in approximately 62% less computational time.

33 ADVANCED PROPULSION SYSTEMS↗

Achilles and GENIE Model Comparisons Between Quasi-Elastic Proton Exclusive Events in DUNE PRISM ND-LAr

In this study, we analyze simulated CCQE proton exclusive events in DUNE s ND-LAr using Achilles and GENIE Monte Carlo event generators. We generate events using muon neutrino exclusive flux files of DUNE s on-axis and 20m off-axis PRISM positions and compare various kinematics of proton exclusive events between the event generators. Events are further analyzed in categories according to the number of final-state protons. This study was motivated by QE s dominance in statistics for off-axis events over other channels in which off-axis events are very useful for characterizing DUNE s flux. Furthermore, our motivation for comparing results between generators is driven by off-axis QE neutrino spectra peaking below 1 GeV. In this region, many low-energy proton events are produced in which models have been shown to have disagreements.

43 PARTICLE ACCELERATORS↗

Collaborative Proposal: Improving understanding of the internal structure and dynamics of deep convection using ARM observations and large eddy simulations

Recent observational and large eddy simulation (LES) modeling studies have nearly unanimously supported the view of deep cumulus convection being composed of a series of quasi-spherical bubbles of buoyant air, known as moist thermals. Despite the prevalence of moist thermals in deep convection, a comprehensive theory for the dynamics of these structures is lacking. Most current conceptual models for cumulus convection are based on canonical scaling theories for dry thermals or plumes; however, there is considerable evidence that the behavior of moist thermals differs markedly from these theories. Furthermore, the theoretical basis for most cumulus parameterizations originates from the plume conceptual model, and therefore these parameterizations are inconsistent with the real structure of moist convection. Motivated by the aforementioned knowledge gaps, this “end-to-end” research effort use theory, observations, numerical simulations, and direct improvements to the Zhang-McFarlane (ZM) convection scheme in the global climate Community Atmosphere Model (CAM) to address the following research questions: What key environmental parameters determine whether or not shallow convection will transition into deep convection, in the context of thermal-like updrafts? What factors regulate the size of thermals within cumulus updrafts? How does vertical wind shear influence thermal behavior, and as a consequence, vertical velocity and mass flux profiles and the shallow-to-deep convective transition? What are the critical processes that determine updraft vertical velocities and their connection to the vertical mass flux profile for thermal-like updrafts? Idealized LES modeling will be used in conjunction with theoretical models for the core properties of thermal-like updrafts to better understand key processes that regulate thermal ascent rates and entrainment properties. Thermal-tracking procedures will be used to characterize the behavior of thermals within the LES, and recently developed direct measures of entrainment and detrainment will be used to quantify entrainment/detrainment rates. Building from these results, we will analyze the structure of moist thermals from hemispheric range-height indicator scans taken during the Atmospheric Radiation Measurement Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign, and from “real case” LES of CACTI events. This combined modeling and observational analysis will provide essential validation for the existing body of research on moist thermal dynamics, which is based primarily on modeling studies. With the insight gained from the aforementioned activities, we will modify the Zhang-McFarlane convection scheme to improve its representation of updraft vertical velocity and entrainment rate profiles. These process-level changes will be tested in the Community Atmosphere Model to assess the impact on global climate simulations.

54 ENVIRONMENTAL SCIENCES↗

The impact of coupled air–sea interaction on extreme East Asian summer monsoon simulation in CMIP5 models

Abstract In this study, the relationship between the ability to simulate air–sea interactions over the western North Pacific (WNP), and to reproduce the extreme East Asian summer monsoon (EASM), were investigated by comparing the performances of several global climate models (GCMs). High ranked in air–sea interaction simulation (HRA) and low ranked in air–sea interaction simulation (LRA) models were selected, according to their performance in simulating relations between sea surface temperature (SST) and precipitation over the WNP, from the ensemble of models that participated in the third and fifth phases of the Coupled Model Intercomparison Project (CMIP3, CMIP5). Compared with CMIP3 models, CMIP5 models exhibited improved simulations of the distinctive air–sea interaction over the WNP, namely, the strong atmospheric forcing on the ocean. Among CMIP5 models, HRA models, which reproduced intrinsic negative correlations between precipitation and SST over the WNP, could simulate the extreme EASM better than LRA models. In particular, HRA models generated a more realistic spatial distribution of the extreme EASM compared with LRA models. The defects of the LRA models resulted from distorted synoptic fields, including underestimated geopotential height and overestimated low‐level wind over the WNP, inducing unrealistic moisture supply and convection due to the exaggerated SST forcing. In contrast, reasonable air–sea interactions represented in HRA models lead to realistic synoptic fields over the WNP, and proper simulation of the extreme EASM.

Kim, Taehyung↗

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↗

Real-Space Local Dynamics of Molten Inorganic Salts Using Van Hove Correlation Function

Molten inorganic salts are attracting resurgent attention because of their unique physicochemical properties, making them promising media for next-generation concentrating solar power systems and molten salt reactors. The dynamics of these highly disordered ionic media is largely studied by theoretical simulations, while the robust experimental techniques capable of observing local dynamics are not well-developed. To provide fundamental insights into the atomic-scale transport properties of molten salts, we report the real-space dynamics of molten magnesium chloride at high temperatures employing the Van Hove correlation function obtained by inelastic neutron scattering. Overall, our results directly depict the distance-dependent dynamics of a molten salt on the picosecond time scale. This study demonstrates the capability of the developed approach to describe the locally correlated- and self-dynamics in molten salts, significantly improving our understanding of the interplay between microscopic structural parameters and their dynamics that ultimately control physical properties of condensed matter in extreme environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigating the Interaction between Persistent Slip Bands and Surface Hard Coatings via Crystal Plasticity Simulations

Fatigue cracks often initiate from the surface extrusion/intrusions formed due to the operation of persistent slip bands (PSBs). Suppression of these surface topographical features by hard surface coatings can significantly extend fatigue lives under lower stress amplitudes (i.e., high cycle fatigue), while cracks initiate early in the coating or in the coating–substrate interface under higher stress amplitudes (i.e., low cycle fatigue), deteriorating the fatigue performance. However, both beneficial and detrimental effects of the coatings appear to be affected by the coating–substrate material combination and coating thickness. A quantitative understanding of the role of these factors in the fatigue performance of materials is still lacking. In this study, crystal plasticity simulations were employed to elucidate the dependence of the coating’s effects on two factors—i.e., the coating thickness and loading amplitudes. The results revealed that the thicker coatings more effectively suppress the operation of the PSBs, but generate higher tensile and shear stresses, normal and parallel to the interfaces, respectively, promoting interfacial delamination. The tensile stresses parallel to the interface within the coating, which favors coating fracture, are not sensitive to the coating thickness.

36 MATERIALS SCIENCE↗

Comparative Analyses of Bioequivalence Assessment Methods for In Vitro Permeation Test Data

ABSTRACT For topical, dermatological drug products, an in vitro option to determine bioequivalence (BE) between test and reference products is recommended. In particular, in vitro permeation test (IVPT) data analysis uses a reference‐scaled approach for two primary endpoints, cumulative penetration amount (AMT) and maximum flux ( J max ), which takes the within donor variability into consideration. In 2022, the Food and Drug Administration (FDA) published a draft IVPT guidance that includes statistical analysis methods for both balanced and unbalanced cases of IVPT study data. This work presents a comprehensive evaluation of various methodologies used to estimate critical parameters essential in assessing BE. Specifically, we investigate the performance of the FDA draft IVPT guidance approach alongside alternative empirical and model‐based methods utilizing mixed‐effects models. Our analyses include both simulated scenarios and real‐world studies. In simulated scenarios, empirical formulas consistently demonstrate robustness in approximating the true model, particularly in effectively addressing treatment–donor interactions. Conversely, the effectiveness of model‐based approaches heavily relies on precise model selection, which significantly influences their results. The research emphasizes the importance of accurate model selection in model‐based BE assessment methodologies. It sheds light on the advantages of empirical formulas, highlighting their reliability compared to model‐based approaches and offers valuable implications for BE assessments. Our findings underscore the significance of robust methodologies and provide essential insights to advance their understanding and application in the assessment of BE, employed in IVPT data analysis.

Leon, Sami↗

Completion of initial reduced-order model for flammable refrigerant dispersal in residential spaces

Environmental regulations aimed at reducing global warming impacts of HVAC&R refrigerants have resulted in the phase-out of chlorofluorocarbons (CFCs) in 2010 and hydrochlorofluorocarbons (HCFCs) by 2030 in developed countries with additional restrictions on the use of hydrofluorocarbons (HFCs) set to take effect in 2036. Many of the remaining alternative refrigerants that have lower global warming potential (GWP) and that are suitable for use in HVAC&R systems (e.g., propane, difluoromethane) are flammable to some degree. Flammable refrigerants introduce new challenges and hazards to property and personal safety related to potential deflagration during system maintenance or due to leakage of the refrigerant accumulating in the conditioned space. Standards have been developed to set maximum charge limits for flammable refrigerants in HVAC&R systems; however, existing safety standards and codes still restrict their use. The bodies that maintain and update these codes need publicly available, science-based information to enable credible guidelines for setting safe charge limits for different flammable refrigerants in different HVAC&R applications. In 2016, the Alliance for Responsible Atmospheric Policy, the Air-Conditioning Heating and Refrigeration Institute (AHRI), ASHRAE, the U.S. Department of Energy (DOE), and the State of California began efforts to develop such information. As part of this effort, Oak Ridge National Laboratory (ORNL) began the current, ongoing project to examine imposed charge limits for flammable refrigerants and identify reasonable adjustments to these limits when found appropriate. Past tasks under this project have included development of experimentally vetted, computational fluid dynamics (CFD) simulation approaches to study the results of leakage of flammable refrigerants from various HVAC&R systems into different types of commercial and residential spaces. In this report, we discuss recent efforts to develop a predictive model of the flammable volume fraction and accumulated refrigerant mass in a single-room residential space resulting from the leak of a flammable refrigerant from a small room air conditioning (RAC) unit. The eventual goal is development of a model suitable for public release which could simulate a range of scenarios. In discussions with the AHRTI Flammable Refrigerant Subcommittee (FRS) at the beginning of this effort, a total of 9 input parameters were chosen for consideration including room area, room/door opening area, ventilation fan flow rate, unit/leak height, leak area, leak rate, total refrigerant charge, refrigerant molecular weight, and state of the unit fan.

42 ENGINEERING↗

Quantifying Properties for a Mechanistic, Predictive Understanding of Aqueous Impact on Aging of Medium and Low Voltage AC and DC Cabling in Nuclear Power Plants

To address the gap in knowledge in understanding degradation in relevant service conditions, this project aimed to develop a mechanistic, predictive model of medium and low voltage cable failure based on the primary environmental degradation parameters of aqueous immersion time, temperature, and the oxidation extent. To do so, we took a two-fold approach toward evaluating degradation as related to the aqueous condition. First, we evaluated the chemical, mechanical, and electrical properties of polymers that comprise the cabling insulation under varied aqueous conditions at different temperatures. Secondly, we evaluated the performance of insulation materials under similar accelerated aging conditions in addition to the dielectric breakdown of the cabling under submersion conditions. Thermal oxidation and immersion of LDPE thin films were done using a Parr vessel and the extent of oxidation was monitored using the carbonyl index from ATR-FTIR spectrum, a ratio quantifying the carbonyl functionality produced on the LDPE film surface. Increasing the temperature, oxidation time, and oxygen pressure increased the water-vapor permeability and the carbonyl content. Thermogravimetric analysis of the films confirmed that there is indeed an initial weight loss of the smaller molecular weight molecules and volatiles. At higher temperatures, there is a secondary mass loss. Up to 70°C, the mass loss curves were similar until 80°C, when the aged films decreased at a much greater rate with respect to increasing temperature. This was seen in the 10, 50, 90% weight loss temperatures. The mechanical properties of sheet PP and HDPE as well as injection molded LDPE, HDPE, HDPP dog bones were plotted with respect to the time spent in the Parr vessel during thermal oxidation. For the sheet PP and HDPE, the UTS and the modulus of elasticity decreased while the elongation increased. The injection-molded LDPE, HDPE, and HDPP dog bones experienced a similar trend where the changes in the properties were within experimental uncertainty. This analysis of LDPE in dry and immersive oxygen-rich environments is an important baseline for future experiments looking at other degradation mechanisms. A predictive aging model for low-density polyethylene insulative cable housings was developed. This model revealed a combination of physical and chemical processes which lead to an increase in permeability and a decrease in strength of the insulator. In creating this predictive model, a gap of knowledge in the literature was revealed in two areas: understanding how the crystallinity of a polymer changes with time and visualizing the predicted pores formed through the insulator which leads to failure. Developing and using an adapted ATR-FTIR method, crystallinity was monitored and compared to bulk crystallinity found via DSC. Inhomogeneity in crystallinity changes were observed but the FTIR method but was found to be less reliable. Pore visualization was found utilizing electrical impedance spectroscopy saturated aged polyethylene films with synthesized citrate-capped gold nanoparticles, and chloroauric acid precursor. results indicate a decreasing impedance of the polyethylene films caused by the transport of ions through the film. SEM imaging was then performed for elemental analysis of the films and counter electrodes for the presence of ions. Chloride and gold ions were detected; however, no nanoparticles were found. This gives an estimated pore size of at least 0.3 nm through the aged film. Cyclic submergence of polyethylene and polypropylene in aqueous solutions of copper sulfate and Harrison’s solution were examined. It was discovered that aging in these mixed conditions and at 90°C did cause a small, yet significant increase to tensile strength for both unaged polyethylene and polypropylene. The PE and PP in cycled and submerged conditions did not have tensile strengths significantly different from dry-aged after 16 weeks. For both solutions, it was determined that capacitance increases with both water tree depth and cable temperature and is relatively independent of changing water tree AR. As for resistance, there is no apparent change between depth percentages of 10 and 80 for both solutions. This is because the water tree has not yet entered the conductor and thus there is no shorting current flowing through the tree yet. Between 80 and 100% water tree depth, it is observed that for both solutions there is a rapid decrease in resistance because the tree is now impacting the conductor shield and allowing shorting current to flow through it. The relationship between resistance and temperature was found to be inversed. As temperature increases, resistance decreases for both water tree solutions with distilled water having the most apparent difference across the range of temperatures simulated. With regards to water tree AR, resistance was found to be relatively independent at depth percentages of 10, 50, and 70, but there was some relationship at 100% depth for both solutions. At this depth, as water tree AR increased, resistance was found to also increase. Through the examination of voltage and electric field distribution plots, it was confirmed for EPR cables that increasing water tree depth leads to increasing distortion. The increase in shorting current at 100% depth for aqueous copper sulfate compared to distilled water was also visualized. Lastly, the relationship between shorting current and temperature for narrow water trees was visualized and validated using the distribution plots. This study has led to a better understanding of the effect of cable temperature, water tree solution, and geometry of the tree region on cable degradation. Furthermore, this study has documented simulation results of water tree degradation in EPR MV cables, which had been previously lacking in information. Future work in this area will add frequency into the mix and examine what effect it has on the rate of cable degradation as a result of water treeing. To compliment the work being carried by UMD on different polymer chemistries, harvested medium voltage (MV) cables were selected for accelerated aging study at ORNL. These medium voltage cables are representative of the vintages and materials that are currently in use at existing nuclear power plants (NPPs). These cables were immersed in water at 90°C and energized for a period of two years at elevated voltage. Partial discharge was measured during this period of time to track potential degradation. Unfortunately, due to the absence of failure, the extent of integration between the UMD techniques and the harvested MV insulation was limited. However, based on the findings and from the UMD research in the previous sections on polyethylene and polypropylene, follow-on characterization for harvested MV cable insulation should focus on sample preparation to take advantage of the UMD techniques to better understand mechanistic degradation in MV cable insulations. This would include: 1.) Utilization of solutions and impedance measurement insulation to track permeability in harvested and accelerated aging insulation samples, 2.) Adaptation of gold nanoparticle synthesis for electrical impedance spectroscopy with supporting SEM characterization to study pore formation in harvested and accelerated aged insulation samples, & 3.) ATR-FTIR crystallinity with multiple temperature dependent harvested MV cable insulation under different air and submerged accelerated aging. In addition, low voltage dielectric spectroscopy and high voltage dissipation factor, tan δ, could be implemented as a condition monitoring tool for harvested MV cable insulation in submerged environments.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Thermal Modeling and Simulation of the Packed-bed Thermal Energy Storage Combined with INL Thermal Energy Distribution System

Dynamic Energy Transport and Integration Lab (DETAIL) at Idaho National Laboratory is to support experimental demonstration and validation research on Nuclear-Renewable Hybrid Energy System [1]. The Thermal Energy Distribution System (TEDS) is a thermal-hydraulic flow loop in DETAIL with its own dedicated control system to support the integration of co-located multiple experimental systems, where a packed-bed thermal energy storage (TES) is installed as a thermal buffer and storage unit. Among various TES options, the packed-bed TES is adopted in TEDS because it offers a low-cost single-tank thermal storage option compared to the traditional two-tank TES. However, since the TES tank is filled with granular fillers having different thermophysical properties from those of TES tank wall, there is a potential thermo-mechanical issue to be carefully addressed like thermal ratcheting which may pose a significant design concern for the packed-bed TES tank. Thermal ratcheting is a thermomechanical process caused by the repeated rearrangement of granular filler inside a TES tank during continuous thermal cycling operation of the packed-bed TES system. If the thermally induced stress exceeds the yield strength of a TES tank wall, catastrophic consequences may happen like rupture of the TES tank. Thus, it is important to understand the phenomenon to ensure the robust operation of the packed-bed TES tank. Given that thermal ratcheting is caused by complex interaction of thermal transport in the porous bed and solid mechanics, the accurate prediction of transient thermal behavior of the packed-bed TES tank, which is the focus of this paper, is critical to the reliable thermal ratcheting analysis. This paper discusses the numerical modeling, simulation, and validation studies that are ongoing at INL to investigate the transient thermal behavior of the packed-bed TES. Of particular concern is the transient thermal process occurring in the packed-bed TES unit that is operated in conjunction with the INL TEDS. The main goal of this research is three-fold: (i) provide preliminary insights into the transient thermal behavior of the TEDS TES tank, (ii) support the thermal measurement and validation plan for TEDS experiment, and (iii) provide transient thermal boundary conditions to support the reliable thermal ratcheting analysis of the TEDS TES tank. For the transient thermal modeling and analysis, a CFD model was developed, and the validity of the modeling approach was examined via comparing the numerical simulation results with the experimental data obtained from various design characteristics of packed-bed TES tanks. Then, the present modeling method was applied for the transient thermal analysis of the TEDS TES tank, and the results are discussed along with the potential improvement of data acquisition strategy for the future TEDS experiments for more precise validation study.

25 ENERGY STORAGE↗

Measured and Simulated Carbon Dynamics in Midwestern U.S. Corn‐Soybean Rotations

Abstract Corn ( Zea mays L.) and soybean ( Glycine max [L.] Merr.) production dominate Midwestern U.S. agriculture and impact the regional carbon and nitrogen cycles. Sustaining soil carbon is important for corn‐soybean production (CS); however, quantifying soil carbon changes requires long‐term field measurements and/or model simulations. In this study, changes in soil organic (SOC), inorganic (SIC), and total (TC) carbon; pH; total nitrogen (TN); and net ecosystem production (NEP) were measured in a conventional corn‐soybean rotation in central Iowa. Soil samples ( n = 42; 0–120 cm depth) were collected from two adjacent fields in 2005 and 2016. Eddy‐flux stations set up in the fields continuously monitored NEP from 2005–2016, and net biome production (NBP) was calculated using yield records. The DayCENT model was used to simulate the effects of conventional management practices on soil carbon and calibrated with field‐measured NEP and SOC. Measured soil TC (0–120 cm) decreased by −14.19 ± 6.25 Mg ha −1 , with highest reductions in SOC and SIC ( p < 0.05) at 0–15 and 90–120 cm, respectively. Measured TN decreased by −0.7 ± 0.29 Mg ha −1 with N‐accumulation at 60–90 cm ( p < 0.05). Eddy‐flux NBP decreased by −13.19 ± 0.05 Mg ha −1 . Soil and eddy‐flux records show a carbon reduction by −1.14 ± 0.63 and −1.20 ± 0.06 Mg ha −1 yr −1 , respectively. The validated DayCENT model suggests that all SOC pools declined. We postulate that conventional CS production has adverse effects on C and N dynamics in Midwestern United States.

Dold, C.↗

Development of in-use engine speed/torque heat maps across multiple heavy-duty commercial vehicle vocations

The U.S. Department of Energy (DOE) established the SuperTruck program with the goal of achieving brake thermal efficiency (BTE) greater than or equal to 55% as demonstrated in an operational heavy-duty (HD) diesel engine at a 65-miles-per-hour (mph) cruise point. Beyond the line-haul application, HD engines operate in a wide range of speed and torque conditions that are unlikely to yield the same efficiency under real-world operation. Thereby, the in-use engine heat maps described in this paper are a valuable tool to illustrate whether the engine-efficiency “sweet spot” matches the most frequent operating conditions. In this study, NREL developed engine heat maps to quantify the important operating points for various vocations using our Fleet DNA database of commercial fleet vehicle operations data. These heat maps clearly show that high-frequency operating points vary significantly according to vehicle vocation, while only a few of them match the sweet spot. Beyond the illustration, engine in-use heat maps can also be leveraged to build up reduced-order engine-efficiency models, needed by many rapid powertrain simulations. As case studies, nine reduced-order models – including line-haul truck, transfer truck, transit bus, transit bus with compressed natural gas (CNG) engine, drayage, refuse pickup, local delivery, utility truck, and school bus with CNG engine – using a trust-region reflective algorithm to fit the on-road data extracted based on the engine in-use heat maps.

33 ADVANCED PROPULSION SYSTEMS↗

A systematic study of proto-neutron star convection in three-dimensional core-collapse supernova simulations

ABSTRACT This paper presents the first systematic study of proto-neutron star (PNS) convection in three dimensions (3D) based on our latest numerical fornax models of core-collapse supernova (CCSN). We confirm that PNS convection commonly occurs, and then quantify the basic physical characteristics of the convection. By virtue of the large number of long-term models, the diversity of PNS convective behaviour emerges. We find that the vigour of PNS convection is not correlated with CCSN dynamics at large radii, but rather with the mass of PNS − heavier masses are associated with stronger PNS convection. We find that PNS convection boosts the luminosities of νμ, ντ, $\bar{\nu }_{\mu }$, and $\bar{\nu }_{\tau }$ neutrinos, while the impact on other species is complex due to a competition of factors. Finally, we assess the consequent impact on CCSN dynamics and the potential for PNS convection to generate pulsar magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗

3D micromechanical simulation of PBX composites

Previous research studied the constitutive response and interface strength of a bonded bicrystal system containing two HMX explosive grains bound with a HTPB polymer binder by performing FE simulations on real geometry imaged with micro Computed Tomography (µCT). The parameters generated from this past study were successfully applied to 2D mesoscale simulations of µCT imaged HMX-HTPB polycrystal samples. Presently, the mechanical response of PBXs containing a Nitro-Plasticized Estane (NPE) binder with similar microstructural geometry to the HMX-HTPB polycrystal samples were studied. Imaging of a HMX-NPE bicrystal sample allowed for simulation of the new material pairing. However, imaging using µCT techniques on the new polycrystal samples was difficult due to the decrease in contrast between the HMX grains and the NPE binder. Mesoscale polycrystal simulations utilizing a single µCT imaged microstructure were performed on a 3D representative volume element (RVE) comparing the delamination behavior and mechanical response of both HMX-HTPB and HMX-NPE PBXs. Furthermore, this technique shows the ability to virtually explore the mechanical response of a range of hypothetical materials that share common microstructural geometric characteristics in addition to studying numerically altered geometries.

36 MATERIALS SCIENCE↗