Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex 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,603 records · Page 89

Continued Validation Studies using the MOOSE Framework for Plasma Simulation with Electromagnetics

Resources and tools for the modeling and simulation of low-temperature plasma (LTP) discharges are increasingly vital to progress in the field in order to properly characterize and study complex source designs and plasma chemistries beyond the scope of traditional diagnostics. Open-source software provides powerful platforms for this work and can enable community-driven LTP R&D. Within the Multiphysics Object-Oriented Simulation Environment (MOOSE) open-source framework [1], capabilities have been demonstrated in the areas of plasma fluids (Zapdos [2]), plasma chemistry (CRANE [3]), and general electromagnetic wave theory (Electromagnetic Library for Kinetics & fluids [ELK] [4]). ELK has since been coupled to Zapdos/CRANE to enable fully coupled electromagnetic plasma simulations. This talk will detail the continued validation efforts and discuss Zapdos-ELK-CRANE code coupling with various-low temperature plasma sources. The impact of fully coupled electromagnetics on process parameters (e.g., temperature and electron/ion energy) versus an electrostatic description will also be discussed. [1] Permann et al., SoftwareX, 11 (2020) 100430. [2] Lindsay et al., J. Phys. D: Appl. Phys. 49 (2016) 235204. [3] Keniley et al., Talk, GEC Session DT3.00002 (2019). [4] Icenhour et al., Talk, GEC Session RR1.00006 (2019). Acknowledgements: INL Graduate Fellowship Program, NSF SI2 Grant 1740300, U.S. Dept. of Energy Office of Science Graduate Student Research Program

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Synthesis and Single Crystals of Refractory Oxides of Lanthanides and Thorium

