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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Development and Verification of a Numerical Library for Solving Global Terrestrial Multiphysics Problems

Abstract Current generation Land Surface Models (LSMs) routinely simulate many nonlinear multiphysics processes. The complexity of future generation LSMs is expected to increase as critical new biophysical and biogeochemical processes are incorporated. Current generation LSMs have several shortcomings including the lack of robust numerical methods to solve discretized equations, monolithic software design that hinders testing of a process representation in isolation from other components, and absence of a flexible coupling framework to solve tightly coupled multiphysics problems. While the LSMs community vigorously evaluates the accuracy of the conceptual model to represent reality (validation), the accuracy of the numerical implementation of the conceptual model (verification) is seldom examined. Method of Manufactured Solutions is a technique for verifying complex codes when analytical solutions are unavailable. In this work, we present a stand‐alone, open source numerical library for solving global terrestrial multiphysics processes that includes a flexible coupling framework. Robust numerical solution for linear and nonlinear equations is provided via the Portable, Extensible Toolkit for Scientific Computation library. We verify the numerical library using Method of Manufactured Solutions for a range of problems comprising single and multiphysics steady state problems that are applied in one or multiple physical domains. This work provides an example of incorporating code verification as an integral part of Land Surface Model development activities.

58 GEOSCIENCES↗
DOE OSTI2015

Accurate modeling and inversion of electrical resistivity data in the presence of metallic infrastructure with known location and dimension

Electrical resistivity tomography (ERT) has been widely used in environmental applications to study processes associated with subsurface contaminants and contaminant remediation. Anthropogenic alterations in subsurface electrical conductivity associated with contamination often originate from highly industrialized areas with significant amounts of buried metallic infrastructure. The deleterious influence of such infrastructure on imaging results generally limits the utility of ERT where it might otherwise prove useful for subsurface investigation and monitoring. In this manuscript we present a method of accurately modeling the effects of buried conductive infrastructure within the forward modeling algorithm, thereby removing them from the inversion results. The method is implemented in parallel using immersed interface boundary conditions, whereby the global solution is reconstructed from a series of well-conditioned partial solutions. Forward modeling accuracy is demonstrated by comparison with analytic solutions. Synthetic imaging examples are used to investigate imaging capabilities within a subsurface containing electrically conductive buried tanks, transfer piping, and well casing, using both well casings and vertical electrode arrays as current sources and potential measurement electrodes. Results show that, although accurate infrastructure modeling removes the dominating influence of buried metallic features, the presence of metallic infrastructure degrades imaging resolution compared to standard ERT imaging. However, accurate imaging results may be obtained if electrodes are appropriately located.

imaging↗
NASA NTRS1974

Energy transfer in volume-reflecting heat shields

An approximate analysis of radiative transfer in highly scattering materials was developed based on the Kubelka-Munk differential equations--a set of two differential equations representing the spatial rate of change of radiative half-fluxes within the scattering media. These approximate solutions of the Kubelka-Munk equations together with analytic solutions for the steady state temperature distribution for two types of boundary conditions are given. These solutions show the influence of back surface reflectance, scattering power, incident radiative flux parameter and boundary conductive flux parameter on overall reflectance and temperature distributions. This radiation field analysis, adapted to spherical geometry, was applied to the evaluation of the thermal performance of teflon and fritted quartz as heat protection materials for entry into the atmosphere of Jupiter.

Weston, K. C.↗
NASA NTRS1977

Analytic theory of orbit contraction

The motion of a satellite in orbit, subject to atmospheric force and the motion of a reentry vehicle are governed by gravitational and aerodynamic forces. This suggests the derivation of a uniform set of equations applicable to both cases. For the case of satellite motion, by a proper transformation and by the method of averaging, a technique appropriate for long duration flight, the classical nonlinear differential equation describing the contraction of the major axis is derived. A rigorous analytic solution is used to integrate this equation with a high degree of accuracy, using Poincare's method of small parameters and Lagrange's expansion to explicitly express the major axis as a function of the eccentricity. The solution is uniformly valid for moderate and small eccentricities. For highly eccentric orbits, the asymptotic equation is derived directly from the general equation. Numerical solutions were generated to display the accuracy of the analytic theory.

