Search NASA⌕ Search

SEARCH · Search NASA

Results for “generalized finite differences”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Redox levels in aqueous solution: Effect of van der Waals interactions and hybrid functionals

We investigate redox levels in aqueous solution using a combination of ab initio molecular dynamics (MD) simulations and thermodynamic integration methods. The molecular dynamics are performed with both the semilocal Perdew-Burke-Ernzerhof functional and a nonlocal functional (rVV10) accounting for van der Waals (vdW) interactions. The band edges are determined through three different schemes, namely, from the energy of the highest occupied and of the lowest unoccupied Kohn-Sham states, from total-energy differences, and from a linear extrapolation of the density of states. It is shown that the latter does not depend on the system size while the former two are subject to significant finite-size effects. For the redox levels, we provide a formulation in analogy to the definition of charge transition levels for defects in crystalline materials. We consider the H{sup +}/H{sub 2} level defining the standard hydrogen electrode, the OH{sup −}/OH{sup ∗} level corresponding to the oxidation of the hydroxyl ion, and the H{sub 2}O/OH{sup ∗} level for the dehydrogenation of water. In spite of the large structural modifications induced in liquid water, vdW interactions do not lead to any significant structural effect on the calculated band gap and band edges. The effect on the redox levels is also small since the solvation properties of ionic species are little affected by vdW interactions. Since the electronic properties are not significantly affected by the underlying structural properties, it is justified to perform hybrid functional calculations on the configurations of our MD simulations. The redox levels calculated as a function of the fraction α of Fock exchange are found to remain constant, reproducing a general behavior previously observed for charge transition levels of defects. Comparison with experimental values shows very good agreement. At variance, the band edges and the band gap evolve linearly with α. For α ≃ 0.40, we achieve a band gap, band-edge positions, and redox levels in overall good agreement with experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

Analysis of subgrid models using direct and large-eddy simulations of isotropic turbulence

Direct and large eddy simulations of forced and decaying isotropic turbulence have been performed using a pseudospectral and a finite-difference code. Subgrid models that include a one-equation subgrid kinetic energy model with and without a stochastic backscatter forcing term and a new scale similarity model have been analyzed in both Fourier space and physical space. The Fourier space analysis showed that the energy transfer across the cutoff wavenumber k(sub c) is dominated by local interaction. The correlation between the exact and the modeled (by a spectral eddy viscosity) nonlinear terms and the subgrid energy transfer in physical space was found to be quite low. In physical space, a similar correlation analysis was carried out using top hat filtering. Results show that the subgrid stress and the energy flux predicted by the subgrid models correlates very well with the exact data. The scale similarity model showed very high correlation for reasonable grid resolution. However, with decrease in grid resolution, the scale similarity model became more uncorrelated, when compared to the kinetic energy subgrid model. The subgrid models were then used for large-eddy simulations for a range of Reynolds number. It was determined that the dissipation was modeled poorly and that the correlation with the exact results was quite low for all the models. In general, for coarse grid resolution, the scale similarity model consistently showed very low correlation while the kinetic energy model showed a relatively higher correlation. These results suggest that to use the scale similarity model relatively fine grid resolution may be required, whereas, the kinetic energy model could be used even in coarse grid.

Menon, S.↗
NASA NTRS1998

A Study of Failure in Small Pressurized Cylindrical Shells Containing a Crack

