Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iteration method”

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,729 records · Page 96

Numerical simulation of an electrothermal deicer pad

A numerical simulation is developed to investigate the removal of ice from composite aircraft blades by means of electrothermal deicing. The model considers one dimensional, unsteady state heat transfer in the composite blade-ice body. The heat conduction equations are approximated by using the Crank-Nicolson finite difference scheme, and the phase change in the ice layer is handled using the Enthalpy method. To solve the system of equations which result, Gauss-Seidel iteration is used. The simulation computes the temperature profile in the composite blade-ice body, as well as the movement of the ice-water interface, as a function of time. This information can be used to evaluate deicer performance. The simulation can also be used to solve a variety of other heat conduction problems involving composite bodies.

Marano, J. J.↗

Predicting Wall Modifications for Adaptive Wind Tunnels

Wall shape changed iteratively until it matches streamlines. FLEXWAL predicts upper and lower wall modifications necessary to remove wall-interference effects in adaptive-wall wind tunnels. FLEXWAL aids in elimination of wall-interference effects on objects tested in typical two-dimensional wind tunnel with rigid sidewalls and flexible, solid floor and ceiling boundaries. Iterative procedure valid for subsonic and transonic test conditions, and convergence of method verified both analytically and experimentally. FLEXWAL written in FORTRAN IV for batch execution.

Everhart, J. L.↗

Investigation to advance prediction techniques of the low-speed aerodynamics of V/STOL aircraft

A computer program, VSAERO, has been applied to a number of V/STOL configurations with a view to advancing prediction techniques for the low-speed aerodynamic characteristics. The program couples a low-order panel method with surface streamline calculation and integral boundary layer procedures. The panel method--which uses piecewise constant source and doublet panels-includes an iterative procedure for wake shape and models boundary layer displacement effect using the source transpiration technique. Certain improvements to a basic vortex tube jet model were installed in the code prior to evaluation. Very promising results were obtained for surface pressures near a jet issuing at 90 deg from a flat plate. A solid core model was used in the initial part of the jet with a simple entrainment model. Preliminary representation of the downstream separation zone significantly improve the correlation. The program accurately predicted the pressure distribution inside the inlet on the Grumman 698-411 design at a range of flight conditions. Furthermore, coupled viscous/potential flow calculations gave very close correlation with experimentally determined operational boundaries dictated by the onset of separation inside the inlet. Experimentally observed degradation of these operational boundaries between nacelle-alone tests and tests on the full configuration were also indicated by the calculation. Application of the program to the General Dynamics STOL fighter design were equally encouraging. Very close agreement was observed between experiment and calculation for the effects of power on pressure distribution, lift and lift curve slope.

Maskew, B.↗

Multigrid method for a vortex breakdown simulation

An inviscid model for a steady axisymmetric flow with swirl was studied. The governing equation is a nonlinear elliptic equation which has more than one solution for a certain range of the swirl parameter. The physically interesting solutions have closed streamlines that look like vortex breakdown (bubble-like solutions). A multigrid method is used to find these solutions. Using an FMG algorithm (nested iteration), the problem is solved in just a few multigrid cycles.

Taasan, S.↗

Computer Program Predicts Turbine-Stage Performance

MTSBL updated version of flow-analysis programs MERIDL and TSONIC coupled to boundary-layer program BLAYER. Method uses quasi-three-dimensional, inviscid, stream-function flow analysis iteratively coupled to calculated losses so changes in losses result in changes in flow distribution. Manner effects both configuration on flow distribution and flow distribution on losses taken into account in prediction of performance of stage. Written in FORTRAN IV.

Boyle, Robert J.↗

A fast, preconditioned conjugate gradient Toeplitz solver