Vinh, N. X.↗
NASA NTRS1990

A defect stream function formulation for compressible turbulent boundary layers

Progress to date on the development of a method for turbulent, wall-bounded flow which uses the defect stream function formulation in the outer layer and an analytic law of the wall and wake formulation in the inner region is reviewed. This two-formulation approach avoids the need to computationally resolve the high-gradient inner layer. One of the most appealing recent developments is the transformation of the compressible governing equation for the defect stream function into a linear, second-order differential equation which has analytic solutions for many problems of practical interest. Numerical and analytic results for incompressible and compressible flows are shown to be in excellent agreement with experimental results. In this paper the two-formulation approach is applied to primitive-variable computations. Excellent comparisons with experiment are presented for two compressible flat plate flows.

Barnwell, Richard W.↗
NASA NTRS1994

Electron transport and energy degradation in the ionosphere: Evaluation of the numerical solution, comparison with laboratory experiments and auroral observations

Auroral electron transport calculations are a critical part of auroral models. We evaluate a numerical solution to the transport and energy degradation problem. The numerical solution is verified by reproducing simplified problems to which analytic solutions exist, internal self-consistency tests, comparison with laboratory experiments of electron beams penetrating a collision chamber, and by comparison with auroral observations, particularly the emission ratio of the N2 second positive to N2(+) first negative emissions. Our numerical solutions agree with range measurements in collision chambers. The calculated N(2)2P to N2(+)1N emission ratio is independent of the spectral characteristics of the incident electrons, and agrees with the value observed in aurora. Using different sets of energy loss cross sections and different functions to describe the energy distribution of secondary electrons that emerge from ionization collisions, we discuss the uncertainties of the solutions to the electron transport equation resulting from the uncertainties of these input parameters.

Lummerzheim, D.↗
NASA NTRS1999

Slow Crack Growth Analysis of Brittle Materials with Finite Thickness Subjected to Constant Stress-Rate Flexural Loading

A two-dimensional, numerical analysis of slow crack growth (SCG) was performed for brittle materials with finite thickness subjected to constant stress-rate ("dynamic fatigue") loading in flexure. The numerical solution showed that the conventional, simple, one-dimensional analytical solution can be used with a maximum error of about 5% in determining the SCG parameters of a brittle material with the conditions of a normalized thickness (a ratio of specimen thickness to initial crack size) T > 3.3 and of a SCG parameter n > 10. The change in crack shape from semicircular to elliptical configurations was significant particularly at both low stress rate and low T, attributed to predominant difference in stress intensity factor along the crack front. The numerical solution of SCG parameters was supported within the experimental range by the data obtained from constant stress-rate flexural testing for soda-lime glass microslides at ambient temperature.

Chio, S. R.↗
NASA NTRS2019

Observations on the Computation of Eigenvalue and Eigenvector Jacobians

Many scientific and engineering problems benefit from analytic expressions for eigenvalue and eigenvector derivatives with respect to the elements of the parent matrix. While there exists extensive literature on the calculation of these derivatives, which take the form of Jacobian matrices, there are a variety of deficiencies that have yet to be addressed — including the need for both left and right eigenvectors, limitations on the matrix structure, and issues with complex eigenvalues and eigenvectors. This work addresses these deficiencies by proposing a new analytic solution for the eigenvalue and eigenvector derivatives. The resulting analytic Jacobian matrices are numerically efficient to compute and are valid for the general complex case. It is further shown that this new general result collapses to previously known relations for the special cases of real symmetric matrices and real diagonal matrices. Finally, the new Jacobian expressions are validated using forward finite differencing and performance is compared with another technique.

Jacobian↗
NASA NTRS1978

