Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mathematics and Computing, Energy - Solar”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Containerless Ripple Turbulence

One of the longest standing unsolved problems in physics relates to the behavior of fluids that are driven far from equilibrium such as occurs when they become turbulent due to fast flow through a grid or tidal motions. In turbulent flows the distribution of vortex energy as a function of the inverse length scale [or wavenumber 'k'] of motion is proportional to 1/k(sup 5/3) which is the celebrated law of Kolmogorov. Although this law gives a good description of the average motion, fluctuations around the average are huge. This stands in contrast with thermally activated motion where large fluctuations around thermal equilibrium are highly unfavorable. The problem of turbulence is the problem of understanding why large fluctuations are so prevalent which is also called the problem of 'intermittency'. Turbulence is a remarkable problem in that its solution sits simultaneously at the forefront of physics, mathematics, engineering and computer science. A recent conference [March 2002] on 'Statistical Hydrodynamics' organized by the Los Alamos Laboratory Center for Nonlinear Studies brought together researchers in all of these fields. Although turbulence is generally thought to be described by the Navier-Stokes Equations of fluid mechanics the solution as well as its existence has eluded researchers for over 100 years. In fact proof of the existence of such a solution qualifies for a 1 M$ millennium prize. As part of our NASA funded research we have proposed building a bridge between vortex turbulence and wave turbulence. The latter occurs when high amplitude waves of various wavelengths are allowed to mutually interact in a fluid. In particular we have proposed measuring the interaction of ripples [capillary waves] that run around on the surface of a fluid sphere suspended in a microgravity environment. The problem of ripple turbulence poses similar mathematical challenges to the problem of vortex turbulence. The waves can have a high amplitude and a strong nonlinear interaction. Furthermore, the steady state distribution of energy again follows a Kolmogorov scaling law; in this case the ripple energy is distributed according to 1/k (sup 7/4). Again, in parallel with vortex turbulence ripple turbulence exhibits intermittency. The problem of ripple turbulence presents an experimental opportunity to generate data in a controlled, benchmarked system. In particular the surface of a sphere is an ideal environment to study ripple turbulence. Waves run around the sphere and interact with each other, and the effect of walls is eliminated. In microgravity this state can be realized for over 2 decades of frequency. Wave turbulence is a physically relevant problem in its own right. It has been studied on the surface of liquid hydrogen and its application to Alfven waves in space is a source of debate. Of course, application of wave turbulence perspectives to ocean waves has been a major success. The experiment which we plan to run in microgravity is conceptually straightforward. Ripples are excited on the surface of a spherical drop of fluid and then their amplitude is recorded with appropriate photography. A key challenge is posed by the need to stably position a 10cm diameter sphere of water in microgravity. Two methods are being developed. Orbitec is using controlled puffs of air from at least 6 independent directions to provided the positioning force. This approach has actually succeeded to position and stabilize a 4cm sphere during a KC 135 segment. Guigne International is using the radiation pressure of high frequency sound. These transducers have been organized into a device in the shape of a dodecahedron. This apparatus 'SPACE DRUMS' has already been approved for use for combustion synthesis experiments on the International Space Station. A key opportunity presented by the ripple turbulence data is its use in driving the development of codes to simulate its properties.

Putterman, Seth↗
NASA NTRS1992

Theory of hydrophobicity: transient cavities in molecular liquids

Observation of the size distribution of transient cavities in computer simulations of water, n-hexane, and n-dodecane under benchtop conditions shows that the sizes of cavities are more sharply defined in liquid water but the most-probable-size cavities are about the same size in each of these liquids. The calculated solvent atomic density in contact with these cavities shows that water applies more force per unit area of cavity surface than do the hydrocarbon liquids. This contact density, or "squeezing" force, reaches a maximum near cavity diameters of 2.4 angstroms. The results for liquid water are compared to the predictions of simple theories and, in addition, to results for a reference simple liquid. The numerical data for water at a range of temperatures are analyzed to extract a surface free energy contribution to the work of formation of atomic-size cavities. Comparison with the liquid-vapor interfacial tensions of the model liquids studied here indicates that the surface free energies extracted for atomic-size cavities cannot be accurately identified with the macroscopic surface tensions of the systems.

NASA Center ARC↗
NASA NTRS1982