A simple factorization is given of an arbitrary hermitian, positive definite matrix in which the factors are well-conditioned, hermitian, and positive definite. In fact, given knowledge of the extreme eigenvalues of the original matrix A, an optimal improvement can be achieved, making the condition numbers of each of the two factors equal to the square root of the condition number of A. This technique is to applied to the solution of hermitian, positive definite Toeplitz systems. Large linear systems with hermitian, positive definite Toeplitz matrices arise in some signal processing applications. A stable fast algorithm is given for solving these systems that is based on the preconditioned conjugate gradient method. The algorithm exploits Toeplitz structure to reduce the cost of an iteration to O(n log n) by applying the fast Fourier Transform to compute matrix-vector products. Matrix factorization is used as a preconditioner.

Pan, Victor↗

Trajectory optimization for real-time guidance. I - Time-varying LQR on a parallel processor

A key algorithmic element of a real-time trajectory optimization hardware/software implementation, the quadratic program (QP) solver element, is presented. The purpose of the effort is to make nonlinear trajectory optimization fast enough to provide real-time commands during guidance of a vehicle such as an aeromaneuvering orbiter. Many methods of nonlinear programming require the solution of a QP at each iteration. In the trajectory optimization case the QP has a special dynamic programming structure, a LQR-like structure. QP algorithm speed is increased by taking advantage of this special structure and by parallel implementation.

Psiaki, Mark L.↗

Theoretical performance of foil journal bearings

A modified forward iteration approach for the coupled solution of foil bearings is presented. The method is used to predict the steady state theoretical performance of a journal type gas bearing constructed from an inextensible shell supported by an elastic foundation. Bending effects are treated as negligible. Finite element methods are used to predict both the foil deflections and the pressure distribution in the gas film.

Carpino, M.↗

Testing of FTS fingers and interface using a passive compliant robot manipulator

This report deals with testing of a pair of robot fingers designed for the Flight Telerobotic Servicer (FTS) to grasp a cylinder type of Orbital Replaceable Unit (ORU) interface. The report first describes the objectives of the study and then the testbed consisting of a Stewart Platform-based manipulator equipped with a passive compliant platform which also serves as a force/torque sensor. Kinematic analysis is then performed to provide a closed-form solution for the force inverse kinematics and iterative solution for the force forward kinematics using the Newton's Raphson Method. Mathematical expressions are then derived to compute force/torques applied to the FTS fingers during the mating/demating with the interface. The report then presents the three parts of the experimental study on the feasibility and characteristics of the fingers. The first part obtains data of forces applied by the fingers to the interface under various misalignments, the second part determines the maximum allowable capture angles for mating, and the third part processes and interprets the obtained force/torque data.

Nguyen, Charles C.↗

Coupled radiation effects in thermochemical nonequilibrium shock-capturing flowfield calculations

Lunar and Mars return conditions are examined using the LAURA flowfield code and the LORAN radiation code to assess the effect of radiative coupling on axisymmetric thermochemical nonequilibrium flows. Coupling of the two codes is achieved iteratively. Special treatment required to couple radiation in a shock-capturing method is discussed. Results indicate that while coupling effects are generally the same as occur in equilibrium flows, under certain conditions radiation can modify the chemical kinetics of a nonequilibrium flow and thus alter relaxation processes. Coupling effects are found to be small for all cases considered, except for a five meter diameter aerobrake returning from Mars at 13.6 km/sec.

Hartung, Lin C.↗

Remote sensing of earth terrain

