SEARCH · Engineering Papers

Results for “Effective temperature model”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,243 records · Page 69

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Impacts of irregularly-distributed acidified brine flow on geo-chemo-mechanical alteration in an artificial shale fracture under differential stress

The efficacy of geological carbon sequestration is reliant on the integrity of the caprock and its resistance to physical and chemical alteration. Caprocks with high abundance of reactive carbonates like calcite are susceptible to acid-promoted dissolution and can result in structural weakening. This work investigates the effect of acidified brine flow through an artificially fractured, high-carbonate (30 % by XRD) shale under differential compressive stress. Cylindrical samples were cut in half vertically and milled to create an artificial fracture with interlocking asperities and open channels. Samples were sheared with a single applied stress in a custom flow cell housed within an industrial CT scanner. Further, either acidic (pH 4) or reservoir-simulated (pH 9.5) brine was flowed through the artificial fracture for 7–8 days under reservoir pressure and room temperature. Model simulations indicate flow mainly occurred in open channels, with limited flow between overlapping asperities. Analysis of fracture surfaces by optical and scanning electron microscopy show increased surface alteration and roughness after exposure to pH 4 versus pH 9.5 brine indicating mineral dissolution/loss, and this effect is greater in areas that receive the highest brine flows. Similarly, surface analysis by scratch testing shows fracture toughness decreases more after exposure to acidic versus reservoir-simulated brine, with the greatest alteration in areas of highest acidic brine flows. Despite weakening, no shear slip occurred. Overall, the results indicate that acidified brine can result in significant physical and geomechanical alteration of irregular fracture surfaces in shale caprock, with greatest effects in preferential flow regions.

58 GEOSCIENCES↗

DOE OSTI2024

Autoregressive distributed lag-based dynamic uniformity modeling and monitoring approaches for superconductor manufacturing

High-temperature superconductors (HTS), known for their high efficiency and low energy loss, have found profound applications across various fields, driving the demand for long, uniformly performing tapes. However, ensuring uniform performance over extended lengths of HTS tapes, often characterized by the consistency of critical current, remains challenging due to fluctuations in growth conditions during manufacturing. Here, to elucidate the mechanisms underlying variations in tape uniformity and enable real-time monitoring of associated parameters, we propose an Autoregressive Distributed Lag (ADL)-based Dynamic Uniformity Modeling and Monitoring (ADUM2) approach. This method integrates uniformity measurement, the identification of critical process parameters and real-time monitoring within the manufacturing process. The ADUM2 approach is applied to the advanced metal organic chemical vapor deposition (A-MOCVD) process, a pilot-scale method for superconductor manufacturing. Our model demonstrates superior performance compared to benchmark methods, accounting for over 80% of the total variance in the data and identifying 13 key process parameters influencing the uniformity of HTS tapes. This study offers significant insights into the high-temperature superconductor manufacturing process and holds the potential to facilitate the production of cost-effective, uniformly performing long superconducting tapes in the future.

autoregressive distributed lag analysis↗

DOE OSTI2026

High-Resolution Regional Atmosphere–Ocean–Wave Coupled Simulations of Hurricane Henri (2021)

To explore the integrated effects of ocean and ocean surface wave related physical processes on TC simulations, a set of three model simulations is performed. * In experiment 'A', the event is modeled using the Weather Research Forecasting (WRF) model alone with prescribed Sea Surface Temperature (SST) at 6-hour intervals. * In experiment ‘AO,’ WRF is coupled with the Finite Volume Community Ocean Model (FVCOM), enabling variable exchange between atmosphere and ocean, but without considering ocean surface wave-related physical processes. * In experiment ‘AOW’, WRF, FVCOM, and Simulating WAves Nearshore (SWAN) exchange variables with each other every hour through the OASIS3-MCT Coupler to allow direct and indirect atmosphere-ocean-wave interactions. * Observational data are also included in this dataset (Dropsonde, HRD-Radar, NDBC_wave). All simulations are initialized at 18:00 UTC on August 19, 2021, within a domain encompassing the western North Atlantic Ocean. The atmospheric domain features a horizontal resolution of 3 km. The ocean domain, which covers a substantial portion of the WRF ocean domain, employs an unstructured triangular grid with resolutions ranging from 3 km near the coast to 9 km in the open ocean, effectively resolving the complex coastline of the U.S. Northeast Coast. Initial and boundary conditions for the atmosphere model are obtained from the 6-hourly 0.25° NCEP (National Centers for Environmental Prediction) Global Forecast System (GFS; NCEP, 2015) data.