The deformation in the vicinity of axial cracks in thin pressurized cylinders is examined using small experimental The deformation in the vicinity of axial cracks in thin pressurized cylinders is examined using small experimental models. The loading applied was either symmetric or unsymmetric about the crack plane, the latter being caused by structural constraints such as stringers. The objective was two fold - one, to provide the experimental results which will allow computer modeling techniques to be evaluated for deformations that are significantly different from that experienced by flat plates, and the other to examine the deformations and conditions associated with the onset of crack kinking which often precedes crack curving. The stresses which control crack growth in a cylindrical geometry depend on conditions introduced by the axial bulging, which is an integral part of this type of failure. For the symmetric geometry, both the hoop and radial strain just ahead off the crack, r = a, were measured and these results compared with those obtained from a variety of structural analysis codes, in particular STAGS [1], ABAQUS and ANSYS. In addition to these measurements, the pressures at the onset of stable and unstable crack growth were obtained and the corresponding crack deformations measured as the pressures were increased to failure. For the unsymmetric cases, measurements were taken of the crack kinking angle, and the displacements in the vicinity of the crack. In general, the strains ahead of the crack showed good agreement between the three computer codes and between the codes and the experiments. In the case of crack behavior, it was determined that modeling stable tearing with a crack-tip opening displacement fracture criterion could be successfully combined with the finite-element analysis techniques as used in structural analysis codes. The analytic results obtained in this study were very compatible with the experimental observations of crack growth. Measured crack kinking angles also showed good agreement with theories based on the maximum principle stress criterion.

Barwell, Craig A.↗
DOE OSTI2025

Impact of representative ground motion level on seismic PSA with the boundary between overestimation and underestimation

One commonly used approach in seismic probabilistic safety assessment (PSA) is the discrete method. This method follows the standard PSA framework and can be applied to various models, such as multi-unit models, while reducing computational costs using standard software. However, due to the inability to subdivide intervals infinitely, the discrete method approximates with a finite number of subintervals. In practice, different numbers of subintervals are applied, and the representative ground motion level is selected based on expert judgment. When employing a smaller number of subintervals, it is important to take caution to prevent underestimation. This study analyzes the impact of the representative ground motion level on seismic risk. It confirms that underestimation can occur with a small number of subintervals depending on the representative ground motion level. This study also proposes a method for determining the boundary of underestimation and overestimation. The method is demonstrated through examples, providing a mathematical foundation for selecting appropriate representative ground motion levels. By avoiding underestimation, this research helps prevent the oversight of significant risk contributors and enhances the understanding of seismic risk.

99 - GENERAL AND MISCELLANEOUS↗
NASA NTRS1995

The scheme machine: A case study in progress in design derivation at system levels

The Scheme Machine is one of several design projects of the Digital Design Derivation group at Indiana University. It differs from the other projects in its focus on issues of system design and its connection to surrounding research in programming language semantics, compiler construction, and programming methodology underway at Indiana and elsewhere. The genesis of the project dates to the early 1980's, when digital design derivation research branched from the surrounding research effort in programming languages. Both branches have continued to develop in parallel, with this particular project serving as a bridge. However, by 1990 there remained little real interaction between the branches and recently we have undertaken to reintegrate them. On the software side, researchers have refined a mathematically rigorous (but not mechanized) treatment starting with the fully abstract semantic definition of Scheme and resulting in an efficient implementation consisting of a compiler and virtual machine model, the latter typically realized with a general purpose microprocessor. The derivation includes a number of sophisticated factorizations and representations and is also deep example of the underlying engineering methodology. The hardware research has created a mechanized algebra supporting the tedious and massive transformations often seen at lower levels of design. This work has progressed to the point that large scale devices, such as processors, can be derived from first-order finite state machine specifications. This is roughly where the language oriented research stops; thus, together, the two efforts establish a thread from the highest levels of abstract specification to detailed digital implementation. The Scheme Machine project challenges hardware derivation research in several ways, although the individual components of the system are of a similar scale to those we have worked with before. The machine has a custom dual-ported memory to support garbage collection. It consists of four tightly coupled processes--processor, collector, allocator, memory--with a very non-trivial synchronization relationship. Finally, there are deep issues of representation for the run-time objects of a symbolic processing language. The research centers on verification through integrated formal reasoning systems, but is also involved with modeling and prototyping environments. Since the derivation algebra is basd on an executable modeling language, there is opportunity to incorporate design animation in the design process. We are looking for ways to move smoothly and incrementally from executable specifications into hardware realization. For example, we can run the garbage collector specification, a Scheme program, directly against the physical memory prototype, and similarly, the instruction processor model against the heap implementation.

