Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mathematical physics”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

JCZ3: An Equation of State Based on the Exponential-Six Intermolecular Potential – Update: Formulation

This report presents a comprehensive analysis of the JCZ3 equation of state (EOS) as implemented in the TIGER code, focusing on its thermodynamic framework, mathematical formulations, and implications for modeling gas phase behavior under varying conditions. Overall, the report affirms the TIGER code's foundational robustness and its potential as a valuable tool for simulating high-pressure, high-temperature gas mixtures, while also highlighting areas for further refinement and validation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2004

Macromolecular Crystallization in Microgravity

The key concepts that attracted crystal growers, macromolecular or solid state, to microgravity research is that density difference fluid flows and sedimentation of the growing crystals are greatly reduced. Thus, defects and flaws in the crystals can be reduced, even eliminated, and crystal volume can be increased. Macromolecular crystallography differs from the field of crystalline semiconductors. For the latter, crystals are harnessed for their electrical behaviors. A crystal of a biological macromolecule is used instead for diffraction experiments (X-ray or neutron) to determine the three-dimensional structure of the macromolecule. The better the internal order of the crystal of a biological macromolecule then the more molecular structure detail that can be extracted. This structural information that enables an understanding of how the molecule functions. This knowledge is changing the biological and chemical sciences with major potential in understanding disease pathologies. Macromolecular structural crystallography in general is a remarkable field where physics, biology, chemistry, and mathematics meet to enable insight to the basic fundamentals of life. In this review, we examine the use of microgravity as an environment to grow macromolecular crystals. We describe the crystallization procedures used on the ground, how the resulting crystals are studied and the knowledge obtained from those crystals. We address the features desired in an ordered crystal and the techniques used to evaluate those features in detail. We then introduce the microgravity environment, the techniques to access that environment, and the theory and evidence behind the use of microgravity for crystallization experiments. We describe how ground-based laboratory techniques have been adapted to microgravity flights and look at some of the methods used to analyze the resulting data. Several case studies illustrate the physical crystal quality improvements and the macromolecular structural advances. Finally, limitations and alternatives to microgravity and future directions for this research are covered.

Snell, Edward H.↗
DOE OSTI2026

Discovering Strong Gravitational Lenses in the Dark Energy Survey with Interactive Machine Learning and Crowd-sourced Inspection with Space Warps

We conduct a search for strong gravitational lenses in the Dark Energy Survey (DES) Year 6 imaging data. We implement a pre-trained Vision Transformer (ViT) for our machine learning (ML) architecture and adopt interactive machine learning to construct a training sample with multiple classes to address common types of false positives. Our ML model reduces ∼236 million DES cutout images to 22,564 targets of interest, including ∼85% of previously reported galaxy–galaxy lens candidates discovered in DES. These targets were visually inspected by citizen scientists, who ruled out ∼90% as false positives. Of the remaining 2618 candidates, 149 were expert-classified as “definite” lenses and 516 as “probable” lenses, for a total of 665 systems, with 147 of these candidates being newly identified. Additionally, we trained a second ViT to find double-source plane lens systems, finding at least one double-source system. Our main ViT excels at identifying galaxy–galaxy lenses, consistently assigning high scores to candidates with high expert assessments. The top 800 ViT-scored images include ∼100 of our “definite” lens candidates. This selection is an order of magnitude higher in purity than previous convolutional neural-network-based lens searches and demonstrates the feasibility of applying our methodology for discovering large samples of lenses in future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Accelerometer Modeling in the State-Space

Generically, accelerometers are sensors that measure the accelerating force applied to an object. Their operation is based on a simple mechanical system and D’Alembert’s force principle along with Newton’s laws of motion. This device is typically mounted to the body of the test object under investigation that is subjected to an accelerating force. The displacement of the seismic or “proof” mass due to the acceleration is sensed by an embedded component from which the acceleration is determined. Interestingly enough, the resulting the resulting mathematical relations governing the dynamical motion is based on the underlying physics, sensor materials, geometry and fabrication. Mathematically, the various designs satisfy the underlying Newtonian equations of motion, but their type and performance differ significantly.

42 ENGINEERING↗
DOE OSTI2026

An Introduction to Variational Calculus

