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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Mathematical modelling of the concave front in the adjacent high explosive detonation problem

This study presents an analysis of the transition-zone in adjacent high explosive (HE) detonation problems which uses a $D, 𝜅, \dot{D}$ relationship, where $D$ is the detonation front-normal velocity, 𝜅 is the detonation front curvature and $\dot{D}$ is the time derivative of detonation front-normal velocity. Our approach extends the traditional $(D, 𝜅)$ model to accurately predict the behaviour of both diverging and converging detonation shock fronts. Our findings affirm that a hyperbolic type of front evolution equation, enhanced with wave acceleration, provides a robust framework for modelling complex shock front dynamics in HE materials. This approach not only captures the natural effects of straightness and boundary slope jumps in the transition-zone but also bridges the gap between mathematical predictions and experimental observations, offering insights into the behaviour of both diverging and converging detonation propagations in a homogeneous HE.

acceleration↗
DOE OSTI2026

Mathematical Models and Numerical Methods for High-Fidelity Simulation of Ignition of Reactive Mixtures by Nanosecond Plasma Discharges in Realistic Configurations

We present a newly developed framework for the numerical simulation of ignition of reactive mixtures using single or repeated nanosecond discharge pulses. The framework builds upon the AMReX library, using the existing compressible solver PeleC and low-Mach solver PeleLMeX and allowing for adaptive mesh refinement, complex geometries, and execution on next-generation high-performance computing (HPC) systems. High-fidelity elementary models are adopted for weakly-ionised plasma discharges with significant energy deposition, consistent with nanosecond discharge pulses, and then implemented in the solver. The treatment of non-thermal electrons and charged species, thermodynamics of non-equilbrium species, plasma kinetics, limiting time scales, and boundary conditions for charged species are discussed and addressed for computational efficiency. The framework is demonstrated for three relevant applications: single and multi-pulse discharges in air, single pulse ignition of an ethylene/air mixture, and a three-dimensional plasma discharge in air with temperature stratification. The successful application of the framework demonstrates the feasibility of high-fidelity simulation of ignition of air/hydrocarbon mixtures in three-dimensions with multiple discharge pulses.

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

Analysis of tumor-immune functional responses in a mathematical model of neoantigen cancer vaccines

Cancer neoantigen vaccines have emerged as a promising approach to stimulating the immune system to fight cancer. We propose a simple model including key elements of cancer-immune interactions and conduct a phase plane analysis to understand the immunological mechanisms of cancer neoantigen vaccines. Analytical results are obtained for two widely used functional forms that represent the killing rate of tumor cells by immune cells: the law of mass action (LMA) and the dePillis-Radunskaya Law (LPR). Using the LMA, our results reveal that a slowly growing tumor can escape the immune surveillance and that there is a unique periodic solution. The LPR offers richer dynamics, in which tumor elimination and uncontrolled tumor growth are both present. We show that tumor elimination requires sufficient number of initial activated T cells in relationship to the malignant cells, which lends support to using the neoantigen cancer vaccine as an adjuvant therapy after the primary tumor is surgically removed or treated using radiotherapy. We also derive a sufficient condition for uncontrolled tumor growth under the assumption of the LPR. Here, the juxtaposition of analyses with these two different choices for the killing rate function highlights their importance on model behavior and biological implications, by which we hope to spur further theoretical and experimental work to understand mechanisms underlying different functional forms for the killing rate.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Evaluation of the rheological and electrical percolation of high‐density polyethylene/carbon black composites using mathematical models