Creep crack-growth: A new path-independent integral (T sub c), and computational studies

The development of valid creep fracture criteria is considered. Two path-independent integral parameters which show some degree of promise are the C* and (Delta T)sub c integrals. The mathematical aspects of these parameters are reviewed by deriving generalized vector forms of the parameters using conservation laws which are valid for arbitrary, three dimensional, cracked bodies with crack surface tractions (or applied displacements), body forces, inertial effects, and large deformations. Two principal conclusions are that (Delta T)sub c has an energy rate interpretation whereas C* does not. The development and application of fracture criteria often involves the solution of boundary/initial value problems associated with deformation and stresses. The finite element method is used for this purpose. An efficient, small displacement, infinitesimal strain, displacement based finite element model is specialized to two dimensional plane stress and plane strain and to power law creep constitutive relations. A mesh shifting/remeshing procedure is used for simulating crack growth. The model is implemented with the quartz-point node technique and also with specially developed, conforming, crack-tip singularity elements which provide for the r to the n-(1+n) power strain singularity associated with the HRR crack-tip field. Comparisons are made with a variety of analytical solutions and alternate numerical solutions for a number of problems.

Stonesifer, R. B.↗
NASA NTRS1992

Computer Program For Variable-Conductance Heat Pipes

VCHPDA provides accurate mathematical models of transient as well as steady-state performance of variable-conductance heat pipes over wide range of operating conditions. Applies to heat pipes with either cold, wicked or hot, nonwicked gas reservoirs and uses ideal-gas law and "flat-front" (negligible vapor diffusion) gas theory. Calculates length of gas-blocked region and temperature of vapor in active portion of heat pipe by solving set of nonlinear equations for conservation of energy and mass. Written in FORTRAN 77.

Antoniuk, D.↗
NASA NTRS1993

Space Shuttle Main Engine performance analysis

For a number of years, NASA has relied primarily upon periodically updated versions of Rocketdyne's power balance model (PBM) to provide space shuttle main engine (SSME) steady-state performance prediction. A recent computational study indicated that PBM predictions do not satisfy fundamental energy conservation principles. More recently, SSME test results provided by the Technology Test Bed (TTB) program have indicated significant discrepancies between PBM flow and temperature predictions and TTB observations. Results of these investigations have diminished confidence in the predictions provided by PBM, and motivated the development of new computational tools for supporting SSME performance analysis. A multivariate least squares regression algorithm was developed and implemented during this effort in order to efficiently characterize TTB data. This procedure, called the 'gains model,' was used to approximate the variation of SSME performance parameters such as flow rate, pressure, temperature, speed, and assorted hardware characteristics in terms of six assumed independent influences. These six influences were engine power level, mixture ratio, fuel inlet pressure and temperature, and oxidizer inlet pressure and temperature. A BFGS optimization algorithm provided the base procedure for determining regression coefficients for both linear and full quadratic approximations of parameter variation. Statistical information relative to data deviation from regression derived relations was also computed. A new strategy for integrating test data with theoretical performance prediction was also investigated. The current integration procedure employed by PBM treats test data as pristine and adjusts hardware characteristics in a heuristic manner to achieve engine balance. Within PBM, this integration procedure is called 'data reduction.' By contrast, the new data integration procedure, termed 'reconciliation,' uses mathematical optimization techniques, and requires both measurement and balance uncertainty estimates. The reconciler attempts to select operational parameters that minimize the difference between theoretical prediction and observation. Selected values are further constrained to fall within measurement uncertainty limits and to satisfy fundamental physical relations (mass conservation, energy conservation, pressure drop relations, etc.) within uncertainty estimates for all SSME subsystems. The parameter selection problem described above is a traditional nonlinear programming problem. The reconciler employs a mixed penalty method to determine optimum values of SSME operating parameters associated with this problem formulation.

Santi, L. Michael↗
NASA NTRS2002

Computational Model of the Modulation of Gene Expression Following DNA Damage

