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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 595 records · Page 33

All sources NASA NTRS DOE OSTI

DOE OSTI2024

SAM Code Enhancements for Modeling of Liquid Metal-Cooled Fast Reactor Concepts

The SAM code is under development and supported by DOE-NE’s Nuclear Energy Advanced Modeling and Simulation program as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. These advanced reactor concepts incorporate novel and improved approaches to achieve safety and economic feasibility. This report summarizes two major efforts in addressing the modeling gaps in SAM for liquid-metal-cooled fast reactor (LMFR) applications, i.e. thermal mixing and stratification phenomena in large pools and corrosion-oxidation of components in flowing lead. A new one-dimensional model for thermal mixing and stratification effects in large pools and enclosures is developed and implemented. Thermal mixing and stratification occur when fluid enters a pool at a temperature different than the bulk fluid itself, a scenario often encountered during transients in pool-type systems. These phenomena are critical for the safety of reactors, impacting phenomena like natural circulation, which is essential for passive cooling. The improved model in SAM addresses limitations of state-of-the-art approaches by combining one-dimensional (1D) channels, representing the coolant jet flow, with lumped-parameter zero-dimensional (0D) pools, representing the rest of coolant in the tank. Energy exchange between the 1D jet and the 0D pools is based on heat transfer correlations calibrated against 3D simulations. It is verified that this model can handle various flow configurations, including hot jets in colder pools, cold jets in hotter pools, and the presence of features like ceilings, free surfaces, and obstacles. Additionally, validation against experimental data demonstrates the ability of the model to capture mixing and stratification effects in a wide range of conditions. The flexibility and improved accuracy of the new model make it a valuable tool for reactor safety analysis, allowing for the simulation of different geometries encountered in advanced reactors. A system-level corrosion modeling capability is developed and implemented in SAM to support Lead Fast Reactor (LFR) development. Although the initial focus of this capability will be on LFR application, this can later be expanded to include other liquid metals such as Lead-Bismuth Eutectic (LBE) and PbLi. This report summarizes the common corrosion mitigation strategies and outlines the progress on implementing and validating a corrosion-oxidation model in SAM. Verification and validation of the corrosion-oxidation portion of the model was performed using analytical solution and measured data from samples tested in the non-isothermal pumped lead loop at IPPE Obninsk. The iron transport and corrosion/precipitation portion of the model was assessed using an analytical model and measured corrosion depths from a natural convection lead loop experiment performed at CEA. It is demonstrated that the model implemented in SAM performed well in these assessments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS1974

Three-dimensional elastic stress and displacement analysis of finite geometry solids containing cracks

The line method of analysis is applied to the Navier-Cauchy equations of elastic equilibrium to calculate the displacement distributions in various bodies containing cracks. The application of this method to these equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. When decoupling the equations and their boundary conditions is not possible, the use of a successive approximation procedure permits the analytical solution of the resulting ordinary differential equations. The results obtained show a considerable potential for using this method in the three-dimensional analysis of finite geometry solids and suggest a possible extension of this technique to nonlinear material behavior.

Gyekenyesi, J. P.↗

NASA NTRS1975

Three-dimensional elastic stress and displacement analysis of finite geometry solids containing cracks

Gyekenyesi, J. P.↗

DOE OSTI2023

Generalized analytical and numerical modeling of optical second harmonic generation in anisotropic crystals and complex heterostructures using #SHAARP package

