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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

A Hydro-mechanical Model and Analytical Solutions for Geomechanical Modeling of Carbon Dioxide Geological Sequestration

We present a hydro-mechanical model for geological sequestration of carbon dioxide. The model considers the poroelastic effects by taking into account the coupling between the geomechanical response and the fluid flow in greater detail. The simplified hydro-mechanical model includes the geomechanical part that relies on the linear elasticity, while the fluid flow is based on the Darcy’s law. Two parts were coupled using the standard linear poroelasticity. Analytical solutions for pressure field were obtained for a typical geological sequestration scenario. The model predicts the temporal and spatial variation of pressure field and effects of permeability and elastic modulus of formation on the fluid pressure distribution.

sequestration, geomechanics, Multi-scale, poro-ela↗
DOE OSTI2025

Integrating Analytical Solutions and U-Net Model for Predicting Groundwater Contaminant Plumes in Pump-and-Treat Systems

Pump-and-treat (P&T) is a common technique for groundwater remediation involving the extraction and treatment of contaminated water above ground. Optimizing the design and operation of the P&T well network is essential for maximizing the system’s effectiveness and efficiency. However, this optimization often necessitates many model evaluations, leading to computationally demanding tasks. This study introduces a novel approach that integrates analytical solutions for groundwater dynamics with the U-Net (Ronneberger et al., 2015) deep learning framework to predict groundwater contaminant plume migration under dynamic pumping conditions. By incorporating the Thiem equation (Thiem, 1906) into the input preprocessing, the U-Net model transforms sparse well data into a continuous spatial field that captures the hydraulic impacts of pumping activities. This integration enables the model to leverage both deep learning capabilities and classical physics-based groundwater theories, enhancing prediction accuracy and computational efficiency. These advancements can facilitate rapid, large-scale evaluations of P&T optimization simulations, allowing for timely and effective decision-making in well placement and system management. We demonstrate the model's robust performance across both simplified transient 2D models and a more complex 3D heterogeneous site model at the 200 West P&T facility at the Hanford Site. The U-Net-based model offers substantial computational advantages, reducing simulation times significantly compared to full physics-based models and providing a powerful tool for rapid site evaluation and P&T system optimization, such as evaluating alternative P&T well network designs. Our findings highlight the potential of advanced machine learning models to significantly enhance the efficiency and sustainability of groundwater remediation efforts, offering a novel application of U-Net architecture in environmental science.

Pump-and-treat↗
DOE OSTI2022

Leakage From Coexisting Geologic Forcing and Injection‐Induced Pressurization: A Semi‐Analytical Solution for Multilayered Aquifers With Multiple Wells

Abstract Abnormal fluid pressures (above or below hydrostatic pressure) can develop and persist in sedimentary basins. The common occurrence of abnormal pressures may cause challenges for project permitting of geological carbon sequestration (GCS), particularly in reservoirs with pre‐injection overpressure. The leaky wells that may exist in some sedimentary basins can provide flow paths between deep brine aquifers and shallower freshwater aquifers. Pre‐injection relative overpressures can cause brine leakage through leaky wells even before any injection occurs. The tendency for flow through leaky wells is coupled with the process of pressure dissipation that occurs through aquitards. Specifically, with non‐zero permeability, aquitards can dissipate pressure over large areal extents and thereby reduce leakage rates through leaky wells. This study presents development of a semi‐analytical solution for hydraulic head and brine leakage in multilayered aquifer–aquitard systems with geologic pressure forcing. The geologic forcing that causes abnormal pressures in the multilayered system can coexist with any number of injection, extraction, and leaky wells that also affect fluid pressure. The semi‐analytical model is applied to explore how leakage through leaky wells varies as functions of pressurization rate, along with aquitard and leaky well properties in an overpressured multilayered system. The results show that although injection‐induced pressures can dissipate rapidly through suitably permeable aquitards, coexisting geologic forcing may create sustained rates of brine leakage into freshwater aquifers through leaky wells. In GCS, a very low‐permeability aquitard with high capillary entry pressure to free‐phase CO 2 is desired to serve as the caprock to prevent leakage of CO 2 from the storage reservoir. Nevertheless, the results from this study show that the brine leakage impact to shallow freshwater aquifers through leaky wells might substantially decrease with increasing aquitard permeability values, as long as small aquitard permeability and high capillary entry pressure serve to prevent CO 2 leakage.

Cihan, A.↗
DOE OSTI2018

Verification of Modelica-Based Models with Analytical Solutions for Tritium Diffusion

