Search NASA⌕ Search

SEARCH · Search NASA

Results for “splitting methods”

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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

On split inclusion problem and fixed point problem for multi-valued mappings

Many of the iterative schemes for solving split inclusion and fixed point problems involve step-sizes that depend on the norm of a bounded linear operator. The implementation of such algorithms are usually difficult to handle. This is because they require the computation of the operator norm. In this paper, we propose an algorithm involving a step-size selected in such a way that its implementation does not require the computation or an estimate of some spectral radius. Using our algorithm we proved strong convergence theorem for split inclusion problem and fixed point problem for multi-valued quasi-nonexpansive mappings in real Hilbert spaces. Our result generalizes some important and recent results in the literature. Some applications of our main result to game theory and variational inequality problem are also presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Crystal growth and quantum oscillations in the topological chiral semimetal CoSi

Here, we survey the electrical transport properties of single-crystalline, topological chiral semimetal CoSi grown via different methods. High-quality CoSi single crystals were found in the growth from a tellurium solution. The sample's high carrier mobility enables us to observe quantum oscillations (QOs) in its thermoelectrical signals. Our analysis of QOs reveals two spherical Fermi surfaces around the R point in its Brillouin zone corner. The extracted Berry phases of these electron orbits are consistent with the -2 chiral charge as reported in density functional theory (DFT) calculations. A detailed analysis of the QOs reveals that the spin-orbit-coupling-induced band splitting is less than 2 meV near the Fermi level, one order of magnitude smaller than our DFT calculation. We also report a large phonon-drag-induced Nernst effect in CoSi at intermediate temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Orbital-selective Mott phase and non-Fermi liquid in FePS 3

The layered metal phosphorous trisulfide FePS 3 is reported to be a Mott insulator at ambient conditions and to undergo structural and insulator-metal phase transitions under pressure. However, the character of the resulting metallic states has not been understood clearly so far. Here, we theoretically study the phase transitions of FePS 3 using first-principles methods based on density functional theory and embedded dynamical mean field theory. We find that the Mott transition in FePS 3 can be orbital selective, with t 2g states undergoing a correlation-induced insulator-to-metal transition while e g states remain gapped. We show that this orbital-selective Mott phase, which occurs only when nonhydrostatic pressure is used, is a bad metal (or non-Fermi liquid) with large fluctuating moments due to Hund's coupling. Further application of pressure increases the crystal-field splitting and converts the system to a conventional Fermi liquid with low-spin configurations dominant. Finally, our results show that FePS 3 is an example of a system that realizes an orbital-selective Mott phase, allowing tuning between correlated and uncorrelated metallic properties in an accessible pressure range (≤ 18 GPa).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1991

A multi-block multigrid method for the solution of the Euler and Navier-Stokes equations for three-dimensional flows

A multi-block multigrid method for the solution of the three-dimensional Euler and Navier-Stokes equations is presented. The basic flow solver is a cell-vertex method which employs central-difference spatial approximations and Runge-Kutta time stepping. The use of local time stepping, implicit residual smoothing, multigrid techniques, and variable-coefficient numerical smoothing results in an efficient and robust scheme. The multi-block strategy places the block loop within the Runge-Kutta loop such that accuracy and convergence are not affected by block boundaries. This has been verified by comparing the results of one- and two-block calculations in which the two-block grid is generated by splitting the one-block grid. Results are presented for both Euler and Navier-Stokes computations of wings and wing-fuselage combinations.

Atkins, H. L.↗
NASA NTRS1994

Idealized textile composites for experimental/analytical correlation

Textile composites are fiber reinforced materials produced by weaving, braiding, knitting, or stitching. These materials offer possible reductions in manufacturing costs compared to conventional laminated composites. Thus, they are attractive candidate materials for aircraft structures. To date, numerous experimental studies have been performed to characterize the mechanical performance of specific textile architectures. Since many materials and architectures are of interest, there is a need for analytical models to predict the mechanical properties of a specific textile composite material. Models of varying sophistication have been proposed based on mechanics of materials, classical laminated plate theory, and the finite element method. These modeling approaches assume an idealized textile architecture and generally consider a single unit cell. Due to randomness of the textile architectures produced using conventional processing techniques, experimental data obtained has been of limited use for verifying the accuracy of these analytical approaches. This research is focused on fabricating woven textile composites with highly aligned and accurately placed fiber tows that closely represent the idealized architectures assumed in analytical models. These idealized textile composites have been fabricated with three types of layer nesting configurations: stacked, diagonal, and split-span. Compression testing results have identified strength variations as a function of nesting. Moire interferometry experiments are being used to determine localized deformations for detailed correlation with model predictions.

