Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mathematics and Computing, Energy”

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 55 records · Page 3

A Hierarchical Approach to Fracture Mechanics

Recent research conducted under NASA LaRC's Creativity and Innovation Program has led to the development of an initial approach for a hierarchical fracture mechanics. This methodology unites failure mechanisms occurring at different length scales and provides a framework for a physics-based theory of fracture. At the nanoscale, parametric molecular dynamic simulations are used to compute the energy associated with atomic level failure mechanisms. This information is used in a mesoscale percolation model of defect coalescence to obtain statistics of fracture paths and energies through Monte Carlo simulations. The mathematical structure of predicted crack paths is described using concepts of fractal geometry. The non-integer fractal dimension relates geometric and energy measures between meso- and macroscales. For illustration, a fractal-based continuum strain energy release rate is derived for inter- and transgranular fracture in polycrystalline metals.

Saether, Erik↗

Space proton transport in one dimension

An approximate evaluation procedure is derived for a second-order theory of coupled nucleon transport in one dimension. An analytical solution with a simplified interaction model is used to determine quadrature parameters to minimize truncation error. Effects of the improved method on transport solutions with the BRYNTRN data base are evaluated. Comparisons with Monte Carlo benchmarks are given. Using different shield materials, the computational procedure is used to study the physics of space protons. A transition effect occurs in tissue near the shield interface and is most important in shields of high atomic number.

NASA Discipline Number 45-10↗

Overall Traveling-Wave-Tube Efficiency Improved By Optimized Multistage Depressed Collector Design

