Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational physics applications”

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 OSTI2021

GeV-scale neutrinos: interactions with mesons and DUNE sensitivity

The simplest extension of the SM to account for the observed neutrino masses and mixings is the addition of at least two singlet fermions (or right-handed neutrinos). If their masses lie at or below the GeV scale, such new fermions would be produced in meson decays. Similarly, provided they are sufficiently heavy, their decay channels may involve mesons in the final state. Although the couplings between mesons and heavy neutrinos have been computed previously, significant discrepancies can be found in the literature. The aim of this paper is to clarify such discrepancies and provide consistent expressions for all relevant effective operators involving mesons with masses up to 2 GeV. Moreover, the effective Lagrangians obtained for both the Dirac and Majorana scenarios are made publicly available as FeynRules models so that fully differential event distributions can be easily simulated. As an application of our setup, we numerically compute the expected sensitivity of the DUNE near detector to these heavy neutral leptons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Quantifying Volume Change in Porous Electrodes via the Multi-Species, Multi-Reaction Model

Automotive manufacturers are working to improve individual cell and overall pack design by increasing their performance, durability, and range, while reducing cost; and active material volume change is one of the more complex aspects that needs to be considered during this process. As the time from initial design to manufacture of electric vehicles is decreased, design work that used to rely solely on testing needs to be supplemented or replaced by virtual methods. As electrochemical engineers drive battery and system design using model-based methods, the need for coupled electrochemical/mechanical models that take into account the active material change utilizing physics based or semi-empirical approaches is necessary. In this study, we illustrated the applicability of a mechano-electrochemical coupled modeling method considering the multi-species, multi-reaction model as popularized by Verbrugge and Baker. To do this, validation tests were conducted using a computer-controlled press apparatus that can control the press displacement and press force with precision. The coupled MSMR volume change model was developed and its applicability to graphite and NMC cells was illustrated. The increased accuracy of the model considering the coupled MSMR volume change approach shows in the importance of accounting for individual gallery volume change behavior on cell level predictions.

25 ENERGY STORAGE↗
DOE OSTI2021

CARDINAL

Cardinal is a code that automates the workflow for iteratively solving neutral particle transport via the OpenMC code and computational fluid dynamics via the NekRS code within the Multiphysics Object Oriented Simulation Environment (MOOSE), a finite element framework for solving general partial differential equations. Cardinal establishes a mapping between the Monte Carlo geometry representation used in OpenMC and the spectral element mesh used in NekRS to other physics applications built upon the MOOSE framework, allowing high-resolution particle transport and fluid dynamics physics feedback to the MOOSE application ecosystem. Heat generation rates as solved by OpenMC are used as heat sources in MOOSE, while the temperature and density fields from NekRS are used as conjugate heat transfer boundary conditions and/or source fields in MOOSE. OpenMC, Nek5000, and MOOSE are open source, community-developed simulation codes.

NOVAK, APRILJEAN RODGERS↗
DOE OSTI2021

Coupling of regional geophysics and local soil-structure models in the EQSIM fault-to-structure earthquake simulation framework

Accurate understanding and quantification of the risk to critical infrastructure posed by future large earthquakes continues to be a very challenging problem. Earthquake phenomena are quite complex and traditional approaches to predicting ground motions for future earthquake events have historically been empirically based whereby measured ground motion data from historical earthquakes are homogenized into a common data set and the ground motions for future postulated earthquakes are probabilistically derived based on the historical observations. This procedure has recognized significant limitations, principally due to the fact that earthquake ground motions tend to be dictated by the particular earthquake fault rupture and geologic conditions at a given site and are thus very site-specific. Historical earthquakes recorded at different locations are often only marginally representative. There has been strong and increasing interest in utilizing large-scale, physics-based regional simulations to advance the ability to accurately predict ground motions and associated infrastructure response. However, the computational requirements for simulations at frequencies of engineering interest have proven a major barrier to employing regional scale simulations. In a U.S. Department of Energy Exascale Computing Initiative project, the EQSIM application development is underway to create a framework for fault-to-structure simulations. This framework is being prepared to exploit emerging exascale platforms in order to overcome computational limitations. This article presents the essential methodology and computational workflow employed in EQSIM to couple regional-scale geophysics models with local soil-structure models to achieve a fully integrated, complete fault-to-structure simulation framework. Here, the computational workflow, accuracy and performance of the coupling methodology are illustrated through example fault-to-structure simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Study of recirculating liquid fuel in a 1D critical stationary system

