Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mathematical model”

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

Reduced order models for Lagrangian hydrodynamics

It is reported, as a mathematical model of high-speed flow and shock wave propagation in a complex multimaterial setting, Lagrangian hydrodynamics is characterized by moving meshes, advection-dominated solutions, and moving shock fronts with sharp gradients. These challenges hinder the existing projection-based model reduction schemes from being practical. We develop several variations of projection-based reduced order model techniques for Lagrangian hydrodynamics by introducing three different reduced bases for position, velocity, and energy fields. A time-windowing approach is also developed to address the challenge imposed by the advection-dominated solutions. Lagrangian hydrodynamics is formulated as a nonlinear problem, which requires a proper hyper-reduction technique. Therefore, we apply the over-sampling DEIM and SNS approaches to reduce the complexity due to the nonlinear terms. Finally, we also present both a posteriori and a priori error bounds associated with our reduced order model. We compare the performance of the spatial and time-windowing reduced order modeling approaches in terms of accuracy and speed-up with respect to the corresponding full order model for several numerical examples, namely Sedov blast, Gresho vortices, Taylor-Green vortices, and triple-point problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Development and Evaluation of Occupancy-Aware Model Predictive Control for Residential Building Energy Efficiency and Occupant Comfort

The residential sector accounts for 25% of global primary energy consumption. Two methods have previously been proposed to reduce residential energy use associated with the provision of occupant thermal comfort: 1. Occupancy-based HVAC control, operating systems only during confirmed occupancy, and 2. model predictive control (MPC), harnessing a mathematical model and forecasts to find optimal operating strategies. Previous studies estimate the average energy savings of the two methods individually in the range of 21% and 16%, respectively. The research presented herein was carried out to evaluate the energy savings potential in residential buildings by combining both approaches across different climates, house vintages, and occupancy patterns. Occupancy and eight different physical modalities (e.g. CO2 and VOC) data were collected from five homes for time periods of 4–9 weeks. Collected data sets were used to train occupancy prediction models suggested by an extensive literature survey of occupancy model types. The trained prediction models were combined with MPC and detailed EnergyPlus building simulation models to evaluate residential building performance in terms of annual energy savings and thermal comfort, along with discomfort exceedance metrics. Multiple home types and regions were analyzed to understand regional and climate-dependent potential. Based on actual field data, the occupancy models had a prediction inaccuracy between 8% and 35% across the investigated homes. Average occupancy for the collected data ranged from 56% to 86%, a typical range reported in the literature. Building simulations were conducted for three control scenarios: conventional thermostatic control, occupancy-based, and occupancy-based MPC. The results indicate that all advanced strategies improve upon the conventional control, with some scenarios cutting energy use in half with only occasional incurrence of discomfort. The findings indicate that occupancy-aware model predictive residential building control has the potential to drastically reduce energy use and associated emissions while maintaining occupant comfort for both new and existing buildings.

ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATION↗
DOE OSTI2023

Catalyst: Fast and flexible modeling of reaction networks

We introduce Catalyst.jl, a flexible and feature-filled Julia library for modeling and high-performance simulation of chemical reaction networks (CRNs). Catalyst supports simulating stochastic chemical kinetics (jump process), chemical Langevin equation (stochastic differential equation), and reaction rate equation (ordinary differential equation) representations for CRNs. Through comprehensive benchmarks, we demonstrate that Catalyst simulation runtimes are often one to two orders of magnitude faster than other popular tools. More broadly, Catalyst acts as both a domain-specific language and an intermediate representation for symbolically encoding CRN models as Julia-native objects. This enables a pipeline of symbolically specifying, analyzing, and modifying CRNs; converting Catalyst models to symbolic representations of concrete mathematical models; and generating compiled code for numerical solvers. Leveraging ModelingToolkit.jl and Symbolics.jl, Catalyst models can be analyzed, simplified, and compiled into optimized representations for use in numerical solvers. Finally, we demonstrate Catalyst’s broad extensibility and composability by highlighting how it can compose with a variety of Julia libraries, and how existing open-source biological modeling projects have extended its intermediate representation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Rich Dynamics of a General Producer–Grazer Interaction Model under Shared Multiple Resource Limitations

