Search NASA⌕ Search

SEARCH · Search NASA

Results for “analytic solution”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Providing plant data analytics through a seamless digital environment - 211

As technology continues to evolve and become more integrated into a worker's daily routine in the commercial nuclear power industry the need for easy access to data becomes a priority. Not only does the need for data increase but the amount of data collected increases. In most cases the data is collected and stored in various software applications, many of which are legacy systems, which do not offer any other option to access the data except through the application's user interface. Furthermore the data get grouped in 'silos' according to work function and not necessarily by subject. Hence, in order to access all the information needed for a particular task or analysis one may have to access multiple applications to gather all the data needed. The industry and the research community have identified the need for a digital architecture and more importantly the need for a Seamless Digital Environment (SDE). An SDE provides a means to access multiple applications, gather the data points needed, conduct the analysis requested, and present the result to the user with minimal or no effort by the user. The study researches the potential approaches to building an analytics solution for equipment reliability, on a small scale, focusing on either a single piece of equipment or a single system. The analytics solution will consist of a data integration layer, predictive and machine learning layer and the user interface layer that will display the output of the analysis in a straight forward, easy to consume manner. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1976

Approximate solutions of radiative transfer in dusty nebulae. II - Hydrogen and helium

The ionization structure of hydrogen and helium in dusty nebulae is discussed. General equations for pure hydrogen nebulae have been modified to include helium. An approximate analytical solution is presented for Stroemgren radii of hydrogen and helium in the absence of dust. It is shown that these results, with simple modifications, are also applicable to dusty nebulae where the effective absorption cross section of dust grains varies slowly with frequency in the range from 1000 to 200 A. No analytical solutions are possible if this cross section varies rapidly with frequency. In this case, however, simple coupled differential equations are derived which can easily be solved numerically. Approximate analytical expressions are given for evaluating the variation of the fraction of ionizing radiation absorbed by dust and the ratio of the volume emission measures of He II to H II regions with the spectrum of the ionizing source, helium abundance, and absorption properties of dust. The effects of dust on the He III zone are discussed. The present results are restricted to spherically symmetric nebulae, but nonuniform gas and dust distributions and clumpiness can be taken into account by the general results.

Dana, R. A.↗
DOE OSTI2025

Development, Verification, and Validation of an OpenFOAM-Based Solver for Modeling Inertial Fusion Energy Chambers

Our work seeks to introduce a computational tool tailored to the physics of inertial fusion energy chambers, in particular, those concepts based on thick liquid walls. In this approach, the structural materials are protected by several neutron mean-free-paths of renewable liquid and thus will be able to survive much longer than un-shielded walls, with virtually all structures lasting for the life of the plant and enabling the use of commercially available and qualified materials. The OpenFOAM-based solver named rhoCentralFoam has been used as a starting point. rhoCentralFoam belongs to the standard OpenFOAM solver toolset. It is a high-speed, explicit compressible flow solver with shock-capturing capability. While the main features have been retained, the solver had to be restructured to make use of tabular data for equations of states, a necessary addition to model the complex thermo-physical properties of ionized gasses. This entailed the need to change the independent state variables used by the solver, resulting in a new thermodynamic library and slightly different solution algorithm. Moreover, a radiation heat transfer model based on the P-1 approximation was added to the solver. The solver is verified against an analytical solution from the Sedov-Taylor-Neumann test problem to showcase the ability of the hydrodynamic solvers to handle strong shocks, whereas the P-1 model was verified using a simple one-dimensional problem with an analytical solution. Additionally, a validation case involving shock-wave propagation through jet array is presented, and the results are compared with experimental data from the open literature. Lastly, in order to showcase the utility of the solver for practical cases, we applied the refined solver to two representative scenarios: gas venting within the HYLIFE-II chamber and the compression of the gas following the partial ablation of the liquid wall.

Chamber dynamics↗
NASA NTRS2003

Finite Element Modeling of a Cylindrical Contact Using Hertzian Assumptions

