Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fluids”

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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quantifying transport and electrocatalytic reaction processes in a gastight rotating cylinder electrode reactor via integration of Computational Fluid Dynamics modeling and experiments

Understanding the complexity of the multiple processes of mass, momentum, charge, and heat transport, and how these affect reaction kinetics at the electrode/electrolyte interface is one of the major challenges in the field of energy and catalysis. The rapid and rational scale-up of electrocatalytic systems to industrial scales require a detailed understanding of nonlinear transport-reaction processes, accessible only through the building of multi-physics models that capture with high fidelity the complexity of real-world devices. The gastight rotating cylinder electrode (RCE) reactor is a promising lab-scale tool that can decouple transport from intrinsic kinetics to generate data for first-principle models useful in the design of industrial, electrochemical reactors. Computational Fluid Dynamics (CFD) studies have previously been used to investigate the bulk flow in RCE reactors for simple corrosion and electroplating processes. However, the quantification of changes in local concentration within the viscous layer where catalysis takes place requires capturing the correct flow conditions inside the hydrodynamic boundary layer near the surface of the electrode. Further, this requires simulations with spatial resolution in the nm and μm scale and temporal resolutions between ms and s scales that are similar to the timescales for reactions on the electrode surface. In this study, experimental electrocatalysis is combined with CFD modeling to elucidate and parameterize the hydrodynamics in a gastight RCE reactor. CFD simulations of the electrochemical ferricyanide reduction reaction under mass transport limited conditions are used to evaluate the validity of the CFD model parameters by comparing calculated dimensionless mass transport descriptors to dimensionless correlations obtained experimentally. Justifications for assumptions and details of the simulation methods used in this study are presented to provide a detailed understanding of the effect that each model parameter has on the ability to accurately simulate electrocatalysis in RCE systems. The simulation methodology reported here is a first step towards the development of multi-scale models for the study of transport dependent electrocatalytic processes, such as the electrochemical transformation of CO 2 to fuels and chemicals.

42 ENGINEERING↗
DOE OSTI2019

Development of an Ultra-Low Head Siphon Hydro Turbine Using Computational Fluid Dynamics

Structural simplification and efficient hydraulic performance are key to effectively utilizing ultra-low head water power resources and reducing hydroelectric unit costs. In this study, a computational fluid dynamics (CFD) method was used to predict hydraulic performance of an axial turbine at the Gaoliangjian power station. CFD results agreed well with field test data. Using the same numerical method, a new siphon turbine was then designed based on the original distributor and turbine runner equipment. This study investigated the effect of different siphon outlet passage geometry parameters, runner blade shapes, and distributors on the hydraulic performance of the siphon turbine. The maximum hydraulic efficiency increased to 87.9% under a head of 2.87 m. Finally, the hydraulic performance of the turbine was compared for four different distributor designs after adding an intake sump at the turbine entrance. The bell-shaped distributor with four guide vanes resulted in the highest power output at the lowest head. Therefore, the siphon turbine is a good option for energy conversion in ultra-low water head settings.

Zhou, Daqing↗
DOE OSTI2022

Boron isotopes in blue diamond record seawater-derived fluids in the lower mantle

Boron is a quintessential crustal element but its conspicuous presence in diamond – a mantle mineral – raises questions about potential subduction pathways for boron and other volatiles. It has been a long-standing goal to characterize the isotopic composition of boron in blue, boron-bearing (Type IIb) diamonds to reveal its origin. Mineral inclusions indicate that Type IIb diamonds crystallize at transition zone to lower mantle depths, meaning that if the boron is subducted it would trace a pathway of volatile elements into the deep mantle. Here, using off-line laser ablation sampling, we present the first boron isotope compositions, along with trace element contents and carbon isotope compositions, of a suite of blue diamonds mainly sourced from the Cullinan diamond mine in South Africa. The ten analyzed blue diamonds have a wide range in δ 11 B, between -9.2 ± 2‰ and -0.5 ± 2‰, compared to the more restricted range for mid-ocean ridge basalts (-7.1 ± 0.9‰). Carbon isotope values for the blue diamonds range between -20.6 and -1.8‰, with a mode at -17‰, significantly more negative than the main mantle mode at -5‰. Combined, the boron and carbon isotope compositions require fluid input from subducted oceanic lithosphere down into deep mantle source regions of blue diamonds. This finding highlights a deep subduction pathway for volatiles to stoke the source regions of deeply-derived magmas.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Toward an Understanding of Surface Layer Formation, Growth, and Transformation at the Glass–Fluid Interface

