Search NASA⌕ Search

SEARCH · Search NASA

Results for “CASE 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 991 records · Page 55

Convective flows in enclosures with vertical temperature or concentration gradients

The transport process in the fluid phase during the growth of a crystal has a profound influence on the structure and quality of the solid phase. In vertical growth techniques the fluid phase is often subjected to vertical temperature and concentration gradients. The main objective is to obtain more experimental data on convective flows in enclosures with vertical temperature or concentration gradients. Among actual crystal systems the parameters vary widely. The parametric ranges studied for mass transfer are mainly dictated by the electrochemical system employed to impose concentration gradients. Temperature or concentration difference are maintained between two horizontal end walls. The other walls are kept insulated. Experimental measurements and observations were made of the heat transfer or mass transfer, flow patterns, and the mean and fluctuating temperature distribution. The method used to visualize the flow pattern in the thermal cases is an electrochemical pH-indicator method. Laser shadowgraphs are employed to visualize flow patterns in the solutal cases.

Wang, L. W.↗

Fast computation of synthetic seismograms within a medium containing remote localized perturbations: a numerical solution to the scattering problem

We derive a fast discrete solution to the scattering problem. This solution allows us to compute accurate synthetic seismograms or waveforms for arbitrary locations of sources and receivers within a medium containing localized perturbations. The key to efficiency is that wave propagation modelling does not need to be carried out in the entire volume that encompasses the sources and the receivers but only within the subvolume containing the perturbations or scatterers. The proposed solution has important applications, for example, it permits the imaging of remote targets located in regions where no sources or receivers are present. Our solution relies on domain decomposition: within a small volume that contains the scatterers, wave propagation is modelled numerically, while in the surrounding volume, where the medium isn't perturbed, the response is obtained through wavefield extrapolation. The originality of this work is the derivation of discrete formulas for representation theorems and Kirchhof–Helmholtz integrals that naturally adapt to the numerical scheme employed for modelling wave propagation. Our solution applies, for example, to finite difference methods or finite/spectral elements methods. The synthetic seismograms obtained with our solution can be considered ‘exact’ as the total numerical error is comparable to that of the method employed for modelling wave propagation. We detail a basic implementation of our solution in the acoustic case using the finite difference method and present numerical examples that demonstrate the accuracy of the method. We show that ignoring some terms accounting for higher order scattering effects in our solution has a limited effect on the computed seismograms and significantly reduces the computational effort. Finally, we show that our solution can be used to compute localized sensitivity kernels and we discuss applications to target oriented imaging. Extension to the elastic case is straightforward and summarized in a dedicated section.

97 MATHEMATICS AND COMPUTING↗

Development of a BISON validation case for the TRISO transient irradiations in NSRR using effective heat capacity methods

The current tristructural isotropic (TRISO) fuel assessment and validation database in BISON primarily covers steady- state irradiation and high-temperature furnace testing. Transient assessment cases are potentially needed to support U.S. industry efforts in designing and deploying commercial reactors using TRISO fuels. Historical transient tests in- volving TRISO fuels used highly conservative conditions compared to the typical high-temperature gas-cooled reactor accident scenarios. Despite this, modeling historical transient tests is fundamental for evaluating BISON’s predictive capabilities, adapting material properties for high-temperature and high-particle-power regimes, and developing a sys- tematic validation approach for TRISO transient applications. This work developed a 1D model of transient experiments carried out at the Nuclear Safety Research Reactor using BISON. BISON’s predictions of energy deposition, UO 2 melting onset, and molten volume fractions are compared against experimental measurements. Melting was modeled using an effective specific heat capacity model for UO 2 . We found that BISON’s predictions are in reasonable agreement with experimental data for low-energy-deposition cases, and that BISON overpredicts melting at higher energy depositions. We also discuss the potential causes of discrepancies between the simulated and measured results and propose ways to further develop the model. Although these simulations used conservative conditions compared to those expected for actual TRISO-fueled reactors, they extended the range of conditions reflected in the data in the existing BISON database for TRISO fuels.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Coupled hydromechanical and electromagnetic disturbances in unsaturated porous materials

