Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mathematical modelling”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Ion exchange behavior of astatine and bismuth

Astatine and bismuth sorption on several ion exchange resins from nitric acid media has been studied. This work covers commercially available resins Dowex 50 × 4, MP thiol, Dowex 1 × 8, and TEVA. One of the main advantages of using ion exchangers for astatine separation and purification is the absence of any organic media in the effluent in comparison with extraction chromatography resins. The behavior of the above-mentioned metals was investigated in up to 4 M HNO 3 solutions. The determined distribution coefficients are greater than 1 for all the resins studied, reaching 10 4 for MP thiol and TEVA resins. The interaction of nitric acid with an anion exchanger in Cl-form results in an exchange reaction, where AtO + chloro complexes can form. To understand the astatine behavior under these conditions, stability constants of AtOCl and AtOCl 2 – complexes have been reevaluated and discrepancy in the literature values has been eliminated. For each resin a thermodynamic model has been developed to suggest a possible mechanism of astatine sorption. Literature and new experimental data on bismuth sorption by studying resins in nitric acid media have been reviewed and a mathematical model to describe its behavior has been suggested. Furthermore, a ratio of corresponding fit functions of At and Bi assigned to the same resin and acidity has been used to estimate the separation factors of these elements.

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

Computational Physics Overview [Slides]

Computational physics is an important part of the overall investment in National Security Science at LANL. Computational physics is the study and implementation of numerical analysis to solve problems in physics for which a quantitative theory already exists. Historically, computational physics was the first application of modern computers in science. There are three key elements to computational physics: mathematical models of physical phenomena and conservation equations, computer codes that implement these models, and computer platforms that execute the code instructions and manipulate the data.

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

Heat pump water heater enhanced with phase change materials thermal energy storage: Modeling study

A promising solution to improve the first hour rating (FHR) of a heat pump water heater (HPWH) involves employing a secondary tank which contains phase change material (PCM) capsules. To better understand the influence of PCM thermal storage on the HPWH operational performance, a dynamic model was developed to simulate and analyze the behavior of a newly developed HPWH technology that incorporates PCM storage into a standard HPWH to optimize key parameters such as the uniform energy factor and FHR. Mathematical models of several key components of the proposed HPWH-PCM integrated thermal energy storage (TES) system, e.g., water heater tank, PCM TES tank, evaporator, compressor, and expansion valve, have been elaborated. Also, a model-based control co-simulation platform was developed to integrate a embed PCM storage HPWH dynamic model with a control model for better supporting control design, analysis, verification, and validation. The model accuracy has been validated through comparing simulation results with lab test results, with a mean average percentage error of <5.5% for most of the selected performance variables. In addition, using the developed co-simulation platform, the demand response control strategy was studied to evaluate the load flexibility of the combined HPWH-PCM storage system by shifting system power usage to outside of the 3.5 h peak load period.

25 ENERGY STORAGE↗
DOE OSTI2021

Targeting a cell surface vitamin D receptor on tumor-associated macrophages in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is an aggressive tumor with limited treatment options and poor prognosis. We applied the in vivo phage display technology to isolate peptides homing to the immunosuppressive cellular microenvironment of TNBC as a strategy for non-malignant target discovery. We identified a cyclic peptide (CSSTRESAC) that specifically binds to a vitamin D receptor, protein disulfide-isomerase A3 (PDIA3) expressed on the cell surface of tumor-associated macrophages (TAM), and targets breast cancer in syngeneic TNBC, non-TNBC xenograft, and transgenic mouse models. Systemic administration of CSSTRESAC to TNBC-bearing mice shifted the cytokine profile toward an antitumor immune response and delayed tumor growth. Moreover, CSSTRESAC enabled ligand-directed theranostic delivery to tumors and a mathematical model confirmed our experimental findings. Finally, in silico analysis showed PDIA3-expressing TAM in TNBC patients. This work uncovers a functional interplay between a cell surface vitamin D receptor in TAM and antitumor immune response that could be therapeutically exploited.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Using PyBioNetFit to leverage qualitative and quantitative data in biological model parameterization and uncertainty quantification