Silicate glass is a metastable and durable solid that has application to a number of energy and environmental challenges (e.g., microelectronics, fiber optics, and nuclear waste storage). If allowed to react with water over time silicate glass develops an altered layer at the solid-fluid interface. In this study, we used borosilicate glass (LAWB45) as model material to develop a robust understanding of altered layer formation (i.e., amorphous hydrated surface layer and crystalline reaction products). Experiments were conducted at high surface area-to-volume ratio (~200,000 m-1) and 90 °C in the pressurized unsaturated flow (PUF) apparatus for 1.5-years to facilitate the formation of thick altered layers and allow for the effluent solution chemistry to be monitored continuously. A variety of microscopy techniques were used to characterized reacted grains and suggest the average altered layer thickness is 13.2 ± 8.3 μm with the hydrated and clay layer representing 74.8% and 25.2% of the total altered layer, respectively. This estimate is within the experimental error of the value estimated from the B release rate data (~10 ±1 μm/yr) over the 1.5-year duration. PeakForce® quantitative nanomechanical mapping results suggest the hydrated layer has a modulus that ranges between ~20 to 40 GPa, which is in the range of porous silica that contains from ~20 to ~50% porosity, yet significantly lower than dense silica (~70 to 80 GPa). Scanning transmission electron microscopy (STEM) images confirm the presence of pores and an analysis of a higher resolution image of a region provides a qualitative estimate of ≥ 22% porosity in this layer with variations in the hydrated layer in void volume with increasing distance from the unaltered glass. Chemical composition analyses, based on a combination of time-of-flight secondary-ion mass spectrometry (ToF-SIMS) and scanning electron microscopy with X-ray energy dispersive spectroscopy (EDS) and STEM-EDS, clearly show that the altered layer is mainly composed of Al, H, Si, and O with the clay layer being enriched in Li, Zn, Fe, and Mg. The amorphous hydrated layer is enriched in Ca, H, and Zr with a minor amount of K. Furthermore, ToF-SIMS results also suggest the B profile is anti-correlated with the H profile in the hydrated layer. Our selected-area electron diffraction results suggest the structure of the hydrated layer closely resembles opal-AG (amorphous gel-like) with an average crystallite size of ~0.7 nm which is smaller than the critical nucleus for silica nanoparticles (i.e., 1.4 to 3 nm). These results suggest the hydrated layer is more consistent with a polymeric gel rather than a colloidal gel and is comprised of molecular units (<1 nm in size) that result from the difficult to hydrolyze bonds, such as Si—O—Zr units, during the glass corrosion process. The size of individual particles or molecular units is a function of formation conditions (e.g., pH, ionic strength, nano-confinement, solute composition) in the hydrated layer.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling gas migration through clay-based buffer material using coupled multiphase fluid flow and geomechanics with stress-dependent gas permeability

A model for gas migration through clay-based buffer material is developed for modeling gas generation and migration associated with deep geologic nuclear waste disposal. The model is based on a multiphase fluid flow and geomechanics simulator that is adapted to consider enhanced gas flow when gas pressure is high enough to approach the confining stress magnitude. A key feature in the model is a direct coupling between gas permeability and stress, through a non-linear stress-dependent permeability function. The model was first tested and calibrated by modelling two different laboratory gas migration tests on Wyoming (MX-80) bentonite samples. The calibrated model was then applied to model gas migration through a bentonite buffer of a large-scale gas injection test (Lasgit) conducted at the Äspö Hard Rock Laboratory in Sweden. Observed preferential gas migration along interfaces (between compacted blocks and along the canister surface) required explicit representation of such interfaces in the model. The model with the stress-dependent gas permeability accurately captured observed experimental responses in terms of gas breakthrough time, peak gas pressure, and cumulative gas flow rates. The calibrated model was finally applied to simulate migration of hydrogen gas generated within a breached nuclear waste canister over 10,000 years, involving migration of much larger gas volumes. For the considered gas generation rate and host rock properties, the generated gas could migrate through the bentonite buffer and released into the surrounding host rock at a maximum gas pressure somewhat higher than the initial total stress, though a significant amount of hydrogen remained within the buffer. This modelling sets the stage for further detailed analysis of the impact of hydrogen gas generation on the long-term performance of nuclear waste repositories.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

