Results for “analytic solution” — Page 25

DOE OSTI2015

Steady state in a gas of inelastic rough spheres heated by a uniform stochastic force

We study here the steady state attained in a granular gas of inelastic rough spheres that is subject to a spatially uniform random volume force. The stochastic force has the form of the so-called white noise and acts by adding impulse to the particle translational velocities. We work out an analytical solution of the corresponding velocity distribution function from a Sonine polynomial expansion that displays energy non-equipartition between the translational and rotational modes, translational and rotational kurtoses, and translational-rotational velocity correlations. By comparison with a numerical solution of the Boltzmann kinetic equation (by means of the direct simulation Monte Carlo method), we show that our analytical solution provides a good description that is quantitatively very accurate in certain ranges of inelasticity and roughness. We also find three important features that make the forced granular gas steady state very different from the homogeneous cooling state (attained by an unforced granular gas). First, the marginal velocity distributions are always close to a Maxwellian. Second, there is a continuous transition to the purely smooth limit (where the effects of particle rotations are ignored). And third, the angular translational-rotational velocity correlations show a preference for a quasiperpendicular mutual orientation (which is called “lifted-tennis-ball” behavior)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

Estimated Fluid Force and Damping Characteristics of a Thin Film Damper Comparison Between Closed-Form Solutions and Numerical Analysis

Analytical solutions of thin film dampers are useful for determining critical speeds and stability of rotor systems. Most thin film dampers in use are of short axial length, and closed-form solutions to the Reynolds equations exist for estimating pressure, forces, and damping for these types of dampers. This article compares the fluid film forces and damping estimated by the short film bearing model form of the Reynolds equations to the calculated forces and damping of a transient computational fluid dynamic simulation. For this comparison, the fluid was assumed to be incompressible, laminar, and isoviscous. The fluid film forces and damping are calculated from integrating the pressure distribution over the surface of the damper due to small amplitude motions about a steady state static off-center circular orbit. In this case, no cavitation is assumed, and the journal has no angular velocity, so direct stiffness cannot be calculated from the closed-form solution. Radial clearance, journal length, and journal eccentricity have a significant effect on fluid force and damping within a thin film damper. Fluid density does not affect fluid force or damping substantially, while fluid viscosity does. Both the closed-form solutions and computational fluid dynamics simulation compare well with each other and reflect these trends.

Cook, Jason↗

DOE OSTI2015

PARTICLE ACCELERATION AT LOW CORONAL COMPRESSION REGIONS AND SHOCKS

We present a study on particle acceleration in the low corona associated with the expansion and acceleration of coronal mass ejections (CMEs). Because CME expansion regions low in the corona are effective accelerators over a finite spatial region, we show that there is a rigidity regime where particles effectively diffuse away and escape from the acceleration sites using analytic solutions to the Parker transport equation. This leads to the formation of broken power-law distributions. Based on our analytic solutions, we find a natural ordering of the break energy and second power-law slope (above the break energy) as a function of the scattering characteristics. These relations provide testable predictions for the particle acceleration from low in the corona. Our initial analysis of solar energetic particle observations suggests a range of shock compression ratios and rigidity dependencies that give rise to the solar energetic particle (SEP) events studied. The wide range of characteristics inferred suggests competing mechanisms at work in SEP acceleration. Thus, CME expansion and acceleration in the low corona may naturally give rise to rapid particle acceleration and broken power-law distributions in large SEP events.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Close-contact melting on an isothermal surface with the inclusion of non-Newtonian effects

