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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Attached Jet Velocity Profiles in Mixing Tanks

Free jets have been studied in detail over much of the last century, but the theory for offset and attached jets remains incomplete. Attached jets differ from free jets in that they lose momentum to nearby surfaces, attenuating their velocities. The velocity profiles of free circular jets are nearly Gaussian, with quantitative mathematical descriptions derived from first principles by Goertler and Tollmien (Rajaratnam, 1976). In contrast, mathematical descriptions of three-dimensional attached jets from circular nozzles remain much less mature. Agelin-Chaab and Tachie (2011) used particle imaging velocimetry of a three-dimensional attached jet to show that the scaled velocity decays with scaled distance from the nozzle with a power law exponent between -1.15 and -1.21, which is larger in magnitude than that of a free jet. However, quantitative analytical expressions for the velocity profiles of attached jets similar to those of free jets remain elusive. This paper addresses this critical gap. Here we evaluate the velocity profiles of three-dimensional offset jets emerging from circular nozzles that that become attached jets. These jets lose momentum due to interactions with nearby surfaces and are important to evaluating flows in mixing vessels and to suspending solids and trapped gases in radioactive waste tanks. Despite the importance of attached jets, prior insight has been purely experimental, limited to overly simplistic analytical models, or restricted to computationally expensive computational fluid dynamics case studies. We compare the expression of Verhoff (1963) to experimental results to find reasonable quantitative agreement. As stated by Agelin-Chaab and Tachie (2011), “detailed velocity measurements of 3D offset jets are rare.” Such remains the case. This study adds to the literature by providing information at two additional Reynolds numbers (1.43 106 and 1.87 106) and evaluating simple but accurate expressions for velocity profiles. These Reynolds numbers are an order of magnitude higher than prior reports. The semi-empirical stream wise velocity profile perpendicular to the surface proposed by Verhoff (1963) is in approximate agreement with these velocity profiles, which is surprising because these attached jets are three-dimensional instead of two-dimensional as evaluated by Verhoff. However, additional work is necessary to fully describe these profiles quantitatively.

jet, attached jet, offset jet↗
DOE OSTI2022

Critical Shear Stress for Erosion Under Laminar Jet Flow

Particle erosion is an important process in both natural and manmade environments. For example, in natural environments particle or sediment erosion rates determine the course of rivers, depth of streams, and fate of deltas. The interplay between erosion rates and flow patterns can have dramatic impact on shorelines. Where hydrodynamic jet flows or rip currents change direction, sediments may accumulate or disperse altering the shape of coastlines and the ecosystems and economies that rely on them. In man-made systems, jet flows may be used to mobilize particles from surfaces or within industrial scale mixing tanks. For example, in the nuclear industry, impinging jet mixers may be used to mobilize and suspend sludges at the bottom of waste tanks. In each of these examples, the Shields diagram is useful to determine the conditions under which erosion occurs. This diagram plots the scaled critical shear stress for erosion as a function of the particle size as embedded in particle Reynolds numbers, Archimedes numbers, or their equivalents. Although the data, that carries a large spread, is largely in agreement with mathematical models, for the smallest particles, model predictions diverge from each other and from experimental data. Indeed, there is substantial uncertainty about the Shields diagrams at relatively small particle sizes and flow rates. Here we evaluate critical shear stresses for erosion under laminar jet flow. To date critical shear stresses for erosion, as embodied in Shields diagrams, have largely been developed for turbulent flow both experimental and mathematically. However, curves on the Shields diagrams diverge among various models at lower particle Reynolds numbers. Indeed, how critical shear stresses for erosion develop under fully laminar jet flow conditions (not simply laminar boundary layers under turbulent flow) remains unclear. Here we address this gap, which is particularly important for mobilization of the smallest particles. We find the Shields parameter to be inversely proportional to the particle Reynolds numbers under laminar jet flow conditions similar to turbulent conditions.

critical shear stresses, erosion rates, effective ↗
DOE OSTI2022