Several studies have been conducted to investigate the physics of liquid fuel reactors, showing also applications with molten salts. A liquid nuclear fuel implies changes in the neutron balance equation to take into account the precursors' displacement and the emission of delayed neutrons in a different position than at the original fission. This requires to upgrade the computer codes normally used to calculate nuclear reactors using only solid fuel. In this work, we revisit a simple problem with liquid fuel, which is proposed for the verification of the numerical solutions obtained by advanced computer codes. This problem studies criticality with constant coefficients, thus neglecting thermal feedback. We elaborate on the analytical solution of the problem, deriving also a generalized eigenvalue problem by finite-volume integration over the cells of a discretized mesh to study the evolution of the dominance ratio with fuel velocity. Finally, we investigate the influence of the fuel velocity on the reactivity of the system. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Geant4 Event Biasing and Fast Simulation

Geant4 offers advanced event biasing techniques to significantly accelerate simulations involving rare events. Various biasing methods, such as leading particle selection, cross-section biasing, radioactive decay enhancement, and bremsstrahlung splitting, enable efficient event sampling, though they require careful handling. Additionally, Geant4 provides a Fast Simulation Interface, allowing the replacement of standard processes in specific region and for selected particles, enabling faster execution or external code integration. Applications of fast simulation include electromagnetic shower modeling in calorimeters, machine learning inference, and offloading tasks to specialized hardware like GPUs, making Geant4 a powerful tool for computationally demanding simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Snowmass 2021 Whitepaper: Proton Structure at the Precision Frontier

An overwhelming number of theoretical predictions for hadron colliders require parton distribution functions (PDFs), which are an important ingredient of theory infrastructure for the next generation of high-energy experiments. This whitepaper summarizes the status and future prospects for determination of high-precision PDFs applicable in a wide range of energies and experiments, in particular in precision tests of the Standard Model and in new physics searches at the high-luminosity Large Hadron Collider and Electron–Ion Collider. We discuss the envisioned advancements in experimental measurements, QCD theory, global analysis methodology, and computing that are necessary to bring unpolarized PDFs in the nucleon to the N2LO and N3LO accuracy in the QCD coupling strength. Special attention is given to the new tasks that emerge in the era of the precision PDF analysis, such as those focusing on the robust control of systematic factors both in experimental measurements and theoretical computations. Various synergies between experimental and theoretical studies of the hadron structure are explored, including opportunities for studying PDFs for nuclear and meson targets, PDFs with electroweak contributions or dependence on the transverse momentum, for incisive comparisons between phenomenological models for the PDFs and computations on discrete lattice, and for cross-fertilization with machine learning/AI approaches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1972

Generalized simulation technique for turbojet engine system analysis

A nonlinear analog simulation of a turbojet engine was developed. The purpose of the study was to establish simulation techniques applicable to propulsion system dynamics and controls research. A schematic model was derived from a physical description of a J85-13 turbojet engine. Basic conservation equations were applied to each component along with their individual performance characteristics to derive a mathematical representation. The simulation was mechanized on an analog computer. The simulation was verified in both steady-state and dynamic modes by comparing analytical results with experimental data obtained from tests performed at the Lewis Research Center with a J85-13 engine. In addition, comparison was also made with performance data obtained from the engine manufacturer. The comparisons established the validity of the simulation technique.

Seldner, K.↗
NASA NTRS1990

Cumulative reports and publications through December 31, 1989

A complete list of reports from the Institute for Computer Applications in Science and Engineering (ICASE) is presented. The major categories of the current ICASE research program are: numerical methods, with particular emphasis on the development and analysis of basic numerical algorithms; control and parameter identification problems, with emphasis on effectual numerical methods; computational problems in engineering and the physical sciences, particularly fluid dynamics, acoustics, structural analysis, and chemistry; computer systems and software, especially vector and parallel computers, microcomputers, and data management. Since ICASE reports are intended to be preprints of articles that will appear in journals or conference proceedings, the published reference is included when it is available.

Source record↗
NASA NTRS1991

Computational fluid dynamics for hypersonic airbreathing aircraft

A development status evaluation is presented for CFD methods applicable to fuselage-integrated scramjet powerplant incorporating hypersonic vehicles; these methods are critically important due to the unavailability of experimental facilities for such elevated Mach number/high-enthalphy conditions. Advancements are required in algorithm robustness and speed, geometric flexibility, and the inclusion of more complete flow physics. The most serious deficiencies lie in turbulence modeling, the lack of complete transition-prediction methods, and combustion modeling.

Thomas, James L.↗
NASA NTRS1993

