SEARCH · Search NASA

Results for “coupling coefficient”

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,135 records · Page 63

All sources NASA NTRS DOE OSTI

DOE OSTI2018

An adaptive least-squares global sensitivity method and application to a plasma-coupled combustion prediction with parametric correlation

We introduce an efficient non-intrusive surrogate-based methodology for global sensitivity analysis and uncertainty quantification. Modified covariance-based sensitivity indices (mCov-SI) are defined for outputs that reflect correlated effects. The overall approach is applied to simulations of a complex plasma-coupled combustion system with disparate uncertain parameters in sub-models for chemical kinetics and a laser-induced breakdown ignition seed. The surrogate is based on an Analysis of Variance (ANOVA) expansion, such as widely used in statistics, with orthogonal polynomials representing the ANOVA subspaces and a polynomial dimensional decomposition (PDD) representing its multi-dimensional components. The coefficients of the PDD expansion are obtained using a least-squares regression, which both avoids the direct computation of high-dimensional integrals and affords an attractive flexibility in choosing sampling points. This facilitates importance sampling using a Bayesian calibrated posterior distribution, which is fast and thus particularly advantageous in common practical cases, such as our large-scale demonstration, for which the asymptotic convergence properties of polynomial expansions cannot be realized due to computation expense. Effort, instead, is focused on efficient finite-resolution sampling. Standard covariance-based sensitivity indices (Cov-SI) are employed to account for correlation of the uncertain parameters. Magnitude of Cov-SI is unfortunately unbounded, which can produce extremely large indices that limit their utility. Alternatively, mCov-SI are then proposed in order to bound this magnitude $\in$ [0, 1]. The polynomial expansion is coupled with an adaptive ANOVA strategy to provide an accurate surrogate as the union of several low-dimensional spaces, avoiding the typical computational cost of a high-dimensional expansion. It is also adaptively simplified according to the relative contribution of the different polynomials to the total variance. Finally, the approach is demonstrated for a laser-induced turbulent combustion simulation model, which includes parameters with correlated effects.

42 ENGINEERING↗

DOE OSTI2021

Enhanced relaxed physical factorization preconditioner for coupled poromechanics

The relaxed physical factorization (RPF) preconditioner is a recent algorithm allowing for the efficient and robust solution to the block linear systems arising from the three-field displacement-velocity-pressure formulation of coupled poromechanics. For its application, however, it is necessary to invert blocks with the algebraic form C^ = (C + βFF T ), where C is a symmetric positive definite matrix, FF T a rank-deficient term, and β a real non-negative coefficient. The inversion of C^, performed in an inexact way, can become unstable for large values of β, as it usually occurs at some stages of a full poromechanical simulation. In this work, we propose a family of algebraic techniques to stabilize the inexact solve with C^. This strategy can prove useful in other problems as well where such an issue might arise, such as augmented Lagrangian preconditioning techniques for Navier-Stokes or incompressible elasticity. First, we introduce an iterative scheme obtained by a natural splitting of matrix C^. Second, we develop a technique based on the use of a proper projection operator annihilating the near-kernel modes of C^. Both approaches give rise to a novel class of preconditioners denoted as Enhanced RPF (ERPF). Furthermore, effectiveness and robustness of the proposed algorithms are demonstrated in both theoretical benchmarks and real-world large-size applications, outperforming the native RPF preconditioner.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Investigating Electron Conductivity Regimes in the Bacterial Cytochrome Wire OmcS

The anaerobic bacterium Geobacter sulfurreducens produces extracellular, electronically conductive cytochrome polymer wires that are conductive over micron length scales. Structure models from cryo-electron microscopy data show OmcS wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. However, the mechanism by which this heme arrangement supports long-range electronic conduction remains unknown. Structure models from cryo-electron microscopy data show these wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. Existing computational models using static heme redox potentials and coupling energies fail to explain experimental observations, predicting conductances 10,000 to 100,000 times lower than measured values. Here, we investigate how dynamic disorder affects site energies, interheme coupling, and long-range electronic conductivity within these cytochrome wires. We introduce an approach to extract charge carrier site information directly from Kohn–Sham density functional theory, without employing projector schemes, and show that site and coupling energies are highly sensitive to changes in interheme geometry and the surrounding electrostatic environment. Unlike models that incorporate dynamic disorder as a thermally averaged quantity, our quantum charge carrier model incorporates proxies for dynamic disorder through decoherence corrections, yielding predicted diffusion coefficient closer to what is expected from experiment and comparable with other organic-based electronic materials. Based on these simulations, we propose that the instantaneous fluctuations of the local electrostatic environment can transiently lift energy degeneracies and delocalize charge carriers. Furthermore, these studies reveal how incorporating dynamic fluctuations associated with the environment resolves the discrepancy between theory and experiment in microbial cytochrome wires and highlight design principles for bioinspired, heme-based conductive materials.

