Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct numerical simulation”

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,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Predictive and Neural Predictive Control of Uncertain Systems

Accomplishments and future work are:(1) Stability analysis: the work completed includes characterization of stability of receding horizon-based MPC in the setting of LQ paradigm. The current work-in-progress includes analyzing local as well as global stability of the closed-loop system under various nonlinearities; for example, actuator nonlinearities; sensor nonlinearities, and other plant nonlinearities. Actuator nonlinearities include three major types of nonlineaxities: saturation, dead-zone, and (0, 00) sector. (2) Robustness analysis: It is shown that receding horizon parameters such as input and output horizon lengths have direct effect on the robustness of the system. (3) Code development: A matlab code has been developed which can simulate various MPC formulations. The current effort is to generalize the code to include ability to handle all plant types and all MPC types. (4) Improved predictor: It is shown that MPC design using better predictors that can minimize prediction errors. It is shown analytically and numerically that Smith predictor can provide closed-loop stability under GPC operation for plants with dead times where standard optimal predictor fails. (5) Neural network predictors: When neural network is used as predictor it can be shown that neural network predicts the plant output within some finite error bound under certain conditions. Our preliminary study shows that with proper choice of update laws and network architectures such bound can be obtained. However, much work needs to be done to obtain a similar result in general case.

Kelkar, Atul G.↗
NASA NTRS2004

An Inverse Interpolation Method Utilizing In-Flight Strain Measurements for Determining Loads and Structural Response of Aerospace Vehicles

An important and challenging technology aimed at the next generation of aerospace vehicles is that of structural health monitoring. The key problem is to determine accurately, reliably, and in real time the applied loads, stresses, and displacements experienced in flight, with such data establishing an information database for structural health monitoring. The present effort is aimed at developing a finite element-based methodology involving an inverse formulation that employs measured surface strains to recover the applied loads, stresses, and displacements in an aerospace vehicle in real time. The computational procedure uses a standard finite element model (i.e., "direct analysis") of a given airframe, with the subsequent application of the inverse interpolation approach. The inverse interpolation formulation is based on a parametric approximation of the loading and is further constructed through a least-squares minimization of calculated and measured strains. This procedure results in the governing system of linear algebraic equations, providing the unknown coefficients that accurately define the load approximation. Numerical simulations are carried out for problems involving various levels of structural approximation. These include plate-loading examples and an aircraft wing box. Accuracy and computational efficiency of the proposed method are discussed in detail. The experimental validation of the methodology by way of structural testing of an aircraft wing is also discussed.

Shkarayev, S.↗
NASA NTRS1995

Gear Crack Propagation Investigation

Reduced weight is a major design goal in aircraft power transmissions. Some gear designs incorporate thin rims to help meet this goal. Thin rims, however, may lead to bending fatigue cracks. These cracks may propagate through a gear tooth or into the gear rim. A crack that propagates through a tooth would probably not be catastrophic, and ample warning of a failure could be possible. On the other hand, a crack that propagates through the rim would be catastrophic. Such cracks could lead to disengagement of a rotor or propeller from an engine, loss of an aircraft, and fatalities. To help create and validate tools for the gear designer, the NASA Lewis Research Center performed in-house analytical and experimental studies to investigate the effect of rim thickness on gear-tooth crack propagation. Our goal was to determine whether cracks grew through gear teeth (benign failure mode) or through gear rims (catastrophic failure mode) for various rim thicknesses. In addition, we investigated the effect of rim thickness on crack propagation life. A finite-element-based computer program simulated gear-tooth crack propagation. The analysis used principles of linear elastic fracture mechanics, and quarter-point, triangular elements were used at the crack tip to represent the stress singularity. The program had an automated crack propagation option in which cracks were grown numerically via an automated remeshing scheme. Crack-tip stress-intensity factors were estimated to determine crack-propagation direction. Also, various fatigue crack growth models were used to estimate crack-propagation life. Experiments were performed in Lewis' Spur Gear Fatigue Rig to validate predicted crack propagation results. Gears with various backup ratios were tested to validate crack-path predictions. Also, test gears were installed with special crack-propagation gages in the tooth fillet region to measure bending-fatigue crack growth. From both predictions and tests, gears with backup ratios (rim thickness divided by tooth height) of 3.3 and 1.0 produced tooth fractures, whereas a backup ratio of 0.3 produced rim fractures. For a backup ratio of 0.5, the experiments produced rim fractures and the predictions produced both rim and tooth fractures, depending on the initial geometry of the crack. Good correlation between predicted and measured crack growth was achieved when the fatigue crack-closure concept was introduced into the analysis. As the gear rim thickness decreased, the compressive cyclic stress in the gear-tooth fillet region increased. This retarded crack growth and increased the number of crack-propagation cycles to failure.

