Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Advanced Modeling and Simulation”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Energy Efficient Mobility Systems 2018 Annual Progress Report

Energy Efficient Mobility Systems (EEMS) Program activities during FY 2018 focused on analytical research to understand the impacts that new mobility technologies and services will have at the vehicle, traveler, and overall transportation system-level. This research included the development of vehicle and transportation system simulation models and tools to evaluate the complex interactions among the various actors within the mobility landscape, analysis of empirical data to characterize which solutions may provide the largest benefits, and development of new control systems and algorithms that use vehicle connectivity and automation to improve the performance and efficiency of individual vehicles as well as the overall traffic system This document presents a brief overview of the EEMS Program and documents progress and results for projects within four of the five EEMS activity areas: (1) the SMART (Systems and Modeling for Accelerated Research in Transportation) Mobility Lab Consortium, (2) High Performance Computing and Big Data Solutions for Mobility Data, (3) Advanced R&D Projects conducted by industry and academia, and (4) Core Modeling, Simulation, and Evaluation, Similarly, the remaining EEMS activity area – (5) Living Labs (managed under VTO’s Technology Integration Program). Each of the individual progress reports provide a project overview and highlights of the technical results.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1994

Ice-atmosphere interactions in the central Arctic: Remotely-sensed and simulated ice concentration and motion

The response of the Beaufort Sea ice pack to the passage of strong low pressure systems is studied using SAR (synthetic aperture radar), AVHRR (advanced very high resolution radiometer), and SSM/I (special sensor microwave/imager) data combined with model simulations. The SAR, AVHRR, and modeled concentrations concur generally in showing a 1 to 5 percent decrease in ice fraction during the passage of the lows. SSM/I derived concentrations appear to underestimate overall concentration and overestimate the opening of the pack relative to the other data sets. Ice motion derived from SAR and AVHRR and model simulations using two ice rheologies agree reasonably well in direction and magnitude, although the cavitating fluid approximation tends to overestimate ice motion compared to a viscous plastic rheology. Ice motions from SAR, AVHRR, and buoy observations can be merged to yield uniform grids particularly well suited to synoptic scale studies and model validation.

Maslanik, J. A.↗
NASA NTRS1994

Particle kinetic simulation of high altitude hypervelocity flight

Rarefied flows about hypersonic vehicles entering the upper atmosphere or through nozzles expanding into a near vacuum may only be simulated accurately with a direct simulation Monte Carlo (DSMC) method. Under this grant, researchers enhanced the models employed in the DSMC method and performed simulations in support of existing NASA projects or missions. DSMC models were developed and validated for simulating rotational, vibrational, and chemical relaxation in high-temperature flows, including effects of quantized anharmonic oscillators and temperature-dependent relaxation rates. State-of-the-art advancements were made in simulating coupled vibration - dissociation - recombination for post-shock flows. Models were also developed to compute vehicle surface temperatures directly in the code rather than requiring isothermal estimates. These codes were instrumental in simulating aerobraking of NASA's Magellan spacecraft during orbital maneuvers to assess heat transfer and aerodynamic properties of the delicate satellite. NASA also depended upon simulations of entry of the Galileo probe into the atmosphere of Jupiter to provide drag and flow field information essential for accurate interpretation of an onboard experiment. Finally, the codes have been used extensively to simulate expanding nozzle flows in low-power thrusters in support of propulsion activities at NASA-Lewis. Detailed comparisons between continuum calculations and DSMC results helped to quantify the limitations of continuum CFD codes in rarefied applications.

Boyd, Iain↗
DOE OSTI2025

Alquimia v1.0: a generic interface to biogeochemical codes – a tool for interoperable development, prototyping and benchmarking for multiphysics simulators

