SEARCH · Engineering Papers

Results for “multiple 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 811 records · Page 45

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Multiple relaxation and inhomogeneous broadening in resonance enhanced Raman scattering - Application to tunable infrared generation

The solutions for the imaginary susceptibility of the Raman field transition with arbitrary relaxation rates and field strengths are examined for differing sets of relaxation rates with emphasis on alkali metal vapors which have spontaneous emission dominated relaxation. The model is further expanded to include Doppler broadening and used to predict the peak gain as a function of detuning for a frequency doubled alexandrite laser-pumped cesium vapor gain cell.

Ryan, J. C.↗

NASA NTRS2017

Parametric Evaluation of Interstellar Exploration Mission Concepts

One persistent difficulty in evaluating the myriad advanced propulsion concepts proposed over the last 60 years is a true apples to apples comparison of the expected gain in performance. This analysis is complicated by numerous factors including, multiple missions of interest to the advanced propulsion community, the lack of a credible closed form solution to 'medium thrust' trajectories, and lack of detailed design data for most proposed concepts that lend credibility to engine performance estimates. This paper describes a process on how to make fair comparisons of different propulsion concepts for multiple missions over a wide range of performance values. The figure below illustrates this process. This paper describes in detail the process and outlines the status so far in compiling the required data. Parametric data for several missions are calculated and plotted against specific power-specific impulse scatter plots of expected propulsion system performance. The overlay between required performance as defined by the trajectory parametrics and expected performance as defined in the literature for major categories of propulsion systems clearly defines which propulsion systems are the most apt for a given mission. The application of the Buckingham Pi theorem to general parameters for interstellar exploration ( mission time, mass, specific impulse, specific power, distance, propulsion source energy/mass, etc.) yields a number of dimensionless variables. The relationships of these variables can then be explored before application to a particular mission. Like in the fields of fluid mechanics and heat transfer, the use of the Buckingham Pi theorem results in new variables to make apples to apples comparisons.

Adams, Robert B.↗

NASA NTRS1976

Elliptic integral solutions to a class of space flight optimization problems

This paper is initially concerned with the minimum-time, exoatmospheric flight of a rocket with constant thrust acceleration magnitude, as in the cases of nuclear and solar electric propulsion. Gravitational acceleration is assumed to be a constant scalar multiple of the radius vector, plus a correction term which is a given function of time. The solution to the state equations is obtained in terms of elliptic integrals. A method is presented for the solution of the two-point boundary-condition problem associated with orbital transfer. At most, the latter method requires iteration upon final time, angle of injection, and two other parameters which are bounded. An example problem is provided which involves a rocket with very low thrust and a spiraling trajectory of many revolutions, but an altitude change of only several hundred miles above the earth. Finally, the original elliptic integral solution is extended to a larger class of low and intermediate thrust problems with constant thrust magnitude, mass decreasing with time, and an inverse square gravitational force.

Andrus, J. F.↗

NASA NTRS1978

Aerodynamic preliminary analysis system. Part 2: User's manual and program description

A comprehensive aerodynamic analysis program based on linearized potential theory is described. The solution treats thickness and attitude problems at subsonic and supersonic speeds. Three dimensional configurations with or without jet flaps having multiple nonplanar surfaces of arbitrary planform and open or closed slender bodies or noncircular contour are analyzed. Longitudinal and lateral-directional static and rotary derivative solutions are generated. The analysis is implemented on a time sharing system in conjunction with an input tablet digitizer and an interactive graphics input/output display and editing terminal to maximize its responsiveness to the preliminary analysis problem. Nominal case computation time of 45 CPU seconds on the CDC 175 for a 200 panel simulation indicates the program provides an efficient analysis for systematically performing various aerodynamic configuration tradeoff and evaluation studies.

Divan, P.↗

NASA NTRS1987

A generalized procedure for constructing an upwind-based TVD scheme

A generalized formulation for constructing second- and higher-order accurate TVD (total variation diminishing) schemes is presented. A given scheme is made TVD by limiting antidiffusive flux differences with some nonlinear functions, so-called limiters. The general idea of the formulation and its mathematical proof of Harten's TVD conditions is shown by applying the Lax-Wendroff method to a scalar nonlinear equation and constant-coefficient system of conservation laws. For the system of equations, several definitions are derived for the argument used in the limiter function and present their performance to numerical experiments. Then the formulation is formally extended to the nonlinear system of equations. It is demonstrated that use of the present procedure allows easy conversion of existing central or upwind, and second- or higher-order differencing schemes so as to preserve monotonicity and to yield physically admissible solutions. The formulation is simple mathematically as well as numerically; neither matrix-vector multiplication nor Riemann solver is required. Roughly twice as much computational effort is needed as compared to conventional scheme. Although the notion of TVD is based on the initial value problem, application to the steady Euler equations of the formulation is also made. Numerical examples including various ranges of problems show both time- and spatial-accuracy in comparison with exact solutions.