Data generated in studies of cellular regulatory systems are often qualitative. For example, measurements of signaling readouts in the presence and absence of mutations may reveal a rank ordering of responses across conditions but not the precise extents of mutation-induced differences. Qualitative data are often ignored by mathematical modelers or are considered in an ad hoc manner, as in the study of Kocieniewski and Lipniacki (2013) [Phys Biol 10: 035006], which was focused on the roles of MEK isoforms in ERK activation. In this earlier study, model parameter values were tuned manually to obtain consistency with a combination of qualitative and quantitative data. This approach is not reproducible, nor does it provide insights into parametric or prediction uncertainties. Here, starting from the same data and the same ordinary differential equation (ODE) model structure, we generate formalized statements of qualitative observations, making these observations more reusable, and we improve the model parameterization procedure by applying a systematic and automated approach enabled by the software package PyBioNetFit. We also demonstrate uncertainty quantification (UQ), which was absent in the original study. Our results show that PyBioNetFit enables qualitative data to be leveraged, together with quantitative data, in parameterization of systems biology models and facilitates UQ. These capabilities are important for reliable estimation of model parameters and model analyses in studies of cellular regulatory systems and reproducibility.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Modeling and control of cascaded bridgeless multilevel rectifier under unbalanced load conditions

The goal of this project is to model and control a novel unidirectional cascaded multilevel bridgeless rectifier as an active front end in medium and high voltage applications. This topology has many advantages over a conventional cascaded H-bridge rectifier, such as lower implementation cost, higher reliability, and greater flexibility with similar power quality. The complete design process of the proposed converter is developed step by step in order to meet all the desired objectives. The steady-state mathematical model is used to develop a method for the voltage balancing of dc cells. Power factor analysis is discussed to mathematically derive requirements for the number of partially controlled and fully controlled H-bridges in the proposed H-bridge converter. Power loss, efficiency, and cost comparison studies between the traditional cascaded H-Bridge converter and the proposed bridgeless converter demonstrate the advantages. After exploring various well-established control methods, a novel control strategy is proposed to achieve dc voltage balancing, fast and robust grid synchronization, power factor correction, and elimination of zero crossing current distortion under both balanced and unbalanced load conditions. The converter can also be used for reactive power compensation in a grid tied power system if a sufficient number of fully controlled H-bridge modules are included. Processor-In-the-Loop (PIL) simulation has been the utilized to validate the performance of discrete control structure. Simulation and experimental results validate the models and control method.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Response to ‘Stochastic and deterministic interpretation of pool models’

Here, we concur with Azizi-Rad et al. (2021) that it is vital to critically evaluate and compare different soil carbon models, and we welcome the opportunity to further describe the unique contribution of the PROMISE model (Waring et al., 2020) to this literature. The PROMISE framework does share many features with established biogeochemical models, as our original manuscript highlighted in Table 1, and our work builds upon model innovations developed by many different groups, including that of Azizi-Rad and colleagues. Yet, the PROMISE framework is distinctive due to where it places mechanistic emphasis, and how these mechanisms are formalized in the mathematical model structure.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A PDE Model of Breast Tumor Progression in MMTV-PyMT Mice

The evolution of breast tumors greatly depends on the interaction network among different cell types, including immune cells and cancer cells in the tumor. This study takes advantage of newly collected rich spatio-temporal mouse data to develop a data-driven mathematical model of breast tumors that considers cells’ location and key interactions in the tumor. The results show that cancer cells have a minor presence in the area with the most overall immune cells, and the number of activated immune cells in the tumor is depleted over time when there is no influx of immune cells. Interestingly, in the case of the influx of immune cells, the highest concentrations of both T cells and cancer cells are in the boundary of the tumor, as we use the Robin boundary condition to model the influx of immune cells. In other words, the influx of immune cells causes a dominant outward advection for cancer cells. We also investigate the effect of cells’ diffusion and immune cells’ influx rates in the dynamics of cells in the tumor micro-environment. Sensitivity analyses indicate that cancer cells and adipocytes’ diffusion rates are the most sensitive parameters, followed by influx and diffusion rates of cytotoxic T cells, implying that targeting them is a possible treatment strategy for breast cancer.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Maximum bound principle preserving integrating factor Runge–Kutta methods for semilinear parabolic equations

