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,333 records · Page 74

Conjugate Compressible Fluid Flow and Heat Transfer in Ducts

A computational approach to modeling transient, compressible fluid flow with heat transfer in long, narrow ducts is presented. The primary application of the model is for analyzing fluid flow and heat transfer in solid propellant rocket motor nozzle joints during motor start-up, but the approach is relevant to a wide range of analyses involving rapid pressurization and filling of ducts. Fluid flow is modeled through solution of the spatially one-dimensional, transient Euler equations. Source terms are included in the governing equations to account for the effects of wall friction and heat transfer. The equation solver is fully-implicit, thus providing greater flexibility than an explicit solver. This approach allows for resolution of pressure wave effects on the flow as well as for fast calculation of the steady-state solution when a quasi-steady approach is sufficient. Solution of the one-dimensional Euler equations with source terms significantly reduces computational run times compared to general purpose computational fluid dynamics packages solving the Navier-Stokes equations with resolved boundary layers. In addition, conjugate heat transfer is more readily implemented using the approach described in this paper than with most general purpose computational fluid dynamics packages. The compressible flow code has been integrated with a transient heat transfer solver to analyze heat transfer between the fluid and surrounding structure. Conjugate fluid flow and heat transfer solutions are presented. The author is unaware of any previous work available in the open literature which uses the same approach described in this paper.

Cross, M. F.↗

Using Markov Models of Fault Growth Physics and Environmental Stresses to Optimize Control Actions

A generalized Markov chain representation of fault dynamics is presented for the case that available modeling of fault growth physics and future environmental stresses can be represented by two independent stochastic process models. A contrived but representatively challenging example will be presented and analyzed, in which uncertainty in the modeling of fault growth physics is represented by a uniformly distributed dice throwing process, and a discrete random walk is used to represent uncertain modeling of future exogenous loading demands to be placed on the system. A finite horizon dynamic programming algorithm is used to solve for an optimal control policy over a finite time window for the case that stochastic models representing physics of failure and future environmental stresses are known, and the states of both stochastic processes are observable by implemented control routines. The fundamental limitations of optimization performed in the presence of uncertain modeling information are examined by comparing the outcomes obtained from simulations of an optimizing control policy with the outcomes that would be achievable if all modeling uncertainties were removed from the system.

Bole, Brian↗

A Summary of the Space-Time Conservation Element and Solution Element (CESE) Method

The space-time Conservation Element and Solution Element (CESE) method for solving conservation laws is examined for its development motivation and design requirements. The characteristics of the resulting scheme are discussed. The discretization of the Euler equations is presented to show readers how to construct a scheme based on the CESE method. The differences and similarities between the CESE method and other traditional methods are discussed. The strengths and weaknesses of the method are also addressed.

Fluid Dynamics↗

Unexpected Nonlinear Effects in Superconducting Transition-Edge Sensors

When a normal metal transitions into the superconducting state the DC resistance drops from a finite value to zero over some finite transition width in temperature, current, and magnetic field. Superconducting transition-edge sensors (TESs) operate within this transition region and uses resistive changes to measure deposited thermal energy. This resistive transition is not perfectly smooth and a wide range of TES designs and materials show sub-structure in the resistive transition (as seen in smooth nonmonotonic behavior, jump discontinuities, and hysteresis in the devices current-voltage relation and derivatives of the resistance with respect to temperature, bias current, and magnetic field). TES technology has advanced to the point where for many applications this structure is the limiting factor in performance and optimization consists of finding operating points away from these structures. For example, operating at or near this structure can lead to nonlinearity in the detectors response and gain scale, limit the spectral range of the detector by limiting the usable resistive range, and degrade energy resolution. The origin of much of this substructure is unknown. This presentation investigates a number of possible sources in turn. First we model the TES as a superconducting weak-link and solve for the characteristic differential equations current and voltage time dependence. We find:(1) measured DC biased current-voltage relationship is the time-average of a much higher frequency limit cycle solution.(2) We calculate the fundamental frequency and estimate the power radiated from the TES treating the bias leads as an antennae.(3) The solution for a set of circuit parameters becomes multivalued leading to current transitions between levels.(4)The circuit parameters can change the measure resistance and mask the true critical current. As a consequence the TES resistance surface is not just a function of temperature, current, and magnetic field but is also a function of the circuit elements (such as shunt resistor, SQUID inductance, and capacitor values). In other words, same device measured in different electrical circuits will have a different resistive surface in temperature, current, and magnetic field. Next we consider that at the transition temperature of a superconductor both the magnetic penetration depth and coherence length are divergent. As a consequence these important characteristic length scales are changing with operating point. We present measurements on devices showing commensurate behavior between these characteristic lengths and the length scale of added normal metal structures. Reordering of proximity vortices leads to discontinuities and irreversibility of the current-voltage curves. Last we consider a weak-link TES including both thermal activated resistance effects and the effect of the magnetic penetration depth being a function of temperature and magnetic field. We derive its impact on the resistive transition surface and the important device parameters a and b.