Variational Calculus is an advanced topic in mathematics. This document exists to compile the basics of Variational Calculus in a readable and digestible manner for application by engineers. This document is a living document and will be updated as time moves forward. It is the hope that by reading this document, the basics of Variational Calculus can be readily applied to any system/function of interest. This does not serve as a replacement for textbooks or other learned sources but should act as a companion piece to better serve the reader in times of confusion and distress.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Bayesian Monte Carlo Evaluation Framework for Imperfect Data and Models [Abstract]

Nuclear data evaluation methods conventionally make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate for minimization of a cost function (even for non-linear models); and that both the model (of, e.g., neutron cross section) and experimental data (including their covariance data) are perfect. These assumptions are inherent to the well-known generalized linear least squares (GLLS) minimization method commonly used for evaluations of resolved resonance region (RRR) neutron cross sections. However, these assumptions are often not justified due to the presence of non-normal PDFs, non-linear models (e.g. R -matrix formalism), and inherent imperfections in data and models (e.g. discrepant data sets, discrepancies between the previous evaluation and newly measured data, or imperfect covariance data). We remove the said assumptions in a mathematical framework of Bayes’ theorem, and implement it using the Metropolis-Hastings Monte Carlo method. Parameters of a new kind are introduced to parameterize inherent imperfections, e.g. , any discrepancies between the theoretical model and measured data. These new parameters enable evaluators to quantify their expert judgement about any discrepancies or imperfections in a reproducible manner. We demonstrate the framework with an ongoing evaluation of 233 U in the eV region using the ENDF-B/VIII library and transmission data measured by Guber, et al. , and compare the posterior parameters to those obtained by conventional evaluation methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1986

Physics through the 1990s: Elementary-particle physics

The volume begins with a non-mathematical discussion of the motivation behind, and basic ideas of, elementary-particle physics theory and experiment. The progress over the past two decades with the quark model and unification of the electromagnetic and weak interactions is reviewed. Existing theoretical problems in the field, such as the origin of mass and the unification of the fundamental forces, are detailed, along with experimental programs to test the new theories. Accelerators, instrumentation, and detectors are described for both current and future facilities. Interactions with other areas of both theoretical and applied physics are presented. The sociology of the field is examined regarding the education of graduate students, the organization necessary in large-scale experiments, and the decision-making process involved in high-cost experiments. Finally, conclusions and recommendations for maintaining US excellence in theory and experiment are given. Appendices list both current and planned accelerators, and present statistical data on the US elementary-particle physics program. A glossary is included.

Source record↗
DOE OSTI2025

Diffusion Codes: Self-Correction from Small(er)-Set Expansion with Tunable Non-locality

Optimal constructions of classical LDPC codes can be obtained by choosing the Tanner graph uniformly at random among biregular graphs. We introduce a class of codes that we call ``diffusion codes'', defined by placing each edge connecting bits and checks on some graph, and acting on that graph with a random SWAP network. By tuning the depth of the SWAP network, we can tune a tradeoff between the amount of randomness -- and hence the optimality of code parameters -- and locality with respect to the underlying graph. For diffusion codes defined on the cycle graph, if the SWAP network has depth $\sim Tn$ with $T> n^{2β}$ for arbitrary $β>0$, then we prove that almost surely the Tanner graph is a lossless ``smaller set'' vertex expander for small sets up size $δ\sim \sqrt T \sim n^β$, with bounded bit and check degree. At the same time, the geometric size of the largest stabilizer is bounded by $\sqrt T$ in graph distance. We argue, based on physical intuition, that this result should hold more generally on arbitrary graphs. By taking hypergraph products of these classical codes we obtain quantum LDPC codes defined on the torus with smaller-set boundary and co-boundary expansion and the same expansion/locality tradeoffs as for the classical codes. These codes are self-correcting and admit single-shot decoding, while having the geometric size of the stabilizer growing as an arbitrarily small power law. Our proof technique establishes mixing of a random SWAP network on small subsystems at times scaling with only the subsystem size, which may be of independent interest.

Combinatorics (math.CO)↗
DOE OSTI2023

Bayesian Monte Carlo Evaluation of Imperfect (n, 233 U) Data and Model