Adams, Daniel O.↗
DOE OSTI2024

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

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

Optimizing accuracy and efficacy in data-driven materials discovery for the solar production of hydrogen

The production of hydrogen fuels, via water splitting, is of practical relevance for meeting global energy needs and mitigating the environmental consequences of fossil-fuel-based transportation. Water photoelectrolysis has been proposed as a viable approach for generating hydrogen, provided that stable and inexpensive photocatalysts with conversion efficiencies over 10% can be discovered, synthesized at scale, and successfully deployed. While a number of first-principles studies have focused on the data-driven discovery of photocatalysts, in the absence of systematic experimental validation, the success rate of these predictions may be limited. We address this problem by developing a screening procedure with co-validation between experiment and theory to expedite the synthesis, characterization, and testing of the computationally predicted, most desirable materials. Here, starting with 70 150 compounds in the Materials Project database, the proposed protocol yielded 71 candidate photocatalysts, 11 of which were synthesized as single-phase materials. Experiments confirmed hydrogen generation and favorable band alignment for 6 of the 11 compounds, with the most promising ones belonging to the families of alkali and alkaline-earth indates and orthoplumbates. This study shows the accuracy of a nonempirical, Hubbard-corrected density-functional theory method to predict band gaps and band offsets at a fraction of the computational cost of hybrid functionals, and outlines an effective strategy to identify photocatalysts for solar hydrogen generation.

08 HYDROGEN↗
DOE OSTI2018

Generalized, Energy-conserving Numerical Simulations of Particles in General Relativity. I. Time-like and Null Geodesics

The numerical integration of particle trajectories in curved spacetimes is fundamental for obtaining realistic models of the particle dynamics around massive compact objects such as black holes and neutron stars. Generalized algorithms capable of handling generic metrics are needed for studies of both standard (Schwarzschild and Kerr metrics) and nonstandard (e.g., Schwarzschild metric plus nonclassical perturbations or multiple black hole metrics) spacetimes. The most commonly employed explicit numerical schemes (e.g., Runge–Kutta) are incapable of producing highly accurate results at critical points, e.g., in the regions close to the event horizon where gravity causes extreme curvature of the spacetime, at an acceptable computational cost. In this work we describe a generalized algorithm for the numerical integration of time-like (massive particles) and null (photons) geodesics in any given 3 + 1 split spacetime. We introduce a new, exactly energy-conserving implicit integration scheme based on the preservation of the underlying Hamiltonian, and we compare its properties with a standard fourth-order Runge–Kutta explicit scheme and an implicit midpoint scheme. We test the numerical performance of the three schemes against analytic solutions of particle and photon orbits in Schwarzschild and Kerr spacetimes. We additionally prove the versatility of our framework in handling more exotic metrics, such as Morris–Thorne wormholes and quantum-perturbed Schwarzschild black holes. The generalized method is also discussed in the perspective of future extensions to more complex particle dynamics, e.g., the addition of the Lorentz force acting on charged particles, which allows for test particle diagnostics in GRMHD simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2005

Turbulent Vortex-Flow Simulation Over a 65 deg Sharp and Blunt Leading-Edge Delta Wing at Subsonic Speeds

Turbulent thin-layer, Reynolds-Averaged Navier-Stokes solutions, based on a multi-block structured grid, are presented for a 65 deg delta wing having either a sharp leading edge (SLE) or blunt leading edge (BLE) geometry. The primary objective of the study is to assess the prediction capability of the method for simulating the leading-edge flow separation and the ensuing vortex flow characteristics. Computational results are obtained for two angles of attack of approximately 13 and 20 deg, at free-stream Mach number of 0.40 and Reynolds number of 6 million based on the wing mean aerodynamic chord. The effects of two turbulence models of Baldwin-Lomax with Degani-Schiff (BL/DS) and the Spalart-Allmaras (SA) on the numerical results are also discussed. The computations also explore the effects of two numerical flux-splitting schemes, i.e., flux difference splitting (fds) and flux vector splitting (fvs), on the solution development and convergence characteristics. The resulting trends in solution sensitivity to grid resolution for the selected leading-edge geometries, angles of attack, turbulence models and flux splitting schemes are also presented. The validity of the numerical results is evaluated against a unique set of experimental wind-tunnel data that was obtained in the National Transonic Facility at the NASA Langley Research Center.