Source record↗
NASA NTRS2008

Decentralized Control of Sound Radiation using a High-Authority/Low-Authority Control Strategy with Anisotropic Actuators

This paper describes a combined control strategy designed to reduce sound radiation from stiffened aircraft-style panels. The control architecture uses robust active damping in addition to high-authority linear quadratic Gaussian (LQG) control. Active damping is achieved using direct velocity feedback with triangularly shaped anisotropic actuators and point velocity sensors. While active damping is simple and robust, stability is guaranteed at the expense of performance. Therefore the approach is often referred to as low-authority control. In contrast, LQG control strategies can achieve substantial reductions in sound radiation. Unfortunately, the unmodeled interaction between neighboring control units can destabilize decentralized control systems. Numerical simulations show that combining active damping and decentralized LQG control can be beneficial. In particular, augmenting the in-bandwidth damping supplements the performance of the LQG control strategy and reduces the destabilizing interaction between neighboring control units.

Schiller, Noah H.↗
NASA NTRS2011

Current and Future Impact Risks from Small Debris to Operational Satellites

The collision between Iridium 33 and Cosmos 2251 in 2009 signaled the potential onset of the collision cascade effect, commonly known as the "Kessler Syndrome", in the low Earth orbit (LEO) region. Recent numerical simulations have shown that the 10 cm and larger debris population in LEO will continue to increase even with a good implementation of the commonly-adopted mitigation measures. This increase is driven by collisions involving large and massive intacts, i.e., rocket bodies and spacecraft. Therefore, active debris removal (ADR) of large and massive intacts with high collision probabilities has been argued as a direct and effective means to remediate the environment in LEO. The major risk for operational satellites in the environment, however, comes from impacts with debris just above the threshold of the protection shields. In general, these are debris in the millimeter to centimeter size regime. Although impacts by these objects are insufficient to lead to catastrophic breakup of the entire vehicle, the damage is certainly severe enough to cause critical failure of the key instruments or the entire payload. The focus of this paper is to estimate the impact risks from 5 mm and 1 cm debris to active payloads in LEO (1) in the current environment and (2) in the future environment based on different projection scenarios, including ADR. The goal of the study is to quantify the benefits of ADR in reducing debris impact risks to operational satellites.

Liou, Jer-Chyi↗
NASA NTRS2015

Scaling Relation for Occulter Manufacturing Errors

An external occulter is a spacecraft own along the line-of-sight of a space telescope to suppress starlight and enable high-contrast direct imaging of exoplanets. The shape of an external occulter must be specially designed to optimally suppress starlight and deviations from the ideal shape due to manufacturing errors can result loss of suppression in the shadow. Due to the long separation distances and large dimensions involved for a space occulter, laboratory testing is conducted with scaled versions of occulters etched on silicon wafers. Using numerical simulations for a flight Fresnel occulter design, we show how the suppression performance of an occulter mask scales with the available propagation distance for expected random manufacturing defects along the edge of the occulter petal. We derive an analytical model for predicting performance due to such manufacturing defects across the petal edges of an occulter mask and compare this with the numerical simulations. We discuss the scaling of an extended occulter test-bed.

Sirbu, Dan↗
NASA NTRS2015

Aeroacoustic Simulations of a Nose Landing Gear Using FUN3D on Pointwise Unstructured Grids

Numerical simulations have been performed for a partially-dressed, cavity-closed (PDCC) nose landing gear configuration that was tested in the University of Florida's open-jet acoustic facility known as the UFAFF. The unstructured-grid flow solver FUN3D is used to compute the unsteady flow field for this configuration. Mixed-element grids generated using the Pointwise(TradeMark) grid generation software are used for these simulations. Particular care is taken to ensure quality cells and proper resolution in critical areas of interest in an effort to minimize errors introduced by numerical artifacts. A hybrid Reynolds-averaged Navier-Stokes/large eddy simulation (RANS/LES) turbulence model is used for these simulations. Solutions are also presented for a wall function model coupled to the standard turbulence model. Time-averaged and instantaneous solutions obtained on these Pointwise grids are compared with the measured data and previous numerical solutions. The resulting CFD solutions are used as input to a Ffowcs Williams-Hawkings noise propagation code to compute the farfield noise levels in the flyover and sideline directions. The computed noise levels compare well with previous CFD solutions and experimental data.