[1] A theory of cross-coupled flow equations in unsaturated soils is necessary to predict (1) electroosmotic flow with application to electroremediation and agriculture, (2) the electroseismic and the seismoelectric effects to develop new geophysical methods to characterize the vadose zone, and (3) the streaming current, which can be used to investigate remotely ground water flow in unsaturated conditions in the capillary water regime. To develop such a theory, the cross-coupled generalized Darcy and Ohm constitutive equations of transport are extended to unsaturated conditions. This model accounts for inertial effects and for the polarization of porous materials. Rather than using the zeta potential, like in conventional theories for the saturated case, the key parameter used here is the quasi-static volumetric charge density of the pore space, which can be directly computed from the quasi-static permeability. The apparent permeability entering Darcy’s law is also frequency dependent with a critical relaxation time that is, in turn, dependent on saturation. A decrease of saturation increases the associated relaxation frequency. The final form of the equations couples the Maxwell equations and a simplified form of two-fluid phases Biot theory accounting for water saturation. A generalized expression of the Richard equation is derived, accounting for the effect of the vibration of the skeleton during the passage of seismic waves and the electrical field. A new expression is obtained for the effective stress tensor. The model is tested against experimental data regarding the saturation and frequency dependence of the streaming potential coupling coefficient. The model is also adapted for two-phase flow conditions and a numerical application is shown for water flooding of a nonaqueous phase liquid (NAPL, oil) contaminated aquifer. Seismoelectric conversions are mostly taking place at the NAPL (oil)/water encroachment front and can be therefore used to remotely track the position of this front. This is not the case for other geophysical methods.

58 GEOSCIENCES↗

SU-E-J-184: Volumetric Indices to Aid Definition of Respiratory Gates with Particular Reference to Lung Stereotactic Body Radiotherapy (SBRT)

Purpose: Conventional definition of respiratory gates relies on tumor motion determination in limited planes. We are proposing a new method to define the RPM gates in a consistent manner ensuring that the tumor motion is restricted to 4 mm in 3D for lung SBRT patients. The method relies on studying the ratio of volumes obtained by GTVnn intersection with GTV50+2mm margin and GTVnn where GTVnn is the GTV volume in phase nn (=0,10,30.) while GTV50+2mm is a pseudo structure created by adding an isotropic margin of 2mm to GTV50. If for any phase nn, above ratio equals 1, it ensures that the tumor motion is ≤ ±2 mm in 3D from GTV50. Methods: This method was tested for 50 patients (14-Central, 36-peripheral) to determine the RPM gates which were then compared with the gates used clinically. The minimum cut-off value of the above coefficient for its inclusion of a phase in RPM gate was taken as 0.97 for central and 0.95 for peripheral tumors. Results: 15 (30%) of the patients did not require any change in the RPM gates w.r.t. gates defined using conventional motion assessment methods. In 15(30%) cases, the RPM gates could have been smaller while in remaining 20 patients, gates could have been larger. 5(/14) patient’s central tumors and 10 (/36) peripheral tumors did not need any gate change. 8(/50) patients could have RPM gate change of 30% while 10(/50) could have a gate change of up to 20%. 10, 20 & 30% RPM gate change could have happened for 11, 10 & 9 patients, respectively. Conclusion: Proposed volumetric indices based method allows a consistent, scientific and objective method to decide optimal RPM gates which is free from any inter or intra person variability and satisfies the tumor motion limits as defined by AAPM TG-76 in totality.

60 APPLIED LIFE SCIENCES↗

Proceedings of the 2004 Workshop on CFD Validation of Synthetic Jets and Turbulent Separation Control