Superconducting Detectors↗

Rapid Assessment of a Lunar Surface Laboratory Module

Johnson Space Center is piloting an innovation team known as the Forge to perform rapid turnaround, concurrent engineering studies. The Forge is intended to be able to pull together a team of strategically selected domain experts to innovatively solve a specified problem within a short period of time. The Forge’s inaugural project being the development of a Lunar Surface Laboratory Module concept. The lab module is intended as an option for the Artemis program’s Lunar Architecture Team. The lab module may allow for increased science productivity and return on investment than a surface architecture with only the habitat. This paper will discuss the configuration of the Lunar Laboratory module, including the pressure vessel, its placement relative to the Surface Habitat and Pressurized Rover, and internal architecture. In order to determine the suitability of such a design, a human factors evaluation is conducted. Historically, a human factors evaluation of a habitable volume could require a month or longer, including time needed to prepare a test plan, conduct the evaluation, compile the data, and produce a report. The entire Forge effort will only last for five weeks, forcing the creation of a rapid assessment process. This rapid assessment includes three evaluation studies, all conducted and processed in the space of a single day, to determine the acceptability of the resulting design. The first evaluation is a card sort evaluation of the laboratory instruments. A brief discussion of the card sort evaluation is provided. This card sort is focused on the relative positioning of the laboratory instruments, determining whether a given instrument should be adjacent to, or separated from, another. It will use human factors personnel and both life and physical scientists as subject matter experts. The second evaluation is a science productivity tabletop evaluation. This evaluation focuses on the selection of instruments and the number of instruments, attempting to assess the amount of work the crew will be able to use the facility to do. It uses life and physical science personnel as subject matter experts. The third evaluation is a CAD or Virtual Reality human factors walkthrough evaluation. This evaluation focuses on the layout of the lab module and its connectivity to the basecamp’s pressurized elements. If possible, astronauts will be used as subject matter experts. Each of these evaluations will focus primarily on objective data with limited subjective comments, reducing the time required for processing responses. The results of each evaluation and their implications will be discussed. Based on the evaluation data, an overall acceptability of the Lunar Surface Laboratory Module produced by the Forge study can be determined. Lessons learned and recommendations for future use of the human factors rapid assessment process in early design trades will also be provided.

Lunar Laboratory↗

A Rapid Assessment of a Lunar Surface Laboratory Module

