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.

SourceResearch connectionFrom yearThrough year

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

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS1996

Electron Flow to a Satellite at High Positive Potential

The Tethered Satellite System (TSS) is designed to deploy a 1.6 m diameter spherical satellite a distance of 20 km above the space shuttle orbiter on an insulated conducting tether. Because of the passage of the conducting tether through the earth's magnetic field, an emf is generated producing a positive satellite potential of about 5000 V. Electron flow under the influence of this high positive potential is the focus of the present analysis. The ionospheric parameters at TSS orbit altitude are; thermal velocity of electrons, 1.9 x 10(exp 5) M/S, thermal velocity of the ions, 1.1 x 10(exp 3) m/s, velocity of the satellite 8 x 10(exp 3) m/s. The electrons, with a Debye length, lambda(D) = 0.49 cm, spiral about the earth's magnetic field lines (0.4 Gauss) with a radius of about 3 cm and the ions spiral with a radius of 5 m. Under these conditions, the electron thermal energy, kT is 0.17 eV. The TSS satellite radius, r(p) is 163 Debye lengths. There is an extensive literature on the interaction of satellites with the near-earth ionospheric plasma. The space charge limitation to the electron current collected by a sphere at positive electrical potential was calculated by Langmuir and Blodgett (1924). Parker and Murphy (1967) recognized the importance of the influence of the earth's magnetic field and used the guiding center approximation to calculate the electron current collected by a positive charged satellite. More recently Ma and Schunk (1989) have calculated the time dependent flow of electrons to a spherical satellite at positive potential utilizing numerical methods and Sheldon (1994) used similar methods to solve this problem for the steady state. In order to analyze some of the phenomena that occurred in the ionosphere during the TSS flights, it would be useful to have analytic expressions for these electron flows. The governing equations are very complex and an exact analytical solution is not likely. An approximate analytical solution is feasible however, and the results of one attempt are presented herein.

Sheldon, John W.↗

NASA NTRS2008

Guidelines and Parameter Selection for the Simulation of Progressive Delamination

Turon s methodology for determining optimal analysis parameters for the simulation of progressive delamination is reviewed. Recommended procedures for determining analysis parameters for efficient delamination growth predictions using the Abaqus/Standard cohesive element and relatively coarse meshes are provided for single and mixed-mode loading. The Abaqus cohesive element, COH3D8, and a user-defined cohesive element are used to develop finite element models of the double cantilever beam specimen, the end-notched flexure specimen, and the mixed-mode bending specimen to simulate progressive delamination growth in Mode I, Mode II, and mixed-mode fracture, respectively. The predicted responses are compared with their analytical solutions. The results show that for single-mode fracture, the predicted responses obtained with the Abaqus cohesive element correlate well with the analytical solutions. For mixed-mode fracture, it was found that the response predicted using COH3D8 elements depends on the damage evolution criterion that is used. The energy-based criterion overpredicts the peak loads and load-deflection response. The results predicted using a tabulated form of the BK criterion correlate well with the analytical solution and with the results predicted with the user-written element.

Song, Kyongchan↗

DOE OSTI2012

Improvement of three-field based safety analysis code, SPACE, through verification and validation

To verify and validate the state-of-the-art safety analysis code, SPACE, features of three field governing equations were arranged and the comparative analysis with both analytic solution and experimental benchmarks were performed. With the world wide streamline on the safety analysis, SPACE uses two-fluid, three-field governing equations which aim to implement multi-dimensional and multi-scaled analysis. Through the Verification and Validation (V and V), developing code based on the three-field governing equation have not only been confirmed but also improved. V and V activities on SPACE have been conducted as to interfacial drag model in vertical turbulent flow. The verification of the interfacial drag model was based on cell size sensitivity test for confirming feasible calculation range of cell size. Furthermore, through the validation with analytic solution, the modification of the interfacial drag model in vertical turbulent annulus flow was carried out and the nominated model was incorporated. Modified interfacial drag model, Asali's correlation, yields improved result than previous interfacial drag model. The validation with experimental benchmarks had also been performed on the nearly horizontal counter current limitation. To confirm the possibility of the further application of the CCFL model in case with horizontal pipe, validation with nearly horizontal pipe experiment was carried out. V and V activities of SPACE with analytic solution and experimental benchmarks allow us to confirm/identify the tailorability for safety analysis. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2014

Asymptotic solution for heat convection-radiation equation