17 WIND ENERGY↗

DOE OSTI2026

High-Resolution Regional Atmosphere–Ocean–Wave Coupled Simulations of Hurricane Henri (2021)

17 WIND ENERGY↗

DOE OSTI2022

Graphite Oxidation Behavior: NRC Graphite Behavior Model

Graphite microstructure and oxidation Graphite crystal structure and basal planes Reactive surface area sites Zig-zag and Arm-chair Pore microstructure effects Interconnectivity Reactivity and oxygen transport Oxygen-graphitereaction Not C-O reaction! Temperature effects Oxygen diffusion to graphite Diffusion into the microstructure High temperature vs low temperature diffusion Oxidation: It’s a gradient of factors

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Dynamic modelling and control strategy of a temperature-driven metal hydride cooling system for buildings

A temperature-driven coupled metal hydride (MH) based thermal energy storage (TES) system can allow to shave and shift the peak energy demand in buildings. The high energy density and long-term (seasonal) energy storage capability are its major advantages over other energy storage methods. The dynamic nature of the MH operation, however, requires controlled hydrogen transfer between the coupled MHs at a rate needed to meet the building's transient load. While temperature-driven MH systems are studied in the literature, their application in buildings and control are scarcely reported. Here, this paper presents a control-based dynamic modeling of the temperature-driven coupled MH-TES system for building cooling applications. The dynamic model is developed in MATLAB(R) Simulink environment, considering the thermodynamic and kinetic behaviors of the MH systems. Based on a preliminary analysis of a property database of over 337 hydrides, we select around 1600 MH pairs suitable for building cooling applications. Each of these MH pairs is studied for their performance using the dynamic model, and among all, Zr 0.76 Ti 0.24 Ni 1.16 Mn 0.63 V 0.14 Fe 0.18 -Ti 0.85 Zr 0.15 Cr 1.2 Mn 0.8 MH pair showed fast dynamics along with high coefficient of performance (COP) of 0.71. A parametric investigation is performed on this MH pair to understand the effect of operating temperatures. Finally, three proportional-integral (PI) feedback controllers are investigated to regulate the temperature, pressure and mass exchange between the coupled MH pairs. The developed PI controller is sufficiently capable of rejecting the signal noise from the hydrogen flow and internal heat exchange processes with root mean square error of 5.78 W between reference and actual cooling load.

08 HYDROGEN↗

DOE OSTI2020

Ion-Exchange Modeling of Crystalline Silicotitanate for Cesium Removal - 20283

The Tank Closure Cesium Removal (TCCR) system is a Savannah River Site (SRS) demonstration 'at-tank' process designed to remove {sup 137}Cs from the high-level aqueous tank waste so that the decontaminated solution can be disposed as low-level waste. Cesium is removed by ion exchange (IX) columns using engineered IONSIV{sup TM} R9120-B form of the Crystalline Silicotitanate (CST) media. The TCCR system is deployed at Tanks 10 and Tank 11 in the SRS H Tank Farm. Water is added to the salt-cake in Tank 10 H to dissolve it. The dissolved salt solution waste is pumped out of Tank 10H (feed tank), through filters and IX columns. The decontaminated salt solution is transferred to Tank 11 (receipt tank), and on to Tank 50H for final disposal in the Saltstone Production Facility. The current TCCR can accommodate lead-lag (two-column) or lead-lag-guard (three-column) configurations to optimize media utilization and achieve the target decontamination. To assist the TCCR operations, a parametric study was conducted to evaluate the impact of different parameters (e.g., column configurations (single column, two columns or three columns in series), waste characteristics, operating temperature, process flow rate, CST average particle size) on the IX column performance including CST bed utilization. The initial results indicate that the IX column performance is improved at slower process flow rate, at lower operating temperature, and with smaller CST average particle size. Multi-column configurations are recommended, because the single-column configuration does not utilize CST bed effectively. This paper demonstrates the versatility of the ion exchange modeling to evaluate the effects of CST characteristics and operational parameters on IX column performances. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Medium-Scale Methanol Pool Fire Model Validation