Vatsa, Veer N.↗
NASA NTRS2017

Conceptual Design of an Electric Sail Technology Demonstration Mission Spacecraft

There is great interest in examining the outer planets of our solar system and Heliopause region (edge of Solar System) and beyond regions of interstellar space by both the Planetary and Heliophysics communities. These needs are well docu-mented in the recent National Academy of Sciences Decadal Surveys. There is significant interest in developing revolutionary propulsion techniques that will enable such Heliopause scientific missions to be completed within 10 to15 years of the launch date. One such enabling propulsion technique commonly known as Electric Sail (E-Sail) propulsion employs positively charged bare wire tethers that extend radially outward from a rotating spacecraft spinning at a rate of one revolution per hour. Around the positively charged bare-wire tethers, a Debye Sheath is created once positive voltage is applied. This sheath stands off of the bare wire tether at a sheath diameter that is proportional to the voltage in the wire coupled with the flux density of solar wind ions within the solar system (or the location of spacecraft in the solar system. The protons that are expended from the sun (solar wind) at 400 to 800 km/sec are electrostatically repelled away from these positively charged Debye sheaths and propulsive thrust is produced via the resulting momentum transfer. The amount of thrust produced is directly proportional to the total wire length. The Marshall Space Flight Center (MSFC) Electric Sail team is currently funded via a two year Phase II NASA Innovative Advanced Concepts (NIAC) awarded in July 2015. The team's current activities are: 1) Developing a Particle in Cell (PIC) numeric engineering model from the experimental data collected at MSFC's Solar Wind Facility on the interaction between simulated solar wind interaction with a charged bare wire that can be applied to a variety of missions, 2) The development of the necessary tether deployers and tethers to enable successful de-ployment of multiple, multi km length bare tethers, 3) Controllability of the space-craft via a voltage bias to steer itself through the solar system to destinations of discovery. These activities once demonstrated analytically, will require a technology demonstration mission (TDM) around the year2020 to demonstrate that all systems work together seamlessly before a Heliophysics Electrostatic Rapid Transit System (HERTS) mission could be initiated. A notional TDM spacecraft that meets the requirements of such a mission will be showcased in this paper.

Wiegmann, Bruce M.↗
NASA NTRS2023

Procedure Automation Rating Matrix

The National Aeronautics and Space Administration (NASA) Advanced Air Mobility (AAM) National Campaign (NC) is researching the means by which future Urban Air Mobility (UAM) aircraft will operate safely in an integrated and scalable airspace architecture. Consistent with this objective, the NASA NC Airspace Procedures team designed a matrix to evaluate UAM instrument flight procedure design, flyability and interoperability of candidate departure, enroute, and approach architectures in live flight or simulation. The Procedure Automation Rating Matrix (PARM) is a multi-dimensional rating scale designed to provide direct feedback from test pilots and operators to airspace procedure designers developing airspace constructs for the integration and scalability of AAM operations in the National Airspace System (NAS). The PARM is assessed using a hierarchical decision tree that guides the operator through a ten-point alpha-numeric rating scale initiated either with or without the use of automation.

National Campaign↗
NASA NTRS1972

Navigation strategy and filter design for solar electric missions

Methods which have been proposed to improve the navigation accuracy for the low-thrust space vehicle include modifications to the standard Sequential- and Batch-type orbit determination procedures and the use of inertial measuring units (IMU) which measures directly the acceleration applied to the vehicle. The navigation accuracy obtained using one of the more promising modifications to the orbit determination procedures is compared with a combined IMU-Standard. The unknown accelerations are approximated as both first-order and second-order Gauss-Markov processes. The comparison is based on numerical results obtained in a study of the navigation requirements of a numerically simulated 152-day low-thrust mission to the asteroid Eros. The results obtained in the simulation indicate that the DMC algorithm will yield a significant improvement over the navigation accuracies achieved with previous estimation algorithms. In addition, the DMC algorithms will yield better navigation accuracies than the IMU-Standard Orbit Determination algorithm, except for extremely precise IMU measurements, i.e., gyroplatform alignment .01 deg and accelerometer signal-to-noise ratio .07. Unless these accuracies are achieved, the IMU navigation accuracies are generally unacceptable.

