Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model Based Systems Engineering”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Sensor System and Observer Algorithm Co-Design For Modern Internal Combustion Engine Air Management Based on H2 Optimization

This paper outlines a novel sensor selection and observer design algorithm for linear time-invariant systems with both process and measurement noise based on H 2 optimization to optimize the tradeoff between the observer error and the number of required sensors. The optimization problem is relaxed to a sequence of convex optimization problems that minimize the cost function consisting of the H 2 norm of the observer error and the weighted l 1 norm of the observer gain. An LMI formulation allows for efficient solution via semi-definite programing. The approach is applied here, for the first time, to a turbo-charged spark-ignited engine using exhaust gas circulation to determine the optimal sensor sets for real-time intake manifold burnt gas mass fraction estimation. Simulation with the candidate estimator embedded in a high fidelity engine GT-Power model demonstrates that the optimal sensor sets selected using this algorithm have the best H 2 estimation performance. Sensor redundancy is also analyzed based on the algorithm results. This algorithm is applicable for any type of modern internal combustion engines to reduce system design time and experimental efforts typically required for selecting optimal sensor sets.

Zhang, Xu↗
DOE OSTI2019

Making Coal Relevant for Small Scale Applications: Modular Gasification for Syngas/Engine CHP Applications in Challenging Environments

This Front-End Engineering Design (FEED) study is funded by the US Department of Energy (DOE), and managed by the National Energy Technology Laboratory (NETL) through a Cooperative Agreement awarded to the University of Alaska Fairbanks (UAF) to perform a technical and cost study for the design and installation of a Hamilton-Maurer International (HMI) gasifier for generating a clean syngas for firing in reciprocating engines. An existing diesel engine generator (DEG) and a commercially available spark ignition engine are used to generate electricity with the low-Btu syngas. Other than the HMI gasifier, only commercially available technologies were assumed for this study. An important objective of this design study was to estimate capital and operating costs of small modular coal and coal/biomass fired generating units to provide wind regulation services. Although this study is based on retrofitting an existing coal fired power plant with the HMI gasifier, sufficient data was developed to characterize greenfield units. The original concept was to provide clean syngas and tars/oils from the gasification process to the existing Colt-Pielstick PC2.6B diesel engine generator (9.6 MWe) installed at the UAF Atkinson Power Plant. The intent was to convert the diesel engine to operate in a dual fuel mode, firing both syngas and tars/oils. As the project progressed, Fairbanks Morse Engine (FME), current licensee for the Colt-Pielstick engine line, was not able to develop a modification package for the engine to allow for syngas firing. Although this type of modification is available for firing natural gas, the capability does not extend to syngas. Western Energy Systems, distributor for Innio Jenbacher, was able to offer the Jenbacher Model JMS620 engine. After review of available options, the project team selected a design configuration using two Jenbacher engines rated at 1965 kWe each to fire the syngas, with the oils/tars directed to the Colt Pielstick diesel to be blended with ultra-low sulfur diesel oil (ULSD). The Jenbacher engines are expected to operate at continuous base load, with some load-following ability, while the existing diesel engine will operate in grid support service, with loads varying from a nominal 25% to 100%. An annual capacity factor for the diesel engine is projected to be between 50% and 60%. The overall gasifier and gas clean-up system is sized to produce enough syngas to adequately fuel three Jenbacher Engines for a nominal power output of 5.9 MWe, but this study assumes that only two Jenbacher engines will be installed due to space constraints in the Atkinson Power Plant, and because of expected air emission permit limitations. There is still potential that the excess syngas capacity equivalent to about 2 MWe could be used as supplemental fuel for the DEG. The exhaust from the Jenbacher engines is passed through heat recovery units to generate low pressure steam, thereby providing a significant portion of the steam required for operating the gasifier. The overall performance of the HMI gasifier/engine configuration is summarized as follows: 1) Overall electric efficiency, Jenbacher and diesel engines, based on gasifier fuel production (syngas and pyrolysis liquids) is 35.2% LHV, 2) Overall cogenerating efficiency (electricity plus steam), based on gasifier fuel production (syngas and pyrolysis liquids) is 54.8% LHV, 3) Overall electric efficiency, Jenbacher and diesel engines, based on coal input, is 29.6% LHV, and 4) Overall cogenerating efficiency, electricity plus steam, based on coal input is 46.0%. This project is estimated to have a Bare Erected Cost (BEC) of $\$$ 33.4 million and a Total Plant Cost (TPC) of $\$$ 45.7 million. Based on the original scope, sending both the clean syngas and the oils and tars to the existing diesel engine, the BEC and TPC values would have been $\$$ 23.3 million and $\$$ 32.1 million, respectively. A variable Cost of Electricity (COE) is calculated to reflect only those costs that will be incurred by incremental operation of the system for the HMI gasifier/Jenbacher engine combination. Operation of the diesel engine is determined separately. The variable COE values are: Case 1: 100% Usibelli Coal: (a) HMI gasifier/Jenbacher Engines -- 8.3 cents/kWh; (b) Diesel Engine -- 20.2 cents/kWh, and Case 2: 80% Usibelli Coal + 20% wood chips (proportion based on Btu input, LHV basis): (a) HMI Gasifier/Jenbacher -- 13.7 cents/kWh; (b) Diesel Engine -- 20.2 cents/kWh. The incremental cost of operating the gasifier/Jenbacher system increases with co-firing of wood chips due to the high cost of this material, over five times as costly as Usibelli coal delivered to the Atkinson Power building. The diesel cost is constant as the oils and tars are priced at the same value as the ULSD for the purposes of fueling the diesel engine. Section 8 of this report provides details of the O&M costs. The overall power system is configured to be modular and scalable to provide a wide range of generating capacity using the same components. This allows for reduction in cost, risk, schedule time to engineer and construct, and provides some degree of repeatability for air permitting.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Graph-Based Modeling for the Detection and Tracking of Sarin-Surrogate-Induced Neurotoxicity Using a Human-Relevant, In-Vitro Brain Model

