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

THCM Coupled Model for Hydrate-Bearing Sediments: Data Analysis and Design of New Field Experiments (Marine and Permafrost Settings)

Gas hydrates are solid compounds made of water molecules clustered around low molecular weight gas molecules such as methane, hydrogen, and carbon dioxide. Methane hydrates form under pressure (P) and temperature (T) conditions that are common in sub-permafrost layers and in deep marine sediments. Stability conditions constrain the occurrence of gas hydrates to submarine sediments and permafrost regions. The amount of technically recoverable methane trapped in gas hydrate may exceed 104tcf. Gas hydrates are a potential energy resource, can contribute to climate change, and can cause large-scale seafloor instabilities. In addition, hydrate formation can be used for CO2 sequestration (also through CO2-CH4 replacement), and efficient geological storage seals. The experimental study of hydrate bearing sediments has been hindered by the very low solubility of methane in water (lab testing), and inherent sampling difficulties associated with depressurization and thermal changes during core extraction. This situation has prompted more decisive developments in numerical modeling in order to advance the current understanding of hydrate bearing sediments, and to investigate/optimize production strategies and implications. The goals of this research has been to addresses the complex thermo-hydro-chemo-mechanical THCM coupled phenomena in hydrate-bearing sediments, using a truly coupled numerical model that incorporates sound and proven constitutive relations, satisfies fundamental conservation principles. Analytical solutions aimed at verifying the proposed code have been proposed as well. These tools will allow to better analyze available data and to further enhance the current understanding of hydrate bearing sediments in view of future field experiments and the development of production technology.

58 GEOSCIENCES↗
DOE OSTI2016

Benchmark for Free-Gas Elastic Scattering with the Heavy-Gas Model and Piecewise-Constant Cross Sections

Accurate Monte Carlo calculations involving thermal neutrons require that codes implement the thermal scattering physics correctly. The heavy-gas model, which is valid for large atomic mass, with constant cross sections in energy has an analytic solution for the thermal neutron spectrum that can be readily used to verify the implementation of the free-gas thermal scattering routines. Previous work has shown that MCNP6.1 correctly computes the free-gas thermal flux with constant cross sections. This summary extends the verification to piecewise-constant cross sections in energy that can better approximate physically realistic cross sections. The obtained solutions for the piecewise-constant cases agree with the respective thermal spectra obtained by MCNP6.1, further verifying that it is accurately simulating the free-gas model. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Observability of localized magnetoplasmons in quantum dots: Scrutinizing the eligibility of far-infrared, Raman, and electron-energy-loss spectroscopies

We investigate a one-component, quasi-zero dimensional, quantum plasma exposed to a parabolic potential and an applied magnetic field in the symmetric gauge. If the size of such a system as can be realized in the semiconducting quantum dots is on the order of the de-Broglie wavelength, the electronic and optical properties become highly tunable. Then the quantum size effects challenge the observation of many-particle phenomena such as the magneto-optical absorption, Raman intensity, and electron energy-loss spectrum. An exact analytical solution of the problem leads us to infer that these many-particle phenomena are, in fact, dictated by the generalized Kohn’s theorem (GKT) in the long-wavelength limit. Maneuvering the confinement and/or the magnetic field furnishes the resonance energies capable of being explored with the FIR, Raman, and/or electron-energy-loss spectroscopy. This implies that either of these probes is competent in observing the localized magnetoplasmons in the system. As an application of the rigorous analytical diagnosis of the system, we have presented various pertinent single-particle, such as Fock-Darwin spectrum, Fermi energy, zigzag excitation spectrum, and magneto-optical transitions, and the many-particle phenomena, such as magneto-optical absorption, Raman intensity, and electron energy-loss probability. In the latter, the energy position of the resonance peaks is observed to be independent of the electron-electron interactions and hence of the number of electrons in the quantum dot in compliance with the GKT. It is found that both confinement potential and magnetic field play a decisive role in influencing the aforementioned many-particle phenomena. Specifically, increasing (decreasing) the strength of the confining potential is found to be analogous to shrinking (expanding) the size of the quantum dots and results into a blue (red) shift in the respective spectra. Intensifying the magnetic field has two-fold effects in the resonance spectra associated with all three spectroscopies considered here: the lower resonance peak observes a red shift, whereas the higher one experiences a blue shift. This is a unique and intriguing behavior observed in the quantum dots with complete confinement. A deeper insight into the physics of the quantum dots is paving the way for their implementation in such diverse fields as quantum computing and medical imaging.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2011

