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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Approximate solutions for diffusive fracture-matrix transfer: Application to storage of dissolved CO 2 in fractured rocks

Analytical solutions with infinite exponential series are available to calculate the rate of diffusive transfer between low-permeability blocks and high-permeability zones in the subsurface. Truncation of these series is often employed by neglecting the early-time regime. Here in this paper, we present unified-form approximate solutions in which the early-time and the late-time solutions are continuous at a switchover time. The early-time solutions are based on three-term polynomial functions in terms of square root of dimensionless time, with the first coefficient dependent only on the dimensionless area-to-volume ratio. The last two coefficients are either determined analytically for isotropic blocks (e.g., spheres and slabs) or obtained by fitting the exact solutions, and they solely depend on the aspect ratios for rectangular columns and parallelepipeds. For the late-time solutions, only the leading exponential term is needed for isotropic blocks, while a few additional exponential terms are needed for highly anisotropic rectangular blocks. The optimal switchover time is between 0.157 and 0.229, with highest relative approximation error less than 0.2%. The solutions are used to demonstrate the storage of dissolved CO 2 in fractured reservoirs with low-permeability matrix blocks of single and multiple shapes and sizes. These approximate solutions are building blocks for development of analytical and numerical tools for hydraulic, solute, and thermal diffusion processes in low-permeability matrix blocks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

Nonlinear theory of the compressible gas flows over a nonuniform boundary in the gravitational field in the shallow-water approximation

A set of equations is derived for the motion of a compressible ideal gas over a nonuniform boundary in the gravitational field in the shallow-water approximation. Classical simple waves are shown not to be the solutions to this set of equations. Generalized simple waves are found to exist only in the case of a linear underlying-surface profile. All continuous and discontinuous solutions are obtained in an explicit form for the case of the boundary in the form of an inclined plane, and an analytical solution is found for the problem of the decay of an arbitrary discontinuity. This solution consists of four wave configurations. Necessary and sufficient conditions are determined for the existence of each configuration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Exact transient photon correlation with arbitrary laser pulses

We present a full quantum theory to study the transient evolution of photon pairs. We introduce a method which gives exact time-dependent solutions of the coupled quantum Langevin equations for a multilevel quantum particle driven by arbitrary time-dependent laser fields. The analytical solutions are used to develop a numerical code for computing exact time evolution of the two-photon correlation function. We analyze the effects of laser pulses sequence, pulse duration, chirping, and initial internal quantum states on the nonclassicality of the photon correlation through the violation of the Cauchy-Schwarz inequality. The results provide a promising possibility of controlling the generation of highly correlated photon pairs using tailored short laser pulses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Annular spherically focused ring transducers for improved single-beam acoustical tweezers

The use of ultrasonic transducers with a central hollow is suggested for improved single-beam acoustical tweezers applications. Within the framework of the Fresnel-Kirchhoff parabolic approximation, a closed-form partial-wave series expansion (PWSE) for the incident velocity potential (or pressure) field is derived for an annular spherically focused ring (asfr) with uniform vibration across its surface in spherical coordinates. The Rayleigh-Sommerfeld diffraction integral and the addition theorems for the Legendre and spherical wave functions are used to obtain the PWSE assuming a weakly focused beam (with a focusing angle α ≤ 20°). The PWSE allows evaluating the incident field from the finite asfr in 3D. Moreover, the obtained solution allows computing efficiently the acoustic scattering and radiation force on a sphere centered on the beam's axis of wave propagation. The analytical solution is valid for wavelengths largely exceeding the radius of the asfr and when the viscosity of the surrounding fluid can be neglected. Numerical predictions for the beam-forming, scattering, and axial time-averaged radiation force are performed with particular emphasis on the asfr thickness, the axial distance separating the sphere from the center of the transducer, the (non-dimensional) size of the transducer, as well as the sphere's elastic properties without restriction to the long- (i.e., Rayleigh) or the short-wavelength (i.e., ray acoustics) regimes. Potential applications of the present solution are in beam-forming design, particle tweezing, and manipulation due to negative forces using ultrasonic asfr transducers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2010

Simulation of guided modes (eigenmodes) and synthesis of a thin-film generalised waveguide Luneburg lens in the zero-order vector approximation