Organophosphorus (OP) nerve agents are a chemical threat to the United States, to the civilian population (e.g., pesticides) and historically weaponized (e.g., sarin) as chemical warfare agents. The unprecedented, accelerated process from “bench-to-bedside” during the SARSCov2 pandemic has made it clear that technology and tools need to be readily available for immediate response. Advances in human organ tissue mimetic systems are a promising technology to evaluate the human-relevant response in vitro for basic and applied research and drug screening. In particular, current brain microphysiological systems (MPS) have the capability to monitor and detect changes in engineered human neural circuit activity. However, current data analytics approaches for these systems lack the granularity to functionally detect and distinguish the different mechanisms that occur in the brain following neurotoxicity, injury, and disease. The goal of this project was to advance the computational analytical capabilities of the brain MPS to detect functional changes in neural circuit structure at different stages of Sarin surrogate-induced neurotoxicity. We developed graph-based models to (1) identify the composition of the neural circuit structure; (2) detect and monitor how this structure changes following sarin-induced neurotoxicity; and (3) evaluate the analytical pipeline using known/promising oxime reactivators. Through experiments on the bMPS where in vitro neuronal cultures were exposed to a sarin surrogate, we demonstrated the capabilities of our computational pipeline to identify different responses in the functional networks of brain cells exposed to low and high concentrations of the nerve agent. We identified a biphasic response of human neural network activity following exposure to a sarin-surrogate that had not been reported in the literature before. The graph-based models and software developed in this project can be used for future studies that leverage the brain MPS technology, such as treatment efficacy assessment.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Reinforcement learning for adaptive maintenance policy optimization under imperfect knowledge of the system degradation model and partial observability of system states