Alquimia v1.0 is a generic interface to geochemical solvers that facilitates development of multiphysics simulators by enabling code coupling, prototyping and benchmarking. The interface enforces the function arguments and their types for setting up, solving, serving up output data and carrying out other common auxiliary tasks while providing a set of structures for data transfer between the multiphysics code driving the simulation and the geochemical solver. Alquimia relies on a single-cell approach that permits operator splitting coupling and parallel computation. We describe the implementation in Alquimia of two widely used open-source codes that perform geochemical calculations: PFLOTRAN and CrunchFlow. We then exemplify its use for the implementation and simulation of reactive transport in porous media by two open-source flow and transport simulators: Amanzi and ParFlow. We also demonstrate its use for the simulation of coupled processes in novel multiphysics applications including the effect of multiphase flow on reaction rates at the pore scale with OpenFOAM, the role of complex biogeochemical processes in land surface models such as the E3SM Land Model (ELM) and the impact of surface–subsurface hydrological interactions on hydrogeochemical export from watersheds with the Advanced Terrestrial Simulator (ATS). These applications make it apparent that the availability of a well-defined yet flexible interface has the potential to improve the software development workflow, freeing up resources to focus on advances in process models and mechanistic understanding of coupled problems.

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

Convection-permitting hindcasting of diurnal variation of Mei-yu rainfall over East China with a global variable-resolution model

Two sets of 48-hr hindcast experiments for Mei-yu season using non-hydrostatic global variable-resolution model (MPAS-Atmosphere) with regional refinement at 4km resolution (V4km) and uniform 60km resolution (U60km) are first-time conducted to attempt at simulating the Mei-yu rainfall in 2015 over East China. The analysis focuses on the seasonal average of diurnal variation of precipitation amount (PA), intensity (PI), and frequency (PF). Both simulations reasonably reproduce the spatial distribution of PA. V4km is more skillful in simulating the spatial distributions and magnitudes of PF and PI. The diurnal cycle of Mei-yu rainfall shows a major early morning peak and a minor afternoon peak, contributed by precipitation during two sub-periods with distinct synoptic circulations. With strong Mei-yu in the first sub-period, the diurnal variation of PA is controlled by nocturnal southwesterly jet. V4km overestimates the morning peak mainly due to its bias in simulating boundary layer inertial oscillation. Although differences in PA between the simulations are small, U60km overestimates PF and underestimates PI. With weak Mei-yu in the second sub-period, the diurnal variation of PA is controlled by both synoptic circulation and local convection. At both resolutions, deviations in the large-scale circulation modulated by a few typhoons lead to positive biases in the morning peak of PA. After removing the typhoon impacts, V4km captures the observed diurnal cycle of PA well, while U60km significantly underestimates PA and PI particularly in the afternoon. Future studies focusing on advancing modeling of southwesterly jet and typhoons may further improve convection permitting simulation of Mei-yu rainfall.

Xu, Mingyue↗
DOE OSTI2017

Generation of Subsurface Voids, Incubation Effect, and Formation of Nanoparticles in Short Pulse Laser Interactions with Bulk Metal Targets in Liquid: Molecular Dynamics Study