Liou, Meng-Sing↗

NASA NTRS1984

Flow in a differentially rotated cylindrical drop at moderate Reynolds number

Galerkin finite-element approximations are combined with computer-implemented perturbation methods for tracking families of solutions to calculate the steady axisymmetric flows in a differentially rotated cylindrical drop as a function of Reynolds number Re, drop aspect ratio and the rotation ratio between the two end disks. The flows for Reynolds numbers below 100 are primarily viscous and reasonably described by an asymptotic analysis. When the disks are exactly counter-rotated, multiple steady flows are calculated that bifurcate to higher values of Re from the expected solution with two identical secondary cells stacked symmetrically about the axial midplane. The new flows have two cells of different size and are stable beyond the critical value Re sub c. The slope of the locus of Re sub c for drops with aspect ratio up to 3 disagrees with the result for two disks of infinite radius computed assuming the similarity form of the velocity field. Changing the rotation ratio for exact counter-rotation ruptures the junction of the multiple flow fields into two separated flow families.

Harriott, G. M.↗

DOE OSTI2025

Resilient Operation of Networked Community Microgrids with High Solar Penetration

This project, funded by the US Department of Energy’s Solar Energy Technologies Office (SETO), focused on the operation of microgrids as a coordinated network. The primary objective, which was successfully achieved, was to develop both control strategies and hardware solutions to support the resilient and efficient operation of networked microgrids with high solar penetration. The work was structured around the following four main tasks: • Development of distributed and scalable optimization algorithms for AC-coupled networked microgrids. • Design and implementation of a novel DC interconnection hardware to enable precise power exchange between microgrids. • Laboratory operational validation of the developed technologies using 480 V testbeds and commercially available hardware. • Field operational validation of the complete solution in Adjuntas, Puerto Rico, interconnecting two kW-scale, split-phase microgrids of Casa Pueblo’s microgrids. This project addressed multiple technical challenges across the domains of optimization, control, hardware interconnection, and protection. One of its key contributions was delivering tangible, real-world solutions for networking microgrids. In contrast to purely theoretical or simulation-based work, this project included full-scale hardware operational validation both in the lab and in the field. The work conducted as part of this project—in collaboration with the University of Puerto Rico; the University of Tennessee, Knoxville; the University of Central Florida; and Casa Pueblo—has advanced the state of the art in networked microgrids. Key contributions include the development of distributed control strategies, practical solutions for real-world implementation challenges, and the introduction of a novel DC interlink approach for microgrid interconnection. The project featured both laboratory and field validation using commercial off-the-shelf components. The field deployment successfully validated that a group of microgrids can operate in a coordinated manner, enabling precise power flow between systems and mutual support during extreme events. This project resulted in 15 journal publications and 15 conference papers; 5 graduate students and 15 undergraduate students were supported. The codes of distributed optimization and forecasting were made open-source through OSTI.gov for distributed optimization and forecasting. All the publications are available in the ORNL-hosted project landing page. The DC interlink with state-of-charge balancing control was operationally validated in Adjuntas by interconnecting two real-world, 240 V split-phase microgrids. To the best knowledge of the team, this represents the first operational validation of AC microgrids interconnected via DC-interlinks. As a culmination of this project, a follow-on grant was awarded to support the technology transfer of the distributed optimization framework to a commercial microgrid controller, Stellar Edge, developed by the California-based company New Sun Road.

14 SOLAR ENERGY↗

NASA NTRS2004

Recent Developments in OVERGRID, OVERFLOW-2 and Chimera Grid Tools Scripts

OVERGRID and OVERFLOW-2 feature easy to use multiple-body dynamics. The new features of OVERGRID include a preliminary chemistry interface, standard atmosphere and mass properties calculators, a simple unsteady solution viewer, and a debris tracking interface. Script library development in Chimera Grid Tools has applications in turbopump grid generation. This viewgraph presentation profiles multiple component dynamics, validation test cases for a sphere, cylinder, and oscillating airfoil, and debris analysis.

