Search NASA⌕ Search

SEARCH · Search NASA

Results for “MATHEMATICAL SOLUTIONS”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Heat pipe dynamic behavior

The vapor flow in a heat pipe was mathematically modeled and the equations governing the transient behavior of the core were solved numerically. The modeled vapor flow is transient, axisymmetric (or two-dimensional) compressible viscous flow in a closed chamber. The two methods of solution are described. The more promising method failed (a mixed Galerkin finite difference method) whereas a more common finite difference method was successful. Preliminary results are presented showing that multi-dimensional flows need to be treated. A model of the liquid phase of a high temperature heat pipe was developed. The model is intended to be coupled to a vapor phase model for the complete solution of the heat pipe problem. The mathematical equations are formulated consistent with physical processes while allowing a computationally efficient solution. The model simulates time dependent characteristics of concern to the liquid phase including input phase change, output heat fluxes, liquid temperatures, container temperatures, liquid velocities, and liquid pressure. Preliminary results were obtained for two heat pipe startup cases. The heat pipe studied used lithium as the working fluid and an annular wick configuration. Recommendations for implementation based on the results obtained are presented. Experimental studies were initiated using a rectangular heat pipe. Both twin beam laser holography and laser Doppler anemometry were investigated. Preliminary experiments were completed and results are reported.

Issacci, F.↗
DOE OSTI2022

The “Fresnel Equations” for Diffuse radiation on Inclined photovoltaic Surfaces (FEDIS)

Here the well-known Fresnel equations solve for the reflection and transmission of light for precise incident angles. The transmission of diffuse radiation incident on a planar or domed surface is often needed for real-world applications. Due to the complexity of the Fresnel equations, the analytical solution of the integration has hitherto been unobtainable over the last centuries. Therefore, this problem was numerically solved by integrating the angular transmittances in space often leading to substantial computing burden and bias in the results. To efficiently estimate the solar energy resource for a glass-covered photovoltaic (PV) module, we derive an analytical solution of diffuse transmission based on the rigorous integration of an alternate form of the Fresnel equations. The approach leads to a simple yet accurate relative transmittance model that reconciles the solar energy sensed by pyranometers and PV panels. With limited and clearly stated approximations, the complex mathematical derivation resulted in an elegant solution. An experiment using 1-year of data at the National Renewable Energy Laboratory's (NREL's) Solar Radiation Research Laboratory (SRRL) shows that the new model dramatically decreases the disparity between the solar radiation measurements by a Kipp and Zonen CM Pyranometer 22 (CMP22) and an IMT reference cell on a 1-axis tracking system. The solution in this paper can be widely used in scientific and engineering research, development, and applications wherever the Fresnel equations are used.

14 SOLAR ENERGY↗
DOE OSTI2022

FEDIS (The “Fresnel Equations” for Diffuse radiation on Inclined photovoltaic Surfaces) [SWR-22-64]

The well-known Fresnel equations solve for the reflection and transmission of light for precise incident angles. The transmission of diffuse radiation incident on a planar or domed surface is often needed for real-world applications. Due to the complexity of the Fresnel equations, the analytical solution of the integration has hitherto been unobtainable over the last centuries. Therefore, this problem was numerically solved by integrating the angular transmittances in space often leading to substantial computing burden and bias in the results. To efficiently estimate the solar energy resource for a glass-covered photovoltaic (PV) module, we derive an analytical solution of diffuse transmission based on the rigorous integration of an alternate form of the Fresnel equations. The approach leads to a simple yet accurate relative transmittance model that reconciles the solar energy sensed by pyranometers and PV panels. With limited and clearly stated approximations, the complex mathematical derivation resulted in an elegant solution.

Xie, Yu↗
NASA NTRS2014

Comptonization in Ultra-Strong Magnetic Fields: Numerical Solution to the Radiative Transfer Problem