Here, tritium transport in metal and molten salt fluids combined with diffusion through high-temperature structural materials is an important phenomenon in both magnetic confinement fusion (MCF) and molten salt reactor (MSR) applications. For MCF, tritium is desirable to capture for fusion fuel. For MSRs, uncaptured tritium potentially can be released to the environment. In either application, quantifying the time- and space-dependent tritium concentration in the working fluid(s) and structural components is necessary.Whereas capability exists specifically for calculating tritium transport in such systems (e.g., using TMAP for fusion reactors), it is desirable to unify the calculation of tritium transport with other system variables such as dynamic fluid and structure temperature combined with control systems such as those that might be found in a system code. Some capability for radioactive trace substance transport exists in thermal-hydraulic systems codes (e.g., RELAP5-3D); however, this capability is not coupled to species diffusion through solids. Combined calculations of tritium transport and thermal-hydraulic solution have been demonstrated with TRIDENT but only for a specific type of MSR.Researchers at Oak Ridge National Laboratory have developed a set of Modelica-based dynamic system modeling tools called TRANsient Simulation Framework Of Reconfigurable Models (TRANSFORM) that were used previously to model advanced fission reactors and associated systems. In this system, the augmented TRANSFORM library includes dynamically coupled fluid and solid trace substance transport and diffusion. Results from simulations are compared against analytical solutions for verification.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

Analytical Solution and Physics of a Propellant Damping Device

NASA design teams have been investigating options for "detuning" Ares I to prevent oscillations originating in the vehicle solid-rocket main stage from synching up with the natural resonance of the rest of the vehicle. An experimental work started at NASA MSFC center in 2008 using a damping device showed great promise in damping the vibration level of an 8 resonant tank. However, the mechanisms of the vibration damping were not well understood and there were many unknowns such as the physics, scalability, technology readiness level (TRL), and applicability for the Ares I vehicle. The objectives of this study are to understand the physics of intriguing slosh damping observed in the experiments, to further validate a Computational Fluid Dynamics (CFD) software in propellant sloshing against experiments with water, and to study the applicability and efficiency of the slosh damper to a full scale propellant tank and to cryogenic fluids. First a 2D fluid-structure interaction model is built to model the system resonance of liquid sloshing and structure vibration. A damper is then added into the above model to simulate experimentally observed system damping phenomena. Qualitative agreement is found. An analytical solution is then derived from the Newtonian dynamics for the thrust oscillation damper frequency, and a slave mass concept is introduced in deriving the damper and tank interaction dynamics. The paper will elucidate the fundamental physics behind the LOX damper success from the derivation of the above analytical equation of the lumped Newtonian dynamics. Discussion of simulation results using high fidelity multi-phase, multi-physics, fully coupled CFD structure interaction model will show why the LOX damper is unique and superior compared to other proposed mitigation techniques.

Yang, H. Q.↗
DOE OSTI2017

Phase field model of fluid-driven fracture in elastic media: Immersed-fracture formulation and validation with analytical solutions

Propagation of fluid-driven fractures plays an important role in natural and engineering processes, including transport of magma in the lithosphere, geologic sequestration of carbon dioxide, and oil and gas recovery from low-permeability formations, among many others. The simulation of fracture propagation poses a computational challenge as a result of the complex physics of fracture and the need to capture disparate length scales. Phase field models represent fractures as a diffuse interface and enjoy the advantage that fracture nucleation, propagation, branching, or twisting can be simulated without ad hoc computational strategies like remeshing or local enrichment of the solution space. Here we propose a new quasi-static phase field formulation for modeling fluid-driven fracturing in elastic media at small strains. The approach fully couples the fluid flow in the fracture (described via the Reynolds lubrication approximation) and the deformation of the surrounding medium. The flow is solved on a lower dimensionality mesh immersed in the elastic medium. This approach leads to accurate coupling of both physics. We assessed the performance of the model extensively by comparing results for the evolution of fracture length, aperture, and fracture fluid pressure against analytical solutions under different fracture propagation regimes. Thus, the excellent performance of the numerical model in all regimes builds confidence in the applicability of phase field approaches to simulate fluid-driven fracture.

58 GEOSCIENCES↗
NASA NTRS2015

An Analytical Solution for Yaw Maneuver Optimization on the International Space Station and Other Orbiting Space Vehicles

This paper presents a new method for optimizing yaw maneuvers, which are the most common large maneuvers on the International Space Station (ISS). The goal of the maneuver optimization is to find a maneuver trajectory with minimal torques acting on the vehicle during the maneuver. Therefore, the thruster firings necessary to perform the maneuver are minimized. Reduction of thruster firings saves propellant and decreases structural loads and contamination of the vehicle critical elements, thus saving the service life of the thrusters and the vehicle itself. Equations describing the pitch and roll motion needed to counteract the major torques during a yaw maneuver are obtained. Also, a yaw rate profile is suggested. In the obtained optimized case, the torques are significantly reduced. The proposed approximate analytical solution does not require extensive computer resources and, therefore, can be implemented using software onboard the ISS. As a result, the maneuver execution will be automatic. This is one of the major benefits of the simplified solution presented in this paper with respect to existing computational approaches. The suggested maneuver optimization method can be used not only for the ISS, but for other space vehicles as well.