Ghaffari, Farhad↗
DOE OSTI2023

Plasma assisted distributed chemical production

The present development is a process to produce commodity chemicals such as methanol and syngas using an integrated plasma catalysis technology. The method comprises providing a fixed or fluidized bed reactor having a microwave plasma flame and a catalyst bed with a catalyst, wherein the catalyst is an alloyed bimetallic nanowire. In the process, the plasma flame fluidizes the catalyst thereby producing a more effective catalyst than the non-fluidized catalyst. It is anticipated that the reactor can have a throughput capacity of up to 30 Lpm/kW and can be effective for the conversion of CO 2 , CH 4 , air, water, and combinations thereof, through reactions such as pure CO 2 splitting, reverse water gas shift (RWGS) for CO production, methanol synthesis, and plasma reforming of methane, thereby making a system that would be attractive for small GTL units.

Sunkara, Mahendra↗
DOE OSTI2020

Impact of Laboratory-Accelerated Aging Methods to Study Alkali–Silica Reaction and Reinforcement Corrosion on the Properties of Concrete

This study focuses on two separate investigations of the main aging mechanisms: alkali–silica reactivity (ASR) and the corrosion of reinforcing steel (rebar) concrete, both of which may result in a premature failure to meet the serviceability or strength requirements of a concrete structure. However, these processes occur very slowly, spanning decades. The impact of direct chemical additives to fresh concrete to accelerate ASR and the corrosion of reinforcing steel on the fresh and hardened properties of the ensuing material are investigated to inform the potential use of chemicals in large-scale studies. The deterioration of reinforced concrete (RC) is determined by means of expansion, cracking, bulk diffusivity and surface resistivity measurements, and compressive, split tensile and flexural strength tests. The results indicate that the addition of sodium hydroxide and calcium chloride can effectively accelerate the crack formation and propagation in concrete due to ASR and the corrosion of rebar, respectively. The ASR-induced cracks maintained a constant crack width from 0.05 mm to 0.1 mm over the measurement period regardless of the intensity of aging acceleration. Adding 4% chloride by weight of cement for accelerating rebar corrosion resulted in an average crack that was 82% larger than in the case of ASR accelerated with the addition of sodium hydroxide. The addition of alkali resulted in an increase in early-age (7-day) strength. At a total alkali loading of 2.98 kg/m3, 3.84 kg/m3 and 5.57 kg/m3, the 28-day compressive strength of concrete decreased by 3%, 10% and 24%, respectively. Similarly, a higher early-age strength and a lower later-age strength was observed for the concrete in the presence of corrosive calcium chloride. The results from this research are expected to inform future studies on the long-term performance of RC structures under accelerated ASR and corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An investigation of the multi-mode Richtmyer-Meshkov instability at a gas/HE interface using Pagosa

Here, in this work, we present a hydrocode Pagosa and explore the Richtmyer-Meshkov Instability (RMI) at an air/high explosive (HE) interface for the first time that is important but has not received much attention yet in the high explosive safety field. Thus, the presented Pagosa can be expected to predict the whole deflagration-to-detonation transition (DDT) process in future. In Pagosa, spatial discretization is implemented on cubic staggered grids by computing different variables at the vertex and the cell center, respectively, a special operator-splitting technique is employed to reduce the computational cost, and an artificial viscosity is added to handle the discontinuous shock waves in our simulations. To quantitatively evaluate the capability of Pagosa to solve these kinds of instabilities, the single mode Rayleigh-Taylor instability (RTI) and the multimode RMI at an air/SF 6 interface are first performed, respectively. The Pagosa results are compared with the related numerical solutions in the existing references and the experimental result. Moreover, a theoretical derivation of growth of RTI is also provided based on our numerical method. Subsequently, we explore the multi-mode RMI at an air/HE interface as well as the effects of several factors using Pagosa. Numerical results show that Pagosa is a powerful toolset to generate the right structures and the amplitude of RTI and RMI at an air/SF 6 interface. The solid HE can be penetrated by a strong shock wave and forms RMI deformations. The RMI at an air/HE interface behaves very different than at an air/SF 6 interface, periodic, decreased oscillation is observed due to material character, and is very sensitive to the initial simulation settings, that is, a tiny change in physical quantities will lead to a remarkable RMI structure, which is also observed in a shock bubble interaction. The findings in this work are significant, and will present a new insight for the high explosive field.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Origin of the magnetic and orbital ordering in α-Sr 2 CrO 4