A large class of semilinear parabolic equations satisfy the maximum bound principle (MBP) in the sense that the time-dependent solution preserves for any time a uniform pointwise bound imposed by its initial and boundary conditions. Here, the MBP plays a crucial role in understanding the physical meaning and the well-posedness of the mathematical model. Investigation on numerical algorithms with preservation of the MBP has attracted increasingly attentions in recent years, especially for the temporal discretizations, since the violation of MBP may lead to nonphysical solutions or even blow-ups of the algorithms. In this paper, we study high-order MBP-preserving time integration schemes by means of the integrating factor Runge-Kutta (IFRK) method. Beginning with the space-discrete system of semilinear parabolic equations, we present the IFRK method in general form and derive the sufficient conditions for the method to preserve the MBP. In particular, we show that the classic four-stage, fourth-order IFRK scheme is MBP preserving for some typical semilinear systems although not strong stability preserving, which can be instantly applied to the Allen-Cahn type of equations. To our best knowledge, this is the first time to present a fourth-order linear numerical method preserving the MBP. In addition, convergence of these numerical schemes is proved theoretically and verified numerically, as well as their efficiency by simulations of 2D and 3D long-time evolutional behaviors. Numerical experiments are also carried out for a model which is not a typical gradient flow as the Allen-Cahn type of equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A comprehensive modeling method for proton exchange membrane electrolyzer development

Hydrogen attracts significant interests as an effective energy carrier that can be derived from renewable sources. Hydrogen production using a proton-exchange membrane (PEM) electrolyzer can efficiently convert renewable power via water splitting in wide scales—from large, centralized generation to on-site production. Mathematical models with multiple scales and fidelities facilitate the continuing improvements of PEM electrolyzer development to improve performance, cost, and reliability. The model scopes and methods are presented in this paper, which also introduces a comprehensive PEM electrolysis modeling tool based on computational fluid dynamics (CFD) software, ANSYS/Fluent. The modeling tool incorporates electrochemical model of a PEM electrolysis cell to simulate the performance of coupled thermal-fluid, species transport, and electrochemical processes in a product-scale cell or stack by leveraging the powerful meshing generation and CFD solver of ANSYS/Fluent. The thermal-fluid modeling includes liquid water/gas two-phase flow and simulates a PEM electrolysis cell by using Fluent user-defined functions as add-on modules accounting for PEM-specific species transport and electrochemical processes. The modeling outcomes expediate PEM electrolyzer scaling up from basic material development and laboratory testing.

08 HYDROGEN↗
DOE OSTI2022

Design and application of a kinetic model of lipid metabolism in Saccharomyces cerevisiae

Lipid biosynthesis plays a vital role in living cells and has been increasingly engineered to overproduce various lipid-based chemicals. However, owing to the tightly constrained and interconnected nature of lipid biosynthesis, both understanding and engineering of lipid metabolism remain challenging, even with the help of mathematical models. Here we report the development of a kinetic metabolic model of lipid metabolism in Saccharomyces cerevisiae that integrates fatty acid biosynthesis, glycerophospholipid metabolism, sphingolipid metabolism, storage lipids, lumped sterol synthesis, and the synthesis and transport of relevant target-chemicals, such as fatty acids and fatty alcohols. The model was trained on lipidomic data of a reference S. cerevisiae strain, single knockout mutants, and lipid overproduction strains reported in literature. The model was used to design mutants for fatty alcohol overproduction and the lipidomic analysis of the resultant mutant strains coupled with model-guided hypothesis led to discovery of a futile cycle in the triacylglycerol biosynthesis pathway. In addition, the model was used to explain successful and unsuccessful mutant designs in metabolic engineering literature. Thus, this kinetic model of lipid metabolism can not only enable the discovery of new phenomenon in lipid metabolism but also the engineering of mutant strains for overproduction of lipids.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Vibration Studies for the CLARA Interferometer