Dobrinskaya, Tatiana↗
NASA NTRS2004

Analytical Solution for the Critical Velocity of Pushing/Engulfment Transition

The distribution of ceramic particles in a metal matrix composite material depends primarily on the interaction of the particles with the solid/liquid interface during the solidification process. A numerical model that describes the evolution of the shape of the solid/liquid interface in the proximity of a foreign particle will presented in this paper. The model accounts for the influence of the temperature gradient and the Gibbs-Thomson and disjoining pressure effects. It shows that for the systems characterized by k(sub p) < k(sub L) the disjoining pressure causes the interface curvature to change its sign in the close-contact particle/interface region. It also shows that the increase of the temperature gradient diminishes the effect of the disjoining pressure. The analysis of the numerical results obtained for a large range of processing conditions and materials parameters has led to the development of an analytical solution for the critical velocity of pushing/engulfinent transition. The theoretical results will be discussed and compared with the experimental measurements performed under microgravity conditions.

Catalina, Adrian V.↗
DOE OSTI2014

Analytical solutions to matrix diffusion problems

We report an analytical method to solve in a few cases of practical interest the equations which have traditionally been proposed for the matrix diffusion problem. In matrix diffusion, elements dissolved in ground water can penetrate the porous rock surronuding the advective flow paths. In the context of radioactive waste repositories this phenomenon provides a mechanism by which the area of rock surface in contact with advecting elements is greatly enhanced, and can thus be an important delay mechanism. The cases solved are relevant for laboratory as well for in situ experiments. Solutions are given as integral representations well suited for easy numerical solution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1975

Analytic solution for low thrust out of ecliptic mission design

Description of an analytic trajectory solution for low thrust out of ecliptic mission design based on the two-variable asymptotic expansion technique. A practical steering program utilizing normal thrusting symmetrical about the nodes is assumed. A two-term expansion yields an essentially exact solution of the problem which is in a form suitable for evaluation on small desk-top computers. It is shown that for a circular base orbit at 1 AU, a 15-kw solar electric system with a specific impulse of 3000 sec is capable of delivering an 1100-kg net payload to a heliographic latitude of about 63 deg.

Flandro, G. A.↗
NASA NTRS1964

Analytical Solutions for Non-Linear Differential Equations with the Help of a Digital Computer

A technique was developed with the help of a digital computer for analytic (algebraic) solutions of autonomous and nonautonomous equations. Two operational transform techniques have been programmed for the solution of these equations. Only relatively simple nonlinear differential equations have been considered. In the cases considered it has been possible to assimilate the secular terms into the solutions. For cases where f(t) is not a bounded function, a direct series solution is developed which can be shown to be an analytic function. All solutions have been checked against results obtained by numerical integration for given initial conditions and constants. It is evident that certain nonlinear differential equations can be solved with the help of a digital computer.

Cromwell, P. C.↗
DOE OSTI2014

SU-C-9A-04: Alternative Analytic Solution to the Paralyzable Detector Model to Calculate Deadtime and Deadtime Loss

Purpose: Some common methods to solve for deadtime are (1) dual-source method, which assumes two equal activities; (2) model fitting, which requires multiple acquisitions as source decays; and (3) lossless model, which assumes no deadtime loss at low count rates. We propose a new analytic alternative solution to calculate deadtime for paralyzable gamma camera. Methods: Deadtime T can be calculated analytically from two distinct observed count rates M1 and M2 when the ratio of the true count rates alpha=N2/N1 is known. Alpha can be measured as a ratio of two measured activities using dose calibrators or via radioactive decay. Knowledge of alpha creates a system with 2 equations and 2 unknowns, i.e., T and N1. To verify the validity of the proposed method, projections of a non-uniform phantom (4GBq 99mTc) were acquired in using Siemens SymbiaS multiple times over 48 hours. Each projection has >100kcts. The deadtime for each projection was calculated by fitting the data to a paralyzable model and also by using the proposed 2-acquisition method. The two estimates of deadtime were compared using the Bland-Altmann method. In addition, the dependency of uncertainty in T on uncertainty in alpha was investigated for several imaging conditions. Results: The results strongly suggest that the 2-acquisition method is equivalent to the fitting method. The Bland-Altman analysis yielded mean difference in deadtime estimate of ∼0.076us (95%CI: -0.049us, 0.103us) between the 2-acquisition and model fitting methods. The 95% limits of agreement were calculated to be -0.104 to 0.256us. The uncertainty in deadtime calculated using the proposed method is highly dependent on the uncertainty in the ratio alpha. Conclusion: The 2-acquisition method was found to be equivalent to the parameter fitting method. The proposed method offers a simpler and more practical way to analytically solve for a paralyzable detector deadtime, especially during physics testing.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2017