The turbine blades in the high-pressure fuel turbopump/alternate turbopump (HPFTP/AT) are subjected to hot gases rapidly flowing around them. This flow excites vibrations in the blades. Naturally, one has to worry about resonance, so a damping device was added to dissipate some energy from the system. The foundation is now laid for a very complex problem. The damper is in contact with the blade, so now there are contact stresses (both normal and tangential) to contend with. Since these stresses can be very high, it is not all that difficult to yield the material. Friction is another non-linearity and the blade is made out of a Nickel-based single-crystal superalloy that is orthotropic. A few approaches exist to solve such a problem and computer models, using contact elements, have been built with friction, plasticity, etc. These models are quite cumbersome and require many hours to solve just one load case and material orientation. A simpler approach is required. Ideally, the model should be simplified so the analysis can be conducted faster. When working with contact problems determining the contact patch and the stresses in the material are the main concerns. Closed-form solutions for non-conforming bodies, developed by Hertz, made out of isotropic materials are readily available. More involved solutions for 3-D cases using different materials are also available. The question is this: can Hertzian1 solutions be applied, or superimposed, to more complicated problems-like those involving anisotropic materials? That is the point of the investigation here. If these results agree with the more complicated computer models, then the analytical solutions can be used in lieu of the numerical solutions that take a very long time to process. As time goes on, the analytical solution will eventually have to include things like friction and plasticity. The models in this report use no contact elements and are essentially an applied load problem using Hertzian assumptions to determine the contact patch dimensions.

Knudsen, Erik↗
DOE OSTI2016

Guided Formation of 3D Helical Mesostructures by Mechanical Buckling: Analytical Modeling and Experimental Validation

3D helical mesostructures are attractive for applications in a broad range of microsystem technologies due to their mechanical and electromagnetic properties as stretchable interconnects, radio frequency antennas, and others. Controlled compressive buckling of 2D serpentine-shaped ribbons provides a strategy to formation of such structures in wide ranging classes of materials (from soft polymers to brittle inorganic semiconductors) and length scales (from nanometer to centimeter), with an ability for automated, parallel assembly over large areas. The underlying relations between the helical configurations and fabrication parameters require a relevant theory as the basis of design for practical applications. In this work, an analytic model of compressive buckling in serpentine microstructures is presented based on the minimization of total strain energy that results from various forms of spatially dependent deformations. Experiments at micro- and millimeter scales, together with finite element analyses, have been exploited to examine the validity of developed model. The theoretical analyses shed light on general scaling laws in terms of three groups of fabrication parameters (related to loading, material, and 2D geometry), including a negligible effect of material parameters and a square root dependence of primary displacements on the compressive strain. Furthermore, analytic solutions were obtained for the key physical quantities (e.g., displacement, curvature and maximum strain). A demonstrative example illustrates how to leverage the analytic solutions in choosing the various design parameters, such that brittle fracture or plastic yield can be avoided in the assembly process.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Validation of chemistry models employed in a particle simulation method

The chemistry models employed in a statistical particle simulation method, as implemented in the Intel iPSC/860 multiprocessor computer, are validated and applied. Chemical relaxation of five-species air in these reservoirs involves 34 simultaneous dissociation, recombination, and atomic-exchange reactions. The reaction rates employed in the analytic solutions are obtained from Arrhenius experimental correlations as functions of temperature for adiabatic gas reservoirs in thermal equilibrium. Favorable agreement with the analytic solutions validates the simulation when applied to relaxation of O2 toward equilibrium in reservoirs dominated by dissociation and recombination, respectively, and when applied to relaxation of air in the temperature range 5000 to 30,000 K. A flow of O2 over a circular cylinder at high Mach number is simulated to demonstrate application of the method to multidimensional reactive flows.

Haas, Brian L.↗
DOE OSTI2017

Diffusion Synthetic Acceleration for High Order S{sub N} Transport on Meshes with Curved Surfaces