Optical second harmonic generation (SHG) is a nonlinear optical effect widely used for nonlinear optical microscopy and laser frequency conversion. The closed-form analytical solution of the nonlinear optical responses is essential for evaluating the optical responses of new materials whose optical properties are unknown a priori. Many approximations have therefore been employed in the existing analytical approaches, such as slowly varying approximation, weak reflection of the nonlinear polarization, transparent medium, high crystallographic symmetry, Kleinman symmetry, easy crystal orientation along a high-symmetry direction, phase matching conditions and negligible interference among nonlinear waves, which may lead to large errors in the reported material properties. To avoid these approximations, here we have developed an open-source package named Second Harmonic Analysis of Anisotropic Rotational Polarimetry (#SHAARP) for single interface (si) and in multilayers (ml) for homogeneous crystals. The reliability and accuracy are established by experimentally benchmarking with both the SHG polarimetry and Maker fringes predicted from the package using standard materials. SHAARP.si and SHAARP.ml are available through GitHub https://github.com/Rui-Zu/SHAARP and https://github.com/bzw133/SHAARP.ml, respectively.

complex systems↗

NASA NTRS1990

The method of lines in analyzing solids containing cracks

A semi-numerical method is reviewed for solving a set of coupled partial differential equations subject to mixed and possibly coupled boundary conditions. The line method of analysis is applied to the Navier-Cauchy equations of elastic and elastoplastic equilibrium to calculate the displacement distributions in various, simple geometry bodies containing cracks. The application of this method to the appropriate field equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. When decoupling of the equations and their boundary conditions is not possible, the use of a successive approximation procedure permits the analytical solution of the resulting ordinary differential equations. The use of this method is illustrated by reviewing and presenting selected solutions of mixed boundary value problems in three dimensional fracture mechanics. These solutions are of great importance in fracture toughness testing, where accurate stress and displacement distributions are required for the calculation of certain fracture parameters. Computations obtained for typical flawed specimens include that for elastic as well as elastoplastic response. Problems in both Cartesian and cylindrical coordinate systems are included. Results are summarized for a finite geometry rectangular bar with a central through-the-thickness or rectangular surface crack under remote uniaxial tension. In addition, stress and displacement distributions are reviewed for finite circular bars with embedded penny-shaped cracks, and rods with external annular or ring cracks under opening mode tension. The results obtained show that the method of lines presents a systematic approach to the solution of some three-dimensional mechanics problems with arbitrary boundary conditions. The advantage of this method over other numerical solutions is that good results are obtained even from the use of a relatively coarse grid.

Gyekenyesi, John P.↗

NASA NTRS1975

Thermomechanical coupling in fatigue fracture of viscoelastic materials

Subcritical crack growth in a linear viscoelastic material subjected to cyclic loading is investigated starting with the thermodynamic power balance. Physically, it is supposed that the subcritical slow crack growth is due to local weakening of the material in the neighborhood of the crack tip. Under fatigue loading it is assumed that all energy dissipation goes into heat and that this heat build-up is the dominant irreversible process governing crack growth. The cycle averaged temperature distribution around the crack tip is obtained from local application of the first law of thermodynamics (i.e., conservation of energy). The analysis leads to the solution of a nonlinear integro-differential equation for crack length as a function of time which is coupled with the local energy equation. A regular perturbation technique is used to obtain an analytical solution which compares very well with experimental results.

Hufferd, W. L.↗

NASA NTRS1974

Finite elasto-plastic deformation. I - Theory and numerical examples

It is demonstrated that the problem of elasto-plastic finite deformation is governed by a quasi-linear model irrespective of deformation magnitude. This feature follows from the adoption of a rate viewpoint toward finite deformation analysis in an Eulerian reference frame. Objectivity of the formulation is preserved by introduction of a frame-invariant stress rate. Equations for piece-wise linear incremental finite element analysis are developed by application of the Galerkin method to the instantaneously linear governing differential equations of the quasi-linear model. Numerical solution capability has been established for problems of plane strain and plane stress. The accuracy of the numerical analysis is demonstrated by consideration of a number of problems of homogeneous finite deformation admitting comparative analytic solution. It is shown that accurate, objective numerical solutions can be obtained for problems involving dimensional changes of an order of magnitude and rotations of a full radian.

Osias, J. R.↗

NASA NTRS1976

Numerical modeling of pollutant transport using a Lagrangian marker particle technique

A derivation and code were developed for the three-dimensional mass transport equation, using a particle-in-cell solution technique, to solve coastal zone waste discharge problems where particles are a major component of the waste. Improvements in the particle movement techniques are suggested and typical examples illustrated. Preliminary model comparisons with analytic solutions for an instantaneous point release in a uniform flow show good results in resolving the waste motion. The findings to date indicate that this computational model will provide a useful technique to study the motion of sediment, dredged spoils, and other particulate waste commonly deposited in coastal waters.

Spaulding, M.↗

NASA NTRS1953

Analytical Derivation and Experimental Evaluation of Short-Bearing Approximation for Full Journal Bearing

An approximate analytical solution including the effect of end leakage from the oil film of short plain bearings is presented because of the importance of endwise flow in sleeve bearings of the short lengths commonly used. The analytical approximation is supported by experimental data, resulting in charts which facilitate analysis of short plain bearings. The analytical approximation includes the endwise flow and that part of the circumferential flow which is related to surface velocity and film thickness but neglects the effect of film pressure on the circumferential flow. In practical use, this approximation applies best to bearings having a length-diameter ratio up to 1, and the effects of elastic deflection, inlet oil pressure, and changes of clearance with temperature minimize the relative importance of the neglected term. The analytical approximation was found to be an extension of a little-known pressure-distribution function originally proposed by Michell and Cardullo.

Dubois, George B↗

NASA NTRS1995

Combining Comparison Functions and Finite Element Approximations in CFD

In a variety of potential flow applications, the modal element method has been shown to significantly reduce the numerical grid, employ a more precise grid termination boundary condition, and give theoretical insight to the flow physics. The method employs eigenfunctions to replace the numerical grid over significant portions of the flow field. Generally, a numerical grid is employed around obstacles with complex geometry while eigenfunctions are applied to regions in the flow field where the boundary conditions can easily be satisfied. To handle a wider class of computational fluid dynamics (CFD) problems, the present paper extends the modal element to include function approximations which do not satisfy the governing differential equation. To accomplish this task, a double modal series approximation and weighted residual constraints are developed to force the comparison functions to satisfy the governing differential equation and to interface properly with the finite element solution. As an example, the method is applied to the problem of potential flow in a channel with two-dimensional cylindrical like obstacles. The calculated flow fields are in excellent agreement with exact analytical solutions.

Baumeister, Kenneth J.↗

NASA NTRS1997

Rotor Vortex Filaments: Living on the Slipstream's Edge

The purpose of this paper is to gain a better understanding of rotor wake evolution in hover and axial flow by deriving an analytical solution for the time dependent behavior of vortex filament circulation and core size. This solution is applicable only for vortex filaments in the rotor far-wake. A primarily inviscid vortex/shear layer interaction (where the slipstream boundary is modeled as a shear layer) has been identified in this analytical treatment. This vortex/shear layer interaction results in decreasing, vortex filament circulation and core size with time. The inviscid vortex/shear layer interaction is shown, in a first-order treatment, to be of greater magnitude than viscous diffusion effects. The rate of contraction, and ultimate collapse, of the vortex filament core is found to be directly proportional to the rotor inflow velocity. This new insight into vortex filament decay promises to help reconcile several disparate observations made in the literature and will, hopefully, promote new advances in theoretical modeling of rotor wakes.

Young, Larry A.↗

NASA NTRS2005

Capillary Driven Flows Along Differentially Wetted Interior Corners

Closed-form analytic solutions useful for the design of capillary flows in a variety of containers possessing interior corners were recently collected and reviewed. Low-g drop tower and aircraft experiments performed at NASA to date show excellent agreement between theory and experiment for perfectly wetting fluids. The analytical expressions are general in terms of contact angle, but do not account for variations in contact angle between the various surfaces within the system. Such conditions may be desirable for capillary containment or to compute the behavior of capillary corner flows in containers consisting of different materials with widely varying wetting characteristics. A simple coordinate rotation is employed to recast the governing system of equations for flows in containers with interior corners with differing contact angles on the faces of the corner. The result is that a large number of capillary driven corner flows may be predicted with only slightly modified geometric functions dependent on corner angle and the two (or more) contact angles of the system. A numerical solution is employed to verify the new problem formulation. The benchmarked computations support the use of the existing theoretical approach to geometries with variable wettability. Simple experiments to confirm the theoretical findings are recommended. Favorable agreement between such experiments and the present theory may argue well for the extension of the analytic results to predict fluid performance in future large length scale capillary fluid systems for spacecraft as well as for small scale capillary systems on Earth.

Golliher, Eric L.↗

DOE OSTI2021

Risk-based area of review estimation in overpressured reservoirs to support injection well storage facility permit requirements for CO 2 storage projects

This paper by the Energy & Environmental Research Center presents a workflow and modeling approach for delineating a risk-based area of review (AOR) to support a U.S. Environmental Protection Agency (EPA) Class VI permit for a carbon dioxide (CO 2 ) storage project. The approach combines semianalytical solutions for estimating formation fluid leakage through a hypothetical leaky wellbore with the results of numerical reservoir simulations to define the AOR. The modeling utilizes 1) semianalytical solutions from the peer-reviewed literature for formation fluid leakage through abandoned wellbores by Raven (1990) and Avci (1994), 2) a FORTRAN model compiled and described in Cihan et al. (2011, 2012) called ASLMA (Analytical Solution for Leakage in Multilayered Aquifers), and 3) a computational framework for estimating a risk-based AOR first proposed by Oldenburg et al. (2014, 2016). Therefore, the approach builds upon well-established research and underlying hydrogeological principles that have been upheld for nearly three decades. Moreover, the ASLMA model has been broadly applied to an array of storage projects. The work presented herein extends these earlier works using a custom wrapper written in the software environment, R (R Core Team, 2020), which was developed to perform multiple runs of the ASLMA model using given ranges for one or more input parameters. In addition, the current work simulates the pressure buildup within the storage reservoir in response to CO 2 injection using a compositional simulator to better accommodate the temporospatial evolution of pressure buildup within the storage reservoir that is more accurately modeled using a heterogeneous geologic model and a compositional simulator that accounts for the multiphase interactions. The workflow is demonstrated using a case study for a 180,000-metric-ton-per-year storage project located in the PCOR (Plains CO 2 Reduction) Partnership region. For the storage project evaluated here, under the scenario where the leaky wellbore is open to a saline aquifer (thief zone) between the overlying seal (cap rock) and the underground sources of drinking water (USDW), the risk-based AOR essentially collapses to the areal extent of the CO 2 plume in the storage reservoir because the pressure buildup in the storage reservoir beyond the CO 2 plume is insufficient to drive formation fluids up a hypothetical leaky wellbore into the USDW. However, even under the conservative assumption that the leaky wellbore is not open to a thief zone, beyond the areal extent of the CO 2 plume, the incremental leakage is less than 400 m 3 over 20 years, which represents ~0.0001% or less of the total volume of water contained within the USDW rock volume. As discussed in the text, the threshold criterion for defining the risk-based AOR is site-specific and should be informed by the results of the sensitivity analysis and available site characterization data. The approach outlined in this paper is designed to be protective of USDWs and, therefore, comply with the Safe Drinking Water Act requirements and provisions for the U.S. EPA Class VI Underground Injection Control (UIC) Program (Class VI Rule) and North Dakota Administrative Code Chapter 43-05-01.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1973

