Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactive Flow”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Tracking Diverse Minerals, Hungry Organisms, and Dangerous Contaminants Using Reactive Transport Models

Beneath our feet is a fascinating world of flowing water, cosmopolitan microbes, and complex mineral assemblages. Yet we see none of it from above. Our quest to investigate these complex subsurface interactions has led to the development of reactive transport models. These are computer algorithms that allow us to explore, in a virtual way, the natural dynamics of Earth's systems and our anthropogenic impact on those systems. Here, we explain the concepts behind reactive transport models—which include the transport of aqueous species and the descriptions of biogeochemical reactions involving solutes, surfaces and microorganisms—and introduce to reactive transport applications in terrestrial and marine environments.

58 GEOSCIENCES↗
DOE OSTI2024

Fluid inertia controls mineral precipitation and clogging in pore to network-scale flows

Mineral precipitation caused by fluid mixing presents complex control and predictability challenges in a variety of natural and engineering processes, including carbon mineralization, geothermal energy, and microfluidics. Precipitation dynamics, particularly under the influence of fluid flow, remain poorly understood. Combining microfluidic experiments and three-dimensional reactive transport simulations, we demonstrate that fluid inertia controls mineral precipitation and clogging at flow intersections, even in laminar flows. We observe distinct precipitation regimes as a function of Reynolds number (Re). At low Reynolds numbers (Re < 10), precipitates form a thin, dense layer along the mixing interface, which shuts precipitation off, while at high Reynolds numbers (Re > 50), strong three-dimensional flows significantly enhance precipitation over the entire intersection, resulting in rapid clogging. When injection rates from two inlets are uneven, flow symmetry-breaking leads to unexpected flow bifurcation phenomena, which result in enhanced concurrent precipitation in both downstream channels. Finally, we extend our findings to rough channel networks and demonstrate that the identified inertial effects on precipitation at the intersection scale are also present and even more dramatic at the network scale. This study sheds light on the fundamental mechanisms underlying mixing-induced mineral precipitation and provides a framework for designing and optimizing processes involving mineral precipitation.

Science & Technology - Other Topics↗
DOE OSTI2018

Field Demonstration of the Combined Use of Thermally-Degrading and Cation-Exchanging Tracers to Predict Thermal Drawdown in a Geothermal Reservoir

This paper presents the results and interpretation of a reactive tracer test in a geothermal reservoir involving a cation-exchanging tracer to interrogate flow pathway surface area to volume ratios and a thermally-degrading tracer to interrogate average flow pathway temperatures. The test interpretations, which are based on comparing the reactive tracer responses to that of a nonreactive tracer, are used to predict when thermal drawdown will occur at the production well due to injecting cooler, spent geothermal water into an injection well. It is shown that the cation- exchanging tracer provides constraints on the timing of thermal drawdown, and the thermally-degrading tracer provides constraints on the magnitude of thermal drawdown. Thus, the combined use of the two reactive tracers provides a better prediction of thermal drawdown than either one alone. Uncertainties associated with the method, including assumptions made about flow pathway geometry, are discussed. It is also shown that interrogating average temperatures along flow pathways using a thermally-degrading tracer at different times can be translated into an estimate of when thermal drawdown will occur without direct knowledge of flow pathway geometry.

15 GEOTHERMAL ENERGY↗
DOE OSTI2018

Molecular admixtures and impurities in atmospheric pressure plasma jets

