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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Uncertainty Quantification for the Reliability of the Analytical Analysis for the Simplified Model of CO2 Geological Sequestration

A hydro-mechanical model with analytical solutions including pressure evolution and geomechanical deformation for geological CO2 injection and sequestration were introduced in our previous work. However, the reliability and accuracy of the hydro-mechanical model and the companion analytical solution are uncertain because of the assumptions and simplifications in the analytical model, though it was validated by a few example cases. This study introduce the method to efficiently measure the accuracy of the analytical model, and specify the acceptable input parameters range that can guarantee the accuracy and reliability of the analytical solution. A coupled hydro-geomechanical subsurface transport simulator STOMP was adopted as a reference to justify the reliability of the hydro-mechanical model and the analytical solution. A quasi-Monte Carlo sampling method was applied to efficiently sample the input parameter space.

analytical analysis↗
DOE OSTI2016

Modeling crossed-beam energy transfer for inertial confinement fusion

We developed a numerical code that describes both the energy transfer occurring when two or more laser beams overlap in a weakly non-homogeneous plasma, and the beam energy losses associated with the electron-ion collisions. The numerical solutions are validated with both the exact analytical solutions in homogeneous plasmas, and with new approximate analytical solutions in non-homogeneous plasmas that include the aforementioned inverse bremsstrahlung effect. Comparisons with kinetic particle-in-cell simulations are satisfactory, provided the acoustic wave-breaking limit and the self-focusing regime are not reached. An application of the Cross-Beam Energy Transfer model is shown for a typical case of indirect-drive implosion in a gold hohlraum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Method and apparatus for simultaneous spectroelectrochemical analysis

An apparatus and method of simultaneous spectroelectrochemical analysis is disclosed. A transparent surface is provided. An analyte solution on the transparent surface is contacted with a working electrode and at least one other electrode. Light from a light source is focused on either a surface of the working electrode or the analyte solution. The light reflected from either the surface of the working electrode or the analyte solution is detected. The potential of the working electrode is adjusted, and spectroscopic changes of the analyte solution that occur with changes in thermodynamic potentials are monitored.

47 OTHER INSTRUMENTATION↗
DOE OSTI2014

Drifting solutions with elliptic symmetry for the compressible Navier-Stokes equations with density-dependent viscosity

In this paper, we investigate the analytical solutions of the compressible Navier-Stokes equations with dependent-density viscosity. By using the characteristic method, we successfully obtain a class of drifting solutions with elliptic symmetry for the Navier-Stokes model wherein the velocity components are governed by a generalized Emden dynamical system. In particular, when the viscosity variables are taken the same as Yuen [M. W. Yuen, “Analytical solutions to the Navier-Stokes equations,” J. Math. Phys. 49, 113102 (2008)], our solutions constitute a generalization of that obtained by Yuen. Interestingly, numerical simulations show that the analytical solutions can be used to explain the drifting phenomena of the propagation wave like Tsunamis in oceans.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Design Guidelines for Predicting Stress in Cemented Doublets Undergoing Temperature Change

This work explores quick predictive methods for calculating potentially risky stresses in cemented doublets underdoing temperature change that agree well with finite element analysis. We also provide guidelines for avoiding stress concentrations. Cemented doublets are ubiquitous in the optics world, but it is not always clear if they will perform as intended in a specific thermal environment due to excessive stress. The stress in the optic can cause debonding or even fracture. Yoder addressed this problem using analytical results from Chen and Nelson but also expressed some skepticism as to the accuracy and usefulness of those results. Closed-form analytical solutions are extremely useful when doing trade-off studies to see if a particular geometry or glass choice will cause stress problems. The problem of bonded materials undergoing temperature change is of great interest in the electronics industry, and improved analytical solutions are now available that address the issues raised by Yoder, and they have been validated using finite element analysis (FEA). But the newer solutions still have limitations. The analytical solutions presented here cannot predict stresses due to stress concentrations or singularities at the edge of the bondline—this still has to be done using FEA—although some advanced methods are available. Also, the existing analytical methods assume a flat interface between adjoining elements; they do not account for the effect of curvature. We find that most doublet failures start at the edges and involve debonding of the glass. Glass fracture is rare for doublets, accounting for less than 10% of the failures. This paper provides the reader with a general introduction to a large and dense body of literature in order to demystify the subject of temperature-induced glass bonding failures.

