Search NASA⌕ Search

SEARCH · Search NASA

Results for “method of characteristics”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

The nonsinglet spin-dependent structure function evolution by Laplace and characteristics methods

We evaluate the non-singlet spin-dependent structure function g{sub 1}{sup NS} at leading order (LO) and next-to-leading order (NLO) by using the Laplace-transform technique and method of characteristics and also obtain its first moment at NLO. The polarized non-singlet structure function results are compared with the data from HERMES (A. Airapetian et al., Phys. Rev. D 75, 012007 (2007)) and E143 (K. Abe et al. (E143 Collab.), Phys. Rev. D 58, 112003 (1998)) at LO and NLO analyses and the first-moment the result at NLO is compared with the result of the NLO GRSV2000 fit. Considering the solution, this method is valid at low- and large-x regions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics

In this work, we present a conservative numerical method for radiation magnetohydrodynamics with frequency-dependent full transport in stationary spacetimes. This method is stable and accurate for both large and small optical depths and radiation pressures. The radiation stress–energy tensor is evolved in flux-conservative form, and closed with a swarm of samples that each transport a multigroup representation of the invariant specific intensity along a null geodesic. In each zone, the enclosed samples are used to efficiently construct a Delaunay triangulation of the unit sphere in the comoving frame, which in turn is used to calculate the Eddington tensor, average source terms, and adaptively refine the sample swarm. Furthermore, radiation four-forces are evaluated in the moment sector in a semi-implicit fashion. The radiative transfer equation is solved in invariant form deterministically for each sample. Since each sample carries a discrete representation of the full spectrum, the cost of evaluating the transport operator is independent of the number of frequency groups, representing a significant reduction of algorithmic complexity for transport in frequency-dependent problems. The major approximation we make in this work is performing scattering in an angle-averaged way. Local adaptivity in samples also makes this scheme more amenable to nonuniform meshes than a traditional Monte Carlo method. We describe the method and present results on a suite of test problems. We find that Method of Characteristics Moment Closure converges at least as ~N -1 , rather than the canonical Monte Carlo N -1/2 , where N is the number of samples per zone.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A time-dependent neutron transport method of characteristics formulation with time derivative propagation

A new time-dependent Method of Characteristics (MOC) formulation for nuclear reactor kinetics was developed utilizing angular flux time-derivative propagation. This method avoids the requirement of storing the angular flux at previous points in time to represent a discretized time derivative; instead, an equation for the angular flux time derivative along 1D spatial characteristics is derived and solved concurrently with the 1D transport characteristic equation. This approach allows the angular flux time derivative to be recast principally in terms of the neutron source time derivatives, which are approximated to high-order accuracy using the backward differentiation formula (BDF). This approach, called Source Derivative Propagation (SDP), drastically reduces the memory requirements of time-dependent MOC relative to methods that require storing the angular flux. An SDP method was developed for 2D and 3D applications and implemented in the computer code DeCART in 2D. DeCART was used to model two reactor transient benchmarks: a modified TWIGL problem and a C5G7 transient. The SDP method accurately and efficiently replicated the solution of the conventional time-dependent MOC method using two orders of magnitude less memory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

The Random Ray Method Versus Multigroup Monte Carlo: The Method of Characteristics in OpenMC and SCONE

The Random Ray Method (TRRM) is a recently developed approach to solving neutral particle transport problems based on the Method of Characteristics. While the method previously has been implemented only in closed-source or limited-functionality codes, this work describes its implementation in two open-source Monte Carlo codes: OpenMC and SCONE. The random ray implementations required small modifications to the existing Multigroup Monte Carlo (MGMC) solvers, offering a rare venue for redundant, fine-grained, "apples-to-apples" speed and accuracy comparisons between transport methods. To this end, TRRM and MGMC solvers are evaluated against each other using each code's native capabilities on reactor eigenvalue problems with different degrees of energy discretization. On the C5G7 benchmark (featuring only seven energy groups), TRRM achieves a maximum pin power error comparable to or lower than that of MGMC for a given run time. On a problem with 69 energy groups, MGMC is found to scale more efficiently, obtaining a lower pin power error for a given run time. However, the defining difference between the two transport methods is found to be their vastly different uncertainty distributions. Specifically, TRRM is found to maintain similar levels of accuracy and uncertainty throughout the simulation domain whereas MGMC can exhibit orders-of-magnitude greater errors in areas of the problem that feature low neutron flux. For instance, TRRM provided an up to 373 times speed advantage compared with MGMC for computing the flux in low-flux regions in the moderator surrounding the C5G7 core.

42 ENGINEERING↗
DOE OSTI2023

