Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Complex media”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Strontium Speciation in Relevant Tank Waste Components Examined by Electrospray Ionization Mass Spectrometry

The identification of chemical species formed in complex nuclear waste is crucial for the development and employment of advanced separations technologies to remediate the Hanford site by processing tank waste. The current Tank Side Cesium Removal (TSCR) process deployed at Hanford utilizes crystalline silicotitanate (CST) ion exchange (IX) media to aid in the separation of low-activity waste for proper treatment and disposal. The inorganic IX media is highly selective for Cs but has been shown to also remove Sr from caustic simulants and small-scale IX processing of Hanford tank waste.(Fiskum, Rovira et al. 2019, Fiskum, Campbell et al. 2021, Westesen, Campbell et al. 2022) Quantitative Sr removal has not been observed in all tank waste supernates tested; thus, to better understand Sr removal and effectively predict processing behavior through TSCR, it is necessary to first investigate Sr speciation in tank waste. This work utilized electrospray ionization mass spectrometry (ESI-MS) to identify ionic Sr complexes that form in the presence of NO 3 –, NO 2 –, OH–, and Cl–. Although our results show that NO 3 –, NO 2 –, and OH– are competitive for Sr 2+ binding, previous data from IX studies indicate that [SrOH] + is not the dominant species of concern in tank waste processing schemes.(Fiskum, Campbell and Trang-Le 2020) Our results show that the [Sr(NO3)]+ species and the [Sr(NO2)] + species form in considerable abundances, which may affect the ability to separate Sr using CST in nuclear waste separation processes.

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

Approaches for the Simulation of Coupled Processes in Evolving Fractured Porous Media Enabled by Exascale Computing

Models have historically represented fractured porous media with continuum descriptions that characterize the media using bulk parameters. The impact of small-scale features is not captured in these models, although they may be controlling the performance of subsurface applications. Pore-scale models can simulate processes in small-scale features by representing the pore space geometry explicitly but are computationally expensive for large domains. The alternative multiscale approach entails the combination of pore-scale and continuum-scale descriptions in a single framework. We use Chombo-Crunch, a computational capability that discretizes complex geometries with an adaptive, embedded boundary method to contrast these two approaches. Chombo-Crunch takes advantage of recent computational performance and memory bandwidth improvements resulting from the emergence of exascale computing resources. These combined improvements enable the efficient simulation of reactive transport in fractured media with a high degree of fidelity and the ability to capture the control small-scale processes exert on the overall medium evolution.

42 ENGINEERING↗
NASA NTRS1964

Applications of the method of Monte Carlo to problems in thermal radiation

A summary of the work involving the Monte Carlo method in the solution of problems in thermal radiation transfer is presented, which indicates general methods previously used for solving problems in which radiation is coupled with other modes of energy transfer. Previous work involving radiation in absorbing-emitting media is included. An example is outlined to indicate the use of the Monte Carlo method in the design of a space radiator. Suggestions are given for solution of a complex case incorporating the effects of coupled conduction, convection and radiation, wavelength dependent and selective surfaces, nonisothermal conditions, and strongly directional or nondiffuse emitting and reflecting surfaces. A discussion is given of the factors that may affect convergence, running time, and accuracy of the Monte Carlo solutions and of the advantages and disadvantages of this approach for practical problems. Also discussed are the case of programming for complex problems and the probable machine time requirements of the method.

THERMAL RADIATION↗
DOE OSTI2026

Finite-element boundary-integral simulation of thin wires and inhomogeneous penetrable bodies in subsurface multilayered anisotropic media