In this work, medium scale (30 cm diameter) methanol pool fires were simulated using the latest fire modeling suite implemented in Sierra/Fuego, a low Mach number multiphysics reacting flow code. The sensitivity of model outputs to various model parameters was studied with the objective of providing model validation. This work also assesses model performance relative to other recently published large eddy simulations (LES) of the same validation case. Two pool surface boundary conditions were simulated. The first was a prescribed fuel mass flux and the second used an algorithm to predict mass flux based on a mass and energy balance at the fuel surface. Gray gas radiation model parameters (absorption coefficients and gas radiation sources) were varied to assess radiant heat losses to the surroundings and pool surface. The radiation model was calibrated by comparing the simulated radiant fraction of the plume to experimental data. The effects of mesh resolution were also quantified starting with a grid resolution representative of engineering type fire calculations and then uniformly refining that mesh in the plume region. Simulation data were compared to experimental data collected at the University of Waterloo and the National Institute of Standards and Technology (NIST). Validation data included plume temperature, radial and axial velocities, velocity temperature turbulent correlations, velocity velocity turbulent correlations, radiant and convective heat fluxes to the pool surface, and plume radiant fraction. Additional analyses were performed in the pool boundary layer to assess simulated flame anchoring and the effect on convective heat fluxes. This work assesses the capability of the latest Fuego physics and chemistry model suite and provides additional insight into pool fire modeling for nonluminous, non-sooting flames.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

LLNL FESP Theory Highlights: November 2024

This talk presented a refined understanding of plasma-facing thermionic emitter cathodes by examining the limitations of the conventional space charge-limited model, which overlooks ion trapping dynamics within the virtual cathode. While conventional theory distinguishes between temperature-limited and space charge-limited regimes based on electron current constraints, recent findings reveal that trapped ions can form a quasi-neutral "second plasma," expanding upstream and interacting with the primary plasma. These effects are amplified in non-planar geometries, such as cylindrical or spherical sheaths, where trapped ions create a positive feedback loop with the upstream plasma. To address these phenomena, a new "Aid-and-Compete" model was developed, offering insights into current enhancement and sheath mode transitions relevant for plasma devices, including emissive cathodes in electric propulsion, fusion applications, thin film deposition systems, and plasma diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Multiscale modeling of hydrogenolysis of ethane and propane on Ru(0001): Implications for plastics recycling

Plastic waste presents an environmental threat. Chemical recycling via hydrogenolysis can convert plastic waste into waxes, lubricants, and fuels. Among catalysts, Ru stands out for its superior activity and selectivity. The chemistry of light alkane hydrogenolysis can help understanding plastics deconstruction. Here, we perform first-principles calculations, develop descriptor-based relations, and conduct microkinetic modeling and analysis on ethane and propane. Predictions are in excellent agreement with experimental data. We identify a similar cracking pattern for both hydrocarbons entailing a deeply dehydrogenated species with the removal of four hydrogen atoms: CHCH*+* → 2CH* for ethane and CH 3 CCH*+2* → CH 3 C* + CH* for propane. We find that the rate-determining step is the C-C cracking for ethane and the first dehydrogenation from the terminal carbon (CH 3 CH 2 CH 3 *+* → CH 3 CH 2 CH 2 *+H*) for propane. The vinyl species CH 2 CH* produced from propane cracking is responsible for whether a single or multiple cracking events occur and effectively controls the selectivity. Specifically, ethane formation in propane hydrogenolysis is suppressed at elevated temperatures due to over-cracking via multiple (two here) cracking events being preferred over hydrogenation and desorption of ethane from the catalyst. Our workflow and models provide a baseline for future studies on heavier hydrocarbons. Insights into recent experimental studies of polyethylene over Ru-based catalysts are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Data-Enabled Predictive Control for Building HVAC Systems