We consider the radiative transfer problem in a plane-parallel slab of thermal electrons in the presence of an ultra-strong magnetic field (B approximately greater than B(sub c) approx. = 4.4 x 10(exp 13) G). Under these conditions, the magnetic field behaves like a birefringent medium for the propagating photons, and the electromagnetic radiation is split into two polarization modes, ordinary and extraordinary, that have different cross-sections. When the optical depth of the slab is large, the ordinary-mode photons are strongly Comptonized and the photon field is dominated by an isotropic component. Aims. The radiative transfer problem in strong magnetic fields presents many mathematical issues and analytical or numerical solutions can be obtained only under some given approximations. We investigate this problem both from the analytical and numerical point of view, provide a test of the previous analytical estimates, and extend these results with numerical techniques. Methods. We consider here the case of low temperature black-body photons propagating in a sub-relativistic temperature plasma, which allows us to deal with a semi-Fokker-Planck approximation of the radiative transfer equation. The problem can then be treated with the variable separation method, and we use a numerical technique to find solutions to the eigenvalue problem in the case of a singular kernel of the space operator. The singularity of the space kernel is the result of the strong angular dependence of the electron cross-section in the presence of a strong magnetic field. Results. We provide the numerical solution obtained for eigenvalues and eigenfunctions of the space operator, and the emerging Comptonization spectrum of the ordinary-mode photons for any eigenvalue of the space equation and for energies significantly lesser than the cyclotron energy, which is on the order of MeV for the intensity of the magnetic field here considered. Conclusions. We derived the specific intensity of the ordinary photons, under the approximation of large angle and large optical depth. These assumptions allow the equation to be treated using a diffusion-like approximation.

acceleration of particles↗
DOE OSTI2023

Scaling similarities and quasinormal modes of D0 black hole solutions

We study the gravity solution dual to the D0 brane quantum mechanics, or BFSS matrix model, in the ’t Hooft limit. The classical physics described by this gravity solution is invariant under a scaling transformation, which changes the action with a specific critical exponent, sometimes called the hyperscaling violating exponent. We present an argument for this critical exponent from the matrix model side, which leads to an explanation for the peculiar temperature dependence of the entropy in this theory, S ∝ T 9/5 . We also present a similar argument for all other Dp-brane geometries. We then compute the black hole quasinormal modes. This involves perturbing the finite temperature geometry. These perturbations can be easily obtained by a mathematical trick where we view the solution as the dimensional reduction of an Ad S2+9/5 × S 8 geometry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1950

Subsonic Gas Flow Past A Wing Profile

The use of the linearized equations of Chaplygin to calculate the subsonic flow of a gas permits solving the problem of the flow about a wing profile for absence and presence of circulation. The solution is obtained in a practical convenient form that permits finding all the required magnitudes for the gas flow (lift, lift moment velocity distribution over the profile, and critical Mach number). This solution is not expressed in simple closed form; for a certain simplifying assumption, however, the equations of Chaplygin can be reduced to equations with constant coefficients, and solutions are obtained by using only the mathematical apparatus of the theory of functions of a complex variable. The method for simplifying the equations was pointed out by Chaplygin himself. These applied similar equations to the solution of the flow problem and obtained a solution for the case of the absence of circulation.

Christianovich, S. A.↗
DOE OSTI2021

Elucidating the Role of Hydroxide Electrolyte on Anion-Exchange-Membrane Water Electrolyzer Performance

Many solid-state devices, especially those requiring anion conduction, often add a supporting electrolyte to enable efficient operation. The prototypical case is that of anion-exchange-membrane water electrolyzers (AEMWEs), where addition of an alkali metal solution improves performance. However, the specific mechanism of this performance improvement is currently unknown. This work investigates the functionality of the alkali metal solution in AEMWEs using experiments and mathematical models. The results show that additional hydroxide plays a key role not only in ohmic resistance of the membrane and catalyst layer but also in the reaction kinetics. The modeling suggests that the added liquid electrolyte creates an additional electrochemical interface with the electrocatalyst that provides ion-transport pathways and distributes product gas bubbles; the total effective electrochemical active surface area in the cell with 1 M KOH is 5 times higher than that of the cell with DI water. In the cell with 1 M KOH, more than 80% of the reaction current is associate with the liquid electrolyte. These results indicate the importance of high pH of electrolyte and catalyst/electrolyte interface in AEMWEs. The understanding of the functionality of the alkali metal solution presented in this study should help guide the design and optimization of AEMWEs.

25 ENERGY STORAGE↗
DOE OSTI2021

Probabilistic Deep Learning Approach to Automate the Interpretation of Multi-phase Diffraction Spectra

Autonomous synthesis and characterization of inorganic materials requires the automatic and accurate analysis of X-ray diffraction spectra. For this task, we designed a probabilistic deep learning algorithm to identify complex multi-phase mixtures. At the core of this algorithm lies an ensemble convolutional neural network trained on simulated diffraction spectra, which are systematically augmented with physics-informed perturbations to account for artifacts that can arise during experimental sample preparation and synthesis. In this work, larger perturbations associated with off-stoichiometry are also captured by supplementing the training set with hypothetical solid solutions. Spectra containing mixtures of materials are analyzed with a newly developed branching algorithm that utilizes the probabilistic nature of the neural network to explore suspected mixtures and identify the set of phases that maximize confidence in the prediction. Our model is benchmarked on simulated and experimentally measured diffraction spectra, showing exceptional performance with accuracies exceeding those given by previously reported methods based on profile matching and deep learning. We envision that the algorithm presented here may be integrated in experimental workflows to facilitate the high-throughput and autonomous discovery of inorganic materials.

