Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mathematics”

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 973 records · Page 54

Dose optimization of an adjuvanted peptide-based personalized neoantigen melanoma vaccine

The advancements in next-generation sequencing have made it possible to effectively detect somatic mutations, which has led to the development of personalized neoantigen cancer vaccines that are tailored to the unique variants found in a patient’s cancer. These vaccines can provide significant clinical benefit by leveraging the patient’s immune response to eliminate malignant cells. However, determining the optimal vaccine dose for each patient is a challenge due to the heterogeneity of tumors. To address this challenge, we formulate a mathematical dose optimization problem based on a previous mathematical model that encompasses the immune response cascade produced by the vaccine in a patient. We propose an optimization approach to identify the optimal personalized vaccine doses, considering a fixed vaccination schedule, while simultaneously minimizing the overall number of tumor and activated T cells. To validate our approach, we perform in silico experiments on six real-world clinical trial patients with advanced melanoma. We compare the results of applying an optimal vaccine dose to those of a suboptimal dose (the dose used in the clinical trial and its deviations). Our simulations reveal that an optimal vaccine regimen of higher initial doses and lower final doses may lead to a reduction in tumor size for certain patients. Our mathematical dose optimization offers a promising approach to determining an optimal vaccine dose for each patient and improving clinical outcomes.

59 BASIC BIOLOGICAL SCIENCES↗

ASCR@40: Highlights and Impacts of ASCR’s Programs

A report compiled by the ASCAC Subcommittee on the 40-year history of ASCR for the U.S. Department of Energy’s Office of Advanced Scientific Computing Research. The Office of Advanced Scientific Computing Research (ASCR) sits within the Office of Science in the Department of Energy (DOE). Per their web pages, “the mission of the ASCR program is to discover, develop, and deploy computational and networking capabilities to analyze, model, simulate, and predict complex phenomena important to the DOE.” This succinct statement encompasses a wide range of responsibilities for computing and networking facilities; for procuring, deploying, and operating high performance computing, networking, and storage resources; for basic research in mathematics and computer science; for developing and sustaining a large body of software; and for partnering with organizations across the Office of Science and beyond. While its mission statement may seem very contemporary, the roots of ASCR are quite deep—long predating the creation of DOE. Applied mathematics and advanced computing were both elements of the Theoretical Division of the Manhattan Project. In the early 1950s, the Manhattan Project scientist and mathematician John von Neumann, then a commissioner for the AEC (Atomic Energy Commission), advocated for the creation of a Mathematics program to support the continued development and applications of digital computing. Los Alamos National Laboratory (LANL) scientist John Pasta created such a program to fund researchers at universities and AEC laboratories. Under several organizational name changes, this program has persisted ever since, and would eventually grow to become ASCR.

97 MATHEMATICS AND COMPUTING↗

ASCR@40: Four Decades of Department of Energy Leadership in Advanced Scientific Computing Research

Throughout its long history, the Office of Advanced Scientific Computing Research (ASCR) has built the critical technologies to ensure U.S. leadership in energy science and national security. It has made its parent agency, the Department of Energy (DOE) and its Office of Science, the world’s recognized leader in computational science. ASCR’s stated mission is “to discover, develop, and deploy computational and networking capabilities to analyze, model, simulate, and predict complex phenomena important to the DOE.” To accomplish this goal, ASCR oversees a large complex of computing and networking facilities and is responsible for procuring, deploying and operating high-performance computing (HPC), networking and storage resources; conducting basic research in mathematics and computer science; developing and sustaining a large body of software; and collaborating with other Office of Science programs, academia and industry. ASCR’s computational science leadership has a long history, predating even DOE’s inception. Applied mathematics and advanced computing were both elements of the Manhattan Project’s Theoretical Division. In the 1950s, DOE’s predecessor, the Atomic Energy Commission, created a mathematics program to develop and apply digital computing by supporting researchers at universities and AEC laboratories. Several organizational and name changes later, this program would grow and become ASCR.

97 MATHEMATICS AND COMPUTING↗

Autonomous Discovery in Science and Engineering

