Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical method”

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 OSTI2020

Improvements to a class of hybrid methods for radiation transport: Nyström reconstruction and defect correction methods

In this study, two modifications are introduced for improving the accuracy, versatility, and robustness of a class of hybrid methods for radiation transport. In general, such methods are constructed by splitting the radiative flux into collided and uncollided components to which low- and high-resolution angular approximations are applied, respectively. In this work we focus on discrete ordinates discretizations of high and low order. The first modification we introduce changes the way in which the collided component is mapped into the uncollided component at the end of each time step in a simulation. The new mapping is a Nyström-type reconstruction that is applicable to arbitrary discrete ordinates quadratures, is guaranteed to preserve positivity of the solution provided that all ordinate weights are positive, is significantly more accurate than previous methods, and can be readily extended to other discretizations such as moment methods, finite element methods, and diffusion approximations. The second modification leverages integral deferred correction (IDC) to iteratively correct for the splitting error introduced by the inconsistency in angular discretization between the collided and uncollided components, in addition to improving the accuracy of the low-order temporal error that is treated by traditional IDC methods. Numerical tests in one- and two-dimensional geometries are used to demonstrate the increased accuracy and efficiency of the proposed modifications. It is found that the two techniques combined yield methods with solution accuracy and memory requirements comparable to that of monolithic discrete ordinates methods while reducing runtime by as much as a factor of between two and ten, depending on the problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Pronghorn: A Multidimensional Coarse Mesh Application for Advanced Reactor Thermal-Hydraulics

This paper presents an overview of Pronghorn, a multiscale thermal-hydraulic (T/H) application developed by Idaho National Laboratory and the University of California, Berkeley. Pronghorn, built on the open-source finite element Multiphysics Object-Oriented Simulation Environment (MOOSE), leverages state-of-the-art physical models, numerical methods, and nonlinear solvers to deliver fast-running advanced reactor T/H simulation capabilities within a modern software engineering environment. This work summarizes the physical models, multiphysics and multiscale coupling, and numerical discretization in Pronghorn with emphasis on our initial target application to pebble bed reactors (PBRs). A diverse set of applications are shown to depressurized natural circulation in the SANA experiments, forced convection in the Pebble Bed Modular Reactor, three-dimensional (3-D)/one-dimensional coupling of Pronghorn and RELAP-7 systems T/H for loop analysis in the High Temperature Reactor Power Module, and forced convection in the Mark-1 Pebble Bed Fluoride-Salt-Cooled High-Temperature Reactor. A multiphysics coupling of Pronghorn, RELAP-7, and Griffin deterministic neutronics for a gas-cooled PBR demonstrates the capability of the MOOSE framework for reactor design calculations. These applications highlight the verification and validation underlying Pronghorn’s software development while emphasizing features that improve upon capabilities offered by legacy tools in areas such as 3-D unstructured meshing, physics modeling, and multiphysics coupling.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Algorithm to calculate off-plane magnetic field from an on-plane field map

An algorithm is presented to calculate the field components [B x (x,y,z), B y (x,y,z), B z (x,y,z)] of the magnetic field at a point (x,y,z) in space, from the knowledge of the components [B x (x,y = 0,z), B y (x,y = 0,z), B z (x,y = 0,z)] on a ”reference plane”, which is normal to the y-axis at y=0. The algorithm, which is a general one and is not restricted to fields with mid-plane symmetry is based on the Taylor series expansion of the magnetic field components at any point in space in terms of the distance (y) of the point from the reference plane. The coefficients of the Taylor series expansion are expressed in terms of the on-plane field components and their partial derivatives with respect to spatial coordinates (x,z). The field components are usually generated from magnetic field measurements on a rectangular grid on the plane. The required partial derivatives of the field components on the plane can be computed with two methods; numerically, or by fitting a polynomial function to the experimentally measured values of the field components and subsequently take the required partial derivatives of the polynomial function. Each of these two method which calculates the partial derivatives of the field components on the plane, is applied to calculate the off-plane field components of the magnetic field generated by a modified Halbach magnet [1] of inner radius of 4.4 cm and it is found that each method provides accurate results with a relative error <1% deviation from the ”ideal” fields if the y distance of the point from the reference plane is smaller than 2 cm. This algorithm is part of the RAYTRACE computer code [2] and has been applied [3] on a dipole magnet with median plane symmetry.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

Visualization of the Mode Shapes of Pressure Oscillation in a Cylindrical Cavity