36 MATERIALS SCIENCE↗
DOE OSTI2015

A Finite Element Model for Simulation of Carbon Dioxide Sequestration

We present a hydro-mechanical model, followed by stress, deformation, and shear-slip failure analysis for geological sequestration of carbon dioxide (CO2). The model considers the poroelastic effects by taking into account of the two-way coupling between the geomechanical response and the fluid flow process. Analytical solutions for pressure and deformation fields were derived for a typical geological sequestration scenario in our previous work. A finite element approach is introduced here for numerically solving the hydro-mechanical model with arbitrary boundary conditions. The numerical approach was built on an open-source finite element code Elmer, and results were compared to the analytical solutions. The shear-slip failure analysis was presented based on the numerical results, where the potential failure zone is identified. Information is relevant to the prediction of the maximum sustainable injection rate or pressure. The effects of caprock permeability on the fluid pressure, deformation, stress, and the shear-slip failure zone were also quantitatively studied. It was shown that a larger permeability in caprock and base rock leads to a larger uplift but a smaller shear-slip failure zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A mass-transfer particle-tracking method for simulating transport with discontinuous diffusion coefficients

The problem of a spatially discontinuous diffusion coefficient (D(x)) is one that may be encountered in hydrogeologic systems due to natural geological features or as a consequence of numerical discretization of flow properties. To date, mass-transfer particle-tracking (MTPT) methods, a family of Lagrangian methods in which diffusion is jointly simulated by random walk and diffusive mass transfers, have been unable to solve this problem. This manuscript presents a new mass-transfer (MT) algorithm that enables MTPT methods to accurately solve the problem of discontinuous D(x). To achieve this, we derive a semi-analytical solution to the discontinuous D(x) problem by employing a predictor-corrector approach, and we use this semi-analytical solution as the weighting function in a reformulated MT algorithm. Furthermore, this semi-analytical solution is generalized for cases with multiple 1D interfaces as well as for 2D cases, including a 2 × 2 tiling of 4 subdomains that corresponds to a numerically-generated diffusion field. The solutions generated by this new mass-transfer algorithm closely agree with an analytical 1D solution or, in more complicated cases, trusted numerical results, demonstrating the success of our proposed approach.

42 ENGINEERING↗
DOE OSTI2018

Semi-numerical solution to a fractal telegraphic dual-porosity fluid flow model

We present a semi-numerical solution to a fractal telegraphic dual-porosity fluid flow model. It combines Laplace transform and finite difference schemes. The Laplace transform handles the time variable whereas the finite difference method deals with the spatial coordinate. This semi-numerical scheme is not restricted by space discretization and allows the computation of a solution at any time without compromising numerical stability or the mass conservation principle. Our formulation results in a non-analytically-solvable second-order differential equation whose spatial discretization yields a tridiagonal linear algebraic system. Moreover, we describe comparisons between semi-numerical and semi-analytical solutions for particular cases. Results agree well with those from semi-analytic solutions. Furthermore, we expose a parametric analysis from the coupled model in order to show the effects of relevant parameters on pressure profiles and flow rates for the case where neither analytic nor semi-analytic solutions are available.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Review of analytical studies on seismic fluid-structure interaction of base-supported cylindrical tanks

Analytical solutions for seismic fluid-structure interaction (FSI) of tanks are used by analysts and engineers for a number of purposes, including preliminary sizing, risk assessment, and verification of numerical models. Analytical solutions for seismic FSI of tanks were first developed in the 1930s and have been extended and/or modified in the following nine decades to accommodate different seismic inputs and boundary conditions. This report reviews published analytical solutions for rigid and flexible, base-supported, cylindrical tanks subjected to unidirectional horizontal motion of a small amplitude. The studies parsed the fluid-structure (FSI) response into an impulsive and a convective component, and addressed frequencies of lateral motions of the tank (i.e., impulsive frequencies), hydrodynamic pressures, frequencies of waves (i.e., convective frequencies), wave heights, and reactions at the base. The solutions were categorized into exact and approximate solutions based on the methodologies used. Herein, FSI responses are calculated and compared for a range of tank dimensions using the analytical solutions from the different studies. The responses presented here are normalized to be unitless, and can be used for tanks with different dimensions and mechanical properties, and subjected to different input motions. A sample steel water tank is analyzed using the normalized solutions for two earthquake ground motions with different frequency contents. Ignoring tank flexibility may significantly underestimate seismic responses of the tank and its contained fluid. The approximate solutions, which have been widely used for decades, may not be sufficiently accurate in the modern era. Calculation errors in the reviewed studies are identified and corrected, and differences are quantified for the sample tank.