With the prevailing presence of drilling wells near the subsurface in mature oil and gas fields, the application of electromagnetic methods can be particularly challenging where the electromagnetic field is affected by the steel casing. In the past decades, borehole-to-surface and crosswell electromagnetic methods have been utilized for monitoring of reservoir and underground CO 2 storage. This paper presents a unified finite-element boundary-integral (FEBI) method capable of simultaneously modeling the complex electromagnetic interactions between thin metallic wires (representing steel casings) with 3D trajectory and arbitrary 3D inhomogeneous penetrable bodies (such as CO 2 plumes or hydrocarbon reservoirs) within anisotropic multilayered subsurface environments. Unlike existing approaches that treat these components separately or require dense discretization, or are limited to vertical wells, our unified formulation preserves flexible electromagnetic coupling while delivering improved computational efficiency. Assuming the background formation is multilayered anisotropic media, the surface integral equation method is applied to model the thin wires and boundaries of the inhomogeneous bodies. Meanwhile, the finite element method is applied to model the volume of inhomogeneous bodies. Here, the performance of the proposed FEBI method is assessed through comparison with reference numerical results and its practical significance is demonstrated through CO 2 plume monitoring scenarios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Spectroscopic Insight on Neodymium Solvation in Lithium Borohydride-Supported Electrolyte

Borohydride-based electrolytes have recently emerged as promising media for the electrodeposition of electropositive metals, including rare earth (RE) elements. While the presence of supporting alkali metal cations and RE counteranions provides essential electrochemical conductivity for achieving fast metal electrodeposition, interactions between the host ligand and solvated neodymium (Nd) complexes remain unclear. This study provides insights into the coordination structure of a concentrated and directly solvated Nd salt in a lithium borohydride-supported electrolyte. Our spectroscopic results indicate that the RE coordination environment is significantly influenced by the solvation mechanism, which can vary between metathesis and complexation pathways, primarily dictated by stoichiometric factors. Under dilute conditions, nearly complete metathesis of anions leads to a high coordination number for the host ligand (borohydride), consistent with the previously reported solvated Nd speciation in chlorine-free electrolytes. In contrast, concentrated dissolution of the Nd salt in the supported electrolyte is dominated by a complexation pathway featuring a Li-ion-paired complex with a low coordination number of the host ligand. Density functional theory (DFT) calculations indicated that the observed blue shift in the borohydride vibration was the result of an increase in electron density drawn into the terminal B–H interbond region from the hydride as the coordination changed from Li to Nd. In conjunction with DFT results, vibrational analyses allowed correlation of the experimental shifts associated with changes in Nd ligation and coordination spheres, further consolidating the prevalence of highly chloride-coordinated species under concentrated conditions. In conclusion, the outcomes of this work illuminate the distinctive and heterogeneous coordination structures that the electroactive RE species can adopt at high concentrations in lithium borohydride-supported electrolytes, as a key step to comprehend the reported metal electrodeposition performance in these media.

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

Elucidating the Transition of 3D Morphological Evolution of Binary Alloys in Molten Salts with Metal Ion Additives

Molten salts serve as effective high-temperature heat transfer fluids and thermal storage media used in a wide range of energy generation and storage facilities, including concentrated solar power plants, molten salt reactors and high-temperature batteries. However, at the salt–metal interfaces, a complex interplay of charge-transfer reactions involving various metal ions, generated either as fission products or through corrosion of structural materials, takes place. Simultaneously, there is a mass transport of ions or atoms within the molten salt and the parent alloys. The precise physical and chemical mechanisms leading to the diverse morphological changes in these materials remain unclear. Here, to address this knowledge gap, this work employed a combination of synchrotron X-ray nanotomography and electron microscopy to study the morphological and chemical evolution of Ni-20Cr in molten KCl-MgCl 2 , while considering the influence of metal ions (Ni 2+ , Ce 3+ , and Eu 3+ ) and variations in salt composition. Our research suggests that the interplay between interfacial diffusivity and reactivity determines the morphological evolution. The summary of the associated mass transport and reaction processes presented in this work is a step forward toward achieving a fundamental comprehension of the interactions between molten salts and alloys. Overall, the findings offer valuable insights for predicting the diverse chemical and structural alterations experienced by alloys in molten salt environments, thus aiding in the development of protective strategies for future applications involving molten salts.

36 - MATERIALS SCIENCE↗
NASA NTRS1994

An electromagnetic noncontacting sensor for thickness measurement in a dispersive medium