Our work describes a novel experimental method to visualize the mode shapes of pressure oscillation in a cylindrical cavity. Acoustic resonance in a cavity is a grand old problem that has been under investigation (using both analytical and numerical methods) for more than a century. In this article, a novel method based on high speed imaging of combustion chemiluminescence was presented to visualize the mode shapes of pressure oscillation in a cylindrical cavity. By generating high-temperature combustion gases and strong pressure waves simultaneously in a cylindrical cavity, the pressure oscillation can be inferred due to the chemiluminescence emissions of the combustion products. We can then visualized the mode shapes by reconstructing the images based on the amplitudes of the luminosity spectrum at the corresponding resonant frequencies. Up to 11 resonant mode shapes were clearly visualized, each matching very well with the analytical solutions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2012

Solving PDEs with Intrepid

Intrepid is a Trilinos package for advanced discretizations of Partial Differential Equations (PDEs). The package provides a comprehensive set of tools for local, cell-based construction of a wide range of numerical methods for PDEs. This paper describes the mathematical ideas and software design principles incorporated in the package. We also provide representative examples showcasing the use of Intrepid both in the context of numerical PDEs and the more general context of data analysis.

Bochev, P.↗
DOE OSTI2013

Modern methods for calculating photoionization and electron-impact ionization of two-electron atoms and molecules

A number of recently developed methods for calculating multifold differential cross sections for photoionization and electron impact ionization of atoms and molecules with two active electrons are reviewed. The methods are based on unique approaches to the calculation of three-body Coulomb wave functions. The exterior scaling method and the driven Schrödinger equation formalism are considered. The effectiveness of the time-dependent approaches to the scattering problem, such as the paraxial approximation and the time-dependent scaling, is demonstrated. A novel numerical method is formulated, which has been developed by the authors to solve the six-dimensional Schrödinger equation for an atom with two active electrons on the basis of the Chang-Fano transformation and the discrete variable representation. The threshold behavior manifested by the angular distributions of the two-electron photoionization of a negative hydrogen ion and a helium atom and the multifold differential cross sections for the electronimpact ionization of hydrogen and nitrogen molecules are analyzed on the basis of numerical simulations. The Wannier law for the angular distribution of double ionization is demonstrated to be incorrect even at very low energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

A Scalable Reduced‐Order Model for the Steady Navier–Stokes Equations

Scaling up new scientific technologies from laboratory to industry often involves demonstrating performance on a larger scale. Computer simulations can accelerate design and predictions in the deployment process, though traditional numerical methods are computationally intractable even for intermediate pilot plant scales. Recently, the component reduced order modeling method has been developed to tackle this challenge by combining projection reduced order modeling and discontinuous Galerkin domain decomposition. However, while many scientific or engineering applications involve nonlinear physics, this method has only been demonstrated for various linear systems. In this work, the component reduced order modeling method is extended to steady Navier–Stokes flow, with application to general nonlinear physics in view. The large‐scale, global domain is decomposed into a combination of small‐scale unit component. Linear subspaces for flow velocity and pressure are identified via proper orthogonal decomposition over sample snapshots collected from each small‐scale unit component. Velocity bases are augmented with a pressure supremizer to satisfy the inf–sup condition for stable pressure prediction. Two different nonlinear reduced order modeling methods are employed and compared for efficient evaluation of nonlinear advection: A third‐order tensor projection operator and the empirical quadrature procedure. The proposed method is demonstrated on the flow over arrays of five different unit objects, achieving a 23‐fold speedup with less than 4% relative error in domains up to 256 times larger than the unit components. Furthermore, a numerical experiment with the pressure supremizer strongly indicates the need for a supremizer for stable pressure prediction. A comparison between the tensorial approach and the empirical quadrature procedure revealed a slight advantage of the empirical quadrature procedure. The framework is compared with an alternating Schwarz‐based reduced‐order approach, demonstrating improved efficiency and robustness for the DG‐based global solver while retaining flexibility for sub‐scale iterative solvers. The method is further extended to a coupled advection–diffusion and Navier–Stokes system, illustrating its applicability to multi‐physics problems and its potential for more general, inter‐coupled nonlinear systems.

42 ENGINEERING↗
DOE OSTI2019

Simulating bubble dynamics in a buoyant system

Summary Multiphase flows are critical components of many physical systems; however, numerical models of multiphase flows with large parameter gradients can be challenging. Here, two different numerical methods, volume of fluid (VOF) and smoothed particle hydrodynamics (SPH), are used to model the buoyant rise of isolated gas bubbles through quiescent fluids for a range of Bond and Reynolds numbers. The VOF is an Eulerian grid–based method, whereas the SPH is Lagrangian and mesh free. Each method has unique strengths and weaknesses, and a comparison of the two approaches as applied to multiphase phenomena has not previously been performed. The VOF and SPH simulations are compared, verified, and validated. Results using two‐dimensional VOF and SPH simulations are similar to each other and are able to reproduce numerical benchmarks and experimental results for sufficiently large Morton and Reynolds numbers. It is also shown that at low Reynolds numbers, the two methods, SPH and VOF, diverge in the transient regime of the bubble rise. Regimes that require simulations capable of representing three‐dimensional drag are identified as well as regimes in which results from VOF and SPH diverge.