Johnson, Steven D.↗
NASA NTRS2013

Spurious Behavior of Shock-Capturing Methods: Problems Containing Stiff Source Terms and Discontinuities

The goal of this paper is to relate numerical dissipations that are inherited in high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities. For pointwise evaluation of the source term, previous studies indicated that the phenomenon of wrong propagation speed of discontinuities is connected with the smearing of the discontinuity caused by the discretization of the advection term. The smearing introduces a nonequilibrium state into the calculation. Thus as soon as a nonequilibrium value is introduced in this manner, the source term turns on and immediately restores equilibrium, while at the same time shifting the discontinuity to a cell boundary. The present study is to show that the degree of wrong propagation speed of discontinuities is highly dependent on the accuracy of the numerical method. The manner in which the smearing of discontinuities is contained by the numerical method and the overall amount of numerical dissipation being employed play major roles. Moreover, employing finite time steps and grid spacings that are below the standard Courant-Friedrich-Levy (CFL) limit on shockcapturing methods for compressible Euler and Navier-Stokes equations containing stiff reacting source terms and discontinuities reveals surprising counter-intuitive results. Unlike non-reacting flows, for stiff reactions with discontinuities, employing a time step and grid spacing that are below the CFL limit (based on the homogeneous part or non-reacting part of the governing equations) does not guarantee a correct solution of the chosen governing equations. Instead, depending on the numerical method, time step and grid spacing, the numerical simulation may lead to (a) the correct solution (within the truncation error of the scheme), (b) a divergent solution, (c) a wrong propagation speed of discontinuities solution or (d) other spurious solutions that are solutions of the discretized counterparts but are not solutions of the governing equations. The present investigation for three very different stiff system cases confirms some of the findings of Lafon & Yee (1996) and LeVeque & Yee (1990) for a model scalar PDE. The findings might shed some light on the reported difficulties in numerical combustion and problems with stiff nonlinear (homogeneous) source terms and discontinuities in general.

Yee, Helen M. C.↗
DOE OSTI2019

A new proposal for a high-performance 4-channel demultiplexer based on 2D photonic crystal using three cascaded ring resonators for applications in advanced optical systems

In this paper, a high-performance 4-channel demultiplexer (DEMUX) based on 2D photonic crystal using three cascaded ring resonators is proposed. The proposed DEMUX is composed of square lattice of silicon dielectric rods with a refractive index of 3.46 suspended in air. For the proposed device, the lattice constant a = 590 nm and the filling factor r/a = 0.2. In this study, the plane wave expansion and finite-difference time-domain methods are used respectively to characterize the photonic band gap and to investigate the optical behaviour of the structure. It is demonstrated that by choosing appropriate values for the radius of the inner rods, scattering rods, coupling rods and the coupling length between the bus and the drop waveguides, the desired wavelengths can effectively be separated. The resonant wavelengths of the proposed DEMUX are 1528.7 nm, 1535.3 nm, 1531 nm, and 1533 nm. The mean quality factor, spectral bandwidth and the transmission efficiency of our device are computed to be respectively 7358.5, 0.217 nm and 99.25%. The proposed DEMUX has a mean channel spacing, minimum and maximum values of crosstalk of about 2.2 nm, − 46.68 dB and − 9.79 dB, respectively. Furthermore, the footprint of the device is 294.4 μm{sup 2}, which is suitable for photonic integrated circuits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Predicting effective thermal and elastic properties of cementitious composites containing polydispersed hollow and core-shell micro-particles