A three-dimensional pressure flux-split RNS application to sub/supersonic flow in inlets and ducts

The reduced Navier-Stokes (RNS) formulation is combined with a pressure based flux split procedure for the computational analysis of three dimensional flow. Only lowest order diffusion, that is required to satisfy the no-slip boundary condition is retained in this approximation. The governing equations are appropriately differenced such that the physical boundary conditions are closely coupled and combined to form a closed discrete system. No numerical or characteristic type boundary conditions are required. A global pressure relaxation procedure, i.e. multi-sweep PNS, is considered. A sparse matrix direct solver is applied to the crossplane solution. At the outflow boundary only the pressure or pressure gradient is prescribed. The resulting formulation has been applied to the solution of three-dimensional subsonic/supersonic internal flows in inlets and ducts.

Khosla, P. K.↗
NASA NTRS2011

Circuit for Communication Over Power Lines

Many distributed systems share common sensors and instruments along with a common power line supplying current to the system. A communication technique and circuit has been developed that allows for the simple inclusion of an instrument, sensor, or actuator node within any system containing a common power bus. Wherever power is available, a node can be added, which can then draw power for itself, its associated sensors, and actuators from the power bus all while communicating with other nodes on the power bus. The technique modulates a DC power bus through capacitive coupling using on-off keying (OOK), and receives and demodulates the signal from the DC power bus through the same capacitive coupling. The circuit acts as serial modem for the physical power line communication. The circuit and technique can be made of commercially available components or included in an application specific integrated circuit (ASIC) design, which allows for the circuit to be included in current designs with additional circuitry or embedded into new designs. This device and technique moves computational, sensing, and actuation abilities closer to the source, and allows for the networking of multiple similar nodes to each other and to a central processor. This technique also allows for reconfigurable systems by adding or removing nodes at any time. It can do so using nothing more than the in situ power wiring of the system.

Krasowski, Michael J.↗
DOE OSTI2025

Development for Integrated System-Level Analysis Capabilities in SAM for Molten Salt Reactors

In recent years, there has been renewed interest in Molten Salt Reactors (MSRs) for their potential advantages compared to reactors that rely on solid fuel. In response to such interest, many methods and codes have been developed to capture the unique features of MSRs. Among them, the System Analysis Module (SAM) is a modern system analysis tool that provides fast-running, modest-fidelity, whole-plant transient analysis capabilities, essential for fast-turnaround design scoping and engineering analyses of advanced reactor concepts. For liquid-fuel MSRs, the complex physics and chemistry involved in MSR operation—such as reactor kinetics, fluid flow, heat transfer, and salt composition dynamics—pose significant challenges for system-level modeling. Specific modeling capabilities are needed for system-level transient simulation. This paper presents recent advancements in SAM capability enhancements for system-level modeling of MSRs, focusing on improved simulation fidelity, computational efficiency, and multi-physics integration. Key enhancements include the development of species transport, Delayed Neutron Precursor (DNP) drift, modified Point Kinetics Equations (PKE), decay heat modeling, key fission product behavior, salt corrosion, and thermal-hydraulic coupling, as well as code robustness and performance enhancements for MSR applications. The code enhancement allows for better predictive accuracy in safety analysis, transient behavior, and operational optimization, thus supporting the design and licensing of next-generation MSRs. Results from case studies are presented to demonstrate the benefits of these enhancements in accurately capturing key reactor transient behaviors.

Hu, Rui (ORCID:0000000237712920)↗
DOE OSTI2022

On the theory of nonhomogeneous nonequilibrium superconductivity in 2D systems with massless fermions

Here we analyze static and nonequilibrium superconducting properties of a 2D relativistic-like model system with local electron-electron interaction, Rashba spin-orbit interaction αR in presence of time-dependent in-plane magnetic field H(t). It is shown that similar to the 2D case with ordinary massive quasiparticle dispersion ε(k)~|k| 2 at large fields, such a system demonstrates a nonhomogeneous superconducting stripe phase with the order parameter Δ(r)=Δ(0)cos(2[μ B B×r]n/$\hbarυ_F$) (B is the magnetic induction, υF is the Fermi velocity, n is the normal to the plane, μ B is the Bohr magneton, and αR$\ll$υF) where the stripes are oriented along the B direction. In the considered system, the inter-stripe period L and the magnitude of the magnetic field B are related by a universal relation BL=$\hbarυ_F$/μ B ≃0.714∙10 -4 Tm. Contrary to the case of massive quasiparticles, where the condition α R ~υ F can be, in principle, satisfied by increasing α R or by charge doping (Fermi velocity decreasing), in a relativistic-like system, where υF is doping-independent and one-two orders of magnitude larger than typical Fermi velocity in the “standard” 2D systems, the stripe phase can be the ground state at a rather low doping level. We also analyzed the nonequilibrium properties of the system with a focus on the melting of the stripe order (when the magnetic field is quenched to a lower value) and stripe dynamics (when the field is rotated by 90° degrees) and found several notable results. In particular, it was shown that the stripe domains melt according to law R~1$\sqrt{t}$ at initial times, while at longer times they shrink exponentially. In the case of the flipped magnetic field, the stripe orientation gradually turns from x- to y-direction, and the intermediate “crossed-stripe” phase takes place during times of order of picoseconds. Such a crossed phase is built of periodic superconducting bubbles that potentially may have applications in modern ultrafast superconducting technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1997

