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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Benchmark solutions for the galactic ion transport equations with spatial and energy coupling

In order to anticipate future space shielding requirements, NASA has initiated an effort to formulate computational methods to simulate radiation effects in space. As part of the program, numerical transport algorithms have been developed for the deterministic Boltzman equation describing galactic cosmic ray (GCR) interactions with matter. It thus becomes necessary to assess the accuracy of proposed deterministic algorithms. For this reason, analytical benchmark solutions to mathematically tractable galactic cosmic ray equations have recently been obtained. Even though these problems involve simplifying assumptions of the associated physics, they still contain the essential features of the basic transport processes. The solutions obtained are features of the basic transport processes. The solutions obtained are compared to results from numerical algorithms in order to ensure proper coding and to provide a measure of the accuracy of the numerical methods used in the algorithm. For the first time, mathematical methods have been applied to the galactic ion transport (GIT) equations in the straight ahead approximation with constant nuclear properties. The approach utilizes a Laplace transforms inversion yielding a closed form benchmark solution which is also computationally efficient.

Ganapol, Barry D.↗

DOE OSTI2020

Numerical methods for nonlocal and fractional models

Partial differential equations (PDEs) are used, with huge success, to model phenomena arising across all scientific and engineering disciplines. However, across an equally wide swath, there exist situations in which PDE models fail to adequately model observed phenomena or are not the best available model for that purpose. On the other hand, in many situations, nonlocal models that account for interaction occurring at a distance have been shown to more faithfully and effectively model observed phenomena that involve possible singularities and other anomalies. In this article, we consider a generic nonlocal model, beginning with a short review of its definition, the properties of its solution, its mathematical analysis, and specific concrete examples. We then provide extensive discussions about numerical methods, including finite element, finite difference, and spectral methods, for determining approximate solutions of the nonlocal models considered. In that discussion, we pay particular attention to a special class of nonlocal models that are the most widely studied in the literature, namely those involving fractional derivatives. The article ends with brief considerations of several modeling and algorithmic extensions which serve to show the wide applicability of nonlocal modeling.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1993

Using emergent order to shape a space society

A fast-growing movement in the scientific community is reshaping the way that we view the world around us. The short-hand name for this movement is 'chaos'. Chaos is a science of the global, nonlinear nature of systems. The center of this set of ideas is that simple, deterministic systems can breed complexity. Systems as complex as the human body, ecology, the mind or a human society. While it is true that simple laws can breed complexity, the other side is that complex systems can breed order. It is the latter that I will focus on in this paper. In the past, nonlinear was nearly synonymous with unsolvable because no general analytic solutions exist. Mathematically, an essential difference exists between linear and nonlinear systems. For linear systems, you just break up the complicated system into many simple pieces and patch together the separated solutions for each piece to form a solution to the full problem. In contrast, solutions to a nonlinear system cannot be added to form a new solution. The system must be treated in its full complexity. While it is true that no general analytical approach exists for reducing a complex system such as a society, it can be modeled. The technical involves a mathematical construct called phase space. In this space stable structures can appear which I use as analogies for the stable structures that appear in a complex system such as an ecology, the mind or a society. The common denominator in all of these systems is that they rely on a process called feedback loops. Feedback loops link the microscopic (individual) parts to the macroscopic (global) parts. The key, then, in shaping a space society, is in effectively using feedback loops. This paper will illustrate how one can model a space society by using methods that chaoticists have developed over the last hundred years. And I will show that common threads exist in the modeling of biological, economical, philosophical, and sociological systems.

Graps, Amara L.↗

DOE OSTI2020

Numerical methods for nonlocal and fractional models

Partial differential equations (PDEs) are used with huge success to model phenomena across all scientific and engineering disciplines. However, across an equally wide swath, there exist situations in which PDEs fail to adequately model observed phenomena, or are not the best available model for that purpose. On the other hand, in many situations, nonlocal models that account for interaction occurring at a distance have been shown to more faithfully and effectively model observed phenomena that involve possible singularities and other anomalies. Here, we consider a generic nonlocal model, beginning with a short review of its definition, the properties of its solution, its mathematical analysis and of specific concrete examples. We then provide extensive discussions about numerical methods, including finite element, finite difference and spectral methods, for determining approximate solutions of the nonlocal models considered. In that discussion, we pay particular attention to a special class of nonlocal models that are the most widely studied in the literature, namely those involving fractional derivatives. The article ends with brief considerations of several modelling and algorithmic extensions, which serve to show the wide applicability of nonlocal modelling.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Using Complex Probability Amplitudes to Simulate Solute Transport in Composite Porous Media