In this paper, we employ a new approximate analytical method called the optimal homotopy asymptotic method (OHAM) to solve steady state heat transfer problem in slabs. The heat transfer problem is modeled using nonlinear two-point boundary value problem. Using OHAM, we obtained the approximate analytical solution for dimensionless temperature with different values of a parameter ε. Further, the OHAM results for dimensionless temperature have been presented graphically and in tabular form. Comparison has been provided with existing results from the use of homotopy perturbation method, perturbation method and numerical method. For numerical results, we used Runge-Kutta Fehlberg fourth-fifth order method. It was found that OHAM produces better approximate analytical solutions than those which are obtained by homotopy perturbation and perturbation methods, in the sense of closer agreement with results obtained from the use of Runge-Kutta Fehlberg fourth-fifth order method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

A Forward Time Stepping Heat Conduction Model for Spot Melt Additive Manufacturing

Solidification dynamics are crucial for determining microstructure development in additively manufactured parts. Multiphysics models based on finite element or finite volume methods may help gain insight for complicated phenomena such as fluid flow, keyholing, and porosity but are too computationally expensive to use for simulating actual builds. Recent analytic and semi-analytic solutions for moving heat sources in a semi-infinite three-dimensional space provide a way to accurately estimate the solidification conditions for entire builds. The downside to these methods is that, unlike finite element or finite volume methods, they cannot use the temperature distribution of the previous timesteps to march the solution forward in time. This paper provides the mathematical formulation and implementation of a forward time stepping (FTS) approach to an existing semi-analytic solution. The speed and accuracy of the two methods are then compared for various scan patterns. The result is that, for spot melts, the forward time-stepping model provides improvements in both speed and accuracy. This is especially true for longer simulations, where the simulation can be orders of magnitude faster. Finally, the longest simulation analyzed in this paper was roughly 30× faster when using the forward time-stepping model versus the straightforward implementation of the semi-analytic solution.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Viscous Gubser flow with conserved charges to benchmark fluid simulations

We present semi-analytical solutions for the evolution of both the temperature and chemical potentials for viscous Gubser flow with conserved charges. Such a solution can be especially useful in testing numerical codes intended to simulate relativistic fluids with large chemical potentials. The freeze-out hypersurface profiles for constant energy density are calculated, along with the corresponding normal vectors, and presented as a new unit test for numerical codes. We also compare the influence of the equation of state on the semi-analytical solutions. We benchmark the newly developed smoothed particle hydrodynamics code ccake that includes both shear viscosity and three conserved charges. Here, the numerical solutions are in excellent agreement with the semi-analytical solution and are also able to accurately reproduce the hypersurface at freeze-out.

Hydrodynamic models↗

DOE OSTI2017

Sedimentary Geothermal Feasibility Study: October 2016

The objective of this project is to analyze the feasibility of commercial geothermal projects using numerical reservoir simulation, considering a sedimentary reservoir with low permeability that requires productivity enhancement. A commercial thermal reservoir simulator (STARS, from Computer Modeling Group, CMG) is used in this work for numerical modeling. In the first stage of this project (FY14), a hypothetical numerical reservoir model was developed, and validated against an analytical solution. The following model parameters were considered to obtain an acceptable match between the numerical and analytical solutions: grid block size, time step and reservoir areal dimensions; the latter related to boundary effects on the numerical solution. Systematic model runs showed that insufficient grid sizing generates numerical dispersion that causes the numerical model to underestimate the thermal breakthrough time compared to the analytic model. As grid sizing is decreased, the model results converge on a solution. Likewise, insufficient reservoir model area introduces boundary effects in the numerical solution that cause the model results to differ from the analytical solution.

15 GEOTHERMAL ENERGY↗

DOE OSTI2014

DEFLECTION OF A HETEROGENEOUS WIDE-BEAM UNDER UNIFORM PRESSURE LOAD

Oregon State University (OSU) and the Idaho National Laboratory (INL) are currently collaborating on a test program which entails hydro-mechanical testing of a generic plate type fuel element, or generic test plate assembly (GTPA), for the purpose of qualitatively demonstrating mechanical integrity of uranium-molybdenum monolithic plates as compared to that of uranium aluminum dispersion, and aluminum fuel plates onset by hydraulic forces. This test program supports ongoing work conducted for/by the Global Threat Reduction Initiative (GTRI) Fuels Development Program. This study’s focus supports the ongoing collaborative effort by detailing the derivation of an analytic solution for deflection of a heterogeneous plate under a uniform, distributed load in order to predict the deflection of test plates in the GTPA. The resulting analytical solutions for three specific boundary condition sets are then presented against several test cases of a homogeneous plate. In all test cases considered, the results for both homogeneous and heterogeneous plates are numerically identical to one another, demonstrating correct derivation of the heterogeneous solution. Two additional problems are presents herein that provide a representative deflection profile for the plates under consideration within the GTPA. Furthermore, qualitative observations are made about the influence of a more-rigid internal fuel-meat region and its influence on the overall deflection profile of a plate. Present work is being directed to experimentally confirm the analytical solution’s results using select materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1975