Bioinorganic chemistry↗

DOE OSTI2017

Transport modeling of convection dominated helicon discharges in Proto-MPEX with the B2.5-Eirene code

Data-constrained interpretative analyses of plasma transport in convection dominated helicon discharges in the Proto-MPEX linear device, and predictive calculations with additional Electron Cyclotron Heating/Electron Bernstein Wave (ECH/EBW) heating, are reported. The B2.5-Eirene code, in which the multi-fluid plasma code B2.5 is coupled to the kinetic Monte Carlo neutrals code Eirene, is used to fit double Langmuir probe measurements and fast camera data in front of a stainless-steel target. The absorbed helicon and ECH power (11 kW) and spatially constant anomalous transport coefficients that are deduced from fitting of the probe and optical data are additionally used for predictive simulations of complete axial distributions of the densities, temperatures, plasma flow velocities, particle and energy fluxes, and possible effects of alternate fueling and pumping scenarios. The somewhat hollow electron density and temperature radial profiles from the probe data suggest that Trivelpiece-Gould wave absorption is the dominant helicon electron heating source in the discharges analyzed here. There is no external ion heating, but the corresponding calculated ion temperature radial profile is not hollow. Rather it reflects ion heating by the electron-ion equilibration terms in the energy balance equations and ion radial transport resulting from the hollow density profile. With the absorbed power and the transport model deduced from fitting the sheath limited discharge data, calculated conduction limited higher recycling conditions were produced by reducing the pumping and increasing the gas fueling rate, resulting in an approximate doubling of the target ion flux and reduction of the target heat flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Transient Modeling and Simulation of a Generic Stable Salt Reactor

A SAM system-level model of a generic stable salt reactor has been developed to investigate thermal-hydraulic behavior and safety performance under steady and transient conditions. The model integrates information generated from a reactor physics analysis using PROTEUS and PERSENT, and a computation fluid dynamics (CFD) analysis using STAR-CCM+. A loose, iterative coupling scheme between PROTEUS and SAM is implemented to calculate the equilibrium power and temperature distributions in the steady-state critical core condition. The converged steady-state model is then used in PERSENT to calculate the four reactivity feedback temperature coefficients (Doppler, fuel density, coolant density, and core radial expansion) and kinetic parameters that are needed in SAM to model the temperature feedback effects in transient simulations. Within the fully enclosed liquid fuel pins, natural convection is the dominant heat transfer mechanism. The STAR-CCM+ model of the fuel pin considers conjugate heat transfer from the liquid fuel salt to the pin cladding and external reactor coolant. The CFD results of the axial and radial temperature profiles are used to empirically determine an effective fuel salt thermal conductivity in the SAM fuel pin model so that the temperatures predicted by the SAM model match as closely as possible the CFD results. In the central region of the fuel pin, the effective thermal conductivity is as high as similar to 60 times the physical fuel salt thermal conductivity. The whole-plant SAM model is then used to simulate an unprotected station blackout transient. The results of this simulation showed that the large negative fuel axial expansion reactivity feedback reduces fission power to similar to 2.4% nominal power. The core is cooled by natural circulation, which removes heat in the core to the emergency heat removal system, and ultimately, to the ambient. However, peak fuel salt and cladding temperatures can potentially reach as high as 1500 K, albeit briefly, if the shutdown mechanism fails to operate.

stable salt reactor; transient simulations; system↗

NASA NTRS1990

An annular gas seal analysis using empirical entrance and exit region friction factors

Wall shear stress results from stationary-rotor flow tests of five annular gas seals are used to develop entrance and exit region friction factor models. The friction factor models are used in a bulk-flow seal analysis which predicts leakage and rotor-dynamic coefficients. The predictions of the analysis are compared to experimental results and to the predictions of Nelson's analysis (1985). The comparisons are for smooth-rotor seals with smooth and honeycomb-stators. The present analysis predicts the destabilizing cross-coupled stiffness of a seal better than Nelson's analysis. Both analyses predict direct damping well and direct stiffness poorly.

Elrod, D. A.↗

NASA NTRS1996

High Temperature Degradation Mechanisms in Polymer Matrix Composites