Maintenance policy optimization usually is faced with challenges that arise from an imperfect knowledge of system degradation models and from the partial observability of system degradation states. Here, this paper proposes a reinforcement learning method to address these two challenges for a class of maintenance problems with Markov degradation processes. The reinforcement learning approach consists of a learning component and a planning component. Using sequentially collected observations, at each step of decision-making the learning component improves the knowledge of system degradation in terms of the probability distributions of the transition rates based on sequential Bayesian inference. Using the updated transition rates, at each step of decision-making the maintenance policy optimization problem is then formulated as a partially observable Markov decision problem, and the planning component computes the optimal maintenance policy that maximizes the expected cumulative reward. The proposed method is illustrated using a numerical example with repair and inspection maintenance actions. The result shows that as more observations are collected, the learning component progressively learns the true system degradation process, and the planning component adjusts the optimal maintenance policy accordingly as well, which leads to increased reward.

42 ENGINEERING↗
DOE OSTI2015

Experimental Systems-Biology Approaches for Clostridia-Based Bioenergy Production

This is the final project report for project "Experimental Systems-Biology Approaches for Clostridia-Based Bioenergy Production" for the funding period of 9/1/12 to 2/28/2015 (three years with a 6-month no-cost extension) OVERVIEW AND PROJECT GOALS The bottleneck of achieving higher rates and titers of toxic metabolites (such as solvents and carboxylic acids that can used as biofuels or biofuel precursors) can be overcome by engineering the stress response system. Thus, understanding and modeling the response of cells to toxic metabolites is a problem of great fundamental and practical significance. In this project, our goal is to dissect at the molecular systems level and build models (conceptual and quantitative) for the stress response of C. acetobutylicum (Cac) to its two toxic metabolites: butanol (BuOH) and butyrate (BA). Transcriptional (RNAseq and microarray based), proteomic and fluxomic data and their analysis are key requirements for this goal. Transcriptional data from mid-exponential cultures of Cac under 4 different levels of BuOH and BA stress was obtained using both microarrays (Papoutsakis group) and deep sequencing (RNAseq; Meyers and Papoutsakis groups). These two sets of data do not only serve to validate each other, but are also used for identification of stress-induced changes in transcript levels, small regulatory RNAs, & in transcriptional start sites. Quantitative proteomic data (Lee group), collected using the iTRAQ technology, are essential for understanding of protein levels and turnover under stress and the various protein-protein interactions that orchestrate the stress response. Metabolic flux changes (Antoniewicz group) of core pathways, which provide important information on the re-allocation of energy and carbon resources under metabolite stress, were examined using 13C-labelled chemicals. Omics data are integrated at different levels and scales. At the metabolic-pathway level, omics data are integrated into a 2nd generation genome-scale model (GSM) (Maranas group). Omics data are also integrated using bioinformatics (Wu and Huang group), whereby regulatory details of gene and protein expression, protein-protein interactions and metabolic flux regulation are incorporated. The PI (Papoutsakis) facilitated project integration through monthly meeting and reports, conference calls, and collaborative manuscript preparation. The five groups collaborated extensively and made a large number of presentations in national and international meetings. It has also published several papers, with several more in the preparation stage. Several PhD, MS and postdoctoral students were trained as part of this collaborative and interdisciplinary project.

09 BIOMASS FUELS↗
DOE OSTI2025

Modelling pulsed field magnetization of iron-based bulk superconductors