A study on computational fluid dynamics modeling of a refrigerated container for COVID-19 vaccine distribution with experimental validation

A key issue with the distribution of vaccines to prevent COVID-19 is the temperature level required during transport, storage, and distribution. Typical refrigerated transport containers can provide a temperature-controlled environment down to -30 °C. However, the Pfizer vaccine must be carefully transported and stored under a lower temperature between -80 °C and -60 °C. One way to provide the required temperature is to pack the vaccine vials into small packages containing dry ice. Dry ice sublimates from a solid to a gas, which limits the allowable transport duration. This can be mitigated by transporting in a -30 °C refrigerated container. Moreover, because the dry ice will sublimate and thereby release CO 2 gas into the transport container, monitoring the CO 2 concentration within the refrigerated container is also essential. In the present work, a 3D computational fluid dynamics model was developed based on a commercially available refrigerated container and validated with experimental data. The airflow, temperature distribution, and CO 2 concentration within the container were obtained from the simulations. The modeling results can provide guidance on preparing experimental setups, thus saving time and lowering cost, and also provide insight into safety precautions needed to avoid hazardous conditions associated with the release of CO 2 during vaccine distribution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

Source Term Modeling for Evaluating the Potential Impacts to Groundwater of Fluids Escaping from a Depleted Oil Reservoir Used for Carbon Sequestration

In recent years depleted oil reservoirs have received special interest as carbon storage reservoirs because of their potential to offset costs through collaboration with enhanced oil recovery projects. Modeling is currently being conducted to evaluate potential risks to groundwater associated with leakage of fluids from depleted oil reservoirs used for storage of CO2. Modeling results reported here focused on understanding how toxic organic compounds found in oil will distribute between the various phases within a storage reservoir after introduction of CO2, understanding the migration potential of these compounds, and assessing potential groundwater impacts should leakage occur. Two model scenarios were conducted to evaluate how organic components in oil will distribute among the phases of interest (oil, CO2, and brine). The first case consisted of 50 wt.% oil and 50 wt.% water; the second case was 90 wt.% CO2 and 10 wt.% oil. Several key organic compounds were selected for special attention in this study based upon their occurrence in oil at significant concentrations, relative toxicity, or because they can serve as surrogate compounds for other more highly toxic compounds for which required input data are not available. The organic contaminants of interest (COI) selected for this study were benzene, toluene, naphthalene, phenanthrene, and anthracene. Partitioning of organic compounds between crude oil and supercritical CO2 was modeled using the Peng-Robinson equation of state over temperature and pressure conditions that represent the entire subsurface system (from those relevant to deep geologic carbon storage environments to near surface conditions). Results indicate that for a typical set of oil reservoir conditions (75°C, and 21,520 kPa) negligible amounts of the COI dissolve into the aqueous phase. When CO2 is introduced into the reservoir such that the final composition of the reservoir is 90 wt.% CO2 and 10 wt.% oil, a significant fraction of the oil dissolves into the vapor phase. As the vapor phase moves up through the stratigraphic column, pressures and temperatures decrease, resulting in significant condensation of oil components. The heaviest organic components condense early in this process (at higher pressures and temperatures), while the lighter components tend to remain in the vapor phase until much lower pressures and temperatures are reached. Based on the model assumptions, the final concentrations of COI to reach an aquifer at 1,520 kPa and 25°C were quite significant for benzene and toluene, whereas the concentrations of polynuclear aromatic hydrocarbons that reach the aquifer were very small. This work demonstrates a methodology that can provide COI source term concentrations in CO2 leaking from a reservoir and entering an overlying aquifer for use in risk assessments.

carbon sequestration, scCO2, crude oil, benzene, B↗
DOE OSTI2019

Validating compositional fluid flow simulations using 4D seismic interpretation and vice versa in the SECARB Early Test—A critical review

