Results for “numerical calculations: variational” — Page 9

DOE OSTI2023

AGR-5/6/7 Thermal Model with Non-uniform Gas Gaps

Fuel compact temperatures are a crucial factor in assessing the irradiation performance of tri-structural isotropic fuel particles. In the absence of direct measurement, fuel compact temperatures were calculated using a three-dimensional finite element thermal model, which is subject to simulation uncertainty. The most dominant factor in the uncertainty of calculated fuel temperatures is the gas gap uncertainty due to the nub-to-shell clearance caused by a design error of AGR-5/6/7 capsules. The thermal model was revised to examine the most probable graphite offset position for six different days during the irradiation for Capsules 1 and 2. The analysis varied the offset distance and azimuthal direction at both the top and bottom of the holder. The best-fit offset was estimated based on the minimum root mean square error of the residuals (measured minus calculated) for the operational thermocouples (TCs). From these results, the following conclusions were made: (1) The holder offsets led to slightly lower average temperatures but wider temperature variations (lower minimum and higher peak fuel temperatures) for both Capsule 1 and Capsule 2. (2) During earlier cycles (162A–164B), when numerous TCs were still operational, the best-fit offset distance varied over a specific range for both the top and bottom ([0.002–0.0035 in.] for Capsule 1 and [0.003-0.004 in] for Capsule 2). In contrast, the offset azimuthal direction varied widely, especially for the offset at the bottom of the Capsule 1 holder. This is because holder movement was somewhat constrained at the top of Capsule 1 by the TC leads running through the capsule head and into the holder and by the through tubes in Capsule 2, but the Capsule 1 bottom did not have this type of constraint. (3) During later cycles, when all TCs failed, applying the maximum possible offset of 0.006 in. to the northwest direction for both the top and bottom resulted in a calculated peak fuel temperature of 1557? in Capsule 1 (i.e., a 135? increase from 1422? with zero offset on September 20, 2019 (166A)); the maximum offset of 0.0068 in. to the south for both top and bottom resulted in a calculated peak fuel temperature of 1110°C in Capsule 2 on April 20, 2020 (i.e., a 116? increase from 994? with zero offset (168A)). High peak fuel temperatures in Capsule 1 during Cycle 166A could be the cause of massive particle failure near the end of this cycle. (4) Even though the highest temperature at the tip of Type-N TCs, such as TC-1-7, slightly exceeded 1000?, the temperature along the TC wire reached as high as 1335? assuming an offset of 0.006 in. in the northwest of Capsule 1 holder near the end of Cycle 166A. This temperature significantly exceeds the temperature threshold at which TC degradation is expected to occur, ultimately contributing to considerable particle failures in Capsule 1. For eight Type-N TCs in Capsule 2, the peak TC line temperature was much lower (i.e., 1029°C for TC-2-5), assuming maximum offset of 00068 in. to the south during cycle 168A.

advanced gas reactors↗

DOE OSTI2023

AGR-5/6/7 Thermal Model with Non-uniform Gas Gaps

