Search NASA⌕ Search

SEARCH · Search NASA

Results for “analytic solution”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 19 records

Geodesically complete analytic solutions for a cyclic universe

We present analytic solutions to a class of cosmological models described by a canonical scalar field minimally coupled to gravity and experiencing self interactions through a hyperbolic potential. Using models and methods inspired by 2T-physics, we show how analytic solutions can be obtained in flat/open/closed Friedmann-Robertson-Walker universes. Among the analytic solutions, there are many interesting geodesically complete cyclic solutions in which the universe bounces at either zero or finite sizes. When geodesic completeness is imposed, it restricts models and their parameters to a certain parameter subspace, including some quantization conditions on initial conditions in the case of zero-size bounces, but no conditions on initial conditions for the case of finite-size bounces. We will explain the theoretical origin of our model from the point of view of 2T-gravity as well as from the point of view of the colliding branes scenario in the context of M-theory. We will indicate how to associate solutions of the quantum Wheeler-deWitt equation with our classical analytic solutions, mention some physical aspects of the cyclic solutions, and outline future directions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Unsteady fluid flow in a slightly curved pipe: A comparative study of a matched asymptotic expansions solution with a single analytical solution

The present work is motivated by the fact that blood flow in the aorta and the main arteries is governed by large finite values of the Womersley number α and for such values of α there is not any analytical solution in the literature. The existing numerical solutions, although accurate, give limited information about the factors that affect the flow, whereas an analytical approach has an advantage in that it can provide physical insight to the flow mechanism. Having this in mind, we seek analytical solution to the equations of the fluid flow driven by a sinusoidal pressure gradient in a slightly curved pipe of circular cross section when the Womersley number varies from small finite to infinite values. Initially the equations of motion are expanded in terms of the curvature ratio δ and the resulting linearized equations are solved analytically in two ways. In the first, we match the solution for the main core to that for the Stokes boundary layer. This solution is valid for very large values of α. In the second, we derive a straightforward single solution valid to the entire flow region and for 8 ≤ α < ∞, a range which includes the values of α that refer to the physiological flows. Each solution contains expressions for the axial velocity, the stream function, and the wall stresses and is compared to the analogous forms presented in other studies. The two solutions give identical results to each other regarding the axial flow but differ in the secondary flow and the circumferential wall stress, due to the approximations employed in the matched asymptotic expansion process. The results on the stream function from the second solution are in agreement with analogous results from other numerical solutions. The second solution predicts that the atherosclerotic plaques may develop in any location around the cross section of the aortic wall unlike to the prescribed locations predicted by the first solution. In addition, it gives circumferential wall stresses augmented by approximately 100% with respect to the matched asymptotic expansions, a factor that may contribute jointly with other pathological factors to the faster aging of the arterial system and the possible malfunction of the aorta.

42 ENGINEERING↗

Comparing the MELCOR Aerosol Deposition Model with Exact Analytical Solutions

Exact analytical solutions are presented for the evolution of the aerosol particle mass density function in a control volume for particle deposition due to gravitational settling, thermophoresis, and diffusion. The solutions are for arbitrary initial mass density functions and are applied for an initial lognormal density function. Integration of these solutions provides the suspended mass in the control volume as a function of time. These solutions serve as an exact benchmark to assess the accuracy of numerical methods. For the numerical algorithm used in MELCOR, excellent agreement is obtained for gravitational settling, diffusive deposition, and thermophoretic deposition for the suspended aerosol mass. In all cases, the default number of discrete particle size bins of 10 is shown to converge, with hardly any advantage to using 20 size bins.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Analytical solutions of the planar cyclic voltammetry process for two soluble species with equal diffusivities and fast electron transfer using the method of eigenfunction expansions

Cyclic voltammetry is a powerful tool that is used for characterizing electrochemical processes. Models of cyclic voltammetry take into account the mass transport of species and the kinetics at the electrode surface. Analytical solutions of these models are not well-known due to the complexity of the boundary conditions. In this study we present closed form analytical solutions of the planar voltammetry model for two soluble species with fast electron transfer and equal diffusivities using the eigenfunction expansion method. Our solution methodology does not incorporate Laplace transforms and yields good agreement with the numerical solution. This solution method can be extended to cases that are more general and may be useful for benchmarking purposes.

36 MATERIALS SCIENCE↗

New analytical solutions for dust acoustic solitary and periodic waves in an unmagnetized dusty plasma with kappa distributed electrons and ions