Propagation of a guided mode (eigenmode) through an integrated optical generalised waveguide Luneburg lens is numerically simulated for the first time in terms of the previously obtained analytical solution of the vector electrodynamic problem in a smoothly irregular four-layer integrated optical 3D waveguide. The dispersion relation for a four-layer continuously irregular integrated optical 3D waveguide is calculated within the approximations of the asymptotic method of comparison waveguides and the method of adiabatic modes, in particular, taking into account the shift of the propagation constants of quasi-TE and quasi-TM modes. A generalised waveguide Luneburg lens with a full aperture is synthesised in the zero-order approximation. The results of numerical simulation demonstrate, on the one hand, a very good coincidence of the solution to the stated problem obtained in the approximation of the method of comparison waveguides with the previous results, and, on the other hand, advantages of our method: more rigorous solution of the problem, more complete consideration of its physical peculiarities, and higher accuracy of calculations. Another undoubted advantage of the analytical method proposed here is that it can be used to analyse similar structures fabricated of dielectrics, magnetics, and metamaterials, including nonlinear materials, in a wide range of electromagnetic wavelengths. (fibre optics)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Vortex patterns in quasi-two-dimensional flows of a viscous rotating fluid

A modified von Karman problem that describes steady vortex flow in a rotating thin viscous fluid layer is solved. An analysis of the effect of bottom friction on the behavior of cyclonic and anticyclonic vortices at arbitrary values of the Rossby number is presented. Several anticyclonic flow patterns are examined. An approximate analytical solution obtained for steady flows is compared with numerical computations of a time-dependent problem. Experimental results on cyclonic and anticyclonic vortices in multiple-vortex quasi-turbulent flow are presented, and their interpretation based on the solution of the model problem is given.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2011

Phase-field modeling of two-dimensional solute precipitation/dissolution: Solid fingers and diffusion-limited precipitation

Two-dimensional dendritic growth due to solute precipitation was simulated using a phase-field model reported earlier [Z. Xu and P. Meakin, J. Chem. Phys. 129, 014705 (2008)]. It was shown that diffusion-limited precipitation due to the chemical reaction at the solid–liquid interface has similarities with diffusion-limited aggregation (DLA). The diffusion-limited precipitation is attained by setting the chemical reaction rate much larger compared to the solute diffusion to eliminate the effect of the interface growth kinetics. The phase-field simulation results were in reasonable agreement with the analytical solutions. The fractal solid fingers can be formed in the diffusion-limited precipitation and have a fractal dimension measured df = 1.68, close to 1.64, the fractal dimensionality of large square lattice DLA clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2011

Phase-field modeling of two-dimensional solute precipitation/dissolution: Solid fingers and diffusion-limited precipitation

Two-dimensional dendritic growth due to solute precipitation was simulated using a phase-field model reported earlier [Z. Xu and P. Meakin, J. Chem. Phys. 129, 014705 (2008)]. It was shown that diffusion-limited precipitation due to the chemical reaction at the solid-liquid interface posses similarities with diffusion-limited aggregation (DLA). The diffusion-limited precipitation is attained by setting the chemical reaction rate much larger compared to the solute diffusion to eliminate the effect of the interface growth kinetics. The phase-field simulation results were in reasonable agreement with the analytical solutions. The fractal solid fingers can be formed in the diffusion-limited precipitation and have a fractal dimension measured , close to 1.64, the fractal dimensionality of large square lattice diffusion-limited aggregation (DLA) clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2012

Exact self-similar solutions for the magnetized Noh Z pinch problem

A self-similar solution is derived for a radially imploding cylindrical plasma with an embedded, azimuthal magnetic field. The plasma stagnates through a strong, outward propagating shock wave of constant velocity. This analysis is an extension of the classic Noh gasdynamics problem to its ideal magnetohydrodynamics (MHD) counterpart. The present exact solution is especially suitable as a test for MHD codes designed to simulate linear Z pinches. To demonstrate the application of the new solution to code verification, simulation results from the cylindrical R-Z version of Mach2 and the 3D Cartesian code Athena are compared against the analytic solution. Alternative routines from the default ones in Athena lead to significant improvement of the results, thereby demonstrating the utility of the self-similar solution for verification.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

SAFR/V1 (EVKLID/V2 Integral Code Module) Aided Simulation of Melt Movement Along the Surface of a Fuel Element in a Fast Reactor During a Serious Accident