36 MATERIALS SCIENCE↗
NASA NTRS1990

Convection-driven tectonics on Venus

An analysis is presented of convective stress coupling to an elastic lithosphere as applied to Venus. Theoretical solutions are introduced for the response of a mathematically thick elastic plate overlying a Newtonian viscous medium with an exponential depth dependence of viscosity, and a Green's function solution is obtained for the viscous flow driven by a harmonic density distribution at a specified depth. An elastic-plastic analysis is carried out for the deformation of a model Venus lithosphere. The results predict that dynamic uplift of Venusian topography must be accompanied by extensive brittle failure and viscous flow in the lithosphere.

Phillips, Roger J.↗
DOE OSTI2021

Analysis of a new multispecies tumor growth model coupling 3D phase-fields with a 1D vascular network

In this study, we present and analyze a mathematical model for tumor growth incorporating ECM erosion, interstitial flow, and the effect of vascular flow and nutrient transport. The model is of phase-field or diffused-interface type in which multiple phases of cell species and other constituents are separated by smooth evolving interfaces. The model involves a mesoscale version of Darcy’s law to capture the flow mechanism in the tissue matrix. Modeling flow and transport processes in the vasculature supplying the healthy and cancerous tissue, one-dimensional (1D) equations are considered. Since the models governing the transport and flow processes are defined together with cell species models on a three-dimensional (3D) domain, we obtain a 3D–1D coupled model.

3D–1D coupled blood flow models↗
DOE OSTI2021

Synchronous and Concurrent Multidomain Computing Method for Cloud Computing Platforms

We present a numerical method for synchronous and concurrent solution of transient elastodynamics problem where the computational domain is divided into subdomains that may reside on separate computational platforms. Here, this work employs the variational multiscale discontinuous Galerkin (VMDG) method to develop interdomain transmission conditions for transient problems. The fine-scale modeling concept leads to variationally consistent coupling terms at the common interfaces. The method admits a large class of time discretization schemes, and decoupling of the solution for each subdomain is achieved by selecting any explicit algorithm. Numerical tests with a manufactured solution problem show optimal convergence rates. The energy history in a free vibration problem is in agreement with that of the solution from a monolithic computational domain.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Novel deformation function creating or destroying any number of even kink solutions

Here, we present a one-parameter family of deformation functions f(Φ) which have novel properties. Firstly, the deformation function is its own inverse. We show that a class of potentials remains invariant under this deformation. Further, when applied to a certain class of kink bearing potentials, one obtains potentials which are not bounded from below. Besides, we show that there is a wide class of potentials having power law kink tails which are connected by this deformation but the corresponding kink tails of the two potentials have different asymptotic behavior. Finally, we show that when this deformation function is applied to an appropriate one-parameter family of potentials having two kink solutions, it creates new potentials with an arbitrary even number (2n) of kink solutions. Conversely, by starting from this potential with 2n kink solutions the deformation function annihilates 2n-2 kinks and we get back the potential which we started with having only two kink solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Direct Extraction of Uranium Oxides with N,N -Di(2-ethylhexyl)isobutyramide

The direct extraction of uranium from voloxidized nuclear fuel into an organic solvent offers several potential advantages over conventional hydrometallurgical reprocessing, including reducing the reprocessing plant footprint, providing an initial degree of decontamination from fission products, and minimizing the amount of secondary waste from nitric acid. In this work, the direct extraction of uranium oxides into 1.5 M N,N-di(2-ethylhexyl)isobutyramide (DEHiBA) in n-dodecane is examined. UV–vis spectra and distribution ratios of HNO 2 in 1.5 M DEHiBA as well as the equilibrium organic phase H 2 O concentrations in HNO 3 -loaded 1.5 M DEHiBA are also reported. Hypothesized reaction stoichiometries for the direct extraction of uranium from UO 2 , α-U 3 O 8 , and ε-UO 3 are verified through analysis of organic-phase U, HNO 2 , and HNO 3 concentrations after dissolution. Water generated by the dissolution results in the formation of a separate aqueous phase, which will need to be accounted for in future flowsheet design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1966

Computer program determines chemical composition of physical system at equilibrium

FORTRAN 4 digital computer program calculates equilibrium composition of complex, multiphase chemical systems. This is a free energy minimization method with solution of the problem reduced to mathematical operations, without concern for the chemistry involved. Also certain thermodynamic properties are determined as byproducts of the main calculations.

Kwong, S. S.↗