Application of Hamilton's law of varying action

The law of varying action enunciated by Hamilton in 1834-1835 permits the direct analytical solution of the problems of mechanics, both stationary and nonstationary, without consideration of force equilibrium and the theory of differential equations associated therewith. It has not been possible to obtain direct analytical solutions to nonstationary systems through the use of energy theory, which has been limited for 140 years to the principle of least action and to Hamilton's principle. It is shown here that Hamilton's law permits the direct analytical solution to nonstationary, initial value systems in the mechanics of solids without any knowledge or use of the theory of differential equations. Solutions are demonstrated for nonconservative, nonstationary particle motion, both linear and nonlinear.

Bailey, C. D.↗

NASA NTRS2001

Electronic Structure and Properties of Deformed Carbon Nanotubes

A theoretical framework based on Huckel tight-binding model has been formulated to analyze the electronic structure of carbon nanotubes under uniform deformation. The model successfully quantifies the dispersion relation, density of states and bandgap change of nanotubes under uniform stretching, compression, torsion and bending. Our analysis shows that the shifting of the Fermi point away from the Brillouin zone vertices is the key reason for these changes. As a result of this shifting, the electronic structure of deformed carbon nanotubes varies dramatically depending on their chirality and deformation mode. Treating the Fermi point as a function of strain and tube chirality, the analytical solution preserves the concise form of undeformed carbon nanotubes. It predicts the shifting, merging and splitting of the Van Hove singularities in the density of states and the zigzag pattern of bandgap change under strains. Four orbital tight-binding simulations of carbon nanotubes under uniform stretching, compression, torsion and bending have been performed to verify the analytical solution. Extension to more complex systems are being performed to relate this analytical solution to the spectroscopic characterization, device performance and proposed quantum structures induced by the deformation. The limitations of this model will also be discussed.

Yang, Liu↗

NASA NTRS1984

Oscillator response to nonstationary excitation

Analytical solutions are presented regarding probability density distributions of various response parameters of a lightly damped oscillator. The oscillator is subjected to a broad-band stochastic excitation which possesses a time-variant power spectrum. The analytical solutions are derived by utilizing appropriate Fokker-Planck equations which govern Markovian approximations of the response parameters considered. The reliability of the approximate analytical solution is tested by using pertinent data generated by a digital Monte Carlo study.

Spanos, P.-T. D.↗

DOE OSTI2021

Verification of numerical models for seismic fluid-structure interaction analysis of internal components in liquid-filled advanced reactors

Earthquake shaking of a liquid-filled advanced reactor induces fluid-structure interaction (FSI) between the reactor vessel (tank), its internal components, and the contained liquid. Seismic design, qualification, and risk assessment of such reactors must consider fluid-structure responses, their geometries and support conditions, and three-directional seismic inputs, all of which require numerical simulations. Numerical models used for analysis of safety-related nuclear equipment must be verified and validated. Here, this paper verifies numerical models of submerged components using analytical solutions. Seismic FSI analysis of the numerical models is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Prior analytical solutions are reworked in this paper, and calculation errors are identified and corrected. The solutions address frequencies of two concentric cylindrical pipes filled with liquid: the inner pipe is a submerged component, and the outer pipe contains the liquid. A unitless frequency coefficient is defined and presented for pipes of different materials, filled with different liquids, and with a range of dimensions for application to advanced reactors. The numerical models are verified here by comparing the lateral frequencies of submerged components with those calculated using the corrected analytical solutions. Recommendations for verification of numerical models of internal components in advanced reactors for seismic FSI analysis are provided. Although the unitless frequency coefficients and verification procedures are developed for application to advanced nuclear reactors, they are broadly applicable to FSI analysis of submerged components in liquid-filled vessels such as storage tanks, boilers, and steam generators.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2018

Numerical similarity solution for a variable coefficient K(m, n) equation

A technique for finding numerical similarity solutions to an initial boundary value problem (IBVP) for generalized K(m, n) equations is described. The equation under consideration is nonlinear and has variable coefficients. The original problem is transformed with the aid of Lie symmetries to an initial value problem (IVP) for a nonlinear third-order ordinary differential equation. The existence and uniqueness of the solution are examined, and the problem is consequently solved with the aid of a finite-difference scheme for various values of the governing parameters. In lieu of an exact symbolic solution, the scheme is validated by comparing the numerical solutions with the approximate analytic solutions obtained with the aid of the method of successive approximations in their region of validity. The accuracy, efficiency, and consistency of the scheme are demonstrated. Numerical solutions to the original initial boundary value problem are constructed for selected parameter values with the aid of the transforms. The qualitative behavior of the solutions as a function of the governing parameters is analyzed, and it is found that the examined IBVPs for generalized K(m, n) equations with variable coefficients that are functions of time, do not admit solitary wave or compacton solutions.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1986