This paper describes a general purpose imaging technology developed by the U.S. Bureau of Mines (USBM) that, when fully implemented, will solve the general problem of 'seeing into the earth.' A first-generation radar coal thickness sensor, the RCTS-1, has been developed and field-tested in both underground and highwall mines. The noncontacting electromagnetic technique uses spatial modulation created by moving a simple sensor antenna in a direction along each axis to be measured while the complex reflection coefficient is measured at multiple frequencies over a two-to-one bandwidth. The antenna motion imparts spatial modulation to the data that enables signal processing to solve the problems of media, target, and antenna dispersion. Knowledge of the dielectric constant of the media is not necessary because the electrical properties of the media are determined automatically along with the distance to the target and thickness of each layer of the target. The sensor was developed as a navigation guidance sensor to accurately detect the coal/noncoal interface required for the USBM computer-assisted mining machine program. Other mining applications include the location of rock fractures, water-filled voids, and abandoned gas wells. These hazards can be detected in advance of the mining operation. This initiating technology is being expanded into a full three-dimensional (3-D) imaging system that will have applications in both the underground and surface environment.

Chufo, Robert L.↗
DOE OSTI2024

High-Temperature Gas-Cooled Pebble-Bed Reactors Running In And Transient Modeling Capabilities Demonstration

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor (HTGR) applications. The first part of this effort focuses on the running-in verification of Griffin's multiphysics capabilities, specifically for simulating the evolution of Pebble Bed reactor cores from startup to equilibrium. In the absence of validation data, code-to-code comparisons are conducted with Kugelpy, showing good agreement for key quantities like maximum power density and fresh core k-eigenvalue predictions. However, discrepancies in equilibrium core predictions suggest potential issues with cross sections, underscoring the need for further refinement and evaluation. The HTTF system analysis code benchmark involves RELAP5-3D, SAM, and GAMMA+ to assess their predictive capabilities for HTTF behavior under both normal operation and pressurized conduction cooldown (PCC) transient conditions. While there is good agreement in predicting major parameters such as coolant temperature, solid temperature, and flow distribution, discrepancies in transient behavior highlight differences in modeling approaches, nodalizations, and heat transfer models. The HTTF lower plenum CFD benchmark employs nekRS to simulate flow mixing phenomena, successfully capturing relevant flow physics and demonstrating mesh independence in complex geometries. Preliminary results suggest a relatively uniform temperature field but significant unsteadiness in the flow, requiring time-averaging analyses. The GPBR200 system analysis code benchmark uses SAM's core channel and porous media models, incorporating an RCCS loop for decay heat removal. During steady-state and transient conditions, including protected de-pressurized and pressurized loss of forced cooling (DLOFC and PLOFC), both models show good agreement in predicting temperature profiles and key parameters. Notably, while the core channel model underpredicts convective heat transfer effects, both models maintain temperatures well below the TRISO fuel safety limit. These benchmarking efforts collectively enhance the predictive capabilities of the tools used in HTGR design and safety analysis, guiding developments to improve their accuracy and applicability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRSDate not supplied

Compass Concurrent Engineering Lessons Learned in Remote and Hybrid Environments

The Compass Team at NASA’s Glenn Research Center (GRC) is a concurrent engineering team which specializes in conceptual spacecraft mission designs. Detailed descriptions of the team, its history, and its operating model can be found in [1] and [2]. During the COVID-19 pandemic, the team was required to move to remote (virtual) operation from their in-person model for approximately 22 months. As the restrictions began to lift and team members were able to return in-person to the Compass Lab, the team moved into a hybrid mode of operation, with some participants still tying in remotely some or all of the time. This paper discusses many of the lessons learned from these experiences, highlighting improvements, outstanding challenges and the tools and methodologies used to address both. Throughout this discussion the terms “in-person”, “hybrid”, “virtual” and “remote” will be used. For the purposes of this paper, “in-person” will be understood to mean when team members are interacting simultaneously, physically within the Compass Lab. “Remote” or “virtual” will refer to when interactions are happening between people who are not co-located using only technology to interface. “Hybrid” will refer to when two or more participants are physically located in the Compass Lab and one or more participant(s) is participating remotely. Media richness is described as “a medium’s ability to communicate effectively based of four factors. They are the capacity for immediate feedback, the number of cues and channels it utilizes, the degree of personalization it affords, and its ability to communicate using natural language” [3]. This theory will be referenced and discussed in multiple of the following sections due to its relevance when choosing how to operate in remote and hybrid modes, as well as in making tool selections. The key to selecting the best mode of communication lays in how complicated the discussion is and the level of ambiguity involved. Not all conversations or interactions require media rich mediums. For example, providing information about which there is little to no ambiguity can easily be done in less rich methods - such as email. A conversation including high levels of ambiguity and/or complex information is better suited to a richer medium, such as in-person or a video call with shared screens.