The present study deals with a theoretical investigation of a close-contact melting (CCM) process involving a vertical cylinder on a horizontal isothermal surface, where the liquid phase is a non-Newtonian viscoplastic fluid that behaves according to the Bingham model. Accordingly, a new approach is formulated based on the thin layer approximation and different quasi-steady process assumptions. By analytical derivation, an algebraic equation that relates the molten layer thickness and the solid bulk height is developed. The problem is then solved numerically, coupled with another equation for the melting rate. The new model shows that as the yield stress increases the melting rate decreases and the molten layer thickness increases. It is found that under certain conditions, the model can be reduced to a form that allows an analytical solution. The approximate model predicts an exponential dependence of both the melt fraction and the molten layer thickness. Comparison between the numerical and analytical solutions shows that the analytical approximation provides an excellent estimation for sufficiently large values of the yield stress. Dimensional analysis, which is supported by the analytical model results, reveals the dimensionless groups that govern the problem. For the general case, the melt fraction is a function of two dimensionless groups. For the analytical approximation, it is shown that the melt fraction is governed by a single dimensionless group and that the molten layer thickness is governed by two dimensionless groups.

Mechanics↗

DOE OSTI2019

Analysis of the Spectral Stability of the Generalized Runge–Kutta Methods Applied to Initial-Boundary-Value Problems for Equations of the Parabolic Type. II. Implicit Methods

We consider specific realizations of different implicit generalized Runge–Kutta methods as applied to the numerical integration with respect to time of initial-boundary-value problems for the second-order parabolic equations and investigate their spectral stability. It is shown that all implicit generalized Runge–Kutta methods are unconditionally spectrally stable but some of them have the conditional property of monotonicity of the numerical solution with respect to time. The functions of spectral stability of the implicit generalized Runge–Kutta methods are rational. We compare the analytic solution of the nonstationary one-dimensional problem of heat conduction with the numerical solutions of this problem obtained by different implicit generalized Runge–Kutta methods. It is shown that, in this case, the application of the one-stage Radau methods with subsequent discretization of the problem with respect to the space variable leads to the classical forward finite difference scheme (Laasonen scheme), whereas the use of the one-stage Gauss–Legendre method leads to a six-point symmetric scheme (Crank–Nicolson scheme). It is shown that diagonally implicit generalized Nørsett and Burrage methods are realized in almost the same way as the one-stage Radau and Gauss–Legendre methods but their accuracy in the time step is 10–1000 times higher. On the basis of comparison of the numerical and analytic solutions, we conclude that, in order to get practically suitable numerical solutions without any restrictions on the time step, it is reasonable to use one- and three-stage generalized Radau methods or two- and four-stage Lobatto IIIC methods. All other explicit and implicit generalized Runge–Kutta methods require certain restrictions imposed on the time step.

42 ENGINEERING↗

DOE OSTI2016

Regarding on the prototype solutions for the nonlinear fractional-order biological population model

In this study, we have submitted to literature a method newly extended which is called as Improved Bernoulli sub-equation function method based on the Bernoulli Sub-ODE method. The proposed analytical scheme has been expressed with steps. We have obtained some new analytical solutions to the nonlinear fractional-order biological population model by using this technique. Two and three dimensional surfaces of analytical solutions have been drawn by wolfram Mathematica 9. Finally, a conclusion has been submitted by mentioning important acquisitions founded in this study.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2013

Alfvén wave collisions, the fundamental building block of plasma turbulence. I. Asymptotic solution

The nonlinear interaction between counterpropagating Alfvén waves is the physical mechanism underlying the cascade of energy to small scales in astrophysical plasma turbulence. Beginning with the equations for incompressible MHD, an asymptotic analytical solution for the nonlinear evolution of these Alfvén wave collisions is derived in the weakly nonlinear limit. The resulting qualitative picture of nonlinear energy transfer due to this mechanism involves two steps: first, the primary counterpropagating Alfvén waves interact to generate an inherently nonlinear, purely magnetic secondary fluctuation with no parallel variation; second, the two primary waves each interact with this secondary fluctuation to transfer energy secularly to two tertiary Alfvén waves. These tertiary modes are linear Alfvén waves with the same parallel wavenumber as the primary waves, indicating the lack of a parallel cascade. The amplitude of these tertiary modes increases linearly with time due to the coherent nature of the resonant four-wave interaction responsible for the nonlinear energy transfer. The implications of this analytical solution for turbulence in astrophysical plasmas are discussed. The solution presented here provides valuable intuition about the nonlinear interactions underlying magnetized plasma turbulence, in support of an experimental program to verify in the laboratory the nature of this fundamental building block of astrophysical plasma turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2012