Tapley, B. D.↗
NASA NTRS1991

Supercritical transition in plane channel flow with spatially periodic perturbations

Laboratory experiments and numerical simulations have been conducted for plane channel flow with a streamwise-periodic array of cylinders. The primary transition in this open flow occurs as a convective rather than absolute instability and leads to traveling-wave packets, which advect out of the system. The ordered secondary state is characteristic of closed flows, in contrast with other open flows where the primary transition often leads directly to turbulence.

Schatz, Michael F.↗
NASA NTRS1989

PRIME: A bottom-up approach to probabilistic rule development

PRIME is a system to be used by an intelligent machine to allow it to operate in an abstract but uncertain (or stochastic) environment. It maintains a model of the effects of the machine's actions in the form of a rule base, which is induced from experience. This bottom-up approach to rule development allows the model to adapt to changes in the environment. Each rule consists of a condition under which the rule is active, an action, the effect of the action on the environment, and an estimate of the probability of this effect occurring. The effect probabilities are used to model the uncertainty in the environment, permitting multiple possible effects for a single action under a particular set of conditions. The objective of the intelligent machine is to satisfy user-specified goals with maximum probability of success. PRIME fulfills this requirement in two ways: it continuously updates the rule base with the most recent information, to ensure the validity of the model; and it generates plans which have the maximum probability of achieving the goals, based on the probability estimates in the rule base. PRIME is composed of three main processes: exploration, generalization, and planning. In exploration, the machine executes various randomly chosen actions, observes the effects on the environment, and updates the rule base accordingly. This process is used to develop the rule base in simulation, as well as to supplement the current knowledge during normal operation. Generalization is the procedure used to induce general rules from experience, which is encoded in the form of specific rules. These general rules extend the machine's knowledge to situations which have not been encountered yet, thereby increasing the capability of the machine to plan effectively. Planning is the process of constructing an optimal sequence of actions to satisfy a goal, using the rule base to predict the effects of these actions and to determine the probability of success of the plan. The rule representation and many other data structures were specifically chosen to maximize the efficiency of these processes. A simulated environment was designed to test the performance of PRIME. The results of experimentation were largely negative. The main problem was that the domain coverage of the rules was inadequate for the number of rules stored in the rule base, due to redundancies in general rules and numerous rules covering ineffective actions. It was determined that a more efficient generalization, and some form of goal-directed exploration, are necessary in order to solve most of the current deficiencies in PRIME.

Miller, Scott A.↗
NASA NTRS2001

Effects of Gravity on Triple Flame Propagation and Stability

Numerical simulations of 2-D triple flames under gravity force have been implemented to identify the effects of gravity on triple flame structure and propagation properties and to understand the mechanisms of instabilities resulting from both heat release and buoyancy effects. A wide range of gravity conditions, heat release and mixing widths for a scalar mixing layer are computed for downward-propagating (in the same direction with the gravity vector) and upward-propagating (in the opposite direction of the gravity vector) triple flames.

Chen, J.-Y.↗
NASA NTRS2002

Development of Reduced-Order Models for Aeroelastic and Flutter Prediction Using the CFL3Dv6.0 Code

A reduced-order model (ROM) is developed for aeroelastic analysis using the CFL3D version 6.0 computational fluid dynamics (CFD) code, recently developed at the NASA Langley Research Center. This latest version of the flow solver includes a deforming mesh capability, a modal structural definition for nonlinear aeroelastic analyses, and a parallelization capability that provides a significant increase in computational efficiency. Flutter results for the AGARD 445.6 Wing computed using CFL3D v6.0 are presented, including discussion of associated computational costs. Modal impulse responses of the unsteady aerodynamic system are then computed using the CFL3Dv6 code and transformed into state-space form. Important numerical issues associated with the computation of the impulse responses are presented. The unsteady aerodynamic state-space ROM is then combined with a state-space model of the structure to create an aeroelastic simulation using the MATLAB/SIMULINK environment. The MATLAB/SIMULINK ROM is used to rapidly compute aeroelastic transients including flutter. The ROM shows excellent agreement with the aeroelastic analyses computed using the CFL3Dv6.0 code directly.