Model predictive control is widely used as a control technology for the computation of optimal control inputs of building heating, ventilating, and air conditioning (HVAC) systems. However, both the benefits and widespread adoption of model predictive control (MPC) are hindered by the effort of model creation, calibration, and accuracy of the predictions. In this paper, we apply the data-enabled predictive control (DeePC) algorithm for designing controls for building HVAC systems. The algorithm solely depends on input/output data from the system to predict future state trajectories without the need for system identification. The algorithm relies on the idea that a vector space of all input–output trajectories of a discrete-time linear time-invariant (LTI) system is spanned by time-shifts of a single measured trajectory, given the input signal is persistently exciting. Closed-loop simulations using EnergyPlus are performed to demonstrate the approach. The simulated building modeled in EnergyPlus is a modified commercial large office prototype building served by an air handling unit-variable air volume HVAC system. Temperature setpoints of zones are used as control variables to minimize the HVAC energy cost of the building considering a time-of-use electricity rate structure. Furthermore, sensitivity analysis is conducted to gain insights into the effect of parameter tuning on DeePC performance. Simulation results are used to illustrate the performance of the algorithm and compare the algorithm with model-based MPC and occupancy-based setpoint controller. Overall, DeePC achieves similar performance compared to MPC for lower engineering effort.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

On the Thermodynamic Condition for Adsorption Azeotropes

Adsorption azeotropy is a common phenomenon in mixed-gas adsorption equilibria. Since the first literature report of binary adsorption azeotropes in 1933, numerous studies investigated the thermodynamic conditions for adsorption azeotropes but failed to reach definitive conclusions. Based on the generalized Langmuir isotherm model for multicomponent adsorption equilibria which takes into account the vacant site as part of the adsorbed phase, this study presents the thermodynamic condition for adsorption azeotropes as derived from the generalized Langmuir isotherm to be the equality of the ratios of adsorbed phase activity coefficient γi and adsorption equilibrium constant $K$$^{o}_{i}$ for the two adsorbates in the binary 1–2 adsorption system, i.e., γ 1 /$K$$^{o}_{1}$ = γ 2 /$K$$^{o}_{2}$. This adsorption azeotropic condition is analogous to the vapor–liquid equilibrium azeotropic condition, i.e., the equality of the products of liquid phase activity coefficient and saturation vapor pressure Pisat for the two components, i.e., γ 1 $P$$^{sat}_{1}$ = γ 2 $P$$^{sat}_{2}$. We validated the thermodynamic condition for adsorption azeotropes with 14 azeotrope-forming adsorption systems in this study. As a result, we investigated the effects of the pressure, temperature, and adsorbed phase nonideality on azeotrope formation.

42 ENGINEERING↗

DOE OSTI2025

Experimental and modeling study of the autoignition behavior of a saturated heterocycle: Pyrrolidine