Abstract Bulk superconductors can be used as super-strength quasi-permanent magnets capable of providing magnetic flux densities considerably superior to conventional permanent magnets. This makes them attractive for several engineering applications that rely on strong magnetic fields like rotating machines, NMR/MRI and magnetic drug delivery systems. Recently, the authors reported a record trapped magnetic field in an iron-based bulk superconductor: 2.83 T was trapped in potassium-doped barium iron arsenide (Ba, K)Fe 2 As 2 (or Ba122) at 5 K. Of particular significance is that the strength and temporal stability of this magnetic field exceeds the requirements of MRI machines, indicating iron-based bulks can now perform at levels demanded by engineering applications. One crucial challenge for their practical use, however, is the need to apply and remove an external magnetic field to magnetize them. Pulsed field magnetization (PFM) shows great promise as a practical method of magnetizing bulks, but the process generates heat in the bulk that is detrimental to its superconducting performance and ability to act as a super-strength magnet. In this paper, coupled electromagnetic–thermal numerical models are used to simulate the PFM of iron-based bulk superconductors. Here we focus on the recent-record-breaking, fine-grain polycrystalline K-doped Ba122 bulks. The impact that the specific J c ( B ) characteristics and thermal properties of the Ba122 material—all of which have been experimentally measured from state-of-the-art samples—have on the magnetic flux dynamics and thermal behaviour during PFM, including the final trapped field, is investigated. We show that because the thermal properties are similar to those of REBa 2 Cu 3 O 7 −δ bulks, a similar response to pulsed fields is obtained. A maximum trapped field of ∼0.81 T (∼43.3% of the maximum trapped field capability under ideal, field-cooling conditions) was simulated at 5 K, with a magnetization efficiency of ∼54%. The modelling framework provides a fast and flexible tool for optimising the practical PFM process at different operating temperatures to maximise the trapped field in state-of-the-art Ba122 bulks and to guide the design of future experiments.

bulk superconductors↗
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 OSTI2015

On the impact of temperatures up to 200°C in clay repositories with bentonite engineer barrier systems: A study with coupled thermal, hydrological, chemical, and mechanical modeling

One of the most important design variables for a geological nuclear waste repository is the temperature limit up to which the engineered barrier system (EBS) and the natural geologic environment can be exposed. Up to now, almost all design concepts that involve bentonite-backfilled emplacement tunnels have chosen a maximum allowable temperature of about 100°C. Such a choice is largely based on the consideration that in clay-based materials illitization and the associated mechanical changes in the bentonite (and perhaps the clay host rock) could affect the barrier attributes of the EBS. However, existing experimental and modeling studies on the occurrence of illitization and related performance impacts are not conclusive, in part because the relevant couplings between the thermal, hydrological, chemical, and mechanical (THMC) processes have not been fully represented in the models. This work presents a fully coupled THMC simulation of a nuclear waste repository in a clay formation with a bentonite-backfilled EBS for 1000 years. Two scenarios were simulated for comparison: a case in which the temperature in the bentonite near the waste canister can reach about 200°C and a case in which the temperature in the bentonite near the waste canister peaks at about 100°C. The model simulations demonstrate some degree of illitization in both the bentonite buffer and the surrounding clay formation. Other chemical alterations include the dissolution of K-feldspar and calcite, and precipitation of quartz, chlorite, and kaolinite. In general, illitization in the bentonite and the clay formation is enhanced at higher temperature. However, the quantity of illitization is affected by many chemical factors and therefore varies a great deal. The most important chemical factors are the concentration of K in the pore water as well as the abundance and dissolution rate of K-feldspar; less important ones are the concentration of sodium and the quartz precipitation rate. In our modeling scenarios, the calculated decrease in smectite volume fraction in bentonite ranges from 1 to 8% of the initial volume fraction of smectite in the 100°C scenario and 1–27% in the 200°C scenario. Chemical changes in the 200°C scenario could also lead to a reduction in swelling stress up to 15–18% whereas those in the 100°C scenario result in about 14–15% reduction in swelling stress for the base case scenario. Model results also show that the 200°C scenario results in a much higher total stress than the 100°C scenario, mostly due to thermal pressurization. While cautions should be taken regarding the model results due to some limitations in the models, the modeling work is illustrative in light of the relative importance of different processes occurring in EBS bentonite and clay formation at higher than 100°C conditions, and could be of greater use when site specific data are available.

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

Floquet Engineering of Interactions and Entanglement in Periodically Driven Rydberg Chains