The purpose of the CLARA experiment is to study the nature of undulator radiation emitted by single electrons circulating in the IOTA storage ring. The classical and quantum properties of the radiation are investigated by measuring its coherence length, intensity fluctuations and time correlations. The experiment took data in IOTA Run 4 (2022-2023). The key component of the apparatus is a Mach-Zehnder interferometer (MZI), in which the optical length of one of the arms can be precisely controlled. For some measurements, the fine regulation of the arm length must be smaller than the radiation wavelength. For this reason, the apparatus is particularly sensitive to mechanical vibrations. In this note, we model and measure the effect of vibrations on the performance of the MZI under various conditions. Several improvements of the setup were implemented to minimize systematic distortions of the observed interference patterns. In this report, we present a mathematical model of the effect of mechanical vibrations on observed detector and coincidence rates. We also describe the measurements that were made to estimate the magnitude and spectra of rate fluctuations and their sources in the CLARA MZI under various conditions. Finally, we estimate the magnitude of arm length fluctuations and we deduce the sensitivity to coincidence-rate variations in our apparatus.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Contributions of biofilm-induced flow heterogeneities to solute retention and anomalous transport features in porous media

We report that microbial biofilms are ubiquitous within porous media and the dynamics of their growth influence surface and subsurface flow patterns which impacts the physical properties of porous media and large-scale transport of solutes. Here, a two-dimensional pore-scale numerical model was used to evaluate the impact of biofilm-induced flow heterogeneities on conservative transport. Our study integrates experimental biofilm images of Paenibacillus 300A strain in a microfluidic device packed with cylindrical grains in a hexagonal distribution, with mathematical modeling. Biofilm is represented as a synthetic porous structure with locally varying physical properties that honors the impact of biofilm on the porous medium. We find that biofilm plays a major role in shaping the observed conservative transport dynamics by enhancing anomalous characteristics. More specifically, when biofilm is present, the pore structure in our geometry becomes more spatially correlated. We observe intermittent behavior in the Lagrangian velocities that switches between fast transport periods and long trapping events. Our results suggest that intermittency enhances solute spreading in breakthrough curves which exhibit extreme anomalous slope at intermediate times and very marked late solute arrival due to solute retention. The efficiency of solute retention by the biofilm is controlled by a transport regime which can extend the tailing in the concentration breakthrough curves. These results indicate that solute retention by the biofilm exerts a strong control on conservative solute transport at pore-scale, a role that to date has not received enough attention.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Modelling and simulation of brinicle formation

Below the Arctic sea ice, under the right conditions, a flux of icy brine flows down into the sea. The icy brine has a much lower fusion point and is denser than normal seawater. As a result, it sinks while freezing everything around it, forming an ice channel called a brinicle (also known as ice stalactite). In this paper, we develop a mathematical model for this phenomenon, assuming cylindrical symmetry. The fluid is considered to be viscous and quasi-stationary. The heat and salt transport are weakly coupled to the fluid motion and are modelled with the corresponding conservation equations, accounting for diffusive and convective effects. Finite-element discretization is employed to solve the coupled system of partial differential equations. We find that the model can capture the general behaviour of the physical system and generate brinicle-like structures while also recovering dendrite composition, which is a physically expected feature aligned with previous experimental results. This represents, to our knowledge, the first complete model proposed that captures the global structure of the physical phenomenon even though it has some discrepancies, such as brine accumulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Models of SIV rebound after treatment interruption that involve multiple reactivation events