Experiments are conducted in both rapid compression machine (RCM) and shock tube (ST) to better quantify autoignition behavior (e.g., ignition delay, heat release) and understand heteroatomic effects in heterocyclic compounds, which are important reference components for the combustion of biomass-derived liquid fuels. These tests focus on the nitrogen-containing, five-membered saturated ring, pyrrolidine, at diluted conditions covering pressures of 20 and 50 bar, temperatures of 720–1450 K and a range of stoichiometries (ϕ = 0.5–2). A chemical kinetic model is developed and coupled to an existing combustion kinetics framework describing key nitrogen containing intermediates (e.g. pyrrole, ammonia and NOx). H-abstraction reactions by OH, H, CH 3 and HO 2 , are determined using ab-initio transition state theory methods, while analogies to cyclopentane are adopted for many other reactions, such as ring-opening. Here, the autoignition measurements reveal the lack of negative temperature coefficient (NTC) behavior and low-temperature chemistry for pyrrolidine, as opposed to its saturated hydrocarbon analogue, cyclopentane. Interestingly, at the lowest temperatures (T < 750 K), the reactivity of cyclopentane is greater than pyrrolidine, while at higher temperatures, pyrrolidine becomes more reactive. Agreement between the experimental measurements and the model is good, and it is found that H-abstraction reactions by HO 2 and ensuing chemistry play key roles in controlling the reactivity of this cyclic amine. Most of the fluxes, i.e., >70 %, are predicted to move through 1- or 2-pyrroline (C 4 H 7 N) and then the cyclic C 4 H 6 N radical, at both lower and higher temperatures, to form either CH 2 CHCHCHNH via ring-opening or pyrrole via β-scission. It appears that the ring opens more easily at lower temperature whereas the C–H β-scission dominates at higher temperature and lower pressure, such that the reaction of the fuel radical intermediate carrying an unpaired electron on the nitrogen atom with HO 2 is the next most notable in promoting oxidation. When comparing pyrrolidine and cyclopentane, which exhibits distinct pathways in different temperature regimes, the pyrrolidine pathways and sensitivity analysis align more closely to the high temperature case of cyclopentane where the important role of HO 2 radicals is seen to provide chain branching through HO 2 reaction with the fuel, accompanied by H 2 O 2 formation and decomposition to OH. The formation of 5-membered diene rings and ring opening reactions are also found to be highly relevant. Of particular note, it is found that there is little influence of small molecule nitrogen-chemistry, e.g., NH 2 , HCN, NO/NO 2 on the reactivity of the pyrrolidine mixtures investigated here where no recirculated combustion gases are included.

Autoignition↗

DOE OSTI2024

Scalable and compact magnetocaloric heat pump technology

Magnetocaloric heat pumping (MCHP) promises to be more efficient than traditional vapor compression while also eliminating the deleterious effects of gaseous refrigerants. While MCHP devices have shown the temperature spans and efficiencies needed for different heating and cooling applications, they struggle to become commercially viable due to their large size and mass, and resultant high cost. This paper evaluates a baseline MCHP device and explores methods to boost its system power density (SPD). The key components of the baseline system are the gadolinium packed-particle bed active magnetic regenerator (AMR) and a magnetic source composed of permanent magnets and high permeability magnetic steel. To enhance the SPD, the paper evaluates maximizing the AMR volume, opting for first-order magnetocaloric materials, optimizing the magnet and AMR geometry, and reducing the size of magnets and magnetic steel parts. At larger thermal powers, increasing the AMR diameter and the number of magnetic poles were evaluated. Using finite element models, solid models, and estimates of magnetocaloric material performance, thermal powers ranging from 37 W to 44 kW at a nominal 10 K temperature span were projected, and SPD was estimated to improve from 6 W/kg to 81 W/kg. Neglecting end effects, an upper limit of 114 W/g is estimated. Compared to SPD of off-the-shelf compressors with similar environment temperatures, MCHP power density using gadolinium is competitive up to roughly 200 W of cooling power. This is extended to 1 kW when using LaFeSi alloys and up to 3 kW in the limiting case. In conclusion, these results indicate that the performance and mass of MCHP can match that of compressors, which is a critical step toward cost-competitive magnetocaloric technology.

42 ENGINEERING↗

DOE OSTI2026

Development of a BISON validation case for the TRISO transient irradiations in NSRR using effective heat capacity methods