We report on strengths and weaknesses of 4D seismic interpretation as a technique for monitoring carbon dioxide that was injected as part of a large scale test associated with commercial enhanced oil recovery (EOR) at Cranfield Field, Mississippi, USA. The goals of the monitoring effort are 1) to make measurements to verify that the CO 2 is contained in the reservoir according to operational designs and model predictions, and 2) that if there are deviations, to provide data which can be used to update the earth models and determine if any mitigation is needed. Our work uses a compositional numerical simulation to model CO 2 flow in the reservoir and compares the results with estimates of changes in seismic properties between a pre-injection and a survey after more than 2 million metric tons injected. The complicated physics of the problem in Cranfield field present challenges to seismic interpretations. These findings reflect partial agreement between the results of the numerical simulation and the seismic interpretations. Possible causes of discrepancies among the fluid flow model and multiple interpretations of the same seismic data sets performed with different work flows illuminates the types of uncertainties that should be considered to achieve the goals of monitoring.

4D seismic↗
DOE OSTI2022

Experimental measurements of fluid flow in an 84-pin hexagonal rod bundle with spacer grid for a gas-cooled fast modular reactor

A 50 MW e helium-cooled fast modular reactor (FMR) is under development, and as part of the Department of Energy Integrated Research Project (IRP), Texas A&M University is conducting the thermal-hydraulic characterization of its baseline core configuration. Here, we experimentally investigated the axial and cross flow fields characteristics in the hexagonal fuel rod bundle composed of 84 rods, a central rod, and spacer grids at a Reynolds number of 12,000. A fully transparent experimental facility resembling of one unit of the fuel assembly was constructed. Time-resolved particle image velocimetry (TR-PIV) measurements were performed to characterize hydraulic behavior downstream the spacer grid. Velocity measurements were conducted in the axial and radial direction of the rod bundle. From the PIV velocity vector fields, the full-field flow statistics were computed for the mean velocity, vorticity, and Reynolds stresses. An analysis of the energy decay downstream of the spacer grid was conducted with calculations of secondary-flow intensities, turbulent kinetic energy, power spectrum, and spatial-temporal velocity cross correlations. The vorticity field was obtained and pairs of counter-rotating vortexes were identified using a 3D reconstruction of several measurement planes. The data from this experimental campaign will be used to validate numerical models based on Computational Fluid Dynamics and other codes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Analytical solution for partial heating on the exterior of the pipe with application to measuring fluid flow rate

An analytical model using Green’s functions for partial external heating of a pipe is developed, which results in an exact mathematical solution for the radial and axial temperature distribution in the pipe wall. Partial heating consists of a constant heat flux function imposed over a small section of the exterior of a pipe and for a limited time duration. The solution comprises steady-state and transient parts, and an algebraic identity is used to decrease the number of summation terms in the slowly-converging steady- state part. Intrinsic verification principles are used to verify the solution. As an example application, this transient solution is applied toward the development of a simple, noninvasive method for in-field measurement of the flow rate in pipes. Furthermore, to simulate this application, a pulse of energy is imposed to the wall of the pipe, and the developed mathematical solution is used to find the flow rate of the fluid inside the pipe. An optimal experiment is designed to find the best measurement location and time.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Prediction of PbLi fluid flow and temperature field in a thermal convection loop for qualification of fusion materials

Two modeling tools, a thermohydraulics code and a computational model in COMSOL Multiphysics, have been developed, tested and then applied to the analysis of fluid flow and heat transfer in a thermal convection loop (TCL). Such a device has recently been used to experimentally evaluate corrosion compatibility of APMT (Fe–21Cr–5Al–3Mo) steel with high-temperature molten eutectic lead-lithium (PbLi) alloy at Oak Ridge National Lab, TN, USA. Here, the 1-D thermohydraulics code allows for rapid calculations of the loop parameters as a function of the applied heat flux. The supplemental COMSOL finite-element computations provide detailed 3-D velocity and temperature field data but are much more time consuming. Both modelling tools demonstrate an excellent agreement in the computed circulation velocity as well as maximum and minimum temperatures. The performed TCL analysis focuses on flow development in the “hot” and “cold” legs, formation of Dean vortices in corner regions, blocking effect of the immersed samples, and radiative and convective cooling effects under the experiment-relevant conditions for Prandtl number Pr = 0.015, Grashof numbers Gr ~10 7 , and Reynolds numbers Re ~10 4 .

36 MATERIALS SCIENCE↗
DOE OSTI2025

Designing slippery rough surfaces to enhance dropwise condensation of low surface tension fluid