Polymer matrix composites are increasingly used in demanding structural applications in which they may be exposed to harsh environments. The durability of such materials is a major concern, potentially limiting both the integrity of the structures and their useful lifetimes. The goal of the current investigation is to develop a mechanism-based model of the chemical degradation which occurs, such that given the external chemical environment and temperatures throughout the laminate, laminate geometry, and ply and/or constituent material properties, we can calculate the concentration of diffusing substances and extent of chemical degradation as functions of time and position throughout the laminate. This objective is met through the development and use of analytical models, coupled to an analysis-driven experimental program which offers both quantitative and qualitative information on the degradation mechanism. Preliminary analyses using a coupled diffusion/reaction model are used to gain insight into the physics of the degradation mechanisms and to identify crucial material parameters. An experimental program is defined based on the results of the preliminary analysis which allows the determination of the necessary material coefficients. Thermogravimetric analyses are carried out in nitrogen, air, and oxygen to provide quantitative information on thermal and oxidative reactions. Powdered samples are used to eliminate diffusion effects. Tests in both inert and oxidative environments allow the separation of thermal and oxidative contributions to specimen mass loss. The concentration dependency of the oxidative reactions is determined from the tests in pure oxygen. Short term isothermal tests at different temperatures are carried out on neat resin and unidirectional macroscopic specimens to identify diffusion effects. Mass loss, specimen shrinkage, the formation of degraded surface layers and surface cracking are recorded as functions of exposure time. Geometry effects in the neat resin, and anisotropic diffusion effects in the composites, are identified through the use of specimens with different aspect ratios. The data is used with the model to determine reaction coefficients and effective diffusion coefficients. The empirical and analytical correlations confirm the preliminary model results which suggest that mass loss at lower temperatures is dominated by oxidative reactions and that these reaction are limited by diffusion of oxygen from the surface. The mechanism-based model is able to successfully capture the basic physics of the degradation phenomena under a wide range of test conditions. The analysis-based test design is successful in separating out oxidative, thermal, and diffusion effects to allow the determination of material coefficients. This success confirms the basic picture of the process; however, a more complete understanding of some aspects of the physics are required before truly predictive capability can be achieved.

Cunningham, Ronan A.↗

NASA NTRS1997

High Temperature Degradation Mechanisms in Polymer Matrix Composites

Polymer matrix composites are increasingly used in demanding structural applications in which they may be exposed to harsh environments. The durability of such materials is a major concern, potentially limiting both the integrity of the structures and their useful lifetimes. The goal of the current investigation is to develop a mechanism-based model of the chemical degradation which occurs, such that given the external chemical environment and temperatures throughout the laminate, laminate geometry, and ply and/or constituent material properties, we can calculate the concentration of diffusing substances and extent of chemical degradation as functions of time and position throughout the laminate. This objective is met through the development and use of analytical models, coupled to an analysis-driven experimental program which offers both quantitative and qualitative information on the degradation mechanism. Preliminary analyses using coupled diffusion/reaction model are used to gain insight into the physics of the degradation mechanisms and to identify crucial material parameters. An experimental program is defined based on the results of the preliminary analysis which allows the determination of the necessary material coefficients. Thermogravimetric analyses are carried out in nitrogen, air, and oxygen to provide quantitative information on thermal and oxidative reactions. Powdered samples are used to eliminate diffusion effects. Tests in both inert and oxidative environments allow the separation of thermal and oxidative contributions to specimen mass loss. The concentration dependency of the oxidative reactions is determined from the tests in pure oxygen. Short term isothermal tests at different temperatures are carried out on neat resin and unidirectional macroscopic specimens to identify diffusion effects. Mass loss, specimen shrinkage, the formation of degraded surface layers and surface cracking are recorded as functions of exposure time. Geometry effects in the neat resin, and anisotropic diffusion effects in the composites, are identified through the use of specimens with different aspect ratios. The data is used with the model to determine reaction coefficients and effective diffusion coefficients. The empirical and analytical correlations confirm the preliminary model results which suggest that mass loss at lower temperatures is dominated by oxidative reactions and that these reaction are limited by diffusion of oxygen from the surface. The mechanism-based model is able to successfully capture the basic physics of the degradation phenomena under a wide range of test conditions. The analysis-based test design is successful in separating out oxidative, thermal, and diffusion effects to allow the determination of material coefficients. This success confirms the basic picture of the process; however, a more complete understanding of some aspects of the physics are required before truly predictive capability can be achieved.

Cunningham, Ronan A.↗

