Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical 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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

i-SPin: an integrator for multicomponent Schrödinger-Poisson systems with self-interactions

Abstract We provide an algorithm and a publicly available code to numerically evolve multicomponent Schrödinger-Poisson (SP) systems with a SO( n ) symmetry, including attractive or repulsive self-interactions in addition to gravity. Focusing on the case where the SP system represents the non-relativistic limit of a massive vector field, non-gravitational self-interactions (in particular spin-spin interactions) introduce complexities related to mass and spin conservation which are not present in purely gravitational systems. We address them with an analytical solution for the `kick' step in the algorithm, where we are able to decouple the multicomponent system completely. Equipped with this analytical solution, the full field evolution is second order accurate, preserves spin and mass to machine precision, and is reversible. Our algorithm allows for an expanding universe relevant for cosmology, and the inclusion of external potentials relevant for laboratory settings.

Astronomy & Astrophysics↗
NASA NTRS1971

Further studies on the dimensionless parameters associated with the 'in vivo' transport of heat within biological tissue.

Dimensionless parameters appearing in the analysis of transport of heat within living biological tissues are presented and discussed. Dimensional analysis as well as some analytical solutions have been employed to obtain the physically and physiologically most significant dimensionless parameters associated with these types of problems. Some analytical solutions are presented for both steady state and transient cases to illustrate the significance of the derived parameters.-

Shitzer, A.↗
NASA NTRS1991

The free compressible viscous vortex

The present study investigates the effects of compressibility on free (unsteady) viscous heat-conducting vortices. Analytical solutions are found in the limit of large but finite Reynolds number and small but finite Mach number. It is shown that the spreading of the vortex causes a radial flow. This flow is given by the solution of an ordinary differential equation, which gives the dependence of the radial velocity on the tangential velocity, density, and temperature profiles of the vortex. Estimates of the radial velocity found by solving this equation are found to be in good agreement with numerical solutions of the full equations. The equations for the viscous evolution are expanded in powers of Mach number to obtain detailed analytical solutions. It is shown that swirling axisymmetric compressible flows generate negative radial velocities far from the vortex core owing to viscous effects, regardless of the initial distributions of vorticity, density, and entropy.

Colonius, Tim↗
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↗
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.↗
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 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 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↗
DOE OSTI2021

Discrete Green’s functions and spectral graph theory for computationally efficient thermal modeling

Here, this work concerns solutions of the heat equation with the spectral graph method, for which the temperature is defined at discrete points in the domain and the spatial relationship among the points is described by a graph. The heat equation on the graph is solved using matrix techniques involving the eigenvectors and eigenvalues of the Laplacian matrix. The spectral graph approach precludes the computationally intensive meshing and numerous time-integration steps of the finite element method. In the present work, the spectral graph method is extended to include heat loss at the boundaries with a generalized boundary condition, and physics-based edge weights are introduced which simplify the calibration process. From this approach a discrete Green’s function is defined which allows for solutions under a variety of heating conditions including: space-varying initial conditions; time-and-space varying internal heating; and, time-and-space-varying heating at boundaries of type 1 (Dirichlet), type 2 (Neumann) and type 3 (Robin). Results are provided for benchmark heat transfer problems in one spatial dimension and in three spatial dimensions, and verification is provided by comparison with exact analytical solutions and finite difference solutions. The spectral graph method converges within 0.4% error of the analytical solution. The practical utility of the approach is demonstrated by thermal simulation of a multilayer additive manufacturing process. The spectral graph results are compared to experimentally-obtained temperature data for two metal parts, with error less than 5% of the experimental measurements, with computation time less than one minute on a desktop computer.

36 MATERIALS SCIENCE↗
DOE OSTI2022

CO 2 zonal injection rate allocation and plume extent evaluation through wellbore temperature analysis

Temperature analysis during a pause in injection operations, known as warmback analysis, has been used in the petroleum industry for evaluating the injection conformance and estimating the location of the flooded front in applications, such as waterflooding oil reservoirs. Here in this work, methods are introduced to extend the application of temperature warmback analysis to estimate the zonal CO 2 injection rate and zonal CO 2 plume extent during geologic CO 2 storage in a saline aquifer. First, novel analytical solutions are developed to model transient temperature in the aquifer during the injection and subsequent shut-in periods considering two-phase flow (gaseous CO 2 and aqueous brine) conditions in the aquifer. The solution involves a discretization of the aquifer into regions; the energy and mass conservation equations for the regions are solved simultaneously considering appropriate boundary conditions at the interfaces. Two solutions techniques are presented: multi-region and three-region solutions. Inverse models are developed accordingly to evaluate the injection profile and estimate the extent of the plume front in the reservoir during the injection period. The multi-region solution results in an inversion approach that requires regression analysis. However, the three-region formulation results in a simple graphical technique for inverse modeling. The analytical solutions are validated against a thermally coupled reservoir simulation tool using different synthetic cases for CO 2 injection in deep saline aquifers. The results of the developed solutions provide a good match with numerical results during forward and inverse modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Model Package Report: Composite Analysis Solid Waste Release Model (CASWR Model)

This document describes the implementation of a solid waste form release model in GoldSim for the Hanford Site Composite Analysis (CA) Update. This Composite Analysis Solid Waste Release model (CASWR model) was designed to generate deterministic radionuclide release rates for Hanford’s Central Plateau solid waste disposal sites using single realizations of release model coefficients. Five generalized waste form types are used for the conceptual model of waste release: surplus reactor block, cement, soil-debris, grouted residual waste, and ancillary equipment. The surplus reactor block waste form consists of radionuclide leaching from surplus production graphite reactor core blocks via unspecified processes (White et al., 1984 as cited in PNNL-15965). The cement waste form represents solidified wastes whose permeability is much lower than that of the surrounding soil. The soil-debris waste form type is defined as unconsolidated waste mixed with soil material. Tanks and canyon complexes comprise the grouted residual waste form such that their solid waste will be grouted and capped with a surface barrier at the completion of their cleanup. Finally, the ancillary equipment waste form constitutes contaminant releases from ancillary and auxiliary waste form residues associated with tank farms at closure. Individual sub-models are implemented to numerically represent a respective waste form within the CASWR Model: Surplus Reactor Block Sub-model, Cement Sub-model, SoilDebris Sub-model, Grouted Residual Waste Sub-model, and Ancillary Equipment Submodel. Advection is assumed to be the primary transport process governing the release of radionuclides in the Ancillary Equipment and Soil-Debris Sub-models. Diffusion is assumed to be the primary release process in the Grouted Residual Waste and Cement Sub-models. An unspecified zero-order release process is considered in the Surplus Reactor Block Sub-model due to the lack of information regarding actual processes involved in irradiated graphite leaching. The Surplus Reactor Block, Cement, and SoilDebris Sub-models were compared against analytical solutions (PNNL-11800, Composite Analysis for Low-Level Waste Disposal in the 200 Area Plateau of the Hanford Site). The agreement between results of these analytical solutions and the corresponding waste form models verified their correct implementation in GoldSim. A 1-D modeling abstraction approach for the Grouted Residual Waste and Ancillary Equipment Sub-models was adopted from existing Performance Assessment (PA) models. Despite the simplifications made in these sub-models, they were found to be appropriate representations of the waste forms, similar to what was used in the Waste Management Area C PA model (RPP-ENV-58782, Performance Assessment of Waste Management Area C, Hanford Site, Washington, Rev. 0). A sensitivity analysis was conducted to identify the most influential parameters in each waste form sub-model. Suggestions for considering pH-dependent and redox-dependent release mechanisms are formulated through the development of a conditional constant approach in GoldSim.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1982

On the off-axis tensile 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 moire interferometry 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 moire displacements throughout the specimen. A 15 degree fiber orientation and coupon aspect ratios of 5 and 15 are presented.

Nemeth, M. P.↗
NASA NTRS1973

Preliminary numerical analysis of improved gas chromatograph model

A mathematical model for the gas chromatograph was developed which incorporates the heretofore neglected transport mechanisms of intraparticle diffusion and rates of adsorption. Because a closed-form analytical solution to the model does not appear realizable, techniques for the numerical solution of the model equations are being investigated. Criteria were developed for using a finite terminal boundary condition in place of an infinite boundary condition used in analytical solution techniques. The class of weighted residual methods known as orthogonal collocation is presently being investigated and appears promising.

Woodrow, P. T.↗