Dust acoustic solitary waves and periodic waves in an unmagnetized dusty plasma with kappa distributed electrons and ions are investigated through non-perturbative approach. Basic equations are reduced to a system of ordinary differential equations involving electrostatic potential. After that by applying the bifurcation theory of planar dynamical systems to this system, we have studied the existence of solitary wave solutions and periodic wave solutions. New analytical solutions for the above waves are derived depending on the parametric space. Regarding the solitary and periodic wave solutions, the combined effects of temperature ratio (σ) of ions and electrons, spectral index (κ) and density ratio (p) are studied on characteristics of dust acoustic (DA) solitary waves and periodic waves. The spectral index (κ), density ratio (p) of ions and electrons and temperature ratio (σ) significantly influence the characteristics of dust acoustic solitary and periodic structures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Clausius–Clapeyron Scaling of CAPE from Analytical Solutions to RCE

By deriving analytical solutions to radiative-convective equilibrium (RCE), it is shown mathematically that convective available potential energy (CAPE) exhibits Clausius-Clapeyron (CC) scaling over a wide range of surface temperatures up to 310 K. Above 310 K, CAPE deviates from CC scaling and even decreases with warming at very high surface temperatures. At the surface temperature of the current tropics, the analytical solutions predict that CAPE increases at a rate of about 6%-7% per kelvin of surface warming. The analytical solutions also provide insight on how the tropopause height and stratospheric humidity change with warming. Changes in the tropopause height exhibit CC scaling, with the tropopause rising by about 400 m per kelvin of surface warming at current tropical temperatures and by about 1-2 km K -1 at surface temperatures in the range of 320-340 K. The specific humidity of the stratosphere exhibits super-CC scaling at temperatures moderately warmer than the current tropics. With a surface temperature of the current tropics, the stratospheric specific humidity increases by about 6% per kelvin of surface warming, but the rate of increase is as high as 30% K -1 at warmer surface temperatures.

58 GEOSCIENCES↗

JOVIAN STRATOSPHERE AS A CHEMICAL TRANSPORT SYSTEM: BENCHMARK ANALYTICAL SOLUTIONS

We systematically investigated the solvable analytical benchmark cases in both one- and two-dimensional (1D and 2D) chemical-advective-diffusive systems. We use the stratosphere of Jupiter as an example but the results can be applied to other planetary atmospheres and exoplanetary atmospheres. In the 1D system, we show that CH{sub 4} and C{sub 2}H{sub 6} are mainly in diffusive equilibrium, and the C{sub 2}H{sub 2} profile can be approximated by modified Bessel functions. In the 2D system in the meridional plane, analytical solutions for two typical circulation patterns are derived. Simple tracer transport modeling demonstrates that the distribution of a short-lived species (such as C{sub 2}H{sub 2}) is dominated by the local chemical sources and sinks, while that of a long-lived species (such as C{sub 2}H{sub 6}) is significantly influenced by the circulation pattern. We find that an equator-to-pole circulation could qualitatively explain the Cassini observations, but a pure diffusive transport process could not. For slowly rotating planets like the close-in extrasolar planets, the interaction between the advection by the zonal wind and chemistry might cause a phase lag between the final tracer distribution and the original source distribution. The numerical simulation results from the 2D Caltech/JPL chemistry-transport model agree well with the analytical solutions for various cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analytical solution for two-phase flow in a wellbore using the drift-flux model

This paper presents analytical solutions for steady-state, compressible two-phase flow through a wellbore under isothermal conditions using the drift flux conceptual model. Although only applicable to highly idealized systems, the analytical solutions are useful for verifying numerical simulation capabilities that can handle much more complicated systems, and can be used in their own right for gaining insight about two-phase flow processes in wells. The analytical solutions are obtained by solving the mixture momentum equation of steady-state, two-phase flow with an assumption that the two phases are immiscible. These analytical solutions describe the steady-state behavior of two-phase flow in the wellbore, including profiles of phase saturation, phase velocities, and pressure gradients, as affected by the total mass flow rate, phase mass fraction, and drift velocity (i.e., the slip between two phases). Close matching between the analytical solutions and numerical solutions for a hypothetical CO{sub 2} leakage problem as well as to field data from a CO{sub 2} production well indicates that the analytical solution is capable of capturing the major features of steady-state two-phase flow through an open wellbore, and that the related assumptions and simplifications are justified for many actual systems. In addition, we demonstrate the utility of the analytical solution to evaluate how the bottomhole pressure in a well in which CO{sub 2} is leaking upward responds to the mass flow rate of CO{sub 2}-water mixture.

42 ENGINEERING↗

Analytic Solution for Time-Dependent Neutron Transport with Delayed Neutrons via Singular Eigenfunction Expansion

