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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Relative permeabilities for two-phase flow through wellbore cement fractures

Multiple fluids are likely to exist in fractures and flow paths associated with leaky wellbores, including liquids (e.g., crude oil) and gases (e.g., gas exsolved from liquid). These fluids occupy and move through different portions of the pore spaces within the fractures depending on many factors, including fluid properties, fracture size, and the amount of the different fluids. Upward leakage of any phase, through the fracture, can contaminate water-bearing formations, create hazardous surface conditions, and compromise the functionality of the wellbore. Early signs of wellbore leaks may be expressed by anomalous pressure behavior at surface monitoring points on cavern storage wells. These pressure anomalies are difficult to interpret, necessitating knowledge of the factors that affect the multiphase flow in fractures and porous media. These parameters are critical to modeling multiphase flow in fractures. This insight can guide further diagnosis and maximize leak remediation. Here, our study focuses on the relationship of the liquid–gas relative permeabilities for representative variable-aperture wellbore cement fracture. To obtain the relative permeability of each phase, two-phase flow tests were conducted where both fluids were flowing simultaneously through a fractured wellbore cement specimen under a range of factors, namely (1) aperture size, (2) capillary numbers, and (3) viscosity ratio. The flow experiments were conducted under a range of confining stresses and flow velocities, using nitrogen gas and silicone oils (of different viscosities) in a specially designed pressure vessel. The sum of gas and oil relative permeabilities were found to be less than one under all conditions, which indicates that the presence of one phase affects the permeability of the other phase, and vice versa. Since the gas phase flow conditions include a significant inertial flow component in addition to viscous flow, the inertial flow coefficients at different saturation states are presented. The factors affecting the relationship between the relative permeabilities are discussed in detail. A new mathematical model for estimating the relative permeability of wellbore cement fracture is presented and experimentally validated.

58 GEOSCIENCES↗
DOE OSTI2025

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗
DOE OSTI2026

An information-matching approach to optimal experimental design and active learning

The efficacy of mathematical models heavily depends on the quality of the training data, yet collecting sufficient data is often expensive and challenging. Many modeling applications require inferring parameters only as a means to predict other quantities of interest (QoI). Because models often contain many unidentifiable (sloppy) parameters, QoIs often depend on a relatively small number of parameter combinations. Therefore, we introduce an information-matching criterion based on the Fisher information matrix to select the most informative training data from a candidate pool. This method ensures that the selected data contain sufficient information to learn only those parameters that are needed to constrain downstream QoIs. It is formulated as a convex optimization problem, making it scalable to large models and datasets. Here, we demonstrate the effectiveness of this approach across various modeling problems in diverse scientific fields, including power systems and underwater acoustics. Finally, we use information-matching as a query function within an active learning (AL) loop for materials science applications. In all these applications, we find that a relatively small set of optimal training data can provide the necessary information for achieving precise predictions. These results are encouraging for diverse future applications, particularly AL in large machine-learning models.

Materials science↗
DOE OSTI2024

Deep Nonparametric Estimation of Operators between Infinite Dimensional Spaces

Learning operators between infinitely dimensional spaces is an important learning task arising in machine learning, imaging science, mathematical modeling and simulations, etc. This paper studies the nonparametric estimation of Lipschitz operators using deep neural networks. Non-asymptotic upper bounds are derived for the generalization error of the empirical risk minimizer over a properly chosen network class. Under the assumption that the target operator exhibits a low dimensional structure, our error bounds decay as the training sample size increases, with an attractive fast rate depending on the intrinsic dimension in our estimation. Our assumptions cover most scenarios in real applications and our results give rise to fast rates by exploiting low dimensional structures of data in operator estimation. We also investigate the influence of network structures (e.g., network width, depth, and sparsity) on the generalization error of the neural network estimator and propose a general suggestion on the choice of network structures to maximize the learning efficiency quantitatively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A Convex Data-Driven Approach for Nonlinear Control Synthesis