Neutral atom arrays driven into Rydberg states constitute a promising approach for realizing programmable quantum systems. Enabled by strong interactions associated with Rydberg blockade, they allow for simulation of complex spin models and quantum dynamics. We introduce a new Floquet engineering technique for systems in the blockade regime that provides control over novel forms of interactions and entanglement dynamics in such systems. Our approach is based on time-dependent control of Rydberg laser detuning and leverages perturbations around periodic many-body trajectories as resources for operator spreading. These time-evolved operators are utilized as a basis for engineering interactions in the effective Hamiltonian describing the stroboscopic evolution. As an example, we show how our method can be used to engineer strong spin exchange, consistent with the blockade, in a one-dimensional chain, enabling the exploration of gapless Luttinger liquid phases. In addition, we demonstrate that combining gapless excitations with Rydberg blockade can lead to dynamic generation of large-scale multipartite entanglement. Experimental feasibility and possible generalizations are discussed.

Floquet systems↗
DOE OSTI2025

Towards philosophical reasoning with agentic LLMs: Socratic method for scientific assistance

As large language models (LLMs) become central tools in science, improving their reasoning capabilities is critical for meaningful and trustworthy applications. We introduce a Socratic agent for scientific reasoning, implemented through a structured system prompt that guides LLMs via classical principles of inquiry. Unlike typical prompt engineering or retrieval-based methods, our approach leverages definition, analogy, hypothesis elimination, and other Socratic techniques to generate more coherent, critical, and domain-aware responses. We evaluate the agent across diverse scientific domains and benchmark it on the abstraction and reasoning corpus challenge dataset, achieving 97.15% under a fixed prompting protocol and without fine-tuning or external tools. Expert evaluation shows improved reasoning depth, clarity, and adaptability over conventional LLM outputs, suggesting that structured prompting rooted in philosophical reasoning can improve the scientific utility of language models.

LLM reasoning↗
DOE OSTI2025

Probing the Kitaev honeycomb model on a neutral-atom quantum computer

Quantum simulations of many-body systems are among the most promising applications of quantum computers. In particular, models based on strongly correlated fermions are central to our understanding of quantum chemistry and materials problems, and can lead to exotic, topological phases of matter. However, owing to the non-local nature of fermions, such models are challenging to simulate with qubit devices. Here we realize a digital quantum simulation architecture for two-dimensional fermionic systems based on reconfigurable atom arrays. We utilize a fermion-to-qubit mapping based on Kitaev’s model on a honeycomb lattice, in which fermionic statistics are encoded using long-range entangled states. We prepare these states efficiently using measurement and feedforward, realize subsequent fermionic evolution through Floquet engineering with tunable entangling gates interspersed with atom rearrangement, and improve results with built-in error detection. Leveraging this fermion description of the Kitaev spin model, we efficiently prepare topological states across its complex phase diagram and verify the non-Abelian spin-liquid phase by evaluating an odd Chern number. We further explore this two-dimensional fermion system by realizing tunable dynamics and directly probing fermion exchange statistics. Finally, we simulate strong interactions and study the dynamics of the Fermi–Hubbard model on a square lattice. These results pave the way for digital quantum simulations of complex fermionic systems for materials science, chemistry and high-energy physics.

atomic and molecular physics↗
DOE OSTI2011

Patient-Specific 3D Pretreatment and Potential 3D Online Dose Verification of Monte Carlo-Calculated IMRT Prostate Treatment Plans

