SEARCH · Engineering Papers

Results for “complex fluid”

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 775 records · Page 43

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Probing multiscale dissolution dynamics in natural rocks through microfluidics and compositional analysis

Mineral dissolution significantly impacts many geological systems. Carbon released by diagenesis, carbon sequestration, and acid injection are examples where geochemical reactions, fluid flow, and solute transport are strongly coupled. The complexity in these systems involves interplay between various mechanisms that operate at timescales ranging from microseconds to years. Current experimental techniques characterize dissolution processes using static images that are acquired with long measurement times and/or low spatial resolution. These limitations prevent direct observation of how dissolution reactions progress within an intact rock with spatially heterogeneous mineralogy and morphology. We utilize microfluidic cells embedded with thin rock samples to visualize dissolution with significant temporal resolution (100 ms) in a large observation window (3 × 3 mm). Here we injected acidic fluid into eight shale samples ranging from 8 to 86 wt % carbonate. The pre- and postreaction microstructures are characterized at the scale of pores (0.1 to 1 µm) and fractures (1 to 1,000 µm). We observe that nonreactive particle exposure, fracture morphology, and loss of rock strength are strongly dependent on both the relative volume of reactive grains and their distribution. Time-resolved images of the rock unveil the spatiotemporal dynamics of dissolution, including two-phase flow effects in real time and illustrate the changes in the fracture interface across the range of compositions. Moreover, the dynamical data provide an approach for characterizing reactivity parameters of natural heterogeneous samples when porous media effects are not negligible. The platform and workflow provide real-time characterization of geochemical reactions and inform various subsurface engineering processes.

58 GEOSCIENCES↗

DOE OSTI2019

Simulation of discrete cracks driven by nearly incompressible fluid via 2D combined finite-discrete element method

Here, we propose a novel hydraulic solver in order to simulate key mechanisms that control fluid-driven cracks in the framework of the combined finite-discrete element method (FDEM). The main innovative aspect of the present work is the independence of the fluid's critical time step size with the fracture opening. This advantage is extremely important because it means that very fine meshes can be used around areas of interest, such as boreholes, without penalizing the computational cost as fractures propagate (ie, open) and the fluid flows through them. This is a great advantage over other recently introduced approaches that exhibit a dependency of the time step in the form of Δt crit ∝ (l/a) 2 where l is the element size and a is the fracture aperture. This paper presents a series of benchmark cases for the proposed solver. The rationale adopted by the authors was to benchmark and validate the implementation of the hydraulic solver in an incremental fashion, starting from the simplest cases and building in complexity. The calculations shown in this work clearly demonstrate that the proposed approach is able to reproduce analytical results for fluid flow through a single crack. The results presented in this paper also demonstrate that the new approach is robust enough to deal with complex fracture patterns and complex geometries; the obtained fluid-driven fracture patterns in the vicinity of a borehole certainly stand to the scrutiny of human visual perception.

58 GEOSCIENCES↗

NASA NTRS1998

Overview of Sensitivity Analysis and Shape Optimization for Complex Aerodynamic Configurations

This paper presents a brief overview of some of the more recent advances in steady aerodynamic shape-design sensitivity analysis and optimization, based on advanced computational fluid dynamics. The focus here is on those methods particularly well- suited to the study of geometrically complex configurations and their potentially complex associated flow physics. When nonlinear state equations are considered in the optimization process, difficulties are found in the application of sensitivity analysis. Some techniques for circumventing such difficulties are currently being explored and are included here. Attention is directed to methods that utilize automatic differentiation to obtain aerodynamic sensitivity derivatives for both complex configurations and complex flow physics. Various examples of shape-design sensitivity analysis for unstructured-grid computational fluid dynamics algorithms are demonstrated for different formulations of the sensitivity equations. Finally, the use of advanced, unstructured-grid computational fluid dynamics in multidisciplinary analyses and multidisciplinary sensitivity analyses within future optimization processes is recommended and encouraged.

Newman, Perry A.↗

DOE OSTI2024