A more complete understanding of reactive chemistry generated by atmospheric pressure plasma jets (APPJs) is critical to many emerging medical, agricultural, and water treatment applications. Adding molecular gases to the noble working gas which flows through the jet is a common method to tailor the resulting production of reactive oxygen and nitrogen species (RONS). In this paper, results are discussed from a computational investigation of the consequences of H 2 O and O 2 admixtures on the reactive chemistry of He APPJs flowing into humid air. This investigation, performed with a 2-dimensional plasma hydrodynamics model, addresses the RONS that are initially produced and the evolution of that chemistry on longer time scales. Without an admixture, the impurities in 99.999% pure helium are a major source of RONS. The addition of H 2 O decreases the production of reactive nitrogen species (RNS) and increases the production of reactive oxygen species (ROS). The addition of O 2 significantly decreases the production of RNS, as well as hydrogen-containing ROS, but increases the production of ROS without hydrogen. This selectivity comes from the lower ionization energy of O 2 compared to N 2 and H 2 O, which then allows for charge exchange reactions. These charge exchange reactions change the RONS which are produced in the afterglow by dissociative recombination. The consequences of impurities were also examined. Humid air impurities as low as 10 ppm in the helium can account for 79%-98% of the production of most RONS in the absence of an intentional admixture. The degree to which the impurities affect the RONS production depends on the electrode configuration and can be reduced by molecular admixtures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Application of NEAMS Codes to Capture MSR Phenomena

This report documents work completed in FY21 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program’s Molten Salt Reactor (MSR) Application Drivers activity at Argonne. The common focus was on identifying the modeling and simulation functional requirements for designing, licensing, and operating MSRs and applying those capabilities already developed in NEAMS tools to example problems of interest. The four main parts of this report each focuses on a specific area of simulation physics as it relates to MSR phenomena: fuel evolution, chemistry, computational fluid dynamics (CFD), and systems analysis. In terms of fuel evolution, which includes depletion, decay, on-line separations, and transmutation, the current state of computational capabilities for modeling this behavior in liquid-fueled molten salt reactor is discussed. Some of the functional requirements to accomplish the various applications of MSR fuel depletion modeling are highlighted, followed by a summary of recent approaches and code development activities. The chemistry functional requirements were discussed in the context of several applications of high importance for MSRs, such as corrosion, salt chemistry, and species transport. Each of these types of chemistry modeling have considerable impact on various aspects of reactor applications, including informing on reactor designs, improving operational efficiencies, analyzing safety and reactivity concerns, and estimating the mechanistic source term of the reactor. A brief overview is also provided on code development activities ongoing under NEAMS relevant to chemistry modeling of MSRs. In terms of CFD applications, the state-of-the-art spectral element code, Nek5000 was used to model the fluid dynamics within a full core of the Molten Salt Fast Reactor (MSFR) concept designed as part of the Euratom EVOL project. This concept was selected as the challenge problem because of its similar features to several U.S. industry concepts. The goal was to model some of the fast MSR design challenges, including potential large internal re-circulations, the need of accurate tracking of delayed neutron precursors (DNP), and the lack of relevant thermal-hydraulics models/correlations, etc. Therefore, a series of CFD models were created for the MSFR core cavity using a k – τ model two-equation model for the turbulence. These first full core results demonstrated that any potential recirculation zones could be properly identified with the current NEAMS CFD capabilities. The development of these models will also set the stage for future testing of Nek5000’s functionalities to model other MSR thermal-fluid phenomena. Lastly, the validation of SAM against experimental data from the Molten Salt Reactor Experiment, which started in FY19, continues with the inclusion of modeling the reactivity insertion tests. This involved using SAM and its point kinetics model for flowing fuel salt to recreate the time dependent power changes and response after positive reactivity insertions at the 1, 5, and 8 MWt power levels. Through these exercises, several code modifications were suggested to the SAM development team and accommodated to enable closer agreement with solid technical and physical justifications. These include adding a moderator reactivity feedback coefficient as an available input and modifying the solution approach for the point kinetics model. Additional SAM development suggestions for flowing fuel MSRs include adding the capability to allow the moderator power change proportionally with the reactor power and enabling specification of the power and DNP distributions separately.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Implementation of a Model Predictive Control Strategy to Regulate Temperature Inside Plug-Flow Solar Reactor With Countercurrent Flow