An important way to verify the proper functioning of a computer code, and its underlying numerical schemes, is by comparing numerical results generated with the code to an analytic solution. To this end, we present an analytic solution for time-dependent neutron transport with delayed neutrons, which can be used to verify computer codes for modeling reactor kinetics phenomena. To develop our analytic solution, we use the singular eigenfunction expansion method. Ganapol et al. have also developed a solution to this problem, but by using the multiple collision technique. Our approach has certain advantages with respect to this previous work, in particular that it yields a closed form (i.e., in terms of an integral rather than an infinite series) explicit expression for the angular flux (i.e., instead of just the scalar flux or the angular flux at a single spatial point). In addition, Su and Olson have presented an analytic solution for a very similar problem: non-equilibrium radiative transfer. Their approach involves a Laplace transform, which they invert by showing there are no singularities in the right complex plane of the Laplace transform parameter and then simply choosing the Bromwich contour as the imaginary axis. To extend this previous work to our case, we would further need to exclude supercritical problems in order to use the same contour. The resulting analytic solution would then be of limited interest. In contrast, our analytic solution is completely general with respect to sub, delayed super, or prompt supercriticality. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Time evolution of the QED vacuum in a uniform electric field: Complete analytic solution by spinorial decomposition

Exact analytical solutions are presented for the time evolution of the density of pairs produced in the QED vacuum by a time-independent, uniform electric field. The mathematical tool used here to describe the pair production is the Dirac-Heisenberg-Wigner function introduced before [Phys. Rev. D 44, 1825 (1991)]. The initial value problem for this function is solved by decomposing the solution into a product of spinors. The equations for spinors are much simpler and are solved analytically. These calculations are nonperturbative since pair production is due to quantum-mechanical tunneling and the explicit solutions clearly exhibit their nonanalytic behavior.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Analytical solutions of the Arrhenius-Semenov problem for constant volume burn

Analytical solutions to the Semenov thermal ignition problem for constant volume burn governed by Arrhenius reaction kinetics are derived. Specifically, an approximate analytical solution technique for the Arrhenius-Semenov differential equation is derived for reaction orders n ϵ R> 0 and exact solutions are also constructed for reaction orders n ϵ N : n ≤ 3. The approximation technique relies on expansion of the respective nondominant terms in the differential equation at the lower and upper bounds of the reaction progress variable in order to create a pair of integrable series. The two integrated series are then connected to create a single continuous analytical solution. Excellent agreement is observed between the analytical approximation and solutions obtained numerically. The presented approximation constitutes a simple and robust strategy for solving the Arrhenius-Semenov problem analytically.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

New analytical solutions to the two-phase water faucet problem

Here, the one-dimensional water faucet problem is one of the classical benchmark problems originally proposed by Ransom to study the two-fluid two-phase flow model. With certain simplifications, such as massless gas phase and no wall and interfacial frictions, analytical solutions had been previously obtained for the transient liquid velocity and void fraction distribution. The water faucet problem and its analytical solutions have been widely used for the purposes of code assessment, benchmark and numerical verifications. In our previous study, the Ransom’s solutions were used for the mesh convergence study of a high-resolution spatial discretization scheme. It was found that, at the steady state, an anticipated second-order spatial accuracy could not be achieved, when compared to the existing Ransom’s analytical solutions. A further investigation showed that the existing analytical solutions do not actually satisfy the commonly used two-fluid single-pressure two-phase flow equations. In this work, we present a new set of analytical solutions of the water faucet problem at the steady state, considering the gas phase density’s effect on pressure distribution. This new set of analytical solutions are used for mesh convergence studies, from which anticipated second-order of accuracy is achieved for the 2nd order spatial discretization scheme. In addition, extended Ransom’s transient solutions for the gas phase velocity and pressure are derived, with the assumption of decoupled liquid and gas pressures. Numerical verifications on the extended Ransom’s solutions are also presented.

97 MATHEMATICS AND COMPUTING↗

Analytic solutions of an unclassified artifact

This report provides the technical detail for analytic solutions for the inner and outer profiles of the unclassified CMM Test Artifact (LANL Part Number 157Y-700373, 5/03/2001) in terms of radius and polar angle. Furthermore, analytic solutions are derived for the legacy Sheffield measurement hardware, also in terms of radius and polar angle, using part coordinates, i.e., relative to the analytic profile solutions obtained. The purpose of this work is to determine the exact solution for the “cosine correction” term inherent to measurement with the Sheffield hardware. The cosine correction is required in order to interpret the actual measurements taken by the hardware in terms of an actual part definition, or “knot-point spline definition,” that typically accompanies a component drawing. Specifically, there are two portions of the problem: first an analytic solution must be obtained for any point on the part, e.g., given the radii and the straight lines that define the part, it is required to find an exact solution for the inner and outer profile for any arbitrary polar angle. Next, the problem of the inspection of this part must be solved, i.e., given an arbitrary sphere (representing the inspection hardware) that comes in contact with the part (inner and outer profiles) at any arbitrary polar angle, it is required to determine the exact location of that intersection. This is trivial for the case of concentric circles. In the present case, however, the spherical portion of the profiles is offset from the defined center of the part, making the analysis nontrivial. Here, a simultaneous solution of the part profiles and the sphere was obtained.

42 ENGINEERING↗