Fuel compact temperatures are a crucial factor in assessing the irradiation performance of tri-structural isotropic fuel particles. In the absence of direct measurement, fuel compact temperatures were calculated using a three-dimensional finite element thermal model, which is subject to simulation uncertainty. The most dominant factor in the uncertainty of calculated fuel temperatures is the gas gap uncertainty due to the nub-to-shell clearance caused by a design error of AGR-5/6/7 capsules. The thermal model was revised to examine the most probable graphite offset position for six different days during the irradiation for Capsules 1 and 2. The analysis varied the offset distance and azimuthal direction at both the top and bottom of the holder. The best-fit offset was estimated based on the minimum root mean square error of the residuals (measured minus calculated) for the operational thermocouples (TCs). From these results, the following conclusions were made: (1) The holder offsets led to slightly lower average temperatures but wider temperature variations (lower minimum and higher peak fuel temperatures) for both Capsule 1 and Capsule 2. (2) During earlier cycles (162A–164B), when numerous TCs were still operational, the best-fit offset distance varied over a specific range for both the top and bottom ([0.002–0.0035 in.] for Capsule 1 and [0.003-0.004 in] for Capsule 2). In contrast, the offset azimuthal direction varied widely, especially for the offset at the bottom of the Capsule 1 holder. This is because holder movement was somewhat constrained at the top of Capsule 1 by the TC leads running through the capsule head and into the holder and by the through tubes in Capsule 2, but the Capsule 1 bottom did not have this type of constraint. (3) During later cycles, when all TCs failed, applying the maximum possible offset of 0.006 in. to the northwest direction for both the top and bottom resulted in a calculated peak fuel temperature of 1557? in Capsule 1 (i.e., a 135? increase from 1422? with zero offset on September 20, 2019 (166A)); the maximum offset of 0.0068 in. to the south for both top and bottom resulted in a calculated peak fuel temperature of 1110°C in Capsule 2 on April 20, 2020 (i.e., a 116? increase from 994? with zero offset (168A)). High peak fuel temperatures in Capsule 1 during Cycle 166A could be the cause of massive particle failure near the end of this cycle. (4) Even though the highest temperature at the tip of Type-N TCs, such as TC-1-7, slightly exceeded 1000?, the temperature along the TC wire reached as high as 1335? assuming an offset of 0.006 in. in the northwest of Capsule 1 holder near the end of Cycle 166A. This temperature significantly exceeds the temperature threshold at which TC degradation is expected to occur, ultimately contributing to considerable particle failures in Capsule 1. For eight Type-N TCs in Capsule 2, the peak TC line temperature was much lower (i.e., 1029°C for TC-2-5), assuming maximum offset of 00068 in. to the south during cycle 168A.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

AGR-5/6/7 Thermal Model with Non-uniform Gas Gaps

Fuel compact temperatures are a crucial factor in assessing the irradiation performance of tri-structural isotropic fuel particles. In the absence of direct measurement, fuel compact temperatures were calculated using a three-dimensional finite element thermal model, which is subject to simulation uncertainty. The most dominant factor in the uncertainty of calculated fuel temperatures is the gas gap uncertainty due to the nub-to-shell clearance caused by a design error of AGR-5/6/7 capsules. The thermal model was revised to examine the most probable graphite offset position for six different days during the irradiation for Capsules 1 and 2. The analysis varied the offset distance and azimuthal direction at both the top and bottom of the holder. The best-fit offset was estimated based on the minimum root mean square error of the residuals (measured minus calculated) for the operational thermocouples (TCs). From these results, the following conclusions were made: (1) The holder offsets led to slightly lower average temperatures but wider temperature variations (lower minimum and higher peak fuel temperatures) for both Capsule 1 and Capsule 2. (2) During earlier cycles (162A–164B), when numerous TCs were still operational, the best-fit offset distance varied over a specific range for both the top and bottom ([0.002–0.0035 in.] for Capsule 1 and [0.003-0.004 in] for Capsule 2). In contrast, the offset azimuthal direction varied widely, especially for the offset at the bottom of the Capsule 1 holder. This is because holder movement was somewhat constrained at the top of Capsule 1 by the TC leads running through the capsule head and into the holder and by the through tubes in Capsule 2, but the Capsule 1 bottom did not have this type of constraint. (3) During later cycles, when all TCs failed, applying the maximum possible offset of 0.006 in. to the northwest direction for both the top and bottom resulted in a calculated peak fuel temperature of 1557? in Capsule 1 (i.e., a 135? increase from 1422? with zero offset on September 20, 2019 (166A)); the maximum offset of 0.0068 in. to the south for both top and bottom resulted in a calculated peak fuel temperature of 1110°C in Capsule 2 on April 20, 2020 (i.e., a 116? increase from 994? with zero offset (168A)). High peak fuel temperatures in Capsule 1 during Cycle 166A could be the cause of massive particle failure near the end of this cycle. (4) Even though the highest temperature at the tip of Type-N TCs, such as TC-1-7, slightly exceeded 1000?, the temperature along the TC wire reached as high as 1335? assuming an offset of 0.006 in. in the northwest of Capsule 1 holder near the end of Cycle 166A. This temperature significantly exceeds the temperature threshold at which TC degradation is expected to occur, ultimately contributing to considerable particle failures in Capsule 1. For eight Type-N TCs in Capsule 2, the peak TC line temperature was much lower (i.e., 1029°C for TC-2-5), assuming maximum offset of 00068 in. to the south during cycle 168A.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1974