Thermal radiation transport is integral to multiphysics applications such as inertial confinement fusion and astrophysics. Materials in these problems can be exceptionally optically thick and diffusive. A transport solver must be capable of accurate and efficient solutions to such problems. Previously, Woods et al. established a proof-of-concept for high order (HO) discontinuous Galerkin finite element (DGFEM) transport on meshes with curved surfaces. They numerically demonstrated the spatial convergence is order p + 1. Remarkably, these convergence rates are relatively invariant to the amount of curvature of the spatial mesh. We present here an extension of this work that demonstrates a rapidly-convergent iterative solver for the FEM transport equations on meshes with curved sides. We use the modified interior penalty (MIP) diffusion synthetic acceleration (DSA) established by Wang and Ragusa. The partially consistent MIP equations stem from the symmetric interior penalty (IP) derivation of the diffusion equation. The IP method is not stable for optically thick cells so it was combined with the diffusion conforming form (DCF), a discretized diffusion equation derived from the discretized S N equations that is not stable for intermediate optical thicknesses. The MIP equations adapt a 'switch' to the IP method in optically thinner regions and the DCF method in optically thicker regions. The literature does not reveal any Fourier analyses performed for iterative transport on meshes with curved surfaces. Wang and Ragusa performed a Fourier analysis for orthogonal meshes and our results appear to be as expected. Previous literature has reported on some of the sensitivities that affect the spectral radius of DSA schemes: polynomial order, and scattering ratios. Presently, we focus on the sensitivity of the spectral radius to polynomial order, scattering ratio, the constant C, and mesh curvature for varying cell thicknesses. We also numerically study the behavior of the spectral radius to verify our DSA implementation. A sequence of problems is designed to exercise the acceleration. As the problems become more optically thick and diffusive, the accelerated solution continues to converge toward the analytic solution. We employ the Modular Finite Element Methods (MFEM) package, developed at Lawrence Livermore National Laboratory, to generate the linear algebraic system of equations for each angle in the discrete ordinates quadrature set, and use UMFPack to solve this system of equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1989

The evolution of energetic particles and the emitted radiation in solar flares

The evolution of accelerated particle distributions in a magnetized plasma and the resulting radiation are calculated, and the results are applied to solar flares. To study the radiation on timescales of order the particle lifetimes, the evolution of the particle distribution is determined by the use of the Fokker-Planck equation including Coulomb collisions and magnetic mirroring. Analytic solution to the equations are obtained for limiting cases such as homogeneous injection in a homogeneous plasma, and for small pitch angle. These analytic solutions are then used to place constraints on flare parameters such as density, loop length, and the injection timescale for very short implusive solar flares. For general particle distributions in arbitrary magnetic field and background density, the equation is solved numerically. The relative timing of microwaves and X-rays during individual flares is investigated. A number of possible sources for excessive microwave flux are discussed including a flattening in the electron spectrum above hard X-ray energies, thermal synchrotron emission, and trapping of electron by converging magnetic fields. Over shorter timescales, the Fokker-Planck equation is solved numerically to calculate the temporal evolution of microwaves and X-rays from nonthermal thick target models. It is shown that magnetic trapping will not account for the observed correlation of microwaves of approximately 0.15 seconds behind X-rays in flares with rapid time variation, and thus higher energy electrons must be accelerated later than lower energy electrons.

Lu, Edward Tsang↗
NASA NTRS1989

Interactions between long and synoptic-scale waves. I - Instability of a nonzonal flow. II - Growth rate of long waves

The interactions between synoptic-scale and long waves were investigated analytically. First, the influence of the long wave on the synoptic-scale wave is examined. The structure of a synoptic-scale wave growing on a wavy basic scale was analyzed under an assumption that the synoptic-scale waves have the structure of the most unstable normal modes. The derived analytical solution, which is simple and is amenable to physical interpretation, can be understood in terms of eddies and their local growth rate. The analytical solution is then used to examine growth of a long wave in the presence of parameterized synoptic-scale waves. Two possibly unstable solutions were found; one is a modification of the linear long wave, and the other a strongly nonlinear solution. In both cases, the synoptic-scale wave increases the growth rate of the long wave.

Ebisuzaki, Wesley↗
DOE OSTI2017

Analytical general solutions for static wormholes in f ( R , T ) gravity

Originally proposed as a tool for teaching the general theory of relativity, wormholes are today approached in many different ways and are seeing as an efficient alternative for interstellar and time travel. Attempts to achieve observational signatures of wormholes have been growing as the subject has become more and more popular. In this article we investigate some f ( R , T ) theoretical predictions for static wormholes, i.e., wormholes whose throat radius can be considered a constant. Since the T -dependence in f ( R , T ) gravity is due to the consideration of quantum effects, a further investigation of wormholes in such a theory is well motivated. We obtain the energy conditions of static wormholes in f ( R , T ) gravity and apply an analytical approach to find their physical and geometrical solutions. We highlight that our results are in agreement with previous solutions and assumptions presented in the literature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Numerical Study on Crossflow Printed Circuit Heat Exchanger for Advanced Small Modular Reactors