NASA NTRS2001

Determining CDOM Absorption Spectra in Diverse Aquatic Environments Using a Multiple Pathlength, Liquid Core Waveguide System

We evaluated the accuracy, sensitivity and precision of a multiple pathlength, liquid core waveguide (MPLCW) system for measuring colored dissolved organic matter (CDOM) absorption in the UV-visible spectral range (370-700 nm). The MPLCW has four optical paths (2.0, 9.8, 49.3, and 204 cm) coupled to a single Teflon AF sample cell. Water samples were obtained from inland, coastal and ocean waters ranging in salinity from 0 to 36 PSU. Reference solutions for the MPLCW were made having a refractive index of the sample. CDOM absorption coefficients, aCDOM, and the slope of the log-linearized absorption spectra, S, were compared with values obtained using a dual-beam spectrophotometer. Absorption of phenol red secondary standards measured by the MPLCW at 558 nm were highly correlated with spectrophotometer values and showed a linear response across all four pathlengths. Values of aCDOM measured using the MPLCW were virtually identical to spectrophotometer values over a wide range of concentrations. The dynamic range of aCDOM for MPLCW measurements was 0.002 - 231.5 m-1. At low CDOM concentrations spectrophotometric aCDOM were slightly greater than MPLCW values and showed larger fluctuations at longer wavelengths due to limitations in instrument precision. In contrast, MPLCW spectra followed an exponential to 600 nm for all samples.

Miller, Richard L.↗

NASA NTRS2000

Determining CDOM Absorption Spectra in Diverse Coastal Environments Using a Multiple Pathlength, Liquid Core Waveguide System. Measuring the Absorption of CDOM in the Field Using a Multiple Pathlength Liquid Waveguide System

We evaluated the accuracy, sensitivity and precision of a multiple pathlength, liquid core waveguide (MPLCW) system for measuring colored dissolved organic matter (CDOM) absorption in the UV-visible spectral range (370-700 nm). The MPLCW has four optical paths (2.0, 9.8, 49.3, and 204 cm) coupled to a single Teflon AF sample cell. Water samples were obtained from inland, coastal and ocean waters ranging in salinity from 0 to 36 PSU. Reference solutions for the MPLCW were made having a refractive index of the sample. CDOM absorption coefficients, a(sub CDOM), and the slope of the log-linearized absorption spectra, S, were compared with values obtained using a dual-beam spectrophotometer. Absorption of phenol red secondary standards measured by the MPLCW at 558 nm were highly correlated with spectrophotometer values (r > 0.99) and showed a linear response across all four pathlengths. Values of a(sub CDOM) measured using the MPLCW were virtually identical to spectrophotometer values over a wide range of concentrations. The dynamic range of a(sub CDOM) for MPLCW measurements was 0.002 - 231.5/m. At low CDOM concentrations (a(sub 370) < 0.1/m) spectrophotometric a(sub CDOM) were slightly greater than MPLCW values and showed larger fluctuations at longer wavelengths due to limitations in instrument precision. In contrast, MPLCW spectra followed an exponential to 600 nm for all samples. The maximum deviation in replicate MPLCW spectra was less than 0.001 absorbance units. The portability, sampling, and optical characteristics of a MPLCW system provide significant enhancements for routine CDOM absorption measurements in a broad range of natural waters.

Miller, Richard L.↗

NASA NTRS2006

Using CO2:CO Correlations to Improve Inverse Analyses of Carbon Fluxes

Observed correlations between atmospheric concentrations of CO2 and CO represent potentially powerful information for improving CO2 surface flux estimates through coupled CO2-CO inverse analyses. We explore the value of these correlations in improving estimates of regional CO2 fluxes in east Asia by using aircraft observations of CO2 and CO from the TRACE-P campaign over the NW Pacific in March 2001. Our inverse model uses regional CO2 and CO surface fluxes as the state vector, separating biospheric and combustion contributions to CO2. CO2-CO error correlation coefficients are included in the inversion as off-diagonal entries in the a priori and observation error covariance matrices. We derive error correlations in a priori combustion source estimates of CO2 and CO by propagating error estimates of fuel consumption rates and emission factors. However, we find that these correlations are weak because CO source uncertainties are mostly determined by emission factors. Observed correlations between atmospheric CO2 and CO concentrations imply corresponding error correlations in the chemical transport model used as the forward model for the inversion. These error correlations in excess of 0.7, as derived from the TRACE-P data, enable a coupled CO2-CO inversion to achieve significant improvement over a CO2-only inversion for quantifying regional fluxes of CO2.