High linear energy transfer (LET) radiation, such as heavy ions or neutrons, has an increased biological effectiveness compared to X rays for gene mutation, genomic instability, and carcinogenesis. In the traditional paradigm, mutations or chromosomal aberrations are causative of late effects. However, in recent years experimental evidence has demonstrated the important role of the description of the modification of gene expression by radiation in understanding the mechanisms of radiation action. In this report, approaches are discussed to the mathematical description of mRNA and protein expression kinetics following DNA damage. Several hypotheses for models of radiation modulation of protein expression are discussed including possible non-linear processes that evolve from the linear dose responses that follow the initial DNA damage produced by radiation.

NASA Center JSC↗
NASA NTRS1989

Slot clock recovery in optical PPM communication systems with avalanche photodiode photodetectors

Slot timing recovery in a direct-detection optical PPM communication system can be achieved by processing the photodetector output waveform with a nonlinear device whose output forms the input to a phase-locked loop. The choice of a simple transition detector as the nonlinearity is shown to give satisfactory synchronization performance. The rms phase error of the recovered slot clock and the effect of slot timing jitter on the bit error probability were directly measured. The experimental system consisted of an AlGaAs laser diode (wavelength = 834 nm) and a silicon avalanche photodiode photodetector. The system used Q = 4 PPM signaling and operated at a source data rate of 25 Mbits/s. The mathematical model developed to compute the rms phase error of the recovered clock is shown to be in good agreement with results of actual measurements of phase errors. The use of the recovered slot clock in the receiver resulted in no significant degradation in receiver sensitivity compared to a system with perfect slot timing. The system achieved a bit error probability of 10 to the -6th at a received optical signal energy of 55 detected photons per information bit.

Davidson, Frederic M.↗
NASA NTRS1990

Computational experience with a three-dimensional rotary engine combustion model

A new computer code was developed to analyze the chemically reactive flow and spray combustion processes occurring inside a stratified-charge rotary engine. Mathematical and numerical details of the new code were recently described by the present authors. The results are presented of limited, initial computational trials as a first step in a long-term assessment/validation process. The engine configuration studied was chosen to approximate existing rotary engine flow visualization and hot firing test rigs. Typical results include: (1) pressure and temperature histories, (2) torque generated by the nonuniform pressure distribution within the chamber, (3) energy release rates, and (4) various flow-related phenomena. These are discussed and compared with other predictions reported in the literature. The adequacy or need for improvement in the spray/combustion models and the need for incorporating an appropriate turbulence model are also discussed.

Raju, M. S.↗
NASA NTRS1974

A system for spacecraft attitude control and energy storage

A conceptual design for a double-gimbal reaction-wheel energy-wheel device which has three-axis attitude control and electrical energy storage capability is given. A mathematical model for the three-axis gyroscope (TAG) was developed, and a system of multiple units is proposed for attitude control and energy storage for a class of spacecraft. Control laws were derived to provide the required attitude-control torques and energy transfer while minimizing functions of TAG gimbal angles, gimbal rates, reaction-wheel speeds, and energy-wheel speed differences. A control law is also presented for a magnetic torquer desaturation system. A computer simulation of a three-TAG system for an orbiting telescope was used to evaluate the concept. The results of the study indicate that all control and power requirements can be satisfied by using the TAG concept.

Shaughnessy, J. D.↗
NASA NTRS1961

THEORETICAL STUDY OF ATOMIC AND MOLECULAR GASES AND THEIR REACTIONS IN THE UPPER ATMOSPHERE. SUMMARY OF WORK DONE IN THE PERIOD NOV. 1, 1960 TO APR. 30, 1961

In the previous report it was stated that the major mathematical problems of computations were nearly solved and an attempt is being made to start application of "Tables of Integrals" to the computation of the different states of atoms and molecules. In this respect I like to report that the computation of a number of states of atoms and molecules has been carried out. The results are so far confined only to the equilibrium distance of molecules. We are extending the calculation to other distances so as to compute potential energy curves of different states of molecules. We are also developing Thomas Fermi Statistical model of the atom to verify our results of more accurate calculation. The results so far have been very encouraging and I hope to report some of these results in the next paper.

Gas↗
NASA NTRS1981

Oblique-incidence secondary emission from charged dielectronics