Abstract Solar-driven thermochemical energy storage systems are proven to be promising energy carriers (solar fuels) to utilize solar energy by using reactive solid-state pellets. However, the production of solar fuel requires a quasi-steady-state process temperature, which represents the main challenge due to the transient nature of solar power. In this work, an adaptive model predictive controller (MPC) is presented to regulate the temperature inside a tubular solar reactor to produce solid-state solar fuel for long-term thermal storage systems. The solar reactor system consists of a vertical tube heated circumferentially over a segment of its length by concentrated solar power, and the reactive pellets (MgMn2O4) are fed from the top end and flow downwards through the heated tube. A countercurrent flowing gas supplied from the lower end interacts with flowing pellets to reduce it thermochemically at a temperature range of 1000—1500 °C. A low-order physical model was developed to simulate the dynamics of the solar reactor including the reaction kinetics, and the proposed model was validated numerically by using a 7-kW electric furnace. The numerical model then was utilized to design the MPC controller, where the control system consists of an MPC code linked to an adaptive system identification code that updates system parameters online to ensure system robustness against external disturbances (sudden change in the flow inside the reactor), model mismatches, and uncertainty. The MPC controller parameters are tuned to enhance the system performance with minimum steady-state error and overshoot. The controller is tested to track different temperature ranges between 500 °C and 1400 °C with different particles/gas mass flowrates and ramping temperature profiles. Results show that the MPC controller successfully regulated the reactor temperature within ± 1 °C of its setpoint and maintained robust performance with minimum input effort when subjected to sudden changes in the amount of flowing media and the presence of chemical reaction.

Engineering↗
NASA NTRS1971

Selective Catalytic Combustion Sensors for Reactive Organic Analysis

Sensors involving a vanadia-alumina catalyst bed-thermocouple assembly satisfy requirements for simple, reproducible and rapid continuous analysis or reactive organics. Responses generally increase with temperature to 400 C and increase to a maximum with flow rate/catalyst volume. Selectivity decreases with temperature. Response time decreases with flow rate and increases with catalyst volume. At chosen optimum conditions calculated response which is additive and linear agrees better with photochemical reactivity than other methods for various automotive sources, and response to vehicle exhaust is insensitive to flow rate. Application to measurement of total reactive organics in vehicle exhaust as well as for gas chromatography detection illustrate utility. The approach appears generally applicable to high thermal effect reactions involving first order kinetics.

Innes, W. B.↗
DOE OSTI2023

Experimental investigation of power transient flow boiling

To advance the mechanistic understanding of power transients in flow boiling, experimental investigations are performed to capture the thermal response of four different kinds of cladding materials, including FeCrAls, to Fuchs reactivity initiated accident transients and linear ramp power transients. The transition boiling regime is phenomenologically reflected in power-transient flow boiling, which is different from steady-state boiling. In light of this, the critical heat flux that is usually featured with temperature overshooting does not indicate thermal safety margin of cladding materials as conservatively as the maximum heat flux, which is greater than critical heat flux in the power transient boiling curve. This is physically attributed to the thermal energy deposition on the cladding wall. In addition, the thermal energy released from the cladding wall results in different boiling heat transfer coefficient for the decreasing and increasing power stages respectively during the Fuchs power transient. The maximum heat flux difference gap, which appears due to changes in cladding materials and transient time scales, is appreciable under an intermediate heat convection regime. However, this difference gap is gradually weakened by the progressive increasing of mass flux and/or inlet subcooling due to the enhanced dominance of heat convection over heat conduction. Finally, the small difference gap of maximum heat flux is found to be insignificant under the weak heat convection because of power transient induced annular flow instability.

42 ENGINEERING↗
DOE OSTI2025

Electrified vapour deposition at ultrahigh temperature and atmospheric pressure for nanomaterials synthesis