Crack opening calculation in phase-field modeling of fluid-filled fracture: A robust and efficient strain-based method

The phase-field method has become popular for the numerical modeling of fluid-filled fractures, thanks to its ability to represent complex fracture geometry without algorithms. However, the algorithm-free representation of fracture geometry poses a significant challenge in calculating the crack opening (aperture) of phase-field fracture, which governs the fracture permeability and hence the overall hydromechanical behavior. Although several approaches have been devised to compute the crack opening of phase-field fracture, they require a sophisticated algorithm for post-processing the phase-field values or an additional parameter sensitive to the element size and alignment. Here, we develop a novel method for calculating the crack opening of fluid-filled phase-field fracture, which enables one to obtain the crack opening without additional algorithms or parameters. Here we transform the displacement-jump-based kinematics of a fracture into a continuous strain-based version, insert it into a force balance equation on the fracture, and apply the phase-field approximation. Through this procedure, we obtain a simple equation for the crack opening which can be calculated with quantities at individual material points. We verify the proposed method with analytical and numerical solutions obtained based on discrete representations of fractures, demonstrating its capability to calculate the crack opening regardless of the element size or alignment.

58 GEOSCIENCES↗

NASA NTRS2015

Unsteady Three-Dimensional Simulation of a Shear Coaxial GO2/GH2 Rocket Injector with RANS and Hybrid-RAN-LES/DES Using Flamelet Models

Historically, the analysis and design of liquid rocket engines (LREs) has relied on full‐scale testing and one‐dimensional empirical tools. The testing is extremely expensive and the one‐dimensional tools are not designed to capture the highly complex, and multi‐dimensional features that are inherent to LREs. Recent advances in computational fluid dynamics (CFD) tools have made it possible to predict liquid rocket engine performance, stability, to assess the effect of complex flow features, and to evaluate injector‐driven thermal environments, to mitigate the cost of testing. Extensive efforts to verify and validate these CFD tools have been conducted, to provide confidence for using them during the design cycle. Previous validation efforts have documented comparisons of predicted heat flux thermal environments with test data for a single element gaseous oxygen (GO2) and gaseous hydrogen (GH2) injector. The most notable validation effort was a comprehensive validation effort conducted by Tucker et al. [1], in which a number of different groups modeled a GO2/GH2 single element configuration by Pal et al [2]. The tools used for this validation comparison employed a range of algorithms, from both steady and unsteady Reynolds Averaged Navier‐Stokes (U/RANS) calculations, large‐eddy simulations (LES), detached eddy simulations (DES), and various combinations. A more recent effort by Thakur et al. [3] focused on using a state‐of‐the‐art CFD simulation tool, Loci/STREAM, on a two‐dimensional grid. Loci/STREAM was chosen because it has a unique, very efficient flamelet parameterization of combustion reactions that are too computationally expensive to simulate with conventional finite‐rate chemistry calculations. The current effort focuses on further advancement of validation efforts, again using the Loci/STREAM tool with the flamelet parameterization, but this time with a three‐dimensional grid. Comparisons to the Pal et al. heat flux data will be made for both RANS and Hybrid RANSLES/ Detached Eddy simulations (DES). Computation costs will be reported, along with comparison of accuracy and cost to much less expensive two‐dimensional RANS simulations of the same geometry.

Westra, Doug G.↗

DOE OSTI2023

Intermediate Temperature Proton‐Conducting Solid Oxide Electrolysis Cells with Improved Performance and Durability