Thixotropic spectra and Ashby-style charts for thixotropy

There is no universal model for thixotropy, and comparing thixotropic effects between different fluids is a subtle yet challenging problem. We introduce a generalized (model-insensitive) framework for comparing thixotropic properties based on thixotropic spectra. A superposition of exponential stress modes distributed over thixotropic time scales is used to quantify buildup and breakdown times and mode strengths in response to step-change input. This mathematical framework is tested with several experimental step-shear rate data on colloidal suspensions. Low-dimensional metrics based on moments of the distribution reveal characteristic average thixotropic properties, which are visualized on Ashby-style diagrams. This method outlines a framework for describing thixotropy across a diverse range of microstructures, supporting scientific studies as well as material selection for engineering design applications.

Mechanics↗
DOE OSTI2024

Direct comparison of using a Z -transformation instead of the traditional $Q^2$ for extraction of the proton radius from $e-p$ scattering data

A discrepancy in the determination of the proton’s charge radius, r p , between muonic hydrogen spectroscopy versus classic atomic spectroscopy and electron scattering data has become known as the proton radius puzzle. Extractions of r p from electron scattering data require determination of the slope of the proton’s charge form factor, $G$$^{p}_{E}$, in the limit of Q 2 → 0 through fitting and extrapolation. Some works have presented the Z-transformation fitting technique as the best choice for this type of extraction due to the true functional form of $G$$^{p}_{E}$ being mathematically guaranteed to exist within the parameter-space of the fit function. In this work, we test this claim by examining the mathematical bias and variances introduced by this technique as compared to the more traditional Q 2 fits using statistically sampled $G$$^{p}_{E}$ parameterizations with known input radii. Furthermore, our tests conclude that the quality of the Z-transformation technique depends on the range of data used. In the case of new experiments, the fit function and technique should be selected in advance by generating realistic pseudodata and assessing the power of different techniques.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Design Principles to Govern Electrode Fabrication for the Lithium Trivanadate Cathode

A full depth of discharge mathematical model for the lithium trivanadate cathode, considering lithiation of the layered α-phase, phase change, and lithiation of the rock-salt like β-phase at lower potentials, is developed. The coupled electrode-scale and crystal-scale model is fit to electrochemical data, and additionally validated with operando EDXRD. There is good agreement between the simulated and measured spatial variation of the volume fraction of the β-phase. This mathematical model is used to guide electrode fabrication, accounting for both ionic and electronic transport effects. Values of three design parameters—electrode thickness, porosity, and volume fraction of conductor—are identified, and the sensitivity of the energy density to these design parameters is quantified. The model is also used to investigate electrode design to create electrodes that deliver the maximum achievable energy density under the constraint that the α to β-phase transition is avoided, since phase change has been demonstrated to reduce cycle life. The energy density sacrificed to avoid phase change decreases at higher discharge rates, but the target values for electrode fabrication remain the same as those when optimizing the electrode for the full depth of discharge.

25 ENERGY STORAGE↗
DOE OSTI2021

An Analysis of the Impact of Particle Growth on Transport Losses in Polymer-Electrolyte Fuel Cells

Voltage cycling causes catalyst nanoparticles in polymer-electrolyte fuel cells to grow. The concomitant loss of interfacial area results in larger kinetic, transport, and possibly ohmic overpotentials. This paper uses recently published experimental data and mathematical models to investigate the evolutions of transport and ohmic resistances to platinum nanoparticles located on the surface and inside the micropores of carbon black supports. Resistance to oxygen transport rises as surface area declines primarily because the flux to each remaining larger platinum particle increases. The path lengths governing oxygen diffusion to surface and buried platinum sites also increase as nanoparticles grow. Platinum nanoparticles on the surface become relatively less favorable as voltage cycling proceeds because they grow faster than platinum in micropores. Because voltage cycling causes total interfacial area to decline, and the fraction located inside micropores to increase, nanoscale ohmic losses increase as a catalyst layer decays. The practical importance of this effect is difficult to discern because proton conductivity in carbon micropores is not well characterized.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Simplified Universal Equations for Ionic Conductivity and Transference Number