Micro- and nano-sized hollow and core-shell particles (CSPs) have attracted tremendous interests in developing functional cementitious composites and concretes. In this paper, a general method is developed to predict the thermal and elastic properties of cementitious composites containing core-shell and hollow micro-particles. The model follows a two-stage homogenization process – CSP inclusions together with their surrounding interfacial transition zone (ITZ) are first treated as equivalent solid particles, and then the homogenized properties of the composite system are obtained using numerical (i.e., finite element) or analytical (Mori-Tanaka) approaches. Compared with other homogenization methods, the proposed model is easy-to-use and versatile. The numerical model, validated by experimental results, is then used to guide the design of functional cementitious composites. The results show that CSP particle size (or relative shell thickness), volume fraction, shell property, and the size of the ITZ have significant impacts on the effective thermal and elastic properties of cementitious composites; whereas the particle size distribution pattern has relatively minor influence. In addition, the effective property of the composite is more sensitive to the particle property change when CSP inclusions are much stiffer or more thermally conductive than the matrix. Lastly, a scoping study is conducted to elucidate the relationship between the effective thermal conductivity and effective elastic moduli for cementitious composites containing CSP additives with different shell materials, particle sizes, and volume concentrations etc.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Orbital Chern Insulator and Quantum Phase Diagram of a Kagome Electron System with Half-Filled Flat Bands

In this work, we study the quantum phase diagram of electrons on kagome lattice with half-filled lowest flat bands by considering the antiferromagnetic Heisenberg interaction J, and short-range Coulomb interaction V. In the weak J regime, we identify a fully spin-polarized phase. The presence of finite V drives a spontaneous chiral current, which makes the system an orbital Chern insulator by contributing an orbital magnetization. Such an out-of-plane orbital magnetization allows the presence of a Chern insulating phase independent of the spin orientation in contrast to the spin-orbit coupling induced Chern insulator that disappears with inplane ferromagnetism constrained by symmetry. Such a symmetry difference provides a criterion to distinguish the physical origin of topological responses in kagome systems. The orbital Chern insulator is robust against small coupling J. By further increasing J, we find that the ferromagnetic topological phase is suppressed, which first becomes partially polarized and then enters a nonmagnetic phase with spin and charge nematicity. The frustrated flat band allows the spin and Coulomb interaction to play an essential role in determining the quantum phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

UQpy: A general purpose Python package and development environment for uncertainty quantification

In this paper, we present the UQpy software toolbox, an open-source Python package for general uncertainty quantification (UQ) in mathematical and physical systems. The software serves as both a user-ready toolbox that includes many of the latest methods for UQ in computational modeling and a convenient development environment for Python programmers advancing the field of UQ. The paper presents an introduction to the software's architecture and existing capabilities, divided in the code in a set of modules centered around different UQ tasks such as sampling methods, generation of random processes and random fields, probabilistic inverse modeling, reliability analysis, surrogate modeling, and active learning. The paper also highlights the importance of the RunModel module, which is used to drive simulations in the uncertainty analyses performed in UQpy. This module conveniently allows the user to define computational models directly in Python, or to run simulations from a third-party software in serial or in parallel. To illustrate the various capabilities, two examples are tracked throughout the paper and analyzed repeatedly for various UQ tasks. The first is a Python model solving a nonlinear structural dynamics problem, used to illustrate UQpy's capabilities in sampling and forward propagation of high dimensional random vectors (stochastic processes), and probabilistic inference. The second model is a third-party Abaqus finite element model solving the thermomechanical response of a beam structure. This example is used to illustrate UQpy's capabilities in variance reduction sampling techniques, reliability analysis, surrogate modeling and active learning techniques.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Zero-pairing and zero-temperature limits of finite-temperature Hartree-Fock-Bogoliubov theory