In remote sensing, the encountered geophysical media such as agricultural canopy, forest, snow, or ice are inhomogeneous and contain scatters in a random manner. Furthermore, weather conditions such as fog, mist, or snow cover can intervene the electromagnetic observation of the remotely sensed media. In the modelling of such media accounting for the weather effects, a multi-layer random medium model has been developed. The scattering effects of the random media are described by three-dimensional correlation functions with variances and correlation lengths corresponding to the fluctuation strengths and the physical geometry of the inhomogeneities, respectively. With proper consideration of the dyadic Green's function and its singularities, the strong fluctuation theory is used to calculate the effective permittivities which account for the modification of the wave speed and attenuation in the presence of the scatters. The distorted Born approximation is then applied to obtain the correlations of the scattered fields. From the correlation of the scattered field, calculated is the complete set of scattering coefficients for polarimetric radar observation or brightness temperature in passive radiometer applications. In the remote sensing of terrestrial ecosystems, the development of microwave remote sensing technology and the potential of SAR to measure vegetation structure and biomass have increased effort to conduct experimental and theoretical researches on the interactions between microwave and vegetation canopies. The overall objective is to develop inversion algorithms to retrieve biophysical parameters from radar data. In this perspective, theoretical models and experimental data are methodically interconnected in the following manner: Due to the complexity of the interactions involved, all theoretical models have limited domains of validity; the proposed solution is to use theoretical models, which is validated by experiments, to establish the region in which the radar response is most sensitive to the parameters of interest; theoretically simulated data will be used to generate simple invertible models over the region. For applications to the remote sensing of sea ice, the developed theoretical models need to be tested with experimental measurements. With measured ground truth such as ice thickness, temperature, salinity, and structure, input parameters to the theoretical models can be obtained to calculate the polarimetric scattering coefficients for radars or brightness temperature for radiometers and then compare theoretical results with experimental data. Validated models will play an important role in the interpretation and classification of ice in monitoring global ice cover from space borne remote sensors in the future. We present an inversion algorithm based on a recently developed inversion method referred to as the Renormalized Source-Type Integral Equation approach. The objective of this method is to overcome some of the limitations and difficulties of the iterative Born technique. It recasts the inversion, which is nonlinear in nature, in terms of the solution of a set of linear equations; however, the final inversion equation is still nonlinear. The derived inversion equation is an exact equation which sums up the iterative Neuman (or Born) series in a closed form and, thus, is a valid representation even in the case when the Born series diverges; hence, the name Renormalized Source-Type Integral Equation Approach.

Yueh, Herng-Aung↗

Experimental Investigation of Jet Impingement Heat Transfer Using Thermochromic Liquid Crystals

Jet impingement cooling of a hypersonic airfoil leading edge is experimentally investigated using thermochromic liquid crystals (TLCS) to measure surface temperature. The experiment uses computer data acquisition with digital imaging of the TLCs to determine heat transfer coefficients during a transient experiment. The data reduction relies on analysis of a coupled transient conduction - convection heat transfer problem that characterizes the experiment. The recovery temperature of the jet is accounted for by running two experiments with different heating rates, thereby generating a second equation that is used to solve for the recovery temperature. The resulting solution requires a complicated numerical iteration that is handled by a computer. Because the computational data reduction method is complex, special attention is paid to error assessment. The error analysis considers random and systematic errors generated by the instrumentation along with errors generated by the approximate nature of the numerical methods. Results of the error analysis show that the experimentally determined heat transfer coefficients are accurate to within 15%. The error analysis also shows that the recovery temperature data may be in error by more than 50%. The results show that the recovery temperature data is only reliable when the recovery temperature of the jet is greater than 5 C, i.e. the jet velocity is in excess of 100 m/s. Parameters that were investigated include nozzle width, distance from the nozzle exit to the airfoil surface, and jet velocity. Heat transfer data is presented in graphical and tabular forms. An engineering analysis of hypersonic airfoil leading edge cooling is performed using the results from these experiments. Several suggestions for the improvement of the experimental technique are discussed.

Dempsey, Brian Paul↗

Coupled Radiation Effects in Thermochemical Nonequilibrium Shock-Capturing Flowfield Calculations

Lunar and Mars return conditions are examined using the LAURA flow field code and the LORAN radiation code to assess the effect of radiative coupling on axisymmetric thermochemical nonequilibrium flows. Coupling of the two codes is achieved iteratively. Special treatment required to couple radiation in a shock-capturing method is discussed. Results indicate that while coupling effects are generally the same as occur in equilibrium flows, under certain conditions radiation can modify the chemical kinetics of a nonequilibrium flow and thus alter relaxation processes. Coupling effects are found to be small for all cases considered, except for a five meter diameter aerobrake returning from Mars at 13.6 kilometers per second.