Chan, William M.↗

NASA NTRS1975

Use of numerically generated body-fitted coordinate systems for solution of the Navier-Stokes equations

A procedure for numerical solution of the time-dependent, two-dimensional incompressible Navier-Stokes equations that can treat the unsteady laminar flow about bodies of arbitrary shape, such as two-dimensional airfoils, multiple airfoils, and submerged hydrofoils, as naturally as it can deal with the flow about simple bodies. The solution is based on a method of automatic numerical generation of a general curvilinear coordinate system with coordinate lines coincident with all boundaries of a general multiconnected region containing any number of arbitrarily shaped bodies. The curvilinear coordinates are generated as the solution of two elliptical partial differential equations with Dirichlet boundary conditions, one coordinate being specified to be constant on each of the boundaries, and a distribution of the other being specified along the boundaries. The solution compares excellently with the Blasius boundary layer solution for the flow past a semiinfinite flat plate.

Thompson, J. F.↗

NASA NTRS2017

Low-Thrust Trajectory Optimization with Simplified SQP Algorithm

The problem of low-thrust trajectory optimization in highly perturbed dynamics is a stressing case for many optimization tools. Highly nonlinear dynamics and continuous thrust are each, separately, non-trivial problems in the field of optimal control, and when combined, the problem is even more difficult. This paper de-scribes a fast, robust method to design a trajectory in the CRTBP (circular restricted three body problem), beginning with no or very little knowledge of the system. The approach is inspired by the SQP (sequential quadratic programming) algorithm, in which a general nonlinear programming problem is solved via a sequence of quadratic problems. A few key simplifications make the algorithm presented fast and robust to initial guess: a quadratic cost function, neglecting the line search step when the solution is known to be far away, judicious use of end-point constraints, and mesh refinement on multiple shooting with fixed-step integration.In comparison to the traditional approach of plugging the problem into a “black-box” NLP solver, the methods shown converge even when given no knowledge of the solution at all. It was found that the only piece of information that the user needs to provide is a rough guess for the time of flight, as the transfer time guess will dictate which set of local solutions the algorithm could converge on. This robustness to initial guess is a compelling feature, as three-body orbit transfers are challenging to design with intuition alone. Of course, if a high-quality initial guess is available, the methods shown are still valid.We have shown that endpoints can be efficiently constrained to lie on 3-body repeating orbits, and that time of flight can be optimized as well. When optimizing the endpoints, we must make a trade between converging quickly on sub-optimal endpoints or converging more slowly on end-points that are arbitrarily close to optimal. It is easy for the mission design engineer to adjust this trade based on the problem at hand.The biggest limitation to the algorithm at this point is that multi-revolution transfers (greater than 2 revolutions) do not work nearly as well. This restriction comes in because the relationship between node 1 and node N becomes increasingly nonlinear as the angular distance grows. Trans-fers with more than about 1.5 complete revolutions generally require the line search to improve convergence. Future work includes: Comparison of this algorithm with other established tools; improvements to how multiple-revolution transfers are handled; parallelization of the Jacobian computation; in-creased efficiency for the line search; and optimization of many more trajectories between a variety of 3-body orbits.

Parrish, Nathan L.↗

NASA NTRS2021

Multi-Frequency NMM3D Simulations of Wave Propagation in Vegetation for Remote Sensing of Soil Moisture

To investigate the feasibility of using multi-frequencies to fur- ther improve the soil moisture retrieval, the recently devel- oped hybrid method is used to perform full-wave scattering simulations of a wheat field at three different frequencies of L-, S- and C- bands. In the hybrid method, the multiple scattering within a single wheat plant is first captured using the T-matrix of vec- tor cylindrical waves calculated from the full-wave solutions of HFSS. In the second step, the scattering among different plants is considered using the Foldy-Lax equations of multiple scattering theory (MST) and the extracted T-matrix. The transmission of microwaves through wheat field at the L-, S- and C-bands are calculated using the hybrid method. Results show that: (1) the transmission obtained from the full-wave simulations is much larger than those computed from the radiative transfer equations (RTE) model, (2) the hybrid method transmission has a weaker frequency dependence than those of the RTE, and (3) the attenuation caused by the vegetation layer would saturate with an increase in frequency.

Yueh, Simon↗

NASA NTRS2003

Development of Numerical Codes for Modeling Electromagnetic Behavior at High Frequencies Near Large Objects