Isothermal Bondi Accretion in Jaffe and Hernquist Galaxies with a Central Black Hole: Fully Analytical Solutions

One of the most active fields of research of modern-day astrophysics is that of massive black hole formation and coevolution with the host galaxy. In these investigations, ranging from cosmological simulations, to semi-analytical modeling, to observational studies, the Bondi solution for accretion on a central point-mass is widely adopted. In this work we generalize the classical Bondi accretion theory to take into account the effects of the gravitational potential of the host galaxy, and of radiation pressure in the optically thin limit. Then, we present the fully analytical solution, in terms of the Lambert–Euler W -function, for isothermal accretion in Jaffe and Hernquist galaxies with a central black hole. The flow structure is found to be sensitive to the shape of the mass profile of the host galaxy. These results and the formulae that are provided, most importantly, the one for the critical accretion parameter, allow for a direct evaluation of all flow properties, and are then useful for the abovementioned studies. As an application, we examine the departure from the true mass accretion rate of estimates obtained using the gas properties at various distances from the black hole, under the hypothesis of classical Bondi accretion. An overestimate is obtained from regions close to the black hole, and an underestimate outside a few Bondi radii; the exact position of the transition between the two kinds of departure depends on the galaxy model.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1988

Analytical solutions of a class of multidimensional Fokker-Planck equations

Analytical representations of the solutions of Kolmogorov's forward (i.e., Fokker-Planck) and backward partial differential equations are obtained by employing Girsanov's (1960) transformation. The equations considered are restricted to those with the drift vector being the gradient of a function and the diffusion matrix being nonsingular and dependent upon time only.

Zhang, Weijian↗
DOE OSTI2011

Analytical solution of the radiative transfer equation for infinite-space fluence

This Brief Report presents the derivation of analytical expressions for the fluence of the steady state radiative transfer equation in an infinitely extended and anisotropically scattering medium in arbitrary dimensions for different source types. The fluence, which is composed of an infinite sum of diffusion-like Green's functions, was compared to the Monte Carlo method. Within the stochastic nature of the Monte Carlo simulations, an exact agreement was found in the steady state and time domains. It is shown that the use of low-order approximations is sufficient for many relevant cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Two-level system coupled to phonons: Full analytical solution

In this study, we propose an analytical procedure to fully solve a two-level system coupled to phonons. Instead of using the common formulation in terms of linear and quadratic system-phonon couplings, we introduce different phonons depending on the system electronic level. We use this approach to recover known results for the linear-coupling limit in a simple way. More importantly, we derive results for the quadratic coupling induced by a phonon frequency change, a problem considered up to now as not analytically solvable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1984

Approximate analytic solutions for singular non-linear oscillators

Mickens (1981, 1984) has considered analytic techniques for obtaining approximate solutions to one-dimensional nonlinear oscillatory systems x(double-dot) + x = lambda f(x, x/dot/, lambda) where lambda is a small positive parameter and f is a nonlinear polynomial function of its arguments. However, in certain cases there is interest in the analysis of physical systems for which the nonlinear function f(x, x/dot/, lambda) is singular for finite values of x or x(dot). The present investigation is concerned with the use of existing approximate analytic schemes to obtain solutions to singular nonlinear oscillatory differential equations.

Bota, K. B.↗
NASA NTRS1984

An analytic solution for the orbital perturbations of the Venus Radar Mapper due to gravitational harmonics

Hill's variational equations are solved analytically for the orbital perturbations of a spacecraft nominally in an elliptic orbit around a non-spherical body. The rotation of the central planet about its spin-axis is not considered in the analysis. The perturbations are restricted to the planetary gravitational harmonics only. An extremely simple algorithm is derived to transform the spherical harmonic potentials to the orbital coordinate system, and the resulting accelerations are shown to be simply trigonometric functions of the true anomaly. With the principal matrix solution for the differential equations of the adjoint system given in closed form, the orthogonality of the trigonometric functions makes it possible to obtain an analytic solution for the non-homogeneous problem, at intervals of 2 pi in true anomaly. The solution for orbital perturbations can be extended over several revolutions by applying well-known results from Floquet's theory. The technique is demonstrated with results presented on the spacecraft periapsis altitude for the forthcoming Venus Radar Mapper Mission.

Vijayaraghavan, A.↗