The Center for Advanced Mathematics for Energy Research Applications (CAMERA) is an integrated, cross-disciplinary center aimed at inventing, developing, and delivering the fundamental new mathematics required to capitalize on experimental investigations at scientific facilities. Jointly funded by the Office of Advanced Scientific Computing Research (ASCR) and the Office of Basic Energy Sciences (BES) within the US Department of Energy’s Office of Science, CAMERA identifies areas in experimental science that can be aided by new mathematical insights, develops the needed algorithmic tools, and delivers them as user-friendly software to the experimental community.

97 MATHEMATICS AND COMPUTING↗

Discovering new governing equations using ML

A hallmark of the scientific process since the time of Newton has been the derivation of mathematical equations meant to capture relationships between observables. As the field of mathematical modeling evolved, practitioners specifically emphasized mathematical formulations that were predictive, generalizable, and interpretable. Machine learning’s ability to interrogate complex processes is particularly useful for the analysis of highly heterogeneous, anisotropic materials where idealized descriptions often fail. As we move into this new era, we anticipate the need to leverage machine learning to aid scientists in extracting meaningful, but yet sometimes elusive, relationships between observed quantities.

97 MATHEMATICS AND COMPUTING↗

ExaSGD: 2022 Kernel Thrust Activities

The Kernel Thrust milestone ADSE22-407 covers the development of device-capable optimization algorithms and solvers technologies required by the ExaSGD project’s software stack in order to solve security-constrained alternating current optimal power flow (SC-ACOPF) problems on emerging exascale architectures. To this extent, in FY22 the main objective of the Kernel Thrust was (i) provide sparse optimization solver that runs efficiently on hardware accelerator devices (i.e., NVIDIA and AMD GPUs) to perform intra-node computations, (ii) strengthen the reliability and increase the performance of the mixed-dense sparse (MDS) solver of HiOp for deployment on the FY22 target architectures, Summit and Crusher, and (iii) increase performance by improving the mathematical algorithm and refining the parallel MPI-based implementation of the coarse-grain parallel solver HiOp-PriDec for capabilities deployment on the FY22 target architectures, Summit and Crusher. This document presents the developments and contributions done by the Kernels Thrust Team in FY22 toward completion of the above-mentioned objectives. These contributions progressed along four main development (sub)thrusts: (1) Design and implementation of a sparse optimization solver for use on hardware accelerators; (2) Improvement of the mathematical algorithm and of the parallel implementation of HiOp-PriDec to ensure readiness and efficient coarse-grain parallelism for FY23 target exascale machine; and (3) Support Software and Application Development Thrusts of the exaSGD project in their deployment of the project’s software stack on AMD- and NVIDIA-based architectures. The development of the sparse optimization solver (thrust 1 above) was new in FY22 and resulted in a new sparse solver in HiOp (available as of version 0.6). The second development thrust was a continuation of the efforts from FY21 and improved the mathematical algorithm and the communication strategy of the HiOp-PriDec solver. The last developement thrust is a large collaborative effort. Namely, the project’s teams from multiple labs (LLNL, PNNL, ORNL, and NREL) performed large-scale demonstration of the ExaSGD software stack, namely the optimization solvers of HiOp interfaced with the modeling front-end ExaGO and the stochastic sampler PowerScenarios. These demonstration efforts solved large-scale instances of the SC-ACOPF challenge problem of medium network sizes (10, 000-bus system) and large number of contingencies on Summit (NVIDIA accelerators) and Crusher (AMD accelerators) systems at ORNL.

97 MATHEMATICS AND COMPUTING↗

Quantification of Type I Interferon Inhibition by Viral Proteins: Ebola Virus as a Case Study

Type I interferons (IFNs) are cytokines with both antiviral properties and protective roles in innate immune responses to viral infection. They induce an antiviral cellular state and link innate and adaptive immune responses. Yet, viruses have evolved different strategies to inhibit such host responses. One of them is the existence of viral proteins which subvert type I IFN responses to allow quick and successful viral replication, thus, sustaining the infection within a host. We propose mathematical models to characterise the intra-cellular mechanisms involved in viral protein antagonism of type I IFN responses, and compare three different molecular inhibition strategies. We study the Ebola viral protein, VP35, with this mathematical approach. Approximate Bayesian computation sequential Monte Carlo, together with experimental data and the mathematical models proposed, are used to perform model calibration, as well as model selection of the different hypotheses considered. Finally, we assess if model parameters are identifiable and discuss how such identifiability can be improved with new experimental data.