Conventional nuclear data evaluation methods using generalized linear least squares make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate to minimize the cost function (even for nonlinear models); the model (e.g., of neutron cross section) and experimental data (including covariance data) are without defect and prior PDFs of parameters and measured data are known perfectly. Neglect of covariance between model parameters and measured data in conventional evaluations contributes to imperfections. These assumptions are inherent to the generalized linear least squares minimization method commonly used for resolved resonance region neutron cross section evaluations but are often not justified due to the presence of non-normal PDFs, nonlinear models (e.g., R-matrix formalism), and inherent imperfections in data and models (e.g., imperfect covariance data). Here, these assumptions are removed in a mathematical framework of Bayes’ theorem, which is implemented using the Metropolis-Hastings Monte Carlo method. Most importantly, new parameters are introduced to parameterize discrepancies between the theoretical model and measured data to quantify judgement about discrepancies or imperfections in a reproducible manner. An evaluation of 233U in the eV region using the ENDF-B/VIII.0 library and transmission data (Guber et al.) is presented, and posterior parameters are compared to those obtained by conventional evaluation methods. This example illustrates the effects of removing the most harmful assumption: that of model-data perfection.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Molecular NMR shieldings, J -couplings, and magnetizabilities from numeric atom-centered orbital based density-functional calculations

This paper reports and benchmarks a new implementation of nuclear magnetic resonance shieldings, magnetizabilities, and J-couplings for molecules within semilocal density functional theory, based on numeric atom-centered orbital (NAO) basis sets. NAO basis sets are attractive for the calculation of these nuclear magnetic resonance (NMR) parameters because NAOs provide accurate atomic orbital representations especially near the nucleus, enabling high-quality results at modest computational cost. Moreover, NAOs are readily adaptable for linear scaling methods, enabling efficient calculations of large systems. Here, the paper has five main parts: (1) It reviews the formalism of density functional calculations of NMR parameters in one comprehensive text to make the mathematical background available in a self-contained way. (2) The paper quantifies the attainable precision of NAO basis sets for shieldings in comparison to specialized Gaussian basis sets, showing similar performance for similar basis set size. (3) The paper quantifies the precision of calculated magnetizabilities, where the NAO basis sets appear to outperform several established Gaussian basis sets of similar size. (4) The paper quantifies the precision of computed J-couplings, for which a group of customized NAO basis sets achieves precision of ~Hz for smaller basis set sizes than some established Gaussian basis sets. (5) The paper demonstrates that the implementation is applicable to systems beyond 1000 atoms in size.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

A Perspective on Sustainable Computational Chemistry Software Development and Integration

The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.

36 MATERIALS SCIENCE↗
DOE OSTI2021

2019 Computing Sciences Strategic Plan

Computing has transformed nearly every aspect of scientific inquiry — across disciplines and across scales — from the behavior of subatomic particles to the formation of structures in the early universe, from the assembly of the human genome to the evolution of earth systems. Over the past two decades, computing has become an integral part of how Berkeley Lab is “Bringing Science Solutions to the World.” Advances in computing and mathematics have been key, with new mathematical models of complex physical phenomena, new methods for analyzing complex data, new algorithms for accuracy and scaling and sophisticated software systems that encapsulate these techniques into open, reusable tools. The performance of NERSC computers and the ESnet network have grown by several orders of magnitude, along with our understanding of how to map scientific computations and workflows onto these systems. From research to facility operations, the passion, talent and dedication of the Computing Sciences Area staff has been the cornerstone of our success. The plan outlined in this document describes the next step in a journey to expand the influence and impact of our efforts, building an increasingly connected global enterprise for science that places more powerful instruments in the hands of scientists, along with more powerful methods and tools for modeling, analysis and prediction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Multilevel Monte Carlo methods for the Grad-Shafranov free boundary problem

The equilibrium configuration of a plasma in an axially symmetric reactor is described mathematically by a free boundary problem associated with the celebrated Grad-Shafranov equation. The presence of uncertainty in the model parameters introduces the need to quantify the variability in the predictions. This is often done by computing a large number of model solutions on a computational grid for an ensemble of parameter values and then obtaining estimates for the statistical properties of solutions. In this study, we explore the savings that can be obtained using multilevel Monte Carlo methods, which reduce costs by performing the bulk of the computations on a sequence of spatial grids that are coarser than the one that would typically be used for a simple Monte Carlo simulation. We examine this approach using both a set of uniformly refined grids and a set of adaptively refined grids guided by a discrete error estimator. Numerical experiments show that multilevel methods dramatically reduce the cost of simulation, with cost reductions typically on the order of 60 or more and possibly as large as 200. Furthermore, adaptive griding results in more accurate computation of geometric quantities such as x-points associated with the model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1997

