Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solve time”

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 1,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Interaction of the sonic boom with atmospheric turbulence

Theoretical research has been carried out to study the effect of free-stream turbulence on sonic boom pressure fields. A new transonic small-disturbance model to analyze the interactions of random disturbances with a weak shock has been developed. The model equation has an extended form of the classic small-disturbance equation for unsteady transonic aerodynamics. An alternative approach shows that the pressure field may be described by an equation that has an extended form of the classic nonlinear acoustics equation that describes the propagation of sound beams with narrow angular spectrum. The model shows that diffraction effects, nonlinear steepening effects, focusing and caustic effects and random induced vorticity fluctuations interact simultaneously to determine the development of the shock wave in space and time and the pressure field behind it. A finite-difference algorithm to solve the mixed-type elliptic-hyperbolic flows around the shock wave has also been developed. Numerical calculations of shock wave interactions with various deterministic and random fluctuations will be presented in a future report.

Rusak, Z.↗
NASA NTRS1993

Passive microwave remote sensing of thin sea ice using principal component analysis

The possibility of developing a satellite-based sea-ice-concentration algorithm which solves for the presence of thinner ice is explored on the basis of time sequences of surface-based measurements of passive MW emission from growing saline ice reported by Wensnahan et al. (1993). It is shown that two classes of thinner ice can be distinguished from mixtures of open water, first-year ice, and multiyear ice. Newly formed ice appears to be optically thick at 37 and 90 GHz and has a relatively dry surface. The thin ice spectrum occurs when the ice is greater than 4 cm thick and appears to result from the accumulation of brine at the ice surface. Thin ice has a relatively stable spectrum characterized by high brightness temperatures, a near-zero spectral gradient at vertical polarization, and a large difference between vertical and horizontal polarizations. It is concluded that thin ice can be detected using satellite data.

Wensnahan, Mark↗
NASA NTRS1993

Optimal convolution SOR acceleration of waveform relaxation with application to semiconductor device simulation

In this paper we describe a novel generalized SOR (successive overrelaxation) algorithm for accelerating the convergence of the dynamic iteration method known as waveform relaxation. A new convolution SOR algorithm is presented, along with a theorem for determining the optimal convolution SOR parameter. Both analytic and experimental results are given to demonstrate that the convergence of the convolution SOR algorithm is substantially faster than that of the more obvious frequency-independent waveform SOR algorithm. Finally, to demonstrate the general applicability of this new method, it is used to solve the differential-algebraic system generated by spatial discretization of the time-dependent semiconductor device equations.

Reichelt, Mark↗
NASA NTRS1994

Interaction of the sonic boom with atmospheric turbulence

Theoretical research was carried out to study the effect of free-stream turbulence on sonic boom pressure fields. A new transonic small-disturbance model to analyze the interactions of random disturbances with a weak shock was developed. The model equation has an extended form of the classic small-disturbance equation for unsteady transonic aerodynamics. An alternative approach shows that the pressure field may be described by an equation that has an extended form of the classic nonlinear acoustics equation that describes the propagation of sound beams with narrow angular spectrum. The model shows that diffraction effects, nonlinear steepening effects, focusing and caustic effects and random induced vorticity fluctuations interact simultaneously to determine the development of the shock wave in space and time and the pressure field behind it. A finite-difference algorithm to solve the mixed type elliptic-hyperbolic flows around the shock wave was also developed. Numerical calculations of shock wave interactions with various deterministic and random fluctuations will be presented in a future report.

Rusak, Zvi↗
NASA NTRS1994

Sonic boom interaction with turbulence

A recently developed transonic small-disturbance model is used to analyze the interactions of random disturbances with a weak shock. The model equation has an extended form of the classic small-disturbance equation for unsteady transonic aerodynamics. It shows that diffraction effects, nonlinear steepening effects, focusing and caustic effects and random induced vorticity fluctuations interact simultaneously to determine the development of the shock wave in space and time and the pressure field behind it. A finite-difference algorithm to solve the mixed-type elliptic hyperbolic flows around the shock wave is presented. Numerical calculations of shock wave interactions with various deterministic vorticity and temperature disturbances result in complicate shock wave structures and describe peaked as well as rounded pressure signatures behind the shock front, as were recorded in experiments of sonic booms running through atmospheric turbulence.