We consider a class of nonlinear control synthesis problems where the underlying mathematical models are not explicitly known. We propose a data-driven approach to stabilize the systems when only sample trajectories of the dynamics are accessible. Our method is built on the density-function-based stability certificate that is the dual to the Lyapunov function for dynamic systems. Unlike Lyapunov-based methods, density functions lead to a convex formulation for a joint search of the control strategy and the stability certificate. This type of convex problem can be solved efficiently using the machinery of the sum of squares (SOS). For the data-driven part, we exploit the fact that the duality results in the stability theory can be understood through the lens of Perron–Frobenius and Koopman operators. This allows us to use data-driven methods to approximate these operators and combine them with the SOS techniques to establish a convex formulation of control synthesis. The efficacy of the proposed approach is demonstrated through several examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A Hierarchical Optimization Method for Electric Vertical Takeoff and Landing Aircraft Network Design

Electric vertical takeoff and landing aircraft (eVTOLs) are expected to serve urban air mobility in a station-to-station configuration, which makes the optimal network design of eVTOL stations a critical question to explore. Existing approaches often face limitations, such as the inability to interact station locations with demand or difficulty in finding the optimal solution for large study regions. Here, this paper first proposes a mathematical model to generate optimal eVTOL station locations while considering associated potential eVTOL demand, and then proposes a heuristic algorithm, Hierarchical Optimization MEthod (HOME), to efficiently solve the model. With a case study of Southern California, HOME was compared to 1) directly solving the original integer linear programming-based network design problem, and 2) employing the widely used genetic algorithm. Results suggest that HOME can find optimal solutions with limited computational resources. The proposed framework powered by HOME provides a computationally efficient way to support urban air mobility planning.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Deep nonparametric estimation of operators between infinite dimensional spaces

Learning operators between infinitely dimensional spaces is an important learning task arising in machine learning, imaging science, mathematical modeling and simulations, etc. This paper studies the nonparametric estimation of Lipschitz operators using deep neural networks. Non-asymptotic upper bounds are derived for the generalization error of the empirical risk minimizer over a properly chosen network class. Under the assumption that the target operator exhibits a low dimensional structure, our error bounds decay as the training sample size increases, with an attractive fast rate depending on the intrinsic dimension in our estimation. Our assumptions cover most scenarios in real applications and our results give rise to fast rates by exploiting low dimensional structures of data in operator estimation. We also investigate the influence of network structures (e.g., network width, depth, and sparsity) on the generalization error of the neural network estimator and propose a general suggestion on the choice of network structures to maximize the learning efficiency quantitatively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Stochastic Learning Approach for Binary Optimization: Application to Bayesian Optimal Design of Experiments

Here, we present a novel stochastic approach to binary optimization suited for optimal experimental design (OED) for Bayesian inverse problems governed by mathematical models such as partial differential equations. The OED utility function, namely, the regularized optimality criterion, is cast into a stochastic objective function in the form of an expectation over a multivariate Bernoulli distribution. The probabilistic objective is then solved by using a stochastic optimization routine to find an optimal observational policy. This formulation (a) is generally applicable to binary optimization problems with soft constraints and is ideal for OED and sensor placement problems; (b) does not require differentiability of the original objective function (e.g., a utility function in OED applications) with respect to the design variable, and thus it enables direct employment of sparsity-enforcing penalty functions such as $\ell_0$, without needing to utilize a continuation procedure or apply a rounding technique; (c) exhibits much lower computational cost than traditional gradient-based relaxation approaches; and (d) can be applied to both linear and nonlinear OED problems with proper choice of the utility function. The proposed approach is analyzed from an optimization perspective with detailed convergence analysis of the optimization approach and is also analyzed from a machine learning perspective with correspondence to policy gradient reinforcement learning. The approach is demonstrated numerically by using an idealized two-dimensional Bayesian linear inverse problem and validated by extensive numerical experiments carried out for sensor placement in a parameter identification setup.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Construction of a Code Verification Matrix for Heat Conduction With Finite Element Code Applications