Experimental measurements and computer simulation of secondary electron emission coefficients for FEP-Teflon, for normal and oblique incidence in the presence of a normal electric field are reported. Knowledge of the electrostatic environment surrounding the specimen calculation of particle trajectories are considered. A simulation using a conformal mapping, a Green's integral, and a trajectory generator provides the necessary mathematical support for the measurements which have been made with normal fields of 1.5 and 2.7 kV/mm. When incidence is normal and energy exceeds the critical energy, the coefficient is given by V.58 and for oblique incidence this expression may be divided by the cosine of the angle. The parameter V sub 0 is a function of normal field. Measurements for values of V sub f are presented.

Robinson, J. W.↗
NASA NTRS1996

Dynamic control of photosynthetic photon flux for lettuce production in CELSS

A new dynamic control of photosynthetic photon flux (PPF) was tested using lettuce canopies growing in the Minitron II plant-growth/canopy gas-exchange system. Canopy photosynthetic rates (Pn) were measured in real time and fedback for further environment control. Pn can be manipulated by changing PPF, which is a good environmental parameter for dynamic control of crop production in a Controlled Ecological Life-Support Systems CELSS. Decision making that combines empirical mathematical models with rule sets developed from recent experimental data was tested. With comparable yield indices and potential for energy savings, dynamic control strategies will contribute greatly to the sustainability of space-deployed CELSS.

NASA Discipline Life Support Systems↗
NASA NTRS1991

One-way nesting for a primitive equation ocean model

Prognostic numerical models for atmospheric and oceanic circulations require initial fields, boundary conditions, and forcing functions in addition to a consistent set of partial differential equations, including a state relation and equations expressing conservation of mass, momentum, and energy. Depending on the horizontal domain to be modeled, the horizontal boundary conditions are either physically obvious or extremely difficult to specify consistently. If the entire atmosphere is modeled, periodic horizontal boundary conditions are appropriate. On the other hand, the physical horizontal boundaries on the entire ocean are solid walls. Obviously, the normal velocity at a solid wall is zero while the specification of the tangential velocity depends on the mathematical treatment of the horizontal viscous terms. Limitations imposed by computer capacity and cost, as well as research interests, have led to the use of limited area models to study flows in the atmosphere and ocean. The limited area models do not have physical horizontal boundaries, merely numerical ones. Correctly determining these open boundary conditions for limited-area numerical models has both intrigued and frustrated numerical modelers for decades. One common approach is to use the closed or solid wall boundary conditions for a limited-area model. The argument given for this approach is that the boundary conditions affect flow near the walls but that none of these effects are propagated into the interior. Therefore, one chooses a big enough domain that the central region of interest is not corrupted by the boundary flow. Research in progress to model the North Atlantic circulation vividly illustrates the pitfalls of this approach. Two model runs are compared: (1) the southern boundary at 20S between latitudes 0 and 40W is artificially closed; and (2) the same boundary is specified as open with an inward transport of 15 Sv (determined from a global model with the same physics) uniformly spread across the boundary. A comparison of both runs is presented.

Blake, D. W.↗
NASA NTRS1998

Analyzing Space-Based Interferometric Measurements of Stars and Network Measurements of Gamma-Ray Bursts

Since the announcement of the discovery of sources of bursts of gamma-ray radiation in 1973, hundreds more reports of such bursts have now been published. Numerous artificial satellites have been equipped with gamma-ray detectors including the very successful Compton Gamma Ray Observatory BATSE instrument. Unfortunately, we have made no progress in identifying the source(s) of this high energy radiation. We suspected that this was a consequence of the method used to define gamma-ray burst source "error boxes." An alternative procedure to compute gamma-ray burst source positions, with a purely physical underpinning, was proposed in 1988 by Taff. Since then we have also made significant progress in understanding the analytical nature of the triangulation problem and in computing actual gamma-ray burst positions and their corresponding error boxes. For the former, we can now mathematically illustrate the crucial role of the area occupied by the detectors, while for the latter, the Atteia et al. (1987) catalog has been completely re-reduced. There are very few discrepancies in locations between our results and those of the customary "time difference of arrival" procedure. Thus, we have numerically demonstrated that the end result, for the positions, of these two very different-looking procedures is the same. Finally, for the first time, we provide a sample of realistic "error boxes" whose non-simple shapes vividly portray the difficulty of burst source localization.

Taff, L. G.↗
NASA NTRS2007

Displaying and Analyzing Antenna Radiation Patterns