Technical understanding: A high-performance and stable intermediate temperature proton-conducting solid oxide electrolysis cell (PCEC) will add significantly to the hydrogen roadmap of EERE. However, the low electrocatalytic activity of the air electrode and the instability of the electrolyte are technical barriers. In this research, we utilized computational fluid dynamic modeling and materials/microstructure engineering to investigate the complex reaction chain and to improve the performance and reliability of the protonic ceramic electrochemical cells. An active and stable Pr 1.8 Ba 0.2 NiO 4+δ -BaZr 0.1 Ce 0.7 Y 0.1 Yb 0.1 O 3-δ (BZCYYb) scaffold conformal coating electrode was rationally designed where the robust Pr 1.8 Ba 0.2 NiO 4+δ functions both as a catalytic active and a protective layer. The PCECs with porosity-optimized scaffold conformal coating electrodes demonstrate 1.46 Acm -2 under 1.3 V at 600°C, which is among the highest in recent years. The final optimization by computational fluid dynamic modeling and conformal coating technique was applied to evolve the PCEC air electrode from Pr 2 NiO 4+δ to Pr 1.8 Ba 0.2 NiO 4+δ -BaZr 0.1 Ce 0.7 Y 0.1 Yb 0.1 O 3-δ scaffold conformal coating. The electrolysis current was boosted from 0.445A/cm 2 to 1.46A/cm 2 at 600°C and 1.3V. A superior low degradation rate of 1mV/kh was achieved in a 5000-hour of operation. Technical effectiveness and economic feasibility: The effectiveness of this TCO conformal coating on electrolyte backbone as air electrode has been proven by long-term demonstration. Based on the highly conductive BZCYYb electrolyte and the TCO steam electrode PNO, this project rationally took advantage of the existent materials to structure a PCEC device with high TRL and low risk. The solution infiltration method used to structure the conformal coating is a standard method used by the solid oxide cell community. Very minimal adoption is needed in the current SOC fabrication process. It is a cost-effective, scalable method for manufacturing, and brings up remarkable impact on the current technical status of PCEC. Benefit to the public: The findings in the research showed the industry that the severe challenge in PCEC can be overcome in a practical way. A various energy-hydrogen ecosystem is not beyond reach. The success of this technology provides sustainable and affordable H 2 to society. It will stimulate investment from the governmental and private sectors in response to the device manufacture and deployment, and the related infrastructure construction. It will generate high-pay domestic jobs and environmental benefits to the general public in the United States.

08 HYDROGEN↗

DOE OSTI2020

Wind Farm Simulation and Layout Optimization in Complex Terrain: Preprint

This work reports on incorporating complex terrain into wind farm simulations for the purpose of layout optimization. Adding complex terrain boundary conditions to NREL's medium fidelity computational fluid dynamics model, WindSE, produces significant separation, flow curvature, and speedup effects that would otherwise be difficult to capture with lower-fidelity models or a flat-terrain assumption. These flow features, in turn, can significantly impact the optimal turbine array layout. We demonstrate the impact of complex terrain on flow in both an idealized and real-world setting, and discuss modifications to the code that enable gradient-based optimization using terrain-aware adjoint gradients. Through several optimization case studies, we show that the layout optimization process takes advantage of speedup effects on terrain high points, and leverages flow curvature effects that modify wake trajectories. This yields substantial power improvements over gridded layouts, and hints at future research directions in simulation and optimization for wake trajectories in complex terrain.

17 WIND ENERGY↗

NASA NTRS1991

Analysis and representation of complex structures in separated flows

We discuss our recent work on extraction and visualization of topological information in separated fluid flow data sets. As with scene analysis, an abstract representation of a large data set can greatly facilitate the understanding of complex, high-level structures. When studying flow topology, such a representation can be produced by locating and characterizing critical points in the velocity field and generating the associated stream surfaces. In 3D flows, the surface topology serves as the starting point. The 2D tangential velocity field near the surface of the body is examined for critical points. The tangential velocity field is integrated out along the principal directions of certain classes of critical points to produce curves depicting the topology of the flow near the body. The points and curves are linked to form a skeleton representing the 2D vector field topology. This skeleton provides a basis for analyzing the 3D structures associated with the flow separation. The points along the separation curves in the skeleton are used to start tangent curve integrations. Integration origins are successively refined to produce stream surfaces. The map of the global topology is completed by generating those stream surfaces associated with 3D critical points.

Helman, James↗

DOE OSTI2018

Lagrangian particle method for compressible fluid dynamics