Organism growth is often determined by multiple resources interdependently. However, growth models based on the Droop cell quota framework have historically been built using threshold formulations, which means they intrinsically involve single-resource limitations. In addition, it is a daunting task to study the global dynamics of these models mathematically, since they employ minimum functions that are non-smooth (not differentiable). To provide an approach to encompass interactions of multiple resources, we propose a multiple-resource limitation growth function based on the Droop cell quota concept and incorporate it into an existing producer–grazer model. The formulation of the producer’s growth rate is based on cell growth process time-tracking, while the grazer’s growth rate is constructed based on optimal limiting nutrient allocation in cell transcription and translation phases. We show that the proposed model captures a wide range of experimental observations, such as the paradox of enrichment, the paradox of energy enrichment, and the paradox of nutrient enrichment. Together, our proposed formulation and the existing threshold formulation provide bounds on the expected growth of an organism. Moreover, the proposed model is mathematically more tractable, since it does not use the minimum functions as in other stoichiometric models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Predictive Models and High-Performance Computing as Tools to Accelerate the Scaling-Up of New Bio-Based Fuels Workshop: Summary Report

This report summarizes the results of a virtual workshop sponsored by the Bioenergy Technologies Office held on June 9–11, 2020. The workshop discussed best practices for utilizing mathematical modeling tools across multiple scales to reduce technology uncertainty and accelerate scaling-up of biorefinery/chemical production equipment and optimize operations.

09 BIOMASS FUELS↗
DOE OSTI2020

Evaluation of temperature changes in a cell culture plate well with absorbing walls containing a 3D model of tumour spheroid under irradiation in the optical range

The temperature conditions for photodynamic therapy of the culture of tumour cells organised in 3D spheroids are studied. The spheroids are placed in the cell culture medium in the well of a special cell culture plate with absorbing walls. A mathematical model is presented that describes the temperature distribution in a well with a cell culture plate heated by optical radiation. It is shown that the temperature in the well volume is not uniformly distributed, and the temperature on the well surface can differ significantly from the temperature in its volume. (laser medicine)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Modeling hepatitis C virus kinetics during liver transplantation highlights the role of the liver in virus clearance

While the liver, specifically hepatocytes, are widely accepted as the main source of hepatitis C virus (HCV) production, the role of the liver/hepatocytes in clearance of circulating HCV remains unknown. Frequent HCV kinetic data were recorded and mathematically modeled from five liver transplant patients throughout the anhepatic (absence of liver) phase and for 4 hr post-reperfusion. During the anhepatic phase, HCV remained at pre-anhepatic levels (n = 3) or declined (n = 2) with t 1/2 ~1 hr. Immediately post-reperfusion, virus declined in a biphasic manner in four patients consisting of a rapid decline (t 1/2 = 5 min) followed by a slower decline (t 1/2 = 67 min). Consistent with the majority of patients in the anhepatic phase, when we monitored HCV clearance at 37°C from culture medium in the absence/presence of chronically infected hepatoma cells that were inhibited from secreting HCV, the HCV t 1/2 in cell culture was longer in the absence of chronically HCV-infected cells. The results suggest that the liver plays a major role in the clearance of circulating HCV and that hepatocytes may be involved.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Ergodic Lagrangian dynamics in a superhero universe

We present a fictional scenario that, while undeniably whimsical, provides the foundation for a unique exercise in extended problem solving, physics analysis, and quantitative model development. Starting with the foundational premise of the Wild Cards shared-world superhero universe, we demonstrate how a variety of concepts appropriate to the advanced undergraduate level—ergodicity, functional analysis, Lagrangian mechanics, and the ever-important simplifying approximation—can be combined into a rich, coherent mathematical model. The goal of this case study is to develop a useful pedagogical exercise in exploring an open-ended research question that presents, at first glance, no clear path forward. Being both eclectic and lengthy, this exercise offers a unique way for students to apply their core physics and mathematics education. It is perhaps best used within a senior honors seminar or within a brief (e.g., January term) elective class.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Modeling suggests SARS-CoV-2 rebound after nirmatrelvir-ritonavir treatment is driven by target cell preservation coupled with incomplete viral clearance