Optimum propellant usage for reaction jet systems of space vehicles

The on-off type control for reaction jet systems is proven to be the optimal fuel scheme. However, due to the nonlinear characteristics of this type control, no direct method of solution is known for this optimal process being applied to the attitude control of space vehicles. This paper will discuss the optimization of fuel usage for an attitude control of space vehicles. A computation technique is developed for the calculation of optimal control law. The method of calculus of variations is applied to the estimation of the changes of performance index as well as terminal constraints. Thus an algorithm is obtained by the steepest descent method. A numerical example is given in the paper.

Liu, T. C.↗

NASA NTRS1984

Dynamic analysis of evolutive conservative systems. Discussion of eigenmode crossings

After an analysis of the close connection between the symmetries of a dynamical system and the multiplicity of its vibrational natural frequencies, it is proved by variational arguments that for a system of invariable symmetry the eigenfrequencies associated with the eigenmodes of a given symmetry type do not cross, in general, during the evolution of this system. The theory is implemented by some numerical calculations applied to the analysis of the evolution of the axisymmetric hydroelastic modes of the Ariane launch vehicle during burning of the first stage.

Morand, H. J. P.↗

NASA NTRS1994

Polarimetric method of estimation of vertical aerosol distribution in application to observations of ozone and NO2

The degree of polarization of skylight at the zenith during twilight depends on the aerosol content in the atmosphere. The long-term observations at the high-mountain research station 'Kislovodsk' (North Caucasus) have shown that the variation of the degree of polarization after the eruption of the El Chichon volcano can serve as the effective parameter characterizing the vertical aerosol stratification in the atmosphere. The results of the measurements are confirmed by the numerical calculations. The algorithm of the retrieval of the vertical aerosol distribution on the base of the measurements of the degree of polarization is proposed. This method can be applied for the increasing of the precision of O3, NO2, and other gas content measurements.

Elansky, Nikolay F.↗

NASA NTRS1995

Effect of mesh distortion on the accuracy of transverse shear stresses and their sensitivity coefficients in multilayered composites

A study is made of the effect of mesh distortion on the accuracy of transverse shear stresses and their first-order and second-order sensitivity coefficients in multilayered composite panels subjected to mechanical and thermal loads. The panels are discretized by using a two-field degenerate solid element, with the fundamental unknowns consisting of both displacement and strain components, and the displacement components having a linear variation throughout the thickness of the laminate. A two-step computational procedure is used for evaluating the transverse shear stresses. In the first step, the in-plane stresses in the different layers are calculated at the numerical quadrature points for each element. In the second step, the transverse shear stresses are evaluated by using piecewise integration, in the thickness direction, of the three-dimensional equilibrium equations. The same procedure is used for evaluating the sensitivity coefficients of transverse shear stresses. Numerical results are presented showing no noticeable degradation in the accuracy of the in-plane stresses and their sensitivity coefficients with mesh distortion. However, such degradation is observed for the transverse shear stresses and their sensitivity coefficients. The standard of comparison is taken to be the exact solution of the three-dimensional thermoelasticity equations of the panel.

Noor, Ahmed K.↗

NASA NTRS2005

Experimental Study on Interactions Between H Atoms and Organic Haze

In Titan s atmosphere composed of N2 and CH4, irradiations of both solar ultraviolet light and charged particles induce active chemical reactions. In the processes of these reactions, a large amount of hydrogen (H) atoms are expected to be formed by dissociation of CH4 and other hydrocarbons [e.g., 1, 2]. Theoretical models suggest that these active H atoms need to be converted to stable hydrogen molecules (H2) efficiently to maintain unsaturated hydrocarbons and organic haze in Titan s atmosphere [e.g., 1]. Furthermore, molecular hydrogen is an important greenhouse effect gas in Titan s atmosphere, and small variation in its abundance strongly affects Titan s surface temperature [3]. Thus, the formation of H2 molecules from H atoms is a key reaction for both the atmospheric chemistry and the surface environment of Titan. Although several numerical calculations have been conducted to investigate the atmospheric chemistry of Titan with hypothesized recombination reactions of H atoms, such as catalytic scheme of C4H2 [e.g., 1, 2], it is still unclear what chemical reaction is responsible for the conversion of H atoms to H2 molecules in Titan s atmosphere.