Enhancing low surface tension liquid condensation is critical for achieving high energy efficiency and reducing the size of thermal energy systems. Extensive research has focused on promoting dropwise condensation of these liquids using state-of-the-art coatings on plain surfaces. However, maintaining dropwise condensation with low surface tension fluids is challenging due to rivulet formation, resulting in wetted tails that transition to filmwise condensation at elevated heat fluxes. Here, to address this issue, we uncover the role of surface structures and surface chemistry in the dropwise condensation of low surface tension ethanol on slippery rough surfaces (SRS). High-performance dropwise condensation has been achieved on slippery microchannels grafted with perfluoropolyether. The SRS uniquely facilitates rapid lateral droplet removal, enabling faster directional droplet shedding without rivulet formation. The resulting higher droplet removal frequency on SRS leads to heat transfer coefficients 100 % and 500 % higher than conventional dropwise and filmwise condensation on plain surfaces, respectively. Our findings uncover the pivotal role of rapid droplet removal through slippery microchannels in sustaining dropwise condensation of low surface tension liquids. This study introduces a new paradigm for promoting dropwise condensation using engineered SRS that incorporates surface structure and surface chemistry. The work will provide fundamental design guidelines to design efficient and compact condensers that use costly or flammable low global warming potential refrigerants in future refrigeration systems.

Dropwise condensation↗
DOE OSTI2019

Synchrotron radiography characterization of the liquid core dynamics in a canonical two-fluid coaxial atomizer

The liquid core of a canonical two-fluid coaxial atomizer has been characterized using synchrotron X-rays. The high energy photons allow for high-speed imaging of attenuation through the dense liquid-gas jet core, resolving the internal structures that include entrapped air bubbles and the formation of liquid ligaments and bags. When the gas-to-liquid momentum ratio increases, the liquid core transitions from an intact column, where primary break-up happens several liquid diameters downstream, to a hollow crown with a downstream span comparable to the liquid diameter that disintegrates by shedding ligaments from its rim. At high gas momentum ratios (limited by the sonic velocity at the gas nozzle exit), this crown suffers partial dewetting and, when angular momentum is added to the gas, it dewets on a large section of the liquid injection needle circumference. This partial crown exhibits azimuthal motions along the circumference, on timescales much longer than the relevant flow timescales. Here, the crown attachment to the liquid needle presents a bi-stable nature. The dramatic changes of the liquid core morphology, as the gas momentum and swirl ratios vary, have a strong impact on the gas-liquid boundary layers, which control the liquid break-up mechanisms and the resulting spray characteristics, such as droplet size distributions and the droplet volume fraction across the spray.

42 ENGINEERING↗
DOE OSTI2020

Efficient bi-directional coupling of 3D computational fluid-particle dynamics and 1D Multiple Path Particle Dosimetry lung models for multiscale modeling of aerosol dosimetry

The development of predictive aerosol dosimetry models has been a major focus of environmental toxicology and pharmaceutical health research for decades. One-dimensional (1D) models successfully predict overall deposition averages but fail to accurately predict local deposition. Computational fluid-particle dynamics (CFPD) models provide site-specific predictions but at a computational cost that prohibits whole lung predictions. Thus, there is a need for developing multiscale strategies to provide a realistic subject-specific picture of the fate of inhaled aerosol in the lungs. CT-based 3D/CFPD models of the large airways were bidirectionally coupled with individualized 1D Navier-Stokes airflow and particle transport based upon the widely used Multiple Path Particle Dosimetry Model (MPPD). Distribution of airflows among lobes was adjusted by measured lobar volume changes observed in CT images between FRC and FRC + 1.5 L. Additionally, as a test of the effectiveness of the coupling procedures, deposition modeling of previous 1 µm aerosol exposure studies was performed. The complete coupled model was run for 3 breaths, with the computation-intense portion being the 3D CFD Lagrangian particle tracking calculation. The average deposition per breath was 11% in the combined multiscale model with site-specific doses available in the CFPD portion of the model and airway- or region-specific deposition available for the MPPD portion. In conclusion, the key methods developed in this study enable predictions of ventilation heterogeneities and aerosol deposition across the lungs that are not captured by 3D or 1D models alone. Overall, these methods can be used as the foundation for multi-scale modeling of the full respiratory system.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Delivery of Vegetable Oil Suspensions in a Shear Thinning Fluid for Enhanced Bioremediation