Zvi Rusak↗
NASA NTRS1994

Magnetometer-only attitude and rate determination for a gyro-less spacecraft

Attitude determination algorithms that requires only the earth's magnetic field will be useful for contingency conditions. One way to determine attitude is to use the time derivative of the magnetic field as the second vector in the attitude determination process. When no gyros are available, however, attitude determination becomes difficult because the rates must be propagated via integration of Euler's equation, which in turn requires knowledge of the initial rates. The spacecraft state to be determined must then include not only the attitude but also rates. This paper describes a magnetometer-only attitude determination scheme with no a priori knowledge of the spacecraft state, which uses a deterministic algorithm to initialize an extended Kalman filter. The deterministic algorithm uses Euler's equation to relate the time derivatives of the magnetic field in the reference and body frames and solves the resultant transcendental equations for the coarse attitude and rates. An important feature of the filter is that its state vector also includes corrections to the propagated rates, thus enabling it to generate highly accurate solutions. The method was tested using in-flight data from the Solar, Anomalous, and Magnetospheric Particles Explorer (SAMPEX), a Small Explorer spacecraft. SAMPEX data using several eclipse periods were used to simulate conditions that may exist during the failure of the on-board digital sun sensor. The combined algorithm has been found effective, yielding accuracies of 1.5 deg in attitude (within even nominal mission requirements) and 0.01 degree per second (deg/sec) in the rates.

Natanson, G. A.↗
NASA NTRS1995

On the breakup of viscous liquid threads

A one-dimensional model evolution equation is used to describe the nonlinear dynamics that can lead to the breakup of a cylindrical thread of Newtonian fluid when capillary forces drive the motion. The model is derived from the Stokes equations by use of rational asymptotic expansions and under a slender jet approximation. The equations are solved numerically and the jet radius is found to vanish after a finite time yielding breakup. The slender jet approximation is valid throughout the evolution leading to pinching. The model admits self-similar pinching solutions which yield symmetric shapes at breakup. These solutions are shown to be the ones selected by the initial boundary value problem, for general initial conditions. Further more, the terminal state of the model equation is shown to be identical to that predicted by a theory which looks for singular pinching solutions directly from the Stokes equations without invoking the slender jet approximation throughout the evolution. It is shown quantitatively, therefore, that the one-dimensional model gives a consistent terminal state with the jet shape being locally symmetric at breakup. The asymptotic expansion scheme is also extended to include unsteady and inerticial forces in the momentum equations to derive an evolution system modelling the breakup of Navier-Stokes jets. The model is employed in extensive simulations to compute breakup times for different initial conditions; satellite drop formation is also supported by the model and the dependence of satellite drop volumes on initial conditions is studied.

Papageorgiou, Demetrios T.↗
NASA NTRS1995

More About Pressure Probes For Turbulent Flows

Array of probes yields data on static and total pressures and cross-stream turbulence. Is combination of separate pressure-probe arrays described in "Measuring Streamwise Momentum and Cross-Stream Turbulence" (ARC-11934), and "Probe Rakes To Measure Static Pressure and Turbulence" (ARC-12973). Equations that express directional responses of various probes in array solved together to eliminate unknown quantities; makes it possible to process digitized time-averaged probe readings together to obtain time-averaged static, stagnation, and cross-stream-turbulence pressures.

Rossow, Vernon J.↗
NASA NTRS1999

Use of Boundary Layer Transition Detection to Validate Full-Scale Flight Performance Predictions