When establishing the pedigree of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model. During this process, numerical results on various meshes are systematically compared to a reference analytic solution. The selection of analytic solutions can be a laborious process, as it is difficult to establish adequate code confidence without performing redundant work. In this report we address this issue by applying a physics-based process that establishes a set of reference problems. In this process, code simulation options are categorized and systematically tested, which ensures that gaps in testing are easily identified and addressed. The resulting problems are primarily intended for code verification analysis but may also be useful for comparison to other simulation codes, troubleshooting activities, or training exercises. The process is used to select fifteen code verification problems relevant for the one-dimensional steady-state heat conduction equation. These problems are applicable to a wide variety of simulation tools, but, in this work, a demonstration is performed using the finite element-based nuclear fuel performance code BISON. Convergence to the analytic solution at the theoretical rate is quantified for a selection of the problems, which establishes a baseline pedigree for the code. Not only can this standard set of conduction solutions be used for verification of other codes, but also the physics-based process for selecting problems can be utilized to quantify and expand testing for any simulation tool.

42 ENGINEERING↗
DOE OSTI2023

Mechanical Actuation via Homeomorphic Transformations of Topological Solitons within Polymer Coatings

Abstract Topological solitons are currently under investigation for their exotic properties, especially in nonlinear physics, optics, and material sciences. However, challenges of robust generation and limited stability over time have hindered their practical uses. To address this issue, an approach is developed to form structured arrays of solitons in films of polymerizable liquid crystals. Their complex molecular architecture is preserved by in situ photopolymerization forming a stable liquid crystal network. Most excitingly, their properties are advanced to include responsiveness functions. When thermally actuated, these topological solitons mediate the reconfiguration of surface topographies. Complex shape changes occur depending on the intrinsic complex spatial distribution of the director, which may even lead to full shape inversion and topographical changes as high as ≈40% of the initial thickness. Conversely, the shape changes provide information on the initial director profile, which is consistent with the mathematical model. The soliton‐containing polymer coatings are applicable in multiple domains, ranging from tunable optics to haptics, and from shape‐coupled sensing systems to temperature‐coupled heat management.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Particle–Continuum Coupling and its Scaling Regimes: Theory and Applications

Abstract This work is motivated by the goal of designing simulation software for technical devices that, at their functional core, rely on atomistic‐scale processes embedded in a larger‐scale fluid environment. The core of the problem is the conceptual and technical approach for coupling particle and continuum representations of a fluid. The state of the art for key aspects including physical modeling, mathematical formalization, computational implementation, and applications, is discussed and organized in a consistent picture across the relevant physical regimes.

Delle Site, Luigi↗
DOE OSTI2026

Bemnifosbuvir: An HCV NS5B Inhibitor With Multiple Modes of Action

Bemnifosbuvir (BEM) is a potent, pan-genotypic inhibitor targeting the hepatitis C virus (HCV) NS5B polymerase. Its antiviral activity was evaluated in an ascending dose phase I clinical trial involving 30 patients treated once a day for 7 days. After treatment initiation, plasma HCV RNA declined in a biphasic manner with a mean reduction of 2.3 log 10 IU/mL after 24 hours and 4.4 log10 IU/mL by day 7 for the highest dose. Alanine aminotransferase (ALT) also normalized in most patients. We employed a multiscale mathematical model fitted to the HCV RNA and ALT dynamics to quantify the antiviral activity and evaluate the modes of action of BEM. We found that models in which BEM only acted as a typical HCV RNA polymerase inhibitor and reduced the intracellular production of HCV RNA did not fit the data as well as models in which BEM had multiple modes of action, including suppressing viral assembly and secretion and enhancing intracellular HCV RNA degradation. BEM's effectiveness in inhibiting intracellular HCV RNA production increased with dose (150 mg/day: 88.2%, 300 mg/day: 98.8%, 600 mg/day: 99.5%), while inhibition of viral assembly and release was ~95% effective regardless of dose. We observed a dose-dependent enhancement in the degradation of intracellular HCV RNA, with degradation rates 1.5-fold higher in patients receiving 300 mg/day and 2.7-fold higher in those receiving 600 mg/day than in patients receiving 150 mg/day. No significant differences in antiviral activity were detected between HCV genotypes 1b and 3 or between patients with and without compensated cirrhosis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Microwave‐assisted hydrodistillation of the essential oil from Algerian Pimpinella anisum seeds