Radiative heating of interstellar grains falling toward the solar nebula: 1-D diffusion calculations

As the dense molecular cloud that was the precursor of our Solar System was collapsing to form a protosun and the surrounding solar-nebula accretion disk, infalling interstellar grains were heated much more effectively by radiation from the forming protosun than by radiation from the disk's accretion shock. Accordingly, we have estimated the temperatures experienced by these infalling grains using radiative diffusion calculations whose sole energy source is radiation from the protosun. Although the calculations are 1-dimensional, they make use of 2-D, cylindrically symmetric models of the density structure of a collapsing, rotating cloud. The temperature calculations also utilize recent models for the composition and radiative properties of interstellar grains (Pollack et al. 1994. Astrophys. J. 421, 615-639), thereby allowing us to estimate which grain species might have survived, intact, to the disk accretion shock and what accretion rates and molecular-cloud rotation rates aid that survival. Not surprisingly, we find that the large uncertainties in the free parameter values allow a wide range of grain-survival results: (1) For physically plausible high accretion rates or low rotation rates (which produce small accretion disks), all of the infalling grain species, even the refractory silicates and iron, will vaporize in the protosun's radiation field before reaching the disk accretion shock. (2) For equally plausible low accretion rates or high rotation rates (which produce large accretion disks), all non-ice species, even volatile organics, will survive intact to the disk accretion shock. These grain-survival conclusions are subject to several limitations which need to be addressed by future, more sophisticated radiative-transfer models. Nevertheless, our results can serve as useful inputs to models of the processing that interstellar grains undergo at the solar nebula's accretion shock, and thus help address the broader question of interstellar inheritance in the solar nebula and present Solar System. These results may also help constrain the size of the accretion disk: for example, if we require that the calculations produce partial survival of organic grains into the solar nebula, we infer that some material entered the disk intact at distances comparable to or greater than a few AU. Intriguingly, this is comparable to the heliocentric distance that separates the C-rich outer parts of the current Solar System from the C-poor inner regions.

NASA Center ARC↗
DOE OSTI2025

Inference of phase field fracture models

The phase field approach to modeling fracture uses a diffuse damage field to represent cracks. This representation mollifies singularities that arise in computations with sharp interface models and some of the resultant difficulties in the mathematical and numerical treatment of fracture. Phase field fracture models have proven effective at representing crack propagation, branching, and merging. Specific formulations, beginning with brittle fracture, have also been shown to converge to classical solutions. Extensions to cover the range of material failure, including ductile and cohesive fracture, lead to an array of possible models. There exists a large body of literature focusing on this class of models and on the impact of model form on the predicted crack evolution. However, there have not been systematic studies into how optimal models may be chosen. Here, we take a first step in this direction by developing formal methods for identification of the best parsimonious model of phase field fracture given full-field data on the damage and deformation fields. We consider some of the main models that have been used for the degradation of elastic response due to damage and its propagation. Our approach builds upon Variational System Identification (VSI), a weak form variant of the Sparse Identification of Nonlinear Dynamics (SINDy). Furthermore, in this first communication we focus on synthetically generated data but we also consider central issues associated with the use of experimental full-field data, such as data sparsity and noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Hybrid large-eddy simulation with adaptive downstream anisotropic eddy viscosity model

Turbulence pervades most flows of engineering interest, and its prediction remains a challenge on both accuracy and cost. One promising predictive approach that reduces cost combines large-eddy simulation (LES) with simulation based on Reynolds averaged Navier–Stokes equations (RANS). This study presents a method to overcome stability and accuracy issues associated with these hybrid LES-RANS methods. The method developed involves extracting the Reynolds stresses from the upstream LES solution and then using that information to convert the downstream RANS model from a scalar eddy-viscosity model to an anisotropic nonlinear eddy-viscosity model. The method developed differs from the downstream tensorial eddy-viscosity model by being independent of the coordinate system. The method developed was evaluated by computing film cooling of a flat plate with the coolant injected through one row of circular holes. Results obtained show instabilities at the LES-to-RANS interface to be eliminated. Also, the method developed yielded solutions that compare reasonably well with those from LES, even though a significant portion of the flow is computed by the adapted anisotropic RANS model instead of LES, which significantly reduced the number of grid points and computational time needed. Since modification of the downstream RANS model is based on information extracted from the upstream LES solution, the method developed is adaptive to the problem being studied.

42 ENGINEERING↗