Silva, Walter A.↗
NASA NTRS2005

Formation Flying of Tethered and Nontethered Spacecraft

A paper discusses the effect of the dynamic interaction taking place within a formation composed of a rigid and a deformable vehicle, and presents the concept of two or more tethered spacecraft flying in formation with one or more separated free-flying spacecraft. Although progress toward formation flight of nontethered spacecraft has already been achieved, the document cites potential advantages of tethering, including less consumption of fuel to maintain formation, very high dynamic stability of a rotating tethered formation, and intrinsically passive gravity-gradient stabilization. The document presents a theoretical analysis of the dynamics of a system comprising one free-flying spacecraft and two tethered spacecraft in orbit, as a prototype of more complex systems. The spacecraft are modeled as rigid bodies and the tether as a mass-less spring with structural viscous damping. Included in the analysis is a study of the feasibility of a centralized control system for maintaining a required formation in low Earth orbit. A numerical simulation of a retargeting maneuver is reported to show that even if the additional internal dynamics of the system caused by flexibility is considered, high pointing precision can be achieved if a fictitious rigid frame is used to track the tethered system, and it should be possible to position the spacecraft with centimeter accuracy and to orient the formation within arc seconds of the desired direction also in the presence of low Earth orbit environmental perturbations. The results of the study demonstrate that the concept is feasible in Earth orbit and point the way to further study of these hybrid tethered and free-flying systems for related applications in orbit around other Solar System bodies.

Quadrelli, Marco B.↗
NASA NTRS1983

Modeling of transient disturbances in coronal-streamer configurations

Numerical simulations of the formation and propagation of mass ejection, loop transients in coronal streamers are discussed. The simulations are accomplished with numerical solutions of the single fluid, ideal MHD equations of motion in the meridional plane. The streamer is produced by simulating the relaxation of an initially radial hydrodynamic flow coupled with a dipole magnetic field. The simulated transient then results from an energy release at the base of the streamer. The legs of the loop transient produced remain essentially stationary while the loop expands mainly in the radial direction with velocities of 400 to 750 km s-1. Once the leading edge of the transient has passed out of the lower corona, the initial streamer configuration is restored after 15 to 24 hours. A second energy release 2 hours later than, and with an energy release identical to, the first does not produce a significant coronal disturbance.

Steinolfson, R. S.↗
NASA NTRS1994

Penetration below a convective zone

Two-dimensional numerical simulations are used to investigate how fully compressible nonlinear convection penetrates into a stably stratified zone beneath a stellar convection zone. Estimates are obtained of the extent of penetration as the relative stability S of the stable to the unstable zone is varied over a broad range. The model deals with a perfect gas possessing a constant dynamic viscosity. The dynamics is dominated by downward-directed plumes which can extend far into the stable material and which can lead to the excitation of a broad spectrum of internal gravity waves in the lower stable zone. The convection is highly time dependent, with the close coupling between the lateral swaying of the plumes and the internal gravity waves they generate serving to modulate the strength of the convection. The depth of penetration delta, determined by the position where the time-averaged kinetic flux has its first zero in the stable layer, is controlled by a balance between the kinetic energy carried into the stable layer by the plumes and the buoyancy braking they experience there. A passive scalar is introduced into the unstable layer to evaluate the transport of chemical species downward. Such a tracer is effectively mixed within a few convective overturning times down to a depth of delta within the stable layer. Analytical estimates based on simple scaling laws are used to interpret the variation of delta with S, showing that it first involves an interval of adiabatic penetration if the local Peclet number of the convection exceeds unity, followed by a further thermal adjustment layer, the depths of each interval scaling in turn as S(exp -1) and S(exp -1/4). These estimates are in accord with the penetration results from the simulations.

Hurlburt, Neal E.↗
NASA NTRS1994

Low-Density Nozzle Flow by the Direct Simulation Monte Carlo and Continuum Methods

Two different approaches, the direct simulation Monte Carlo (DSMC) method based on molecular gasdynamics, and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gasdynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental data, in which density and rotational temperature variations along the centerline and at various locations inside a low-density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data, both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity-temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chang-Hong↗