Abstract Essential oils (EOs) from anise seeds ( Pimpinella anisum ) collected in the Algerian region were extracted by microwave‐assisted hydrodistillation (MAHD) and hydrodistillation (HD). The oil yield (%), chemical composition, energetic consumption and the antimicrobial properties were evaluated of the obtained EO. The yields obtained by both processes were 2.81 ± 0.01% for 90 minutes and 3.30 ± 0.05% for 330 minutes for MAHD and HD, respectively. The determination of the chemical composition of EOs by GC‐MS revealed the presence of a majority compound, trans ‐anethole 81.52% and 91.31% for MAHD and HD, respectively. The other main compounds present in the oils for MAHD and HD were, respectively: methyl chavicol (estragole) (0.63%‐1.39%), cis ‐anethole (0.23%‐0.34%), anisaldehyde (1.1%‐1.57%), anisyl methyl ketone (0.28%‐0.13%), trans ‐methyl isoellgenol (0.31%‐0.17%), γ ‐himachalene (0.5%‐0.36%), trans ‐pseudoisoeugenyl 2‐methylbutyrate (6.53%‐2.53%) and epoxy‐pseudoisoeugenyl 2‐methylbutyrate (2.01%‐0.56%). The EO obtained by MAHD detected other compounds, fenchone 0.41%, limonene 0.68%, cis ‐jasmolacatone extra C 0.41%, trans ‐isolongifolanone 0.2% and α ‐himachalene 0.16%. The energetic consumption by MAHD was significantly less than HD. Theoretical data from mathematical modelling were compared with practical data, and the coefficient of determination was between 0.93 and 0.97. The sensibility of bacterial strains to the essential oil obtained by MAHD is marked relatively to that obtained by HD, and for the sensibility of fungal strains, it was conserved. The MAHD process revealed a good yield of anise seeds EOs while reducing extraction time and energetic consumption. The EO extracted by this process revealed the presence of other compounds while improving antibacterial activity.

Boumahdi, Yamina↗
DOE OSTI2025

Three-Receiver Quantum Broadcast Channels: Classical Communication with Quantum Non-unique Decoding

In network communication, it is common in broadcasting scenarios for there to exist a hierarchy among receivers based on information they decode due, for example, to different physical conditions or premium subscriptions. This hierarchy may result in varied information quality, such as higher-quality video for certain receivers. This is modeled mathematically as a degraded message set, indicating a hierarchy between messages to be decoded by different receivers, where the default quality corresponds to a common message intended for all receivers, a higher quality is represented by a message for a smaller subset of receivers, and so forth. We extend these considerations to quantum communication, exploring three-receiver quantum broadcast channels with two- and three-degraded message sets. Our technical tool involves employing quantum non-unique decoding, a technique we develop by utilizing the simultaneous pinching method. Here, we construct one-shot codes for various scenarios and find achievable rate regions relying on various quantum Rényi mutual information error exponents. Our investigation includes a comprehensive study of pinching across tensor product spaces, presenting our findings as the asymptotic counterpart to our one-shot codes. By employing the non-unique decoding, we also establish a simpler proof to Marton’s inner bound for two-receiver quantum broadcast channels without the need for more involved techniques. Additionally, we derive no-go results and demonstrate their tightness in special cases.