42 ENGINEERING↗
DOE OSTI2015

Fundamental solutions in a half space of two-dimensional hexagonal quasicrystal and their applications

Fundamental phonon-phason field in a half-infinite space of two-dimensional hexagonal quasicrystal is derived, on the basis of general solutions in terms of quasi-harmonic functions, by virtue of the trial-and-error technique. Extended Boussinesq and Cerruti problems are studied. Appropriate potential functions are assumed and corresponding fundamental solutions are explicitly derived in terms of elementary functions. The boundary integral equations governing the contact and crack problems are constructed from the present fundament solutions. The obtained analytical solutions can serve as guidelines for future indentation tests via scanning probe microscopy and atomic force microscopy methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Progress and investigation on lattice Boltzmann modeling of multiple immiscible fluids or components with variable density and viscosity ratios

Lattice Boltzmann models for simulating multiphase flows are relatively new, and much work remains to be done to demonstrate their ability to solve fundamental test cases before they are considered for engineering problems. From this perspective, a hydrodynamic lattice Boltzmann model for simulating immiscible multiphase flows with high density and high viscosity ratios, up to O(1000) and O(100) respectively, is presented and validated against analytical solutions. The method is based on a two phase flow model with operators extended to handle N immiscible fluids. The current approach is O(N) in computational complexity for the number of different gradient approximations. This is a major improvement, considering the O(N{sup 2}) complexity found in most works. A sequence of systematic and essential tests have been conducted to establish milestones that need to be met by the proposed approach (as well as by other methods). First, the method is validated qualitatively by demonstrating its ability to address the spinodal decomposition of immiscible fluids. Second, the model is quantitatively verified for the case of multilayered planar interfaces. Third, the multiphase Laplace law is studied for the case of three fluids. Fourth, a quality index is developed for the three-phase Laplace–Young’s law, which concerns the position of the interfaces between the fluids resulting from the different surface tensions. The current model is compatible with the analytical solution, and is shown to be first order accurate in terms of this quality index. Finally, the multilayered Couette’s flow is studied. In this study, numerical results can recover the analytical solutions for all the selected test cases, as long as unit density ratios are considered. For high density and high viscosity ratios, the analytical solution is recovered for all tests, except that of the multilayered Couette’s flow. Numerical results and a discussion are presented for this unsuccessful test case. It is believed that other LB models may have the same problem in addressing the simulation of multiphase flows with variable density ratios.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1999

A Spherical Earth Solution for TOA Lightning Location Retrieval

With the advent of high-speed digital computers, the method of chi square minimization is a highly practical means for analyzing a wide variety of (otherwise intractable) nonlinear inversion problems in applied mathematical physics. Little thought or effort is required to apply the chi square method to obtain quick and reasonable estimates of a solution, and the method offers a means to assess retrieval errors. Because the method is simple and practical it is sometimes hastily applied to problems that can be solved by formal analytic or quasi-analytic means. Presently, Global Atmospherics Inc. (GAI) finds the minimum of a chi square function to analyze time-of-arrival (TOA) and magnetic bearing data derived from the National Lightning Detection Network (NLDN); ellipsoidal Earth geometry is assumed. An analytic solution to this problem has not yet been attained, but the consideration and solving of less general problem statements might eventually lead to a final solution. In the present study, the problem of retrieving lightning ground-strike location on a spherical Earth surface using a network of 4 or more time-of-arrival (TOA) sensors is considered. It is shown that this problem has an analytic solution and therefore does not require the use of nonlinear estimation theory (such as the chi square method mentioned above). The mathematical robustness of the analytic solution is tested using computer-generated lightning sources and simulated TOA measurement errors. A quasi-analytic extension of the spherical Earth solution for an oblate spheroidal Earth geometry is considered in a related study. The incorporation of magnetic bearing information into these analytic solutions would lead to a general and elegant analytic retrieval scheme that would most likely replace the chi square estimation theory currently employed by Global Atmospherics Inc. (GAI).