Probability amplitudes are fundamental to quantum mechanics and offer robust descriptions of complicated systems, which have allowed physicists to explain behaviors inaccessible to classical physics. This article ponders how some of the same conceptual underpinnings of the mathematics used for modeling quantum systems might be applied to subsurface water resources problems and speculates how these tools could facilitate applications on quantum computers. A probability amplitude-based model for describing advective-dispersive transport in porous media using linear operators is investigated. The proposed complex valued model decomposes spreading into two “sub-continuum partial dispersion” coefficients, and this recovers classical spreading when the sum of these coefficients is the Fickian dispersion coefficient. However, the probability amplitudes have a manyto-one relationship to a probability distribution, so it embeds a level of heterogeneity into seemingly equivalent functions. Two propagators with different sub-continuum coefficients may have the same macroscopic behavior when either is considered in isolation, but when they act on the other the system’s behavior changes. Additionally, differences in the amplitudes cause a reduction in spreading as velocity correlations are disrupted, despite both propagators having identical dispersion coefficients, and this cannot be achieved using classical methods without changing the dispersion coefficient. The main point is that these amplitude-based models offer a way to embed information about the system into the propagators, instead of just “averaging it out” when making an upscaled model.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

New results on asymmetric thick branes

This work deals with the presence and stability of thick brane solutions in the warped five dimensional braneworld scenario with a single extra spatial dimension of infinite extent. We combine two distinct procedures that give rise to new possibilities, allowing that we describe models of asymmetric thick branes, with the asymmetry being controlled by a single real parameter. We illustrate the main results with some distinct models, which show that the method works for both standard and generalized models, and the solutions are gravitationally stable against small perturbations of the metric.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1990

Extension of transonic flow computational concepts in the analysis of cavitated bearings

An analogy between the mathematical modeling of transonic potential flow and the flow in a cavitating bearing is described. Based on the similarities, characteristics of the cavitated region and jump conditions across the film reformation and rupture fronts are developed using the method of weak solutions. The mathematical analogy is extended by utilizing a few computational concepts of transonic flow to numerically model the cavitating bearing. Methods of shock fitting and shock capturing are discussed. Various procedures used in transonic flow computations are adapted to bearing cavitation applications, for example, type differencing, grid transformation, an approximate factorization technique, and Newton's iteration method. These concepts have proved to be successful and have vastly improved the efficiency of numerical modeling of cavitated bearings.

Vijayaraghavan, D.↗

NASA NTRS1990

Extension of transonic flow computational concepts in the analysis of cavitated bearings

Vijayaraghavan, D.↗

NASA NTRS1983

An efficient computational method for characterizing the effects of random surface errors on the average power pattern of reflectors

Based on the works of Ruze (1966) and Vu (1969), a novel mathematical model has been developed to determine efficiently the average power pattern degradations caused by random surface errors. In this model, both nonuniform root mean square (rms) surface errors and nonuniform illumination functions are employed. In addition, the model incorporates the dependence on F/D in the construction of the solution. The mathematical foundation of the model rests on the assumption that in each prescribed annular region of the antenna, the geometrical rms surface value is known. It is shown that closed-form expressions can then be derived, which result in a very efficient computational method for the average power pattern. Detailed parametric studies are performed with these expressions to determine the effects of different random errors and illumination tapers on parameters such as gain loss and sidelobe levels. The results clearly demonstrate that as sidelobe levels decrease, their dependence on the surface rms/wavelength becomes much stronger and, for a specified tolerance level, a considerably smaller rms/wavelength is required to maintain the low sidelobes within the required bounds.

Rahmat-Samii, Y.↗

DOE OSTI2019

The exact discontinuity of a partial wave along the left-hand cut and the exact N∕D method in non-relativistic scattering

We first deduce the analytical continuation in the complex planes of the initial and final three-momenta of the Lippmann–Schwinger equation in coupled or uncoupled partial-wave amplitudes. This result allows us to deduce a master equation whose solution is the exact discontinuity of the on-shell partial-wave amplitudes along the left-hand cut. This equation is always a linear non-singular integral equation whose solution is fixed exclusively by the knowledge of the potential, applicable to either regular or singular potentials. The capability of calculating exactly this discontinuity allows one to settle the exact N∕D method in two-body non-relativistic scattering for coupled and uncoupled waves. We exemplify this new advance in scattering theory by explicitly checking the agreement between the Lippmann–Schwinger equation with the corresponding solutions of the exact N∕D method for some examples that involve regular and singular potentials, either attractive or repulsive.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1998