The papers presented here are from the Langley Research Center Workshop on Computational Fluid Dynamics (CFD) Validation of Synthetic Jets and Turbulent Separation Control (nicknamed "CFDVAL2004"), held March 2004 in Williamsburg, Virginia. The goal of the workshop was to bring together an international group of CFD practitioners to assess the current capabilities of different classes of turbulent flow solution methodologies to predict flow fields induced by synthetic jets and separation control geometries. The workshop consisted of three flow-control test cases of varying complexity, and participants could contribute to any number of the cases. Along with their workshop submissions, each participant included a short write-up describing their method for computing the particular case(s). These write-ups are presented as received from the authors with no editing. Descriptions of each of the test cases and experiments are also included.

Rumsey, Christopher L.↗

Rate-Compatible LDPC Codes with Linear Minimum Distance

A recently developed method of constructing protograph-based low-density parity-check (LDPC) codes provides for low iterative decoding thresholds and minimum distances proportional to block sizes, and can be used for various code rates. A code constructed by this method can have either fixed input block size or fixed output block size and, in either case, provides rate compatibility. The method comprises two submethods: one for fixed input block size and one for fixed output block size. The first mentioned submethod is useful for applications in which there are requirements for rate-compatible codes that have fixed input block sizes. These are codes in which only the numbers of parity bits are allowed to vary. The fixed-output-blocksize submethod is useful for applications in which framing constraints are imposed on the physical layers of affected communication systems. An example of such a system is one that conforms to one of many new wireless-communication standards that involve the use of orthogonal frequency-division modulation

Divsalar, Dariush↗

SU-F-I-41: Calibration-Free Material Decomposition for Dual-Energy CT

Purpose: To eliminate tedious phantom calibration or manually region of interest (ROI) selection as required in dual-energy CT material decomposition, we establish a new projection-domain material decomposition framework with incorporation of energy spectrum. Methods: Similar to the case of dual-energy CT, the integral of the basis material image in our model is expressed as a linear combination of basis functions, which are the polynomials of high- and low-energy raw projection data. To yield the unknown coefficients of the linear combination, the proposed algorithm minimizes the quadratic error between the high- and low-energy raw projection data and the projection calculated using material images. We evaluate the algorithm with an iodine concentration numerical phantom at different dose and iodine concentration levels. The x-ray energy spectra of the high and low energy are estimated using an indirect transmission method. The derived monochromatic images are compared with the high- and low-energy CT images to demonstrate beam hardening artifacts reduction. Quantitative results were measured and compared to the true values. Results: The differences between the true density value used for simulation and that were obtained from the monochromatic images, are 1.8%, 1.3%, 2.3%, and 2.9% for the dose levels from standard dose to 1/8 dose, and are 0.4%, 0.7%, 1.5%, and 1.8% for the four iodine concentration levels from 6 mg/mL to 24 mg/mL. For all of the cases, beam hardening artifacts, especially streaks shown between dense inserts, are almost completely removed in the monochromatic images. Conclusion: The proposed algorithm provides an effective way to yield material images and artifacts-free monochromatic images at different dose levels without the need for phantom calibration or ROI selection. Furthermore, the approach also yields accurate results when the concentration of the iodine concentrate insert is very low, suggesting the algorithm is robust with respect to the low-contrast scenario.

60 APPLIED LIFE SCIENCES↗

Relativistic configuration-interaction and many-body-perturbation-theory calculations of U i hyperfine constants