Nernst-Einstein equation can provide a reasonable estimate of the ionic conductivity of dilute solutions. For concentrated solutions, alternate methods such as Green–Kubo relations and Einstein relations are more suitable to account for ion-ion interactions. Such computations can be expensive for multicomponent systems. Simplified mathematical expressions like the Nernst-Einstein equation do not exist for concentrated multicomponent mixtures. Newman’s treatment of multicomponent concentrated solutions yields a conductivity relation in terms of species concentration and Onsager phenomenological coefficients. However, the estimation of these phenomenological coefficients is not straightforward. Here, mathematical formulations that relate the phenomenological coefficients with the friction coefficients are developed, leading to simplified, ready-to-use expressions of conductivity and transference numbers that can be used for a wide range of ionic mixtures. This approach involves spectral decomposition of the matrix of Onsager phenomenological coefficients. The general analytical expressions for conductivity and transference number are simplified for binary electrolytes, and numerical solutions are provided for ternary and quaternary mixtures with ion dissociation.

Electrochemistry↗
DOE OSTI2020

CSBCOO SPMM: Compressed Sparse Block Coordinate format sparse matrix dense vector multiply benchmark (CSBCOO SPMM) v1.0

The purpose of this microbenchmark is to provide a means to explore different programming methodologies using a simple, but not trivial, mathematical kernel. The kernel is based on the matrix vector product in the MFDn nuclear configuration interaction code's LOBPCG eigensolver. This particular mathematical operation is a key element in the solution of systems of equations and block eigensolvers.

Cook, Brandon↗
DOE OSTI2022

WORM (Write One, Read Many)

WORM (Write One, Run Many) is an easy to use, cross platform, embedded and extensible, functional programming language designed to facilitate the creation of input-decks for computer codes that use standard ASCII text files for input. WORM makes it easy to create generic (yet, complex and powerful) reusable models. Additionally its nature allows for complex calculations and routines to be coded once and easily reused, further simplifying the creation of input decks. WORM (Write One, Run Many) is a powerful and versatile tool designed to improve the efficiency of today’s criticality safety analyst by allowing: + input decks for parametric studies to be created quickly and easily, + calculations and variables to be imbedded into any input deck, thus allowing for meaningful parameter specifications, + problems to be specified using any combination of units, and + complex mathematically defined models to be created. A very simple syntax is employed, and therefore the WORM is easy to learn. A WORM model is essentially a standard input deck with some of its numerical values replaced by WORM code. WORM code may include and evaluate the following mathematical operators and functions: addition, subtraction, multiplication, division, exponentiation, modulus, sine, cosine, tangent, arcsine, arccosine, arctangent, the natural logarithm, logarithm base 10, integer truncation, absolute value, and random number. Several common constants, e.g., pi, e, and Avogadros’s Number (both as 6.022e23 and 0.6022), are predefined in WORM. Additionally, many unit conversion factors are also predefined: millimeters, meters, inches, feet, yards, and mils to centimeters; kilograms, pounds, and ounces to grams; liters, milliliters, gallons, and fluid ounces to cubic centimeters; and angular degrees to radians. For parametric studies, WORM supports various shorthand list specifications: the explicit step size, linear interpolation, and logarithmic interpolation. The list notation sequentially assigns multiple values to a name. WORM creates an input deck for each value of the name. If multiple lists are used, WORM steps through each list individually, i.e., WORM creates input decks corresponding to each and every permutation of the list values. Additionally, a library of standard material definitions and Perl subroutines are included. Any one of these files can be incorporated into the subject model with a simple WORM read command. WORM is completely written in Perl, the Practical Extraction and Reporting Language. Perl is one of the most portable programming languages available today. As such, the WORM works on practically any computer platform.

Sartor, Raymond↗
DOE OSTI2023

EI_MS_ML