The ability of short pulse laser ablation in liquids to produce clean colloidal nanoparticles and unusual surface morphology has been employed in a broad range of practical applications. In this paper, we report the results of large-scale molecular dynamics simulations aimed at revealing the key processes that control the surface morphology and nanoparticle size distributions by pulsed laser ablation in liquids. The simulations of bulk Ag targets irradiated in water are performed with an advanced computational model combining a coarse-grained representation of liquid environment and an atomistic description of laser interaction with metal targets. For the irradiation conditions that correspond to the spallation regime in vacuum, the simulations predict that the water environment can prevent the complete separation of the spalled layer from the target, leading to the formation of large subsurface voids stabilized by rapid cooling and solidification. The subsequent irradiation of the laser-modified surface is found to result in a more efficient ablation and nanoparticle generation, thus suggesting the possibility of the incubation effect in multipulse laser ablation in liquids. The simulations performed at higher laser fluences that correspond to the phase explosion regime in vacuum reveal the accumulation of the ablation plume at the interface with the water environment and the formation of a hot metal layer. The water in contact with the metal layer is brought to the supercritical state and provides an environment suitable for nucleation and growth of small metal nanoparticles from metal atoms emitted from the hot metal layer. The metal layer itself has limited stability and can readily disintegrate into large (tens of nanometers) nanoparticles. The layer disintegration is facilitated by the Rayleigh–Taylor instability of the interface between the higher density metal layer decelerated by the pressure from the lighter supercritical water. The nanoparticles emerging from the layer disintegration are rapidly cooled and solidified due to the interaction with water environment, with a cooling rate of ~2 × 10 12 K/s observed in the simulations. The computational prediction of two distinct mechanisms of nanoparticle formation yielding nanoparticles with different characteristic sizes provides a plausible explanation for the experimental observations of bimodal nanoparticle size distributions in laser ablation in liquids. The ultrahigh cooling and solidification rates suggest the possibility for generation of nanoparticles featuring metastable phases and highly nonequilibrium structures.

Chemistry↗
DOE OSTI2017

Advanced Fluid Reduced Order Models for Compressible Flow.

This report summarizes fiscal year (FY) 2017 progress towards developing and implementing within the SPARC in-house finite volume flow solver advanced fluid reduced order models (ROMs) for compressible captive-carriage flow problems of interest to Sandia National Laboratories for the design and qualification of nuclear weapons components. The proposed projection-based model order reduction (MOR) approach, known as the Proper Orthogonal Decomposition (POD)/Least- Squares Petrov-Galerkin (LSPG) method, can substantially reduce the CPU-time requirement for these simulations, thereby enabling advanced analyses such as uncertainty quantification and de- sign optimization. Following a description of the project objectives and FY17 targets, we overview briefly the POD/LSPG approach to model reduction implemented within SPARC . We then study the viability of these ROMs for long-time predictive simulations in the context of a two-dimensional viscous laminar cavity problem, and describe some FY17 enhancements to the proposed model reduction methodology that led to ROMs with improved predictive capabilities. Also described in this report are some FY17 efforts pursued in parallel to the primary objective of determining whether the ROMs in SPARC are viable for the targeted application. These include the implemen- tation and verification of some higher-order finite volume discretization methods within SPARC (towards using the code to study the viability of ROMs on three-dimensional cavity problems) and a novel structure-preserving constrained POD/LSPG formulation that can improve the accuracy of projection-based reduced order models. We conclude the report by summarizing the key takeaways from our FY17 findings, and providing some perspectives for future work.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

CFD Simulations of Lower Plenum Mixing

Review of model development and validation performed in the Advanced Reactor Technologies (ART) program for thermal mixing at the outlet of High Temperature Gas Reactors (HTGRs). Understanding the mixing that occurs in the lower plenum in an HTGR is necessary to facilitate design improvements and to perform reactor safety analysis. Numerical models are one possible approach to gain a better understanding of mixing in the lower plenum. Given the complexity of the geometry and the intense mixing present, it is important to perform validation of numerical models. Three models have been developed during FY2025: a porous media with Pronghorn, a Reynolds Averaged Navier Stokes (RANS) with STAR-CCM+, and a Large Eddy Simulation (LES) with NekRS. The reference facility is a scaled-down version of the lower plenum of the High Temperature Gas-Cooled Reactor - Pebble-bed Module (HTR-PM) demonstration reactor. Preliminary results of the porous media and the RANS shows general good agreement against experimental benchmark data. Future work will leverage high-fidelity results obtained through LES to guide model selection and improvements to the lower-fidelity models, with particular attention to the Pronghorn porous media.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

WellPINN: Accurate Well Representation for Transient Fluid Pressure Diffusion in Subsurface Reservoirs With Physics‐Informed Neural Networks