Vapour-phase synthesis methods have shown promise for the scalable synthesis of nanomaterials and coatings. However, the vaporization of different precursors for the synthesis of a broad nanomaterial space, particularly at atmospheric pressure, while maintaining compositional and structural control of the final product is challenging. Here we report the generation of an ultrahigh-temperature atomic vapour at atmospheric pressure based on electrified heating, for the growth of multi-elemental nanomaterials and thin films. This process relies on a reactor design whereby solid-state precursors are vaporized within a semi-confined space beneath an electrified heater that can reach ~3,000 K. The proximity of the heater rapidly breaks down the bonds of metal salt precursors and decomposes them into an atomic vapour that expands into a high-temperature (>2,000 K), highly reactive and high-flux vapour (10 21 –10 22 atoms per cm 2 per second) that travels upwards in a directional flow. When mixed with entrained ambient gases, the highly reactive atomic species rapidly nucleate and grow into the desired final products, including alloys, oxides, sulfides and thin films, which can be deposited on a low-temperature substrate. This EVD approach can synthesize a broad range of functional nanomaterials at atmospheric pressure, including single-phase multi-elemental nanomaterials formed under thermodynamically non-equilibrium conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Growth of (1 1 1) and (2 0 0) orientation cubic MgZnO thin films under different oxygen flow rate by PLD method and its difference in element composition and optical absorption characteristics

Under different migration energy of reactive Mg, Zn and O atoms from MgZnO target at different oxygen flow rate, (2 0 0) and (1 1 1) orientations MgZnO thin films with cubic structure were fabricated on fused quartz substrate by PLD method. And MgZnO thin film possesses relatively higher Zn composition and lower Mg composition when deposited more along (1 1 1) orientation. The band gap and UV absorption characteristics of MgZnO thin film do not change completely in accordance with the Mg/Zn atom ratio of MgZnO thin films deposited at different oxygen flow rate, but influenced more by the ratio between Mg and Zn atoms that combined with O atoms in MgZnO crystal lattice and the grain boundary density of MgZnO thin films deposited at different oxygen flow rate.

36 MATERIALS SCIENCE↗
DOE OSTI2017

On the physics-based processes behind production-induced seismicity in natural gas fields

Induced seismicity due to natural gas production is observed at different sites worldwide. Common understanding states that the pressure drop caused by gas production leads to compaction, which affects the stress field in the reservoir and the surrounding rock formations and hence reactivates preexisting faults and induces earthquakes. We show that the multiphase fluid flow involved in natural gas extraction activities should be included. We use a fully coupled fluid flow and geomechanics simulator, which accounts for stress-dependent permeability and linear poroelasticity, to better determine the conditions leading to fault reactivation. In our model setup, gas is produced from a porous reservoir, divided into two compartments that are offset by a normal fault. Results show that fluid flow plays a major role in pore pressure and stress evolution within the fault. Fault strength is significantly reduced due to fluid flow into the fault zone from the neighboring reservoir compartment and other formations. We also analyze scenarios for minimizing seismicity after a period of production, such as (i) well shut-in and (ii) gas reinjection. In the case of well shut-in, a highly stressed fault zone can still be reactivated several decades after production has ceased, although on average the shut-in results in a reduction in seismicity. In the case of gas reinjection, fault reactivation can be avoided if gas is injected directly into the compartment under depletion. However, gas reinjection into a neighboring compartment does not stop the fault from being reactivated.

03 NATURAL GAS↗
DOE OSTI2015

The response of soil carbon storage and microbially mediated carbon turnover to simulated climatic disturbance in a northern peatland forest. Revisiting the concept of soil organic matter recalcitrance

