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.
A two-dimensional finite difference numerical technique is presented to determine the temperature distribution in a solid blade of a radial guide vane. A computer program is written in Fortran IV for IBM 370/165 computer. The computer results obtained from these programs have a similar behavior and trend as those obtained by experimental results.
A two dimensional finite difference numerical technique is presented to determine the temperature distribution of an internal cooled blade of radial turbine guide vanes. A simple convection cooling is assumed inside the guide vane blade. Such cooling has relatively small cooling effectiveness at the leading edge and at the trailing edge. Heat transfer augmentation in these critical areas may be achieved by using impingement jets and film cooling. A computer program is written in FORTRAN IV for IBM 370/165 computer.
Approximate solution for large asteroid distribution with masses near limiting largest mass of population
An exact solution of the problem of designing an airfoil with a prescribed velocity distribution on the suction surface in a given uniform flow of an incompressible perfect fluid is obtained by replacing the boundary of the airfoil by vortices. By this device, a method of solution is developed that is applicable both to isolated airfoils and to airfoils in cascade. The conformal transformation of the designed airfoil into a circle can then be obtained and the velocity distribution at any angle of attack computed. Numerical illustrations of the method are given for the airfoil in cascade.
We introduce NeuroSA, a neuromorphic architecture specifically designed to ensure asymptotic convergence to the ground state of an Ising problem using a Fowler-Nordheim quantum mechanical tunneling based threshold-annealing process. The core component of NeuroSA consists of a pair of asynchronous ON-OFF neurons, which effectively map classical simulated annealing dynamics onto a network of integrate-and-fire neurons. The threshold of each ON-OFF neuron pair is adaptively adjusted by an FN annealer and the resulting spiking dynamics replicates the optimal escape mechanism and convergence of SA, particularly at low-temperatures. To validate the effectiveness of our neuromorphic Ising machine, we systematically solved benchmark combinatorial optimization problems such as MAX-CUT and Max Independent Set. Across multiple runs, NeuroSA consistently generates distribution of solutions that are concentrated around the state-of-the-art results (within 99%) or surpass the current state-of-the-art solutions for Max Independent Set benchmarks. Furthermore, NeuroSA is able to achieve these superior distributions without any graph-specific hyperparameter tuning. For practical illustration, we present results from an implementation of NeuroSA on the SpiNNaker2 platform, highlighting the feasibility of mapping our proposed architecture onto a standard neuromorphic accelerator platform.
A numerical lifting surface formulation, including computed results for planar wing cases is presented. This formulation, referred to as the vortex spline scheme, combines the adaptability to complex shapes offered by paneling schemes with the smoothness and accuracy of loading function methods. The formulation employes a continuous distribution of singularity strength over a set of panels on a paneled wing. The basic distributions are independent, and each satisfied all the continuity conditions required of the final solution. These distributions are overlapped both spanwise and chordwise. Boundary conditions are satisfied in a least square error sense over the surface using a finite summing technique to approximate the integral. The current formulation uses the elementary horseshoe vortex as the basic singularity and is therefore restricted to linearized potential flow. As part of the study, a non planar development was considered, but the numerical evaluation of the lifting surface concept was restricted to planar configurations. Also, a second order sideslip analysis based on an asymptotic expansion was investigated using the singularity spline formulation.
A new numerical formulation with computed results, is presented. This formulation combines the adaptability to complex shapes offered by paneling schemes with the smoothness and accuracy of the loading function methods. The formulation employs a continuous distribution of singularity strength over a set of panels on a paneled wing. The basic distributions are independent, and each satisfies all of the continuity conditions required of the final solution. These distributions are overlapped both spanwise and chordwise (termed 'spline'). Boundary conditions are satisfied in a least square error sense over the surface using a finite summing technique to approximate the integral.
The simplified shear solution method is presented for approximating the through-thickness shear stress distribution within a composite laminate or panel based on laminated beam theory. The method does not consider the solution of a particular boundary value problem; rather it requires only knowledge of the global shear loading, geometry, and material properties of the laminate or panel. It is thus analogous to lamination theory in that ply level stresses can be efficiently determined from global load resultants (as determined, for instance, by finite element analysis) at a given location in a structure and used to evaluate the margin of safety on a ply by ply basis. The simplified shear solution stress distribution is zero at free surfaces, continuous at ply boundaries, and integrates to the applied shear load. Comparisons to existing theories are made for a variety of laminates, and design examples are provided illustrating the use of the method for determining through-thickness shear stress margins in several types of composite panels and in the context of a finite element structural analysis.
Rising penetration of consumer-owned Distribution Grid Resources (DGRs), increasingly managed by third party aggregators and enrolled in grid services and wholesale market programs, can create localized congestion in distribution networks. Managing these constraints is challenging due to a persistent coordination and information gap: utilities are accountable for reliability and have network topology and state visibility, while aggregators control the DGR capability needed to relieve congestion. This survey synthesizes congestion management solutions for distribution systems with high DGR penetration, covering both market-based mechanisms (distribution level markets, locational pricing, flexibility auctions) and non-market-based solutions (network reconfiguration, direct DGR control, demand response, curtailment, etc.). The literature is organized across three decision horizons: long term planning, operational planning, and real-time operation. Special attention is devoted to emerging distribution system operator architectures and coordination frameworks spanning transmission system operators, aggregators, and microgrids. Drawing on recent case studies and implementations, we distill best practices, identify key technical and economic barriers, and outline research directions. The evidence points to a shift toward integrated congestion management that combines market signals with technical controls, enabled by improved monitoring, forecasting, and closed loop control capabilities.
When electronic oscillators are coupled to nearest neighbors to form an array on a hexagonal lattice, the planar phase distributions desired for excitation of a phased array antenna are not steady state solutions of the governing non-linear equations describing the system. Thus the steady state phase distribution deviates from planar. It is shown to be possible to obtain an exact solution for the steady state phase distribution and thus determine the deviation from the desired planar distribution as a function of beam steering angle.
The Least-cost Optimal Distribution Grid Expansion (LODGE) model provides the optimal portfolio of distribution system upgrades—e.g., voltage regulators, feeder reconductoring, transformer upgrades and non-wires alternatives (NWA), such as strategic siting of storage and distributed generation—to interconnect distributed energy resources (DERs) and enable load growth. It can be used to assess grid infrastructure costs and explore policy and regulatory solutions for distribution planning and DER valuation.In 2025, Berkeley Lab conducted three pilot analyses to validate LODGE results with empirical utility data before the model’s first release in 2026. The pilots, done with utilities in Washington, Colorado, and New Mexico, provide examples that illustrate how the model works, what it can do, and the value of the analysis.
As computing technology continues to advance, computational modeling of scientific and engineering problems produces data of increasing complexity: large in size and unstructured in shape. Volume visualization of such data is a challenging problem. This paper proposes a distributed parallel solution that makes ray-casting volume rendering of unstructured-grid data practical. Both the data and the rendering process are distributed among processors. At each processor, ray-casting of local data is performed independent of the other processors. The global image composing processes, which require inter-processor communication, are overlapped with the local ray-casting processes to achieve maximum parallel efficiency. This algorithm differs from previous ones in four ways: it is completely distributed, less view-dependent, reasonably scalable, and flexible. Without using dynamic load balancing, test results on the Intel Paragon using from two to 128 processors show, on average, about 60% parallel efficiency.
Analysis of phase diagrams and cation distributions within crystalline solutions as means of obtaining thermodynamic data on rock forming crystalline solutions is discussed along with some aspects of partitioning of elements in coexisting phases. Crystalline solutions, components in a silicate mineral, and chemical potentials of these components were defined. Examples were given for calculating thermodynamic mixing functions in the CaW04-SrW04, olivine-chloride solution, and orthopyroxene systems.
A viscous-inviscid interaction algorithm is developed for prediction of two-dimensional mean and fluctuating velocity distributions in the wake immediately downstream of an airfoil trailing edge. A composite pressure field is defined, and a Poisson equation solved for transverse pressure variations. A parabolized form of the time-averaged steady Navier-Stokes equations are solved in conjunction with a viscous-augmented two-dimensional inviscid potential flow analysis. A tensor constitutive equation is employed to predict Reynolds stress distributions from solutions of a turbulence kinetic energy two equation closure model. Numerical predictions compared favorably with detailed experimental data for mean and fluctuating velocities, and Reynolds shear stress distributions, in the trailing edge region of a NACA 63-012 airfoil.
The described investigation is concerned with the solution of the non-LTE optically thick transfer equations for hydrogen, carbon, and other constituents to determine semiempirical models for six components of the quiet solar chromosphere. For a given temperature-height distribution, the solution is obtained of the equations of statistical equilibrium, radiative transfer for lines and continua, and hydrostatic equilibrium to find the ionization and excitation conditions for each atomic constituent. The emergent spectrum is calculated, and a trial and error approach is used to adjust the temperature distribution so that the emergent spectrum is in best agreement with the observed one. The relationship between semiempirical models determined in this way and theoretical models based on radiative equilibrium is discussed by Avrett (1977). Harvard Skylab EUV observations are used to determine models for a number of quiet-sun regions.
Experiments were conducted in air, CO2, and CO2-argon gas mixtures to obtain heating distribution data over a 140 deg blunt cone with various corner radii. The effect of corner radius on the heating distribution over the forebody of the cone was included in the investigation. These experiments provide data for validation of two-dimensional axisymmetric and three-dimensional Navier-Stokes solutions. Heating distribution data and measured bow shock wave stand-off distances for 0 deg angle of attack were compared with predicted values using a two-dimensional axisymmetric Navier-Stokes code.
Mixing of two fluids generated by steady and particularly g-jitter acceleration is fundamental towards the understanding of transport phenomena in a microgravity environment. We propose to carry out flight and ground-based experiments to quantify flow fields due to g-jitter type of accelerations using Stereo Imaging Velocimetry (SIV), and measure the concentration field using laser fluorescence. The understanding of the effects of g-jitter on transport phenomena is of great practical interest to the microgravity community and impacts the design of experiments for the Space Shuttle as well as the International Space Station. The aim of our proposed research is to provide quantitative data to the community on the effects of g-jitter on flow fields due to mixing induced by buoyancy forces. The fundamental phenomenon of mixing occurs in a broad range of materials processing encompassing the growth of opto-electronic materials and semiconductors, (by directional freezing and physical vapor transport), to solution and protein crystal growth. In materials processing of these systems, crystal homogeneity, which is affected by the solutal field distribution, is one of the major issues. The understanding of fluid mixing driven by buoyancy forces, besides its importance as a topic in fundamental science, can contribute towards the understanding of how solutal fields behave under various body forces. The body forces of interest are steady acceleration and g-jitter acceleration as in a Space Shuttle environment or the International Space Station. Since control of the body force is important, the flight experiment will be carried out on a tunable microgravity vibration isolation mount, which will permit us to precisely input the desired forcing function to simulate a range of body forces. To that end, we propose to design a flight experiment that can only be carried out under microgravity conditions to fully exploit the effects of various body forces on fluid mixing. Recent flight experiments, by the P.I. through collaboration with the Canadian Space Agency (STS-85, August 1997), aimed at determining the stability of the interface between two miscible liquids inside an enclosure show that a long liquid column (5 cm) under microgravity isolation conditions can be stable, i.e. the interface remains sharp and vertical over a short time scale; thus transport occurs by molecular mass diffusion. On the other hand, when the two liquids were excited from a controlled vibration source (Microgravity Vibration Isolation Mount) two to four mode large amplitude quasi-stationary waves were observed. The data was limited to CCD recording of the dynamics of the interface between the two fluids. We propose to carry out flight experiments to quantify the dynamics of the flow field using Stereo Imaging Velocimetry and measure the concentration field using laser fluorescence. The results will serve as a basis to understand effects of g-jitter on transport phenomena, in this case mass diffusion. As the measurement of the kinematics of the flow field will shed light on the instability mechanism. The research will allow measurement of the flow field in microgravity environment to prove two hypotheses: (1) Maxwell's hypothesis: finite convection always exists in diffusing systems, and (2) Quasi-stationary waves inside a bounded enclosure in a microgravity environment is generated by Kelvin-Helmholtz instability; resonance of the interface which produces incipient mixing is due to Rayleigh-Taylor instability. The first hypothesis can be used as a benchmark experiment to illustrate diffusive mixing. The second hypothesis will lead to the understanding of g-jitter effects on buoyancy driven flow fields which occur in many situations involving materials processing, and other basic fluid physics phenomena. In addition, the second hypothesis will also provide insight in how Rayleigh-Taylor and Kelvin-Helmholtz instabilities propagate concentration fronts during mixing. Measurement of the flow field using SIV is important because it is the flow field which causes instability at the interface between the two fluids. Mixing driven by buoyancy induced flow fields will be addressed both experimentally and computationally. The experimental effort will address the kinematics of mixing: stretching, transport and chaos. Quantification of the mechanisms of mixing will consists of measuring the flow field using the SIV system at Glenn and capturing the dynamics of the interface, to measure mass transport, using a CCD camera. These experiments will be carried out within the framework of Earth's gravity and g-jitter microgravity acceleration as in a Space Shuttle environment or the International Space Station. The g-jitter will be induced and controlled using a tunable vibration isolation platform to isolate against vibration as well as input periodic and random vibration to the system. The parametric range of the microgravity experiment will be extended from the experiments on STS-85 to investigate higher mode quasi-stationary waves (8 to 12), as well as resonance regions which leads to chaos and turbulence. Ground-based experiments will focus on effects of vibration on stably stratified fluid layers in order to scale for possible scenarios in a microgravity environment. These vibrations will be subjected perpendicular to the concentration field on the ground since the parallel case can only be carried out in a microgravity environment. The concept of dynamical similarity will be applied to tune the experiments as closely as possible to a Space Shuttle environment or the International Space Station. The computational effort will take advantage of the Computational Laboratory at Glenn to corroborate the experimental findings with predictions of the dynamics of the flow field using the codes FLUENT (finite difference based) and FIDAP (finite element based). We will investigate two important cases, single-fluid model to address dilute systems with negligible jump in viscosity and the more general two-fluid model which accounts for finite jump in viscosity. Apart from its microgravity relevance, this experiment is well suited to study dynamics in nonlinear systems.
The fluid flow and solute segregation which occur during directional solidification, including effects of gravity and microgravity were studied. The nature of the fluid flow, its effects on the shape of the crystal-metal interface and the resulting distribution of solutes are examined. During solidification of a binary alloy at constant velocity vertically upwards, thermosolutal convection occurs if the solute rejection or preferential incorporation at the crystal-melt interface decreases the density of the melt. Numerical calculations of the solute, temperature, and flow fields are carried out by using previously developed algorithms for Schmidt numbers Sc = 1, 10, 81. The time periodic flow, which occurs over a narrow range of parameters, is investigated.