Concurrent Engineering↗
DOE OSTI2026

Nanorod Pair Complexes Manipulated via Magnetic Casimir Forces

Here, we demonstrate a route to modulate Casmir–Lifshitz interactions between anisotropic nanoparticles and magnetic fluids. By semiclassical quantum electrodynamics, we study ground state dispersion forces for cylindrical dielectric nanorods made of polystyrene (PS) and zinc oxide (ZnO) embedded in toluene-based host media with gold-coated magnetite nanoparticles and also predict magnetic contributions to the fully retarded excited state interaction. The variation in magnetic permeability enables tuning between repulsive and attractive interactions, and measurable magnetic Casimir traps are predicted between a pair of ZnO–PS nanoparticles whose equilibrium position can be modulated over an order of magnitude with a small variation in the size of the magnetite nanoparticle. This provides an alternative magnetic Casimir-effect pathway to reversibly tune quantum electromagnetic forces at the nanoscale for the assembly and enhancement of colloidal stability.

Excited-state interaction↗
DOE OSTI2025

2025 ASMS Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS

Title (20 words): Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS Introduction (120 words): One of the challenges facing wide-spread adoption of nuclear power is the development of efficient separation processes for used nuclear fuel. The molecules in separation processes are subjected to an extreme environment due to the high radiation fields from the used fuel and highly acidic media used for fuel dissolution, which results in significant molecular degradation, leading to reduced process efficiency. These degradation products must be identified and studied so mitigation strategies can be developed to maintain process efficiency. However, complex systems can have many degradation products, complicating identification. Untargeted analysis tools could be used to understand radiation chemistry in complex systems. However, this would necessitate improved understanding of the gas-phase fragmentation mechanisms of fuel cycle molecules like tetraoctyldiglycolamide (TODGA). Methods (120 words): The gas-phase fragmentation of protonated TODGA was investigated using collision-induced dissociation (CID), resonance ejection, and infrared ion spectroscopy (IRIS). CID and resonance ejection experiments were conducted using a Bruker Daltonics (Bremen, Gemany) SolariX XR fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. IRIS spectra of protonated TODGA and its two CID fragmentation products were measured using a modified Bruker amaZon Speed ETD 3D quadrupole ion trap mass spectrometer coupled to the Free Electron Lasers for Infrared eXperiments (FELIX) free electron laser. Measured spectra were compared with density functional theory (DFT) calculations using the Gaussian 16, Revision C.02 software package with the ?B97X-D functional and def2-TZVPP basis sets. Candidate structures were generated using the CREST 3.0 conformational sampling software tool. Preliminary Data (300 words): Collision-induced dissociation of protonated TODGA ([C36H73N2O3]+, m/z=581.562) results two fragment ions, one at m/z=340.285 assigned as [C20H38NO3]+ and the other at m/z=312.290, assigned as [C19H38NO2]+. Based on the assigned formula and the structure of protonated TODGA, the fragment at m/z=340.285 is likely formed from elimination of neutral dioctylamine. Comparison of the IRIS spectrum of m/z=340.285 with DFT predictions suggests it contains a ring structure, and is assigned as N-octyl-N-(6-oxo-1,4-dioxan-2-ylidene)octan-1-aminium. Based on this structure and the structure of protonated TODGA, we hypothesize this fragment formed from elimination of neutral dioctylamine followed by a ring closure mechanism. Comparison of the IRIS spectrum of the fragment at m/z=312.290 with DFT predictions also indicated the presence of a ring structure, assigned as N-(1,3-dioxolan-4-ylidene)-N-octyloctan-1-aminium. This product could be formed from elimination of carbon monoxide from the ring of m/z=340.285 as a sequential fragmentation or formed directly from protonated TODGA via elimination of neutral N,N-dioctylformamide followed by a ring closure. Resonance ejection experiments where m/z=340 was continuously ejected from the IRC cell showed no decrease in intensity of m/z=312.290 across several collision energies, suggesting that the later, direct formation mechanism, dominates. The location of the ionizing proton in protonated TODGA is important for modeling the fragmentation mechanisms. DFT calculations suggested that the position of bands involving the coupled vibrations of the amide C—N and C=O bonds in TODGA are the most sensitive to proton location. Evaluation of the IRIS spectrum of protonated TODGA suggests that the ionizing proton is located between the two amid oxygens. This protonation location was calculated to lie approximately 30 kJ/mol lower in energy than the next lowest energy location, with the proton located solely on one of the amide oxygens. Novel aspect (20 words): Infrared ion spectroscopy combined with resonance ejection experiments and density functional theory to probe the collision-induced dissociation mechanism of tetraoctyldiglycolamide.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2025