A Microsoft Excel Program for Two-Well Tracer Test Data Analysis

Two-well tracer tests are often conducted to investigate subsurface solute transport in the field. Analyzing breakthrough curves in the extraction and monitoring wells using numerical methods is nontrivial due to highly nonuniform flow conditions. We extended, and implemented analytical solutions for the convection-dispersion equation for an injection-extraction well-duplet in a homogeneous confined aquifer under steady state conditions. Functions were provided to calculate the concentrations in the extraction and monitoring wells at any location due to a step or pulse injection. Type curves for a step injection were compared with those calculated by numerically integrating the solution for a pulse injection. The results from the two approaches are similar when the dispersivity is small. As the dispersivity increases, the latter was found to be more accurate but requires more computing time. The code was verified by comparing the results with published type curves and applied to analyze data from the literature. It can be used as a first approximation for two-well tracer test data analysis, and to check accuracy of numerical solutions. The code and example files are publically-available.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1993

An extended laser flash technique for thermal diffusivity measurement of high-temperature materials

Knowledge of thermal diffusivity data for high-temperature materials (solids and liquids) is very important in analyzing a number of processes, among them solidification, crystal growth, and welding. However, reliable thermal diffusivity versus temperature data, particularly those for high-temperature liquids, are still far from complete. The main measurement difficulties are due to the presence of convection and the requirement for a container. Fortunately, the availability of levitation techniques has made it possible to solve the containment problem. Based on the feasibility of the levitation technology, a new laser flash technique which is applicable to both levitated liquid and solid samples is being developed. At this point, the analysis for solid samples is near completion and highlights of the technique are presented here. The levitated solid sample which is assumed to be a sphere is subjected to a very short burst of high power radiant energy. The temperature of the irradiated surface area is elevated and a transient heat transfer process takes place within the sample. This containerless process is a two-dimensional unsteady heat conduction problem. Due to the nonlinearity of the radiative plus convective boundary condition, an analytic solution cannot be obtained. Two options are available at this point. Firstly, the radiation boundary condition can be linearized, which then accommodates a closed-form analytic solution. Comparison of the analytic curves for the temperature rise at different points to the experimentally-measured values will then provide the thermal diffusivity values. Secondly, one may set up an inverse conduction problem whereby experimentally obtained surface temperature history is used as the boundary conditions. The thermal diffusivity can then be elevated by minimizing the difference between the real heat flux boundary condition (radiation plus convection) and the measurements. Status of an experimental study directed at measuring the thermal diffusivity of high-temperature solid samples of pure Nickel and Inconel 718 superalloys are presented. Preliminary measurements showing surface temperature histories are discussed.

Shen, F.↗

DOE OSTI2017

Sierra/SolidMechanics 4.46 Verification Tests Manual

Presented in this document is a small portion of the tests that exist in the Sierra/SolidMechanics (Sierra/SM) verification test suite. Most of these tests are run nightly with the Sierra/SM code suite, and the results of the test are checked versus the correct analytical result. For each of the tests presented in this document, the test setup, a description of the analytic solution, and comparison of the Sierra/SM code results to the analytic solution is provided. Mesh convergence is also checked on a nightly basis for several of these tests. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems. Additional example problems are provided in the Sierra/SM Example Problems Manual. Note, many other verification tests exist in the Sierra/SM test suite, but have not yet been included in this manual.

42 ENGINEERING↗

DOE OSTI2018

Sierra Structural Dynamics Verification Test Manual (4.48)

This document presents tests from the Sierra Structural Mechanics verification test suite. Each of these tests is run nightly with the Sierra/SD code suite and the results of the test checked versus the correct analytic result. For each of the tests presented in this document the test setup, derivation of the analytic solution, and comparison of the Sierra/SD code results to the analytic solution is provided. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems.