Arai, Erin↗
DOE OSTI2024

Using Hydrodynamic Similarity as a Verification Method for Impact Cratering Simulations in the FLAG Hydrocode

Hydrodynamic codes (hydrocodes) are common tools for modeling hypervelocity impacts to provide insight into the physical phenomenon. Hydrocodes can simulate impacts from micrometer to kilometer spatial scales and reach impact velocities difficult to achieve in experimental settings. However, numerical models are approximations, and demonstrating that a numerical method is capable of providing physical results for these models is essential. In this work, we employ a hydrocode verification technique that leverages hydrodynamic similarity, a mathematical property of the conservation equations of fluid mechanics that form the basis for hydrocode models. Using the FLAG hydrocode, we simulate aluminum (Al) and basalt projectiles and targets at spatial scales spanning 7 orders of magnitude (hundreds of micrometers to kilometers). These materials were chosen because Al-6061 is a common material in spacecraft and satellites and basalt is a useful approximation of rocky astronomical bodies. Our results show that hydrodynamic similarity holds for each material model used and across spatial scales. We show that under certain conditions hydrodynamic similarity can apply in the presence of gravity and that similarity does not hold in the presence of strength models. We conclude that the FLAG hydrocode preserves important mathematical properties of fluid dynamics in hypervelocity impacts of Al-6061 and basalt.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

An efficient cost reduction procedure for bounded-control LQR problems

A novel approach has been developed for approximating the solution to the constrained LQR problem, based on updating the final state and costate of a related regular problem, and on slightly shifting the switching times (the instants when the control meets the constraints). The main result is the expression of a suboptimal control in feedback form using the solution of some compatible Riccati equation. The gradient method is applied to reduce the cost via explicit algebraic formula for its partial derivatives with respect to the hidden final state/costate of the related regular problem and to the switching times. The numerical method is termed efficient because it does not involve integrations of states or cost trajectories, and reduces to its minimum the dimension of the unknown parameters at the final condition. All the relevant objects are calculated from a few auxiliary matrices, which are computed only once. The scheme is here applied to two case studies whose optimal solutions are known. The first example is a two-dimensional model of the ‘cheapest stop of a train’ problem. The second one refers to the temperature control of a metallic strip leaving a multi-stand rolling mill, a problem with a high-dimensional state.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Theoretical methods to design and test quantum simulators for the compact Abelian Higgs model

The lattice compact Abelian Higgs model is a nonperturbative regularized formulation of low-energy scalar quantum electrodynamics. In 1+1 dimensions, this model can be quantum simulated using a ladder-shaped optical lattice with Rydberg-dressed atoms. In this setup, one spatial dimension is used to carry the angular momentum of the quantum rotors. One can use truncations corresponding to spin-2 and spin-1 to build local Hilbert spaces associated with the links of the lattice. We argue that ladder-shaped configurable arrays of Rydberg atoms can be used for the same purpose. We make concrete proposals involving two and three Rydberg atoms to build one local spin-1 space (a qutrit). We show that the building blocks of the Hamiltonian calculations are models with one and two spins. We compare target and simulators using perturbative and numerical methods. The two-atom setup provides an easily controllable simulator of the one-spin model while the three-atom setup involves solving nonlinear equations. We discuss approximate methods to couple two spin-1 spaces. The article provides analytical and numerical tools necessary to design and build the proposed simulators with current technology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Correction of steel casing effect for density log using numerical and experimental methods in the slim borehole

Density log is widely applied for a variety of fields such as the petroleum exploration, mineral exploration, and geotechnical survey. The logging condition of density log is normally open holes but there are frequently cased boreholes. The primary calibration curve by slim hole logging manufacturer is normally the calibration curves for the variation of borehole diameter. In this study, we have performed the correction of steel casing effects using numerical and experimental methods. We have performed numerical modeling using the Monte Carlo N-Particle (MCNP) code based on Monte Carlo method, and field experimental method from open and cased hole log. In this study, we used the FDGS (Formation Density Gamma Sonde) for slim borehole with a 100 mCi 137Cs source, three inch borehole and steel casing. The casing effect between numerical and experimental method is well matched.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Using Small-Angle Neutron Scattering to Detect Nanoscopic Lipid Domains