The role of fast and slow dynamics in nonlinear resonant ultrasound spectroscopy of consolidated granular materials

Abstract Elastic nonlinearity observed in consolidated granular media can be attributed to the combination of slow and fast effects, which give rise to hysteresis and relaxation of both modulus and damping after the sample is perturbed. A consequence is a high level of complexity in the measurements of the sample linear and nonlinear elastic parameters. The results of experiments are dependent on the experimental protocol that is adopted to measure the relevant quantities and it is hard to quantify parameters with accuracy and repeatability. Here we focus on examining Nonlinear Resonant Ultrasound Spectroscopy, showing experimentally the role of slow dynamics in the process and quantifying/discussing its influence on the quantification of nonlinearity. We also propose a model to describe the process, which shows that different contributions to nonlinearity (e.g., classical and hysteretic) could be due to physical features (defects) relaxing with different relaxation times.

Science & Technology - Other Topics↗
NASA NTRS1990

Multiple sensor multifrequency eddy current monitor for solidification and growth

A compact cylindrical multisensor eddy current measuring system with integral furnace was develop to monitor II-VI crystal growth to provide interfacial information, solutal segregation, and conductivities of the growth materials. The use of an array of sensors surrounding the furnace element allows one to monitor the volume of interest. Coupling these data with inverse multifrequency analysis allows radial conductivity profiles to be generated at each sensor position. These outputs were incorporated to control the processes within the melt volume. The standard eddy current system functions with materials whose electric conductivities are as low as 2E2 Mhos/m. A need was seen to extend the measurement range to poorly conducting media so the unit was modified to allow measurement of materials conductivities 4 order of magnitude lower and bulk dielectric properties. Typically these included submicron thick films and semiinsulating GaAs. This system was used to monitor complex heat transfer in grey bodies as well as semiconductor and metallic solidification.

Wallace, John↗
DOE OSTI2024

Sorption kinetics and stability of conventional adsorbents for mercury remediation

In-situ remediation of mercury at numerous contaminated sites worldwide is a challenging and costly endeavor due to the persistency of this contaminant. In this study, we evaluated eight commercially available sorbent media ranging from carbon-, clays- and silica-based materials (PBC– Biochar, eSorb – Sorbster, nsPAC – Powdered Activated Carbon, fsPAC – Powdered Activated Carbon with Mackinawite, F300 – Filtrasorb 300, Si-SH – Silica Thiol, eBind – RemBind, Q-Clay – Organoclay PM-199), for their effectiveness in sorbing mercury (Hg 2+ ) and mercury complexed with dissolved organic matter (Hg-DOM). Under the chosen experimental conditions of this study, results showed that in the absence of DOM, the kinetic rates of Hg 2+ sorption onto the evaluated sorbents were in the order of 0.31 min-1 (Si-SH) to 2.98 min -1 (nsPAC), whereas in the presence of DOM, the rates varied from 0.16 min-1 (F300) to 0.95 min -1 (nsPAC). The measured sorption capacity for Hg 2+ in the absence of DOM varied from 3.02 mg/g (Q-Clay) to 35.15 mg/g (Si-SH), whereas in the presence of DOM, calculated partition coefficient (KD) ranged from 69.7 mL/g (Q-Clay) to 41,510 mL/g (Si-SH). Furthermore, kinetic data suggest liquid film diffusion was the rate-limiting steps governing mercury sorption onto the studied media. Overall, the obtained study parameters (kinetics/isotherm) are particularly important in informing robust engineering designs for deployment of the vast majority of evaluated sorbents. Thus, sorbent-based strategies offer viable solutions for cost-effective cleanup of mercury at industrially contaminated sites.

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