Accurate representation of pumping wells is essential for reliable reservoir characterization and simulation of operational scenarios in subsurface flow models. Physics-informed neural networks (PINNs) are emerging as a promising alternative to numerical models for reservoir modeling, offering seamless integration of monitoring data and governing physical equations. However, existing PINN-based studies face major challenges in capturing fluid pressure near wells when using a source/sink term, particularly during the early stages after pumping begins. We address this problem by introducing WellPINN, a workflow in which an initially trained PINN infers fluid pressure across the entire reservoir domain using a large equivalent well radius. This initial PINN solution is then locally refined around the well by a set of subdomain PINNs that are trained for smaller equivalent well radii. Continuity across these subdomain interfaces as well as at the initial condition is ensured by hard-constraining each PINN on its subdomain boundary. Our results demonstrate WellPINN as the first workflow of its kind to focus on accurate inference of fluid pressure from pumping rates throughout the entire injection period, significantly advancing the potential of PINNs for inverse modeling and operational scenario simulations. All data and code for this paper are openly available at https://doi.org/10.20350/DIGITALCSIC/17260.

58 GEOSCIENCES↗
NASA NTRS2021

Hybrid phenology matching model for robust crop phenological retrieval

Crop phenology regulates seasonal agroecosystem carbon, water, and energy exchanges, and is a key component in empirical and process-based crop models for simulating biogeochemical cycles of farmlands, assessing gross and net primary production, and forecasting the crop yield. The advances in phenology matching models provide a feasible means to monitor crop phenological progress using remote sensing observations, with a priori information of reference shapes and reference phenological transition dates. Yet the underlying geometrical scaling assumption of models, together with the challenge in defining phenological references, hinders the applicability of phenology matching in crop phenological studies. The objective of this study is to develop a novel hybrid phenology matching model to robustly retrieve a diverse spectrum of crop phenological stages using satellite time series. The devised hybrid model leverages the complementary strengths of phenometric extraction methods and phenology matching models. It relaxes the geometrical scaling assumption and can characterize key phenological stages of crop cycles, ranging from farming practice-relevant stages (e.g., planted and harvested) to crop development stages (e.g., emerged and mature). To systematically evaluate the influence of phenological references on phenology matching, four representative phenological reference scenarios under varying levels of phenological calibrations in terms of time and space are further designed with publicly accessible phenological information. The results indicate that the hybrid phenology matching model can achieve high accuracies for estimating corn and soybean phenological growth stages in Illinois, particularly with the year- and region-adjusted phenological reference (R-squared higher than 0.9 and RMSE less than 5 days for most phenological stages). The inter-annual and regional phenological patterns characterized by the hybrid model correspond well with those in the crop progress reports (CPRs) from the USDA National Agricultural Statistics Service (NASS). Compared to the benchmark phenology matching model, the hybrid model is more robust to the decreasing levels of phenological reference calibrations, and is particularly advantageous in retrieving crop early phenological stages (e.g., planted and emerged stages) when the phenological reference information is limited. This innovative hybrid phenology matching model, together with CPR-enabled phenological reference calibrations, holds 3 considerable promise in revealing spatio-temporal patterns of crop phenology over extended geographical regions.

Phenology↗
DOE OSTI2021

A modeling framework for designing and evaluating curbside traffic management policies at Dallas-Fort Worth International Airport