Johnson Space Center has piloted an innovation team known as the Forge to perform rapid turnaround, concurrent engineering studies. The Forge is intended to be able to pull together a team of strategically selected domain experts to innovatively solve a specified problem within a short period of time. The Forge’s inaugural project was the development of a Lunar Surface Laboratory Module concept. The lab module is intended as an option for the NASA Lunar Architecture Team and may allow for increased science productivity and return on investment than with only the Surface Habitat. The Forge team developed two concepts, a stationary laboratory and a mobile laboratory. Human factors evaluations have been conducted to determine the acceptability of the resulting design concepts. The human factors evaluations conducted include a card sort evaluation of the laboratories by science discipline and the System Usability Scale (SUS). The card sort focused subject matter expert expectations for the relative positioning of the laboratory functional areas, determining whether a given lab area should be adjacent to, or separated from, another. The SUS focused on layout efficiency. Both used human factors personnel and both life and physical scientists as subject matter experts. Each of these evaluations will focus primarily on objective data with limited subjective comments, reducing the time required for processing responses. The results of the card sort and SUS evaluations and their implications will be discussed. Based on the evaluation data, an overall acceptability of the Lunar Surface Laboratory Module produced by the Forge study can be determined. Lessons learned and recommendations for future use of human factors rapid assessment processes in early design trades will also be provided.

Forge↗

Adaptive Optimization for System Performance and Combined Bernstein Polynomial, Optimal Reciprocal Collision Avoidance, Differential Dynamic Programming for Trajectory Replanning and Collision Avoidance for UAM Vehicles

The emerging urban air mobility (UAM) sector in aerospace is driving development of unconventional multi-modal vehicle configurations and autonomous flight. The combination of multi-modal vehicle dynamics, complex environment, requirements to deal with flight contingencies in an efficient and safe manner, as well as necessity for precise trajectory following and performance, are the driving influence behind adaptive optimization for system performance. We are interested in trajectory optimization algorithm that would system parameter estimation and identifying the optimal switching time between modes of hybrid dynamical systems. This presentation discusses a parameterized optimal control trajectory optimization algorithm that is an extended and generalized version of Differential Dynamic Programming (DDP), titled Parameterized Differential Dynamic Programming (PDDP). DDP is an efficient trajectory optimization algorithm relying on second order approximations of a system’s dynamics and cost function and has recently been applied to optimize systems with time invariant parameters. Experiments are presented applying PDDP to solve model predictive control (MPC) and moving horizon estimation (MHE) tasks simultaneously. In particular, PDDP is used to determine the optimal transition point between flight regimes of a complex urban air mobility (UAM) class vehicle exhibiting multiple phases of flight and to identify and compensate for actuation faults.

optimization↗

Rapid Assessment of a Lunar Surface Laboratory Module

Johnson Space Center has piloted an innovation team known as the Forge to perform rapid turnaround, concurrent engineering studies. The Forge is intended to be able to pull together a team of strategically selected domain experts to innovatively solve a specified problem within a short period of time. The Forge’s inaugural project was the development of a Lunar Surface Laboratory Module concept. The lab module is intended as an option for the NASA Lunar Architecture Team and may allow for increased science productivity and return on investment than with only the Surface Habitat. The Forge team developed two concepts, a stationary laboratory and a mobile laboratory. Human factors evaluations have been conducted to determine the acceptability of the resulting design concepts. The human factors evaluations conducted include a card sort evaluation of the laboratories by science discipline and the System Usability Scale (SUS). The card sort is focused subject matter expert expectations for the relative positioning of the laboratory functional areas, determining whether a given lab area should be adjacent to, or separated from, another. The SUS will focus on layout efficiency. Both will use human factors personnel and both life and physical scientists as subject matter experts. Each of these evaluations will focus primarily on objective data with limited subjective comments, reducing the time required for processing responses. The results of the card sort and SUS evaluations and their implications will be discussed. Based on the evaluation data, an overall acceptability of the Lunar Surface Laboratory Module produced by the Forge study can be determined. Lessons learned and recommendations for future use of human factors rapid assessment processes in early design trades will also be provided.

Human Factors↗

Discretization Error Estimation and Control for Farfield Acoustic Signatures