Recently, variational Hartree-Fock-Bogoliubov (HFB) mean-field equations were shown to possess a mathematically well-defined zero-pairing limit, independently of the closed- or open-shell character of the system under consideration. This limit is non-trivial for open-shell systems such that HFB theory does {\it not} reduce to the Hartree-Fock (HF) formalism in all cases. The present work extends this analysis to finite-temperature HFB (FTHFB) theory by investigating the behavior of this more general formalism in the combined zero-temperature and zero-pairing limits. The zero-pairing and zero-temperature limits of the FTHFB statistical density operator constrained to carry an arbitrary (integer) number of particles A on average is worked out analytically and realized numerically using a two-nucleon interaction. While the FTHFB density operator reduces to the projector corresponding to a pure HF Slater determinant for closed-shell nuclei, the FTHFB formalism does not reduce to the HF theory in all cases in the zero-temperature and zero-pairing limits, i.e. for open-shell nuclei. However, the fact that a nucleus can be of open-shell character in these joint limits is necessarily the result of some symmetry restrictions. Whenever it is the case, the non-trivial description obtained for open-shell systems is shown to depend on the order with which both limits are taken, i.e. the two limits do not commute for these systems. When the zero-temperature limit is performed first, the FTHFB density operator is demoted to a projector corresponding to a pure state made out of a linear combination of a finite number of Slater determinants with different (even) numbers of particles. When the zero-pairing limit is performed first, the FTHFB density operator remains a statistical mixture of a finite number of Slater determinants with both even and odd particle numbers. While the entropy (pairing density) is zero in the first (second) case, it does not vanish in the second (first) case in spite of the temperature (pairing) tending towards zero. The difference between both limits can have striking consequences for the (thermal) expectation values of observables. For instance, the particle-number variance does not vanish in either case and has limiting values that differ by a factor of two in both cases. In conclusion, while in the textbook situation associated with closed-shell nuclei Hartree-Fock-Bogoliubov (finite-temperature Hartree-Fock) theory reduces to Hartree-Fock theory in the zero-pairing (zero-temperature) limit, the present analysis demonstrates that a non trivial and unexpected limit is obtained for this formalism in open-shell systems. This result sheds a new light on certain aspects of this otherwise very well-studied many-body formalism.

20 ≤ A ≤ 38↗
NASA NTRS1994

Discrete-Layer Piezoelectric Plate and Shell Models for Active Tip-Clearance Control

The objectives of this work were to develop computational tools for the analysis of active-sensory composite structures with added or embedded piezoelectric layers. The targeted application for this class of smart composite laminates and the analytical development is the accomplishment of active tip-clearance control in turbomachinery components. Two distinct theories and analytical models were developed and explored under this contract: (1) a discrete-layer plate theory and corresponding computational models, and (2) a three dimensional general discrete-layer element generated in curvilinear coordinates for modeling laminated composite piezoelectric shells. Both models were developed from the complete electromechanical constitutive relations of piezoelectric materials, and incorporate both displacements and potentials as state variables. This report describes the development and results of these models. The discrete-layer theories imply that the displacement field and electrostatic potential through-the-thickness of the laminate are described over an individual layer rather than as a smeared function over the thickness of the entire plate or shell thickness. This is especially crucial for composites with embedded piezoelectric layers, as the actuating and sensing elements within these layers are poorly represented by effective or smeared properties. Linear Lagrange interpolation polynomials were used to describe the through-thickness laminate behavior. Both analytic and finite element approximations were used in the plane or surface of the structure. In this context, theoretical developments are presented for the discrete-layer plate theory, the discrete-layer shell theory, and the formulation of an exact solution for simply-supported piezoelectric plates. Finally, evaluations and results from a number of separate examples are presented for the static and dynamic analysis of the plate geometry. Comparisons between the different approaches are provided when possible, and initial conclusions regarding the accuracy and limitations of these models are given.

Heyliger, P. R.↗
DOE OSTI2015

A simple structure of all circular-air-holes photonic crystal fiber for achieving high birefringence and low confinement loss