The goal of this project was to investigate changes in the structure of dissolved and solid phase organic matter, the production of CO 2 and CH 4 , and the composition of decomposer microbial communities in response to the climatic forcing of environmental processes that determine the balance between carbon gas production versus storage and sequestration in peatlands. Cutting-edge analytical chemistry and next generation sequencing of microbial genes were been applied to habitats at the Marcell Experimental Forest (MEF), where the US DOE’s Oak Ridge National Laboratory and the USDA Forest Service are constructing a large-scale ecosystem study entitled, “Spruce and Peatland Responses Under Climatic and Environmental Change”(SPRUCE). Our study represented a comprehensive characterization of the sources, transformation, and decomposition of organic matter in the S1 bog at MEF. Multiple lines of evidence point to distinct, vertical zones of organic matter transformation: 1) the acrotelm consisting of living mosses, root material, and newly formed litter (0-30 cm), 2) the mesotelm, a mid-depth transition zone (30-75 cm) characterized by labile organic C compounds and intense decomposition, and 3) the underlying catotelm (below 75cm) characterized by refractory organic compounds as well as relatively low decomposition rates. These zones are in part defined by physical changes in hydraulic conductivity and water table depth. O-alkyl-C, which represents the carbohydrate fraction in the peat, was shown to be an excellent proxy for soil decomposition rates. The carbon cycle in deep peat was shown to be fueled by modern carbon sources further indicating that hydrology and surface vegetation play a role in belowground carbon cycling. We provide the first metagenomic study of an ombrotrophic peat bog, with novel insights into microbial specialization and functions in this unique terrestrial ecosystem. Vertical structuring of microbial communities closely paralleled the chemical evolution of peat, with large shifts in microbial populations occurring in the biogeochemical hotspot, the mesotelm, where the highest rates of decomposition were detected. Stable isotope geochemistry and potential rates of methane production paralleled vertical changes in methanogen community composition to indicate a predominance of acetoclastic methanogenesis mediated by the Methanosarcinales in the mesotelm, while hydrogen-utilizing methanogens dominated in the deeper catotelm. Evidence pointed to the availability of phosphorus as well as nitrogen limiting the microbially-mediated turnover of organic carbon at MEF. Prior to initiation of the experimental treatments, our study provided key baseline data for the SPRUCE site on the vertical stratification of peat decomposition, key enzymatic pathways, and microbial taxa containing these pathways. The sensitivity of soil carbon turnover to climate change is strongly linked to recalcitrant carbon stocks and the temperature sensitivity of decomposition is thought to increase with increasing molecular complexity of carbon substrates. This project delivered results on how climate change perturbations impact the microbially-mediated turnover of recalcitrant organic matter in peatland forest soils, both under controlled conditions in the laboratory and at the ecosystem-scale in the field. This project revisited the concept of “recalcitrance” in the regulation of soil carbon turnover using a combination of natural abundance radiocarbon and optical spectroscopic measurements on bulk DOM, and high resolution molecular characterization of DOM. The project elucidated how organic matter reactivity and decomposition will respond to climate change in a both a qualitative (organic matter lability) and quantitiative (increased rates) manner. An Aromaticity Index was developed to represent a more direct and accurate parameter for modeling of DOM reactivity in peatlands. The abundance and community composition of soil microorganisms that mediate C cycling were interrogated with depth in the peat, with season, and in manipulated climate enclosures at unprecedented resolution. Therefore this project delivered strategic new insights on the functioning of peatland ecosystems that collectively store approximately one-third of the world's soil carbon. Furthermore, results from the detailed characterization of DOM lability and microbial community structure/ function will be employed to further develop biogeochemical models to include microbial respiration pathways as well as to track carbon flow with a term that incorporates relative reactivity based on aromaticity index. As it stands now, detailed soil organic matter structure and microbial parameters are not included in Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Reactive Transport Modeling of Anthropogenic Carbon Mineralization in Stacked Columbia River Basalt Reservoirs

Numerical simulation of CO2 storage in basalts and related reactive lithologies requires modeling complex, coupled hydrologic and chemical processes, including multi-phase flow and transport, partitioning of CO2 into the aqueous phase, and chemical interactions with aqueous fluids and rock minerals. We conducted reactive transport simulations of the Wallula pilot-scale CO2 injection into the flow tops of the Grande Ronde Basalt using our PNNL STOMP-CO2 simulator with the ECKEChem reactive module. Our mineralization simulation of the ~1,000 tons of injected CO2 into the interflow zones was based on the hydrologic transport model we previously developed. For this work, the simulations considered geochemical reactions involving the basalt components, precipitates, formation brine, and injected CO2. In our benchmark case, carbonate minerals precipitated, resulting in ~20% of the CO2 being mineralized in 10 years. Increasing the reaction rate of a single primary mineral phase (clinopyroxene) by an order of magnitude resulted in a carbon mineralization reaction extent of ~90% over the same time interval. Based on these initial sensitivity analysis results, it is clear that a thorough understanding of primary mineral dissolution rates is required for accurately predicting long-term fate and transport of injected CO2 into basalt formations. Our reactive transport numerical simulations will be key components of commercial-scale CO2 storage operation permitting, de-risking, and optimization in mafic and ultramafic reservoirs.