At the completion of this program, we can report that we developed a considerable degree of technical improvements in our ability to perform hydrothermal reactions at high temperatures and pressures. We can now routinely perform reactions at 700-750°C and 200 MPa. Currently we are in the process of exploiting this new technology synthesizing a range of exotic new materials investigating relatively poorly understood materials. Our initial efforts focused on the chemistry of rare earth oxides with tetravalent and pentavalent oxides. We recently published a study of the lanthanides with Nb 5+ and Ta 5+ ions, where we grew oxides such as RENdO 4 and RETaO 4 as high quality single crystals. These compounds were targeted as potential hosts for luminescent and scintillation materials, particularly given that they are among the densest oxide hosts and hence have good potential as absorbers for high energy radiation like X-rays and gamma rays. We also isolated a range of unusual new rare earth tantalates with very complex structures. indicating that the chemistry is very sensitive to conditions. We performed some fairly comprehensive examinations of the solid-state chemistry of rare earth ions with various tetravalent metal ions especially Si 4+ , Ge 4+ , Sn 4+ and Ti 4+ . Given the potential role of rare earth silicates in immobilizing radioactive waste elements in long-term storage, and the similarity of our hydrothermal fluids with known geological conditions, this chemistry continues to be relevant. We prepared an extensive series of new lanthanide germanates (e.g. RE 13 Ge 6 O 31 (OH), BaRE 10 (GeO 4 ) 4 O 8 ). and found that there is there is almost no overlap between the chemistry of the rare earth silicates. Stannic oxide (SnO 2 ) is much more refractory and requires higher temperatures and of mineralizer concentrations. One significant result is the growth of RE 2 Sn 2 O 7 pyrochlore single crystals. These are of interest because the rare earth stannate pyrochlores are known to display a wide range of magnetic frustration such as spin ice behavior. We grew high quality single crystals of rare earth germanate and stannate pyrochlores and this led to a collaboration with Professor Kate Ross at Colorado State. Preliminary measurements, indicate that the crystals contain no detectable defects or site disorder. Initial neutron diffraction on single crystals was performed at Oak Ridge, and more detailed experiments involving the Ross group are underway at both NIST and ORNL. This particular chemistry has turned out to be the most potentially significant work on this project and the collaborative effort with Prof. Ross is the topic of a DoE renewal project on quantum materials. Our initial foray into the hydrothermal chemistry of rare earth titanates has also been very promising and a range of cubic and polar ferroic phases of the light rare earths RE 2 Ti 2 O 7 (RE = La - Pr) in the P2 1 phase. We also discovered an interesting new phase Ce 2 Ti 4 O 11 that can have implications in heavy metal immobilization and storage. along with a series of new rare earth titanates (La 5 Ti 4 O 15 (OH) Sm 3 TiO 5 (OH) 3 and Lu 5 Ti 2 O 11 (OH) with exceptionally complex structures. One interesting sidelight has been high temperature hydrothermal chemistry terbium, including the growth of large crystals of TbO(OH). This is not a new compound but it is the first time it has been grown as large single crystals. The Tb atom density is almost as high as that in Tb 2 O 3 and has a very high Verdet constant (ca. 70), making it a very attractive candidate as a Faraday rotator. Unfortunately it is not in a cubic structure but he material is hard, stable, pure and inexpensive, so should still be an attractive Faraday oscillator. We recently received a patent on this material. We also synthesized K 2 Tb(Ge 2 O 7 ) containing stable octahedral Tb 4+ ions, which appears to be the first example of a well-characterized Tb 4+ complex. Given that Tb 4+ has been proposed as a benign surrogate for more treacherous tetravalent ions such as Cf 4+ and Bk 4+ , we think that Tb 4+ silicates can be a particularly useful study for actinide immobilization and related work. We also began reaction studies with rare earths and both ReO 2 and RuO 2 . These resulted in large single crystals of species like RE 5 Ru 2 O 12 , RE 4 Re 2 O 11 , REReO 4 and RE 2 ReO 5 . Several of these samples have already been sent to ORNL for magnetic and neutron diffraction studies.

36 MATERIALS SCIENCE↗

Thermal Modeling and Correlation of the Space Environments Complex Vacuum Chamber and Cryoshroud

From Presentation: Formerly known as the Space Power Facility (SPF) Located at Plum Brook Station in Sandusky, OH -- Houses large-scale test facilities including the Space Simulation Vacuum Chamber -- Largest vacuum chamber in the world -- 100 ft diameter, 122 ft tall -- Cryoshroud operates between -250 °F and +140 °F -- To be used for Orion EM-1 thermal vacuum test -- Thermal Desktop model of vacuum chamber, cryoshroud, mechanical ground support equipment (MGSE) including fluid model of GN2 -- Primary uses: Predict MGSE temperatures; Predict heat load on GN2 system -- Includes 96 TD Temperature Measures that represent test thermocouples -- Steady state correlation with 3 tests: 2 hot; 1 cold -- Two ISTs performed with no MGSE inside the cryoshroud: Insulated test: operated at 170 °F for 37 hours; Uninsulated test: operated at 170 °F for 14 days -- Only cryoshroud thermocouples were installed

Thermal modeling↗

Evaluation and beneficiation of high sulfur and high alkali fly ashes for use as supplementary cementitious materials in concrete

Coal-fuel power plants with semi-dry or dry flue gas desulfurization (FGD) systems produce high sulfur and/or high alkali fly ashes due to comingling of fly ash with FGD products. Such fly ashes do not meet the SO 3 content limit (5.0% max.) of ASTM C618 or are unable to mitigate the alkali-silica reaction. The mineralogy of the sulfur present in these ashes can vary significantly (e.g., CaSO 4 , CaSO 3 , Na2SO 4 ) based on the FGD technology used and this affects the performance of these fly ashes in concrete. Thus, the single SO 3 % limit of ASTM C618 is unable to capture the complexity and performance of fly ash, and this results in elimination of potentially viable pozzolans for concrete. Here this study performs a systematic investigation of the effect of SO 3 type and content in fly ash on various performance parameters of cement-fly ash pastes and mortars, including workability (flow and flow retention), pore fluid pH, setting time, strength development, and potential for deleterious expansion. To better quantify and understand these effects, the study considers both real and doped fly ashes (i.e., a blend of specification-compliant fly ash with target sulfur compounds). The poor performance observed in the case of setting time and pore fluid pH was successfully mitigated using chemical admixtures.

36 MATERIALS SCIENCE↗

The interplay between vapour, liquid, and solid phases in laser powder bed fusion

The capability of producing complex, high performance metal parts on demand has established laser powder bed fusion (LPBF) as a promising additive manufacturing technology, yet deeper understanding of the laser-material interaction is crucial to exploit the potential of the process. By simultaneous in-situ synchrotron x-ray and schlieren imaging, we probe directly the interconnected fluid dynamics of the vapour jet formed by the laser and the depression it produces in the melt pool. The combined imaging shows the formation of a stable plume over stable surface depressions, which becomes chaotic following transition to a full keyhole. We quantify process instability across several parameter sets by analysing keyhole and plume morphologies, and identify a previously unreported threshold of the energy input required for stable line scans. The effect of the powder layer and its impact on process stability is explored. These high-speed visualisations of the fluid mechanics governing LPBF enable us to identify unfavourable process dynamics associated with unwanted porosity, aiding the design of process windows at higher power and speed, and providing the potential for in-process monitoring of process stability.

36 MATERIALS SCIENCE↗

MFIX DEM Enhancement for Industry-Relevant Flows (Final Report)

The overall goal of this two-phase project is to implement performance improvements of the Multiphase Flow with Interphase Exchanges (MFIX) Discrete Element Model (DEM) code that enable a transformative shift for industrial use. Prior to this effort, the largest simulations performed using MFIX are O(10 7 ) particles. This falls short of the O(10 9 ) particle simulations that must be completed on a timescale of days or weeks (vs. months or years) to enable simulations with physically-relevant domain sizes to be incorporated into industrial design cycles within five years. This was accomplished by tailoring best-in-class practices to bear on the unique challenges posed by the MFIX-DEM algorithm and code base. Scientific simulations (e.g., in cosmology, turbulent combustion) routinely use massively parallel computing to update far more particles in short wall clock times. Results from Phase 1 (1.5 years in duration) indicated significant gains in speed were possible for a wide range of benchmark cases. Moreover, a survey sent to >35 companies indicates that the timing is ideal for such an enhanced tool, with >80% of the respondents indicating that DEM is already value-added or will be within the next 5 years, and >70% of the respondents indicating that improved speed is the top computational priority. In Phase 2 (3.5 years in duration), the two major barriers that hinder industry from effectively using multiphase Computational Fluid Dynamics (CFD) to cut costs and improve performance, namely computational overhead and confidence in predictions, continued to be addressed. Regarding the former, the results from Phase 1 to guide the effort, with enhancements focused on an improved time-stepping algorithm and particle sorting. Four target problems of 1 billion particles each and increasing complexity were identified: homogeneous cooling, tumbler with continuous particle size distribution, discharge from a rectangular hopper and a cylindrical riser. Each of these were successfully simulated for relevant time scales (on order of seconds) using less than 24 hours of wall clock time. These represent the first 1-billion particle DEM simulations performed with MFIX, namely using the MFIX-Exa code. This code is currently under development at NETL in collaboration with Lawrence Berkeley National Laboratory. Regarding the second barrier on predictive uncertainty, experiments from Phase 1 (interacting nozzles - hydrodynamics only) and Phase 2 (very small-scale segregation experiments) were used to demonstrate the ability of two simplified approaches to uncertainty quantification (UQ). By limiting the number of particles, UQ based on the simplified treatment was compared to standard UQ, which was shown to have much higher computational demands. Experiments were also performed on a pilot-scale stripper unit to provide validation data for future CFD-DEM simulations and UQ.

20 FOSSIL-FUELED POWER PLANTS↗

The Art of Automation: Translating Electron Microscopy Workflows Into Automated Processes

Acquiring data using a scanning transmission electron microscope (STEM) is a complex, multi-step process. The intricacy of the process depends on the type of sample, composition of the material, desired results of the experiment, resolution requirement and other experimental factors. Each experiment presents unique complications, such as sample drift and contamination, that the microscopist must consider when acquiring data. All these challenges are handled fluidly and expertly by experienced microscopists, but to reach new levels of innovation in material development, including greater reproducibility, throughput, and precision, the automation of these workflows is essential. The initial phase of this work involved translating intuition-based workflows into discrete, programmable steps. Some common key stages in STEM workflows are the initial tuning, scanning the sample for areas of interest, and then acquiring the data. Each stage can be broken further into specific parameter adjustments, such as aberration correction and dwell time optimization, depending on the experiment. When deconstructing various experiments each step was assessed for automation feasibility based on the amount of real time operator decisions. There are steps that lend themselves to automation more readily than others, such as course focusing and sample screening, but there is potential for full automation of all stages with time. As an initial step, an automated montage routine was developed, allowing for the efficient acquisition of large portions of the sample without requiring continuous intervention from the operator. The automation of this small process of the procedure demonstrates the value of this capability. A major challenge in automation arises from discrepancies between commanded, reported and actual stage movements. Using systematic tests, stage movement was quantified. This error can be corrected algorithmically for more accurate workflows in the future. Expanding automation capabilities would result in larger, more efficient data acquisition which allows for more robust statistical analysis. Additionally, this work lays the groundwork for a closed loop system where machine learning algorithms would intake automatically acquired data and make real time decisions. By progressively automating this instrument, this work establishes the foundation for fully automated experimentation in transmission electron microscopy.

97 MATHEMATICS AND COMPUTING↗

Harvesting Energy from Wastewater by Converting Sewage

This project aims were to develop and demonstrate a scalable, integrated process to convert sewage sludge into renewable natural gas (RNG), enabling wastewater treatment plants (WWTPs) to become net energy producers. The system proposal integrates autothermal hydrothermal liquefaction (AT-HTL), supercritical salt precipitation (SCSP), and hydrothermal gasification (HTG), collectively forming the Supercritical Sludge-to-Gas (SC-S2G) platform. Initially, batch hydrothermal liquefaction reactions were used to screen sewage sludge using AT-HTL (later termed RI-HTL) conversion to biocrude, aqueous and char phases compared to hydrothermal liquefaction (HTL). Significant improvement in biocrude yield using peroxide addition at O:C ratio of 0.05 and under conditions of 300°C for 10 minutes gave 57% biocrude yield and 85% fluid carbon yield (biocrude plus aqueous), while minimizing the loss of carbon to char solids (~7%). Hence, RI-HTL was shown to be effective for conversion of real sewage sludge. The corrosion of the alloy reactor tubes or vessels is an important factor when developing a process that includes an oxidant and a chemically complex feed like sewage sludge. We investigated the corrosion rates on metal alloys at 350°C for 240 hours. Corrosion rates of 0.21 and 0.26 mpy for 304L and 316L stainless steel were measured respectively. The corrosion information obtained in this investigation was utilized by PNNL for design, materials sourcing and construction of the pilot scale continuous flow system.

09 BIOMASS FUELS↗

Advances in computational design and analysis of airbreathing propulsion systems

The development of commercial and military aircraft depends, to a large extent, on engine manufacturers being able to achieve significant increases in propulsion capability through improved component aerodynamics, materials, and structures. The recent history of propulsion has been marked by efforts to develop computational techniques that can speed up the propulsion design process and produce superior designs. The availability of powerful supercomputers, such as the NASA Numerical Aerodynamic Simulator, and the potential for even higher performance offered by parallel computer architectures, have opened the door to the use of multi-dimensional simulations to study complex physical phenomena in propulsion systems that have previously defied analysis or experimental observation. An overview of several NASA Lewis research efforts is provided that are contributing toward the long-range goal of a numerical test-cell for the integrated, multidisciplinary design, analysis, and optimization of propulsion systems. Specific examples in Internal Computational Fluid Mechanics, Computational Structural Mechanics, Computational Materials Science, and High Performance Computing are cited and described in terms of current capabilities, technical challenges, and future research directions.

Klineberg, John M.↗

Advances in computational design and analysis of airbreathing propulsion systems

The development of commercial and military aircraft depends, to a large extent, on engine manufacturers being able to achieve significant increases in propulsion capability through improved component aerodynamics, materials, and structures. The recent history of propulsion has been marked by efforts to develop computational techniques that can speed up the propulsion design process and produce superior designs. The availability of powerful supercomputers, such as the NASA Numerical Aerodynamic Simulator, and the potential for even higher performance offered by parallel computer architectures, have opened the door to the use of multi-dimensional simulations to study complex physical phenomena in propulsion systems that have previously defied analysis or experimental observation. An overview of several NASA Lewis research efforts is provided that are contributing toward the long-range goal of a numerical test-cell for the integrated, multidisciplinary design, analysis, and optimization of propulsion systems. Specific examples in Internal Computational Fluid Mechanics, Computational Structural Mechanics, Computational Materials Science, and High Performance Computing are cited and described in terms of current capabilities, technical challenges, and future research directions.

Klineberg, John M.↗

Cytometer on a Chip

A cytometer now under development exploits spatial sorting of sampled cells on a microarray chip followed by use of grating-coupled surface-plasmon-resonance imaging (GCSPRI) to detect the sorted cells. This cytometer on a chip is a prototype of contemplated future miniature cytometers that would be suitable for rapidly identifying pathogens and other cells of interest in both field and laboratory applications and that would be attractive as alternatives to conventional flow cytometers. The basic principle of operation of a conventional flow cytometer requires fluorescent labeling of sampled cells, stringent optical alignment of a laser beam with a narrow orifice, and flow of the cells through the orifice, which is subject to clogging. In contrast, the principle of operation of the present cytometer on a chip does not require fluorescent labeling of cells, stringent optical alignment, or flow through a narrow orifice. The basic principle of operation of the cytometer on a chip also reduces the complexity, mass, and power of the associated laser and detection systems, relative to those needed in conventional flow cytometry. Instead of making cells flow in single file through a narrow flow orifice for sequential interrogation as in conventional flow cytometry, a liquid containing suspended sampled cells is made to flow over the front surface of a microarray chip on which there are many capture spots. Each capture spot is coated with a thin (.50-nm) layer of gold that is, in turn, coated with antibodies that bind to cell-surface molecules characteristic of the cell species of interest. The multiplicity of capture spots makes it possible to perform rapid, massively parallel analysis of a large cell population. The binding of cells to each capture spot gives rise to a minute change in the index of refraction at the surface of the chip. This change in the index of refraction is what is sensed in GCSPRI, as described briefly below. The identities of the various species in a sample of cells is spatially encoded in the chip by the pattern of capture spots. The number of cells of a particular species is determined from the magnitude of the GCSPRI signal from that spot. GCSPRI as used here can be summarized as follows: The cytometer chip is fabricated with a diffraction grating on its front surface. The chip is illuminated with a light emitting diode (LED) from the front. By proper choice of grating parameters and of the wavelength and the angle of incidence of a laser beam, laser light can be made to be coupled into an electromagnetic mode that resonates with surface plasmons and thus couples light into surface plasmons. Coupling of light into a surface plasmon at a given location reduces the amount of incident light reflected from that location. A change in the index of refraction at the surface of a capture spot gives rise to a change in the resonance condition. Depending on the specific design, the change in the index of refraction could manifest itself as a brightening or darkening, a change in the wavelength needed to excite the plasmon at a given angle of incidence, or a change in the angle of incidence needed to excite the plasmon at a given wavelength. Whereas a multiwavelength laser system with multichannel detection would be needed to detect multiple species in conventional flow cytometry, it suffices to use an LED and a single detector channel in the GCSPRI approach: this contributes significantly to reductions in cost, complexity, size, mass, and power. GCSPRI cytometer chips could be made of plastic and could be mass-produced cheaply by use of molding and other methods adopted from the manufacture of digital video disks. These methods are amenable to a high degree of miniaturization: such additional features as fluidic channels, reaction chambers, and fluid-coupling ports could readily be incorporated into the chips, without incurring substantial additional costs.

Fernandez, Salvador M.↗

Cytometer on a Chip

A cytometer now under development exploits spatial sorting of sampled cells on a microarray chip followed by use of grating-coupled surface-plasmon-resonance imaging (GCSPRI) to detect the sorted cells. This cytometer on a chip is a prototype of contemplated future miniature cytometers that would be suitable for rapidly identifying pathogens and other cells of interest in both field and laboratory applications and that would be attractive as alternatives to conventional flow cytometers. The basic principle of operation of a conventional flow cytometer requires fluorescent labeling of sampled cells, stringent optical alignment of a laser beam with a narrow orifice, and flow of the cells through the orifice, which is subject to clogging. In contrast, the principle of operation of the present cytometer on a chip does not require fluorescent labeling of cells, stringent optical alignment, or flow through a narrow orifice. The basic principle of operation of the cytometer on a chip also reduces the complexity, mass, and power of the associated laser and detection systems, relative to those needed in conventional flow cytometry. Instead of making cells flow in single file through a narrow flow orifice for sequential interrogation as in conventional flow cytometry, a liquid containing suspended sampled cells is made to flow over the front surface of a microarray chip on which there are many capture spots. Each capture spot is coated with a thin (approximately 50-nm) layer of gold that is, in turn, coated with antibodies that bind to cell-surface molecules characteristic of one the cell species of interest. The multiplicity of capture spots makes it possible to perform rapid, massively parallel analysis of a large cell population. The binding of cells to each capture spot gives rise to a minute change in the index of refraction at the surface of the chip. This change in the index of refraction is what is sensed in GCSPRI, as described briefly below. The identities of the various species in a sample of cells is spatially encoded in the chip by the pattern of capture spots. The number of cells of a particular species is determined from the magnitude of the GCSPRI signal from that spot. GCSPRI as used here can be summarized as follows: The cytometer chip is fabricated with a diffraction grating on its front surface. The chip is illuminated with a light emitting diode (LED) from the front. By proper choice of grating parameters and of the wavelength and the angle of incidence of a laser beam, laser light can be made to be coupled into an electromagnetic mode that resonates with surface plasmons and thus couples light into surface plasmons. Coupling of light into a surface plasmon at a given location reduces the amount of incident light reflected from that location. A change in the index of refraction at the surface of a capture spot gives rise to a change in the resonance condition. Depending on the specific design, the change in the index of refraction could manifest itself as a brightening or darkening, a change in the wavelength needed to excite the plasmon at a given angle of incidence, or a change in the angle of incidence needed to excite the plasmon at a given wavelength. Whereas a multiwavelength laser system with multichannel detection would be needed to detect multiple species in conventional flow cytometry, it suffices to use an LED and a single detector channel in the GCSPRI approach: this contributes significantly to reductions in cost, complexity, size, mass, and power. GCSPRI cytometer chips could be made of plastic and could be mass-produced cheaply by use of molding and other methods adopted from the manufacture of digital video disks. These methods are amenable to a high degree of miniaturization: such additional features as fluidic channels, reaction chambers, and fluid-coupling ports could readily be incorporated into the chips, without incurring substantial additional costs.

Fernandez, Salvador M.↗

Probabilistic Modeling Of Ocular Biomechanics In VIIP: Risk Stratification

Visual Impairment and Intracranial Pressure (VIIP) syndrome is a major health concern for long-duration space missions. Currently, it is thought that a cephalad fluid shift in microgravity causes elevated intracranial pressure (ICP) that is transmitted along the optic nerve sheath (ONS). We hypothesize that this in turn leads to alteration and remodeling of connective tissue in the posterior eye which impacts vision. Finite element (FE) analysis is a powerful tool for examining the effects of mechanical loads in complex geometries. Our goal is to build a FE analysis framework to understand the response of the lamina cribrosa and optic nerve head to elevations in ICP in VIIP. To simulate the effects of different pressures on tissues in the posterior eye, we developed a geometric model of the posterior eye and optic nerve sheath and used a Latin hypercubepartial rank correlation coef-ficient (LHSPRCC) approach to assess the influence of uncertainty in our input parameters (i.e. pressures and material properties) on the peak strains within the retina, lamina cribrosa and optic nerve. The LHSPRCC approach was repeated for three relevant ICP ranges, corresponding to upright and supine posture on earth, and microgravity [1]. At each ICP condition we used intraocular pressure (IOP) and mean arterial pressure (MAP) measurements of in-flight astronauts provided by Lifetime Surveillance of Astronaut Health Program, NASA Johnson Space Center. The lamina cribrosa, optic nerve, retinal vessel and retina were modeled as linear-elastic materials, while other tissues were modeled as a Mooney-Rivlin solid (representing ground substance, stiffness parameter c1) with embedded collagen fibers (stiffness parameters c3, c4 and c5). Geometry creationmesh generation was done in Gmsh [2], while FEBio was used for all FE simulations [3]. The LHSPRCC approach resulted in correlation coefficients in the range of 1. To assess the relative influence of the uncertainty in an input parameter on the peak strains, we ranked and then normalized these coefficients, considering that normalized values 0.5 implied a substantial influence on the range of the peak strains in the optic nerve head (ONH). IOP and ICP were found to have a major influence on the peak strains in the ONH, as did optic nerve and LC stiffness. Interestingly, the stiffness of the sclera far from the scleral canal did not have a large influence on peak strains in ONH tissues; however, the collagen fiber stiffness in the peripapillary sclera and annular ring both influenced the peak strains within the ONH. We have created a physiologically relevant model that incorporated collagen fibers to study the effects of elevated ICP. Elevated ICP resulted in strains in the optic nerve that are not predicted to occur on earth: the upright or supine conditions. We found that IOP, ICP, lamina cribrosa stiffness and optic nerve stiffness had the highest association with these extreme strains in the ONH. These extreme strains may activate mechanosensitive cells that induce tissue remodeling and are a risk factor for the development of VIIP.

biomechanics↗

Geothermal Heat Recovery Complex: Large-Scale, Deep Direct-Use System in a Low-Temperature Sedimentary Basin (Final Report)

A feasibility study of using deep direct-use (DDU) geothermal energy to heat agricultural research facilities (ARFs) was conducted at the University of Illinois at Urbana-Champaign (U of IL) and its similar application to military facilities in the Illinois Basin (ILB). The geothermal energy system (GES) investigated utilizes low-temperature (30–90°C; 90–190°F) geothermal fluid (i.e., brine) from an extraction well that is part of a deep, two-well (doublet) system that extends to the bottom of the ILB. The geothermal reservoir modeled, the Mt. Simon Sandstone (MSS), is about 1,280 m (4,200 feet) deep and 457 m thick (1,500 feet) beneath the U of IL. The DDU GES surface infrastructure includes heat exchangers connected in-parallel to pipelines carrying the geothermal fluid and fresh cold and hot water. Analysis of the GES indicated that the MSS can provide a baseload of 2 MMBtu/hr to heat the ARFs by extracting 954 m³/d (6,000 barrels/day [bbl/d]) of geothermal fluid that has a temperature of 44–46 °C (111–115 °F). In addition to analyzing the levelized cost of heat (LCOH) and life cycle costs, the environmental effects of the DDU GES were evaluated, including reduced greenhouse gas (GHG) emissions and water consumption. Multiple system designs were evaluated and then ranked based on their maximum heating performance, energy efficiency, and cost recovery. This feasibility study identified the key components of the fully-integrated DDU technology that can be implemented, both technically and economically. The results and information from this study provides end-users and policy makers with guidance for additional research on the specific components of DDU technology such that its widespread use can provide an uninterruptible energy source, increase resilience from extreme weather conditions, reduce U.S. dependency on fossil fuels, and reduce greenhouse gas (GHG) emissions. The site-specific part of this study gives U of IL administrators a realistic and pragmatic assessment of the financial resources necessary to add a DDU GES in the MSS to the campus’ energy portfolio.

15 GEOTHERMAL ENERGY↗

Grid generation for a complex aircraft configuration

The procedure used to create a grid around the F/A-18 fighter aircraft is presented. This work was done for the NASA High Alpha Technology Program. As part of this program, LeRC is numerically and experimentally investigating the flow in the F/A-18 inlet duct at high angles of attack. A grid was needed which could be used to calculate both the external and internal flow around the F/A-18. The grid had to be compatible with the computational fluid dynamics (CFD) codes PARC3D and CFL3D. The programs used to create this grid were I3GVIRGO and GRIDGEN. A surface definition used to create the grid was obtained from McDonnell Aircraft Company (MCAIR) and was composed of numerous files each containing a point definition of a portion of the aircraft. These files were read into the geometry manipulation program I3GVIRGO, where they were modified and grouped into smaller GRIDGEN database files. Next, the block outlines and boundary conditions were specified in the GRIDBLOCK program. The GRIDGEN2D program was used to create the surface grid on the block faces, and GRIDGEN3D was used to create the full 3-D grid.

Bruns, Jim↗

Near-continuum, hypersonic oxygen flow over a double cone simulated by direct simulation Monte Carlo informed from quantum chemistry

A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.

Mechanics↗

Terrain-Influenced Winds and Fire-Fire Interactions in Wildland Fire Simulations [Dissertation]

Ensemble-based approaches to prescribed fire planning cannot be supported by computational fluid dynamics based models like FIRETEC and the Wildland-Urban Interface Fire Dynamics Simulator (WFDS) because they are too computationally expensive and cannot leverage large eddy simulation approaches like CAWFE and WRF-SFIRE because they have too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, that is typical for instance in the Western United States. This dissertation describes a variety of improvements made to QUIC-Fire and its various incorporated algorithms in an effort to make it a viable tool in simulating wildland and prescribed fires on terrain. The modifications made to QUIC-Fire are described in three chapters. The first chapter describes the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. The terraininfluenced winds it generates are analyzed and compared. In particular, this chapter presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captures the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. The model agrees well with measurements on the upwind side of the peak, but overestimates speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field is obtained in 6 s, making the solver appealing for coupled fire-atmosphere simulations. The Socorro Mountain is highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect is addressed in this work). In the second chapter the implementation of the terrain-following version of QUIC-URB into QUIC-Fire, and the necessary changes needed to include terrain are described. No changes to the underlying fire spread algorithm are made other than what is required to correctly account for the inclusion of terrain. Previously published FIRETEC results that use five different topographies that share the same centerline profile are compared to simulation results from the modified QUIC-Fire that use the same topographies and fuels. QUIC-Fire results show overall similar behaviors in terms of how the topographies affect fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all topographies are over-predicted when compared to FIRETEC. Lateral fire growth shows similar trends with FIRETEC between topographies but does not capture the increase in spread due to a diagonal interface between grassland and forested fuel region of the domain. These results suggest that there are three algorithms within QUIC-Fire that could use improvement: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB. Lastly, the third chapter describes a modification to the initial guess used for the calculation of the QUIC-URB mass-conserved wind solution during fire simulations. The modification is aimed at improving fire-fire interactions in QUIC-Fire simulations. The modification consists of using the solution from the previous timestep as the starting point for the calculation of the solution for the next timestep. Fire-fire interactions is greatly improved by the change but a new source of error is introduced. Due to how plumes are modelled in QUIC-Fire the new solution contains errors where gaps in the plume structure are present. However, these errors are mostly limited to the upper atmosphere, where they do not affect fire behavior at the surface, and their magnitude isn’t significant enough to discount the amount of new fire phenomenology now captured in QUIC-Fire with the change.