Finite deformation of elasto-plastic solids

A theoretical basis is established for analysis of finite deformation of metals. The observation that finite deformation of such elastoplastic materials may be viewed as a process rather than an event leads to derivation of a complete initial and boundary value problem distinguished by its quasilinear nature. This feature of the formulation motivates adoption of an incremental approach to numerical problem solving. Numerical solution capability is established for problems of plane stress and plane strain. The validity of the theory and numerical analysis is demonstrated by consideration of a number of problems of homogeneous finite deformation for which analytic solutions are available. Subsequently the analysis is employed for the investigation of necking in flat metal tensile bars. The results of this investigation provide the first full numerical solutions for tensile necking in plane stress and plane strain. In addition a basis is provided for assessment of the validity of stress-strain relations inferred from tensile test data.

Osias, J. R.↗

NASA NTRS1984

Implementation of uniform perturbation method for potential flow past axisymmetric and two-dimensional bodies

The aerodynamic characteristics of potential flow past an axisymmetric slender body and a thin airfoil are calculated using a uniform perturbation analysis method. The method is based on the superposition of potentials of point singularities distributed inside the body. The strength distribution satisfies a linear integral equation by enforcing the flow tangency condition on the surface of the body. The complete uniform asymptotic expansion of its solution is obtained with respect to the slenderness ratio by modifying and adapting an existing technique. Results calculated by the perturbation analysis method are compared with the existing surface singularity panel method and some available analytical solutions for a number of cases under identical conditions. From these comparisons, it is found that the perturbation analysis method can provide quite accurate results for bodies with small slenderness ratio. The present method is much simpler and requires less memory and computation time than existing surface singularity panel methods of comparable accuracy.