Various fluids such as water, gases (helium), molten salts (FLiNaK, FLiBe) and liquid metal (sodium) are used as a coolant of advanced small modular reactors (SMRs). The printed circuit heat exchanger (PCHE) has been adopted as the intermediate and/or secondary heat exchanger of SMR systems because this heat exchanger is compact and effective. The size and cost of PCHE can be changed by the coolant type of each SMR. In this study, the crossflow PCHE analysis code for advanced small modular reactor has been developed for the thermal design and cost estimation of the heat exchanger. The analytical solution of single pass, both unmixed fluids crossflow heat exchanger model was employed to calculate a two dimensional temperature profile of a crossflow PCHE. The analytical solution of crossflow heat exchanger was simply implemented by using built in function of the MATLAB program. The effect of fluid property uncertainty on the calculation results was evaluated. In addition, the effect of heat transfer correlations on the calculated temperature profile was analyzed by taking into account possible combinations of primary and secondary coolants in the SMR systems. Size and cost of heat exchanger were evaluated for the given temperature requirement of each SMR.

42 ENGINEERING↗
DOE OSTI2015

Efficient simulation and model reformulation of two-dimensional electrochemical thermal behavior of lithium-ion batteries

Lithium-ion batteries are an important technology to facilitate efficient energy storage and enable a shift from petroleum based energy to more environmentally benign sources. Such systems can be utilized most efficiently if good understanding of performance can be achieved for a range of operating conditions. Mathematical models can be useful to predict battery behavior to allow for optimization of design and control. An analytical solution is ideally preferred to solve the equations of a mathematical model, as it eliminates the error that arises when using numerical techniques and is usually computationally cheap. An analytical solution provides insight into the behavior of the system and also explicitly shows the effects of different parameters on the behavior. However, most engineering models, including the majority of battery models, cannot be solved analytically due to non-linearities in the equations and state dependent transport and kinetic parameters. The numerical method used to solve the system of equations describing a battery operation can have a significant impact on the computational cost of the simulation. In this paper, a model reformulation of the porous electrode pseudo three dimensional (P3D) which significantly reduces the computational cost of lithium ion battery simulation, while maintaining high accuracy, is discussed. This reformulation enables the use of the P3D model into applications that would otherwise be too computationally expensive to justify its use, such as online control, optimization, and parameter estimation. Furthermore, the P3D model has proven to be robust enough to allow for the inclusion of additional physical phenomena as understanding improves. In this study, the reformulated model is used to allow for more complicated physical phenomena to be considered for study, including thermal effects.

25 ENERGY STORAGE↗
DOE OSTI2021

Analytical modeling for redox flow battery design

Deeper market penetration of redox flow batteries requires optimization of the cell performance. Though important for optimization, detailed analytical solutions have not been developed for electrolyte flow, mass and charge transport, and reaction kinetics within redox flow batteries. To this end, here we present analytical solutions to active species concentration and over-potential based on advection-diffusion transport for ions and Bulter-Volmer model for interface reaction kinetics. The solutions were validated with results from a finite element model. These solutions were then applied to investigate the relationship between over-potential and state of charge, current density, reaction rate constant, flow velocity, diffusivity, total active species concentration, and electrode structure. Explicit formulas were identified for minimum activation over-potential and limiting current density as well as their dependence on electrolyte properties, operation conditions, and electrode structure. With our new mathematical formulas, this work provides a theoretical framework for flow battery design.

25 ENERGY STORAGE↗
DOE OSTI2017

Spontaneous imbibition of water and determination of effective contact angles in the Eagle Ford Shale Formation using neutron imaging