Abstract In this work, conductive polymer composites (CPCs) of bio‐based polyethylene (BioPe) containing different concentrations of carbon black (CB) were developed. By using oscillatory rheology analysis, a Newtonian plateau was observed in BioPe, and all BioPe/CB composites had a behavior of a pseudo‐solid and that composites with volume fractions ranging from 0.24 to 0.56 presented higher viscosity, storage, and loss modulus. This suggests the formation of a percolated network and by using the power‐law models, it was observed that the electrical percolation threshold was higher than the rheological percolation threshold. The electrical conductivity was measured using the four‐point probe method and a sigmoid model was used to predict the CPCs' electrical conductivity percolation threshold. The results indicated that the four‐point probe method presented satisfactory results according to the calculated standard deviations and voltage–current characteristics for each round of measurements considering the same ranging as used in rheology analysis. The analytical model used showed a coefficient of determination ( R 2 ) higher than 95%, allowing the prediction of the electrical conductivity of the CPC and the percolation threshold as a function of the volumetric fraction of the CB.

da SiIva, Moacy P.↗
DOE OSTI2021

Mathematical modeling of novel porous transport layer architectures for proton exchange membrane electrolysis cells

Thin foil based porous transport layers (PTLs) that contain highly structured pore arrays have shown promise as anode PTLs in proton exchange membrane electrolysis cells. These novel PTLs, fabricated with advanced manufacturing techniques, produce thin, tunable, multifunctional layers with reduced flow and interfacial resistances and high thermal and electric conductivities. To further optimize their design, it is important to understand their fundamental impact on the transport of protons, electrons, and liquid/vapor mixtures in the electrode. In this work, we develop a two-dimensional multiphysics model to simulate the coupled electrochemistry and multiphase transport in an electrolysis cell operated with the novel PTL architecture. The results show that larger pores improve access of water to the anode catalyst layer, which is beneficial for both the oxygen evolution reaction and membrane hydration. Larger pore sizes also improve oxygen gas transport from the catalyst layer, because generated oxygen gas is forced to travel in-plane through the anode catalyst layer until it reaches a pore opening that is connected to a channel. The discussed results confirm that the proposed thin foil based PTLs are fundamentally different from conventional PTLs, such as felts or layered meshes. The model developed in this work also provides generalizable insight into fundamental PEMEC phenomena, such as the competition between liquid and gas phase transport, membrane hydration and water management, and nonuniform electrochemical reactions, which are processes relevant to all PEMEC designs.

25 ENERGY STORAGE↗
DOE OSTI2020

A mathematical model of solid-state dewetting of barium thin films on W (112)

Solid state dewetting occurs when a thin solid film is heated. The temperature of dewetting depends on the thickness of the film; dewetting can be observed in the range of 1/3 to 1/2 of the bulk melting temperature. While remaining solid, the film behaves in a manner similar to liquids dewetting and agglomerating to forming islands or droplets. One of the possible mechanisms is the conversion of a metastable thin film geometry into a more stable form. Heating the metastable film gives the film atoms higher mobility, and the films retract, dewetting the surface. This atomic motion can be restricted due to surface anisotropy. Herein we present in situ emission microscopy observations of barium thin films deposited onto W (112) by thermal evaporation. From the modeling viewpoint, the evolution of the film in this system could be divided in four stages: (i) the nucleation and growth of the thin film as a simply connected region; (ii) formation of droplets/islands/hillocks; (iii) nucleation of holes; (iv) evolution of the components of the disconnected film to their equilibrium state. We present a continuum model that is qualitatively consistent with the evolution of the film observed at the initial stage of the experiment and discuss the later stages of the evolution of surface structures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A mathematical model of asynchronous data flow in parallel computers *

Abstract We present a simplified model of data flow on processors in a high-performance computing framework involving computations necessitating inter-processor communications. From this ordinary differential model, we take its asymptotic limit, resulting in a model which treats the computer as a continuum of processors and data flow as an Eulerian fluid governed by a conservation law. We derive a Hamilton–Jacobi equation associated with this conservation law for which the existence and uniqueness of solutions can be proven. We then present the results of numerical experiments for both discrete and continuum models; these show a qualitative agreement between the two and the effect of variations in the computing environment’s processing capabilities on the progress of the modelled computation.

Barnard, Richard C.↗
DOE OSTI2020

Mathematical modeling of hepatitis C RNA replication, exosome secretion and virus release