Sekine, Y.↗

NASA NTRS1996

Multigrid Acceleration of Time-Accurate Navier-Stokes Calculations

A numerical scheme to solve the unsteady Navier-Stokes equations is described. The scheme is fully implicit in time and is unconditionally stable (at least for first- and second-order discretizations of the physical time derivatives). With unconditional stability, the choice of the time step is based on the physical phenomena to be resolved rather than limited by numerical stability. This is especially important for high Reynolds number viscous flows, where the spatial variation of grid cell size can be as much as six orders of magnitude. A multigrid-multiblock, steady-state, three-dimensional Navier-Stokes solver, TLNS3D, was modified to iteratively invert the equations at each physical time step. The implementation of this procedure in TLNS3D is discussed. The implications of applying several popular turbulence models to unsteady flow are also considered. Numerical results are presented to show the application of the scheme to various two-dimensional turbulent flows. The results of a three-dimensional laminar flow calculation are also given.

Melson, N. Duane↗

NASA NTRS2016

The Variations of Neutron Component of Lunar Radiation Background from LEND LRO Observations

Lunar neutron flux data measured by the Lunar Exploration Neutron Detector (LEND) on board NASA's Lunar Reconnaissance Orbiter (LRO) were analyzed for the period 2009-2014.We have re-evaluated the instrument's collimation capability and re-estimated the neutron counting rate measured in the Field of View (FOV) of the LEND collimated detectors, and found it to be 1.070.1counts per second. We derived the spectral density of the neutron flux for various lunar regions using our comprehensive numerical model of orbital measurements. This model takes into account the location of the LEND instrument onboard LRO to calculate the surface leakage neutron flux and its propagation to the instrument detectors. Based on this we have determined the lunar neutron flux at the surface to be approx. 2 neutrons/ [sq cm/ sec] in the epithermal energy range, 0.4e V to 1keV. We have also found variations of the lunar neutron leakage flux with amplitude as large as a factor of two, by using multi-year observations to explore variations in the Galactic Cosmic Ray (GCR) flux during the 23rd-24th solar cycles.

Litvak, M. L.↗

NASA NTRS1983

Aircraft aerodynamic prediction method for V/STOL transition including flow separation

A numerical procedure was developed for the aerodynamic force and moment analysis of V/STOL aircraft operating in the transition regime between hover and conventional forward flight. The trajectories, cross sectional area variations, and mass entrainment rates of the jets are calculated by the Adler-Baron Jet-in-Crossflow Program. The inviscid effects of the interaction between the jets and airframe on the aerodynamic properties are determined by use of the MCAIR 3-D Subsonic properties are determined by use of the MCAIR 3-D Subsonic Potential Flow Program, a surface panel method. In addition, the MCAIR 3-D Geometry influence Coefficient Program is used to calculate a matrix of partial derivatives that represent the rate of change of the inviscid aerodynamic properties with respect to arbitrary changes in the effective wing shape.

Gilmer, B. R.↗

NASA NTRS1987

Application of turbulence modeling to predict surface heat transfer in stagnation flow region of circular cylinder

A theoretical analysis and numerical calculations for the turbulent flow field and for the effect of free-stream turbulence on the surface heat transfer rate of a stagnation flow are presented. The emphasis is on the modeling of turbulence and its augmentation of surface heat transfer rate. The flow field considered is the region near the forward stagnation point of a circular cylinder in a uniform turbulent mean flow. The free stream is steady and incompressible with a Reynolds number of the order of 10 to the 5th power and turbulence intensity of less than 5 percent. For this analysis, the flow field is divided into three regions: (1) a uniform free-stream region where the turbulence is homogeneous and isotropic; (2) an external viscid flow region where the turbulence is distorted by the variation of the mean flow velocity; and, (3) an anisotropic turbulent boundary layer region over the cylinder surface. The turbulence modeling techniques used are the kappa-epsilon two-equation model in the external flow region and the time-averaged turbulence transport equation in the boundary layer region. The turbulence double correlations, the mean velocity, and the mean temperature within the boundary layer are solved numerically from the transport equations. The surface heat transfer rate is calculated as functions of the free-stream turbulence longitudinal microlength scale, the turbulence intensity, and the Reynolds number.