59 BASIC BIOLOGICAL SCIENCES↗

The Weak Form Is Stronger Than You Think

The weak form is a ubiquitous, well-studied, and widely-utilized mathematical tool in modern computational and applied mathematics. In this work we provide a survey of both the history and recent developments for several fields in which the weak form can play a critical role. In particular, we highlight several recent advances in weak form versions of equation learning, parameter estimation, and coarse graining, which offer surprising noise robustness, accuracy, and computational efficiency. We note that this manuscript is a companion piece to our October 2024 SIAM News article of the same name. Here we provide more detailed explanations of mathematical developments as well as a more complete list of references. Lastly, we note that the software with which to reproduce the results in this manuscript is also available on our group's GitHub website https://github.com/MathBioCU .

26A33, 35D30, 62FXX, 62JXX, 65L09, 65M32, 68Q32,↗

Computational Math Problems for a Clean Energy Future

Cutting edge computational mathematics are ubiquitous in renewable energy research. Problems in resilient and reliable electric grid operations, infrastructure planning, wind farm yaw control, and more demand sophisticated and scalable computational tools that enable the transition of renewable energy technologies from proof of concept to deployment into our energy system. The mission of the Computational Science Center at NREL is to lead the lab's efforts to solve energy challenges using high-performance computing (HPC), computational science, applied mathematics, scientific data management, visualization, and informatics. In this poster, we provide a short overview of three areas of computational mathematics research at NREL: wind power scenario generation for stochastic grid operations and infrastructure planning, improved rational function approximations for electromagnetic transients codes, and wind farm yaw control using a combination of the Alternating Direction Method of Multipliers (ADMM) and reinforcement learning (RL). Increasing penetrations of renewable energy into power grids motivate the investigation of new approaches to characterizing uncertainty for five-minute economic dispatch problems. Similarly, as the penetration of distributed energy resources on power grids increases, it becomes important to revisit our methods of modelling transient phenomena, i.e. electromagnetic transients programs. Finally, the combination of ADMM and RL for wind farm yaw control presented here can potentially increase the efficiency of the deployed distributed controllers by orders of magnitude.

ADMM↗

PCET‐Driven Reactivity of Neptunyl(VI) Yields Oxo‐Bridged Np(V) and Np(IV) Species

Two unconventional polynuclear complexes of neptunium (Np) featuring mono-mathematical equation -oxo motifs have been accessed by proton-coupled electron transfer (PCET) reactivity involving the dissolution of neptunyl(VI) diacetate dihydrate (NpO 2 (OAc) 2 (H 2 O) 2 ∙ HOAc) in methanol followed by addition of a pentadentate Schiff-base ligand. One complex is a mixed-valent [Np V ,Np IV , Np V ] trimer with two bridging mathematical equation μ 2 -oxos and the other is a [Np V , Np V ] dimer featuring a single mathematical equation μ 2 -oxo. In both complexes the outer Np centers are capped with terminal oxo ligands. Spectroscopic and spectrokinetic studies aimed at elucidating mechanistic details of complex formation in this system show that intermediate multinuclear [Np V O 2 (OAc)] n species form prior to metal chelation by the ligand; electrolysis experiments demonstrate that production of Np(V) gives rise to asynchronous proton transfer that does not occur otherwise (in the Np(VI) state) as well as condensation with loss of H 2 O and formation of the polynuclear complexes. We attribute the oxo-deficient nature of these products, with respect to conventional actinyl ([AnO 2 ] m+ ) species, to the reduction/condensation reaction sequence of PCET.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Some Problems in Density Functional Theory

Though calculations based on density functional theory (DFT) are used remarkably widely in chemistry, physics, materials science, and biomolecular research and though the modern form of DFT has been studied for almost 60 years, some mathematical problems remain. From a physical science perspective, it is far from clear whether those problems are of major import. For context, we provide an outline of the basic structure of DFT as it is presented and used conventionally in physical sciences, note some unresolved mathematical difficulties with those conventional demonstrations, then pose several questions regarding both the time-independent and time-dependent forms of DFT that could benefit from attention in applied mathematics. Furthermore, progress on any of these would aid in development of better approximate functionals and in interpretation of DFT.