In a subset of SARS-CoV-2-infected individuals treated with the antiviral nirmatrelvir-ritonavir, the virus rebounds following treatment. The mechanisms driving this rebound are not well understood. We used a mathematical model to describe the longitudinal viral load dynamics of 51 individuals treated with nirmatrelvir-ritonavir, 20 of whom rebounded. Target cell preservation, either by a robust innate immune response or initiation of N-R near the time of symptom onset, coupled with incomplete viral clearance, appears to be the main factor leading to viral rebound. Moreover, the occurrence of viral rebound is likely influenced by the time of treatment initiation relative to the progression of the infection, with earlier treatments leading to a higher chance of rebound. A comparison with an untreated cohort suggests that early treatments with nirmatrelvir-ritonavir may be associated with a delay in the onset of an adaptive immune response. Nevertheless, our model demonstrates that extending the course of nirmatrelvir-ritonavir treatment to a 10-day regimen may greatly diminish the chance of rebound in people with mild-to-moderate COVID-19 and who are at high risk of progression to severe disease. Altogether, our results suggest that in some individuals, a standard 5-day course of nirmatrelvir-ritonavir starting around the time of symptom onset may not completely eliminate the virus. Thus, after treatment ends, the virus can rebound if an effective adaptive immune response has not fully developed. These findings on the role of target cell preservation and incomplete viral clearance also offer a possible explanation for viral rebounds following other antiviral treatments for SARS-CoV-2.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Parsimonious models of in-host viral dynamics and immune response

Mathematical models of in-host viral dynamics and immune response are a vital tool for patient-specific estimation of the initial viral load, prediction of the course of an infection, etc. The COVID-19 pandemics has given impetus to the development of models with an ever-increasing degree of complexity. We show that one of the most popular models---the Target Cell Limited model---fails the identifiability test, i.e., its parameters cannot be uniquely inferred from readily available data such as viral load measurements. Here, we present a model that is both identifiable and parsimonious according to information criteria. Our model's predictions match both reported observations of COVID-19 patients and predictions of its more complex counterparts.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Posiva's Onkalo{sup R} FST (Full-Scale In Situ System Test) - 20051

The FST - Full-scale In Situ System Test - is a unique full-scale demonstration of disposal system in crystalline host rock. The FST was constructed in ONKALO{sup R} demonstration area at the disposal depth of 420 metres, in Posiva's actual final disposal facility at Olkiluoto, Finland. The demonstration has included design, manufacture, installation and comprehensive instrumentation for monitoring the evolution of the disposal system and its components. The demonstration tunnel and experimental deposition holes were excavated in 2010 in ONKALO to address that they can be excavated to meet the requirements set for them. Two experimental deposition holes, bored in the floor of the tunnel, host each a canister surrounded by bentonite clay buffer rings/segments in the holes. Inside the canisters are heating elements to simulate decay heat emitted by spent fuel. The tunnel part used for the test is about 50 m long and it was filled with bentonite clay backfill after the installation of the canisters and buffer components. The tunnel is closed with a steel-reinforced concrete plug. The test is equipped with 500 sensors to monitor early evolution of the system. To forecast the evolution of the system, a mathematical model has been created. (authors)

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

Modeling HIV-1 Within-Host Dynamics After Passive Infusion of the Broadly Neutralizing Antibody VRC01

VRC01 is a broadly neutralizing antibody that targets the CD4 binding site of HIV-1 gp120. Passive administration of VRC01 in humans has assessed the safety and the effect on plasma viremia of this monoclonal antibody (mAb) in a phase 1 clinical trial. After VRC01 infusion, the plasma viral load in most of the participants was reduced but had particular dynamics not observed during antiretroviral therapy. In this paper, we introduce different mathematical models to explain the observed dynamics and fit them to the plasma viral load data. Based on the fitting results we argue that a model containing reversible Ab binding to virions and clearance of virus-VRC01 complexes by a two-step process that includes (1) saturable capture followed by (2) internalization/degradation by phagocytes, best explains the data. This model predicts that VRC01 may enhance the clearance of Ab-virus complexes, explaining the initial viral decay observed immediately after antibody infusion in some participants. Because Ab-virus complexes are assumed to be unable to infect cells, i.e., contain neutralized virus, the model predicts a longer-term viral decay consistent with that observed in the VRC01 treated participants. By assuming a homogeneous viral population sensitive to VRC01, the model provides good fits to all of the participant data. However, the fits are improved by assuming that there were two populations of virus, one more susceptible to antibody-mediated neutralization than the other.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

A multiscale model of the action of a capsid assembly modulator for the treatment of chronic hepatitis B