Development of Curved-Plate Elements for the Exact Buckling Analysis of Composite Plate Assemblies Including Transverse-Shear Effects

The analytical formulation of curved-plate non-linear equilibrium equations including transverse-shear-deformation effects is presented. The formulation uses the principle of virtual work. A unified set of non-linear strains that contains terms from both physical and tensorial strain measures is used. Linearized, perturbed equilibrium equations (stability equations) that describe the response of the plate just after buckling occurs are then derived after the application of several simplifying assumptions. These equations are then modified to allow the reference surface of the plate to be located at a distance z(sub c) from the centroidal surface. The implementation of the new theory into the VICONOPT exact buckling and vibration analysis and optimum design computer program is described as well. The terms of the plate stiffness matrix using both Classical Plate Theory (CPT) and first-order Shear-Deformation Plate Theory (SDPT) are presented. The necessary steps to include the effects of in-plane transverse and in-plane shear loads in the in-plane stability equations are also outlined. Numerical results are presented using the newly implemented capability. Comparisons of results for several example problems with different loading states are made. Comparisons of analyses using both physical and tensorial strain measures as well as CPT and SDPF are also made. Results comparing the computational effort required by the new analysis to that of the analysis currently in the VICONOPT program are presented. The effects of including terms related to in-plane transverse and in-plane shear loadings in the in-plane stability equations are also examined. Finally, results of a design-optimization study of two different cylindrical shells subject to uniform axial compression are presented.

McGowan, David Michael↗
NASA NTRS2022

Review of Solar Energetic Particle Models

Solar Energetic Particle (SEP) events are interesting from a scientific perspective as they are the product of a broad set of physical processes from the corona out through the extent of the heliosphere, and provide insight into processes of particle acceleration and transport that are widely applicable in astrophysics. From the operations perspective, SEP events pose a radiation hazard for aviation, electronics in space, and human space exploration, in particular for missions outside of the Earth’s protective magnetosphere including to the Moon and Mars. Thus, it is critical to improve the scientific understanding of SEP events and use this understanding to develop and improve SEP forecasting capabilities to support operations. Many SEP models exist or are in development using a wide variety of approaches and with differing goals. These include computationally intensive physics-based models, fast and light empirical models, machine learning-based models, and mixed-model approaches. The aim of this paper is to summarize all of the SEP models currently developed in the scientific community, including a description of model approach, inputs and outputs, free parameters, and any published validations or comparisons with data.

Kathryn Whitman↗
DOE OSTI2024

Phonon screening and dissociation of excitons at finite temperatures from first principles

The properties of excitons, or correlated electron–hole pairs, are of paramount importance to optoelectronic applications of materials. A central component of exciton physics is the electron–hole interaction, which is commonly treated as screened solely by electrons within a material. However, nuclear motion can screen this Coulomb interaction as well, with several recent studies developing model approaches for approximating the phonon screening of excitonic properties. While these model approaches tend to improve agreement with experiment, they rely on several approximations that restrict their applicability to a wide range of materials, and thus far they have neglected the effect of finite temperatures. Here, we develop a fully first-principles, parameter-free approach to compute the temperature-dependent effects of phonon screening within the ab initio GW -Bethe–Salpeter equation framework. We recover previously proposed models of phonon screening as well-defined limits of our general framework, and discuss their validity by comparing them against our first-principles results. We develop an efficient computational workflow and apply it to a diverse set of semiconductors, specifically AlN, CdS, GaN, MgO, and SrTiO 3 . We demonstrate under different physical scenarios how excitons may be screened by multiple polar optical or acoustic phonons, how their binding energies can exhibit strong temperature dependence, and the ultrafast timescales on which they dissociate into free electron–hole pairs.

Science & Technology - Other Topics↗