Relativistic solitary waves with phase modulation embedded in long laser pulses in plasmas

We investigate the existence of nonlinear phase-modulated relativistic solitary waves embedded in an infinitely long circularly polarized electromagnetic wave propagating through a plasma. These states are exact nonlinear solutions of the 1-dimensional Maxwell-fluid model for a cold plasma composed of electrons and ions. The solitary wave, which consists of an electromagnetic wave trapped in a self-generated Langmuir wave, presents a phase modulation when the group velocity V and the phase velocity V{sub ph} of the long circularly polarized electromagnetic wave do not match the condition VV{sub ph} = c{sup 2}. The main properties of the waves as a function of their group velocities, wavevectors, and frequencies are studied, as well as bifurcations of the dynamical system that describes the waves when the parameter controlling the phase modulation changes from zero to a finite value. Such a transition is illustrated in the limit of small amplitude waves where an analytical solution for a grey solitary wave exists. The solutions are interpreted as the stationary state after the collision of a long laser pulse with an isolated solitary wave.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Kinklike structures in an arcsin real scalar dynamics

Highlights: • Kinklike analyses in a real scalar theory with an arcsin dynamic. • First-order formalism implementation and deformation procedure in an arcsin scalar dynamics. In this paper, we analyze kink-like analytical solutions in a real scalar theory with an arcsin dynamics inspired by the arcsin electrodynamics presented in Kruglov (2015). This analysis is done by means of the first-order formalism. This formalism provides a framework where the equations of motion can be simplified by preserving the linear stability of the theory. In this work, the deformation procedure is implemented with the aim of finding exact solutions in systems with generalized dynamics. Along the paper, we explore how the first-order formalism is implemented in the arcsin kinetics and how such a term influences the kink-like solutions. As a part of the result of our paper, we show that the kink-like solutions are similar to the ones obtained in the standard scalar kinetic theory. We also show that the extra parameter, that controls the non-linearities of the model, plays an essential role in the energy densities and stability potentials. These quantities vary according to this parameter. The goals here are to show how the first-order framework is implemented in this arcsin scenario and to present the analytical kink-like solutions that can be found by means of the first-order framework and deformation method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Motion of a virtual cathode in a cylindrical channel with electron beam transport in the “compressed” state

This paper studies the motion of a virtual cathode in a two-section drift tube with the formation and breakup of the “compressed” state of an electron beam. Experimental arrangements to intercept part of the injected current during the voltage pulse and to provide virtual cathode motion toward the collector are proposed. The arrangements were implemented on the SINUS-7 high-current electron accelerator. Theoretical and experimental dependences of the virtual cathode velocity on the injected current and cathode voltage are presented. The experimental data on virtual cathode motion agree with its theoretical model based on analytical solutions of equations assisted by computer simulation with the PIC code KARAT. The results of the work demonstrate the feasibility of controlling the virtual cathode motion which can be used in collective ion acceleration and microwave generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2010

The wavefunction of a gravitating shell

We have calculated a discrete spectrum and found an exact analytical solution in the form of Meixner polynomials for the wavefunction of a thin gravitating shell in the Reissner-Nordstroem geometry. We show that there is no extreme state in the quantum spectrum of the gravitating shell, as in the case of an extreme black hole.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Atomic density functional and diagram of structures in the phase field crystal model

The phase field crystal model provides a continual description of the atomic density over the diffusion time of reactions. We consider a homogeneous structure (liquid) and a perfect periodic crystal, which are constructed from the one-mode approximation of the phase field crystal model. A diagram of 2D structures is constructed from the analytic solutions of the model using atomic density functionals. The diagram predicts equilibrium atomic configurations for transitions from the metastable state and includes the domains of existence of homogeneous, triangular, and striped structures corresponding to a liquid, a body-centered cubic crystal, and a longitudinal cross section of cylindrical tubes. The method developed here is employed for constructing the diagram for the homogeneous liquid phase and the body-centered iron lattice. The expression for the free energy is derived analytically from density functional theory. The specific features of approximating the phase field crystal model are compared with the approximations and conclusions of the weak crystallization and 2D melting theories.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

