Techno-economic uncertainty analysis of wet waste-to-biocrude via hydrothermal liquefaction
Not Available
Li, Shuyun↗
SEARCH · Engineering Papers
Results for “Coupled models”
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.
TChem-atm v1.0
SAND2024-11300O TChem-atm is a software library that was developed to solve complex kinetic models for atmospheric chemistry applications. TChem-atm interface employs a hierarchical parallelism design to exploit the massive parallelism available from modern computing platforms. It also supports gas atmospheric chemistry applications, e.g., the energy exascale earth system model. TChem can be used as a box model or coupled with a climate model to compute the time evolution of gas tracer species. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.
Safta, Cosmin↗
SUHMO: an adaptive mesh refinement SUbglacial Hydrology MOdel v1.0
Abstract. Water flowing under ice sheets and glaciers can have a strong influence on ice dynamics, particularly through pressure changes, suggesting that a comprehensive ice sheet model should include the effect of basal hydrology. Modeling subglacial hydrology remains a challenge, however, mainly due to the range of spatial and temporal scales involved – from subglacial channels to vast subglacial lakes. Additionally, networks of subglacial drainage channels dynamically evolve over time. To address some of these challenges, we have developed an adaptive mesh refinement (AMR) model based on the Chombo software framework. We extend the model proposed by Sommers et al. (2018) with a small but significant change to accommodate the transition from unresolved to resolved flow features. We handle the strong nonlinearities present in the equations by resorting to an efficient nonlinear full approximation scheme multigrid (FAS-MG) algorithm. We outline the details of the algorithm and present convergence analysis results demonstrating its good performance. Additionally, we present results validating our approach, using test cases from the Subglacial Hydrology Model Intercomparison Project (SHMIP) (de Fleurian et al., 2018). We finish by presenting a more complex, 100 km-by-100 km synthetic test case with peaks and valleys that we use to investigate the effective pressure distribution as the number of AMR levels increases. These preliminary results suggest that a minimum spatial resolution is needed to properly capture channel features, but additional work is required to precisely quantify this and its impact on accurately modeling the coupled ice sheet–hydrology system. The efficiency of our approach, relying on localized refinement, is also demonstrated. Future work will include coupling the SUbglacial Hydrology MOdel (SUHMO) with the BISICLES AMR ice sheet model (Cornford et al., 2013), both built on the same numerical framework.
58 GEOSCIENCES↗
Near-Surface Hydrology and Soil Properties Drive Heterogeneity in Permafrost Distribution, Vegetation Dynamics, and Carbon Cycling in a Sub-Arctic Watershed: Modeling Archive
This Modeling Archive is in support of a NGEE-Arctic publication: Shirley et al. (2022) “Near-Surface Hydrology and Soil Properties Drive Heterogeneity in Permafrost Distribution, Vegetation Dynamics, and Carbon Cycling in a Sub-Arctic Watershed". [DOI].The dataset contains outputs from the global sensitivity analysis (GSA) of the “ecosys” model as reported in Shirley et al. (2022). The study showed that discontinuous permafrost environments are characterized by complex feedback loops and strong spatial heterogeneity which is created by variability in near-surface hydrology and soil properties. Additionally, the study demonstrated that missing representation of sub-grid heterogeneity in terrestrial ecosystem models can lead to biased estimates of the high-latitude carbon budget. Included in this dataset are the factor values for each run in the GSA and the model outputs used in this study. Included are two *.csv data files and one *.pdf.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).
54 ENVIRONMENTAL SCIENCES↗
A Machine Learning-Assisted Framework to Control Thermally Anisotropic Building Envelopes in Residential Buildings
To curb the energy consumption of buildings and their related CO2 emissions, Oak Ridge National Laboratory (ORNL) has developed the thermally anisotropic building envelope (TABE) —a multi-layer design comprising insulation materials and metal foils connected to thermal loops. In this study, a machine learning-assisted framework was developed to control the TABE in residential buildings to reduce the computation load for future optimal rule-based control and application. First, a 2D finite element model was established in COMSOL to calculate the hourly heat flux through exterior walls installed with the TABE. Then, TABE wall heat fluxes were simulated for various given indoor and outdoor boundary conditions, as well as thermal loops fluid temperatures and flow rates. Since the finite element simulations are computationally expensive, an artificial neural network (ANN) was then trained to use as a proxy of the finite element (COMSOL) modeling. Finally, the trained ANN model was coupled with the EnergyPlus model to predict the energy consumption of a US Department of Energy prototype single-family house installed with the TABE. An optimal simple rule-based control was determined from predefined rules for a case study. The results demonstrate that the developed machine learning–assisted framework can reduce 99.9% of the computation time while efficiently managing residential building energy for installed TABE walls.
Shen, Zhenglai↗
A Unified Analytical Model for Pressure Solution With Fully Coupled Diffusion and Reaction
Abstract Geophysical models for pressure solution are typically developed for diffusion‐controlled or reaction‐controlled scenarios. We present a unified analytical model that considers fully coupled diffusion and reaction during pressure solution. The model recovers the diffusion‐controlled and reaction‐controlled models in the literature as specific limiting cases. When diffusion and reaction exhibit comparable influences, we validate the proposed model against independent numerical simulations. The proposed model is then employed in interpreting experimental measurements, demonstrating a better agreement compared to previous interpretations.
Wang, Ziyan [Energy Geosciences Division Lawrence ↗
A model to separate conduction and radiation in high temperature laser flash measurements for semi-transparent materials
Laser flash analysis (LFA) is a widely used technique in the scientific and industrial communities for measuring the thermal diffusivity and conductivity of materials. However, traditional fitting methods based solely on heat conduction do not account for radiation heat transfer in transparent and semi-transparent materials, which becomes significant at high temperatures. In this study, we have developed a coupled radiation-conduction model that considers the optical properties of materials to accurately extract thermal diffusivity from LFA measurements. We conducted LFA measurements on fused silica glass samples of varying thicknesses and within the temperature range of 25 to 1000 °C and fitted the data with the model. Since the coupled model requires optical properties as input parameters, which are often unavailable for new materials, we have also proposed a modified transparent model. Compared to the coupled model, results from the modified transparent model generally overestimate thermal conductivity, but the errors are acceptable in silica glass with small thickness. Here, we also investigated the applicability of the modified transparent model with materials of different optical properties.
14 SOLAR ENERGY↗
Superposed periodic kink and pulse solutions of coupled nonlinear equations
Here, we present novel previously unexplored periodic solutions, expressed in terms of Jacobi elliptic functions, for both a coupled Φ 4 model and a coupled nonlinear Schrödinger equation (NLS) model. Remarkably, these solutions can be elegantly reformulated as a linear combination of periodic kinks and antikinks, or as a combination of two periodic kinks or two periodic pulse solutions. However, we also find that for $m=0$ and a specific value of the periodicity (or at a nonzero value of the elliptic modulus $m$) this superposition does not hold. These results demonstrate that the notion of superposed solutions extends to the coupled nonlinear equations as well.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Data Release for Retreat and Regrowth of the Greenland Ice Sheet During the Last Interglacial as Simulated by the CESM2-CISM2 Coupled Climate–Ice Sheet Model
CESM2 and CISM2 data files for figures in "Retreat and Regrowth of the Greenland Ice Sheet During the Last Interglacial as Simulated by the CESM2-CISM2 Coupled Climate–Ice Sheet Model" (Sommers et al., 2021, Paleoceanography and Paleoclimatology)
54 ENVIRONMENTAL SCIENCES↗
DECOVALEX-2023: Task C Final Report
The Full-scale Emplacement (FE) heater experiment at the Mont Terri Underground Rock Laboratory (URL) was designed and conducted by Nagra to replicate an emplacement tunnel of Nagra’s reference repository design at 1:1 scale. Alongside testing the technical feasibility of constructing disposal tunnels, emplacing waste containers in the tunnels and then backfilling them, the main goals of the FE experiment are (1) to obtain a better understanding of the coupled effects of induced thermo-hydro-mechanical (THM) processes that may occur and (2) to validate existing coupled THM models (Müller et al., 2017). A key aspect of ensuring safety for repositories located in low-permeability rock involves minimizing any damage to the rock itself, thereby preserving its integrity and promoting a stable environment Amongst a number of processes that could damage the rock is the increase in pore pressure due to thermal loading caused by heat emitted from the waste. To reduce the potential damage of the rock, it is important to analyse the evolution of heat over time due to the heat load of the containers and assess possible consequences by coupled THM models. The aim of Task C of DECOVALEX-2023 was to build 3D numerical models of the FE experiment, focussing in particular on the heating induced pore pressure change in the Opalinus Clay. Data from a large number of sensors were available from the FE experiment for model comparison. These sensors measured temperature and relative humidity in the bentonite around the heaters, and temperature, pressure and displacement/strain in the surrounding Opalinus clay. Data were available from the start of excavation (April 2012) up to August 2020 for most sensors (more than 5 years from the start of heating in December 2014). To fulfil the overall aim of the task, the work was broken down into a number of steps, starting with simpler models to build confidence in each team’s approach and then moving to more complex models that better represent the FE experiment. Step 0 consisted of 2D benchmark models, gradually increasing the number of processes that are represented from thermal (T) only models in Step 0a, to coupled thermal hydraulic (TH) models in Step 0b with a representation of changing porosity, to coupled thermo-hydro-mechanical (THM) models in Step 0c, where porosity changes are calculated by the mechanical model. A detailed specification of processes, parameters, initial and boundary conditions was provided for this step, with the ambition that all teams would work towards close agreement in their model results, thus building confidence in the model implementations. vi It was not straightforward to achieve agreement between the teams, so additional steps (Step 0b2, 0b3, 0c2, 0c3) were added along with derivation of some analytical solutions against which the models could be compared. The reasons for the differences between teams were investigated and found to be caused primarily by different conceptual model assumptions (including temperature dependence of the thermal expansion of water), different model formulations (including porosity evolution) and differences in modelled domain sizes, boundary conditions and grid discretisation. This demonstrates that comparisons between multiple modelling teams and/or comparison with analytical results and experimental data are highly beneficial in providing an indication of uncertainty in model predictions. At the conclusion of Step 0, almost all teams had achieved a close agreement in model results and those that had not achieved an agreement knew the reason for this. Step 1 moved from 2D models to 3D models of the FE experiment without adding technical features like shotcrete or EDZ, and only considering the heating phase. Initially the 3D model was tightly specified to continue to build confidence in the model implementations (Step 1a). The results of Step 1a were compared to the data from the FE-experiment without the teams seeing the data. The teams were then provided with a sub-set of the data from the FE-experiment and invited to consider how best to use the large dataset for model comparison (Step 1b). Teams were then asked to use the data provided to calibrate their models, only changing material property values rather than adding features or processes to their models (Step 1c). In Step 1, teams were asked to only model the heating phase of the experiment, so pressure in the Opalinus Clay was reported as change in pressure since the initial conditions were specified rather than modelled. The change from 2D to 3D models was accompanied by an increase in the dispersion of results between the teams. Some of this was resolved during the task, but some remained and is potentially due to model discretisation. Calibration of parameters was useful in improving the fit of the models to the data but the remaining differences indicated that the models were missing features or processes. In Step 2, the teams were asked to update their models with additional features and processes as well as calibrating parameters to try and improve the fit of the models to the data. Teams were encouraged to represent ventilation of the open FE tunnel prior to backfilling with heaters and bentonite and in Step 2, the absolute pressure in the Opalinus Clay was compared between the teams. Teams took different approaches, but there was consideration of adding shotcrete and an EDZ into the model, representing stress change during excavation and different approaches to modelling ventilation of the FE tunnel. Overall, the documented results showed a very good agreement for temperature. The results for porewater pressure evolution showed a significant improvement for most teams compared to Step 1c with a good agreement to the measurements for several teams whereas some teams overpredicted the pressure increase and others overpredicted the drainage effect especially for the sensors close to the heater. Step 3 was an opportunity for teams to use the models developed in Step 1 and Step 2 to make predictions about the temperature and pressure changes that will be expected at the FE experiment over the next few years in light of the planned changes in thermal output of the heaters.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Modeling Land-Atmosphere Coupling at Cloud-Resolving Scale Within the Multiple Atmosphere Multiple Land (MAML) Framework in SP-E3SM
Representing subgrid variabilities of land surface processes and their upscaled effects is crucial for global climate modeling. Here, we implement a Multiple Atmosphere Multiple Land (MAML) framework in the superparamaterized version of E3SM (SP32 E3SM) to explicitly simulate the subgrid variabilities of land states and fluxes at cloud-resolving scale and their interactions with atmosphere. Comparing to the standard SP-E3SM in which all the atmospheric columns of the cloud resolving model embedded within the global atmospheric model grid interact with the same land surface (i.e., Multiple Atmosphere Single Land (MASL)), the impact of MAML on the strength of land-atmosphere coupling is limited, partly because the current implementation mainly facilitates one-way coupling between the cloud-resolving model and the land surface model. Despite such limitation, MAML increases the surface latent heat flux at the expense of sensible heat flux, and increases precipitation in India, Amazon, and Central Africa, reducing the model dry bias compared to the standard SP-E3SM. By employing a normalized gross moist stability (NGMS) diagnostic framework, we find that the increase in precipitation minus evaporation (PE) is primarily driven by the change in large-scale moisture convergence, particularly by the increase of water vapor in the lower atmosphere, while the local effect of total surface energy flux plays a minor role in the P-E change. More specifically, MAML changes the surface energy partitioning (evaporative fraction), increases the atmosphere water vapor, and further increases P-E by decreasing the NGMS. Finally, future development in the MAML framework is discussed.
54 ENVIRONMENTAL SCIENCES↗
Uncertainty Quantification in High-Low Dynamic System Coupling using RAVEN and TRANSFORM
This work demonstrates new functionality and applications stemming from the development of high-fidelity to low-fidelity (high-low) coupling for system simulations and to further explore the capabilities of the Risk Analysis Virtual Environment (RAVEN) in the performance of uncertainty quantification in this kind of high-low coupled system models. The work builds from previous work on high-low coupling that utilized COBRA-TF (CTF), the high-fidelity subchannel analysis code, with a low fidelity model built in ORNL’s TRANSFORM, the system analysis code, utilizing the Functional Mock-Up Interface (FMI). Steady-state and transient analysis examples using the high/low coupled models generated from CTF and TRANSFORM/FMI are investigated. The workflows for both steady-state and transient coupled simulations are described. A steady-state parameter sweep and uncertainty analysis of the primary flow rates and reactor power are demonstrated. Likewise, a transient pump trip and power ramp sensitivity studies are also demonstrated. This work elucidates some of the potential benefits and future needs of using RAVEN for high/low system coupling analysis of energy systems. It also shows some of the difficulties that can be encountered in coupling system simulations.
Williams, Wesley↗
Precision Higgs couplings in neutral naturalness models: an effective field theory approach
The Higgs sector in neutral naturalness models provides a portal to the hidden sectors, and thus measurements of Higgs couplings at current and future colliders play a central role in constraining the parameter space of the model. We investigate a class of neutral naturalness models, in which the Higgs boson is a pseudo-Goldstone boson from the universal SO(N)/SO(N - 1) coset structure. Integrating out the radial mode from the spontaneous global symmetry breaking, we obtain various dimension-six operators in the Standard Model effective field theory, and calculate the low energy Higgs effective potential with radiative corrections included. We perform a χ 2 fit to the Higgs coupling precision measurements at current and future colliders and show that the new physics scale could be explored up to 2.3 (2.4) TeV without (with) the Higgs invisible decay channels at future Higgs factories. The limits are comparable to the indirect constraints obtained via electroweak precision measurements.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Prediction of softening kinetics and recrystallization texture in non-isothermally annealed bulged tubes using CPFEM and CA models
A hierarchically coupled cellular automata (CA) model, crystal plasticity finite element method (CPFEM), and thermal finite element (FE) model is developed to predict the softening kinetics of the bulged steel tube during non-isothermal annealing. Through the developed model, the kinetics of softening mechanisms including static recovery (SRV) and static recrystallization (SRX), as well as the recrystallization texture are predicted. Later, the Johnson-Mehl-Avrami-Kohnogorov (JMAK) model based on the predicted SRX data is developed to interpret the recrystallization behavior of the material. To perform this study, diverse experimental tests including tube hydroforming (THF), annealing, uniaxial tensile test, hardness test, as well as microstructure observations through optical microscopy and Electron Backscatter Diffraction (EBSD) tests on steel tube are performed. The obtained experimental data are utilized to calibrate and verify the implemented CPFEM model for simulation of THF process, thermal FE model for prediction of the local temperature over annealing time, and CA algorithm for modeling of the softening kinetics and texture evolution throughout the annealing process. The study shows that the predicted deformation characteristics, softening kinetics, recrystallization texture and temperature profile during non-isothermal annealing are in good agreement with experimental data. During the annealing process, a total of four stages for the kinetics of softening mechanisms is observed: No softening; SRV only; SRV dominant; and SRX dominant. During the progress of SRX, the behavior of recrystallization is abruptly changed, confirming that two different mechanisms are controlling the kinetics of transformation.
36 MATERIALS SCIENCE↗
A 1D coupled transport/cold plasma wave model for parallel ponderomotive density modification near RF actuators
The plasma density profile close to RF actuators plays a vital role in the efficiency of power coupling from the antenna to the core plasma of fusion devices. Nonlinear interactions can be significant when considering the interaction of applied RF power and the tenuous SOL plasma, and the inclusion of this region in computational analysis is essential to providing robust predictive simulations of plasma response to RF application. Of particular interest is the observed density modification during RF operation. Here we present a self consistent description of the nonlinear ponderomotive force in the antenna near field region, via coupling a 1D field aligned transport model with a 1D cold plasma wave model. Results are presented, which show a parallel electric field of 300 Vm -1 drove a density reduction of ~ 20% over 25 RF periods.
Barnett, Rhea↗
Modeling of thermal pressurization in tight claystone using sequential THM coupling: Benchmarking and validation against in-situ heating experiments in COx claystone
We apply thermoporoelasticity and a sequentially coupling technique for modeling thermally-driven coupled Thermo-Hydro-Mechanical (THM) processes in tight claystone. A THM benchmark case with a corresponding analytic solution for thermoporoelasticity under a constant heat loading verifies the model. Thereafter, two in situ heating experiments are simulated for model validation: a smaller-scale heating experiment (TED experiment) and a larger-scale experiment (ALC experiment) in Callovo-Oxfordian (COx) claystone at the Meuse/Haute-Marne underground research laboratory in France. The model exhibits good performance to match the observed temperature and pore pressure evolution for the smaller-scale TED experiment. For the larger-scale ALC experiment, general trends of thermal-pressurization are captured in the modeling, but pressure is underestimated at some monitoring points during cool-down. This indicates that the THM response in the field may be affected by the variability of rock's properties or irreversible or time-dependent mechanical processes that are not included in the current thermoporoelastic model. The main contributions of this work are as follows: (1) we verify and validate the numerical simulator, TOUGH-FLAC, to be a valuable coupled THM modeling tool; (2) prove that the laboratory determined material parameters can be used as reference values for upscaling experiments. However, to better identify and quantify THM processes with modeling of in situ tests, more emphasize should be dedicated to obtaining high-quality mechanical deformation data.
58 GEOSCIENCES↗
Wildfire towers drive firebrand lofting: insights from coupled fire-atmosphere model simulations
Wildfire behavior is shaped by complex fire dynamics, with firebrands playing a critical role in spot fire ignition and fire spread. While previous studies have explored firebrand generation and transport, the specific role of towers and troughs from wildland fires in the lofting of firebrands remains unquantified. This study addresses that gap by using physics-based coupled fire-atmosphere model simulations to examine how wildfire towers (updrafts) and troughs (downdrafts) influence firebrand lofting. Our results show that the majority of firebrands (78.85%) are lofted from towers, where strong updrafts drive long-range transport. In contrast, only 21.15% of firebrands are lofted within troughs, where downdrafts cause most firebrands to fall near the fireline. We also find that firebrand size significantly influences lofting behavior, with smaller particles (1 mm radius) exhibiting the strongest correlation with updraft intensity. These findings highlight the dominant role of wildfire towers in promoting long-distance firebrand dispersal—an essential factor in rapid wildfire growth and wildland-urban interface (WUI) fire risks. By quantifying the relationship between firebrand lofting and fire-induced atmospheric features, this study provides critical insights to improve spot fire modeling, support mitigation planning, and enhance firefighter and WUI community safety in spot fire-prone regions.
54 ENVIRONMENTAL SCIENCES↗
Prediction of the thermal-hydraulic-mechanical response of a geological repository at large scale and sensitivity analyses
To better understand the Thermal-Hydraulic-Mechanical (THM) response of the rock caused by the high-level radioactive waste (HLW) released heat in the Callovo-Oxfordian formation (COx) in the near-field and far-field areas, a series of coupled THM modelling is performed. This study presents a case study based on the French geological repository for HLW (Cigéo project) to better assess modelling of deep geological repositories (DGR) within the DECOVALEX-2019 framework. In this study, the proposed coupled THM model, implemented in the platform of COMSOL, is validated by comparing the modelled THM results for a point heat source in an infinite rock mass with an analytical solution. It is also validated by comparing its thermal component for the Base Case with accurate results calculated using another numerical method. To shorten the calculation time, the DGR is modelled using six cuboid blocks representing 162 HLW cells and only six central HLW cells of interest modelled in detail. This simplification is validated. This model can easily be used to model not only the near-field THM response of the rock but also the far-field THM response. The influence of the boundary conditions applied on the gallery wall, HLW cell walls, and external surfaces of the model and the influence of different vertical dimensions of the geometry are studied. The sensitivity analyses of the THM parameters of COx on the THM response at different locations are performed. Finally, the influence of using a 2-dimensional (2D) model to represent a 3-dimensional (3D) repository is also investigated.
58 GEOSCIENCES↗