Adaptive Angular Quadrature Scheme for a backwards-in-time Method of Characteristics Solution to the Radiative Transfer Equation [Slides]

Radiative transfer/radiation transport are important problems to solve in astrophysics and high energy density physics. Various methods exist to solve radiation transport, such as Monte Carlo (MC), Discrete Ordinates (S N ), Method of Characteristics (MOC), and the spherical harmonics (P N ) method. Method of Characteristics requires “launching” of rays in discrete directions. Unresolved details of angular mesh create ray effects and can miss sources in the domain. Ray effects can lead to unphysical “stepping” in solution and incorrect energy deposition. Adaptive quadrature schemes can be used to detect and mitigate these effects. The Method of Characteristics (MOC) is a common method for solving hyperbolic PDEs in radiation transport and supersonic flow problems. Generally in MOC for radiation transport, virtual particles are tracked from birth to the end of a timestep. This requires interpolation to go from final location to cell averaged or corner values of angular intensity. Backwards-in-Time (BIT) particle tracking avoids this by prescribing the final position of the virtual particle at the cell nodes/corners. Angular intensities are computed at time k + 1 by launching ray back to previous timestep(s), or t = 0. Scheme allows solution to be computed as the characteristic ray is traced backwards in time.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

The Subray Method of Characteristics

The MPACT code is being jointly developed by the University of Michigan and Oak Ridge National Laboratory. It uses the 2-D/1-D method to solve neutron transport problems for reactors. The 2-D/1-D method decomposes the problem into a stack of 2-D planes and uses a high-fidelity transport method to resolve all heterogeneity in each plane. These planes are then coupled axially, using a lower-order solver. With this scheme, three-dimensional (3-D) solutions to the transport equation can be obtained at a much lower cost. The 2-D/1-D method assumes that the materials are axially homogeneous for each 2-D plane. Violation of this assumption requires homogenization, which can significantly reduce the accuracy of the calculation. This work proposes the subray method of characteristics (subray MOC) as a solution to this problem. Subray MOC is a subgrid method that allows local heterogeneities to be directly resolved by method of characteristics while treating the rest of the 2-D plane as axially uniform. This improves the accuracy in the neighborhood of the heterogeneity while minimizing the increase in run time. The method was applied to variations of the C5G7 benchmark problems and compared with a previously developed subgrid method called the subplane collision probabilities (SCP) method. Comparisons were made among results obtained using subray MOC, the SCP method, and no subgrid method. Subray MOC consistently performed best, reducing maximum 3-D power distribution errors from as high as 30% to 2% or less. Moreover, it consistently outperformed the SCP method with run times that were shorter than the reference calculations.

2-D/1-D method↗
DOE OSTI2016

Second-order method of characteristics through local cell refinement

A new approach to the method of characteristics is developed for 2D unstructured triangular mesh that is locally second order accurate with respect to mesh cell size. The second order accuracy is obtained through a local cell refinement and flux reconstruction process. Numerical examples are provided to demonstrate the theoretical error prediction for simple test problems. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Quantifying hydrate solidification front advancing using method of characteristics

Here, we develop a one-dimensional analytical solution based on the method of characteristics to explore hydrate formation from gas injection into brine-saturated sediments within the hydrate stability zone. Our solution includes fully coupled multiphase and multicomponent flow and the associated advective transport in a homogeneous system. Our solution shows that hydrate saturation is controlled by the initial thermodynamic state of the system and changed by the gas fractional flow. Hydrate saturation in gas-rich systems can be estimated by 1 – cl 0 /cl e when Darcy flow dominates, where cl 0 is the initial mass fraction of salt in brine, and cl e is the mass fraction of salt in brine at three-phase (gas, liquid, and hydrate) equilibrium. Hydrate saturation is constant, gas saturation and gas flux decrease, and liquid saturation and liquid flux increase with the distance from the gas inlet to the hydrate solidification front. The total gas and liquid flux is constant from the gas inlet to the hydrate solidification front and decreases abruptly at the hydrate solidification front due to gas inclusion into the hydrate phase. The advancing velocity of the hydrate solidification front decreases with hydrate saturation at a fixed gas inflow rate. This analytical solution illuminates how hydrate is formed by gas injection (methane, CO 2 , ethane, propane) at both the laboratory and field scales.

58 GEOSCIENCES↗
DOE OSTI2017

A new approach to three-dimensional neutron transport solution based on the method of characteristics and linear axial approximation