Hartung, Lin C.↗

Thermal Model of a Current-Carrying Wire in a Vacuum

A computer program implements a thermal model of an insulated wire carrying electric current and surrounded by a vacuum. The model includes the effects of Joule heating, conduction of heat along the wire, and radiation of heat from the outer surface of the insulation on the wire. The model takes account of the temperature dependences of the thermal and electrical properties of the wire, the emissivity of the insulation, and the possibility that not only can temperature vary along the wire but, in addition, the ends of the wire can be thermally grounded at different temperatures. The resulting second-order differential equation for the steady-state temperature as a function of position along the wire is highly nonlinear. The wire is discretized along its length, and the equation is solved numerically by use of an iterative algorithm that utilizes a multidimensional version of the Newton-Raphson method.

Border, James↗

ISS Double-Gimbaled CMG Subsystem Simulation Using the Agile Development Method

This paper presents an evolutionary approach in simulating a cluster of 4 Control Moment Gyros (CMG) on the International Space Station (ISS) using a common sense approach (the agile development method) for concurrent mathematical modeling and simulation of the CMG subsystem. This simulation is part of Training systems for the 21st Century simulator which will provide training for crew members, instructors, and flight controllers. The basic idea of how the CMGs on the space station are used for its non-propulsive attitude control is briefly explained to set up the context for simulating a CMG subsystem. Next different reference frames and the detailed equations of motion (EOM) for multiple double-gimbal variable-speed control moment gyroscopes (DGVs) are presented. Fixing some of the terms in the EOM becomes the special case EOM for ISS's double-gimbaled fixed speed CMGs. CMG simulation development using the agile development method is presented in which customer's requirements and solutions evolve through iterative analysis, design, coding, unit testing and acceptance testing. At the end of the iteration a set of features implemented in that iteration are demonstrated to the flight controllers thus creating a short feedback loop and helping in creating adaptive development cycles. The unified modeling language (UML) tool is used in illustrating the user stories, class designs and sequence diagrams. This incremental development approach of mathematical modeling and simulating the CMG subsystem involved the development team and the customer early on, thus improving the quality of the working CMG system in each iteration and helping the team to accurately predict the cost, schedule and delivery of the software.

Inampudi, Ravi↗

A preliminary review and evaluation of historical unrest and eruptive activity at Misti volcano, Arequipa, Peru

Since the Spanish arrived in Arequipa in the mid-16th century, 11 potential small-volume eruptions or periods of unrest are reported at Misti. Yet, documentation of these events is scant, often based upon secondary or tertiary sources, and has never been assessed using modern volcanological methods. As a result, Misti?s historical record contains significant uncertainty and iterated errors. We assess these events with modern volcanological knowledge and when possible return to primary published and archival sources. Additionally, we scrutinized historical newspapers, photos, and contemporary travelers? accounts to better understand the volcano?s activity since the early-19th century. Misti?s earliest reported activity in 1542CE illustrates the problems associated with the volcano?s historical record. The 1542CE event first appears in the literature in the late-19th century. Despite some uncertainty of its validity in these reports, it remains listed as a potential eruption. Numerous contemporaneous accounts document a major earthquake in Arequipa on 22 January 1582CE, causing major damage and widespread fear of an eruption. Upon investigation of the volcano, such fear was quickly dispelled. We suggest that Misti?s 1542CE event did not occur, rather its origin is the fear induced by the 1582CE tectonic earthquake and an iterated typographic error in the date. Historical documents also sometimes record when Misti was explicitly not active, such as in a letter from a resident of Arequipa published in 1872CE. Events are also missing from Misti?s reported history, such as a period of enhanced degassing in 1971CE that is recorded in newspaper articles and eye-witness accounts. In summary, we consider three of Misti?s reported events to be false, four to be possible but unlikely, three are likely, and the 1985CE degassing crisis certainly occurred. We further uncovered preliminary evidence of up to four additional degassing crises at the volcano during the 19th and 20th centuries.