Koshak, William J,↗
DOE OSTI2014

SU-E-T-517: Analytic Formalism to Compute in Real Time Dose Distributions Delivered by HDR Units

Purpose: Develop an analytical algorithm to compute the dose delivered by Ir-192 dwell positions with high accuracy using the 3-dimensional (3D) dose distribution of an HDR source. Using our analytical function, the dose delivered by an HDR unit as treatment progresses can be determined using the actual delivered temporal and positional data of each individual dwell. Consequently, true delivered dose can be computed when each catheter becomes active. We hypothesize that the knowledge of such analytical formulation will allow developing HDR systems with a real time treatment evaluation tool to avoid mistreatments. Methods: In our analytic formulation, the dose is computed by using the full anisotropic function data of the TG 43 formalism with 3D ellipsoidal function. The discrepancy between the planned dose and the delivered dose is computed using an analytic perturbation method over the initial dose distribution. This methodology speeds up the computation because only changes in dose discrepancies originated by spatial and temporal deviations are computed. A dose difference map at the point of interest is obtained from these functions and this difference can be shown during treatment in real time to examine the treatment accuracy. Results: We determine the analytical solution and a perturbation function for the 3 translational 3 rotational, and 1D temporal errors in source distributions. The analytic formulation is a sequence of simple equations that can be processed in any modern computer in few seconds. Because computations are based in an analytical solution, small deviations of the dose when sub-millimeter positional changes occur can be detected. Conclusions: We formulated an analytical method to compute 4D dose distributions and dose differences based on an analytical solution and perturbations to the original dose. This method is highly accurate and can be.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Lie group analysis of the point-reactor neutron kinetics equations for various reactivity models

This work applies Lie Group Theory to the nonlinear point-reactor neutron kinetics equations for several reactivity feedback and reactivity insertion models. The goal of this work is to connect known analytical solutions to their corresponding Lie groups as well as to obtain new analytical solutions and the circumstances from which they arise. The reactivity models we consider are the Nordheim-Fuchs Model, an arbitrary-order polynomial reactivity insertion model, and the Fuchs Ramp-Insertion Model. In all three models, we generalize the analysis when possible and obtain analytical solutions in every case. For the first two models, we succeed in finding the Lie groups corresponding to classical solutions and find additional solutions that were otherwise unknown. Finally, for the final model, we obtain a reduced-order nonlinear kinetics equation which is not analytically solvable, but by applying a second-order truncation, we find a novel closed-form solution and demonstrate its validity with a numerical example.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

On the first crossing distributions in fractional Brownian motion and the mass function of dark matter haloes

We construct an integral equation for the first crossing distributions for fractional Brownian motion in the case of a constant barrier and we present an exact analytical solution. Additionally we present first crossing distributions derived by simulating paths from fractional Brownian motion. We compare the results of the analytical solutions with both those of simulations and those of some approximated solutions which have been used in the literature. Finally, we present multiplicity functions for dark matter structures resulting from our analytical approach and we compare with those resulting from N-body simulations. We show that the results of analytical solutions are in good agreement with those of path simulations but differ significantly from those derived from approximated solutions. Additionally, multiplicity functions derived from fractional Brownian motion are poor fits of the those which result from N-body simulations. We also present comparisons with other models which are exist in the literature and we discuss different ways of improving the agreement between analytical results and N-body simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Analyzing the Validity of Brazilian Testing Using Digital Image Correlation and Numerical Simulation Techniques

