Search NASA⌕ Search

SEARCH · Search NASA

Results for “1-D model”

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 91 records · Page 5

A numerical simulation of the full two-dimensional electrothermal de-icer pad

The ability to predict the time-temperature history of electrothermal de-icer pads is important in the subsequent design of improved and more efficient versions. These de-icer pads are installed near the surface of aircraft components, for the specific purpose of removing accreted ice. The proposed numerical model can incorporate the full 2-D geometry through a section of a region (i.e., section of an airfoil), that current 1-D numerical codes are unable to do. Thus, the effects of irregular layers, curvature, etc., can now be accounted for in the thermal transients. Each layer in the actual geometry is mapped via a body-fitted coordinate transformation into uniform, rectangular computational grids. The relevant heat transfer equations are transformed and discretized. To model the phase change that might occur in any accreted ice, in an enthalpy formulation the phase change equations are likewise transformed and discretized. The code developed was tested against numerous classical numerical solutions, as well as against experimental de-icing data on a UH1H rotor blade obtained from the NASA Lewis Research Center. The excellent comparisons obtained show that this code can be a useful tool in predicting the performance of current de-icer models, as well as in the designing of future models.

Masiulaniec, Konstanty C.↗

Effect of wave-particle interaction on the outflow of ions at high latitudes

The objective of the research is to use the low-frequency wave spectrum measured by the Plasma Wave Instrument on the DE-1 spacecraft to include the wave-particle interaction (WPI) in the different polar wind models that are available at Utah State University. A Monte Carlo technique was used to simulate the ion diffusion in the velocity space due to scattering by the waves. This enabled us to study the effect of WPI on the magnitude as well as the composition of the outflow of the ionospheric ions. In particular, in the first year the 1-D semi-kinetic code developed by Barakat and Schunk (1983) was modified in order to include the effect of the WPI.

Barakat, Abdallah R.↗

Trajectory optimization based on differential inclusion

A method for generating finite-dimensional approximations to the solutions of optimal control problems is introduced. By employing a description of the dynamical system in terms of its attainable sets in favor of using differential equations, the controls are completely eliminated from the system model. Besides reducing the dimensionality of the discretized problem compared to state-of-the-art collocation methods, this approach also alleviates the search for initial guesses from where standard gradient search methods are able to converge. The mechanics of the new method are illustrated on a simple double integrator problem. The performance of the new algorithm is demonstrated on a 1-D rocket ascent problem ('Goddard Problem') in presence of a dynamic pressure constraint.

Seywald, Hans↗

Trajectory optimization based on differential inclusion