Understanding of fundamental processes and prediction of optimal parameters during the horizontal drilling and hydraulic fracturing process results in economically effective improvement of oil and natural gas extraction. Although, the modern analytical and computational models can capture fracture growth, there is a lack of experimental data on spontaneous imbibition and wettability in oil and gas reservoirs for the validation of further model development. In this work, we used neutron imaging to measure the spontaneous imbibition of water into fractures of Eagle Ford Shale with known geometries and fracture orientations. An analytical solution for a set of nonlinear second-order differential equations was applied to the measured imbibition data to determine effective contact angles. The analytical solution fit the measured imbibition data reasonably well and determined effective contact angles were slightly higher than static contact angles due to effects of in-situ changes in velocity, surface roughness, and heterogeneity of mineral surfaces on the fracture surface. Additionally, small fracture widths may have retarded imbibition and affected model fits, which suggests that average fracture widths are not satisfactory for modeling imbibition in natural systems.

DiStefano, Victoria H. (ORCID:0000000160239716)↗
DOE OSTI2020

Multipole-based cable braid electromagnetic penetration model: magnetic penetration case

The goal of this paper is to present, for the first time, calculations of the magnetic penetration case of a first principles multipole-based cable braid electromagnetic penetration model. As a first test case, a one-dimensional array of perfect electrically conducting wires, for which an analytical solution is known, is investigated: we compare both the self-inductance and the transfer inductance results from our first principles cable braid electromagnetic penetration model to those obtained using the analytical solution. These results are found in good agreement up to a radius to half spacing ratio of about 0.78, demonstrating a robustness needed for many commercial and non-commercial cables. We then analyze a second set of test cases of a square array of wires whose solution is the same as the one-dimensional array result and of a rhomboidal array whose solution can be estimated from Kley's model. As a final test case, we consider two layers of one-dimensional arrays of wires to investigate porpoising effects analytically. We find good agreement with analytical and Kley's results for these geometries, verifying our proposed multipole model. Note that only our multipole model accounts for the full dependence on the actual cable geometry which enables us to model more complicated cable geometries.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Comprehensive Studies on Hot Compaction and Vibration-Assisted Compaction Tests of Aluminum Powder

Abstract Aluminum powder compaction was studied using both test and simulation. Cold compaction, hot compaction, and vibration-assisted (cold) compaction tests were conducted to achieve different density ratios. First, the hot compaction test (at 300 °C, compression pressure 140 MPa) improved about 6% compared with cold compaction under the same compression pressure. Second, although the relative density ratio does not obviously improve at a vibration-assisted (cold) compaction, the strength of the specimens made under vibration loading is much better than those of cold compaction. Additionally, finite element models with well-calibrated Drucker–Prager Cap (DPC) material constitutive model were built in abaqus/standard to simulate the powder compaction process. The results of the finite element model have very good correlations with test results up to the tested range, and this finite element model further predicts the loading conditions needed to achieve the higher density ratios. Two exponential equations of the predicted density ratio were obtained by combining the test data and the simulation results. A new analytical solution was developed to predict the axial pressure versus the density ratio for the powder compaction according to DPC material model. The results between the analytical solution and the simulation model have a very good match.

Engineering↗
NASA NTRS1983

On the off-axis tension test for unidirectional composites

The off axis tensile test was examined experimentally to obtain actual displacement fields over the surface of graphite polyimide coupon specimens. The experimental results were compared with approximate analytical solutions and generated finite element results. An optical method of high sensitivity moireinterferometry was used to determine the actual displacements to high precision. The approximate analytical solution and the finite element results compare very favorably with the measured centerline displacements in the test section, and the finite element displacement fields provide excellent agreement with the moiredisplacements throughout the specimen. A 15 degree fiber orientation and coupon aspect ratios of 5 and 15 are presented. Previously announced in STAR as N82-23553

Herakovich, C. T.↗
NASA NTRS1996

A Model for Axial Magnetic Bearings Including Eddy Currents

This paper presents an analytical method of modelling eddy currents inside axial bearings. The problem is solved by dividing an axial bearing into elementary geometric forms, solving the Maxwell equations for these simplified geometries, defining boundary conditions and combining the geometries. The final result is an analytical solution for the flux, from which the impedance and the force of an axial bearing can be derived. Several impedance measurements have shown that the analytical solution can fit the measured data with a precision of approximately 5%.

Kucera, Ladislav↗