The unambiguous identification of compounds from their electron ionization mass (EI-MS) spectra remains a significant unsolved problem in the field of metabolomics and analytical chemistry as a whole. Typically EI-MS spectra are compared using various mathematical operations that convert the spectral similarity or differences into a distance-like metric that roughly approximates the similarity of any two spectra. A commonly used metric for this is the cosine similarity metric which has values close to one for very similar spectra and a value of zero for very dissimilar spectra; however, no metric is perfect. Due to the prevalence of structurally-similar compounds such as isomers and the prevalence of certain fragmentation patterns across structurally-dissimilar compounds, the unambiguous assignment of EI-MS spectra compounds remains difficult. Frequently, querying an observed EI-MS spectrum against a large database such as the NIST17 library yields multiple possible assignments requiring the end user to distinguish between multiple high scoring hits, or multiple low scoring hits while keeping in mind that the correct hit may not be in the database at all. Although techniques such as orthogonal information from techniques such as chromatography can greatly aid in unambiguous assignment, this also requires more complicated experimental designs and access to more complicated analytical instrumentation. Substructures can be trivially detected and represented as strings using a previously published technique called node coloring from a known chemical structure. However, for experimentally-derived EI-MS spectra this information must be derived from the spectra itself (i.e., because we do not know what compound it represents). To achieve this, the software uses techniques from the field of machine learning and a large training dataset of EI-MS spectra corresponding to known structures annotated with substructure strings, to build models that can predict the presence of a given chemical substructure from an EI-MS spectrum directly.If these predictions are of high-quality (i.e., are unlikely to be false positives), the presence of one or more predicted substructures can be used to constrain the number of possible hits for a query spectrum. Mathematically, this restriction could be expressed in many forms, but the most straight-forward implementation is to weight the cosine similarity of a query spectrum and a plausible database match with a Tanimoto-like coefficient based on the ratio of the number of substructures predicted to the number of substructures present in the potential database hit. Determining which combination of models best reduces assignment ambiguity will be achieved using a combination of manual curation and optimization techniques such as genetic algorithms. This software will perform all the steps necessary to construct said models from a training dataset and evaluate them using a holdout dataset. Various statistical analyses can be performed to determine if this approach does decrease assignment ambiguity. For example, if this approach works, on average, the rank-order of the correct assignment for the holdout set of EI-MS spectra should decrease and the weighted cosine similarities for most of the possible matches in the database should be better than the unweighted cosine similarities. Furthermore, this same pipeline can be used on real experimental data to generate less ambiguous assignments.

Mitchell, Joshua↗
DOE OSTI2023

A python package for analyzing Resilience of Complex Systems (pyRoCS) v.0.0

SAND2024-01040O PyRoCS software synthesizes mathematical equations from several domains—including information theory, ecology, and engineering sciences—to support resilience analysis for complex systems. Resilience is the ability of the complex system being analyzed to withstand, operate through, and recover from a disruption. The complex system can be a physical system such as an electric grid, an organization such as a company, or even a subfunction of an organization. Existing mathematical equations for resilience analysis are found within multiple domains including information theory, biological sciences, and complex systems. This package synthesizes and refactors equations from these various domains to make them more generalizable for application across different types of complex systems relevant for resilience analysis. Users will be able to apply these equations to characterize different components of complex systems based on available data. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Verzi, Stephen↗
DOE OSTI2023

Boost Converter Loss Analysis v0.1

This software helps the user mathematically determine the losses in a boost converter they are designing. They can enter the component parasitic values from their component datasheets, and the software will calculate exactly where the losses are coming from. This can help the designer target and reduce specific sources of loss. It also helps system engineers compare the losses between AC and DC versions of the boost converter. Finally, it establishes a set of mathematical equations to help academics study loss theoretically.

Gerber, Daniel↗
DOE OSTI2026

Fortran mimetic abstraction language (Formal) v0.1.