Cao, Ruoshi↗
DOE OSTI2019

Mechanisms of Methane Hydrate Formation in Geological Systems

Natural gas hydrate is ice-like mixture of gas (mostly methane) and water that is widely found in sediments along the world's continental margins and within and beneath permafrost and glaciers in a near-surface depth interval where the pressure is sufficiently high and temperature sufficiently low for gas hydrate to be stable. We categorize the myriad of geological gas hydrate deposits into five characteristic types. We then review the multiple quantitative models that have proposed to describe the genesis of these deposits and describe how each may have formed. We emphasize the importance of coupling multiphase flow (free gas and liquid water) and multicomponent reactive transport with geological history to describe the dynamical processes of gas hydrate formation and evolution in geological systems. A better insight into the kinetics of methane formation from microbial biogenesis and the processes of multiphase flow at the pore scale will advance our knowledge of how these systems form. By understanding the generation and evolution of gas hydrate through time, we will better decipher the role of gas hydrate in the carbon cycle, its potential to contribute to climate change and geohazards, and how to design optimal strategies for gas production from hydrate reservoirs.

58 GEOSCIENCES↗
NASA NTRS2002

Effects of aging on vasoconstrictor and mechanical properties of rat skeletal muscle arterioles

Exercise capacity and skeletal muscle blood flow during exercise are reduced with advancing age. This reduction in blood flow capacity may be related to increased reactivity of skeletal muscle resistance vessels to vasoconstrictor stimuli. The purpose of this study was to test the hypothesis that aging results in increased vasoconstrictor responses of skeletal muscle resistance arterioles. First-order (1A) arterioles (90-220 microm) from the gastrocnemius and soleus muscles of young (4 mo) and aged (24 mo) Fischer-344 rats were isolated, cannulated, and pressurized via hydrostatic reservoirs. Vasoconstriction in response to increases in norepinephrine (NE; 1 x 10(-9)-1 x 10(-4) M) and KCl (20-100 mM) concentrations and increases in intraluminal pressure (10-130 cmH(2)O) were evaluated in the absence of flow. Responses to NE and KCl were similar in both soleus and gastrocnemius muscle arterioles from young and aged rats. In contrast, active myogenic responses to changes in intraluminal pressure were diminished in soleus and gastrocnemius arterioles from aged rats. To assess whether alterations in the mechanical properties of resistance arterioles underlie altered myogenic responsiveness, passive diameter responses to pressure and mechanical stiffness were evaluated. There was no effect of age on the structural behavior (passive pressure-diameter relationship) or stiffness of arterioles from either the soleus or gastrocnemius muscles. These results suggest that aging does not result in a nonspecific decrease in vasoconstrictor responsiveness of skeletal muscle arterioles. Rather, aging-induced adaptations of vasoreactivity of resistance arterioles appear to be limited to mechanisms that are uniquely involved in the signaling of the myogenic response.

NASA Discipline Musculoskeletal↗
DOE OSTI2018

The influence of Preslip Sealing on the Permeability Evolution of Fractures and Faults

The evolution of permeability on fractures and faults during the full earthquake cycle is shown to be sensitive to sealing during the repose phase. We explore the combined effect of static loading followed by fracture reactivation on permeability evolution via slide-hold-slide experiments. During the hold periods, permeability exhibits a slow but continuous reduction. The permeability decay is consistent with power law compaction of the aperture coupled with cubic law flow. With increasing hold periods, permeability evolves following reactivation from net reduction to net increase with the magnitude of the permeability change dependent on the hold period. This implies that the tight interlocking of asperities during interseismic repose primes the fault for permeability enhancement following reactivation. The inferred mechanism is via shear dilation with the probable involvement of unclogging. This result identifies that preslip sealing during repose is an essential component in the cyclic permeability evolution throughout the seismic cycle.

58 GEOSCIENCES↗
DOE OSTI2020

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