A study into the problem of determining electromagnetic solutions at high frequencies for problems involving complex geometries, large sizes and multiple sources (e.g. antennas) has been initiated. Typical applications include the behavior of antennas (and radiators) installed on complex conducting structures (e.g. ships, aircrafts, etc..) with strong interactions between antennas, the radiation patterns, and electromagnetic signals is of great interest for electromagnetic compatibility control. This includes the overall performance evaluation and control of all on-board radiating systems, electromagnetic interference, and personnel radiation hazards. Electromagnetic computational capability exists at NASA LaRC, and many of the codes developed are based on the Moment Method (MM). However, the MM is computationally intensive, and this places a limit on the size of objects and structures that can be modeled. Here, two approaches are proposed: (i) a current-based hybrid scheme that combines the MM with Physical optics, and (ii) an Alternating Direction Implicit-Finite Difference Time Domain (ADI-FDTD) method. The essence of a hybrid technique is to split the overall scattering surface(s) into two regions: (a) a MM zone (MMZ) which can be used over any part of the given geometry, but is most essential over irregular and "non-smooth" geometries, and (b) a PO sub-region (POSR). Currents induced on the scattering and reflecting surfaces can then be computed in two ways depending on whether the region belonged to the MMZ or was part of the POSR. For the MMZ, the current calculations proceed in terms of basis functions with undetermined coefficients (as in the usual MM method), and the answer obtained by solving a system of linear equations. Over the POSR, conduction is obtained as a superposition of two contributions: (i) currents due to the incident magnetic field, and (ii) currents produced by the mutual induction from conduction within the MMZ. This effectively leads to a reduction in the size of linear equations from N to N - Npo with N being the total number of segments for the entire surface and Npo the number of segments over the POSR. The scheme would be appropriate for relatively large, flat surfaces, and at high frequencies. The ADI-FDTD scheme provides for both transient and steady state analyses. The restrictive Courant-Friedrich-Levy (CFL) condition on the time-step is removed, and so large time steps can be chosen even though the spatial grids are small. This report includes the problem definition, a detailed discussion of both the numerical techniques, and numerical implementations for simple surface geometries. Numerical solutions have been derived for a few simple situations.

Joshi, R. P.↗

DOE OSTI2026

Low-stress creep behavior of multiple salt caverns under cyclic operations

Here, the use of multiple salt caverns for large-scale underground hydrogen storage offers a strategic solution to the growing demand for efficient and flexible energy storage, due to the inherent low permeability, high mechanical strength, and self-healing properties of the salt formation. Despite these advantages, the long-term geomechanical response of salt caverns under cyclic injection and withdrawal remains a significant challenge for ensuring operational stability. This mechanistic study highlights the critical role of low-stress creep, also known as pressure solution creep, in governing deformation of salt caverns during cyclic operations. Unlike conventional creep mechanisms that dominate under high-stress conditions, low-stress creep becomes particularly relevant under the moderate stress changes induced by repeated injection and withdrawal, which provides a previously underexplored control on cavern closure. Our three-dimensional computational model integrates the low-stress creep mechanism into the Munson–Dawson creep equation across a spectrum of operational scenarios including variations in cavern pressure cycles, injection-withdrawal frequencies, and caprock mechanical property. The low-stress creep components accelerate the cavern creep closure, which is amplified with shorter cycle frequency and larger pressure difference. Sensitivity analyses capture the complex interactions between multiple caverns and surround formations at the field scale that allows for a quantitative assessment of deformation patterns and closure rates over time and thereby informs the design and operation of sustainable underground hydrogen storage system.

Cyclic operation↗

NASA NTRS2018

Improving Scientific Payload Communication Over Multiple, Low-Capacity Links

Communication is an important part of every airborne science mission. It provides a means for a large group of ground-based scientists to monitor and adjust experiments aboard aircraft that often cannot support a large number of passengers. In some cases, such as unmanned aircraft, communication becomes critical to a successful mission. Existing communication systems require the use of multiple, low-capacity channels, or links, to handle the capacity requirements of the scientists. Current communication solutions attempt to balance performance and reliability, but ultimately suffer from some critical design flaws that limit the ability of the system to handle additional links. We show how adding resources to the current system results in decreased performance and reliability. We explore the use of a new communication protocol, MultiPath Transmission Control Protocol, to make better use of multiple communication links and propose a total system capable of scaling to any number of channels efficiently and dynamically. The tools developed to support and measure this solution are also detailed.

MPTCP↗

NASA NTRS2009