In order to assess the efficacy of novel HIV-1 treatments leading to a functional cure, the time to viral rebound is frequently used as a surrogate endpoint. The longer the time to viral rebound, the more efficacious the therapy. In support of such an approach, mathematical models serve as a connection between the size of the latent reservoir and the time to HIV-1 rebound after treatment interruption. The simplest of such models assumes that a single successful latent cell reactivation event leads to observable viremia after a period of exponential viral growth. Here we consider a generalization developed by Pinkevych et al. and Hill et al. of this simple model in which multiple reactivation events can occur, each contributing to the exponential growth of the viral load. We formalize and improve the previous derivation of the dynamics predicted by this model, and use the model to estimate relevant biological parameters from SIV rebound data. We confirm a previously described effect of very early antiretroviral therapy (ART) initiation on the rate of recrudescence and the viral load growth rate after treatment interruption. We find that every day ART initiation is delayed results in a 39% increase in the recrudescence rate (95% credible interval: [18%, 62%]), and a 11% decrease of the viral growth rate (95% credible interval: [4%, 20%]). We show that when viral rebound occurs early relative to the viral load doubling time, a model with multiple successful reactivation events fits the data better than a model with only a single successful reactivation event.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Mechanism of astatine and bismuth sorption on extraction chromatography resins from nitric acid media

Systematic study of astatine (At) sorption on extraction chromatography resins from nitric acid media has been performed at Texas A&M University. This work covers commercially available resins such as LN, TK400, SR, TRU, and UTEVA. Acquiring an understanding of At behavior in up to 4 M HNO 3 solutions is a key step to develop a rapid and reliable system for At separation and purification. The determined distribution coefficients are greater than 20 for all the resins studied, reaching 400 for TRU and UTEVA resins. The sorption uptake decreases in the order: TRU ≳ UTEVA >TK400 ≳ SR >LN. For each resin a thermodynamic model has been developed to suggest a possible mechanism of At sorption and corresponding sorption constants have been estimated. Literature data along with new results on bismuth (Bi) sorption by studied resins in nitric acid media have been reviewed and a mathematical model to describe its behavior has been suggested. A ratio of corresponding fit functions of At and Bi assigned to the same resin and acidity has been used to estimate separation factors of these elements. Furthermore, only the TRU resin is not suitable for efficient At isolation from Bi, while the other four resins demonstrate promising results.

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

Engineering tertiary chirality in helical biopolymers

Tertiary chirality describes the handedness of supramolecular assemblies and relies not only on the primary and secondary structures of the building blocks but also on topological driving forces that have been sparsely characterized. Helical biopolymers, especially DNA, have been extensively investigated as they possess intrinsic chirality that determines the optical, mechanical, and physical properties of the ensuing material. Here, we employ the DNA tensegrity triangle as a model system to locate the tipping points in chirality inversion at the tertiary level by X-ray diffraction. We engineer tensegrity triangle crystals with incremental rotational steps between immobile junctions from 3 to 28 base pairs (bp). We construct a mathematical model that accurately predicts and explains the molecular configurations in both this work and previous studies. Our design framework is extendable to other supramolecular assemblies of helical biopolymers and can be used in the design of chiral nanomaterials, optically active molecules, and mesoporous frameworks, all of which are of interest to physical, biological, and chemical nanoscience.

Science & Technology - Other Topics↗
DOE OSTI2024

Deciphering the morphology of transition metal carbonate cathode precursors

The performance and life of Li-ion battery cathode materials are determined by both the composition (crystal structure and transition metal ratio) and the morphology (particle size, size distribution, and surface area). Careful control of these two aspects is the key to long lasting, high-energy batteries that can undergo fast charge. Developing such cathodes requires manipulation of the synthesis conditions, namely the coprecipitation process to develop the precursor and a calcination step to lithiate and convert it to a transition metal oxide. In this paper, we utilize a combination of controlled synthesis, microscopic and spectroscopic characterization, and multi-scale mathematical modeling to shed light on the synthesis of cathode precursors. The complex interplay between the various chemical reactions in the co-precipitation process is studied to provide experimentalists with guidance on achieving composition control during synthesis. Further, the formation of a variety of morphologies of the primary particles and the driving force for agglomeration is mathematically described, for the first time, based on an energy minimization approach. Results suggest that the presence of Ni and/or Co significantly lowers the reaction rate constant compared to Mn, resulting in agglomerated growth in the former and single crystal growth in the latter. Modeling studies are used to provide a phase map describing the synthesis conditions needed to control the secondary particle size and corresponding size distribution. Finally, this paper represents an important step in developing a computationally guided approach to the synthesis of battery cathode materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