The current tristructural isotropic (TRISO) fuel assessment and validation database in BISON primarily covers steady- state irradiation and high-temperature furnace testing. Transient assessment cases are potentially needed to support U.S. industry efforts in designing and deploying commercial reactors using TRISO fuels. Historical transient tests in- volving TRISO fuels used highly conservative conditions compared to the typical high-temperature gas-cooled reactor accident scenarios. Despite this, modeling historical transient tests is fundamental for evaluating BISON’s predictive capabilities, adapting material properties for high-temperature and high-particle-power regimes, and developing a sys- tematic validation approach for TRISO transient applications. This work developed a 1D model of transient experiments carried out at the Nuclear Safety Research Reactor using BISON. BISON’s predictions of energy deposition, UO 2 melting onset, and molten volume fractions are compared against experimental measurements. Melting was modeled using an effective specific heat capacity model for UO 2 . We found that BISON’s predictions are in reasonable agreement with experimental data for low-energy-deposition cases, and that BISON overpredicts melting at higher energy depositions. We also discuss the potential causes of discrepancies between the simulated and measured results and propose ways to further develop the model. Although these simulations used conservative conditions compared to those expected for actual TRISO-fueled reactors, they extended the range of conditions reflected in the data in the existing BISON database for TRISO fuels.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

A New Search for Low-mass Dark Matter and an Examination and Reduction of the Uncertainty due to the Photoelectric Absorption Cross Section using a Cryogenic Silicon Detector with Single-charge Sensitivity

Many astrophysical observations point to the abundant existence of dark matter (DM) in the universe composed of particles beyond the Standard Model. However despite numerous experiments using different techniques, DM particles have yet to be directly observed. The Super Cryogenic Dark Matter Search (SuperCDMS) SNOLAB experiment will employ silicon and germanium crystal detectors operated at temperatures as low as 15~mK to probe DM interactions using phonon and ionization signals. Recently, R\&D; facilities have developed gram-sized high-voltage eV-scale (HVeV) silicon detectors that achieve single-electron-hole-pair resolution. By probing effective DM-electron interactions, these devices can be used for searches of low-mass DM candidates.This dissertation presents a DM search experiment known as HVeV Run~2 that employs a second-generation HVeV detector operated in an above-ground laboratory at Northwestern University (IL, USA). Energy spectra are obtained from a blind analysis with 0.39 and 1.2\,g-days of exposure with the detector biased at 60 and 100\,V, respectively. The 0.93~gram detector achieves a 3\,eV phonon energy resolution, corresponding to a world-leading charge resolution of 3\,\% of a single electron-hole pair for a detector bias of 100\,V. With charge carrier trapping and impact ionization effects incorporated into the DM signal models, the resulting exclusion limits are reported for inelastic DM-electron scattering for DM masses from 0.5--$10^{4}$\,MeV$/c^{2}$; in the mass range from 1.2--50\,eV$/c^{2}$ the limits for dark photon and axion-like particle absorption are reported.Several DM search experiments, including HVeV Run~2, are sensitive to low-mass DM candidates that rely on the temperature-dependent photoelectric absorption cross section of silicon. However discrepancies in the underlying literature data result in dominating systematic uncertainties on the DM exclusion limits. In order to reduce these systematic uncertainties, this dissertation presents a novel method of making a direct, low-temperature measurement of the photoelectric absorption cross section of silicon at energies near the band gap (1.2--2.8\,eV).

Wilson, Matthew James↗

DOE OSTI2022

Cloud adjustments from large-scale smoke–circulation interactions strongly modulate the southeastern Atlantic stratocumulus-to-cumulus transition