Purpose: Fast and reliable comprehensive quality assurance tools are required to validate the safety and accuracy of complex intensity-modulated radiotherapy (IMRT) plans for prostate treatment. In this study, we evaluated the performance of the COMPASS system for both off-line and potential online procedures for the verification of IMRT treatment plans. Methods and Materials: COMPASS has a dedicated beam model and dose engine, it can reconstruct three-dimensional dose distributions on the patient anatomy based on measured fluences using either the MatriXX two-dimensional (2D) array (offline) or a 2D transmission detector (T2D) (online). For benchmarking the COMPASS dose calculation, various dose-volume indices were compared against Monte Carlo-calculated dose distributions for five prostate patient treatment plans. Gamma index evaluation and absolute point dose measurements were also performed in an inhomogeneous pelvis phantom using extended dose range films and ion chamber for five additional treatment plans. Results: MatriXX-based dose reconstruction showed excellent agreement with the ion chamber (<0.5%, except for one treatment plan, which showed 1.5%), film ({approx}100% pixels passing gamma criteria 3%/3 mm) and mean dose-volume indices (<2%). The T2D based dose reconstruction showed good agreement as well with ion chamber (<2%), film ({approx}99% pixels passing gamma criteria 3%/3 mm), and mean dose-volume indices (<5.5%). Conclusion: The COMPASS system qualifies for routine prostate IMRT pretreatment verification with the MatriXX detector and has the potential for on-line verification of treatment delivery using T2D.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2026

Finite-element-based simulations of electrodes for CO 2 cascade reduction reactions

The multielectron reduction of CO 2 to liquid fuels could be a path to scalable energy storage, but reaching this goal requires major advances in catalysis and systems engineering. Cascade catalysis, which couples sequential reactions without isolating intermediates, has emerged as a promising route to enhance selectivity and efficiency in CO 2 reduction (CO 2 R). In this review, we examine how finite-element-based simulations of continuum model [finite element method (FEM)] approaches are being used to analyze and guide CO 2 R cascade systems. We first outline the fundamentals of cascade catalysis and recent advances in catalytic materials (metallic, molecular, and hybrid architectures). We then focus on FEM developments at the electrode and device scales, emphasizing how these models capture transport phenomena, local microenvironments, and geometry-dependent effects. To clarify design principles, we present case studies of cascade electrodes organized in systems without and with integrated semiconductors. We further emphasize the integration of FEM with multiscale frameworks (density functional theory, molecular dynamics, kinetic Monte Carlo) and its role in bridging atomic-level insights with device-level performance. Finally, we identify current limitations and future prospects, including improved boundary conditions, coupling with operando experiments, and machine learning-accelerated model development. Together, these insights provide design principles for next-generation CO 2 R cascade systems for efficient solar fuel production.

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

Update to the Performance Assessment for the Savannah River Site Saltstone Disposal Facility - 20124

In 2019, Savannah River Remediation developed a revision to the performance assessment (PA) on behalf of the U.S. Department of Energy (DOE) Savannah River Operations Office (SR) for the near-surface disposal of low-level waste at the Savannah River Site (SRS) Saltstone Disposal Facility (SDF). Soluble waste from SRS Tank Farms undergoes salt processing to remove cesium and other high activity constituents. The low-activity decontaminated salt solution (DSS) is then immobilized by mixing it into a cementitious waste form known as saltstone. After mixing, the saltstone is poured into leak-tight concrete vaults, known as saltstone disposal units (SDUs), where the waste form cures. By the time of facility closure, the SDF is expected to consist of 15 SDUs with a combined capacity of 1.06 E+09 L (280 Mgal) of cured saltstone. The facility operates under a Disposal Authorization Statement (DAS) from DOE and a permit from the South Carolina Department of Health and Environmental Control (SCDHEC). Since the start of operations in 1990, the SDF has received almost 6.7 E+07 L (18 Mgal) of DSS, resulting in the safe disposal of 2.7 E+16 Bq (7.3 E+05 Ci) of activity. Due to the radioactive decay of short-lived contaminants, the total remaining activity in the disposed waste is estimated to be approximately 1.4 E+16 Bq (3.9 E+05 Ci), as of September 2018. The Disposal Authorization Statement requires a demonstration that the system of engineered and natural features of the disposal facility will limit releases from the facility and be protective of human health and the environment for at least the next 1,000 years. The long-term performance of the facility was evaluated under the requirements of the DoE's Radioactive Waste Management Manual (US DOE Manual 435.1-1). Simulations were performed to demonstrate that the disposal facility would meet performance objectives specified in the manual. The evaluation was based on numerical models that simulate the releases of contaminants from the saltstone waste form. Contaminants were transported through groundwater and air pathways to points of assessment to evaluate compliance. In addition, the potential consequences of an inadvertent human intrusion (IHI) were also evaluated. In accordance with the guidance and recommendations in US DoE's technical standard for DAS, deterministic and probabilistic analyses were performed to demonstrate the SDF system, which includes the engineered cover system, the SDUs, the waste form, and the natural features of the site, provides a reasonable expectation that saltstone disposal will meet performance objectives. Post-closure doses to future members of the public (MOP) were evaluated; dose estimates from the air pathway were well below the 1.0 E-04 Sv/yr (10 mrem/yr) performance objective within 1,000 years, and the calculated doses to the MOP from all exposure pathways, including the groundwater and air pathways, were well below the 2.5 E-04 Sv/yr (25 mrem/yr) performance objective within 1,000 years. Doses following an assumed intrusion event within the facility boundaries were well below the IHI performance objectives, where the acute IHI dose was below the acute IHI dose performance objective of 5.0 E-03 Sv (500 mrem) and the chronic IHI dose was below the chronic IHI dose performance objective of 1.0 E-03 Sv/yr (100 mrem/yr) within 1,000 years. (author)

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