The basic approaches used in the SAFR/V1 module of the integral code EVKLID/V2 to simulate the movement of the melt formed upon melting of fuel elements are presented. The system of mass, energy, and momentum conservation equations used to simulate the movement of melt is presented. Special attention is devoted to methods of numerical approximation of the equations as well as to the solution of problems involving smearing of the solution at the melt boundary. The realized methods of stimulating the motion of melt have been verified on the basis of tests with known analytical solutions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1993

Transport methods and interactions for space radiations

This report presents a brief history leading to the involvement of the Langley Research Center of the National Aeronautics and Space Administration (NASA) in space-radiation physics and protection. Indeed, a relatively complete summary of technical capability as of the summer of 1990 is given. The Boltzmann equations for coupled ionic and neutronic fields are presented and inversion techniques for the Boltzmann operator are discussed. Errors generated by the straight ahead approximation are derived and are shown to be negligible for most problems of space-radiation protection. A decoupling of projectile propagation from the target fields greatly simplifies the Boltzmann equations and allows an analytic solution of the target fragment transport. Analytic and numerical methods of solving the projectile transport equations are discussed. The nuclear physics underlying the coefficients in the Boltzmann equation is discussed. A coupled-channel optical model is found as a consequence of the loose binding of nuclear matter and closure of the nuclear states in high-energy reactions. Transport solutions with the developed data base are used with laboratory experiments to validate both the transport code and the data base. Numerical benchmarks and comparison with Monte Carlo calculations are also used for code validation.

Wilson, John W.↗
DOE OSTI2011

Light bullets in the spatiotemporal nonlinear Schroedinger equation with a variable negative diffraction coefficient

We report approximate analytical solutions to the (3+1)-dimensional spatiotemporal nonlinear Schroedinger equation, with the uniform self-focusing nonlinearity and a variable negative radial diffraction coefficient, in the form of three-dimensional solitons. The model may be realized in artificial optical media, such as left-handed materials and photonic crystals, with the anomalous sign of the group-velocity dispersion (GVD). The same setting may be realized through the interplay of the self-defocusing nonlinearity, normal GVD, and positive variable diffraction. The Hartree approximation is utilized to achieve a suitable separation of variables in the model. Then, an inverse procedure is introduced, with the aim to select a suitable profile of the modulated diffraction coefficient supporting desirable soliton solutions (such as dromions, single- and multilayer rings, and multisoliton clusters). The validity of the analytical approximation and stability of the solutions is tested by means of direct simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2011

Finding a spherically symmetric cosmology from observations in observational coordinates — advantages and challenges

One of the continuing challenges in cosmology has been to determine the large-scale space-time metric from observations with a minimum of assumptions — without, for instance, assuming that the universe is almost Friedmann-Lemaître-Robertson-Walker (FLRW). If we are lucky enough this would be a way of demonstrating that our universe is FLRW, instead of presupposing it or simply showing that the observations are consistent with FLRW. Showing how to do this within the more general spherically symmetric, inhomogeneous space-time framework takes us a long way towards fulfilling this goal. In recent work researchers have shown how this can be done both in the traditional Lemaître-Tolman-Bondi (LTB) 3 + 1 coordinate framework, and in the observational coordinate (OC) framework, in which the radial coordinate y is null (light-like) and measured down the past light cone of the observer. In this paper we investigate the stability of solutions, and the use of data in the OC field equations including their time evolution — i.e. our procedure is not restricted to our past light cone — and compare both approaches with respect to the singularity problem at the maximum of the angular-diameter distance, the stability of solutions, and the use of data in the field equations. We also compare the two approaches with regard to determining the cosmological constant Λ. This allows a more detailed account and assessment of the OC integration procedure, and enables a comparison of the relative advantages of the two equivalent solution frameworks. Both formulations and integration procedures should, in principle, lead to the same results. However, as we show in this paper, the OC procedure manifests certain advantages, particularly in the avoidance of coordinate singularities at the maximum of the angular-diameter distance, and in the stability of the solutions obtained. This particular feature is what allows us to do the best fitting of the data to smooth data functions and the possibility of constructing analytic solutions to the field equations. Smoothed data functions enable us to include properties that data must have within the model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Molecular dynamic simulation of Ar-Kr mixture across a rough walled nanochannel: Velocity and temperature profiles

