SEARCH · Search NASA
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.
Surface coordinates of face-milled spiral-bevel gear teeth computed by method involving numerical solution of governing equations. Needed to generate mathematical models of tooth surfaces for use in finite-element analyses of stresses, strains, and vibrations in meshing spiral-bevel gears.
This paper briefly describes the nanometer-class metrology system used in this testbed to estimate the length and orientation of the science baseline vector, which cannot be measured directly. The focus is on the mathematical inversion problem that results and its solution.
Fractional calculus is a relatively old yet emerging field of mathematics with the widest range of engineering and biomedical applications. Despite being an incredibly powerful tool, it, however, requires promotion in the engineering community. Rheology is undoubtedly one of the fields where fractional calculus has become an integral part of cutting-edge research. There exists extensive literature on the theoretical, experimental, and numerical treatment of various fractional viscoelastic flows in constraint geometries. However, the general theoretical approach that unites several most commonly used models is missing. Here we present exact analytical solutions for fractional viscoelastic flow in a circular pipe. We find velocity profiles and shear stresses for fractional Maxwell, Kelvin–Voigt, Zener, Poynting–Thomson, and Burgers models. The dynamics of these quantities are studied with respect to normalized pipe radius, fractional orders, and elastic moduli ratio. Three different types of behavior are identified: monotonic increase, resonant, and aperiodic oscillations. The models developed are applicable in the widest material range and allow for the alteration of the balance between viscous and elastic properties of the materials.
We analyze the well-posedness of an anisotropic, nonlocal diffusion equation. Establishing an equivalence between weighted and unweighted anisotropic nonlocal diffusion operators in the vein of unified nonlocal vector calculus, we apply our analysis to a class of fractional-order operators and present rigorous estimates for the solution of the corresponding anisotropic anomalous diffusion equation. Furthermore, we extend our analysis to the anisotropic diffusion-advection equation and prove well-posedness for fractional orders s ∊ [0.5, 1). We also present an application of the advection-diffusion equation to anomalous transport of solutes.
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.
Computational aeroelastic analyses typically use a mathematical model for the structural modes of a flexible structure and a nonlinear aerodynamic model that can generate a plurality of unsteady aerodynamic responses based on the structural modes for conditions defining an aerodynamic condition of the flexible structure. In the present invention, a linear state-space model is generated using a single execution of the nonlinear aerodynamic model for all of the structural modes where a family of orthogonal functions is used as the inputs. Then, static and dynamic aeroelastic solutions are generated using computational interaction between the mathematical model and the linear state-space model for a plurality of periodic points in time.
By using the theory of polynomial matrix pencils, we construct the asymptotics of linearly independent solutions of homogeneous singularly perturbed system of linear differential equations of any order m with matrices at higher derivatives degenerating as a small parameter approaches zero. The general case is analyzed. In this case, the limit matrix pencil has several finite and infinite elementary divisors of both the same and different multiplicities. The corresponding asymptotic estimates are presented.
An approach to the solution of matrix problems resulting from integral equations of mathematical physics is presented. Based on the inherent smoothness in such equations, the problem is reformulated using a set of orthogonal basis vectors, leading to an equivalent coefficient problem which can be of lower order without significantly impairing the accuracy of the solution. This approach was evaluated using a two-dimensional Neumann problem describing the inviscid, incompressible flow over an airfoil. Two different kinds of mode functions were investigated, namely eigenfunction series and Fourier series. The method using Fourier series was found preferable. It uses all of the coefficients from a Fast Fourier Transform algorithm in an approximate method which exploits the known structure of the transformed coefficient matrix and very promising results for the flow over a realistic airfoil are obtained. On the basis of the results presented here, an order of magnitude reduction in this computer time can be expected for such problems as compared with the time for a direct matrix solution.
Verification of multi-physics simulation software against problems with known analytic or semi-analytic solutions is an important aspect of research into a wide variety of fields involving the motion of fluids, shock physics and other dynamic material properties. Previous work comparing simulation results against analytic solutions has been ad-hoc, with developers frequently writing their own analytic solvers. This has resulted in a large amount of duplicated effort. The python library ExactPack has been developed as a collection of analytic and semi-analytic solvers to a variety of multi-physics problems, providing a consistent API to a set of well-tested solver implementations.
The use of a point source of an impactor energy and momentum to replace the effects of the impactor is examined. The general framework and notation of the impact cratering problems are described; it is determined that the cratering phenomena are governed by Froude, Cauchy, and Reynolds numbers. The coupling parameter concept is defined mathematically as the measure that governs limit point source solutions. Examples of cases where coupling parameters are used are presented. The relationships of the coupling parameter concept with steady flow and the Z-model of cratering of Maxwell (1973, 1977) are studied. Crater size, ejecta distributions, growth histories, time of formation, melt volume, and shock decay for various scale factors for impact cratering mechanics are calculated, and the applicability of the coupling parameter to the study of cratering mechanics is revealed.
In this study, a time-periodic MGRIT algorithm is proposed as a means to reduce the time-to-solution of numerical algorithms by exploiting the time periodicity inherent to many applications in science and engineering. The time-periodic MGRIT algorithm is applied to a variety of linear and nonlinear single- and multiphysics problems that are periodic-in-time. It is demonstrated that the proposed parallel-in-time algorithm can obtain the same time-periodic steady-state solution as sequential time-stepping. It is shown that the required number of MGRIT iterations can be estimated a priori and that the new MGRIT variant can significantly and consistently reduce the time-to-solution compared to sequential time-stepping, irrespective of the number of dimensions, linear or nonlinear PDE models, single-physics or coupled problems and the employed computing resources. The numerical experiments demonstrate that the time-periodic MGRIT algorithm enables a greater level of parallelism yielding faster turnaround, and thus, facilitating more complex and more realistic problems to be solved.
Abstract We show that local weak solutions to anisotropic elliptic equations with bounded and measurable coefficients, whose prototype is are continuous. Our analysis exploits a particular intrinsic geometry that, within the limiting case on the exponents, allows for a reduction of the oscillation of the solution.
This project concerned the construction, testing and analysis of computational algorithms for solving parameterized and stochastic partial differential equations. The study and understanding of equations of this type is of fundamental importance in numerous engineering and scientific applications. Examples include simulation of plasma dynamics in models of electric propulsion and nuclear fusion, simulation of multiphase flows, such as the flow of water, gas and oil in reservoirs, and structural analysis of the dependence of structures on materials. Parametrization is used in such settings when properties of the models such as viscosity of fluids or electric resistivity of materials are not precisely understood and instead are treated as random variables. The resulting solutions are themselves random, and having such solutions will enable engineers to use probabilistic methods to assess the likelihood of events, for example, whether a pollutant in a liquid will exceed a limit, and to use such analyses to develop ways to ensure positive outcomes. Construction of accurate (high resolution) computational solutions is expensive, requiring significant computer time and computational resources, and there is need to reduce computational cost to make simulation useful and effective. The aim of the project was to construct algorithms to efficiently compute surrogate solutions to parameterized problems to allow for efficient and accurate simulation. The technical approach used focused on two related strategies, based on rank-reduction methods and reduced-order models. These methods construct surrogate solutions of parameter-dependent models by projection or interpolation into low-dimensional approximation spaces. Cost savings are achieved if the low-dimensional spaces can be identified and constructed efficiently and if the resulting low-dimensional algebraic systems can be solved cheaply. Accomplishments include: Theoretical and empirical demonstration of the effectiveness of fast multigrid solution strategies for computing low-rank representations of parameter-dependent solutions to discrete partial differential equations, including the first proof establishing so-called textbook convergence properties for low-rank methods. Development of efficient solution algorithms for solving nonlinear parameter-dependent partial differential equations used in models of fluid dynamics. Developent of efficient algorithms for low-rank representation of solutions of time-dependent simulations of fluid dynamics using multi-dimensional tensor representations of solutions.
Exact solutions of the following theoretical problem are presented: A plane atmosphere is in hydrostatic equilibrium with a uniform gravity. The ideal gas law is assumed. Heat is generated everywhere at a rate proportional to the local density. The atmosphere is maintained in a steady state through cooling by thermal conduction and radiation. This problem is reducible to quadratures for a thermal conductivity which is an arbitrary, but prescribed, function of the temperature, and for a radiative loss which is expressible as the product of the density and an arbitrary, but prescribed, function of the pressure. The analysis is carried out for the case of power law thermal conductivity, and a radiative loss proportional to the square of the density and to the first power of the temperature. The radiative cooling function adopted here has the basic mathematical form for an optically thin medium. The solutions reproduce the macroscopic ordering of a hot 'corona' separated from a 'photosphere' by a layer of temperature minimum. The analytic solutions allow direct illustration of the interplay between steady energy transport and the requirements of hydrostatic equilibrium.
In several recent studies, it has been proposed that the fundamental understanding of penetrant transport in dense polymer membranes occurring via the solution-diffusion model, which has been the generally accepted theoretical framework for describing penetrant transport in such materials for the past several decades, is flawed. An alternate mechanistic framework based on the idea of two-phase flow in a porous medium (i.e., pore-flow) has been broadly advanced instead, with proponents of this approach claiming that the pore-flow theoretical framework provides the necessary mechanistic insight to design novel polymeric membrane materials for emerging applications. In this study, we show experimental results for hydraulic permeation of water that are entirely consistent with the solution-diffusion theory, without modification, for three dense polymeric membranes: crosslinked poly(ethylene glycol diacrylate) (XLPEGDA), Nafion 117 ionomer in the sodium counterion form (Nafion 117-Na), and cellulose acetate (CA). By measuring water flux at transmembrane pressures up to 240 bar, we observe a nonlinear relationship between the transmembrane pressure (TMP) and water flux, J w , for XLPEGDA and Nafion 117-Na, while this relationship is linear for CA. We demonstrate that the behavior of these three materials is described via the solution-diffusion model. According to the solution-diffusion model, flux is, to a good approximation, proportional to the transmembrane concentration difference induced by the pressure difference across the membrane, rather than to TMP itself. Water sorption isotherms are reported for all three materials. They further justify the nonlinear relationship between TMP and J w observed in XLPEGDA and Nafion 117-Na, emphasizing that the nonlinearity in the flux/TMP relationship stems from nonlinearities in the sorption isotherm with pressure. Additionally, the relationship between water flux and TMP can be predicted, a priori, with no adjustable parameters when a predictive model for the diffusion coefficient of water is employed in conjunction with the experimental water sorption isotherms in the solution-diffusion model. Furthermore, our results demonstrate the validity of the solution-diffusion model to describe transport of penetrants in dense polymer membranes, while highlighting the sensitivity of the solution-diffusion model to the many physical and mathematical simplifications commonly applied to the theory in literature.
We present and analyze several definitions of Pareto optimality for multicriteria optimization or decision problems with uncertainty primarily in their objective function values. In comparison to related notions of Pareto robustness, we first provide a full characterization of an alternative efficient set hierarchy that is based on six different ordering relations both with respect to the multiple objectives and a possibly finite, countably infinite or uncountable number of scenarios. We then establish several scalarization results for the generation of the corresponding efficient points using generalized weighted-sum and epsilon-constraint techniques. Finally, we leverage these scalarization results to also derive more general conditions for the existence of efficient points in each of the corresponding optimality classes, under suitable assumptions.
Partial differential equations (PDEs) play a central role in the mathematical analysis and modeling of complex dynamic processes across all corners of science and engineering. Their solution often requires laborious analytical or computational tools, associated with a cost that is markedly amplified when different scenarios need to be investigated, for example, corresponding to different initial or boundary conditions, different inputs, etc. In this work, we introduce physics-informed DeepONets, a deep learning framework for learning the solution operator of arbitrary PDEs, even in the absence of any paired input-output training data. We illustrate the effectiveness of the proposed framework in rapidly predicting the solution of various types of parametric PDEs up to three orders of magnitude faster compared to conventional PDE solvers, setting a previously unexplored paradigm for modeling and simulation of nonlinear and nonequilibrium processes in science and engineering.
The linear Boltzmann equation can be cast in a form mathematically identical to the radiation-transport equation. A multigroup procedure is used to reduce the energy (or velocity) dependence of the transport equation to a series of one-speed problems. Each of these one-speed problems is equivalent to the monochromatic radiative-transfer problem, and existing software is used to solve this problem in slab geometry. The numerical code conserves particles in elastic collisions. Generic examples are provided to illustrate the applicability of this approach. Although this formalism can, in principle, be applied to a variety of test particle or linearized gas dynamics problems, it is particularly well-suited to study the thermalization of suprathermal particles interacting with a background medium when the thermal motion of the background cannot be ignored. Extensions of the formalism to include external forces and spherical geometry are also feasible.