We investigate the utility of adjoint-based error estimates for sonic boom farfield simulations governed by solutions of the augmented Burgers’ equation. Solution of this nonlinear system uses operator splitting with a second-order finite volume discretization in space and second-order Runge-Kutta time marching, while the absorption and molecular relaxation are solved using second-order central differencing. The discretization error in selected ground sonic boom cost functionals is estimated using the method of adjoint-weighted residuals. Key elements of the implementation process are emphasized with details provided on the practical aspects as appliedto the sonic boom farfield propagation. We establish the accuracy of the adjoint solutions usingcomplex step and finite difference approaches, and examine the accuracy of the error estimates using analytical N-wave solutions. We then apply it to a pressure waveform corresponding to the X-59 research aircraft. The investigations demonstrate that the method of adjoint-weighted residuals accurately predicts the level of discretization error present in sonic boom farfield simulations while offering insight into which features of the near field signal are the primary drivers of ground noise metrics. The numerical results indicate that at sampling frequencies as low as50kHz, discretization error in the propagation is under 0.01 dB[A] for realistically complex examples.

CST↗

Design of a Thrust Vector Servoelastic Filter for SLS Block 2 Booster Thrust Vector Control System

This paper presents a method to improve the thrust vector servoelastic (TVSE) response of a launch vehicle. TVSE is the coupling between the Thrust Vector Control (TVC) system nozzle dynamics and the vehicle flexible modes. A poorly phased TVC response, i.e. in phase with the structural velocity, can lead to instability of the flexible modes due to the energy injection from the thrust. For large boosters with heavy nozzles, TVSE is more likely to cause instability. This instability is part of the plant dynamics and may not be readily detected by the flight control system. In this paper, TVSE is treated as a disturbance rejection problem with the inertial force of the flexible motion viewed as an external disturbance to the TVC actuator control system. By effectively rejecting this disturbance, the engine can be decoupled from the flexible dynamics thus improving TVSE stability. The actuator disturbance transfer function is modified by inserting a filter into the servo-loop control system. The filter aims to reduce the magnitude and introduce phase lead of the disturbance transfer function over the TVSE sensitive frequency range while maintaining adequate stability margins of the actuator control system. A prefilter is also designed to improve command tracking performance. The filter and prefilter resulting from this method are usually of low order and can be included to an existing actuator controller if the controller is programmable. This method was demonstrated for a hydraulic Simplex actuator model. Its applications to the Space Launch System (SLS) Block 1 and Block 2 TVC actuator are also described. The Block 2 actuator controller was modified based on the results of this method. The tradeoff between TVSE stability and the servo-loop stability can be balanced by adjusting the amount of phase lead introduced by the TVSE filter. This method provides a useful way to solve the TVSE instability issue without costly and time-consuming hardware solutions. This method is based on the basic control principles of loop shaping and disturbance rejection and has been proven to be effective under proper circumstances.

Thrust vector control↗

Design of a Thruster Vector Servo-Elastic Filter

This presents a method to improve the thruster vector servo elasticity (TVSE) stability. TVSE is the coupling between the TVC system engine dynamics and the vehicle flexible modes. It can lead to instability of the flexible modes due to the energy injection from the thrust. For large boosters with heavy TVC engine, TVSE is more likely to cause instability. This instability is part of the plant dynamics and is not detected by the flight control. In this paper, the TVSE is treated as a disturbance rejection problem with the inertial force of the flexible motion viewed as an external disturbance to the TVC actuator control system. By effectively rejecting this disturbance, the engine can be decoupled from the flexible dynamics thus improving TVSE stability. The actuator disturbance transfer function is modified by inserting a TVSE filter into the control system. The filter aims to reduce the magnitude and introduce phase lead of the disturbance transfer function over the TVSE sensitive frequency range while maintaining adequate stability margins of the actuator control system. A prefilter is also designed for good command tracking. The filter and prefilter resulting from this method are usually of low order and can be included to the existing actuator controller if the controller is programmable. This method was demonstrated for a hydraulic Simplex actuator model. Its applications to SLS B1 and B2 TVC actuator are also described. B2 actuator controller is modified based on this method. The tradeoff between TVSE stability and the control system stability can be balanced by adjusting the amount of phase lead introduced by the TVSE filter. This method provides a useful way to solve the TVSE instability issue without costly and time-consuming a TVC system redesign. This method is based on the basic control principles of loop shaping and disturbance rejection and has been proven to be effective under proper circum-stances.