Full-scale flight performance predictions can be made using CFD or a combination of CFD and analytical skin-friction predictions. However, no matter what method is used to obtain full-scale flight performance predictions knowledge of the boundary layer state is critical. The implementation of CFD codes solving the Navier-Stokes equations to obtain these predictions is still a time consuming, expensive process. In addition, to ultimately obtain accurate performance predictions the transition location must be fixed in the CFD model. An example, using the M2.4-7A geometry, of the change in Navier-Stokes solution with changes in transition and in turbulence model will be shown. Oil flow visualization using the M2.4-7A 4.0% scale model in the 14'x22' wind tunnel shows that fixing transition at 10% x/c in the CFD model best captures the flow physics of the wing flow field. A less costly method of obtaining full-scale performance predictions is the use of non-linear Euler codes or linear CFD codes, such as panel methods, combined with analytical skin-friction predictions. Again, knowledge of the boundary layer state is critical to the accurate determination of full-scale flight performance. Boundary layer transition detection has been performed at 0.3 and 0.9 Mach numbers over an extensive Reynolds number range using the 2.2% scale Reference H model in the NTF. A temperature sensitive paint system was used to determine the boundary layer state for these conditions. Data was obtained for three configurations: the baseline, undeflected flaps configuration; the transonic cruise configuration; and, the high-lift configuration. It was determined that at low Reynolds number conditions, in the 8 to 10 million Reynolds number range, the baseline configuration has extensive regions of laminar flow, in fact significantly more than analytical skin-friction methods predict. This configuration is fully turbulent at about 30 million Reynolds number for both 0.3 and 0.9, Mach numbers. Both the transonic cruise and the high-lift configurations were fully turbulent aft of the leading-edge flap hingeline at all Reynolds numbers.

Hamner, Marvine↗
NASA NTRS2001

Computational Aeroacoustics by the Space-time CE/SE Method