36 MATERIALS SCIENCE↗

Mechanical systems and kinematics

Kinematics is the study of how motion is achieved without reference to the forces that create the motion. This includes the basic structures of all robotic systems, such as links and joints, as well as their mathematical representations. This chapter will provide a brief overview of the types of motion used in additive manufacturing (AM) systems, as well as the mathematical equations that govern them. Rotations are explained using both Euler rotations and quaternions. Combined rotation and translations are shown using homogeneous transform matrices. The robot Jacobian, which relates joint velocities to the end effector translational velocities, is explained with an example calculation. This section will give the reader the basic mathematical foundation in kinematics that is necessary to understand the mechanical underpinnings of path planning for AM systems.

Wang, Peter↗

Assessment the leaching characteristics and long-term leaching behavior of some radionuclides from synthesized zeolite cement matrix

Highlights: • Synthesis and characterization of zeolite cement synthesized from fly ash • Leaching tests were used to assess the release rate of {sup 85}Sr and {sup 60}Co radionuclides. • Simple mathematical models were used to detect leaching mechanisms. • Long-term leaching behavior of {sup 85}Sr and {sup 60}Co from cementitious matrix was assessed. The development of cementitious materials remains a vital goal to produce valuable products with good mechanical, physical and chemical properties suitable for the safe disposal of concentrated radionuclides resulted from the treatment of contaminated solutions. In this paper, zeolite cement with good properties was laboratory synthesized from industrial by-products fly ash and characterized using various characterization techniques. The leaching characteristics of strontium and cobalt radionuclides from the synthesized zeolite cement were scrutinized conforming with the standard leaching methodology of International Atomic Energy Authority (IAEA). A mechanical strength assessment was executed to characterize the extent of immobilization process of the solidified matrices. The cumulative leaching fraction (φ) of the two studied radionuclides was found to be less than 5% in all examined conditions, which implies the applicability of the IAEA-recommended methodology for estimating the diffusion coefficient. The experimental leaching data were regressed nonlinearly to various mathematical kinetic models to assess the controlling leaching mechanism and to determine the leaching parameters. The regression results indicated that strontium ({sup 85}Sr) and cobalt ({sup 60}Co) leaching resulted from two succeeding mechanisms: as first order kinetic reaction, and then diffusion. The calculated values of leachability indices signify that the performance of each of the studied matrices is within an acceptable range. A simplified mathematical model, rooted in the first order reaction and diffusion mechanisms, was simulated to predict the radionuclides leaching rates from zeolite cement matrices. By comparing the synthesized zeolite cement with other cementitious materials, it can be concluded that the synthesized material can be classified as an efficient material suitable to immobilize {sup 85}Sr and {sup 60}Co from radioactive wastes. The acquired findings demonstrated that the studied immobilized waste matrices have acceptable mechanical effectiveness.

36 MATERIALS SCIENCE↗

Optimal PMU design based on sampling model and sensitivity analysis

The precise measurements of the synchrophasor and frequency from phasor measurement units (PMUs) are widely used in power grid applications. With the improvement of the technique, applications always require stable dynamic performance and higher accuracy for the synchrophasor and frequency measurements, which is challenging for PMU development. To evaluate the contribution of PMU hardware to measurement accuracy, this paper proposes a general-purpose sampling model to analyze the measurement error. In the proposed sampling model, a strict mathematical derivation is derived, where its error is purely determined by the parameters of the PMU hardware. The sensitivity analysis is carried out by three methods, including mathematical analysis, computer simulation, and variance-based sensitivity analysis. Through the sensitivity analysis, this paper establishes the systematic formulation and the inclusion of synchrophasor, frequency, and ROCOF. Experimental results based on the real-world testbench involving distribution-level PMUs match the mathematical analysis conclusion, which verifies the correctness of the general-purpose sampling model. Furthermore, a strategy for the optimal PMU design is proposed, which could guide PMU design in the future.

42 ENGINEERING↗