Motivated by recent experimental progress in transition metal oxides with the K 2 NiF 4 structure, in this study we investigate the magnetic and orbital ordering in α-Sr 2 CrO 4 . Using first principles calculations, first we derive a three-orbital Hubbard model, which reproduces the ab initio band structure near the Fermi level. The unique reverse splitting of t 2g orbitals in α-Sr 2 CrO 4 , with the 3d 2 electronic configuration for the Cr 4+ oxidation state, opens up the possibility of orbital ordering in this material. Using real-space Hartree-Fock for multi-orbital systems, we constructed the ground state phase diagram for the two dimensional compound α-Sr 2 CrO 4 . We found stable ferromagnetic, antiferromagnetic, antiferro-orbital, and staggered orbital stripe ordering in robust regions of the phase diagram. Furthermore, using the density matrix renormalization group method for two-leg ladders with the realistic hopping parameters of α-Sr 2 CrO 4 , we explore magnetic and orbital ordering for experimentally relevant interaction parameters. Again, we find a clear signature of antiferromagnetic spin ordering along with antiferro-orbital ordering at moderate to large Hubbard interaction strength. We also explore the orbital-resolved density of states with Lanczos, predicting insulating behavior for the compound α-Sr 2 CrO 4 , in agreement with experiments. Finally, an intuitive understanding of the results is provided based on a hierarchy between orbitals, with d xy driving the spin order, while electronic repulsion and the effective one dimensionality of the movement within the d xz and d yz orbitals driving the orbital order.

36 MATERIALS SCIENCE↗
NASA NTRS2021

Synthesis and Characterization of Fe(III)-Fe(II)- Mg- Al Smectite Solid Solutions and Implications for Planetary Science

This study demonstrates the synergies and limits of multiple measurement types for the detection of smectite chemistry and oxidation state on planetary surfaces to infer past geochemical conditions. Smectite clay minerals are common products of water-rock interactions throughout the solar system, and their detection and characterization provides important clues about geochemical conditions and past environments if sufficient information about their composition can be discerned. Here, we synthesize and report on the spectroscopic properties of a suite of smectite samples that span the intermediate compositional range between Fe(II), Fe(III), Mg, and Al end-member species using bulk chemical analyses, X-ray diffraction, Vis/IR reflectance spectroscopy, UV and green-laser Raman spectroscopy, and Mössbauer spectroscopy. Our data show that smectite composition and the oxidation state of octahedral Fe can be reliably identified in the near infrared on the basis of combination and fundamental metal-OH stretching modes between 2.1–2.9 μm, which vary systematically with chemistry. Smectites dominated by Mg or Fe(III) have spectrally distinct fundamental and combination stretches, whereas Al-rich and Fe(II)-rich smectites have similar fundamental minima near 2.76 μm, but have distinct combination M-OH features between 2.24 and 2.36 μm. We show that with expanded spectral libraries that include intermediate composition smectites and both Fe(III) and Fe(II) oxidation states, more refined characterization of smectites from MIR data is now possible, as the position of the 450 cm–1 absorption shifts systematically with octahedral Fe content, although detailed analysis is best accomplished in concert with other characterization methods. Our data also provide the first Raman spectral libraries of smectite clays as a function of chemistry, and we demonstrate that Raman spectroscopy at multiple excitation wavelengths can qualitatively distinguish smectite clays of different structures and can enhance interpretation by other types of analyses. Our sample set demonstrates how X-ray diffraction can distinguish between dioctahedral and trioctahedral smectites using either the (02,11) or (06,33) peaks, but auxiliary information about chemistry and oxidation state aids in specific identifications. Finally, the temperature-dependent isomer shift and quadrupole splitting in Mössbauer data are insensitive to changes in Fe content but reliability differentiates Fe within the smectite mineral structure.