We propose a simple structure of photonic crystal fibers (PCFs) with high birefringence and low confinement loss based on one rectangular centric ring of smaller circular air holes (CAHs) in the fiber core, and three rings of larger CAHs in the fiber cladding. This simple geometry (using all CAHs with two different air hole sizes) is capable of achieving a flexible control of the birefringence, B = 5.501 × 10{sup −3}, and ultra-low confinement loss, 7.30 × 10{sup −5 }dB/km, at an excitation wavelength of λ = 1550 nm. The birefringence value is ∼5.0 times greater than that obtained for conventional CAH PCF. This simple structure has the added advantage from the view point of easy fabrication, robustness, and cost. A full-vector finite element method combined with anisotropic perfectly matched layers was used to analyze the various fiber structures. We have analyzed four cases of CAH PCFs, focusing on the core asymmetry design as opposed to the conventional approach of CAHs or elliptical air holes on the cladding and core. The robustness against manufacturing inaccuracies of the proposed structure has also been further investigated in this work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

A General Study of Counterflow Diffusion Flames for Supercritical CO2 Combustion

Abstract A counterflow diffusion flame for supercritical CO2 combustion is investigated at various CO2 dilution levels and pressures by accounting for real gas effects into both thermal and transport properties. The UCF 1.1 24-species mechanism is used to account the chemistry. The nature of important nonpremixed combustion characteristics such as Prandtl number, thermal diffusivity, Lewis number, stoichiometric scalar dissipation rate, flame thickness, and Damköhler number are investigated with respect to CO2 dilution and pressure. The results show that the aforementioned parameters are influenced by both dilution and pressure; the dilution effect is more dominant. Further, the result shows that Prandtl number increases with CO2 dilution and at 90% CO2 dilution, the difference between the Prandtl number of the inlet jets and the flame is minimal. Also, the common assumption of unity Lewis number in the theory and modeling of nonpremixed combustion does not hold reasonable for sCO2 applications due to large difference of Lewis number across the flame and the Lewis number on the flame drop significantly with an increase in the CO2 dilution. An interesting relation between Lewis number and CO2 dilution is observed. The Lewis number of species drops by 15% when increasing the CO2 dilution by 30%. Increasing the CO2 dilution increases both the flow and chemical timescales; however, chemical timescale increases faster than the flow timescales. The magnitudes of the Damköhler number signify the need to consider finite rate chemistry for sCO2 applications. Further, the Damköhler numbers at 90% sCO2 dilution are very small; hence, laminar flamelet assumptions in turbulent combustion simulations are not physically correct for this application. Also, it is observed that the Damköhler number drops nonlinearly with increasing CO2 dilution in the oxidizer stream. This is a very important observation for the operation of sCO2 combustors. Further, the flame thickness is found to increase with CO2 dilution and reduce with pressure.

Engineering↗
DOE OSTI2022

3D simulations of spinlike flames in Co/Al multilayers with enhanced conduction losses

Reactive Co/Al multilayers are uniformly structured materials that may be ignited to produce rapid and localized heating. Prior studies varying the bilayer thickness (i.e., sum of two individual layers of Co and Al) have revealed different types of flame morphologies, including: (a) steady/planar, (b) wavy/periodic, and (c) transverse bands, originating in the flame front. These instabilities resemble the “spin waves” first observed in the early studies of solid combustion (i.e., Ti cylinder in a N 2 atmosphere), and are likewise thought to be due to the balance of heat released by reaction and heat conduction forward into the unreacted multilayer. However, the multilayer geometry and three-dimensional (3D) edge effects are relatively unexplored. In this work, a new diffusion-limited reaction model for Co/Al multilayers was implemented in large, novel 3D finite element analysis (FEA) simulations, in order to study the origins of these spinlike flames. This reaction model builds upon previous work by introducing three new phase-dependent property models for: (1) the diffusion coefficient, (2) anisotropic thermal conductivity tensor, and (3) bulk heat capacity, as well as one additional model for the bilayer-dependent heat of reaction. These novel 3D simulations are the first to predict both steady and unsteady flames in Co/Al multilayers. Moreover, two unsteady modes of flame propagation are identified, which depend on the enhanced conduction losses with slower flames, as well as flame propagation around notched edges. Finally, future work will consider the generality of the current modeling approach and also seek to define a more generalized set of stability criteria for additional multilayer systems.

42 ENGINEERING↗
DOE OSTI2016