42 ENGINEERING↗

DOE OSTI2019

Sierra/SolidMechanics 4.54 Verification Tests Manual

Presented in this document is a small portion of the tests that exist in the Sierra/SolidMechanics (Sierra/SM) verification test suite. Most of these tests are run nightly with the Sierra/SM code suite, and the results of the test are checked versus the correct analytical result. For each of the tests presented in this document, the test setup, a description of the analytic solution, and comparison of the Sierra/SM code results to the analytic solution is provided. Mesh convergence is also checked on a nightly basis for several of these tests. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems. Additional example problems are provided in the Sierra/SM Example Problems Manual. Note, many other verification tests exist in the Sierra/SM test suite, but have not yet been included in this manual.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Sierra/SD - Verification Test Manual - 4.56

This document presents tests from the Sierra Structural Mechanics verification test suite. Each of these tests is run nightly with the Sierra/SD code suite and the results of the test checked versus the correct analytic result. For each of the tests presented in this document the test setup, derivation of the analytic solution, and comparison of the Sierra/SD code results to the analytic solution is provided. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Sierra/SolidMechanics 4.58 Verification Tests Manual

Presented in this document is a small portion of the tests that exist in the Sierra/SolidMechanics (Sierra/SM) verification test suite. Most of these tests are run nightly with the Sierra/SM code suite, and the results of the test are checked versus the correct analytical result. For each of the tests presented in this document, the test setup, a description of the analytic solution, and comparison of the Sierra/SM code results to the analytic solution is provided. Mesh convergence is also checked on a nightly basis for several of these tests. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems. Additional example problems are provided in the Sierra/SM Example Problems Manual. Note, many other verification tests exist in the Sierra/SM test suite, but have not yet been included in this manual.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2019

Sierra/SolidMechanics 4.52 Verification Tests Manual

Presented in this document is a small portion of the tests that exist in the Sierra/SolidMechanics (Sierra/SM) verification test suite. Most of these tests are run nightly with the Sierra/SM code suite, and the results of the test are checked versus the correct analytical result. For each of the tests presented in this document, the test setup, a description of the analytic solution, and comparison of the Sierra/SM code results to the analytic solution is provided. Mesh convergence is also checked on a nightly basis for several of these tests. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems. Additional example problems are provided in the Sierra/SM Example Problems Manual. Note, many other verification tests exist in the Sierra/SM test suite, but have not yet been included in this manual.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Sierra/SD – Verification Test Manual – 5.0

This document presents tests from the Sierra Structural Mechanics verification test suite. Each of these tests is run nightly with the Sierra/SD code suite and the results of the test checked versus the correct analytic result. For each of the tests presented in this document the test setup, derivation of the analytic solution, and comparison of the Sierra/SD code results to the analytic solution is provided. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Sierra/SD–Verification Test Manual - 4.58

This document presents tests from the Sierra Structural Mechanics verification test suite. Each of these tests is run nightly with the Sierra/SD code suite and the results of the test checked versus the correct analytic result. For each of the tests presented in this document the test setup, derivation of the analytic solution, and comparison of the Sierra/SD code results to the analytic solution is provided. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Sierra/SolidMechanics 5.2 Verification Tests Manual

Presented in this document is a small portion of the tests that exist in the Sierra/SolidMechanics (Sierra / SM) verification test suite. Most of these tests are run nightly with the Sierra/SM code suite, and the results of the test are checked versus the correct analytical result. For each of the tests presented in this document, the test setup, a description of the analytic solution, and comparison of the Sierra/SM code results to the analytic solution is provided. Mesh convergence is also checked on a nightly basis for several of these tests. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems. Additional example problems are provided in the Sierra/SM Example Problems Manual. Note, many other verification tests exist in the Sierra/SM test suite, but have not yet been included in this manual.

97 MATHEMATICS AND COMPUTING↗