Microfluidic Mixing Technology for a Universal Health Sensor

A highly efficient means of microfluidic mixing has been created for use with the rHEALTH sensor an elliptical mixer and passive curvilinear mixing patterns. The rHEALTH sensor provides rapid, handheld, complete blood count, cell differential counts, electrolyte measurements, and other lab tests based on a reusable, flow-based microfluidic platform. These geometries allow for cleaning in a reusable manner, and also allow for complete mixing of fluid streams. The microfluidic mixing is performed by flowing two streams of fluid into an elliptical or curvilinear design that allows the combination of the flows into one channel. The mixing is accomplished by either chaotic advection around micro - fluidic loops. All components of the microfluidic chip are flow-through, meaning that cleaning solution can be introduced into the chip to flush out cells, plasma proteins, and dye. Tests were performed on multiple chip geometries to show that cleaning is efficient in any flowthrough design. The conclusion from these experiments is that the chip can indeed be flushed out with microliter volumes of solution and biological samples are cleaned readily from the chip with minimal effort. The technology can be applied in real-time health monitoring at patient s bedside or in a doctor s office, and real-time clinical intervention in acute situations. It also can be used for daily measurement of hematocrit for patients on anticoagulant drugs, or to detect acute myocardial damage outside a hospital.

Chan, Eugene Y.↗

DOE OSTI2022

Development and performance of INCONEL® alloy 740H® seam-welded piping

INCONEL® alloy 740H® is an age-hardenable nickel-based alloy approved for pressure vessels and piping within ASME Section I, VIII, and B31.1. Currently, the code applies a weld strength reduction factor (WSRF) of 0.7 to the allowable stresses for longitudinal seam welds in the time dependent creep regime. In this work, a full-scale seam weld was successfully produced using typical industrial practices. The component was solution heat-treated and aged after fabrication to improve the performance of the weldment. Tensile, bend, impact, and low-cycle fatigue tests showed the component met all the specification requirements and elevated temperature properties were within base metal expectations with failures predominately in the weld metal. Long-term creep tests, including large samples more representative of the entire weldment, were fabricated and tested to times in excess of 10,000 h at multiple temperatures. Analysis of the creep data supports a WSRF of 0.9 for the solution annealed + aged weldments in contrast to the WSRF of 0.7 currently applied to welded + aged weldments. Metallurgical analysis shows that the solution annealing causes recrystallization of the weld metal and reduces the chemical compositional and microstructural gradients in the weldments resulting in acceptable short-term performance. In long-term creep, evidence for coarsened zones in the weldments due to discontinuous coarsening reactions were identified as the mechanism leading to accelerated creep damage formation in the weld metal and sample failure.

14 SOLAR ENERGY↗

NASA NTRS1982

Multiple-grid acceleration of Lax-Wendroff algorithms

A technique for accelerating the convergence of a one-step Lax-Wendroff method to a steady-state solution is discussed and its applicability extended to the more general class of two-step Lax-Wendroff methods. Several two-step methods which lead to quite efficient multiple grid algorithms are discussed. Computational results are presented using the full two dimensional Euler equations for both subcritical and shocked supercritical flows. Extensions and generalizations are mentioned.

Johnson, G. M.↗

DOE OSTI2024

Noise robust detection of quantum phase transitions

Quantum computing allows for the manipulation of highly correlated states whose properties quickly go beyond the capacity of any classical method to calculate. Thus one natural problem which could lend itself to quantum advantage is the study of ground-states of condensed matter models, and the transitions between them. However, current levels of hardware noise can require extensive application of error-mitigation techniques to achieve reliable computations. In this work, we use several IBM devices to explore a finite-size spin model with multiple “phaselike” regions characterized by distinct ground-state configurations. Using preoptimized Variational Quantum Eigensolver (VQE) solutions, we demonstrate that in contrast to calculating the energy, where zero-noise extrapolation is required in order to obtain qualitatively accurate yet still unreliable results, calculations of the energy derivative, two-site spin correlation functions, and the fidelity susceptibility yield accurate behavior across multiple regions, even with minimal or no application of error-mitigation approaches. Taken together, these sets of observables could be used to identify level crossings in a simple, noise-robust manner which is agnostic to the method of ground state preparation. This work shows promising potential for near-term application to identifying quantum phase transitions, including avoided crossings and nonadiabatic conical intersections in electronic structure calculations. Published by the American Physical Society 2024

Lively, Kevin (ORCID:0000000320981494)↗