The Fortran mimetic abstraction language ("Formal") is a domain-specific language (DSL) embedded in Fortran 202Y [1]. Formal provides novel software abstractions for simulating phenomena governed by the partial differential equations (PDEs) of vector and tensor calculus. Such equations model an extremely broad set of physical phenomena, ranging from atmospheric winds to light propagation. Formal's data structures and algorithms mimic in form and behavior continuous functions and operators. Formal supports these mathematical constructs using mimetic discretizations that define a discrete calculus satisfying various tensor calculus theorems, thereby ensuring high-fidelity representations of the physics being modeled. [2] Formal 0.1.0 also lays a foundation for the future use of Fortran 202Y type-safe templates to facilitate the formal verification of tensor contractions in computational physics and artificial intelligence [3]. [1] "Fortran 202Y" is Fortran standard committee's informal designation for the next Fortran revision, which will likely be "Fortran 2028". [2] Corbino, J. and Castillo, J. (2020) Journal of Computational and Applied Mathematics, https://doi.org/10.1016/j.cam.2019.06.042. [3] Haveraaen, M., Järvi, J., & Rouson, D. (2019). Reflecting on Generics for Fortran. https://j3-fortran.org/doc/year/19/19-188.pdf.

Rouson, Damian [Lawrence Berkeley National Laborat↗
DOE OSTI2021

How Snow Aggregate Ellipsoid Shape and Orientation Variability Affects Fall Speed and Self-Aggregation Rates

Snow aggregate shapes and orientations have long been known to exhibit substantial variability. Despite this observed variability, most weather and climate prediction models use fixed power-law functions that deterministically map particle size to mass and fall speed. As such, integrated quantities like precipitation and self-aggregation rates currently ignore nonlinear effects resulting from variation in shape and orientation for aggregates of the same size. This study therefore develops an analytic framework that couples an empirically based bivariate distribution of ellipsoid shapes to classical hydrodynamic theory so as to capture an appropriate dispersion of masses, projected areas, and fall speeds for an assumed size distribution. For a fixed aggregate size, shape variations produce approximately ±0.13 m s -1 standard deviation of fall speed which increases the mass flux fall speed dispersion by more than 100% over traditional microphysics models. This increased fall speed dispersion results predominantly from shape-induced mass dispersion whereas orientation and drag dispersion play a lesser role. Shape variations can increase mass- and reflectivity-weighted fall speeds by up to 60% of traditional models whereas self-aggregation rates can increase by a factor of 100 for very small slope parameters. This implies that aggregate shape variations effectively forestall the theorized onset of fall speed distribution narrowing and subsequent quenching of the aggregation process. As a result, it is likely that secondary ice formation is necessary to prevent an ever decreasing slope parameter. The mathematical theory presented in this study is used to develop simple correction factors for snow forecast and climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Ginkgo - A math library designed to accelerate Exascale Computing Project science applications

Large-scale simulations require efficient computation across the entire computing hierarchy. A challenge of the Exascale Computing Project (ECP) was to reconcile highly heterogeneous hardware with the myriad of applications that were required to run on these supercomputers. Mathematical software forms the backbone of almost all scientific applications, providing efficient abstractions and operations that are crucial to harness the performance of computing systems. Ginkgo is one such mathematical software library, nurtured by ECP, providing high-performance, user-friendly, and performance portable interfaces for applications in ECP and beyond. In this paper, we elaborate on Ginkgo’s philosophy of high-performance software that is sustainable, reproducible, and easy to use. We showcase the wide feature set of solvers and preconditioners available in Ginkgo and the central concepts involved in their design. We elaborate on four different ECP software integrations: MFEM, PeleLM + SUNDIALS, XGC, and ExaSGD that use Ginkgo to accelerate their science runs. Performance studies of different problems from these applications highlight the effectiveness of Ginkgo and the benefits incurred by these ECP applications.

Cojean, Terry↗
DOE OSTI2021

Stress analysis using direct-S wavelets produced by a vertical vibrator

We have determined how to extract the azimuth of maximum horizontal stress (SHmax) in deep rocks by doing a simple 360° mathematical rotation of a downgoing direct-S wavelet generated at the baseplate of a surface-based vertical vibrator. We worked with direct-S wavelets that travel through stressed rocks to a deep horizontal vertical seismic profile (VSP) geophone. We find that the azimuth where a polarity reversal occurs in mathematical rotations of this downgoing direct-S wavelet defines the azimuth of the SHmax in the overlying rocks. We tested this direct-S wavelet rotation method for determining the SHmax azimuth at a site in the Illinois Basin using legacy VSP data acquired in 2013. SHmax azimuths indicated by this simple wavelet-rotation method were determined when vertical vibrators were stationed at zero offset, at far offset, and at different azimuths around a VSP receiver well. These VSP-based SHmax azimuths agreed with the azimuth of the SHmax found by traditional minifrac tests in the VSP receiver well. This simple VSP data analysis procedure for detecting the azimuth of maximum horizontal stress has never, to our knowledge, been reported or discussed in the geophysical literature. This technical finding should be of interest to the worldwide geophysical community, especially to people who need to monitor how stress fields shift when fluids are injected into, or extracted from, deep porous reservoirs.

Geochemistry & Geophysics↗
DOE OSTI2022

Analytic error analysis of cross section interpolation methods in nodal diffusion codes - I : Theory

This paper discusses two cross section interpolation methods commonly found in popular nodal codes; the partial derivatives and multiple tables models. The motivation for choosing a model, and thus a case matrix structure, is a trade off between accuracy and computational cost. Due to decades of experience, there are default structures that are sufficient for current light water reactor analysis. However, this is not necessarily the case for advanced reactor designs. Therefore, it is advantageous to understand the sources of error in cross section interpolation models so that the quality of a case matrix may be improved. A mathematical framework for these models is presented in this work that provides a more rigorous connection between the nuclear engineering field's cross section interpolation methods and the broader mathematical field of function approximation. The two cross section models examined in this paper were found to utilize Lagrange interpolation and are a subset of Lagrange tensor products. Classical results of Lagrange polynomial error analysis were then applied to the partial derivative and multiple tables models to derive expressions for the total point-wise error. The analytical results classify the total error into two parts: the model form error and interpolation error. Finally, based on our observations, a better foundation for improving the quality of a case matrix is proposed. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

A quantitative model used to compare within-host SARS-CoV-2, MERS-CoV, and SARS-CoV dynamics provides insights into the pathogenesis and treatment of SARS-CoV-2

The scientific community is focused on developing antiviral therapies to mitigate the impacts of the ongoing novel coronavirus disease 2019 (COVID-19) outbreak. This will be facilitated by improved understanding of viral dynamics within infected hosts. Here, using a mathematical model in combination with published viral load data, we compare within-host viral dynamics of SARS-CoV-2 with analogous dynamics of MERS-CoV and SARS-CoV. Our quantitative analyses using a mathematical model revealed that the within-host reproduction number at symptom onset of SARS-CoV-2 was statistically significantly larger than that of MERS-CoV and similar to that of SARS-CoV. In addition, the time from symptom onset to the viral load peak for SARS-CoV-2 infection was shorter than those of MERS-CoV and SARS-CoV. These findings suggest the difficulty of controlling SARS-CoV-2 infection by antivirals. We further used the viral dynamics model to predict the efficacy of potential antiviral drugs that have different modes of action. The efficacy was measured by the reduction in the viral load area under the curve (AUC). Our results indicate that therapies that block de novo infection or virus production are likely to be effective if and only if initiated before the viral load peak (which appears 2–3 days after symptom onset), but therapies that promote cytotoxicity of infected cells are likely to have effects with less sensitivity to the timing of treatment initiation. Furthermore, combining a therapy that promotes cytotoxicity and one that blocks de novo infection or virus production synergistically reduces the AUC with early treatment. Our unique modeling approach provides insights into the pathogenesis of SARS-CoV-2 and may be useful for development of antiviral therapies.

59 BASIC BIOLOGICAL SCIENCES↗