On principles of emergent organization

After more than a century of concerted effort, physics still lacks basic principles of spontaneous organization. To appreciate why, we first state the problem, outline historical approaches, and survey the present state of the physics of self-organization. This frames the particular challenges arising from mathematical intractability and the resulting need for computational approaches, as well as those arising from a chronic failure to define structure. Then an overview of two modern mathematical formulations of organization—intrinsic computation and evolution operators—lays out a way to overcome these challenges. Together, the vantage point they afford shows how to account for the emergence of structured states via a statistical mechanics of systems arbitrarily far from equilibrium. The result is a constructive path forward to principles of organization that builds on mathematical identification of structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Predicting Characteristics of the Water Cycle From Scaling Relationships

Over multi-decadal time scales, assuming that changes in subsurface water storage are negligible, the continental precipitative water flux, P, can be divided into two principal components, Q (run-off, including soil infiltration and groundwater recharge) and ET (evapotranspiration). Taking into account a broadly applied Budyko's phenomenology to describe the relationship of ET/P as a function of PET/P, where PET is the potential evapotranspiration, we propose a theoretical framework for predicting characteristics of the water cycle from scaling relationships. In this framework, the ecosystem net primary productivity is expressed in terms of soil formation and vegetation growth, which is mathematically optimized with respect to the water partitioning, generating directly the value ET/P. Here, the mathematical optimization is based on the general ecological principle that dominant ecosystems tend to be those that, for any given conditions, maximize conversion of atmospheric carbon to biomass. It is shown that application of the results of mathematical optimization to water-limited ecosystems is possible by applying the optimization only to a vegetation covered portion of the surface. For energy-limited ecosystems, the optimization can be applied only to a portion of the precipitation equal to PET, assuming that the remaining P simply runs off. We use theoretical and actual values of plant root fractal dimensionalities, d f , to predict ranges of ET/P as a function of PET/P for 0 ≤ PET/P ≤ 1 and compare with annual and multi-decadal means of ET/P. By comparing the developed approach with a large amount of data collected from the literature, we demonstrate its successful applications to both water- and energy-limited systems.

54 ENVIRONMENTAL SCIENCES↗

The cluster decomposition of the configurational energy of multicomponent alloys

Abstract The cluster expansion method (CEM) is a widely used lattice-based technique in the study of multicomponent alloys. Despite its prevalent use, a clear understanding of expansion terms is lacking. We present a modern mathematical formalism of the CEM and introduce thecluster decomposition—a unique and basis-independent decomposition for functions of the atomic configuration in a crystal. We identify the cluster decomposition as an invariant ANOVA decomposition; and demonstrate how functional analysis of variance and sensitivity analysis can be used to interpret interactions among species. Furthermore, we show how the mathematical structure of the cluster decomposition enables numerical evaluation that scales with the number of clusters and is independent of the number of species. Overall, our work enables rigorous interpretations of interactions among species, provides opportunities to explore parameter estimation beyond linear regression, introduces a numerical efficient implementation, and enables analysis of cluster expansions based on established mathematical and statistical principles.

Chemistry↗

Individualized empirical baselines for evaluating the energy performance of existing buildings

The evaluation of building energy performance requires a baseline for comparison. Common empirical baselines are usually used for existing buildings since they are fast and convenient. However, the same type of building at the same location will receive the same baseline despite their difference in usage. Individualized baselines by creating building energy models are possible solutions, but it is labor intensive and time-consuming. To fill the gap, this study is to develop individualized empirical baselines for existing buildings in a fast way. First, common empirical baselines are created based on survey data. Then, to get training samples, building energy models for large-scale existing buildings are created and simulated. So finally, based on simulation results, mathematical models to get individualized empirical baselines in a fast way are created. U.S. medium office buildings were used as an example to demonstrate the method. We developed 30 mathematical models for medium office buildings in two vintages (constructed before 1980 and after 1980) and 15 climate zones. The mean absolute percentage errors (MAPE) between the individualized empirical baselines and the modeled baselines for those 30 mathematical models are all lower than 5.5%. An engineer can obtain the individualized empirical baseline for an existing building in a few seconds by using the open-source tool we developed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