Depressed Collector Design The microwave traveling wave tube (TWT) is used widely for space communications and high-power airborne transmitting sources. One of the most important features in designing a TWT is overall efficiency. Yet, overall TWT efficiency is strongly dependent on the efficiency of the electron beam collector, particularly for high values of collector efficiency. For these reasons, the NASA Glenn Research Center developed an optimization algorithm based on simulated annealing to quickly design highly efficient multistage depressed collectors (MDC's). Simulated annealing is a strategy for solving highly nonlinear combinatorial optimization problems. Its major advantage over other methods is its ability to avoid becoming trapped in local minima. Simulated annealing is based on an analogy to statistical thermodynamics, specifically the physical process of annealing: heating a material to a temperature that permits many atomic rearrangements and then cooling it carefully and slowly, until it freezes into a strong, minimum-energy crystalline structure. This minimum energy crystal corresponds to the optimal solution of a mathematical optimization problem. The TWT used as a baseline for optimization was the 32-GHz, 10-W, helical TWT developed for the Cassini mission to Saturn. The method of collector analysis and design used was a 2-1/2-dimensional computational procedure that employs two types of codes, a large signal analysis code and an electron trajectory code. The large signal analysis code produces the spatial, energetic, and temporal distributions of the spent beam entering the MDC. An electron trajectory code uses the resultant data to perform the actual collector analysis. The MDC was optimized for maximum MDC efficiency and minimum final kinetic energy of all collected electrons (to reduce heat transfer). The preceding figure shows the geometric and electrical configuration of an optimized collector with an efficiency of 93.8 percent. The results show the improvement in collector efficiency from 89.7 to 93.8 percent, resulting in an increase of three overall efficiency points. In addition, the time to design a highly efficient MDC was reduced from a month to a few days. All work was done in-house at Glenn for the High Rate Data Delivery Program. Future plans include optimizing the MDC and TWT interaction circuit in tandem to further improve overall TWT efficiency.

Vaden, Karl R.↗

Analysis of aerothermal loads on spherical dome protuberances

Hypersonic flow over spherical dome protuberances was investigated to determine increased pressure and heating loads to the surface. The configuration was mathematically modeled in a time-dependent three-dimensional analysis of the conservation of mass, momentum (Navier-Stokes), and energy equations. A boundary mapping technique was used to obtain a rectangular parallelepiped computational domain, a MacCormack explicit time-split predictor-corrector finite difference algorithm was used to obtain solutions. Results show local pressures and heating rates for domes one-half, one, and two boundary layer thicknesses high were increased by factors on the order of 1.4, 2, and 6, respectively. Flow over the lower dome was everywhere attached while flow over the intermediate dome had small windward and leeside separations. The higher dome had an unsteady windward separation region and a large leeside separation region. Trailing vortices form on all domes with intensity increasing with dome height. Discussion of applying the results to a thermally bowed thermal protection system are presented.

Olsen, G. C.↗

Analysis of aerothermal loads on spherical dome protuberances

Hypersonic flow over spherical dome protuberances was investigated to determine increased pressure and heating loads to the surface. The configuration was mathematically modeled in a time-dependant three-dimensional analysis of the conservation of mass, momentum (Navier-Stokes), and energy equations. A boundary mapping techique was used to obtain a rectangular parallelepiped computational domain, and a MacCormack explicit time-split predictor-corrector finite difference algorithm was used to obtain solutions. Results show local pressures and heating rates for domes one-half, one, and two boundary layer thicknesses high were increased by factors on the order of 1.4, 2, and 6, respectively. However, because lee-side pressure and thermal loads were reduced the two lower height domes did not experience any net increase in total loads. The total loads on the higher dome were increased by twenty-five percent. Flow over the lower dome was everywhere attached while flow over the intermediate dome had small windward and leeside separations. The higher dome had an unsteady windward separation region and a large leeside separation region. Trailing vortices form on all domes with intensity increasing with dome height. Discussions of applying the results to a thermally bowed thermal protection system are presented.

Olsen, G. C.↗

A radiometric model of an earth radiation budget radiometer optical system with diffuse-specular surfaces

The Earth Radiation Budget Experiment (ERBE) is to fly on NASA's Earth Radiation Budget Satellite (ERBS) and on NOAA F and NOAA G. Large spatial scale earth energy budget data will be derived primarily from measurements made by the ERBE nonscanning instrument (ERBE-NS). A description is given of a mathematical model capable of simulating the radiometric response of any of the ERBE-NS earth viewing channels. The model uses a Monte Carlo method to accurately account for directional distributions of emission and reflection from optical surfaces which are neither strictly diffuse nor strictly specular. The model computes radiation exchange factors among optical system components, and determines the distribution in the optical system of energy from an outside source. Attention is also given to an approach for implementing the model and results obtained from the implementation.

Luther, M. R.↗

Notes on aerodynamic forces 1 : rectilinear motion

The study of the motion of perfect fluids is of paramount importance for the understanding of the chief phenomena occurring in the air surrounding an aircraft, and for the numerical determination of their effects. The author recently successfully employed some simple methods for the investigation of the flow of a perfect fluid that have never been mentioned in connection with aeronautical problems. These methods appeal particularly to the engineer who is untrained in performing laborious mathematical computations, as they do away with these and allow one to obtain many interesting results by the mere application of some general and well-known principles of mechanics. Discussed here are the kinetic energy of moving fluids, the momentum of a body in a perfect fluid, two dimensional flow, three dimensional flow, and the distribution of the transverse forces of very elongated surfaces of revolution.

Munk, Max M↗

A Six-Node Curved Triangular Element and a Four-Node Quadrilateral Element for Analysis of Laminated Composite Aerospace Structures

Mathematical development and some computed results are presented for Mindlin plate and shell elements, suitable for analysis of laminated composite and sandwich structures. These elements use the conventional 3 (plate) or 5 (shell) nodal degrees of freedom, have no communicable mechanisms, have no spurious shear energy (no shear locking), have no spurious membrane energy (no membrane locking) and do not require arbitrary reduction of out-of-plane shear moduli or under-integration. Artificial out-of-plane rotational stiffnesses are added at the element level to avoid convergence problems or singularity due to flat spots in shells. This report discusses a 6-node curved triangular element and a 4-node quadrilateral element. Findings show that in regular rectangular meshes, the Martin-Breiner 6-node triangular curved shell (MB6) is approximately equivalent to the conventional 8-node quadrilateral with integration. The 4-node quadrilateral (MB4) has very good accuracy for a 4-node element, and may be preferred in vibration analysis because of narrower bandwidth. The mathematical developments used in these elements, those discussed in the seven appendices, have been applied to elements with 3, 4, 6, and 10 nodes and can be applied to other nodal configurations.

Martin, C. Wayne↗

Liquefaction dynamics of n-octadecane in cylindrical coordinates

The study presents one method of mathematically modeling the effects of gravity-induced free convection upon the melting phenomena of a finite paraffin slab. The theoretical model is based on the numerical computer solution of the solid phase and liquid phase energy equations, coupled with an imposed ideal-viscous flow approximation. An implicit alternating-direction technique is used as the finite-difference approximation for the basic energy equations in the numerical solution. The theoretical predictions are compared to experimental data from an annular test cell, in which n-octadecane, C18H38, was used as the test material. The results show that the model is a good first approximation to the solution of solid-liquid phase change involving gravity-induced free convection.

Bain, R. L.↗

Metabolic Adaptation to Muscle Ischemia

Although all tissues in the body can adapt to varying physiological/pathological conditions, muscle is the most adaptable. To understand the significance of cellular events and their role in controlling metabolic adaptations in complex physiological systems, it is necessary to link cellular and system levels by means of mechanistic computational models. The main objective of this work is to improve understanding of the regulation of energy metabolism during skeletal/cardiac muscle ischemia by combining in vivo experiments and quantitative models of metabolism. Our main focus is to investigate factors affecting lactate metabolism (e.g., NADH/NAD) and the inter-regulation between carbohydrate and fatty acid metabolism during a reduction in regional blood flow. A mechanistic mathematical model of energy metabolism has been developed to link cellular metabolic processes and their control mechanisms to tissue (skeletal muscle) and organ (heart) physiological responses. We applied this model to simulate the relationship between tissue oxygenation, redox state, and lactate metabolism in skeletal muscle. The model was validated using human data from published occlusion studies. Currently, we are investigating the difference in the responses to sudden vs. gradual onset ischemia in swine by combining in vivo experimental studies with computational models of myocardial energy metabolism during normal and ischemic conditions.

Cabrera, Marco E.↗

Mathematical Model of a Regenerative Fuel Cell for System Optimization

This thesis developed a system-level optimization model of a regenerative fuel cell (RFC) system for long-duration, off-world energy storage applications. Prior RFC design studies have typically been limited to reduced parameter sets and simplified constraints due to computational limitations relative to the number of relevant degrees of freedom. As a result, important nonlinear interactions between subsystems have not been fully captured. This work began to address that gap by developing a higher-fidelity, nonlinear optimization framework that incorporates a broader set of design variables and coupled constraints, enabling a multidimensional model that captures the coupled behavior of RFC subsystems and demonstrates the feasibility of applying optimization to such systems. An expanded system-level optimization approach was established that captures interactions between electrochemical performance, structural requirements, and storage design. This enabled a more comprehensive evaluation of trade-offs than conventional formulations. The model integrates four coupled subsystems: a fuel cell, an electrolyzer, reactant gas, and high-pressure storage tanks, and was formulated to accommodate a wide range of mission parameters, including operational time and required output power. It incorporates constraints on available solar array power, reactant mass balance between production and consumption, and pressure-dependent storage requirements. To enable reliable convergence, the optimization problem was reformulated to reduce dimensionality and improve numerical stability, with subsystem models organized for efficient evaluation. Problem dimensionality was reduced by consolidating lower-level design variables into higher-level representative quantities, and subsystem behavior was evaluated within the optimization loop. A multi-start initialization strategy was employed to mitigate sensitivity to local minima and improve solution quality, while nonlinear relationships were solved using robust numerical methods. The results showed that convergence was achieved across a range of required output power values. Specific energy reached a maximum at a critical mission power level, where the electrolyzer power matched the available solar input and operated near its voltage and current density limits. Beyond this point, further increases in required power resulted in less mass-efficient operation, increasing total system mass and reducing overall performance. The developed model represents an advancement in RFC system-level optimization by enabling analysis of a broader and more tightly coupled design space than previous considerations. While convergence behavior and computational cost remain challenges, the methods introduced improve solvability and allow inclusion of additional design variables with minimal loss of physical fidelity. However, the numerical results should not be interpreted as definitive design recommendations, as the model includes simplifying assumptions and omits several higher-order effects. Future work should extend this framework by incorporating additional subsystems and loss mechanisms, such as thermal management, parasitic power consumption, and reactant losses, to improve fidelity and ensure more representative design conclusions.

Electrochemistry↗

Crash simulation of composite and aluminum helicopter fuselages using a finite-element program

This paper describes a mathematical investigation of the crash-impact responses of an all-composite helicopter cockpit section incorporating an energy absorbing concept and one of conventional aluminum construction, using the Grumman DYCAST finite-element nonlinear structural dynamics computer program. The overall objective of this initial investigation was to explore the application of DYCAST to the design of future helicopter structures. The specific purposes are to assess the value of the finite-element code as a crashworthiness design analysis tool, and to compare the responses of the two fuselage sections. The results indicate that the finite-element simulation is a powerful analysis tool, providing sufficient detail to evaluate individual structural components, and that the crashworthiness of an all-composite cockpit section could be enhanced with energy-absorbing concepts

Winter, R.↗

The small community solar thermal power experiment. Parabolic dish technology for industrial process heat application

Aspects of incorporating a thermal energy transport system (ETS) into a field of parabolic dish collectors for industrial process heat (IPH) applications were investigated. Specific objectives are to: (1) verify the mathematical optimization of pipe diameters and insulation thicknesses calculated by a computer code; (2) verify the cost model for pipe network costs using conventional pipe network construction; (3) develop a design and the associated production costs for incorporating risers and downcomers on a low cost concentrator (LCC); (4) investigate the cost reduction of using unconventional pipe construction technology. The pipe network design and costs for a particular IPH application, specifically solar thermally enhanced oil recovery (STEOR) are analyzed. The application involves the hybrid operation of a solar powered steam generator in conjunction with a steam generator using fossil fuels to generate STEOR steam for wells. It is concluded that the STEOR application provides a baseline pipe network geometry used for optimization studies of pipe diameter and insulation thickness, and for development of comparative cost data, and operating parameters for the design of riser/downcomer modifications to the low cost concentrator.

Polzien, R. E.↗

Machine-Learned Committor Functions for Reactive Molecular Dynamics

Reactive molecular dynamics (MD) is a powerful tool for atomistic-scale modeling of a diverse range of chemical processes. However, scaling these simulations to large systems and long times scales remains a challenge because of the complexity of the potential energy function required. The authors previously developed a heuristic approach, called REACTER, that incorporates reactivity in MD simulations in a less general but much more computationally efficient manner. REACTER uses standard, fixed valence force fields as the underlying potentialenergy surface for describing all interatomic interactions but adds a procedure for enforcing user-defined reactions that occur when certain geometric constraints on relative atomic positions are satisfied. Further, these bonding changes can be accepted or rejected with a probability related tothe local thermal energy. This work seeks to generalize this approach by replacing the set of user defined geometric constraints and energetic criteria with a committor function that specifies the probability of a reaction occurring on the basis of the local atomic configuration. The committor function is a useful mathematical tool for modeling rare events but, unfortunately, is very difficult to compute for realistic systems in a general way. This work describes a method for approximating the committor function using a machine learning approach, specifically a deep neural network trained with data from reactive MD and DFT-based dynamics simulations. This network is coupled to the existing REACTER protocol, as implemented in the LAMMPS MD package, and used to make on-the-fly predictions of reaction probabilities without the more extensive user input previously required. The new method is demonstrated using the polymerization of polystyrene as a case study. Although very dependent on the quality and quantity of training data, machine-learned committor functions show promise as a method for incorporating reaction probability from higher level calculations into highly scalable MD simulations.

polymer simulations↗

Synthetic battery cycling

The trend in energy storage is toward systems with high-voltage and high-power outputs. These conditions accentuate the cumulative effects of small differences in characteristics from cell to cell in any large multicell battery. These considerations are of particular concern with a number of emerging electrochemical concepts where convenient overcharge reactions do not exist. 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 will be 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.↗

Calculations enable optimum design of magnetic brake

Mathematical analysis and computations determine optimum magnetic coil configurations for a magnetic brake which controllably decelerates a free falling load to a soft stop. Calculations on unconventionally wound coils determine the required parameters for the desired deceleration with minimum electrical energy supplied to the stationary coil.

Kosmahl, H. G.↗

Phosphoric acid fuel cell power plant system performance model and computer program

A FORTRAN computer program was developed for analyzing the performance of phosphoric acid fuel cell power plant systems. Energy mass and electrochemical analysis in the reformer, the shaft converters, the heat exchangers, and the fuel cell stack were combined to develop a mathematical model for the power plant for both atmospheric and pressurized conditions, and for several commercial fuels.

Alkasab, K. A.↗

Chaos: Understanding and Controlling Laser Instability

In order to characterize the behavior of tunable diode lasers (TDL), the first step in the project involved the redesign of the TDL system here at the University of Tennessee Molecular Systems Laboratory (UTMSL). Having made these changes it was next necessary to optimize the new optical system. This involved the fine adjustments to the optical components, particularly in the monochromator, to minimize the aberrations of coma and astigmatism and to assure that the energy from the beam is focused properly on the detector element. The next step involved the taking of preliminary data. We were then ready for the analysis of the preliminary data. This required the development of computer programs that use mathematical techniques to look for signatures of chaos. Commercial programs were also employed. We discovered some indication of high dimensional chaos, but were hampered by the low sample rate of 200 KSPS (kilosamples/sec) and even more by our sample size of 1024 (1K) data points. These limitations were expected and we added a high speed data acquisition board. We incorporated into the system a computer with a 40 MSPS (million samples/sec) data acquisition board. This board can also capture 64K of data points so that were then able to perform the more accurate tests for chaos. The results were dramatic and compelling, we had demonstrated that the lead salt diode laser had a chaotic frequency output. Having identified the chaotic character in our TDL data, we proceeded to stage two as outlined in our original proposal. This required the use of an Occasional Proportional Feedback (OPF) controller to facilitate the control and stabilization of the TDL system output. The controller was designed and fabricated at GSFC and debugged in our laboratories. After some trial and error efforts, we achieved chaos control of the frequency emissions of the laser. The two publications appended to this introduction detail the entire project and its results.

Blass, William E.↗