Wang, Chi R.↗

NASA NTRS1982

Solar irradiance modulation by active regions from 1969 through 1980

The solar irradiance variations resulting from sunspot deficits and facular excesses in emission have been calculated from 1969 through 1980. Agreement appears to exist between the calculations and the major features seen with the Nimbus 7 cavity pyrheliometer and with both the major and minor features detected by the Solar Maximum Mission ACRIM experiment. The 12-year irradiance variations calculated suggest a larger variance with increased solar activity, and little change in the average irradiance with solar activity. The largest excursions over these 12 years show a 0.4% variation. Removal of the activity influences upon solar irradiance during the numerous rocket experiments observing the solar constant may allow a better value for this quantity to be determined.

Schatten, K. H.↗

NASA NTRS1989

Numerical model of solar dynamic radiator for parametric analysis

Growth power requirements for Space Station Freedom will be met through addition of 25 kW solar dynamic (SD) power modules. Extensive thermal and power cycle modeling capabilities have been developed which are powerful tools in Station design and analysis, but which prove cumbersome and costly for simple component preliminary design studies. In order to aid in refining the SD radiator to the mature design stage, a simple and flexible numerical model was developed. The model simulates heat transfer and fluid flow performance of the radiator and calculates area mass and impact survivability for many combinations of flow tube and panel configurations, fluid and material properties, and environmental and cycle variations.

Rhatigan, Jennifer L.↗

NASA NTRS1992

Chaotic evolution of the solar system

The evolution of the entire planetary system has been numerically integrated for a time span of nearly 100 million years. This calculation confirms that the evolution of the solar system as a whole is chaotic, with a time scale of exponential divergence of about 4 million years. Additional numerical experiments indicate that the Jovian planet subsystem is chaotic, although some small variations in the model can yield quasi-periodic motion. The motion of Pluto is independently and robustly chaotic.

Sussman, Gerald J.↗

NASA NTRS1984

Second- and third-order upwind difference schemes for hyperbolic conservation laws

Second- and third-order two time-level five-point explicit upwind-difference schemes are described for the numerical solution of hyperbolic systems of conservation laws and applied to the Euler equations of inviscid gas dynamics. Nonliner smoothing techniques are used to make the schemes total variation diminishing. In the method both hyperbolicity and conservation properties of the hyperbolic conservation laws are combined in a very natural way by introducing a normalized Jacobian matrix of the hyperbolic system. Entropy satisfying shock transition operators which are consistent with the upwind differencing are locally introduced when transonic shock transition is detected. Schemes thus constructed are suitable for shockcapturing calculations. The stability and the global order of accuracy of the proposed schemes are examined. Numerical experiments for the inviscid Burgers equation and the compressible Euler equations in one and two space dimensions involving various situations of aerodynamic interest are included and compared.

Yang, J. Y.↗

NASA NTRS1972

A method for calculating laminar viscous, compressible flows with small pressure gradients

A method is presented for calculating laminar viscous, compressible flows in which pressure gradients are assumed to be neglibible in energy conservation and in developing property relations, but not in momentum conservation. The equations are developed from Navier-Stokes equations in a noninertial reference frame. The method is an extension of an approximate method developed for calculating combined forced and contained natural convection in a rotating tank. The restriction to small density and temperature variations in the previous method is removed. For illustration, a two-dimensional simulation of the rotating tank problem is formulated, and the numerical computation procedure to be used is indicated.

Martin, E. D.↗

NASA NTRS1980

An adaptive computation mesh for the solution of singular perturbation problems

In singular perturbation problems, control of zone size variation can affect the effort required to obtain accurate, numerical solutions of finite difference equations. The mesh is generated by the solution of potential equations. Numerical results for a singular perturbation problem in two dimensions are presented. The mesh was used in calculations of resistive magnetohydrodynamic flow in two dimensions.

Brackbill, J. U.↗