Combined Numerical/Analytical Perturbation Solutions of the Navier-Stokes Equations for Aerodynamic Ejector/Mixer Nozzle Flows

In spite of rapid advances in both scalar and parallel computational tools, the large number of variables involved in both design and inverse problems make the use of sophisticated fluid flow models impractical, With this restriction, it is concluded that an important family of methods for mathematical/computational development are reduced or approximate fluid flow models. In this study a combined perturbation/numerical modeling methodology is developed which provides a rigorously derived family of solutions. The mathematical model is computationally more efficient than classical boundary layer but provides important two-dimensional information not available using quasi-1-d approaches. An additional strength of the current methodology is its ability to locally predict static pressure fields in a manner analogous to more sophisticated parabolized Navier Stokes (PNS) formulations. To resolve singular behavior, the model utilizes classical analytical solution techniques. Hence, analytical methods have been combined with efficient numerical methods to yield an efficient hybrid fluid flow model. In particular, the main objective of this research has been to develop a system of analytical and numerical ejector/mixer nozzle models, which require minimal empirical input. A computer code, DREA Differential Reduced Ejector/mixer Analysis has been developed with the ability to run sufficiently fast so that it may be used either as a subroutine or called by an design optimization routine. Models are of direct use to the High Speed Civil Transport Program (a joint government/industry project seeking to develop an economically.viable U.S. commercial supersonic transport vehicle) and are currently being adopted by both NASA and industry. Experimental validation of these models is provided by comparison to results obtained from open literature and Limited Exclusive Right Distribution (LERD) sources, as well as dedicated experiments performed at Texas A&M. These experiments have been performed using a hydraulic/gas flow analog. Results of comparisons of DREA computations with experimental data, which include entrainment, thrust, and local profile information, are overall good. Computational time studies indicate that DREA provides considerably more information at a lower computational cost than contemporary ejector nozzle design models. Finally. physical limitations of the method, deviations from experimental data, potential improvements and alternative formulations are described. This report represents closure to the NASA Graduate Researchers Program. Versions of the DREA code and a user's guide may be obtained from the NASA Lewis Research Center.

DeChant, Lawrence Justin↗

DOE OSTI2020

Electromagnetic waves in an optical photonic lattice

Laser light propagation in a photonic lattice consisting of two parallel waveguide arrays is theoretically studied using the coupled mode method, with the interaction of each waveguide with the nearest neighbours and between the waveguides of the arrays being taken into account. Analytical expressions are obtained that make it possible to accurately predict the presence of localisation of light depending on the coupling constants. Particular solutions of a system of coupled waves are found, which describe strongly localised light propagating without transverse diffraction along the entire lattice at certain values of the coupling constants. The emergence of spatially limited transverse diffraction of light is predicted. (integral optics)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Control of the optical properties of gum and dentin tissue of a human tooth at laser spectral lines in the range of 200 – 800 nm

A change in the optical properties of the tissue of the gums and dentin of a human tooth is experimentally investigated in vitro with the application of a 99.5 % glycerol solution. The kinetics of the total transmittance in the spectral range of 200 – 800 nm on characteristic laser lines is studied. It is shown that the effectiveness of optical clearing, expressed as a relative increase in the total transmittance, for gum tissue in the laser radiation range of 200 – 800 nm can vary significantly (by an order of magnitude). In this case, the absolute values of optical transmittance in the UV region are relatively small due to the strong absorption by the endogenous chromophores of the studied tissues, and in the visible and NIR regions, on the contrary, they are significant, which is due to the absence of strongly absorbing chromophores. The effectiveness of dentin clearing in the entire visible and NIR spectral regions is on average 25 % for 300 minutes of immersion of the sample in a 99.5 % glycerol solution. The presented results are necessary for increasing the efficiency of laser therapy and methods for optical monitoring of morphological changes in biological tissues for the purpose of early diagnosis of diseases. (paper)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

Extending the operation range of a phase-sensitive optical time-domain reflectometer by using fibre with chirped Bragg gratings

>A record long external-disturbance detection distance of 140 km without employing optical amplifiers in the fibre line is demonstrated by a phase-sensitive optical time-domain reflectometer (Φ-OTDR) using a fibre section with inscribed chirped fibre Bragg gratings (CFBGs). This solution provides simultaneous operation of Φ-OTDR with a standard single-mode fibre at a distance of up to 100 km and with a CFBG fibre section located from the fibre line input point by at least 140 km. (paper)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

The Improved Orbit Solution and Component Masses of bet LMi