Emerging mobility technologies are changing the transportation system landscape. This is especially evident at airports, such as the Dallas-Fort Worth International Airport (DFW). Without careful analysis, these changes could lead to inefficient and costly airport operations. This paper presents a modeling framework that integrates travel mode encoding, demand projection, and microsimulation to enable airports to develop, simulate, and evaluate curbside traffic managements policies and measure their impact. Here, the framework is utilized to analyze several traffic scenarios and policies for DFW: a baseline scenario which represents DFW traffic pattern as observed in 2018 and projected to 2045, a transit network company (TNC) electrification policy, a TNC queuing policy, a policy that increased transit ridership, a bus-only policy which considers the use of only buses inside DFW, an autonomous vehicle (AV) policy which investigates the impact of autonomous vehicle (AV) adoption on airport operations, and an example COVID-19 scenario which models the impact of the COVID19 pandemic. The simulations’ results demonstrate that: increasing the DFW transit ridership postpones the need for airport curbside expansion the most; encouraging shared-mobility with the bus-only policy produces the most savings in curbside congestion delays; automation and electrification for all passenger vehicle trips to/from DFW generates the most saving in fuel consumption and emissions; and uncontrolled AV adoption incurs the highest increase in fuel consumption, delay, and emissions and could require immediate airport capacity extension. Without policy intervention or investment in additional infrastructure capacity, these results predict the current operations would face significant congestion on high demand days starting as early as 2028. While derived in close partnership with DFW, the methodology presented here can be generalized to any airport.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRSDate not supplied

Machine Learning for NASA Advanced Information Systems

NASA's Advanced Information Systems Technology (AIST) Program is one of several Technology programs managed by the Earth Science Technology Office (ESTO) in the Earth Science Division (ESD). The AIST Program focuses on advanced information systems and novel computer science technologies that will be needed by NASA Earth Science in the next 5 to 10 years. The three main thrusts of the AIST Program deal with Novel Observing Strategies (NOS), Analytic Collaborative Frameworks (ACF) and Earth System Digital Twins (ESDT). For all these thrusts, Machine Learning (ML) is increasingly being used in multiple aspects of Earth science systems, e.g., for onboard autonomy and decision making, for the analysis of massive and diverse datasets as well as more recently for developing surrogate models that will represent one of the main components of future Digital Twins of the Earth. Particularly, ESDT technologies developed by the AIST Program will allow to develop integrated Earth Science frameworks that will mirror the Earth with state-of-the-art models (Earth system models and others), timely and relevant observations, and analytic tools. These information systems will be used for supporting near- and long-term science and policy decisions. ESDT frameworks will build on previously developed AIST capabilities and technologies to integrate interconnected models with continuous streams of observations, data analytics, data assimilation, simulations, advanced visualizations and the ability to conduct "what-if" scenarios. This talk will describe the three thrusts of the AIST Program with a special focus on Machine Learning and how it is being used at all steps of the Earth Science data lifecycle.

Mathematical and Computer Sciences (General)↗
NASA NTRS2009

Scientific Benefits of Space Science Models Archiving at Community Coordinated Modeling Center

The Community Coordinated Modeling Center (CCMC) hosts a set of state-of-the-art space science models ranging from the solar atmosphere to the Earth's upper atmosphere. CCMC provides a web-based Run-on-Request system, by which the interested scientist can request simulations for a broad range of space science problems. To allow the models to be driven by data relevant to particular events CCMC developed a tool that automatically downloads data from data archives and transform them to required formats. CCMC also provides a tailored web-based visualization interface for the model output, as well as the capability to download the simulation output in portable format. CCMC offers a variety of visualization and output analysis tools to aid scientists in interpretation of simulation results. During eight years since the Run-on-request system became available the CCMC archived the results of almost 3000 runs that are covering significant space weather events and time intervals of interest identified by the community. The simulation results archived at CCMC also include a library of general purpose runs with modeled conditions that are used for education and research. Archiving results of simulations performed in support of several Modeling Challenges helps to evaluate the progress in space weather modeling over time. We will highlight the scientific benefits of CCMC space science model archive and discuss plans for further development of advanced methods to interact with simulation results.

Kuznetsova, Maria M.↗
NASA NTRS2008

Arctic Mixed-phase Clouds Simulated by a Cloud-Resolving Model: Comparison with ARM Observations and Sensitivity to Microphysics Parameterizations