Radiant Energy Display and Analysis Software Package (REDAP) is a computer program for processing antenna-radiation- pattern data that have been preprocessed by a data-collection program and stored in a spherical-coordinate format. REDAP is designed specifically for application to data generated in testing of antennas in an anechoic chamber at Johnson Space Center; parts of REDAP may be reusable for processing antenna-test data collected elsewhere. REDAP provides a graphical user interface (GUI) and executes mathematical and plotting routines on a personal computer.

Siekierski, Jim↗
NASA NTRS1981

Creep crack-growth: A new path-independent T sub o and computational studies

Two path independent integral parameters which show some degree of promise as fracture criteria are the C* and delta T sub c integrals. The mathematical aspects of these parameters are reviewed. This is accomplished by deriving generalized vector forms of the parameters using conservation laws which are valid for arbitrary, three dimensional, cracked bodies with crack surface tractions (or applied displacements), body forces, inertial effects and large deformations. Two principal conclusions are that delta T sub c is a valid crack tip parameter during nonsteady as well as steady state creep and that delta T sub c has an energy rate interpretation whereas C* does not. An efficient, small displacement, infinitestimal strain, displacement based finite element model is developed for general elastic/plastic material behavior. For the numerical studies, this model is specialized to two dimensional plane stress and plane strain and to power law creep constitutive relations.

Stonesifer, R. B.↗
NASA NTRS1981

Synthetic battery cycling

The use of interactive computer graphics is suggested as an aid in battery system development. Mathematical representations of simplistic but fully representative functions of many electrochemical concepts of current practical interest will permit battery level charge and discharge phenomena to be analyzed in a qualitative manner prior to the assembly and testing of actual hardware. This technique is a useful addition to the variety of tools available to the battery system designer as he bridges the gap between interesting single cell life test data and reliable energy storage subsystems.

Thaller, L. H.↗
NASA NTRSDate not supplied

Computationally Guided Development of Components for High Energy Density Solid-State Lithium-Sulfur Batteries

All electric vertical take-off and landing vehicles (eVTOL) for urban air mobility (UAM) concepts face numerous challenging technical barriers before their introduction into the consumer marketplace. The most challenging of these technical barriers to overcome is developing an energy storage system capable of meeting the rigorous aerospace safety and performance criteria1. The performance metrics for eVTOL craft, such as specific energy, specific power, and safety, exceed those of electric automobiles by a factor of two to four. Current state-of-the-art (SOA) lithium-ion batteries are incapable of meeting the key performance criteria of energy and safety for eVTOL. Therefore, next generation advanced chemistries and designs must be developed to meet required performance metrics for electric aviation2. Beyond lithium-ion chemistries, such as lithium-sulfur, show promise in their high energy, while limitations exist in their power and cyclability due to low electrical conductivity and high intermediate solubility in organic liquid electrolytes. Several strategies to overcome the low electrical conductivity involve the use of selenium as a dopant in the active sulfur material, along with the incorporation of 2-dimensional electron-conducting holey-graphene to improve the composite cathodes electronic conductivity. Furthermore, combining this chemistry with a solid-electrolyte avoids the components’ dissolution issues3. However, the development of composite solid-state cathodes is non-trivial as several components must be intimately mixed so that the active component has sufficient access to both electrons and lithium ions to undergo full electrochemical conversion. Mathematical modeling of battery components can assist experimental design through a robust and rigorous combination of computational modeling techniques covering multiple length scales. The objective is to leverage modern computational materials methods combined with battery multiphysics tools to develop radically advanced compatible cathode and electrolyte materials, build and test solid state lithium-sulfur cells and packs. A NASA-based cross-organizational team of high-powered experts combined integrated computational predictive modeling, fundamental chemistry analysis, advanced material science, and battery cell development to tackle this very challenging, multidisciplinary problem. This presentation will show a multiscale computational modeling approach that has produced a novel particle dynamics method called Solid Electrolyte Sphere Approximation Model (SESAM). SESAM modeling targets the 1-10 µm scale structures and provides electromechanical and grain interactions for predictive design guidelines for the manufacturing of solid-state components. Parameters such as particle size and volume fraction of the constituent materials were modeled and experimentally fabricated to optimize electrochemical performance through improved microstructure design. Experimental feedback was provided through ionic and electronic conductivity assessment and structural analysis of developed materials and cell components.

battery↗