The cell plasma membrane is a complex system, which is thought to be capable of exhibiting non-random lateral organization. Studies of live cells and model membranes have yielded mechanisms responsible for the formation, growth, and maintenance of nanoscopic heterogeneities, although the existence and mechanisms that give rise to these heterogeneities remain controversial. Small-angle neutron scattering (SANS) is a tool ideally suited to interrogate lateral heterogeneity in model membranes, primarily due to its unique spatial resolution (i.e., "5 100 nm) and its ability to resolve structure with minimal perturbation to the membrane. In this review we examine several methods used to analyze the SANS signal arising from freely suspended unilamellar vesicles containing lateral heterogeneity. Specifically, we discuss an analytical model for a single, round domain on a spherical vesicle. We then discuss a numerical method that uses Monte Carlo simulation to describe systems with multiple domains and/or more complicated morphologies. Also discussed are several model-independent approaches that are sensitive to membrane heterogeneity. The review concludes with several recent applications of SANS to the study of membrane raft mixtures.

Pan, Jianjun↗
DOE OSTI2015

On control of continuous dynamical polysystems in discrete times

This paper considers control systems with a finite number of control parameters, i.e. dynamical polysystems. It is assumed that control switchings can occur only at certain discrete times. Statistical methods in number theory are used to investigate the controllability of these systems. Existence of control switching times is established through analysis of Diophantine equations. The values of control switching times are found through suitable numerical methods. Certain assertions on the controllability of the polysystems of the class under consideration are proved. Some examples are considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Brownian dynamics of confined suspensions of active microrollers

In this study, we develop efficient numerical methods for performing many-body Brownian dynamics simulations of a recently observed fingering instability in an active suspension of colloidal rollers sedimented above a wall [M. Driscoll, B. Delmotte, M. Youssef, S. Sacanna, A. Donev, and P. Chaikin, Nat. Phys. (2016), preprint arXiv:1609.08673. We present a stochastic Adams-Bashforth integrator for the equations of Brownian dynamics, which has the same cost but is more accurate than the widely used Euler-Maruyama scheme, and use a random finite difference to capture the stochastic drift proportional to the divergence of the configuration-dependent mobility matrix. We generate the Brownian increments using a Krylov method and show that for particles confined to remain in the vicinity of a no-slip wall by gravity or active flows, the number of iterations is independent of the number of particles. Our numerical experiments with active rollers show that the thermal fluctuations set the characteristic height of the colloids above the wall, both in the initial condition and the subsequent evolution dominated by active flows. Finally, the characteristic height in turn controls the time scale and wavelength for the development of the fingering instability

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Analytical calculation for the solid angle subtended at an arbitrarily positioned point source by an ellipse - Paper 38

The geometrical solid angle can be used to calibrate radiation detectors or to determine the activity of radioactive sources. In the present work, we introduce a direct analytical method for calculating the solid angle subtended at an arbitrarily positioned point source by a closed contour (ellipse). The validity of the derived analytical expressions was successfully confirmed by the comparison with some published data (Numerical Method). Direct mathematical expressions to calculate the geometrical solid angle subtended by an elliptical disk detector with respect to arbitrarily positioned radiating point sources have been derived. The agreement between the calculated and the published solid angle values is excellent; the discrepancies were (0.0%). This means that the present approach is efficient and sufficiently powerful to evaluate the geometrical efficiency of the elliptical disk detectors. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2018

Evaluation of Neutron-induced Cross Sections and their Related Covariances with Physical Constraints

Nuclear data, along with numerical methods and the associated calculation schemes, continue to play a key role in reactor design, reactor core operating parameters calculations, fuel cycle management and criticality safety calculations. Due to the intensive use of Monte-Carlo calculations reducing numerical biases, the final accuracy of neutronic calculations increasingly depends on the quality of nuclear data used. This paper gives a broad picture of all ingredients treated by nuclear data evaluators during their analyses. After giving an introduction to nuclear data evaluation, we present implications of using the Bayesian inference to obtain evaluated cross sections and related uncertainties. In particular, a focus is made on systematic uncertainties appearing in the analysis of differential measurements as well as advantages and drawbacks one may encounter by analyzing integral experiments. The evaluation work is in general done independently in the resonance and in the continuum energy ranges giving rise to inconsistencies in evaluated files. For future evaluations on the whole energy range, we call attention to two innovative methods used to analyze several nuclear reaction models and impose constraints. Finally, we discuss suggestions for possible improvements in the evaluation process to master the quantification of uncertainties. These are associated with experiments (microscopic and integral), nuclear reaction theories and the Bayesian inference.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