Single-layer mixed-phase stratiform (MPS) Arctic clouds, which formed under conditions of large surface heat flux combined with general subsidence during a subperiod of the Atmospheric Radiation Measurement (ARM) Program Mixed-Phase Arctic Cloud Experiment (M-PACE), are simulated with a cloud resolving model (CRM). The CRM is implemented with either an advanced two-moment (M05) or a commonly used one-moment (L83) bulk microphysics scheme and a state-of-the-art radiative transfer scheme. The CONTROL simulation, that uses the M05 scheme and observed aerosol size distribution and ice nulei (IN) number concentration, reproduces the magnitudes and vertical structures of cloud liquid water content (LWC), total ice water content (IWC), number concentration and effective radius of cloud droplets as suggested by the M-PACE observations. It underestimates ice crystal number concentrations by an order of magnitude and overestimates effective radius of ice crystals by a factor of 2-3. The OneM experiment, that uses the L83 scheme, produces values of liquid water path (LWP) and ice plus snow water path (ISWP) that were about 30% and 4 times, respectively, of those produced by the CONTROL. Its vertical profile of IWC exhibits a bimodal distribution in contrast to the constant distribution of IWC produced in the CONTROL and observations.

Xu, Kuan-Man↗
DOE OSTI2023

Simulation-based characterization of the variability of earthquake risk to buildings in the near-field

Recent advancements in high performance computing platforms and computational workflow for regional-scale simulations are enabling unprecedented modeling of fault-to-structure earthquake processes. Regional simulations resolving ground motions at frequencies relevant to engineered systems are becoming computationally viable and provide a new capability to improve understanding of the geographical distribution and intensity of risk to buildings and critical infrastructure. As computational capabilities advance, it is essential to move beyond illustrative single rupture realizations for scenario earthquake events towards the development of a full suite of rupture realizations that appropriately characterize the range of risk to building systems. The work described in this article investigates the application of a suite of fault rupture realizations with the objective of assessing near-fault, site-specific seismic demand variability for building structures. A representative high-performance regional-scale computational model is utilized to execute ground motion and building response simulations based on 18 kinematic rupture realizations of an M7 strike-slip scenario earthquake. The fault rupture models for the scenario earthquake are created by systematically perturbing the hypocenter location and stochastically generating rupture parameters (slip, rise time, rake angle) to represent a breadth of ground motion intensities resulting from the spatial and temporal variabilities of an earthquake rupture process. The resulting seismic demand variability for three-story (short period) and forty-story (long period) steel moment-resisting frame buildings is characterized in terms of the median and distribution of peak inter-story drift ratio for a range of near-fault sites. The full suite of 18 fault rupture realizations and approximately 280,000 nonlinear dynamic building simulations indicate that the three-story building undergoes higher median seismic demand and significantly greater variability of demand at a given site than the forty-story building, which has important implications for the level of certainty in predicting building performance during an earthquake. The simulations performed provide deeper insight into the relationship between fault rupture parameterization and building response, which is essential information for developing a representative suite of rupture realizations for specific earthquake scenarios.

58 GEOSCIENCES↗
DOE OSTI2017

Combustion energy frontier research center (CEFRC) final report (August 1, 2009 – July 31, 2016)