Thrust vector control↗

On recursive least-squares filtering algorithms and implementations

In many real-time signal processing applications, fast and numerically stable algorithms for solving least-squares problems are necessary and important. In particular, under non-stationary conditions, these algorithms must be able to adapt themselves to reflect the changes in the system and take appropriate adjustments to achieve optimum performances. Among existing algorithms, the QR-decomposition (QRD)-based recursive least-squares (RLS) methods have been shown to be useful and effective for adaptive signal processing. In order to increase the speed of processing and achieve high throughput rate, many algorithms are being vectorized and/or pipelined to facilitate high degrees of parallelism. A time-recursive formulation of RLS filtering employing block QRD will be considered first. Several methods, including a new non-continuous windowing scheme based on selectively rejecting contaminated data, were investigated for adaptive processing. Based on systolic triarrays, many other forms of systolic arrays are shown to be capable of implementing different algorithms. Various updating and downdating systolic algorithms and architectures for RLS filtering are examined and compared in details, which include Householder reflector, Gram-Schmidt procedure, and Givens rotation. A unified approach encompassing existing square-root-free algorithms is also proposed. For the sinusoidal spectrum estimation problem, a judicious method of separating the noise from the signal is of great interest. Various truncated QR methods are proposed for this purpose and compared to the truncated SVD method. Computer simulations provided for detailed comparisons show the effectiveness of these methods. This thesis deals with fundamental issues of numerical stability, computational efficiency, adaptivity, and VLSI implementation for the RLS filtering problems. In all, various new and modified algorithms and architectures are proposed and analyzed; the significance of any of the new method depends crucially on specific application.

Hsieh, Shih-Fu↗

Direct inference of nuclear equation-of-state parameters from gravitational-wave observations

The observation of neutron star mergers with gravitational waves (GWs) has provided a new method to constrain the dense-matter equation of state (EOS) and to better understand its nuclear physics. However, inferring nuclear microphysics from GW observations necessitates the sampling of EOS model parameters that serve as input for each EOS used during the GW data analysis. The sampling of the EOS parameters requires solving the Tolman–Oppenheimer–Volkoff (TOV) equations a large number of times—a process that slows down each likelihood evaluation in the analysis on the order of a few seconds. Here, we employ emulators for the TOV equations built using multilayer perceptron neural networks to enable direct inference of nuclear EOS parameters from GW strain data. Our emulators allow us to rapidly solve the TOV equations, taking in EOS parameters and outputting the associated tidal deformability of a neutron star in only a few tens of milliseconds. We implement these emulators in PyCBC to directly infer the EOS parameters using the event GW170817, providing posteriors on these parameters informed solely by GWs. We benchmark these runs against analyses performed using the full TOV solver and find that the emulators achieve speed ups of nearly two orders of magnitude, with negligible differences in the recovered posteriors. Additionally, we constrain the slope and curvature of the symmetry energy at the 90% upper credible interval to be $L$ sym ≲ 106 MeV and $K$ sym ≲ 26 MeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Evaluation of a transfinite element numerical solution method for nonlinear heat transfer problems

Laplace transform techniques have been widely used to solve linear, transient field problems. A transform-based algorithm enables calculation of the response at selected times of interest without the need for stepping in time as required by conventional time integration schemes. The elimination of time stepping can substantially reduce computer time when transform techniques are implemented in a numerical finite element program. The coupling of transform techniques with spatial discretization techniques such as the finite element method has resulted in what are known as transfinite element methods. Recently attempts have been made to extend the transfinite element method to solve nonlinear, transient field problems. This paper examines the theoretical basis and numerical implementation of one such algorithm, applied to nonlinear heat transfer problems. The problem is linearized and solved by requiring a numerical iteration at selected times of interest. While shown to be acceptable for weakly nonlinear problems, this algorithm is ineffective as a general nonlinear solution method.