The general solution to the classical problem of the finite Euler-Bernoulli beam

An analytical solution is obtained for the problem of free and forced vibrations of a finite Euler-Bernoulli beam with arbitrary (partially fixed) boundary conditions. The effects of linear viscous damping, Winkler foundation, constant axial tension, a concentrated mass, and an arbitrary forcing function are included in the analysis. No restriction is placed on the values of the parameters involved, and the solution presented here contains all cited previous solutions as special cases.

Hussaini, M. Y.↗
DOE OSTI2014

Binding energy and mechanical stability of single- and multi-walled carbon nanotube serpentines

Recently, Geblinger et al. [Nat. Nanotechnol. 3, 195 (2008)] and Machado et al. [Phys. Rev. Lett. 110, 105502 (2013)] reported the experimental and molecular dynamics realization of S-like shaped single-walled carbon nanotubes (CNTs), the so-called CNT serpentines. We reported here results from continuum modeling of the binding energy γ between different single- and multi-walled CNT serpentines and substrates as well as the mechanical stability of the CNT serpentine formation. The critical length for the mechanical stability and adhesion of different CNT serpentines are determined in dependence of E{sub i}I{sub i}, d, and γ, where E{sub i}I{sub i} and d are the CNT bending stiffness and distance of the CNT translation period. Our continuum model is validated by comparing its solution to full-atom molecular dynamics calculations. The derived analytical solutions are of great importance for understanding the interaction mechanism between different single- and multi-walled CNT serpentines and substrates.

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

Modeling elasto-viscoplasticity in a consistent phase field framework

Existing continuum level phase field plasticity theories seek to solve plastic strain by minimizing the shear strain energy. However, rigorously speaking, for thermodynamic consistency it is required to minimize the total strain energy unless there is proof that hydrostatic strain energy is independent of plastic strain which is unfortunately absent. In this work, we extend the phase-field microelasticity theory of Khachaturyan et al. by minimizing the total elastic energy with constraint of incompressibility of plastic strain. We show that the flow rules derived from the Ginzburg-Landau type kinetic equation can be in line with Odqvist's law for viscoplasticity and Prandtl-Reuss theory. Free surfaces (external surfaces or internal cracks/voids) are treated in the model. Deformation caused by a misfitting spherical precipitate in an elasto-plastic matrix is studied by large-scale three-dimensional simulations in four different regimes in terms of the matrix: (a) elasto-perfectly-plastic, (b) elastoplastic with linear hardening, (c) elastoplastic with power-law hardening, and (d) elasto-perfectly-plastic with a free surface. The results are compared with analytical/numerical solutions of Lee et al. for (a-c) and analytical solution derived in this work for (d). Additionally, the J integral of a fixed crack is calculated in the phase-field model and discussed in the context of fracture mechanics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Quantifying local and global mass balance errors in physics-informed neural networks

Physics-informed neural networks (PINN) have recently become attractive for solving partial differential equations (PDEs) that describe physics laws. By including PDE-based loss functions, physics laws such as mass balance are enforced softly in PINN. This paper investigates how mass balance constraints are satisfied when PINN is used to solve the resulting PDEs. We investigate PINN’s ability to solve the 1D saturated groundwater flow equations (diffusion equations) for homogeneous and heterogeneous media and evaluate the local and global mass balance errors. We compare the obtained PINN’s solution and associated mass balance errors against a two-point finite volume numerical method and the corresponding analytical solution. We also evaluate the accuracy of PINN in solving the 1D saturated groundwater flow equation with and without incorporating hydraulic heads as training data. We demonstrate that PINN’s local and global mass balance errors are significant compared to the finite volume approach. Tuning the PINN’s hyperparameters, such as the number of collocation points, training data, hidden layers, nodes, epochs, and learning rate, did not improve the solution accuracy or the mass balance errors compared to the finite volume solution. Mass balance errors could considerably challenge the utility of PINN in applications where ensuring compliance with physical and mathematical properties is crucial.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Complete spectrum of the infinite-U Hubbard ring using group theory