This paper presents the results from a molecular dynamics simulation of mixture of argon and krypton in the Poiseuille flow across a rough walled nanochannel. The roughness effect on liquid nanoflows has recently drawn attention The computational software used for carrying out the molecular dynamics simulations is LAMMPS. The fluid flow takes place between two parallel plates and is bounded by horizontal rough walls in one direction and periodic boundary conditions are imposed in the other two directions. Each fluid atom interacts with other fluid atoms and wall atoms through Leenard-Jones (LJ) potential with a cut off distance of 5.0. To derive the flow a constant force is applied whose value is varied from 0.1 to 0.3 and velocity profiles and temperature profiles are noted for these values of forces. The velocity profile and temperature profiles are also looked at different channel widths of nanochannel and at different densities of mixture. The velocity profile and temperature profile of rough walled nanochannel are compared with that of smooth walled nanochannel and it is concluded that mean velocity increases with increase in channel width, force applied and decrease in density also with introduction of roughness in the walls of nanochannel mean velocity again increases and results also agree with the analytical solution of a Poiseuille flow.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A note on the cylindrical solitary waves in an electron-acoustic plasma with vortex electron distribution

In the present work, we consider the propagation of nonlinear electron-acoustic non-planar waves in a plasma composed of a cold electron fluid, hot electrons obeying a trapped/vortex-like distribution, and stationary ions. The basic nonlinear equations of the above described plasma are re-examined in the cylindrical coordinates through the use reductive perturbation method in the long-wave approximation. The modified cylindrical Korteweg-de Vries equation with fractional power nonlinearity is obtained as the evolution equation. Due to the nature of nonlinearity, which is fractional, this evolution equation cannot be reduced to the conventional Korteweg–de Vries equation. An analytical solution to the evolution equation, by use of the method developed by Demiray [Appl. Math. Comput. 132, 643 (2002); Comput. Math. Appl. 60, 1747 (2010)] and a numerical solution by employing a spectral scheme are presented and the results are depicted in a figure. The numerical results reveal that both solutions are in good agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Geometrodynamics of a thin shell in the Reissner-Nordström metric

The geometrodynamics of a thin, dust, electrically neutral shell in the metric of a charged Reissner-Nordström black hole is considered. The Arnowitt-Deser-Misner formalism is used to construct the thin-shell Hamiltonian. The wave equation is derived. The wave equation is shown to be a second-order homogeneous difference equation. Exact analytical solutions have been found.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Nonlinear relativistic quantum theory of Cherenkov emission of longitudinal Langmuir waves in a plasma

A nonlinear relativistic quantum theory of stimulated Cherenkov emission of longitudinal waves by a relativistic monoenergetic electron beam in a cold isotropic plasma is presented. The theory makes use of a quantum model based on the Klein-Gordon equation. The instability growth rates are obtained in the linear approximation and are shown to go over to the familiar growth rates in the classical limit. The mechanisms for the nonlinear saturation of relativistic Cherenkov beam instabilities are described with allowance for quantum effects, and the corresponding analytic solutions are derived.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Stability analysis and future singularity of the m{sup 2} R □{sup -2} R model of non-local gravity

We analyse the classical stability of the model proposed by Maggiore and Mancarella, where gravity is modified by a term ∼ m{sup 2} R □{sup -2} R to produce the late-time acceleration of the expansion of the universe. Our study takes into account all excitations of the metric that can potentially drive an instability. There are some subtleties in identifying these modes, as a non-local field theory contains dynamical fields which yet do not correspond to degrees of freedom. Since some of them are ghost-like, we clarify the impact of such modes on the stability of the solutions of interest that are the flat space-time and cosmological solutions. We then find that flat space-time is unstable under scalar perturbations, but the instability manifests itself only at cosmological scales, i.e. out of the region of validity of this solution. It is therefore the stability of the FLRW solution which is relevant there, in which case the scalar perturbations are known to be well-behaved by numerical studies. By finding the analytic solution for the late-time behaviour of the scale factor, which leads to a big rip singularity, we argue that the linear perturbations are bounded in the future because of the domination of Hubble friction. In particular, this effect damps the scalar ghost perturbations which were responsible for destabilizing Minkowski space-time. Thus, the model remains phenomenologically viable.

79 ASTRONOMY AND ASTROPHYSICS↗