A new Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is suitable for the simulation of complex free surface / multi-phase flows. The main contributions of our method, which is different from SPH in all other aspects, are (a) significant improvement of approximation of differential operators based on a polynomial fit via weighted least squares approximation and the convergence of prescribed order, (b) a second-order particle-based algorithm that reduces to the first-order upwind method at local extremal points, providing accuracy and long term stability, and (c) more accurate resolution of entropy discontinuities and states at free inter-faces. While the method is consistent and convergent to a prescribed order, the conservation of momentum and energy is not exact and depends on the convergence order . The method is generalizable to coupled hyperbolic-elliptic systems. As a result, numerical verification tests demonstrating the convergence order are presented as well as examples of complex multiphase flows.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Final Scientific/Technical Report: Correlations and Fluctuations in Weakly Collisional Plasma

Plasma is a state of matter that exhibits a very rich range of phenomena. To begin with, plasma is both electrical and mechanical - bringing together theories of particle motion and the electromagnetic field. Furthermore, and especially important for this project, a weakly-collisional plasma, such as is found in high-temperature (fusion energy) experiments on earth and the majority of contexts in space and astrophysics, has many moving parts. For example, sitting in earth’s atmosphere we are immersed in a mechanical wave field (sound), a possibly turbulent fluid motion (wind), and an electromagnetic vector wave field with two polarizations (light). This is already enough to produce a rich range of possibilities. In plasma, the electromagnetic field is coupled to the mechanical motion of the medium because it is ionized. Furthermore, a weakly-collisional plasma supports an infinite number of mechanically independent fluids. Thus, plasmas support an infinite number of independent electromechanical waves. Much has been done to describe plasmas with "reduced models" of various kinds. The goal of this project was to both explore the validity of reduced plasma models that are in use, and to propose and validate new models of plasma motion. The primary means to his end was laboratory experiments employing both electrical probes and laser spectroscopy. Laser spectroscopy enables many techniques which can separate the spectrum of independent fluid motions in the ion phase-space. The choice was to focus on low frequency electrostatic waves because the electron motion is relatively simple, the experiments can be on a spatial scale of a few meters, and all the relevant parameters can be measured with a few lasers systems. No study of this kind had previously been undertaken for the study of plasmas. The validation of theories required that the experimental descriptions be compared with theory and simulation in detail. It was found that even multi-fluid theories leave out a large part of the complexity of plasma motion. Reduced descriptions were found to fail under most circumstances. A new technique was developed that enabled a measurement of the phase-space resolved ion correlation function for the first time. The wide range of plasma dynamics possible became clear through this technique. It was found that collisionless (Vlasov) theory has a large field of application even when the plasma is weakly-collisional. A new approach, the kinetic wave expansion, was proposed, tested and found to be very useful for describing electrostatic ion waves. This project demonstrated a new way of looking at the "degrees-of-freedom" of plasmas and provided significant validation tests of fluid and kinetic plasma descriptions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Evidence for Carbonate Surface Complexation during Forsterite Carbonation in Wet Supercritical Carbon Dioxide

Continental flood basalts are attractive formations for geologic sequestration of carbon dioxide because of their reactive divalent-cation containing silicates, such as forsterite (Mg2SiO4), suitable for long-term trapping of CO2 mineralized as metal carbonates. The goal of this study was to investigate at a molecular level the carbonation products formed during the reaction of forsterite with supercritical CO2 (scCO2) as a function of the concentration of H2O adsorbed to the forsterite surface. Experiments were performed at 50 °C and 90 bar using an in situ IR titration capability, and post-reaction samples were examined by ex situ techniques, including SEM, XPS, FIB-TEM, TGA-MS, and MAS-NMR. Carbonation products and reaction extents varied greatly with adsorbed H2O. We show for the first time evidence of Mg-carbonate surface complexation under wet scCO2 conditions. Carbonate is found to be coordinated to Mg at the forsterite surface in a predominately bidentate fashion at adsorbed H2O concentrations below 27 µmol/m2. Above this concentration and up to 76 µmol/m2, monodentate coordinated complexes become dominant. Beyond a threshold adsorbed H2O concentration of 76 µmol/m2, crystalline carbonates continuously precipitate as magnesite, and the particles that form are hundreds of times larger than the estimated thicknesses of the adsorbed water films of about 7 to 15 Å. At an applied level, these results suggest that mineral carbonation in scCO2 dominated fluids near the wellbore and adjacent to caprocks will be insignificant and limited to surface complexation, unless adsorbed H2O concentrations are high enough to promote crystalline carbonate formation. At a fundamental level, the surface complexes and their dependence on adsorbed H2O concentration give insights regarding forsterite dissolution processes and magnesite nucleation and growth.