A new approach based on the method of characteristics (MOC) is proposed to solve the neutron transport equation. A new three-dimensional (3D) spatial discretization is applied to avoid the instability issue of the transverse leakage iteration of the traditional 2D/1D approach. In this new approach, the axial and radial variables are discretized in two different ways: the linear expansion is performed in the axial direction, then, the 3D solution of the angular flux is transformed to be the planar solution of 2D angular expansion moments, which are solved by the planar MOC sweeping. Based on the boundary and interface continuity conditions, the 2D expansion moment solution is equivalently transformed to be the solution of the axially averaged angular flux. Using the piecewise averaged angular flux at the top and bottom surfaces of 3D meshes, the planes are coupled to give the 3D angular flux distribution. The 3D CMFD linear system is established from the surface net current of every 3D pin-mesh to accelerate the convergence of power iteration. The STREAM code is extended to be capable of handling 3D problems based on the new approach. Several benchmarks are tested to verify its feasibility and accuracy, including the 3D homogeneous benchmarks and heterogeneous benchmarks. The computational sensitivity is discussed. The results show good accuracy in all tests. With the CMFD acceleration, the convergence is stable. In addition, a pin-cell problem with void gap is calculated. This shows the advantage compared to the traditional 2D/1D MOC methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Conjugacy Analysis of the Finite-Difference Equations of the Method of Characteristics on the Basis of a Linear-Algebraic Formulation

A linear-algebraic form of the equations of the method of characteristics, which is used to approximate the neutron transport equation, is obtained. It is shown on the basis of the obtained linear-algebraic form that the discrete form of the conjugate equation differs from the algebraically discrete problem constructed by linear-algebraic transformations of the discrete form of the normal problem. The reason for the discrepancy lies in the approximation of the volumes of the spatial cells in covering the working region by a network of characteristics. It is shown by means of test calculations that when the network of characteristics is refined the solution of the conjugate transport equation converges to the solution of the algebraically conjugate problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Determination of the parameters of a holographic layer from its spectral characteristics

Methods for estimating the main parameters of holographic sensors (refractive index modulation depth and hologram thickness) from transmission spectra in the absence of absorption and light scattering are discussed. The consideration is performed for layers oriented parallel to the holographic layer surface under normal light incidence. Direct numerical solution of the problem of light propagation in a periodic nonabsorbing medium is used to study the reflection and transmission spectra of the holographic layer in a wide range of variation in its thickness and the refractive index modulation depth. A classification of the reflection regimes from the holographic layer is proposed (from weak reflection to the photonic crystal regime). A comparison with the results obtained by the coupled-wave analysis is performed, and the limitations of this method at a significant spectral detuning from resonance and under conditions of strong reflection are revealed. It is shown that the main hologram parameters can be estimated from the experimental transmission spectrum of the phase hologram (in the case of strong reflection) based on the spectral dip parameters. (holographic sensors)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Linear source method of characteristics in HELIOS2

The Studsvik Scandpower code HELIOS2 is a fast lattice transport code that has gained wide use for reactor fuel cycle analysis in commercial, regulatory, and research environments. It contains two transport solvers - Collision Probabilities (CP) and Method of Characteristics (MoC) - with a general geometry ability and generalized depletion for a variety of fuel types, experimental facilities, and isotope production applications. The HELIOS2 lattice code has been extended with the development of a linear source MoC option. The linear source equations are derived within the scope of the HELIOS2 general geometry and non-cyclic ray tracing. The implementation was tested for fine and coarse spatial mesh models of the 2D C5G7 MOX computational benchmark, which provides a pure and extreme test of a transport solver. The linear source yields at least comparable accuracy with reduced computational resource requirements relative to the flat-source approximation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Preliminary study on applying the method of characteristics to pebble-bed HTR - 18637

Neutron transport in high temperature pebble-bed gas-cooled reactor, which is a promising generation IV reactor candidate with cylindrical core, graphite moderator and spherical TRISO coated particles fuel, is of great interest. However, the stochastic method namely Monte Carlo faces with tremendous simulation of particles and the conventional two-step method featured by homogenization only can give the scalar flux. In this paper, method of characteristics is trying to be applied to pebble bed reactor because of its advantages in arbitrary geometry and the ability to obtain solution in neutronically large domains. Some challenges occur when it comes to HTR. First, the amount of calculation and the memory requirement are huge due to the large quantities of spherical fuel and stochastic distribution of pebbles. Second, the arrangement of characteristics and ray tracing process is more complex comparing to LWR with cuboid geometry because the stochastic distribution of the spherical assemblies and the longer neutron mean free path in pebble bed HTR. Third, the treatment of the multiscale sphere fuel need be detail studied to determine the subdivision of flat source region. At present, preliminary research on this subject is undergoing. This paper is concentrated on the modeling of the pebble-bed reactor core by constructive solid geometry (CSG) which allows complex spatial models to be built using boolean operations, characteristics arrangement and the treatment of boundary condition. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