Hepatitis C virus (HCV) causes acute hepatitis C and can lead to life-threatening complications if it becomes chronic. The HCV genome is a single plus strand of RNA. Its intracellular replication is a spatiotemporally coordinated process of RNA translation upon cell infection, RNA synthesis within a replication compartment, and virus particle production. While HCV is mainly transmitted via mature infectious virus particles, it has also been suggested that HCV-infected cells can secrete HCV RNA carrying exosomes that can infect cells in a receptor independent manner. In order to gain insight into these two routes of transmission, we developed a series of intracellular HCV replication models that include HCV RNA secretion and/or virus assembly and release. Fitting our models to in vitro data, in which cells were infected with HCV, suggests that initially most secreted HCV RNA derives from intracellular cytosolic plus-strand RNA, but subsequently secreted HCV RNA derives equally from the cytoplasm and the replication compartments. Furthermore, our model fits to the data suggest that the rate of virus assembly and release is limited by host cell resources. Including the effects of direct acting antivirals in our models, we found that in spite of decreasing intracellular HCV RNA and extracellular virus concentration, low level HCV RNA secretion may continue as long as intracellular RNA is available. This may possibly explain the presence of detectable levels of plasma HCV RNA at the end of treatment even in patients that ultimately attain a sustained virologic response.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

The effect of enrofloxacin on enteric Escherichia coli : Fitting a mathematical model to in vivo data

Antimicrobial drugs administered systemically may cause the emergence and dissemination of antimicrobial resistance among enteric bacteria. To develop logical, research-based recommendations for food animal veterinarians, we must understand how to maximize antimicrobial drug efficacy while minimizing risk of antimicrobial resistance. Our objective is to evaluate the effect of two approved dosing regimens of enrofloxacin (a single high dose or three low doses) on Escherichia coli in cattle. We look specifically at bacteria above and below the epidemiological cutoff (ECOFF), above which the bacteria are likely to have an acquired or mutational resistance to enrofloxacin. We developed a differential equation model for the antimicrobial drug concentrations in plasma and colon, and bacteria populations in the feces. The model was fit to animal data of drug concentrations in the plasma and colon obtained using ultrafiltration probes. Fecal E. coli counts and minimum inhibitory concentrations were measured for the week after receiving the antimicrobial drug. We predict that the antimicrobial susceptibility of the bacteria above the ECOFF pre-treatment strongly affects the composition of the bacteria following treatment. Faster removal of the antimicrobial drugs from the colon throughout the study leads to improved clearance of bacteria above the ECOFF in the low dose regimen. If we assume a fitness cost is associated with bacteria above the ECOFF, the increased fitness costs leads to reduction of bacteria above the ECOFF in the low dose study. These results suggest the initial E. coli susceptibility is a strong indicator of how steers respond to antimicrobial drug treatment.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Mathematical Modeling of Multi-Phenomena Anisotropic Systems: Ejection of Primary Aerosols during the Fast Pyrolysis of Biomass

This study introduces a novel particle model for biomass fast pyrolysis, incorporating an anisotropic cylindrical particle to address mass and energy transport coupled with aerosol ejection, which previous models have overlooked. The main contribution lies in developing a model that considers aerosol generation in anisotropic cylindrical particles for the first time, addressing bubbling dynamics and bursting within the liquid phase. The population balance equation describes bubble dynamics and aerosol formation, capturing phenomena like nucleation, growth, coalescence, and bursting. The model employs the method of moments with bubble volume as an internal variable, substantially reducing computational costs by eliminating dependence on this variable. Results highlight the significant impact of anisotropy and particle size on aerosol ejection: smaller, less elongated particles experience faster heating, quicker conversion, and the increased accumulation of the liquid intermediate phase. Specifically, 1 mm diameter particles yield higher concentrations of metaplast and bio-oil aerosols, exceeding 15%, compared to concentrations below 11% for 3 mm particles. This model provides insights into aerosol structure (volume, surface area), aiding in understanding aerosol reactivity at the reactor scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