A method for generating finite-dimensional approximations to the solutions of optimal control problems is introduced. By employing a description of the dynamical system in terms of its attainable sets in favor of using differential equations, the controls are completely eliminated from the system model. Besides reducing the dimensionality of the discretized problem compared to state-of-the-art collocation methods, this approach also alleviates the search for initial guesses from where standard gradient search methods are able to converge. The mechanics of the new method are illustrated on a simple double integrator problem. The performance of the new algorithm is demonstrated on a 1-D rocket ascent problem (`Goddard Problem') in presence of a dynamic pressure constraint.

Seywald, Hans↗

The 3D rocket combustor acoustics model

The theory and procedures for determining the characteristics of pressure oscillations in rocket engines with prescribed burning rate oscillations are presented. Analyses including radial and hub baffles and absorbers can be performed in one, two, and three dimensions. Pressure and velocity oscillations calculated using this procedure are presented for the SSME to show the influence of baffles and absorbers on the burning rate oscillations required to achieve neutral stability. Comparisons are made between the results obtained utilizing 1-D, 2-D, and 3-D assumptions with regards to capturing the physical phenomena of interest and computational requirements.

Priem, Richard J.↗

A first scramjet study

A variety of related scramjet engine topics are examined. The flow is assumed to be 1-D, the gas is thermally and calorically perfect, and focus is on low hypersonic Mach numbers. The thrust and lift of an exposed half nozzle, which is used on the aerospace plane, is evaluated as well as a fully confined nozzle. A rough estimate of the drag of an aerospace plane is provided. Thermal effects and shock waves are next discussed. A parametric scramjet model is then presented based on the influence coefficient method, which evaluates the dominant scramjet processes. The independent parameters are the ratio of specific heats, a nondimensional heat addition parameter, and four Mach numbers. The total thrust generated by the combustor and nozzle is shown to be independent of the heat release distribution and the combustor exit Mach number, providing thermal choking is avoided. An operating condition for the combustor is found that maximizes the thrust. An alternative condition is explored when this optimum is no longer realistic. This condition provides a favorable pressure gradient and a reasonable area ratio for the combustor. Parametric results based on the model is provided.

Emanuel, George↗

The Navy/NASA Engine Program (NNEP89): A user's manual

An engine simulation computer code called NNEP89 was written to perform 1-D steady state thermodynamic analysis of turbine engine cycles. By using a very flexible method of input, a set of standard components are connected at execution time to simulate almost any turbine engine configuration that the user could imagine. The code was used to simulate a wide range of engine cycles from turboshafts and turboprops to air turborockets and supersonic cruise variable cycle engines. Off design performance is calculated through the use of component performance maps. A chemical equilibrium model is incorporated to adequately predict chemical dissociation as well as model virtually any fuel. NNEP89 is written in standard FORTRAN77 with clear structured programming and extensive internal documentation. The standard FORTRAN77 programming allows it to be installed onto most mainframe computers and workstations without modification. The NNEP89 code was derived from the Navy/NASA Engine program (NNEP). NNEP89 provides many improvements and enhancements to the original NNEP code and incorporates features which make it easier to use for the novice user. This is a comprehensive user's guide for the NNEP89 code.

Plencner, Robert M.↗

High temperature furnace modeling and performance verifications

Analytical, numerical, and experimental studies were performed on two classes of high temperature materials processing sources for their potential use as directional solidification furnaces. The research concentrated on a commercially available high temperature furnace using a zirconia ceramic tube as the heating element and an Arc Furnace based on a tube welder. The first objective was to assemble the zirconia furnace and construct parts needed to successfully perform experiments. The 2nd objective was to evaluate the zirconia furnace performance as a directional solidification furnace element. The 3rd objective was to establish a data base on materials used in the furnace construction, with particular emphasis on emissivities, transmissivities, and absorptivities as functions of wavelength and temperature. A 1-D and 2-D spectral radiation heat transfer model was developed for comparison with standard modeling techniques, and were used to predict wall and crucible temperatures. The 4th objective addressed the development of a SINDA model for the Arc Furnace and was used to design sample holders and to estimate cooling media temperatures for the steady state operation of the furnace. And, the 5th objective addressed the initial performance evaluation of the Arc Furnace and associated equipment for directional solidification. Results of these objectives are presented.

Smith, James E., Jr.↗

Oblique MHD cosmic-ray modified shocks: Two-fluid numerical simulations

We present the first results of time dependent, two-fluid, cosmic-ray (CR) modified, MHD shock simulations. The calculations were carried out with a new numerical code for 1-D ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they appear to have converged to dynamical steady states. We also present results from an initial time dependent simulation which extend the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescale than parallel shocks.

Frank, Adam↗

FD-TD modeling of 2-D dielectric waveguides for propagation and scattering of femtosecond optical solitons

Experimentalists have produced all-optical switches capable of 100-fs responses. To adequately model such switches, nonlinear effects in optical materials (both instantaneous and dispersive) must be included. In principle, the behavior of electromagnetic fields in nonlinear dielectrics can be determined by solving Maxwell's equations subject to the assumption that the electric polarization has a nonlinear relation to the electric field. However, until our previous work, the resulting nonlinear Maxwell's equations have not been solved directly. Rather, approximations have been made that result in a class of generalized nonlinear Schrodinger equations (GNLSE) that solve only for the envelope of the optical pulses. In this paper, we present first-time calculations from the vector nonlinear Maxwell's equations of femtosecond soliton propagation and scattering, including carrier waves, in two-dimensional systems of dielectric waveguides exhibiting the Kerr and Raman quantum effects. We use the finite-difference time-domain (FD-TD) method in an extension of our 1-D work. There, in a fundamental innovation, we treated the linear and nonlinear convolutions for the electric polarization as new dependent variables. By differentiating these convolutions in the time domain, we derived an equivalent system of coupled, nonlinear second-order ODE's. These equations together with Maxwell's equations form the system that is solved to determine the electromagnetic fields in inhomogeneous nonlinear dispersive media. Backstorage in time is limited to only that needed by the time-integration algorithm for the ODE's, rather than that needed to store the time-history of the kernel functions of the convolutions (1000-10,000 time steps). Thus, a 2-D nonlinear optics model from Maxwell's equations is now feasible.

Joseph, Rose↗

Two-photon annihilation of thermal pairs in strong magnetic fields

The annihilation spectrum of pairs with 1-D thermal distributions in the presence of a strong magnetic field is calculated. Numerical analysis of the spectrum are performed for mildly relativistic temperatures and for different angles of emission with respect to field lines. Teragauss magnetic fields are assumed so that conditions are typical of gamma ray burst and pulsar environments. The spectra at each viewing angle reveal asymmetric line profiles that are signatures of the magnetic broadening and red shifting of the line: these asymmetries are more prominent for small viewing angles. Thermal Doppler broadening tends to dominate in the right wing of the line and obscures the magnetic broadening more at high temperatures and smaller viewing angles. This angular dependence of the line asymmetry may prove a valuable diagnostic tool. For low temperatures and magnetic field strengths, useful analytic expressions are presented for the line width, and also for the annihilation spectrum at zero viewing angle. The results presented find application in gamma ray burst and pulsar models, and may prove very helpful in deducing field strengths and temperatures of the emission regions of these objects from line observations made by Compton GRO and future missions.

Baring, Matthew G.↗

Direct simulation of rarefied hypersonic flows

As the capability of the space transportation vehicles (STV's) expand to meet the requirements for future space exploration and utilization, the effects of rarefied hypersonic flows will play a more significant role in defining the aerodynamic and aerothermodynamic performance of STV's. This is particularly true of the low lift/drag aeroassisted STV's where aerobraking occurs at relatively high altitudes and high velocity. Because of the limitations of the continuum description as expressed by the Navier-Stokes equations and the difficulties of solving the Boltzmann equation, the particle of molecular approach has been developed over the last three decades for modeling rarefied gas effects. The direct simulation Monte Carlo (DSMC) method of Bird is the most used method today for simulating rarefied flows. The DSMC method provides a direct physical simulation as opposed to a numerical solution of a set of model equations. This is accomplished by developing phenomenological models of the relevant physical events. The DSMC method accounts for translational, thermal, chemical, and radiative nonequilibrium effects. The general features of the DSMC method, the numerical requirements for obtaining meaningful results, the modeling used to simulate high temperature gas effects, and applications of the method to calculate the flow about an aeroassist flight experiment vehicle (AFE) are reviewed. The AFE simulates a geosynchronous return while entering the Earth's upper atmosphere at approximately 10 km/s. Results obtained using a general 3-D code are presented for the more rarefied portion of the atmospheric encounter (altitudes of 200 to 100 km) emphasizing surface, flowfield, and aerodynamic characteristics of the AFE. Finally, results obtained using axisymmetric and 1-D versions of the code are presented for lower altitude conditions.

Moss, James N.↗

One- and two-dimensional Stirling machine simulation using experimentally generated flow turbulence models

Investigations of one- and two-dimensional (1- or 2-D) simulations of Stirling machines centered around experimental data generated by the U. of Minnesota Mechanical Engineering Test Rig (METR) are covered. This rig was used to investigate oscillating flows about a zero mean with emphasis on laminar/turbulent flow transitions in tubes. The Space Power Demonstrator Engine (SPDE) and in particular, its heater, were the subjects of the simulations. The heater was treated as a 1- or 2-D entity in an otherwise 1-D system. The 2-D flow effects impacted the transient flow predictions in the heater itself but did not have a major impact on overall system performance. Information propagation effects may be a significant issue in the simulation (if not the performance) of high-frequency, high-pressure Stirling machines. This was investigated further by comparing a simulation against an experimentally validated analytic solution for the fluid dynamics of a transmission line. The applicability of the pressure-linking algorithm for compressible flows may be limited by characteristic number (defined as flow path information traverses per cycle); this warrants further study. Lastly the METR was simulated in 1- and 2-D. A two-parameter k-w foldback function turbulence model was developed and tested against a limited set of METR experimental data.

Goldberg, Louis F.↗

Singularity classification as a design tool for multiblock grids

A major stumbling block in interactive design of 3-D multiblock grids is the difficulty of visualizing the design as a whole. One way to make this visualization task easier is to focus, at least in early design stages, on an aspect of the grid which is inherently easy to present graphically, and to conceptualize mentally, namely the nature and location of singularities in the grid. The topological behavior of a multiblock grid design is determined by what happens at its edges and vertices. Only a few of these are in any way exceptional. The exceptional behaviors lie along a singularity graph, which is a 1-D construct embedded in 3-D space. The varieties of singular behavior are limited enough to make useful symbology on a graphics device possible. Furthermore, some forms of block design manipulation that appear appropriate to the early conceptual-modeling phase can be accomplished on this level of abstraction. An overview of a proposed singularity classification scheme and selected examples of corresponding manipulation techniques is presented.

Jones, Alan K.↗

Density measurements in an expanding flow using holographic interferometry

A nonequilibrium expansion of nitrogen through a 2-dimensional nozzle of a reflected shock tunnel is investigated with laser holography, time resolved pressure measurements and emission spectroscopy to characterize the flow field for vibrational relaxation measurements. A 2D computation is used to simulate the nozzle flow field. Synthetic holograms are generated using the computed density profiles and are compared with the experimental holograms. Also, using a quasi 1-D approximation, populations for each vibrational level are computed by (1) solving the master equations and (2) using the Landau-Teller formulation. Computational accuracy of both models is assessed against the experimental data.

Sharma, Surendra P.↗

A simulation study of the performance of the NASA (2,1,6) convolutional code on RFI/burst channels

In an earlier report, the LINKABIT Corporation studied the performance of the (2,1,6) convolutional code on the radio frequency interference (RFI)/burst channel using analytical methods. Using an R(sub 0) analysis, the report concluded that channel interleaving was essential to achieving reliable performance. In this report, Monte Carlo simulation techniques are used to study the performance of the convolutional code on the RFI/burst channel in more depth. The basic system model under consideration is shown. The convolutional code is the NASA standard code with generators g(exp 1) = 1+D(exp 2)+D(exp 3)+D(exp 5)+D(exp 6) and g(exp 2) = 1+D+D(exp 2)+D(exp 3)+D(exp 6) and d(sub free) = 10. The channel interleaver is of the convolutional or periodic type. The binary output of the channel interleaver is transmitted across the channel using binary phase shift keying (BPSK) modulation. The transmitted symbols are corrupted by an RFI/burst channel consisting of a combination of additive white Gaussian noise (AWGN) and RFI pulses. At the receiver, a soft-decision Viterbi decoder with no quantization and variable truncation length is used to decode the deinterleaved sequence.

Perez, Lance C.↗

Directions of Moving Plaids is Biased by Asymmetric Viewing Windows

Directionally selective V1 neurons are tuned to particular spatio-temporal frequencies and respond to local 1-D edge motion. At least some MT neurons however appear to respond to the actual velocity of moving 2D patterns. To better understand how the 1D local motion information available from VI is integrated to derive a 2D velocity signal we investigated human perception of moving plaids, 2-D patterns composed of the sum of two 1-D sinusoidal gratings of different orientations. We measured the effect of the shape of the viewing window on the perceived direction of plaid motion. The plaids were spatially windowed by 2-D Gaussians with unequal standard deviations (sigma 1, sigma 2). Four observers indicated perceived direction by adjusting a pointer. Direction errors were measured as a function of the difference between window orientation and true plaid direction (DELTA THETA) for several grating spatial frequencies (SF = 0.3, 0.6, 1.2 c/d) and window aspect ratios (AR = sigma 1/sigma 2 = 1, 1.4, 2, 4). Observers showed systematic errors in perceived direction (approx. 15 for AR = 4 and SF = 0.6 c/d) that peaked at DELTA THETA approx. 40. The errors increased for increasing aspect ratio and decreased for increasing spatial frequencies (or number of cycles). These results show that despite the unambiguous motion of the plaids, under these conditions human misperceive the motion. These data constrain models of motion integration from V1 to MT and, in particular, are inconsistent with algorithms that use either the Intersection of Constraints rule or cross correlation to compute the perceived direction of motion.

Beutter, B. R.↗

Direction of Moving Plaids is Biased by Asymmetric Viewing Windows

Directionally selective V1 neurons are tuned to particular spatiotemporal frequencies and respond to local 1-D edge motion. At least some MT neurons however appear to respond to the actual velocity of moving 2D patterns (Movshon et al., EBR, 11:117, 1986). To better understand how the 1D local motion information available from V1 is integrated to derive a 2D velocity signal we investigated human perception of moving plaids, 2-D patterns composed of the sum of two 1-D sinusoidal gratings of different orientations. We measured the effect of the shape of the viewing window on the perceived direction of plaid motion. The plaids were spatially windowed by 2-D Gaussians with unequal standard deviations (sigma l, sigma 2). Four observers indicated perceived direction by adjusting a pointer. Direction errors were measured as a function of the difference between window orientation and true plaid direction (delta theta) for several grating spatial frequencies (SF = 0.3, 0.6, 1.2 c/d) and window aspect ratios (AR = sigma 1/sigma 2 = 1, 1.4, 2,4). Observers showed systematic errors in perceived direction (approx. 15 degrees for R = 4 and SF = 0.6 c/d) that peaked at delta theta approximately 40 degrees. The errors increased for increasing aspect ratio and decreased for increasing spatial frequencies (or number of cycles). These results show that despite the unambiguous motion of the plaids, within asymmetric windows, human can systematically misperceive plaid direction. These data constrain models of motion integration within extrastriate cortex and, in particular, are inconsistent with algorithms that use either the Intersection of Constraints rule or cross correlation to compute velocity.

Beutter, Brent Robert↗