The Combustion Energy Frontier Research Center (CEFRC) was established to tackle the single overarching grand challenge of energy sustainability, energy security and global warming: to develop a “validated, predictive, multi-scale, combustion modeling capability to optimize the design and operation of evolving fuels in advanced engines for transportation applications,” as identified in the DOE report on “Basic Energy Needs for Clean and Efficient Combustion of 21st Century Transportation Fuels”. The challenge is particularly daunting since energy conversion efficiencies and exhaust emissions are governed by coupled chemical and transport processes at multiple length scales ranging from electron excitation to molecular rearrangements to nanoscale particulate formation to turbulent fuel/air mixing. To tackle this challenge, the CEFRC assembled a world-class team of 15 principal investigators, with the objectives to: 1) develop and test theoretical models to predict elementary reaction rates, molecule thermalization rates, chemical bond dissociation energies, and nonequilibrium transport properties using quantum chemistry calculations that account for strong electron correlation and multiple electronic potential energy surfaces; 2) develop automated kinetic mechanism generation, reduction, and error control methods for predicting alternative fuel including biofuel oxidation, fuel droplet decomposition, and NOx and particulate formation; 3) validate and improve the predictions of these models by measuring ignition delay times, species profiles, flame structures, burning limits, turbulence-kinetic coupling, and NOx and soot emissions at high-pressures and near-limit conditions, by using advanced experimental diagnostic techniques including multiple laser techniques, molecular beam sampling and synchrotron photoionization, and by conducting the measurements in high-pressure shock tubes, jet-stirred and flow reactors, flame bombs, counterflow flames, and advanced-design rapid compression ignition instruments; and 4) develop a suite of validated petascale high-fidelity simulation and modeling capabilities to understand and predict chemistry-turbulence-radiation coupling for new fuels in new regimes, including the high pressure, low-temperature combustion in advanced engine and turbine designs, and 5) establish a knowledge highway between researchers and engineers in academia, national laboratories, and industry to facilitate the dissemination and exchange of knowledge on national and international levels, and enrich the talent pool and capabilities of the next generation of combustion scientists and engineers. The technical activities of the CEFRC were conducted through three Disciplinary Working Groups – Chemistry Theory, Experiment and Mechanism, and Reacting Flows, which coordinated the Center’s research on the development of combustion chemistry of Foundation Fuels (C0–C4 hydrocarbons), Alcohols, and Biodiesel through three corresponding Mechanism Thrust Groups. Such a two-dimensional coordinated and tightly interwoven research structure has been proven to be highly effective in assuring the interplay between the developments of the fundamentals of combustion science and the utilization of the various categories of fuels. The Center has accomplished the above goals over the five year period (August 1, 2009 – July 31, 2014) with appropriated funding, followed by two additional no-cost-extension (NCE) years. The research results are documented in 230 journal publications, with six legacy review papers on the study of combustion chemistry using shock tubes, flow reactors, rapid compression machines, and flames, on uncertainty quantification, and on theoretical reaction dynamics and chemical modeling of combustion. A robust outreach program complemented these PI-led research activities, consisting of: 1) a roving post-doc program comprised of a corps of Center-appointed, co- or multi-sponsored post-doctoral fellows with rotating assignments to conduct seed projects initiated by at least two PIs, in residence with these sponsoring PIs, to rapidly pursue new and high-risk, high-payoff interdisciplinary ideas; 2) an annual summer school on combustion heavily attended (~200) by senior graduate students and practicing researchers covering advanced topics on chemical kinetics, fluid mechanics, turbulent combustion, engine combustion, new technologies, etc.; 3) a robust open web-site providing Center and community information as well as the lecture videos and notes of the summer school; and 4) widely distributed biannual newsletters.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Coupling Surface Flow with High-performance Subsurface Reactive Flow and Transport Code PFLOTRAN

Water exchange between the surface and subsurface is important for both water resource management and environmental protection. In this paper, we develop coupled surface and subsurface flow simulation capability in a parallel subsurface flow and reactive transport code PFLOTRAN. We sequentially couple the diffusion wave-based surface flow with the subsurface flow governedby the Richards equation in PFLOTRAN. These two flow domains are linked with a boundary condition switching method that ensures continuity of pressure and flux at the surface-subsurface interface. We verify the coupled code against other existing hydrologic models and observation data using a number of numerical experiments. The coupled hydrological model exhibits good performance in strong parallel scaling tests. The new coupled surface and subsurface simulator significantly advance community simulation capability towards improving integrated hydrologic and biogeochemical understanding of complex systems such as watersheds and river corridors. Keywords: Surface flow, Integrated hydrological modeling, Boundary condition switching, Parallel computing

Wu, Runjian↗