Characterizing the mechanical behavior of rocks plays a crucial role to optimize the fracturing process in unconventional reservoirs. However, due to the intrinsic anisotropy and heterogeneity in unconventional resources, fracture process prediction remains the most significant challenge for sustainable and economic hydrocarbon production. During the deformation tracking under compression, deploying conventional methods (strain gauge, extensometer, etc.) is insufficient to measure the deformation since the physical attachment of the device is restricted to the size of the sample, monitoring limited point-wise deformation, producing difficulties in data retrieval, and a tendency to lose track in failure points, etc. Where conventional methods are limited, the application of digital image correlation (DIC) provides detailed and additional information of strain evolution and fracture patterns under loading. DIC is an image-based optical method that records an object with a camera and monitors the random contrast speckle pattern painted on the facing surface of the specimen. To overcome the existing limitations, this paper presents numerical modeling of Brazilian disc tests under quasi-static conditions to understand the full-field deformation behaviors and finally, it is validated by DIC. As the direct tensile test has limitations in sample preparation and test execution, the Brazilian testing principle is commonly used to evaluate indirectly the tensile strength of rocks. The two-dimensional numerical model was built to predict the stress distribution and full-field deformation on Brazilian disc under compression based on the assumptions of a homogenous, isotropic and linear elastic material. The uniaxial compression test was conducted using the DIC technique to determine the elastic properties of Spider Berea sandstone, which were used as inputs for the simulation model. The model was verified by the analytical solution and compared with the digital image correlation. The numerical simulation results showed that the solutions matched reasonably with the analytical solutions where the maximum deviation of stress distribution was obtained as 14.59%. The strain evolution (normal and shear strains) and displacements along the central horizontal and vertical planes were investigated in three distinguishable percentages of peak loads (20%, 40%, and 90%) to understand the deformation behaviors in rock. The simulation results demonstrated that the strain evolution contours consistently matched with DIC generated contours with a reasonable agreement. The changes in displacement along the central horizontal and vertical planes showed that numerical simulation and DIC generated experimental results were repeatable and matched closely. In terms of validation, Brazilian testing to measure the indirect tensile strength of rocks is still an issue of debate. The numerical model of fracture propagation supported by digital image correlation from this study can be used to explain the fracturing process in the homogeneous material and can be extended to non-homogeneous cases by incorporating heterogeneity, which is essential for rock mechanics field applications.

58 GEOSCIENCES↗
NASA NTRS1993

Simulation modeling for long duration spacecraft control systems

The use of simulation is described and it is contrasted to analytical solution techniques for evaluation of analytical reliability models. The role importance sampling plays in simulation of models of this type was also discussed. The simulator tool used for our analysis is described. Finally, the use of the simulator tool was demonstrated by applying it to evaluate the reliability of a fault tolerant hypercube multiprocessor intended for spacecraft designed for long duration missions. The reliability analysis was used to highlight the advantages and disadvantages offered by simulation over analytical solution of Markovian and non-Markovian reliability models.

Boyd, Mark A.↗
DOE OSTI2014

Sci—Fri PM: Topics — 05: Experience with linac simulation software in a teaching environment

Medical linear accelerator education is usually restricted to use of academic textbooks and supervised access to accelerators. To facilitate the learning process, simulation software was developed to reproduce the effect of medical linear accelerator beam adjustments on resulting clinical photon beams. The purpose of this report is to briefly describe the method of operation of the software as well as the initial experience with it in a teaching environment. To first and higher orders, all components of medical linear accelerators can be described by analytical solutions. When appropriate calibrations are applied, these analytical solutions can accurately simulate the performance of all linear accelerator sub-components. Grouped together, an overall medical linear accelerator model can be constructed. Fifteen expressions in total were coded using MATLAB v 7.14. The program was called SIMAC. The SIMAC program was used in an accelerator technology course offered at our institution; 14 delegates attended the course. The professional breakdown of the participants was: 5 physics residents, 3 accelerator technologists, 4 regulators and 1 physics associate. The course consisted of didactic lectures supported by labs using SIMAC. At the conclusion of the course, eight of thirteen delegates were able to successfully perform advanced beam adjustments after two days of theory and use of the linac simulator program. We suggest that this demonstrates good proficiency in understanding of the accelerator physics, which we hope will translate to a better ability to understand real world beam adjustments on a functioning medical linear accelerator.

07 ISOTOPE AND RADIATION SOURCES↗