Neutral uranium (U i) is a very difficult atom for theoretical calculations due to a large number (six) of valence electrons, strong valence-valence and valence-core correlations, high density of states, and relativistic effects. Configuration-interaction many-body perturbation theory (CI-MBPT) can efficiently treat valence-core correlations and relativistic effects, but because the formalism was developed for a Dirac-Hartree-Fock (DHF) starting potential that does not contain valence electrons, quite large CI space is needed to compensate for the + 6 charge of such a potential. Much more efficient is the relativistic configuration-interaction (RCI) approach, which uses a relatively accurate starting DHF potential that includes some valence electrons to make the valence-electron Hamiltonian diagonally dominated for some states. Here we report calculations of U i hyperfine constants of several low-energy states using the RCI method with the starting potential that includes four f valence electrons. With this starting potential, it is possible to use the single-configuration approximation or small basis sets to obtain quite accurate results for hyperfine-structure constants. In fact, by scaling the nuclear magnetic moment, the agreement for five levels was within 5% and a new magnetic moment can be recommended, 0.43(2). When two states are mixed, it is difficult to predict the exact mixing fractions, so in this case we proposed a method in which the mixing fraction is found by optimization of one parameter in the MBPT correction to obtain correct g factors. The resulting hyperfine constants are significantly improved. Finally, the CI-MBPT approach was also tested, with some success for low-energy states, with the limited basis set. The methods investigated here can be further developed to include more extensive data sets to improve accuracy and can be applied to other atoms and for calculations of other properties, for example, relevant to fundamental symmetry tests.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Particle Method for the Multispecies Landau Equation

Abstract The multispecies Landau collision operator describes the two-particle, small scattering angle or grazing collisions in a plasma made up of different species of particles such as electrons and ions. Recently, a structure preserving deterministic particle method (Carrillo et al. in J. Comput. Phys. 7:100066, 2020) has been developed for the single species spatially homogeneous Landau equation. This method relies on a regularization of the Landau collision operator so that an approximate solution, which is a linear combination of Dirac delta distributions, is well-defined. Based on a weak form of the regularized Landau equation, the time dependent locations of the Dirac delta functions satisfy a system of ordinary differential equations. In this work, we extend this particle method to the multispecies case, and examine its conservation of mass, momentum, and energy, and decay of entropy properties. We show that the equilibrium distribution of the regularized multispecies Landau equation is a Maxwellian distribution, and state a critical condition on the regularization parameters that guarantees a species independent equilibrium temperature. A convergence study comparing an exact multispecies Bobylev-Krook-Wu (BKW) solution to the particle solution shows approximately 2nd order accuracy. Important physical properties such as conservation, decay of entropy, and equilibrium distribution of the particle method are demonstrated with several numerical examples.

Mathematics↗

IMEX HDG-DG: A coupled implicit hybridized discontinuous Galerkin and explicit discontinuous Galerkin approach for shallow water systems

We propose IMEX HDG-DG schemes for planar and spherical shallow water systems. Of interest is subcritical flow, where the speed of the gravity wave is faster than that of nonlinear advection. In order to simulate these flows efficiently, we split the governing system into a stiff part describing the gravity wave and a non-stiff part associated with nonlinear advection. The former is discretized implicitly with the HDG method while an explicit Runge-Kutta DG discretization is employed for the latter. The proposed IMEX HDG-DG framework: 1) facilitates high-order solutions both in time and space; 2) avoids overly small time-step sizes; 3) requires only one linear system solve per time stage; 4) relative to DG generates smaller and sparser linear systems while promoting further parallelism; and 5) suppresses the fast modes in the system with a large time-step size. Numerical results for various test cases demonstrate that our methods are beneficial for applications where non-stiff terms are accurately treated while stiff terms are less accurately handled.

42 ENGINEERING↗

OECD/NEA EGMPEBV Activities in Multi-Physics Verification and Validation