In recent years, a new numerical methodology for conservation laws-the Space-Time Conservation Element and Solution Element Method (CE/SE), was developed by Dr. Chang of NASA Glenn Research Center and collaborators. In nature, the new method may be categorized as a finite volume method, where the conservation element (CE) is equivalent to a finite control volume (or cell) and the solution element (SE) can be understood as the cell interface. However, due to its rigorous treatment of the fluxes and geometry, it is different from the existing schemes. The CE/SE scheme features: (1) space and time treated on the same footing, the integral equations of conservation laws are solve( for with second order accuracy, (2) high resolution, low dispersion and low dissipation, (3) novel, truly multi-dimensional, simple but effective non-reflecting boundary condition, (4) effortless implementation of computation, no numerical fix or parameter choice is needed, an( (5) robust enough to cover a wide spectrum of compressible flow: from weak linear acoustic waves to strong, discontinuous waves (shocks) appropriate for linear and nonlinear aeroacoustics. Currently, the CE/SE scheme has been developed to such a stage that a 3-13 unstructured CE/SE Navier-Stokes solver is already available. However, in the present paper, as a general introduction to the CE/SE method, only the 2-D unstructured Euler CE/SE solver is chosen as a prototype and is sketched in Section 2. Then applications of the CE/SE scheme to linear, nonlinear aeroacoustics and airframe noise are depicted in Sections 3, 4, and 5 respectively to demonstrate its robustness and capability.

Loh, Ching Y.↗
NASA NTRS2005

A First and Second Order Moment Approach to Probabilistic Control Synthesis

This paper presents a robust control design methodology based on the estimation of the first two order moments of the random variables and processes that describe the controlled response. Synthesis is performed by solving an multi-objective optimization problem where stability and performance requirements in time- and frequency domains are integrated. The use of the first two order moments allows for the efficient estimation of the cost function thus for a faster synthesis algorithm. While reliability requirements are taken into account by using bounds to failure probabilities, requirements related to undesirable variability are implemented by quantifying the concentration of the random outcome about a deterministic target. The Hammersley Sequence Sampling and the First- and Second-Moment- Second-Order approximations are used to estimate the moments, whose accuracy and associated computational complexity are compared numerically. Examples using output-feedback and full-state feedback with state estimation are used to demonstrate the ideas proposed.

Crespo, Luis G.↗
NASA NTRS2019

Series solution for motion in an arbitrary potential field

The equations of motion for a test particle in an arbitrary potential field is solved as a formal power series. This solution is in terms of the time derivatives of the motion of the test particle. These derivative may be used to construct a Taylor series for the motion of the particle in the vicinity of given initial conditions.

Strange, Nathan↗
NASA NTRSDate not supplied

Full Stall Simulations of a Redesigned Ventilation Fan for the ISS

The concept of a stall is studied rigorously in the aerospace industry. From a design standpoint, instabilities such as stall are undesirable– operation in the stall regime has a tremendous impact on aerodynamic performance as well as structural integrity. In extreme cases, operating in stall conditions can cause failure. Further, stall can cause a loss of lift on aircraft wings or a loss of thrust in aircraft engines. In any case, stall continues to be a topic of interest in the aerospace industry. The present work aims to analyze the stall characteristics of a ventilation fan that was recently designed for the International Space Station (ISS). Although the ventilation fan has a rotor-stator design, this paper considers a rotor-only configuration. The FUN3D Computational Fluid Dynamic (CFD) solver developed by NASA Langley Research Center was used to simulate the operational characteristics of the ventilation fan. FUN3D solves the Unsteady Reynolds-Averaged Naiver-Stokes (URANS) equations using implicit time marching and a dynamic overset grid. The FUN3D solver was originally written for exterior flow fields; however, this work represents an extension of the FUN3D solver to turbomachinery or interior flow fields. The FUN3D results for the rotor-only ventilation fan accurately captured the operational characteristics inherent to compressors– the results shared similar performance trends when compared to the experimental results for the rotor-stator case. A peak adiabatic efficiency of 96% occurred at a MFR of 105.2 CFM, the minimum aerodynamically stable point. A computationally stable stall occurred at a mass flow rate of 43.8 CFM where the adiabatic efficiency dropped to 69%.

CFD↗
NASA NTRS2023

Aviation 2072: Scenario Planning for Wicked Problems

The world is changing rapidly. And with it, so are the potential futures for humanity. Many of these potential futures pose hidden threats, while others offer new opportunities for aviation and the broader aerospace community. To mitigate these risks and capitalize on opportunities, organizations must work to deeply understand these potential future scenarios and uncover any underlying drivers. Organizations will be more equipped and informed when making strategic investment decisions by better understanding these threats, opportunities, and drivers . The future scenarios and their threats and opportunities described in this paper were identified through a series of extensive brainstorming workshops, collectively titled MADNESS (Mapping for Aviation Driven by Needs Emergence and Satiating Society). These workshops engaged NASA civil servants and contractors with diverse backgrounds through a facilitated process of future scenario development and critique. The participants were led through exercises focused on a time horizon spanning from 2022 to 2072 and on problems that aviation might solve and or create during this period. From the broad results, two patterns of critical uncertainty emerged: “availability” (scarcity versus abundance) and “transparency” (openness versus security and privacy). The successful iteration and ideation across the MADNESS workshops also suggests a repeatable mechanism for strategic risk exploration across a variety of NASA and industry stakeholders.

Mapping↗
NASA NTRS2003

Shared Activity Coordination

Interacting agents that interleave planning and execution must reach consensus on their commitments to each other. In domains where agents have varying degrees of interaction and different constraints on communication and computation, agents will require different coordination protocols in order to efficiently reach consensus in real time. We briefly describe a largely unexplored class of real-time, distributed planning problems (inspired by interacting spacecraft missions), new challenges they pose, and a general approach to solving the problems. These problems involve self-interested agents that have infrequent communication but collaborate on joint activities. We describe a Shared Activity Coordination (SHAC) framework that provides a decentralized algorithm for negotiating the scheduling of shared activities in a dynamic environment, a soft, real-time approach to reaching consensus during execution with limited communication, and a foundation for customizing protocols for negotiating planner interactions. We apply SHAC to a realistic simulation of interacting Mars missions and illustrate the simplicity of protocol development.

coordination↗
DOE OSTI2025

A time-parallel multiple-shooting method for large-scale quantum optimal control

Quantum optimal control plays a crucial role in quantum computing by providing the interface between compiler and hardware. Solving the optimal control problem is particularly challenging for multi-qubit gates, due to the exponential growth in computational complexity with the system's dimensionality and the deterioration of optimization convergence. To ameliorate the computational complexity of time-integration, this paper introduces a multiple-shooting approach in which the time domain is divided into multiple windows and the intermediate states at window boundaries are treated as additional optimization variables. Further, this enables parallel computation of state evolution across time-windows, significantly accelerating objective function and gradient evaluations. Since the initial state matrix in each window is only guaranteed to be unitary upon convergence of the optimization algorithm, the conventional gate trace infidelity is replaced by a generalized infidelity that is convex for non-unitary state matrices. Continuity of the state across window boundaries is enforced by equality constraints. A quadratic penalty optimization method is used to solve the constrained optimal control problem, and an efficient adjoint technique is employed to calculate the gradients in each iteration. We demonstrate the effectiveness of the proposed method through numerical experiments on quantum Fourier transform gates in systems with 2, 3, and 4 qubits, noting a speedup of 80x for evaluating the gradient in the 4-qubit case, highlighting the method's potential for optimizing control pulses in multi-qubit quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing is developed using a general buoyancy–shear–drag model. Analytic solutions to the simplified model equations corresponding to each instability are derived separately, which are then used to calibrate the model coefficients to predict specific values of the mixing layer growth parameters and exponents. The buoyancy-shear–drag equations for the bubble and spike mixing layer widths are then solved numerically, and a turbulent diffusivity (or viscosity) is constructed dimensionally as the product of the mixing layer width, h, and its time-derivative, dh/dt. Surrogate turbulent fields are then constructed by multiplying a presumed, approximate self-similar spatial profile by appropriate functions of h and dh/dt. Using several simplifying approximations, the turbulent diffusion equations satisfied by the mean mass fraction and mean shear velocity are solved analytically in a reference frame moving with the mean advection velocity, and it is shown that the mean fields evolve in space and time as expected. The explicit modeling and solution of turbulent transport equations (e.g., for the turbulent kinetic energy and its dissipation rate or a lengthscale) are not required. By using separate equations for the bubble and spike evolution, the model naturally captures the asymmetry induced by an increasing density contrast between the heavy and light fluids. The model also includes molecular dissipation and diffusion, and therefore, can describe transition to fully-developed turbulence if the initial turbulent diffusivity/viscosity is much smaller than the molecular diffusivity/viscosity. This model is applied to constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and variable-density Kelvin–Helmholtz mixing layers to demonstrate its utility. It is shown that the numerical solutions of the model calibrated using specific values of the instability growth parameters and exponents: (1) produces mixing layer widths in agreement with the expected self-similar growth power laws; (2) gives spatiotemporal profiles of turbulent fields that are expected and consistent with previous results; (3) predicts the expected power-law growths and decays of the spatially-integrated turbulent fields, and; (4) gives spatiotemporal profiles of the mean fields that are expected and consistent with previous results. This relatively simple model can be implemented numerically, and is expected to be useful for fundamental as well as astrophysical and inertial confinement fusion applications.

42 ENGINEERING↗
NASA NTRS1953

Transient temperature distributions in simple conducting bodies steadily heated through a laminar boundary layer

An analysis is made of the transient heat-conduction effects in three simple semi-infinite bodies: the flat insulated plate, the conical shell, and the slender solid cone. The bodies are assumed to have constant initial temperatures and, at zero time, to begin to move at a constant speed and zero angle of attack through a homogeneous atmosphere. The heat input is taken as that through a laminar boundary layer. Radiation heat transfer and transverse temperature gradients are assumed to be zero. The appropriate heat-conduction equations are solved by an iteration method, the zeroeth-order terms describing the situation in the limit of small time. The method is presented and the solutions are calculated to three orders which are sufficient to give reasonably accurate results when the forward edge has attained one-half the total temperature rise (nose half-rise time). Flight Mach number and air properties occur as parameters in the result. Approximate expressions for the extent of the conduction region and nose half-rise times as functions of the parameters of the problem are presented. (author)

Parker, Hermon M↗