In situ anaerobic biological processes are widely applied for dechlorination of chlorinated solvents in groundwater. A wide range of organic substrates have been tested and applied to support the dechlorination processes. Vegetable oils are a promising substrate and have been shown to induce effective dechlorination, have limited geochemical impacts, and good longevity. Distribution of vegetable oil in the subsurface, because it is a non-aqueous phase material, has typically been addressed by creating emulsified oil solutions. In this study, inexpensive waste vegetable oils were suspended in a xanthan gum solution, a shear-thinning fluid, as an alternative oil delivery mechanism. The stability, oil droplet size and distribution, and rheological behavior of the oil suspensions that are created in the xanthan solutions were studied in batch experiments. The injectability of the suspensions and oil distribution in porous medium were evaluated in column tests. Numerical modeling of the oil droplet transport and distribution in porous media was conducted to help interpret the column-test data. Batch studies showed that simple mixing of vegetable oil and xanthan solution produced stable suspensions of the oil as micron-size droplets. The mixture rheology retains shear-thinning properties that facilitate improved uniformity of substrate distribution in heterogeneous aquifers. Column tests demonstrated successful injection of the vegetable oil suspension into porous medium. This study provided evidence that vegetable oil suspensions in xanthan are a potential substrate to support in situ anaerobic bioremediation with favorable injection properties.

Zhong, Lirong↗
DOE OSTI2011

Dimension reduction method for ODE fluid models

We develop a dimension reduction method for large size ODE systems obtained from a dis- cretization of partial differential equations of viscous fluid flow of nearly constant density. The method is also applicable to other large size classical particle systems with negligibly small variations of concentration. We propose a new computational closure for mesoscale balance equations based on numerical iterative deconvolution. To illustrate the computa- tional advantages of the proposed reduction method we use it to solve a system of Smoothed Particle Hydrodynamic ODEs describing Poiseuille flows driven by uniform and periodic (in space) body forces. For the Poiseuille flow driven by the uniform force the coarse solution was obtained with the zero-order deconvolution. For the flow driven by the periodic body force, the first-order deconvolution was necessary to obtain a sufficiently accurate solution.

dimension reduction↗
DOE OSTI2013

Continuum-kinetic-microscopic model of lung clearance due to core-annular fluid entrainment

The human lung is protected against aspirated infectious and toxic agents by a thin liquid layer lining the interior of the airways. This airway surface liquid is a bilayer composed of a viscoelastic mucus layer supported by a fluid film known as the periciliary liquid. The viscoelastic behavior of the mucus layer is principally due to long-chain polymers known as mucins. The airway surface liquid is cleared from the lung by ciliary transport, surface tension gradients, and airflow shear forces. This work presents a multiscale model of the effect of airflow shear forces, as exerted by tidal breathing and cough, upon clearance. The composition of the mucus layer is complex and variable in time. To avoid the restrictions imposed by adopting a viscoelastic flow model of limited validity, a multiscale computational model is introduced in which the continuum-level properties of the airway surface liquid are determined by microscopic simulation of long-chain polymers. A bridge between microscopic and continuum levels is constructed through a kinetic-level probability density function describing polymer chain configurations. The overall multiscale framework is especially suited to biological problems due to the flexibility afforded in specifying microscopic constituents, and examining the effects of various constituents upon overall mucus transport at the continuum scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Stage-parallel fully implicit Runge–Kutta solvers for discontinuous Galerkin fluid simulations

In this paper, we develop new techniques for solving the large, coupled linear systems that arise from fully implicit Runge–Kutta methods. This method makes use of the iterative preconditioned GMRES algorithm for solving the linear systems, which has seen success for fluid flow problems and discontinuous Galerkin discretizations. By transforming the resulting linear system of equations, one can obtain a method which is much less computationally expensive than the untransformed formulation, and which compares competitively with other time-integration schemes, such as diagonally implicit Runge–Kutta (DIRK) methods. We develop and test several ILU-based preconditioners effective for these large systems. We additionally employ a parallel-in-time strategy to compute the Runge–Kutta stages simultaneously. Numerical experiments are performed on the Navier–Stokes equations using Euler vortex and 2D and 3D NACA airfoil test cases in serial and in parallel settings. In conclusion, the fully implicit Radau IIA Runge–Kutta methods compare favorably with equal-order DIRK methods in terms of accuracy, number of GMRES iterations, number of matrix–vector multiplications, and wall-clock time, for a wide range of time steps.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