Palmer, Paul I.↗

NASA NTRSDate not supplied

Aeroelastic Analysis of Mach 0.8 Transonic Truss-Braced Wing Aircraft

This paper presents an aeroelastic analysis of the Mach 0.8 Transonic Truss-Braced Wing(TTBW) aircraft jig shape using an in-house developed tool based on VSPAERO. A vortex-lattice model of the Mach 0.8 TTBW model is developed, and a transonic and viscous flow correction method is implemented to account for transonic and viscous flow effects. A correction method for the wing-strut interference aerodynamics is developed and applied to the VSPAERO solver. The Galerkin method is used to calculate the geometry deformation under aerodynamic force. The aero-structural analysis solver VSPAERO coupled to the mode shapes computed by NASTRAN using the Galerkin method provides a rapid aircraft aero-structual analysis. A high-fidelity CFD solver FUN3D is used to verify the results. The aeroelastic simulation results show that the aeroelastic lift coefficient is reduced about0.05∼0.07, drag polar is not affected by aeroelasticity, and the pitching moment is reduced about 30% at Mach0.8 and altitude 40,000 ft.

TTBW Aeroelastic↗

NASA NTRS2008

The Influence of Ablation on Radiative Heating for Earth Entry

Using the coupled ablation and radiation capability recently included in the LAURA flowfield solver, this paper investigates the influence of ablation on the shock-layer radiative heating for Earth entry. The extension of the HARA radiation model, which provides the radiation predictions in LAURA, to treat a gas consisting of the elements C, H, O, and N is discussed. It is shown that the absorption coefficient of air is increased with the introduction of the C and H elements. A simplified shock layer model is studied to show the impact of temperature, as well as the abundance of C and H, on the net absorption or emission from an ablation contaminated boundary layer. It is found that the ablation species reduce the radiative flux in the vacuum ultraviolet, through increased absorption, for all temperatures. However, in the infrared region of the spectrum, the ablation species increase the radiative flux, through strong emission, for temperatures above 3,000 K. Thus, depending on the temperature and abundance of ablation species, the contaminated boundary layer may either provide a net increase or decrease in the radiative flux reaching the wall. To assess the validity of the coupled ablation and radiation LAURA analysis, a previously analyzed Mars-return case (15.24 km/s), which contains significant ablation and radiation coupling, is studied. Exceptional agreement with previous viscous shock-layer results is obtained. A 40% decrease in the radiative flux is predicted for ablation rates equal to 20% of the free-stream mass flux. The Apollo 4 peak-heating case (10.24 km/s) is also studied. For ablation rates up to 3.4% of the free-stream mass flux, the radiative heating is reduced by up to 19%, while the convective heating is reduced by up to 87%. Good agreement with the Apollo 4 radiometer data is obtained by considering absorption in the radiometer cavity. For both the Mars return and the Apollo 4 cases, coupled radiation alone is found to reduce the radiative heating by 30 60% and the convective heating by less than 5%.

Johnston, Christopher O.↗

NASA NTRS1975

Effect of vertical-tail location on the aerodynamic characteristics at subsonic speeds of a close-coupled canard configuration

The effects were studied of various vertical-tail configurations on the longitudinal and lateral directional-stability characteristics of a general research fighter model utilizing wing-body-canard. The study indicates that the addition of the high canard resulted in an increase in total lift at angles of attack above 4 deg with a maximum lift coefficient about twice as large as that for the wing-body configuration. For the wing-body (canard off) configuration, the center-line vertical tail indicates positive vertical-tail effectiveness throughout the test angle-of-attack range; however, for this configuration none of the wing-mounted vertical-tail locations tested resulted in a positive directional-stability increment at the higher angles of attack. For the wing-body-canard configuration several outboard locations of the wing-mounted vertical tails were found.

Huffman, J. K.↗

NASA NTRS2015

CFD Modeling of the Multipurpose Hydrogen Test Bed (MHTB) Self-Pressurization and Spray Bar Mixing Experiments in Normal Gravity: Effect of the Accommodation Coefficient on the Tank Pressure