The Expert Group on Multi-Physics Experimental Data, Benchmarks and Validation (EGMPEBV) was established in 2014 by the Nuclear Energy Agency (NEA) at the Organization for Economic Co-operation and Development (OECD) to establish processes for certification of experimental data and development of benchmark models for validation of Multi-Physics (MP) Modeling and Simulation (M and S) computational systems. The aim of the group is to provide member countries with guidelines and recommendations for validating and improving their novel multi-physics simulations, and with an access to key experimental data. To be used to their full potential, these novel tools will require a more complex suite of validation experiments due to the multiple energy, spatial and temporal scales, as well as the number of physical phenomena being simulated. They will also require accurate measurements of all terms, including coupling terms representing the feedback effects, in situations where the appropriate experimental techniques and facilities might not exist today. However, the ability to conduct validation experiments has for certain applications progressed at a slower pace than for computational methods, or in some cases it has significantly regressed (through shutdown of facilities, or through retirement of experts). In addition, as modelling capabilities are reaching deeper levels in single physics phenomena and in coupled behaviors, appropriate experimental techniques may not yet exist. Validation of multiple physics models requires a wide variety of experimental data, which emphasizes the importance of maximizing the use of historically accumulated data to avoid the significant cost of performing similar experiments today. Moreover, the pedigree and quality of this data as it pertains to modern day simulation methods may be insufficient for rigorous validation of high fidelity modeling and simulation tools. Thus, these data can often exhibit limitations including a lack of measurement accuracy or a lack of knowledge of the measurement uncertainty in the experimental techniques; an absence of measurements for some critical parameters; an obscuring of relationships between individual physics phenomena owing to the output data resolution; or simply an inability to allow for a detailed reinterpretation of the experiments based on the documented information. The organization of the EGMPEBV relies on three Task Forces, focused on 1) Experimental Data Qualification and Benchmark Evaluation; 2) Methods and Standards, and 3) Specific Applications. Task Force 1 is focused on providing better and more accurate experimental datasets, which support validation of high-fidelity multi-physics modelling and simulation (M and S) tools. At the same time it has to consider the incompleteness of past experimental data sets (in terms of data, documentation, or uncertainties) and with the limited number of available multi-physics experimental facilities. Task Force 2 is designed to develop validation methods and guidelines, and uncertainty qualification for the new family of multi-physics, multi-scale codes in the context of emerging demands such as longer fuel cycles and power uprate. Task Force 3 is focused on applications of validation experiments and has the objective to organize multi-physics M and S benchmarks involving validation experiments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Earth rotation information derived from MERIT and POLARIS VLBI observations

Advanced methods involving observations of extraterrestrial objects, such as artificial satellites or quasars, may make it possible to monitor geometrical variations in survey networks of regional, continental, and global scale with a spatial resolution of a few centimeters and a temporal resolution of better than one day. However, in connection with a realization of this potential, it is necessary to account for variations in the orientation of the earth in space. The accuracy provided by the conventional approaches for determining the orientation of the earth is not sufficient in this case. However, an appropriate method could be based on the utilization of independent clock radio interferometry (commonly referred to as VLBI) observations of extragalactic radio sources. The project POLARIS (Polar-motion Analysis by Radio Interferometric Sampling) is concerned with an implementation of this method. The MERIT (Monitor Earth Rotation and Intercompare the Techniques of observation and analysis) observations represent an aid to project POLARIS. Attention is given to details regarding these programs and the results obtained thus far.

Robertson, D. S.↗

Alternative Modal Basis Selection Procedures for Nonlinear Random Response Simulation

Three procedures to guide selection of an efficient modal basis in a nonlinear random response analysis are examined. One method is based only on proper orthogonal decomposition, while the other two additionally involve smooth orthogonal decomposition. Acoustic random response problems are employed to assess the performance of the three modal basis selection approaches. A thermally post-buckled beam exhibiting snap-through behavior, a shallowly curved arch in the auto-parametric response regime and a plate structure are used as numerical test articles. The results of the three reduced-order analyses are compared with the results of the computationally taxing simulation in the physical degrees of freedom. For the cases considered, all three methods are shown to produce modal bases resulting in accurate and computationally efficient reduced-order nonlinear simulations.

Przekop, Adam↗

Solving the Bateman equations in CASMO5 using implicit ode numerical methods for stiff systems