STEADY GENERAL RELATIVISTIC MAGNETOHYDRODYNAMIC INFLOW/OUTFLOW SOLUTION ALONG LARGE-SCALE MAGNETIC FIELDS THAT THREAD A ROTATING BLACK HOLE

General relativistic magnetohydrodynamic (GRMHD) flows along magnetic fields threading a black hole can be divided into inflow and outflow parts, according to the result of the competition between the black hole gravity and magneto-centrifugal forces along the field line. Here we present the first self-consistent, semi-analytical solution for a cold, Poynting flux–dominated (PFD) GRMHD flow, which passes all four critical (inner and outer, Alfvén, and fast magnetosonic) points along a parabolic streamline. By assuming that the dominating (electromagnetic) component of the energy flux per flux tube is conserved at the surface where the inflow and outflow are separated, the outflow part of the solution can be constrained by the inflow part. The semi-analytical method can provide fiducial and complementary solutions for GRMHD simulations around the rotating black hole, given that the black hole spin, global streamline, and magnetizaion (i.e., a mass loading at the inflow/outflow separation) are prescribed. For reference, we demonstrate a self-consistent result with the work by McKinney in a quantitative level.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

A realistic renormalizable supersymmetric E₆ model

A complete realistic model based on the supersymmetric version of E₆ is presented. It consists of three copies of matter 27, and a Higgs sector made of 2×(27+27⁻)+351´+351´⁻ representations. An analytic solution to the equations of motion is found which spontaneously breaks the gauge group into the Standard Model. The light fermion mass matrices are written down explicitly as non-linear functions of three Yukawa matrices. This contribution is based on Ref. [1].

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

DIRECT INTEGRATION OF THE COLLISIONLESS BOLTZMANN EQUATION IN SIX-DIMENSIONAL PHASE SPACE: SELF-GRAVITATING SYSTEMS

We present a scheme for numerical simulations of collisionless self-gravitating systems which directly integrates the Vlasov-Poisson equations in six-dimensional phase space. Using the results from a suite of large-scale numerical simulations, we demonstrate that the present scheme can simulate collisionless self-gravitating systems properly. The integration scheme is based on the positive flux conservation method recently developed in plasma physics. We test the accuracy of our code by performing several test calculations, including the stability of King spheres, the gravitational instability, and the Landau damping. We show that the mass and the energy are accurately conserved for all the test cases we study. The results are in good agreement with linear theory predictions and/or analytic solutions. The distribution function keeps the property of positivity and remains non-oscillatory. The largest simulations are run on 64{sup 6} grids. The computation speed scales well with the number of processors, and thus our code performs efficiently on massively parallel supercomputers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Thermal plasma irradiation under the action of the microwave-frequency-modulated laser beam

The results of the calculation of the thermal irradiation of the laser plasma formed by a powerful laser beam with the microwave-frequency-modulated intensity are presented. The analytical solution has been obtained for the case of the light detonation regime. It has been shown that the modulation of the gasdynamic parameters due to the absorption of the laser radiation leads to the modulation of the spectral and integral brightness observed from the thermal plasma irradiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Element free Galerkin formulation of composite beam with longitudinal slip

Behaviour between two materials in composite beam is assumed partially interact when longitudinal slip at its interfacial surfaces is considered. Commonly analysed by the mesh-based formulation, this study used meshless formulation known as Element Free Galerkin (EFG) method in the beam partial interaction analysis, numerically. As meshless formulation implies that the problem domain is discretised only by nodes, the EFG method is based on Moving Least Square (MLS) approach for shape functions formulation with its weak form is developed using variational method. The essential boundary conditions are enforced by Langrange multipliers. The proposed EFG formulation gives comparable results, after been verified by analytical solution, thus signify its application in partial interaction problems. Based on numerical test results, the Cubic Spline and Quartic Spline weight functions yield better accuracy for the EFG formulation, compares to other proposed weight functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2010