Smoke from southern Africa blankets the southeastern Atlantic Ocean from June to October, producing strong and competing aerosol radiative effects. Smoke effects on the transition between overcast stratocumulus and scattered cumulus clouds are investigated along a Lagrangian (air-mass-following) trajectory in regional climate and large eddy simulation models. Results are compared with observations from three recent field campaigns that took place in August 2017: ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES), CLouds and Aerosol Radiative Impacts and Forcing: Year 2017 (CLARIFY), and Layered Atlantic Smoke Interactions with Clouds (LASIC). The case study is set up around the joint ORACLES–CLARIFY flight that took place near Ascension Island on 18 August 2017. Smoke sampled upstream on an ORACLES flight on 15 August 2017 likely entrained into the marine boundary layer later sampled during the joint flight. The case is first simulated with the WRF-CAM5 regional climate model in three distinct setups: (1) FireOn, in which smoke emissions and any resulting smoke–cloud–radiation interactions are included; (2) FireOff, in which no smoke emissions are included; (3) RadOff, in which smoke emissions and their microphysical effects are included but aerosol does not interact directly with radiation. Over the course of the Lagrangian trajectory, differences in free tropospheric thermodynamic properties between FireOn and FireOff are nearly identical to those between FireOn and RadOff, showing that aerosol–radiation interactions are primarily responsible for the free tropospheric effects. These effects are non-intuitive: in addition to the expected heating within the core of the smoke plume, there is also a “banding” effect of cooler temperature (~1–2 K) and greatly enhanced moisture (>2 g kg -1 ) at the plume top. This banding effect is caused by a vertical displacement of the former continental boundary layer in the free troposphere in the FireOn simulation resulting from anomalous diabatic heating due to smoke absorption of sunlight that manifests primarily as a few hundred meters per day reduction in large-scale subsidence over the ocean. A large eddy simulation (LES) is then forced with free tropospheric fields taken from the outputs for the WRF-CAM5 FireOn and FireOff runs. Cases are run by selectively perturbing one variable (e.g., aerosol number concentration, temperature, moisture, vertical velocity) at a time to better understand the contributions from different indirect (microphysical), “large-scale” semi-direct (above-cloud thermodynamic and subsidence changes), and “local” semi-direct (below-cloud smoke absorption) effects. Despite a more than 5-fold increase in cloud droplet number concentration when including smoke aerosol concentrations, minimal differences in cloud fraction evolution are simulated by the LES when comparing the base case with a perturbed aerosol case with identical thermodynamic and dynamic forcings. A factor of 2 decrease in background free tropospheric aerosol concentrations from the FireOff simulation shifts the cloud evolution from a classical entrainment-driven “deepening–warming” transition to trade cumulus to a precipitation-driven “drizzle-depletion” transition to open cells, however. The thermodynamic and dynamic changes caused by the WRF-simulated large-scale adjustments to smoke diabatic heating strongly influence cloud evolution in terms of both the rate of deepening (especially for changes in the inversion temperature jump and in subsidence) and in cloud fraction on the final day of the simulation (especially for the moisture “banding” effect). Such large-scale semi-direct effects would not have been possible to simulate using a small-domain LES model alone.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Modeling creep behavior in ceramic matrix composites

Here, in this work, a three-dimensional viscoplasticity formulation with progressive damage is developed and used to investigate the complex time-dependent constituent load transfer and progressive damage behavior in ceramic matrix composites (CMCs) subjected to creep. The viscoplasticity formulation is based on Hill's orthotropic plastic potential, an associative flow rule, and the Norton-Bailey creep power law with Arrhenius temperature dependence. A fracture mechanics-informed isotropic matrix damage model is used to account for CMC brittle matrix damage initiation and propagation, in which two scalar damage variables capture the effects of matrix porosity as well as matrix property degradation due to matrix crack initiation and propagation. The Curtin progressive fiber damage model is utilized to simulate progressive fiber failure. The creep-damage formulation is subsequently implemented as a constitutive model in the generalized method of cells (GMC) micromechanics formulation to simulate time-dependent deformation and material damage under creep loading conditions. The developed framework is used to simulate creep of single fiber SiC/SiC microcomposites. Simulation results are in excellent agreement with experimental and numerical data available in the literature.

36 MATERIALS SCIENCE↗