bet LMi is a double-lined visual binary with an orbital period of ∼39 yr. Via a simultaneous fitting to both astrometric and radial velocity measurements, we give a complete and improved orbit solution with high precision. Then, the component masses are precisely determined as 2.98 ± 0.10 M {sub ⊙} and 1.92 ± 0.04 M {sub ⊙} with a relative precision of ∼3%, respectively. The orbital parallax is determined to be 19.6 ± 0.2 mas, which is two times more precise than Hipparcos parallax. With the known apparent magnitudes and magnitude difference of the components, we derive the luminosity of the components as 50.7 ± 1.8 L {sub ⊙} and 9.1 ± 4.1 L {sub ⊙}. The estimated radii of the components are 9.4 ± 0.3 R {sub ⊙} and 3.7 ± 1.5 R {sub ⊙}.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1971

Application of Monte Carlo techniques to optimization of high-energy beam transport in a stochastic environment

An algorithm employing a modified sequential random perturbation, or creeping random search, was applied to the problem of optimizing the parameters of a high-energy beam transport system. The stochastic solution of the mathematical model for first-order magnetic-field expansion allows the inclusion of state-variable constraints, and the inclusion of parameter constraints allowed by the method of algorithm application eliminates the possibility of infeasible solutions. The mathematical model and the algorithm were programmed for a real-time simulation facility; thus, two important features are provided to the beam designer: (1) a strong degree of man-machine communication (even to the extent of bypassing the algorithm and applying analog-matching techniques), and (2) extensive graphics for displaying information concerning both algorithm operation and transport-system behavior. Chromatic aberration was also included in the mathematical model and in the optimization process. Results presented show this method as yielding better solutions (in terms of resolutions) to the particular problem than those of a standard analog program as well as demonstrating flexibility, in terms of elements, constraints, and chromatic aberration, allowed by user interaction with both the algorithm and the stochastic model. Example of slit usage and a limited comparison of predicted results and actual results obtained with a 600 MeV cyclotron are given.

Parrish, R. V.↗

DOE OSTI2020

Relating the depth of the water table to the depth of weathering

Weathering of bedrock creates and occludes permeability, affecting subsurface water flow. Often, weathering intensifies above the water table. On the other hand, weathering can also commence below the water table. To explore relationships between weathering and the water table, a simplified weathering model for an eroding hillslope was formulated that takes into account both saturated and unsaturated subsurface water flow (but does not fully account for changes in dissolved gas chemistry). The phreatic line was calculated using solutions to mathematical treatments for both zones. In the model, the infiltration rate at the hill surface sets both the original and the eventual steady-state position of the water table with respect to the weathering reaction front. Depending upon parameters, the weathering front can locate either above or below the water table at steady state. Erosion also affects the water table position by changing porosity and permeability even when other hydrological conditions (e.g. hydraulic conductivity of parent material, infiltration rate at the surface) do not change. The total porosity in a hill (water storage capacity) was found to increase with infiltration rate (all else held constant). Furthermore, this effect was diminished by increasing the erosion rate. We also show examples of how the infiltration rate affects the position of the water table and how infiltration rate affects weathering advance.

58 GEOSCIENCES↗

DOE OSTI2022

It’s Not What You Take Up, It’s What You Keep: How Discoveries from Diverse Disciplines Directed the Development of the FDG PET/CT Scan

Although imaging glucose metabolism with positron emission tomography combined with X-ray CT (FDG-PET/CT) has become a standard diagnostic modality for the discovery and surveillance of malignant tumors and inflammatory processes, its origins extend back to more than a century of notable discoveries in the fields of inorganic and organic chemistry, nuclear physics, mathematics, biochemistry, solute transport physiology, metabolism, and imaging, accomplished by pioneering and driven investigators, of whom at least ten were recipients of the Nobel Prize. These tangled and diverse roots eventually coalesced into the FDG-PET/CT method, that through its many favorable characteristics inherent in the isotope used ( 18 F), the accurate imaging derived from coincidence detection of positron annihilation radiation combined with computed tomography, and the metabolic trapping of 2-deoxy-2-[ 18 F] fluoro-d-glucose (FDG) in tissues, provides safety, sensitivity, and specificity for tumor and inflammation detection. Here, the authors hope that this article will increase the appreciation among its readers of the insight, creativity, persistence, and drive of the many investigators who made this technique possible. This article is followed by a review of the many applications of FDG-PET/CT to the gastrointestinal tract and hepatobiliary system (Mandelkern in Dig Dis Sci 2022)

2-DG↗