The Bateman equations, which describe the transmutation of nuclides over time as a result of radioactive decay, absorption, and fission, are often numerically stiff. This is especially true if short-lived nuclides are included in the system. This paper describes the use of implicit numerical methods for o D Es applied to the stiff Bateman equations, specifically employing the Backward Differentiation Formulas (BDF) form of the linear multistep method. As is true in other domains, using an implicit method removes or lessens the (sometimes severe) step-length constraints by which explicit methods must abide. To gauge its accuracy and speed, the BDF method is compared to a variety of other solution methods, including Runge-Kutta explicit methods and matrix exponential methods such as the Chebyshev Rational Approximation Method (CRAM). A preliminary test case was chosen as representative of a PWR lattice depletion step and was solved with numerical libraries called from a Python front-end. The Figure of Merit (a combined measure of accuracy and efficiency) for the BDF method was nearly identical to that for CRAM, while explicit methods and other matrix exponential approximations trailed behind. The test case includes 319 nuclides, in which the shortest-lived nuclide is {sup 98}Nb with a half-life of 2.86 seconds. Finally, the BDF and CRAM methods were compared within CASMO5, where CRAM had a FOM about four times better than BDF, although the BDF implementation was not fully optimized. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Alternative Modal Basis Selection Procedures For Reduced-Order Nonlinear Random Response Simulation

Three procedures to guide selection of an efficient modal basis in a nonlinear random response analysis are examined. One method is based only on proper orthogonal decomposition, while the other two additionally involve smooth orthogonal decomposition. Acoustic random response problems are employed to assess the performance of the three modal basis selection approaches. A thermally post-buckled beam exhibiting snap-through behavior, a shallowly curved arch in the auto-parametric response regime and a plate structure are used as numerical test articles. The results of a computationally taxing full-order analysis in physical degrees of freedom are taken as the benchmark for comparison with the results from the three reduced-order analyses. For the cases considered, all three methods are shown to produce modal bases resulting in accurate and computationally efficient reduced-order nonlinear simulations.

Przekop, Adam↗

SULI Technical Report: WIRE (Draft)

Local climate undoubtedly affects facility maintenance requirements; however, knowledge of the link between weather and facility condition over time is grounded in the 1990s and needs to be updated. WIRE aims to quantify the relationship between climate and infrastructure by calculating climate-based service life estimates to be used in BUILDER, a program the Department of Energy (DOE) uses for facility management. Phase One of this project is a case study on window-or-wall-mounted AC units. The independent variables influencing service life are load, maintenance quality, and effective age, with load calculated from three weather variables: temperature, humidity, and solar radiation. Given sufficient data, service life can be calculated as a function of these variables via regression. WIRE is currently collecting data. Preliminary estimated service life projections as a function of load derived from Lawrence Livermore National Laboratory (LLNL) data support the need for updated analysis on the relationship between climate and facility condition. With the groundwork laid, WIRE can generate deliverables immediately upon arrival of data. The methods used for this case study constitute a framework for deriving the relationship between weather variables and service life for other components.

42 ENGINEERING↗

Noise and drift analysis of non-equally spaced timing data

Generally, it is possible to obtain equally spaced timing data from oscillators. The measurement of the drifts and noises affecting oscillators is then performed by using a variance (Allan variance, modified Allan variance, or time variance) or a system of several variances (multivariance method). However, in some cases, several samples, or even several sets of samples, are missing. In the case of millisecond pulsar timing data, for instance, observations are quite irregularly spaced in time. Nevertheless, since some observations are very close together (one minute) and since the timing data sequence is very long (more than ten years), information on both short-term and long-term stability is available. Unfortunately, a direct variance analysis is not possible without interpolating missing data. Different interpolation algorithms (linear interpolation, cubic spline) are used to calculate variances in order to verify that they neither lose information nor add erroneous information. A comparison of the results of the different algorithms is given. Finally, the multivariance method was adapted to the measurement sequence of the millisecond pulsar timing data: the responses of each variance of the system are calculated for each type of noise and drift, with the same missing samples as in the pulsar timing sequence. An estimation of precision, dynamics, and separability of this method is given.

Vernotte, F.↗