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.
Mathematical Solutions Toolset is a collection of five software programs that rapidly solve some common mathematical problems. The programs consist of a set of Microsoft Excel worksheets. The programs provide for entry of input data and display of output data in a user-friendly, menu-driven format, and for automatic execution once the input data has been entered.
The availability of formaldehyde (HCHO) (a proxy for volatile organic compound reactivity) and nitrogen dioxide (NO 2 ) (a proxy for nitrogen oxides) tropospheric columns from ultraviolet–visible (UV–Vis) satellites has motivated many to use their ratios to gain some insights into the near-surface ozone sensitivity. Strong emphasis has been placed on the challenges that come with transforming what is being observed in the tropospheric column to what is actually in the planetary boundary layer (PBL) and near the surface; however, little attention has been paid to other sources of error such as chemistry, spatial representation, and retrieval uncertainties. Here we leverage a wide spectrum of tools and data to quantify those errors carefully. Concerning the chemistry error, a well-characterized box model constrained by more than 500 h of aircraft data from NASA's air quality campaigns is used to simulate the ratio of the chemical loss of HO 2 + RO 2 (LRO x ) to the chemical loss of NO x (LNO x ). Subsequently, we challenge the predictive power of HCHO/NO 2 ratios (FNRs), which are commonly applied in current research, in detecting the underlying ozone regimes by comparing them to LRO x /LNO x . FNRs show a strongly linear (R 2 =0.94) relationship with LRO x /LNO x , but only on the logarithmic scale. Following the baseline (i.e., ln(LRO x /LNO x ) = −1.0 ± 0.2) with the model and mechanism (CB06, r2) used for segregating NO x -sensitive from VOC-sensitive regimes, we observe a broad range of FNR thresholds ranging from 1 to 4. The transitioning ratios strictly follow a Gaussian distribution with a mean and standard deviation of 1.8 and 0.4, respectively. This implies that the FNR has an inherent 20 % standard error (1σ) resulting from not accurately describing the RO x –HO x cycle. We calculate high ozone production rates (PO 3 ) dominated by large HCHO × NO 2 concentration levels, a new proxy for the abundance of ozone precursors. The relationship between PO 3 and HCHO × NO 2 becomes more pronounced when moving towards NO x -sensitive regions due to nonlinear chemistry; our results indicate that there is fruitful information in the HCHO × NO 2 metric that has not been utilized in ozone studies. The vast amount of vertical information on HCHO and NO 2 concentrations from the air quality campaigns enables us to parameterize the vertical shapes of FNRs using a second-order rational function permitting an analytical solution for an altitude adjustment factor to partition the tropospheric columns into the PBL region. We propose a mathematical solution to the spatial representation error based on modeling isotropic semivariograms. Based on summertime-averaged data, the Ozone Monitoring Instrument (OMI) loses 12 % of its spatial information at its native resolution with respect to a high-resolution sensor like the TROPOspheric Monitoring Instrument (TROPOMI) (> 5.5 × 3.5 km 2 ). A pixel with a grid size of 216 km 2 fails at capturing ∼ 65 % of the spatial information in FNRs at a 50 km length scale comparable to the size of a large urban center (e.g., Los Angeles). We ultimately leverage a large suite of in situ and ground-based remote sensing measurements to draw the error distributions of daily TROPOMI and OMI tropospheric NO 2 and HCHO columns. At a 68 % confidence interval (1σ), errors pertaining to daily TROPOMI observations, either HCHO or tropospheric NO 2 columns, should be above 1.2–1.5 × 10 16 molec. cm −2 to attain a 20 %–30 % standard error in the ratio. This level of error is almost non-achievable with the OMI given its large error in HCHO. The satellite column retrieval error is the largest contributor to the total error (40 %–90 %) in the FNRs. Due to a stronger signal in cities, the total relative error (< 50 %) tends to be mild, whereas areas with low vegetation and anthropogenic sources (e.g., the Rocky Mountains) are markedly uncertain (> 100 %). Our study suggests that continuing development in the retrieval algorithm and sensor design and calibration is essential to be able to advance the application of FNRs beyond a qualitative metric.
A study is made of the mathematical solution of the differential equation of motion of a test particle in the equatorial plane of the Kerr gravitational field, using S (Schwarzschild-like) coordinates. A qualitative solution of this equation leads to the conclusion that there can only be 25 different types of orbits. For each value of a, the results are presented in a master diagram for which h and e are the parameters. A master diagram divides the h, e parameter space into regions such that at each point within one of these regions the types of admissible orbits are qualitatively the same. A pictorial representation of the physical orbits in the r, phi plane is also given.
A method is presented for the calculation of elastic stresses in symmetrical disks typical of those of a high-temperature gas turbine. The method is essentially a finite-difference solution of the equilibrium and compatibility equations for elastic stresses in a symmetrical disk. Account can be taken of point-to-point variations in disk thickness, in temperature, in elastic modulus, in coefficient of thermal expansion, in material density, and in Poisson's ratio. No numerical integration or trial-and-error procedures are involved and the computations can be performed in rapid and routine fashion by nontechnical computers with little engineering supervision. Checks on problems for which exact mathematical solutions are known indicate that the method yields results of high accuracy. Illustrative examples are presented to show the manner of treating solid disks, disks with central holes, and disks constructed either of a single material or two or more welded materials. The effect of shrink fitting is taken into account by a very simple device.
With the accelerating penetration of solar energy in energy systems, market operations concerning the solar-associated intermittency are widely discussed. How to correctly model the solar generation in the market and solve the uncertainty-based operation problems call for solutions. Unlike other renewable resources, solar power can be more easily applied to the demand side through behind-the-meter installations in a distributed manner with significant future extensibility. The corresponding operational rules in the markets are largely being revised. Over the past two decades, with the development of deregulated power markets, numerous studies related to uncertainty-based market operations tend to transform the market from a deterministic and centralized framework to a stochastic and decentralized one. As the relevant literature is most recent, this paper reviews their contributions to electricity market operations with solar integration and underlines the research directions of this topic and ongoing industry practices. A thorough review of the electricity markets worldwide with solar energy is provided. A variety of proposed mathematical solutions to the problem are also addressed.
This work describes the development of a micromechanical-based constitutive model accounting for the effect of internal expansion on the residual elasticity of a Hashin composite material. This material is made of spherical inclusions that are subjected to gradual swelling within a quasi-brittle matrix. The main focus of this work is to describe and analyze the material mechanical response, with an additional focus on the internal swelling’s effect on the stress–strain response and residual elasticity. Microstructural features and parameters of major importance for the mechanical responses were identified. The innovative characteristics of the proposed approach are summarized as follows: (1) a full determination of the physics of a complete-damage problem throughout the whole process of inclusion swelling with upscaling techniques, which transfers the microcrack-related properties from the lower scale to upper scale; and (2) an evolution of the mechanical fields and the corresponding residual elasticity for various inclusion swelling levels. Concerning the matrix–inclusion composite, it was hypothesized that only the matrix was susceptible to cracking, with varied degrees of damage, whereas the inclusions behave elastically and the elastic modulus of the expanding inclusions remains constant. It is to emphasize that the gradual swelling of inclusions is modeled by an increasing strain in the current micromechanical-based constitutive model, and the microcracks are represented by a set of randomly oriented penny-shaped microcracks with identical radii (namely, closed cracks). The main contribution of the current research is to establish the exact mathematical solutions for the mechanical fields (stress, strain) caused by the swelling of the inclusions (mechanical loading) and derive the effective residual elasticity of a composite (structural response) subject to internal expansion and quasi-brittle damage. Based on the assumption of closed cracks that could be extended to open cracks in the upcoming work, the results proposed by the present paper help to understand the non-linear mechanical behavior of quasi-brittle materials subject to microcracking and provide a theoretical framework to be used as academical benchmark for numerical simulations.
Models of the electrical double layer (EDL) at electrode/liquid-electrolyte interfaces no longer hold for all-solid-state electrochemistry. Here we show a more general model for the EDL at a solid-state electrochemical interface based on the Poisson–Fermi–Dirac equation. By combining this model with density functional theory predictions, the interconnected electronic and ionic degrees of freedom in all-solid-state batteries, including the electronic band bending and defect concentration variation in the space-charge layer, are captured self-consistently. Along with a general mathematical solution, the EDL structure is presented in various materials that are thermodynamically stable in contact with a lithium metal anode: the solid electrolyte Li 7 La 3 Zr 2 O 12 (LLZO) and the solid interlayer materials LiF, Li 2 O and Li 2 CO 3 . Here, the model further allows design of the optimum interlayer thicknesses to minimize the electrostatic barrier for lithium ion transport at relevant solid-state battery interfaces.
Approximate mathematical solutions to space vehicle flight dynamic response control problem
Mathematical solution for electron continuity equation in decaying plasma afterglow
Mathematical solution of nonlinear transverse combustion instability problems in thin annular liquid propellant rocket chambers
Objectives of study were: investigation of temperature changes caused by welding arc with analysis of temperature distribution; development of system of mathematical statements describing thermal stresses and plastic strains during welding; and development of system of mathematical solutions and computer programs for one-dimensional analysis.
The research is reported concerning the development of a system of mathematical solutions and computer programs for one- and two-dimensional analyses for thermal stresses. Reports presented include: the investigation of thermal stress and buckling of tantalum and columbium sheet; and analysis of two dimensional thermal strains and metal movement during welding.
Mathematical solution permits showing that, for circular bends in hardwalled ducts of rectangular cross sections, specific acoustic impedance depends on both frequency and geometry of bend.
A computer program was developed to calculate the three dimensional, steady, incompressible, inviscid, irrotational flow field at the propeller plane (propeller removed) located upstream of an arbitrary airframe geometry. The program uses a horseshoe vortex of known strength to model the wing. All other airframe surfaces are modeled by a network source panels of unknown strength which is exposed to a uniform free stream and the wing-induced velocity field. By satisfying boundary conditions on each panel (the Neumann problem), relaxed boundary conditions being used on certain panels to simulate inlet inflow, the source strengths are determined. From the known source and wing vortex strengths, the resulting velocity fields on the airframe surface and at the propeller plane are obtained. All program equations are derived in detail, and a brief description of the program structure is presented. A user's manual which fully documents the program is cited. Computer predictions of the flow on the surface of a sphere and at a propeller plane upstream of the sphere are compared with the exact mathematical solutions. Agreement is good, and correct program operation is verified.
Results of studies of wind shear hazard to aircraft operation are summarized. Existing wind shear profiles currently used in computer and flight simulator studies are reviewed. The governing equations of motion for an aircraft are derived incorporating the variable wind effects. Quantitative discussions of the effects of wind shear on aircraft performance are presented. These are followed by a review of mathematical solutions to both the linear and nonlinear forms of the governing equations. Solutions with and without control laws are presented. The application of detailed analysis to develop warning and detection systems based on Doppler radar measuring wind speed along the flight path is given. A number of flight path deterioration parameters are defined and evaluated. Comparison of computer-predicted flight paths with those measured in a manned flight simulator is made. Some proposed airborne and ground-based wind shear hazard warning and detection systems are reviewed. The advantages and disadvantages of both types of systems are discussed.
Computer program developed that provides mathematical solution to design and construction of combustion chambers for jet engines. Improved results in areas of combustor flow fields accomplished by this computerprogram solution, cheaper and quicker than experiments involving real systems for models.
Gravity waves and their associated breaking into turbulence are very important in producing the overall picture of middle atmosphere global dynamics and associated transport. It is shown in this research that MST radars represent a most powerful technique for obtaining the needed parameters for gravity-wave-induced drag and diffusion effects as well as measuring wave accelerations and diffusion directly. A mathematical solution to this problem is that of radiative equilibrium with a balanced thermal wind.
The coil planet centrifuge designed by Ito employs flow of a single liquid phase, through a rotating coiled tube in a centrifugal force field, to provide a separation of particles based on sedimentation rates. Mathematical solutions are derived for the linear differential equations governing particle behavior in the coil planet centrifuge device. These solutions are then applied as the basis of a model for optimizing particle separations.