We present a full analytical solution of the multiconfigurational strongly correlated mixed-valence problem corresponding to the N-Hubbard ring filled with N−1 electrons, and infinite on-site repulsion. While the eigenvalues and the eigenstates of the model are known already, analytical determination of their degeneracy is presented here for the first time. The full solution, including degeneracy count, is achieved for each spin configuration by mapping the Hubbard model into a set of Hückel-annulene problems for rings of variable size. The number and size of these effective Hückel annulenes, both crucial to obtain Hubbard states and their degeneracy, are determined by solving a well-known combinatorial enumeration problem, the necklace problem for N−1 beads and two colors, within each subgroup of the C{sub N−1} permutation group. Symmetry-adapted solution of the necklace enumeration problem is finally achieved by means of the subduction of coset representation technique [S. Fujita, Theor. Chim. Acta 76, 247 (1989)], which provides a general and elegant strategy to solve the one-hole infinite-U Hubbard problem, including degeneracy count, for any ring size. The proposed group theoretical strategy to solve the infinite-U Hubbard problem for N−1 electrons is easily generalized to the case of arbitrary electron count L, by analyzing the permutation group C{sub L} and all its subgroups.

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

An asymptotically compatible meshfree quadrature rule for nonlocal problems with applications to peridynamics

In this paper, we present a meshfree quadrature rule for compactly supported nonlocal integro-differential equations (IDEs) with radial kernels. We apply this rule to develop a meshfree discretization of a peridynamic solid mechanics model that requires no background mesh. Existing discretizations of peridynamic models have been shown to exhibit a lack of asymptotic compatibility to the corresponding linearly elastic local solution. By posing the quadrature rule as an equality constrained least squares problem, we obtain asymptotically compatible convergence by introducing polynomial reproduction constraints. Our approach naturally handles traction-free conditions, surface effects, and damage modeling for both static and dynamic problems. We demonstrate high-order convergence to the local theory by comparing to manufactured solutions and to cases with crack singularities for which an analytic solution is available. We verify the applicability of the approach to realistic problems by reproducing high-velocity impact results from the Kalthoff–Winkler experiments.

42 ENGINEERING↗
DOE OSTI2023

A fast Fourier transform-based solver for elastic micropolar composites

This work presents a spectral micromechanical formulation for obtaining the full-field and homogenized response of elastic micropolar composites. The algorithm relies on a coupled set of convolution integral equations for the micropolar strains, where periodic Green’s operators associated with a linear homogeneous reference medium are convolved with functions of the Cauchy and couple stress fields that encode the material’s heterogeneity, as well as any potential material nonlinearity. Such convolution integral equations take an algebraic form in the reciprocal Fourier space that can be solved iteratively. In this vein, the fast Fourier transform (FFT) algorithm is leveraged to accelerate the numerical solution, resulting in a mesh-free formulation in which the periodic unit cell representing the heterogeneous material can be discretized by a regular grid of pixels in two dimensions (or voxels in three dimensions). For verification, the numerical solutions obtained with the micropolar FFT solver are compared with analytical solutions for a matrix with a dilute circular inclusion subjected to plane strain loading. The developed computational framework is then used to study length-scale effects and effective (micropolar) moduli of composites with various topological configurations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2003

The Effect of Rotating a Faraday Disc Perpendicular to an Applied Magnetic Field Theory and Experiment

A magnetohydrodynamic model that examines the effect of rotating an electrically conducting cylinder with a uniform external magnetic field applied orthogonal to its axis is presented. Noting a simple geometry, it can be classified as a fundamental dynamo problem. For the case of an infinitely long cylinder, an analytical solution is obtained and analyzed in detail. A semi-analytical model was developed that considers a finite cylinder. Experimental data from a spinning brass wheel in the presence of Earth's magnetic field were compared to the proposed theory and found to fit well.

Mazuruk, Konstantin↗