Chronic hepatitis B virus (HBV) infection is strongly associated with increased risk of liver cancer and cirrhosis. While existing treatments effectively inhibit the HBV life cycle, viral rebound frequently occurs following treatment interruption. Consequently, functional cure rates of chronic HBV infection remain low and there is increased interest in a novel treatment modality, capsid assembly modulators (CAMs). Here, we develop a multiscale mathematical model of CAM treatment in chronic HBV infection. By fitting the model to participant data from a phase I trial of the first-generation CAM vebicorvir, we estimate the drug’s dose-dependent effectiveness and identify the physiological mechanisms that drive the observed biphasic decline in HBV DNA and RNA, and mechanistic differences between HBeAg-positive and negative infection. Finally, we demonstrate analytically and numerically that the relative change of HBV RNA more accurately reflects the antiviral effectiveness of a CAM than the relative change in HBV DNA.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Modeling resistance to the broadly neutralizing antibody PGT121 in people living with HIV-1

PGT121 is a broadly neutralizing antibody in clinical development for the treatment and prevention of HIV-1 infection via passive administration. PGT121 targets the HIV-1 V3-glycan and demonstrated potent antiviral activity in a phase I clinical trial. Resistance to PGT121 monotherapy rapidly occurred in the majority of participants in this trial with the sampled rebound viruses being entirely resistant to PGT121 mediated neutralization. However, two individuals experienced long-term ART-free viral suppression following antibody infusion and retained sensitivity to PGT121 upon viral rebound. Here, we develop mathematical models of the HIV-1 dynamics during this phase I clinical trial. We utilize these models to understand the dynamics leading to PGT121 resistance and to identify the mechanisms driving the observed long-term viral control. Our modeling highlights the importance of the relative fitness difference between PGT121 sensitive and resistant subpopulations prior to treatment. Specifically, by fitting our models to data, we identify the treatment-induced competitive advantage of previously existing or newly generated resistant population as a primary driver of resistance. Finally, our modeling emphasizes the high neutralization ability of PGT121 in both participants who exhibited long-term viral control.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

DC Microgrids Under Denial of Service Attacks: Feasibility and Stability Issues

Distributed control has been extensively studied for a DC microgrid to coordinate multiple grid assets where each asset can generate its own control decisions based on communications within its neighborhood. Nevertheless, such a control strategy usually has a strong dependence on active information exchange and needs to always maintain a pre-determined communication graph. This makes it vulnerable to cyber attacks, especially denial of service (DoS) attacks that block communication channels that effectively change the communication graph. In this paper, we study two potential impacts of DoS attacks: the inability to maintain steady-state and dynamical generation-demand balance, i.e. the feasibility and stability issues. We develop a mathematical model to describe the cyber-induced system steady-state and dynamical performance. The model reveals the impacts of DoS attacks and enables formal analysis. Case studies are shown to verify the proposed work.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2020

Quantitative and Predictive Genetic Parts for Plant Synthetic Biology

Plant synthetic biology aims to harness the natural abilities of plants and to turn them to new purposes. A primary goal of plant synthetic biology is to produce predictable and programmable genetic circuits from simple regulatory elements and well-characterized genetic components. The number of available DNA parts for plants is increasing, and the methods for rapid quantitative characterization are being developed, but the field of plant synthetic biology is still in its early stages. We here describe methods used to describe the quantitative properties of genetic components needed for plant synthetic biology. Once the quantitative properties and transfer function of a variety of genetic parts are known, computers can select the optimal components to assemble into functional devices, such as toggle switches and positive feedback circuits. However, while the variety of circuits and traits that can be put into plants are limitless, doing synthetic biology in plants poses unique challenges. Plants are composed of differentiated cells and tissues, each representing potentially unique regulatory or developmental contexts to introduced synthetic genetic circuits. Further, plants have evolved to be highly sensitive to environmental influences, such as light or temperature, any of which can affect the quantitative function of individual parts or whole circuits. Measuring the function of plant components within the context of a plant cell and, ideally, in a living plant, will be essential to using these components in gene circuits with predictable function. Mathematical modeling will be needed to account for the variety of contexts a genetic part will experience in different plant tissues or environments. With such understanding in hand, it may be possible to redesign plant traits to serve human and environmental needs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Zika virus dynamics: Effects of inoculum dose, the innate immune response and viral interference

Experimental Zika virus infection in non-human primates results in acute viral load dynamics that can be well-described by mathematical models. The inoculum dose that would be received in a natural infection setting is likely lower than the experimental infections and how this difference affects the viral dynamics and immune response is unclear. Here we study a dataset of experimental infection of non-human primates with a range of doses of Zika virus. We develop new models of infection incorporating both an innate immune response and viral interference with that response. We find that such a model explains the data better than models with no interaction between virus and the immune response. We also find that larger inoculum doses lead to faster dynamics of infection, but approximately the same total amount of viral production.

59 BASIC BIOLOGICAL SCIENCES↗