58 GEOSCIENCES↗

Numerical simulation using vorticity-vector potential formulation

An accurate and efficient computational method is needed for three-dimensional incompressible viscous flows in engineering applications. On solving the turbulent shear flows directly or using the subgrid scale model, it is indispensable to resolve the small scale fluid motions as well as the large scale motions. From this point of view, the pseudo-spectral method is used so far as the computational method. However, the finite difference or the finite element methods are widely applied for computing the flow with practical importance since these methods are easily applied to the flows with complex geometric configurations. However, there exist several problems in applying the finite difference method to direct and large eddy simulations. Accuracy is one of most important problems. This point was already addressed by the present author on the direct simulations on the instability of the plane Poiseuille flow and also on the transition to turbulence. In order to obtain high efficiency, the multi-grid Poisson solver is combined with the higher-order, accurate finite difference method. The formulation method is also one of the most important problems in applying the finite difference method to the incompressible turbulent flows. The three-dimensional Navier-Stokes equations have been solved so far in the primitive variables formulation. One of the major difficulties of this method is the rigorous satisfaction of the equation of continuity. In general, the staggered grid is used for the satisfaction of the solenoidal condition for the velocity field at the wall boundary. However, the velocity field satisfies the equation of continuity automatically in the vorticity-vector potential formulation. From this point of view, the vorticity-vector potential method was extended to the generalized coordinate system. In the present article, we adopt the vorticity-vector potential formulation, the generalized coordinate system, and the 4th-order accurate difference method as the computational method. We present the computational method and apply the present method to computations of flows in a square cavity at large Reynolds number in order to investigate its effectiveness.

Tokunaga, Hiroshi↗