Water-induced surface ordering facilitating the microcutting of ductile metals

Metal cutting is a crucial process in modern manufacturing. Enhancing the machinability of metals can significantly improve their production efficiency and surface integrity. Coating surface-active media (SAM) on the free surface of the metals before cutting is an easy method to improve machinability, which usually pertains to the category of the renowned Rehbinder effect. However, the existing SAM is usually hazardous and complex materials. Besides, the effect of SAM on the local structure of the metal surface remains unclear. In this study, water is employed as a simple yet often overlooked SAM in the microcutting of copper. Using water as SAM also allows the employment of X-ray absorption fine structure spectroscopy (XAFS) to study the local structure of copper with and without water coating. Results show that water coating on the free surface of copper can significantly reduce the cutting force and chip thickness, and improve the surface finish. Interestingly, removing the water coating enables the recovery of the cutting force, demonstrating a reversible effect. Based on the XAFS results and molecular dynamics simulation, a water-induced surface ordering mechanism is proposed to explain the findings from the microcutting experiments. This mechanism suggests that water molecules can induce surface ordering in copper, resulting in reduced surface energy and fracture toughness of copper, thus enhancing machinability. In conclusion, this work provides valuable insights into the comprehension of the Rehbinder effect and shows that picometer-scale modifications of the surface atom arrangement can considerably alter the deformation mode of metals, paving the way for the development of new manufacturing processes.

36 MATERIALS SCIENCE↗
NASA NTRS1983

Wave scattering and guidance by dielectric waveguides with periodic surfaces

The extended-boundary-condition method is used to solve the wave scattering and guidance by a stratified medium with a periodic surface. First, the transition matrix for two sets of Flocquet waves incident from above and below a periodic surface separating the two media is derived. The results are applied to the diffraction of an incident plane wave by a multilayered medium having a periodic surface. The theoretical results can easily be applied to the guidance problem for which the complex guiding constant of the periodic structure is obtained. Numerical results for both the scattering and the guidance problems are obtained. The theory is compared with experiment for metallic gratings with surface-plasmon excitations. The calculated guiding constants of a thin film with a periodic sinusoidal surface are also compared with other numerical methods.

Chuang, S. L.↗
NASA NTRS1986

Presence, segregation and reactivity of H, C and N dissolved in some refractory oxides

The sources of impurities, particularly carbon, in high melting oxides and silicates are discussed, along with detection and quantification methods. The impurities are important for their effects on bulk material properties through the media of, e.g., surface or grain boundary characteristics. The impurities are usually encountered by the contact of the oxide (refractory) material with volatiles such as H2O and CO2, which become incorporated in the material and form anion complexes with oxygen acting as a covalent bonded ligand. The specific processes undergone by MgO in assimilating C impurities are delineated, using data obtained with X-ray photoelectron spectroscopy, Auger electron spectroscopy, secondary ion mass spectrometry and nuclear reaction profiling. Finally, maintenance of a supersaturated solid solution with C impurities by space charge control is described as a means of offset impurity effects.

Freund, F.↗
NASA NTRS1988

Calculation of the recirculating compressible flow downstream a sudden axisymmetric expansion

Significant progress has been made during the last five years to adapt conventional Navier-Stokes solver for handling nonconservative equations. A primary type of application is to use transport equation turbulence models, but the extension is also possible for describing the transport of nonpassive scalars, such as in reactive media. Among others, combustion and gas dissociation phenomena are topics needing a considerable research effort. An implicit two step scheme based on the well-known MacCormack scheme has been modified to treat compressible turbulent flows on complex geometries. Implicit treatment of nonconservative equations (in the present case a two-equation turbulence model) opens the way to the coupled solution of thermochemical transport equations.

Vandromme, D.↗