Cerro, J. A.↗

Diffusive and radiative effects on vaporization times of drops in film boiling

Diffusive and radiative effects are incorporated into an analysis for the vaporization time of drops in film boiling. The momentum, energy, and continuity equations are solved with some appropriate simplifications so as to obtain a simple closed form solution for the overall film boiling heat transfer coefficient. Next, a theoretical expression for the droplet vaporization time is developed and compared to the measured vaporization times of water droplets vaporizing into air, argon, nitrogen, and helium. The agreement between experiment and theory is good. Under the helium blanket, the diffusive evaporative component is significant in comparison to the film boiling component.

Baumeister, K. J.↗

Multi-Robot Assembly Scheduling for the Lunar Crater Radio Telescope on the Far-Side of the Moon

The Lunar Crater Radio Telescope (LCRT) is a pro- posed ultra-long-wavelength radio telescope to be constructed on the far side of the moon. The proposed telescope will be constructed by deploying a 1km wire mesh in a 3-5km crater using a team of wall-climbing DuAxel robots. In this work, we consider the problem of generating minimum-time assembly sequences for LCRT, using realistic models of travel speed and lighting. Specifically, we pose the assembly sequencing problem as a mixed-integer linear program (MILP), which we solve to global optimality using commercial solvers. We present methods for modeling time-varying travel and assembly times, based on variable lighting conditions (including crater shadowing), and show how such time-varying parameters can be incorporated into the MILP. Finally, we present numerical studies of our method, showing how makespan varies with the number of assembly robots.

Schwager, Mac↗

Network-Wide Traffic Signal Control Using Bilinear System Modeling and Adaptive Optimization

This study proposes a new multi-input multi-output optimal bilinear signal control method in which a bilinear dynamic model approximation is used to capture the nonlinear dynamics of the urban traffic networks. With signal green time splits as the control input and traffic delay changes as the output for each intersections in the network, a bilinear system model was developed, which, on the basis of linear system modeling, takes interactions among traffic delays and signal timing splits into consideration. Based on the bilinear system modeling framework, we conducted two steps in each time interval to derive traffic control strategies: (1) we used the normalized least-squared algorithm to estimate system parameters; and (2) we solved an online optimization problem to obtain the updated traffic control inputs for the signal timing that minimizes future traffic delays. We evaluated the proposed method in a microscopic traffic simulation environment (VISSIM) with a 35-intersection network of Bellevue city in Washington. Two different traffic demand patterns: (1) normal traffic demands; and (2) time-varying traffic demands were simulated to compare the performance of different control strategies. Experimental results show that (1) the proposed bilinear system model can better describe traffic system dynamics than linear-model based methods, such as our previously developed linear-quadratic regulator control; and (2) the proposed method outperforms the state-of-the-art signal control strategies, namely the max-pressure and the self-organizing traffic light control methods. We have also shown that the proposed method is applicable to all other possible network layouts and signal controller phasing structures.

42 ENGINEERING↗

Application of the space-time conservation element and solution element method to two-dimensional advection-diffusion problems

The existing 2-D alpha-mu scheme and alpha-epsilon scheme based on the method of space-time conservation element and solution element, which were constructed for solving the linear 2-D unsteady advection-diffusion equation and unsteady advection equation, respectively, are tested. Also, the alpha-epsilon scheme is modified to become the V-E scheme for solving the nonlinear 2-D inviscid Burgers equation. Numerical solutions of six test problems are presented in comparison with their exact solutions or numerical solutions obtained by traditional finite-difference or finite-element methods. It is demonstrated that the 2-D alpha-mu, alpha-epsilon, and nu-epsilon schemes can be used to obtain numerical results which are more accurate than those based on some of the traditional methods but without using any artificial tuning in the computation. Similar to the previous 1-D test problems, the high accuracy and simplicity features of the space-time conservation element and solution element method have been revealed again in the present 2-D test results.

Wang, Xiao-Yen↗