Electromagnetic bound states in the radiation continuum for periodic double arrays of subwavelength dielectric cylinders

Electromagnetic bound states in the radiation continuum are studied for periodic double arrays of subwavelength dielectric cylinders in TM polarization. They are similar to localized waveguide mode solutions of Maxwell's equations for metal cavities or defects of photonic crystals, but, in contrast to the latter, their spectrum lies in the radiation continuum. The phenomenon is identical to the existence of bound states in the radiation continuum in quantum mechanics, discovered by von Neumann and Wigner. In the formal scattering theory, these states appear as resonances with the vanishing width. For the system studied, the bound states are shown to exist at specific distances between the arrays in the spectral region where one or two diffraction channels are open. Analytic solutions are obtained for all bound states (below the radiation continuum and in it) in the limit of thin cylinders (the cylinder radius is much smaller than the wavelength). The existence of bound states is also established in the spectral region where three and more diffraction channels are open, provided the dielectric constant and radius of the cylinders are fine-tuned. The near field and scattering resonances of the structure are investigated when the distance between the arrays varies in a neighborhood of its critical values at which the bound states are formed. In particular, it is shown that the near field in the scattering process becomes significantly amplified in specific regions of the array as the distance approaches its critical values. The effect may be used to control optical nonlinear effects by varying the distance between the arrays near its critical values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2010

Numerical evolution of multiple black holes with accurate initial data

We present numerical evolutions of three equal-mass black holes using the moving puncture approach. We calculate puncture initial data for three black holes solving the constraint equations by means of a high-order multigrid elliptic solver. Using these initial data, we show the results for three black hole evolutions with sixth-order waveform convergence. We compare results obtained with the BAM and AMSS-NCKU codes with previous results. The approximate analytic solution to the Hamiltonian constraint used in previous simulations of three black holes leads to different dynamics and waveforms. We present some numerical experiments showing the evolution of four black holes and the resulting gravitational waveform.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2010

Lemaitre-Tolman-Bondi cosmological wormhole

We present a new analytical solution of the Einstein field equations describing a wormhole shell of zero thickness joining two Lemaitre-Tolman-Bondi universes, with no radial accretion. The material on the shell satisfies the energy conditions and, at late times, the shell becomes comoving with the dust-dominated cosmic substratum.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

SU-E-T-508: Internal Organ Motion Effect On Radiation Dose to a Point Under Half-Beam Block Match Line

Purpose Half-beam block is a field matching technique frequently used in radiotherapy. With no setup error, a well calibrated linac, and no internal organ motion, two photon fields can be matched seamlessly dosimetry-wise with their central axes passing the match line. However, in actual clinical situations, internal organ motion is often inevitable. This study was conducted to investigate its influence on radiation dose to patient internal points directly under the matching line. Methods A clinical setting is modeled as two half-space (x<0 and x<0) radiation fields that are turned on sequentially with a time gap of integer times of the patient internal organ motion period (T{sub 0}). Our point of interest moves with patient internal organs periodically and evenly in and out of the radiation fields, resulting in an average location at x=0. When the fields are delivered without any motion management, the initial phase of the point’s movement is unknown. Statistical methods are used to compute the expected value ( ) and variance (σ) of the point dose given the uncertainty. Results Analytical solutions are obtained for and s of dose received by a point directly under the match line. is proportional to the total beam-on time (T1), and σ demonstrates previously unknown periodic behavior. /<Dσ has local maximums (minimums) when T{sub 1} is 2n*T{sub 0} ((2n+1)*T{sub 0}) (n is an integer) with an upper envelope scales as 1/T{sub 1}. Worst case scenario the point dose can be zero or double of the expected value. Conclusion We have analytically analyzed the internal organ motion effect on radiation dose received by a patient internal point directly under a half-beam block match line. Our results help us to better understand this phenomenon and facilitate the reduction of point dosimetric uncertainties in our clinical practice.

60 APPLIED LIFE SCIENCES↗