Salek, Farzin [Technical University of Munich (Ger↗
DOE OSTI2022

Transport and mechanical behavior in PEO-LLZO composite electrolytes

Abstract Composite solid electrolytes (CEs), wherein ion-conducting polymer and ceramic/glass is mixed, are promising candidates for all-solid-state batteries due to their promise of acceptable ionic conductivity and mechanical properties compared to their individual constituents. While numerous studies have focused on improving the performance of CEs, it is still unclear what the material targets are that can result in improved macroscopic performance especially in light of the coupled needs for high transport and high mechanical strength in these materials. In this study, a two-dimensional (2D) mathematical model is developed to investigate electrochemical and mechanical characteristics of CEs. The model is compared to CEs consisting of poly-ethylene-oxide (PEO) polymer and lithium lanthanum zirconium oxide (LLZO) ceramic material with examination of the impact of varying LLZO volume fractions. The potential drop at the PEO-LLZO interface is evaluated using the junction potential theory. Using experimental data from the literature, the model estimates the ionic conductivity, effective transference number, and mechanical stiffness of the CEs. While the mechanical stiffness improves with increasing volume fraction of LLZO, the impact on conductivity and transference number depends on interfacial resistance at the interface. Finally, the study reports CE’s potential to enhance Li-ion transport and mechanical properties to inhibit lithium (Li) dendrite growth.

25 ENERGY STORAGE↗
DOE OSTI2022

Physics-informed neural networks (PINNs) for fluid mechanics: a review

Despite the significant progress over the last 50 years in simulating flow problems using numerical discretization of the Navier–Stokes equations (NSE), we still cannot incorporate seamlessly noisy data into existing algorithms, mesh-generation is complex, and we cannot tackle high-dimensional problems governed by parametrized NSE. Moreover, solving inverse flow problems is often prohibitively expensive and requires complex and expensive formulations and new computer codes. Here, we review flow physics-informed learning, integrating seamlessly data and mathematical models, and implement them using physics-informed neural networks (PINNs). Here, we demonstrate the effectiveness of PINNs for inverse problems related to three-dimensional wake flows, supersonic flows, and biomedical flows.

42 ENGINEERING↗
DOE OSTI2022

Thermal Pressure in the Thermal Equation of State for Solid and a Proposed Substitute

Abstract The thermal equation of state (TEOS) for solids is a mathematic model among pressure, temperature and density, and is essential for geophysical, geochemical, and other high pressure–temperature (high P–T) researches. However, in the last few decades, there has been a growing concern about the accuracy of the pressure scales of the calibrants, and efforts have been made to improve it by either introducing a reference standard or building new thermal pressure models. The existing thermal equation of state, P ( V , T ) = P ( V , T 0 ) + P th ( V , T ), consists of an isothermal compression and an isochoric heating, while the thermal pressure is the pressure change in the isochoric heating. In this paper, we demonstrate that, for solids in a soft pressure medium in a diamond anvil cell, the thermal pressure can neither be determined from a single heating process, nor from the thermal pressure of its calibrant. To avoid the thermal pressure, we propose to replace the thermal pressure with a well-known thermal expansion model, and integrate it with the isothermal compression model to yields a Birch–Murnaghan-expansion TEOS model, called VPT TEOS. The predicted pressure of MgO and Au at ambient pressure from Birch–Murnaghan-expansion VPT TEOS model matches the experimental pressure of zero (0) GPa very well, while the pressure prediction from the approximated Anderson PVT TEOS exhibit a big deviation and a wrong trend.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Methodologies and strategies for detecting fiber orientation in polymeric fiber-reinforced composites

In this review paper different methods and techniques for representing fiber orientation in advanced polymer matrix composites are presented. A description of the effect of fiber orientation on processing, fabrication and mechanical properties of a composite material is presented. The paper discusses the mathematical modeling and modeling techniques for fiber structure representation in polymeric composites, with an emphasis on fiber orientation in a composite part. The paper is divided into two sections; namely—destructive techniques and non-destructive techniques. The capabilities and limitations of each technique with respect to fiber orientation detection and measurements are discussed.

composites↗