Satellite attitude motion models for capture and retrieval investigations

The primary purpose of this research is to provide mathematical models which may be used in the investigation of various aspects of the remote capture and retrieval of uncontrolled satellites. Emphasis has been placed on analytical models; however, to verify analytical solutions, numerical integration must be used. Also, for satellites of certain types, numerical integration may be the only practical or perhaps the only possible method of solution. First, to provide a basis for analytical and numerical work, uncontrolled satellites were categorized using criteria based on: (1) orbital motions, (2) external angular momenta, (3) internal angular momenta, (4) physical characteristics, and (5) the stability of their equilibrium states. Several analytical solutions for the attitude motions of satellite models were compiled, checked, corrected in some minor respects and their short-term prediction capabilities were investigated. Single-rigid-body, dual-spin and multi-rotor configurations are treated. To verify the analytical models and to see how the true motion of a satellite which is acted upon by environmental torques differs from its corresponding torque-free motion, a numerical simulation code was developed. This code contains a relatively general satellite model and models for gravity-gradient and aerodynamic torques. The spacecraft physical model for the code and the equations of motion are given. The two environmental torque models are described.

Cochran, John E., Jr.↗

DOE OSTI2011

Quantifying non-Markovianity of continuous-variable Gaussian dynamical maps

We introduce a non-Markovianity measure for continuous-variable open quantum systems based on the idea put forward in H.-P. Breuer et al.[Phys. Rev. Lett. 103, 210401 (2009);], that is, by quantifying the flow of information from the environment back to the open system. Instead of the trace distance we use here the fidelity to assess distinguishability of quantum states. We employ our measure to evaluate non-Markovianity of two paradigmatic Gaussian channels: the purely damping channel and the quantum Brownian motion channel with Ohmic environment. We consider different classes of Gaussian states and look for pairs of states maximizing the backflow of information. For coherent states we find simple analytical solutions, whereas for squeezed states we provide both exact numerical and approximate analytical solutions in the weak coupling limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1992

A simplified method for thermal analysis of a cowl leading edge subject to intense local shock-wave-interference heating

Type IV shock wave interference heating on a blunt body causes extremely intense heating over a very localized region of the body. An analytical solution is presented to a heat transfer problem that approximates the shock wave interference heating of an engine cowl leading edge of the National Aero-Space Plane. The problem uses a simplified geometry to represent the leading edge. An analytical solution is developed that provides a means for approximating maximum temperature differences between the outer and inner surface temperatures of the leading edge. The solution is computationally efficient and, as a result, is well suited for conceptual and preliminary design or trade studies. Transient and steady state analyses are conducted, and results obtained from the analytical solution are compared with results of 2-D thermal finite element analyses over a wide range of design parameters. Isotropic materials as well as laminated composite materials are studied. Results of parametric studies are presented to indicate the effects of the thickness of the cowl leading edge and the width of the region heated by the shock wave interference on the thermal response of the leading edge.

Mcgowan, David M.↗

DOE OSTI2020

Thermo-mechanical modeling and experimental validation for multilayered metallic microstructures

This paper presents thermo-mechanical modeling on multilayered metallic microstructures consisting of two bilayers that enable 3D deployment once being Joule heated. Additionally, a representative design of such microstructures, as well as their working principle is introduced. Furthermore, a mathematical model is derived using Euler–Bernoulli beam theory, and analytical solutions describing out-of-plane displacement of the multilayered microstructures when uniformly heated, which is validated by finite element analysis. Parametric analysis on thermal load and dimensional change of beam shows a good agreement between the analytical solutions and results given by finite element analysis. The phases of the near-equiatomic NiTi layer are analyzed by EDS and XRD, which prove to be austenite during the Joule heating process. The mechanical properties of austenite NiTi are incorporated into the analytical solutions that provides a good estimation of 3D deployment of microstructures made of NiTi and aluminum. The experimental results provide an approximately 7 µm out-of-plane deployment by applying a uniform temperature increase of 132 K, and parametric analysis on the dimension of top aluminum beam offers another promising approach to larger 3D deployment. The proposed mathematical model provides an efficient tool for predicting the out-of-plane deployment of such multilayered microstructures.

36 MATERIALS SCIENCE↗

DOE OSTI2016

SIERRA Low Mach Module: Fuego Verification Manual (V. 4.40)

The verification process is intended to demonstrate that a code is solving the governing equations correctly. Verification problems consist of analytical solutions to the mathematical model equations, derived analytical solutions based on an equivalent source term generation procedure (i.e., manufactured solutions), and established numerical bench mark solutions.

42 ENGINEERING↗