A CFD model for simulating the self-pressurization of a large scale liquid hydrogen storage tank is utilized in this paper to model the MHTB self-pressurization experiment. The kinetics-based Schrage equation is used to account for the evaporative and condensi ng interfacial mass flows in this model. The effect of the accommodation coefficient for calculating the interfacial mass transfer rate on the tank pressure during tank selfpressurization is studied. The values of the accommodation coefficient which were considered in this study vary from 1.0e-3 to 1.0e-1 for the explicit VOF model and from 1.0e-4 to 1.0e-3 for the implicit VOF model. The ullage pressure evolutions are compared against experimental data. A CFD model for controlling pressure in cryogenic storage tanks by spraying cold liquid into the ullage is also presented. The Euler-Lagrange approach is utilized for tracking the spray droplets and for modeling the interaction between the droplets and the continuous phase (ullage). The spray model is coupled with the VOF model by performing particle tracking in the ullage, removing particles from the ullage when they reach the interface, and then adding their contributions to the liquid. Droplet-ullage heat and mass transfer are modeled. The flow, temperature, and interfacial mass flux, as well as droplets trajectories, size distribution and temperatures predicted by the model are presented. The ul lage pressure and vapor temperature evolutions are compared with experimental data obtained from the MHTB spray bar mixing experiment. The effect of the accommodation coefficient for calculating the interfacial and droplet mass transfer rates on the tank pressure during mixing of the vapor using spray is studied. The values used for the accommodation coefficient at the interface vary from 1.0e-5 to 1.0e-2. The droplet accommodation coefficient values vary from 2.0e-6 to 1.0e-4.

cryogenic fluid storage↗

DOE OSTI2023

Two-fluid physical modeling of superconducting resonators in the ARTEMIS framework

In this work, we implement a new London equation module for superconductivity in the GPU-enabled ARTEMIS framework, and couple it to a finite-difference time-domain solver for Maxwell's equations. We apply this two-fluid approach to model a superconducting coplanar waveguide (CPW) resonator. We validate our implementation by verifying that the theoretical skin depth and reflection coefficients can be obtained for several superconductive materials, with different London penetration depths, over a range of frequencies. Our convergence studies show that the algorithm is second-order accurate in both space and time, except at superconducting interfaces where the approach is spatially first-order. In our CPW simulations, we leverage the GPU scalability of our code to compare the two-fluid model to more traditional approaches that approximate superconducting behavior and demonstrate that superconducting physics can show comparable performance to the assumption of quasi-infinite conductivity as measured by the Q-factor.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Constraint energy minimizing generalized multiscale finite element method for multi-continuum Richards equations

In fluid flow simulation, the multi-continuum model is a useful strategy. When the heterogeneity and contrast of coefficients are high, the system becomes multiscale, and some kinds of reduced order methods are demanded. Combining these techniques with nonlinearity, we will consider in this paper a dual-continuum model which is generalized as a multi-continuum model for a coupled system of nonlinear Richards equations as unsaturated flows, in complex heterogeneous fractured porous media; and we will solve it by a novel multiscale approach utilizing the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM). In particular, such a nonlinear system will be discretized in time and then linearized by Picard iteration (whose global convergence is proved theoretically). Subsequently, we tackle the resulting linearized equations by the CEM-GMsFEM and obtain proper offline multiscale basis functions to span the multiscale space (which contains the pressure solution). More specifically, we first introduce two new sources of samples, and the GMsFEM is used over each coarse block to build local auxiliary multiscale basis functions via solving local spectral problems, that are crucial for detecting high-contrast channels. Second, per oversampled coarse region, local multiscale basis functions are created through the CEM as constrainedly minimizing an energy functional. Various numerical tests for our approach reveal that the error converges with the coarse-grid size and that only few oversampling layers as well as basis functions are needed.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Updated Thermofluid Performance of the Simplified Flat Variant of the HEMJ

Our group has recently developed and studied “finger”-type divertors that are a simplified version of the helium-cooled modular divertor with multiple jets (HEMJ) using coupled computational fluid dynamics and thermal stress simulations. Such a simplified geometry could reduce complexity and cost given the large number of fingers required to cover the total divertor target area. Previous experimental studies for this simplified flat design reported lower heat transfer coefficients and higher pressure drops than the HEMJ, contrary to numerical predictions. Subsequent measurements determined that the original test section had significant dimensional variations in the jet exit holes. A new test section was therefore manufactured and tested in the Georgia Tech (GT) helium loop. The experimental results presented here for this test section at maximum heat flux of 7.1 MW/m 2 are in good agreement with numerical predictions. Correlations developed from these experimental data are extrapolated to predict the maximum heat flux that can be accommodated by the flat design and the coolant pumping power requirements under prototypical conditions. Lastly, numerical simulations are used to estimate the sensitivity of the flat design to geometric variations typical of manufacturing tolerances and variations in the gap width.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