Asymptotic analysis of the fluid flow with a pressure-dependent viscosity in a system of thin pipes

We consider the incompressible fluid with a pressure-dependent viscosity flowing through a multiple pipe system. The viscosity–pressure relation is given by the Barus law commonly used in the engineering applications. Assuming that the ratio between pipes thickness and its length is small, we propose a rigorous asymptotic approach based on the concept of the transformed pressure. As a result, we obtain new macroscopic model describing the effective behavior of the fluid in the system. In particular, the generalized version of the Kirchhoff’s law is derived giving the explicit formula for the junction pressure. The error estimate for the asymptotic approximation is also provided. Mathematical analysis presented here can be applied to a general viscosity–pressure relation satisfied by other empiric laws.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Knowledge-matching based computational framework for genome-scale metabolic model refinement

Genome-scale metabolic models (GEMs) are mathematically structured knowledge base reconstructed from annotated genome of different organisms. With the advancement of next-generation sequencing technology, many organisms have had their genomes sequenced. However, obtaining a high-quality GEM is highly time-consuming, even with the introduction of several genome-scale reconstruction tools that offer automated draft network generation and gap filling. It has been recognized that the iterative process of manual curation and refinement is the limiting step of GEM development, and how to expedite the GEM refinement is still an open question. As cellular metabolism is a complex system with very high degree of freedom and redundancy, the principles and techniques developed in process systems engineering can be adapted to expedite GEM refinement. In this paper we present a knowledge-matching based computation framework for GEM refinement, and demonstrate the effectiveness of the proposed solution using the refinement of a GEM for Clostridium tyrobutyricum.

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

Residential Real-time Price Response Simulation

The electric industry is gaining experience with innovative price responsive demand pilots and limited roll-outs to customers. One of these pilots is investigating real-time pricing signals to engage end-use systems and local distributed generation and storage in a distributed optimization process. Attractive aspects about the approach include strong scalability characteristics, simplified interfaces between automation devices, and the adaptability to integrate a wide variety of devices and systems. Experience in this nascent field is revealing a rich array of for engineering decisions and the application of complexity theory. To test the decisions, computer simulations are used to reveal insights about design, demand elasticity, and the limits of response (including consumer fatigue). Agent-based approaches lend themselves well in the simulation to modeling the participation and interaction of each piece of equipment on a distribution feeder. This paper discusses rate design and simulation experiences at the distribution feeder level where consumers and their HVAC systems and water heaters on a feeder receive real-time pricing signals.

99 GENERAL AND MISCELLANEOUS↗