Experimental demonstration of revival of oscillations from death in coupled nonlinear oscillators

We experimentally demonstrate that a processing delay, a finite response time, in the coupling can revoke the stability of the stable steady states, thereby facilitating the revival of oscillations in the same parameter space where the coupled oscillators suffered the quenching of oscillation. This phenomenon of reviving of oscillations is demonstrated using two different prototype electronic circuits. Further, the analytical critical curves corroborate that the spread of the parameter space with stable steady state is diminished continuously by increasing the processing delay. Finally, the death state is completely wiped off above a threshold value by switching the stability of the stable steady state to retrieve sustained oscillations in the same parameter space. The underlying dynamical mechanism responsible for the decrease in the spread of the stable steady states and the eventual reviving of oscillation as a function of the processing delay is explained using analytical results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Numerical Simulation of Critical Factors Controlling Heat Extraction from Geothermal Systems Using a Closed-Loop Heat Exchange Method

Closed-loop heat exchange for geothermal energy production involves injecting working fluid down a well that extends through the geothermal resource over a significant length to absorb heat by conduction through the well pipe. The well then needs to return to the surface for energy recovery and fluid re-injection to complete the cycle. We have carried out mixed convective-conductive fluid-flow modeling using a wellbore flow model for TOUGH2 called T2Well to investigate the critical factors that control closed-loop geothermal energy recovery. T2Well solves a mixed explicit-implicit set of momentum equations for flow in the pipe with full coupling to the implicit three-dimensional integral finite difference equations for Darcy flow in the porous medium. T2Well has the option of modeling conductive heat flow from the porous medium to the pipe by means of a semi-analytical solution which makes the computation very efficient because the porous medium does not have to be discretized. When the fully three-dimensional option is chosen, the porous medium is discretized and heat flow to the pipe is by conduction and convection, depending on reservoir permeability and other factors. Simulations of the closed-loop system for a variety of parameter values have been carried out to elucidate the heat recovery process. To the extent that convection may occur to aid in heat delivery to the pipe, the permeability of the geothermal reservoir, whether natural or stimulated, is an important property in heat extraction. The injection temperature and flow rate of the working fluid strongly control the ultimate energy recovery. Pipe diameter also plays a strong role in heat extraction, but is correlated with flow rate. Similarly, the choice of working fluid plays an important role, with water showing better heat extraction than CO2 for certain flow rates, while the CO2 has higher pressure at the production wellhead which can aid in surface energy recovery. In general, we find complex interactions between the critical factors that will require advanced computational approaches to fully optimize.

Oldenburg, Curtis↗
DOE OSTI2019

Development of a Modeling Approach to Describe Thermal Conductivity of Prototypic SiC-SiC tubes for LWR Fuel Cladding

SiC-SiC composites are a promising material intended to be used as a for fuel cladding for Light Water Reactors (LWRs). This material has a high steam oxidation resistance, high strength and low neutron cross section. Thermal conductivity is one of the most important physical properties for assessing the performance of SiC based fuel cladding. However, there is a significant lack of data for tube materials because measurement is challenging owing to the specimen geometry. This report summarizes the modeling methodology, experimental methodologies and software that will be used to study the thermal conductivity of different SiC-SiC composite tubes. The proposed experiments described here, are directed to create a baseline of thermal conductivity values to validate experimental simulations. In order to obtain the thermal conductivity of these materials the Laser Flash Analysis (LFA) will be implemented to obtain the thermal diffusivity of flat and sections of round shape specimens. Thermal diffusivity experiments will be complemented by Finite Element Method (FEM) simulations and a branch of this methodology known as Image-Based Finite Element Method (IBFEM). The FEM simulations will help to identify the possible limitations of experimental measurements in round shape specimens and verify the results of experiments in flat shape samples. The results of the IBFEM will help to quantify the possible effects of porosity in SiC-SiC composites. The findings of the aforementioned activities will be also applied to SiC-SiC neutron irradiated tube samples in order to capture the damage produced by this ageing mechanism.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