Christopher Harpel↗

Drought Monitoring with Drones: A Hundred Fields at a Time

Drought frequency and severity are likely to increase due to global warming. Droughts already have a substantial negative influence on agriculture and the economy and finding ways to reduce their effects could have a monumental impact. Although NASA already has satellites deployed to collect drought data, these satellites are more for global drought indexing than local. To alleviate droughts on a local level, this paper proposes the use of drones to map soil moisture, plant health, and other drought indicators. The proposed drone design is a fixed-wing UAV equipped with a hyperspectral camera, a LiDAR sensor, and an array of weather sensors. These tools will permit it to reliably capture the necessary data to enhance suggestions on improving drought management practices. The collected data from the drones will be deployable in many ways, including for agricultural and non-agricultural applications. The hyperspectral camera has applications in monitoring the health of crops within a field to direct relief measures to the crops most in need. Thermal and LiDAR imaging can be deployed for locating leaks, predicting shortages of water bodies, and determining a field’s water needs. To implement a drought-monitoring drone, the recommended steps include building a prototype drone design that is equipped with the outlined instruments. The prototype drone could then be deployed to collect training data to guide a neural network that would provide interpretations and predictions from the data for users. After the prototype design is iterated upon, it will be ready for deployment and inform water management methods to serve the world in our battle against drought.

climate change↗

Overlapping Conditions in IMPACT

Introduction The IMPACT 1.0 model assumes independence of 119 different medical conditions. That is, it is assumed that one medical event maps to a single medical condition. However, it is often the case that a single medical event may lead to several different concurrent conditions (e.g., a single accident resulting in multiple fractures, chest/abdominal trauma, and sepsis). There is a significant overlap between many of the conditions in ICL 1.0. For example, depression and anxiety are two separate conditions but often co-exist. This overlap may lead to an overestimation in the incidence calculation. In future iterations of the model, we plan to address progression of medical conditions. Methods There are two possible methods to approach progression of a medical condition. The first is to take any medical condition in ICL, calculate its incidence and determine its outcomes no matter what sequelae occur. For example, calculate the incidence of prostatitis and determine its metrics (e.g., loss of crew life, risk of medical evacuation or task time affected) no matter what the sequelae might be (sepsis, hydronephrosis, renal failure, etc.). The second approach is to take any medical condition in the ICL and calculate the outcomes based on transitions to any other condition. For example, calculate the incidence of nephrolithiasis and determine the incidence of each possible transition (e.g., hydronephrosis, renal failure, UTI, pyelonephritis, sepsis) and then determine its metrics. With perfectly informative evidence, the answers should be the same using either method. The question is which approach best reduces the opportunity for overlapping. Results The team identified 31 conditions that were the most likely to transition to a secondary condition such as sepsis and respiratory failure. It was felt that the first approach, as described above, would be the easiest to implement and would have the greatest reduction in overlapping conditions. The identified conditions were primarily infectious conditions, like UTI and pneumonia (which could lead to both respiratory failure and sepsis), toxic inhalations, cardiac arrest, seizures, and trauma. Once these conditions were identified, the conditions were divided between the four clinicians on the team to identify those that were not consistent with the others in way they were approached. There were 13 conditions that were found to be problematic. Some condition definitions needed to be changed to eliminate references to transitions or to update them in cases where the definitions were out of date. Some conditions required new evidence to correct the loss of crew life, and one condition was felt to have so much overlap that it was recommended for elimination. A second effort required modification of the sepsis incidence calculation to prevent double counting of the conditions, such as prostatitis, that could progress to sepsis. Conclusion The team made significant progress in clarifying the method for dealing with transitions from one medical condition to another. In addition, errors were found and corrected in definitions, incidence, and loss of crew life as well as a condition in which reference was made to a CLiFF that was no longer included. These modifications will provide greater fidelity in the IMPACT model and reduce overlap.

Arian Anderson↗