Environmental Molecular Sciences Laboratory↗

DOE OSTI2015

Modeling the Interaction Between Hydraulic and Natural Fractures Using Dual-Lattice Discrete Element Method

The interaction between hydraulic fractures (HF) and natural fractures (NF) will lead to complex fracture networks due to the branching and merging of natural and hydraulic fractures in unconventional reservoirs. In this paper, a newly developed hydraulic fracturing simulator based on discrete element method is used to predict the generation of complex fracture network in the presence of pre-existing natural fractures. By coupling geomechanics and reservoir flow within a dual lattice system, this simulator can effectively capture the poro-elastic effects and fluid leakoff into the formation. When HFs are intercepting single or multiple NFs, complex mechanisms such as direct crossing, arresting, dilating and branching can be simulated. Based on the model, the effects of injected fluid rate and viscosity, the orientation and permeability of NFs and stress anisotropy on the HF-NF interaction process are investigated. Combined impacts from multiple parameters are also examined in the paper. The numerical results show that large values of stress anisotropy, intercepting angle, injection rate and viscosity will impede the opening of NFs.

02 PETROLEUM↗

DOE OSTI2022

Continued Verification of MOOSE Structural Mechanics Tools for Modeling Core Bowing Phenomena in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY22, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21. Several benchmark verification examples were selected from the IAEA verification and validation report. These examples involve clusters of ducts representative of a sector of a hexagonal reactor core which bow into each other and cause contact and load pad elevations, as well as single ducts subjected to irradiation fields undergoing swelling and subsequent bowing. The MOOSE-based results were compared to both IAEA benchmark participants’ results, analytic equations as available, and NUBOW-3D, a beam model code developed by Argonne National Laboratory. In every case, the MOOSE results agreed with other simulations results, providing additional verification basis of the tools for this particular physics application.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

FY23 Status Report on MOOSE-Based Approaches to Modeling Core Bowing in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY23, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21 and FY22. Several benchmark verification examples were selected from the IAEA verification and validation report with increased number of ducts and more complex contact interaction behavior. These examples involve a full symmetric sector with restraint rings at multiple load pad locations to simulate a limited free-bow restraint system concept, as well as irradiation induced swelling and creep effects in a sector. In addition, improvements to the contact module sideset assignment were assessed and compared with previous MOOSE results to verify the contact behavior. The MOOSE-based results were compared to IAEA benchmark participants’ results. In every case, the MOOSE results agreed with other simulations results for estimating bowing behavior, providing additional verification basis of the tools for this particular physics application. Estimation of contact forces were mostly in agreement, with a few outlier results. A plan was suggested for dealing with the discrepancies with estimating contact force values. In addition, a 1-way coupling demonstration was performed using subchannel analysis code Pronghorn-SC and MOOSE on an ABR-1000-design sodium-cooled fast reactor assembly to evaluate coolant and duct wall temperatures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2019

Extrusion of Complex Surface Meshes Utilizing Face Offsetting and Mean Curvature Smoothing

Grids for three-dimensional Computational Fluid Dynamics problems frequently require a prismatic layer of cells, typically extruded in the off-body direction from a two-dimensional surface mesh to properly resolve boundary layers. When the surface geometry is complex, the extrusion process can fail, resulting in the physical boundary layer being thicker than the prismatic layer, leading to under-resolved boundary layers. To address the shortcomings of existing grid generation tools, a new tool has been developed as part of the Mesh Tools suite, that is capable of extruding complex surface meshes. The new tool uses a face offsetting method to preserve surface curvature and a mean curvature smoothing algorithm to prevent the cells from self-intersecting in concave regions. Testing of the new tool found that not only were complex surface meshes able to be extruded to the desired thickness, but extrusion of simple surface meshes was also improved, due to the face offsetting method.