Wong, T. C.↗

NASA NTRS1998

Magnetohydrodynamic Origin of Jets from Accretion Disks

A review is made of magnetohydrodynamic (MHD) theory and simulation of outflows from disks for different distributions of magnetic field threading the disk. In one limit of a relatively weak, initially diverging magnetic field, both thermal and magnetic pressure gradients act to drive matter to an outflow, while a toroidal magnetic field develops which strongly collimates the outflow. The collimation greatly reduces the field divergence and the mass outflow rate decreases after an initial peak. In a second limit of a strong magnetic field, the initial field configuration was taken with the field strength on the disk decreasing outwards to small values so that collimation was reduced. As a result, a family of stationary solutions was discovered where matter is driven mainly by the strong magnetic pressure gradient force. The collimation in this case depends on the pressure of an external medium. These flows are qualitatively similar to the analytic solutions for magnetically driven outflows. The problem of the opening of a closed field line configuration linking a magnetized star and an accretion disk is also discussed.

Lovelace, R. V. E.↗

DOE OSTI2026

Propagation of partially spatially coherent laser beams in instantaneous Kerr media

The propagation of intense, partially spatially coherent laser beams in a medium with instantaneous third-order susceptibility is studied analytically and numerically. For sufficiently high power relative to that required for nonlinear self-focusing, the propagation initially proceeds in two stages. In the first stage, spatial coherence builds up, and in the second stage, the number of speckles reduces. Once the degree of coherence is sufficiently high, whole-beam self-focusing occurs. The beam power is mostly confined within the initial spot radius. Two analytical approaches for describing the evolution of the beam are presented. The method of moments leads to an analytical solution for the rms spot radius that is in excellent agreement with simulations. This method does not require any knowledge of the field statistics beyond the initial conditions and provides no information about the evolution of the individual speckles. The other approach employs a self-similar solution for the second-order coherence function of the field and assumes that the fourth-order coherence function is factorizable and obeys complex circular Gaussian random statistics. The latter method also leads to an analytical expression for the spot radius, but its predictions for the qualitative evolution of the speckles disagree with wave-optics simulations.

lasers↗

NASA NTRS1987

Reynolds stress closure in jet flows using wave models

Numerical methods were developed that will form the computational part of the turbulence closure scheme. A wave model was developed for the two-dimensional shear layer. This configuration is being used as a test case for the closure schemes. Various numerical schemes were examined to give efficient solutions of the Rayleigh equation for this geometry. These include both spectral and finite difference methods. Secondly, numerical methods are under development to solve the non-separable Rayleigh equation. This solution is required for the closure scheme in more complex geometries. A model problem was used to assist in the algorithm development. Two-dimensional spectral methods and a hybrid spectral/finite difference technique were developed. An analytic solution of the Rayleigh equation for a basic elliptic flow was obtained. This will be used to verify the stability codes developed for arbitrary geometries. Other numerical methods for solving the Rayleigh equation based on the boundary element technique were also examined. These solutions are forming the basis of a model for the shock structure in jets of arbitrary geometry.

Morris, P. J.↗