Valerie K Fox↗
DOE OSTI2022

Ten Years of VIIRS Land Surface Temperature Product Validation

The Visible Infrared Imaging Radiometer Suite (VIIRS) Land Surface Temperature (LST) has been operationally produced for a decade since the Suomi National Polar-orbiting Partnership (SNPP) launched in October 2011. A comprehensive evaluation of its accuracy and precision will be helpful for product users in climate studies and atmospheric models. In this study, the VIIRS LST is validated with ground observations from multiple high-quality radiation networks, including six stations from the Surface Radiation budget (SURFRAD) network, two stations from the Baseline Surface Radiation Network (BSRN), and 13 stations from the Atmospheric Radiation Measurement (ARM) network, to evaluate its performance over various land-cover types. The VNP21A1 LST was validated against the same ground observations as a reference. The results yield a close agreement between the SNPP VIIRS LST and ground LSTs with a bias of -0.4 K and a RMSE of 1.96 K over six SURFRAD sites; a bias of -0.2 K and a RMSE of 1.93 K over two BSRN sites; and a bias of -0.1 K and a RMSE of 1.7 K over the 13 ARM sites. The time series of the LST errors over individual sites indicate seasonal cycles. The data anomaly over the BSRN site in Cabauw and the SURFRAD site in Desert Rock is revealed and discussed in this study. In addition, a method using Landsat-8 data is applied to quantify the heterogeneity level of each ground station and the results provide promising insights. The validation results demonstrate the maturity of the JPSS VIIRS LST products and their readiness for various application studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Highly Ordered Mesoporous Bimetallic Phosphides as Efficient Oxygen Evolution Electrocatalysts

Oxygen evolution from water using earth abundant transition metal based catalysts is of importance for the commercialization of water electrolyzer. Herein, we report a novel hard templating method to synthesize transition metal phosphides with uniform shape and size. By virtue of the structural feature, synergistic effect among metals and the in situ formed active species, the as-prepared phosphides with optimized composition present enhanced electrocatalytic performance towards oxygen evolution reaction in alkaline solution. In details, the most efficient catalyst reaches a current density of 10 mA/cm2 at a potential of 1.511 V vs reversible hydrogen electrode, which is much lower than that of commercial RuO2 catalyst. Our work offers a new strategy to optimize the catalysts for water splitting by controlling the morphology and composition.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2017

Mesoporous Iron Sulfide for Highly Efficient Electrocatalytic Hydrogen Evolution

Here, we report a facile synthetic protocol to pre-pare mesoporous FeS 2 without the aid of hard template as an electrocatalyst for the hydrogen evolution reaction (HER). The mesoporous FeS 2 materials with high surface area were successfully prepared by a sol-gel method follow-ing a sulfurization treatment in an H 2 S atmosphere. A re-markable HER catalytic performance was achieved with a low overpotential of 96 mV at a current density of 10 mA·cm 2 and a Tafel slope of 78 mV per decade under alka-line conditions (pH 13). These theoretical calculations indicate that the excellent catalytic activity of mesoporous FeS 2 is attributed to the exposed (210) facets. The mesoporous FeS 2 material might be a promising alternative to the Pt-based electrocatalysts for water splitting.

36 MATERIALS SCIENCE↗
NASA NTRS1999

A Pseubo-Temporal Multi-Grid Relaxation Scheme for Solving the Parabolized Navier-Stokes Equations

A multi-grid, flux-difference-split, finite-volume code, VULCAN, is presented for solving the elliptic and parabolized form of the equations governing three-dimensional, turbulent, calorically perfect and non-equilibrium chemically reacting flows. The space marching algorithms developed to improve convergence rate and or reduce computational cost are emphasized. The algorithms presented are extensions to the class of implicit pseudo-time iterative, upwind space-marching schemes. A full approximate storage, full multi-grid scheme is also described which is used to accelerate the convergence of a Gauss-Seidel relaxation method. The multi-grid algorithm is shown to significantly improve convergence on high aspect ratio grids.

Morrison, J. H.↗