mesh extrusion↗

DOE OSTI2024

Catalytic resonance theory: Circumfluence of programmable catalytic loops

Chemical reactions on heterogeneous catalyst surfaces exhibit complex networks of elementary reactions with multiple pathways to fluid phase products, sometimes leading to surface reaction loops consisting of a closed cycle reaction pathway. While conventional catalysts at steady state exhibit zero net flux in either direction around a catalytic loop, the loop turnover frequency of three-species surface loops was evaluated in this work via microkinetic modeling to assess the reaction loop behavior resulting from a catalytic surface oscillating between two or more surface energy states. For dynamic heterogeneous catalysts undergoing applied oscillations of surface energy (i.e., programs), surface reaction loops of three species were shown to exhibit non-zero net flow of molecules around the loop, with the extent of loop turnover frequency varying with the applied frequency and amplitude of a sinusoidal or square wave oscillation. Alternatively, some dynamic surface reactions exhibited oscillations only between two surface species or resulted in surfaces covered by a single species at all times. Here, more complex behavior was observed for dynamic surface programs of three distinct electronic states, with the temporal order of applied states controlling the direction of flow of molecules within a three-molecule surface loop. Catalytic loops have the potential to limit overall catalytic reaction rates and use energy in programmable catalysts, while some applications could purposely impose non-zero loop turnover frequency for improved surface reaction control.

Catalysis↗

NASA NTRS2001

B-spline Method in Fluid Dynamics

B-spline functions are bases for piecewise polynomials that possess attractive properties for complex flow simulations : they have compact support, provide a straightforward handling of boundary conditions and grid nonuniformities, and yield numerical schemes with high resolving power, where the order of accuracy is a mere input parameter. This paper reviews the progress made on the development and application of B-spline numerical methods to computational fluid dynamics problems. Basic B-spline approximation properties is investigated, and their relationship with conventional numerical methods is reviewed. Some fundamental developments towards efficient complex geometry spline methods are covered, such as local interpolation methods, fast solution algorithms on cartesian grid, non-conformal block-structured discretization, formulation of spline bases of higher continuity over triangulation, and treatment of pressure oscillations in Navier-Stokes equations. Application of some of these techniques to the computation of viscous incompressible flows is presented.

Botella, Olivier↗

DOE OSTI2017

Verification of the Coupled Neutronics/Thermal-Hydraulics Code TANSY against the MSFR Benchmark

The MSRs is a liquid-fuel reactor in which uranium, fission products and actinides are dissolved in the molten salt. The molten salt reactors are featured by circulating fuel, low pressure and high temperature. The MSRs meets many of the future goals of nuclear technology, in particular for what concerns as improved sustainability and unique characteristics in terms of actinide burning. The fluid nature of fuel in MSRs gives extra flexibility in reactor design, fuel fabrication and recycling. But because of the application of fluid fuel, MSRs have several special characteristics compared with the traditional solid-fuel reactors. Firstly, the transport of the delayed neutron precursors by the fluid flow reduces the contribution of delayed neutrons to the chain reaction.. Secondly, the fact that the fuel is dissolved in the coolant rather than separated from the coolant by the claddings results in a much stronger coupling phenomenon between the neutronics and thermal-hydraulics. Thirdly, the shape of the fuel is determined by the container because of the characteristics of fluid, the unstructured meshes are required to model the complex geometry accurately. However most of the conventional reactor codes were developed for solid-fuel reactor and therefore not capable of taking liquid fuel flow into account. To address all these particular modeling issues of MSRs, a multi-physics code named TANSY has been developed in our laboratory. In this paper, the benchmark problem based on the MSFR design was calculated to verify the code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