Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical techniques”

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,441 records · Page 80

Boundary shape identification problems in two-dimensional domains related to thermal testing of materials

The identification of the geometrical structure of the system boundary for a two-dimensional diffusion system is reported. The domain identification problem treated here is converted into an optimization problem based on a fit-to-data criterion and theoretical convergence results for approximate identification techniques are discussed. Ressults of numerical experiments to demonstrate the efficacy of the theoretical ideas are reported.

Banks, H. T.↗

Spike: Artificial intelligence scheduling for Hubble space telescope

Efficient utilization of spacecraft resources is essential, but the accompanying scheduling problems are often computationally intractable and are difficult to approximate because of the presence of numerous interacting constraints. Artificial intelligence techniques were applied to the scheduling of the NASA/ESA Hubble Space Telescope (HST). This presents a particularly challenging problem since a yearlong observing program can contain some tens of thousands of exposures which are subject to a large number of scientific, operational, spacecraft, and environmental constraints. New techniques were developed for machine reasoning about scheduling constraints and goals, especially in cases where uncertainty is an important scheduling consideration and where resolving conflicts among conflicting preferences is essential. These technique were utilized in a set of workstation based scheduling tools (Spike) for HST. Graphical displays of activities, constraints, and schedules are an important feature of the system. High level scheduling strategies using both rule based and neural network approaches were developed. While the specific constraints implemented are those most relevant to HST, the framework developed is far more general and could easily handle other kinds of scheduling problems. The concept and implementation of the Spike system are described along with some experiments in adapting Spike to other spacecraft scheduling domains.

Johnston, Mark↗

On estimating the phase of a periodic waveform in additive Gaussian noise, part 3

Motivated by advances in signal processing technology that support more complex algorithms, researchers have taken a new look at the problem of estimating the phase and other parameters of a nearly periodic waveform in additive Gaussian noise, based on observation during a given time interval. Parts 1 and 2 are very briefly reviewed. In part 3, the actual performances of some of the highly nonlinear estimation algorithms of parts 1 and 2 are evaluated by numerical simulation using Monte Carlo techniques.

Rauch, L. L.↗

A multi-block multigrid method for the solution of the Euler and Navier-Stokes equations for three-dimensional flows

A multi-block multigrid method for the solution of the three-dimensional Euler and Navier-Stokes equations is presented. The basic flow solver is a cell-vertex method which employs central-difference spatial approximations and Runge-Kutta time stepping. The use of local time stepping, implicit residual smoothing, multigrid techniques, and variable-coefficient numerical smoothing results in an efficient and robust scheme. The multi-block strategy places the block loop within the Runge-Kutta loop such that accuracy and convergence are not affected by block boundaries. This has been verified by comparing the results of one- and two-block calculations in which the two-block grid is generated by splitting the one-block grid. Results are presented for both Euler and Navier-Stokes computations of wings and wing-fuselage combinations.

Atkins, H. L.↗

Propagation characteristics of optical fiber structures with arbitrary shape and index variation

The application of the scalar wave-fast Fourier transform (SW-FFT) technique to the computation of the propagation characteristics of some complex optical fiber structures is presented. The SW-FFT technique is based on the numerical solution of the scalar wave equation by a forward-marching fast Fourier transform method. This solution yields the spatial configuration of the fields as well as its modal characteristics in and around the guiding structure. The following are treated by the SW-FFT method: analysis of coupled optical fibers and computation of their odd and even modes and coupling length; the solution of tapered optical waveguides (transitions) and the study of the effect of the slope of the taper on mode conversion; and the analysis of branching optical fibers and demonstration of their mode-filtering and/or power-dividing properties.

Manshadi, F.↗

A two-dimensional Euler solution for an unbladed jet engine configuration

A two dimensional nonaxisymmetric Euler solution in a geometric representation of a jet engine configuration without blades is presented. The domain, including internal and external flow, is covered with a multiblock grid. To construct the grid, a domain decomposition technique is used to subdivide the domain and smooth grids are dimensioned and placed in each block. The grid contains 44 blocks which cover the external field, the inlet, bypass duct, core duct and nozzle of the nonaxisymmetric engine configuration. The geometry is symmetric about the meanline of the hub, but the grid is not since there is no symmetry condition applied to the grid between the two halves. With a symmetric grid at zero angle of attack, the measures of the solution would cancel exactly. With an asymmetric grid, the solution will not necessarily be symmetric and the lift coefficient will not necessarily be zero. Thus, grid asymmetry can be exploited to verify the resolution of the solution. The solution may be verified on the basis of five theoretical quantities: conservation of mass and energy, deviation of the lift coefficient from zero, deviation of the drag coefficient from zero, deviations from constant entropy, and deviations in the pressure distributions over the symmetric surfaces of the components. This technique is suitable for obtaining numerical solutions in complex geometries and provides a foundation for complete engine throughflow calculations.

Stewart, Mark E. M.↗

Some further observations on droplet combustion characteristics - NASA LeRC-Princeton results

Experimental and numerical studies are reviewed which are designed to examine the effects of droplet/gas motion, product-intermediate absorption, extinction, and sooting on droplet combustion. The experimental work at the NASA-Lewis Research Center involves a 2.2-s droptower for investigating microgravitational effects of droplet combustion over a relatively extended range. The droplet-gas velocities are very low because the spherosymmetrical nature of the major combustion processes produces a quiescent environment. The refined experimental results are combined with numerical modeling based on a technique that is fully transient, comprehensive, and has few empirical simplifications. The combination of techniques improves the present understanding of convection-induced effects, reducing soot formation, and promoting quiescent droplet combustion.

Choi, Mun Y.↗

Rough surface scattering from moving ocean surfaces as an indicator of scattering mechanisms

Ever since Crombie used the Doppler spectrum of RF signals scattered from the ocean surface to provide conclusive evidence for Bragg resonant scattering, scattering from the ocean surface has served as a testing ground for the identification of governing scattering mechanisms in rough surface scattering. The moving ocean surface has the advantage for rough surface scattering that, being a dispersive medium, features with different spatial scales the tendency to travel at different speeds. Thus, by looking at the temporal signature of the scattered field, one can try to isolate the scattering mechanism. This feature has been exploited by various authors. To give a recent example, we cite the work of Plant, who examined the Doppler spectrum of the return signal to infer the presence of two-scale scattering in scattering at medium incidence angles. The problem with restricting oneself to the examination of the Doppler spectrum, as has been done traditionally, is that one restricts the phenomena under investigation to temporally homogeneous phenomena: transient phenomena, which are localized in time, are not localized in frequency space. There is increasing evidence that this type of phenomena, due to specular scattering or breaking waves, may also play an important role in determining the scattering mechanism in ocean-like surfaces. To overcome this problem, we introduce the use of the Wavelet Transform to study the frequency-temporal signatures of the scattered field from moving ocean waves. We calculate this signature using various analytic scattering theories and show that the Wavelet Transform provides a useful tool for separating the different scattering mechanisms operating in scattering in ocean-like surfaces. Next we simulate realistic nonlinear moving ocean waves and calculate the temporal scattered field signature by using the method of moments and the stop-start approximation: the surface is assumed stationary during a scattering event, but moves between incident EM pulses. We compute and compare the Doppler frequency spectrum and the Wavelet Transform of the resulting signal for various types of surfaces and show how different scattering mechanisms may be isolated by this technique. Finally, we compare these numerical results against the analytical predictions in order to obtain the domain of applicability of various scattering theories for moving ocean surfaces which may include wave breaking and specular glints.

Rodriguez, Ernesto↗

High Resolution Euler Solvers Based on the Space-Time Conservation Element and Solution Element Method

The I-D, quasi I-D and 2-D Euler solvers based on the method of space-time conservation element and solution element are used to simulate various flow phenomena including shock waves, Mach stem, contact surface, expansion waves, and their intersections and reflections. Seven test problems are solved to demonstrate the capability of this method for handling unsteady compressible flows in various configurations. Numerical results so obtained are compared with exact solutions and/or numerical solutions obtained by schemes based on other established computational techniques. Comparisons show that the present Euler solvers can generate highly accurate numerical solutions to complex flow problems in a straightforward manner without using any ad hoc techniques in the scheme.

Wang, Xiao-Yen↗

Preconditioning the Helmholtz Equation for Rigid Ducts

An innovative hyperbolic preconditioning technique is developed for the numerical solution of the Helmholtz equation which governs acoustic propagation in ducts. Two pseudo-time parameters are used to produce an explicit iterative finite difference scheme. This scheme eliminates the large matrix storage requirements normally associated with numerical solutions to the Helmholtz equation. The solution procedure is very fast when compared to other transient and steady methods. Optimization and an error analysis of the preconditioning factors are present. For validation, the method is applied to sound propagation in a 2D semi-infinite hard wall duct.

Baumeister, Kenneth J.↗

Optical Sensor/Actuator Locations for Active Structural Acoustic Control

Researchers at NASA Langley Research Center have extensive experience using active structural acoustic control (ASAC) for aircraft interior noise reduction. One aspect of ASAC involves the selection of optimum locations for microphone sensors and force actuators. This paper explains the importance of sensor/actuator selection, reviews optimization techniques, and summarizes experimental and numerical results. Three combinatorial optimization problems are described. Two involve the determination of the number and position of piezoelectric actuators, and the other involves the determination of the number and location of the sensors. For each case, a solution method is suggested, and typical results are examined. The first case, a simplified problem with simulated data, is used to illustrate the method. The second and third cases are more representative of the potential of the method and use measured data. The three case studies and laboratory test results establish the usefulness of the numerical methods.

Piezoelectric actuators↗

Experiment 13: The Study of Dopant Segregation Behavior During the Growth of GaAs in Microgravity on USML-2

An investigation into the segregation behavior of selenium doped gallium arsenide (Se/GaAs) during directional solidification in the microgravity environment was conducted using the Crystal Growth Furnace (CGF) aboard the second United States Microgravity Laboratory (USML-2). Two crystals were successfully processed on USML-2, which lasted from October 20 to November 7, 1995. The first sample was processed for 67 hours, 45 minutes (MET 5/04:53:45-8/00:23:50) and included 19 hours of growth at 0.5 microns/sec which yielded 3.42 cm of sample length, and 5 hours of growth at 1.5 microns/sec which yielded 2.7 cm of sample. During the second experiment, the furnace temperature was adjusted to move the melt-solid interface position towards the hot end of the furnace. The second sample was processed for 50 hours, 10 minutes (MET 8/18:48:49-10/21:58:54) and included 11 hours of growth at 0.5 microns/sec which yielded 1.98 cm of sample, and 1 hour, 25 minutes of growth at 5.0 microns/sec which yielded 2.6 cm of sample. This sample provides an order of magnitude change in growth rate and reproduces one of the growth rates used during USML-1. In contrast to the results from USML-1, no voids were present in either crystal grown on USML-2. The absence of voids in either sample indicates that growth rate changes alone were not responsible for the formation of voids found in the crystals grown on USML-1. Sections of the ground-based and flight crystals grown on USML-2 were cut and polished. All of the interface demarcation lines expected from the current pulse interface demarcation (CPID) system have been identified. These measurements have been analyzed for interface positions, interface shapes, and growth rates. Using a newly developed technique, based on experimental and numerical results, the seeding interface reproducibility from run to run was <= 2.5 mm. The seeding interface position could be controllably moved, with respect to the furnace zones, by adjusting the control set points of the heating zones. The interface shapes flattened slightly as the interface position moved closer to the hot zone but was always an unfavorable concave into the solid shape. The growth rate was found to equal the furnace translation rate, after a 2 -hour transient, for growth rates <= 1.0 microns/sec. Segregation measurements for the ground-based crystals are indicative of complete mixing behavior, as expected. Segregation measurements of the flight crystals are still in progress.

Matthiesen, David H.↗

Synthesis of Optimal Constant-Gain Positive-Real Controllers for Passive Systems

This paper presents synthesis methods for the design of constant-gain positive real controllers for passive systems. The results presented in this paper, in conjunction with the previous work by the authors on passification of non-passive systems, offer a useful synthesis tool for the design of passivity-based robust controllers for non-passive systems as well. Two synthesis approaches are given for minimizing an LQ-type performance index, resulting in optimal controller gains. Two separate algorithms, one for each of these approaches, are given. The synthesis techniques are demonstrated using two numerical examples: control of a flexible structure and longitudinal control of a fighter aircraft.

Mao, Y.↗

Analysis of Cirrus Cloud Microphysical Data

The First International Satellite Cloud Climatology Regional Experiment (FIRE) program has the goal of improving our capabilities to understand, model and detect the properties of climatically-important clouds. This is being undertaken through a three-pronged effort of modeling, long-term observations and short-term intensive field studies. Through examination of satellite and other data it is apparent that stratus and cirrus cloud types have the greatest impact on climate due to their radiative effects and ubiquitous nature. As a result, the FIRE program has developed two paths of investigation, each having its own subset of research objectives and measurement programs. The work conducted under this grant was directed toward furthering our understanding of cirrus cloud systems. While it is known that cirrus are climatically important, the magnitude and even sign of the impact is unclear. Cirrus clouds affect the transfer of radiation according to their physical depth and location in the atmosphere and their microphysical composition. However, significant uncertainties still exist in how cirrus clouds form and how they are maintained, what their physical properties are and how they can be parameterized in numerical models. Better remote sensing techniques for monitoring cirrus cloud systems and improved modeling of radiative transfer through ice particles are also needed. A critical element in resolving these issues is a better understanding of cirrus cloud microphysical properties and how they vary. The focus of the research to be conducted under this grant was th data collected in situ by the University of North Dakota Citation aircraft. The goals of this research were to add to the body of knowledge of cirrus cloud microphysics, particularly at the small end of the size spectrum; and analyze the spatial variation of cirrus clouds.

Poellot, Michael R.↗

Decoupling the Role of Particle Inertia and Gravity on Particle Dispersion

Particle dispersion and the influence that particle momentum exchange has on the properties of a turbulent carrier flow in micro-gravity environments challenge present understanding and predictive schemes. The objective of this effort has been to develop and assess high-fidelity simulation tools for predicting particle transport within micro-gravity environments suspended in turbulent flows. The computational technique is based on Direct Numerical Simulation (DNS) of the incompressible Navier-Stokes equations. The particular focus of the present work is on the class of dilute flows in which particle volume fractions and inter-particle collisions are negligible. Particle motion is assumed to be governed by drag with particle relaxation times ranging from the Kolmogorov scale to the Eulerian timescale of the turbulence and particle mass loadings up to one. The velocity field was made statistically stationary by forcing the low wavenumbers of the flow. The calculations were performed using 96(exp 3) collocation points and the Taylor-scale Reynolds number for the stationary flow was 62. The effect of particles on the turbulence was included in the Navier-Stokes equations using the point-force approximation in which 96(exp 3) particles were used in the calculations. DNS results show that particles increasingly dissipate fluid kinetic energy with increased loading, with the reduction in kinetic energy being relatively independent of the particle relaxation time. Viscous dissipation in the fluid decreases with increased loading and is larger for particles with smaller relaxation times. Fluid energy spectra show that there is a non-uniform distortion of the turbulence with a relative increase in small-scale energy. The non-uniform distortion significantly affects the transport of the dissipation rate, with the production and destruction of dissipation exhibiting completely different behaviors. The spectrum of the fluid-particle energy exchange rate shows that the fluid drags particles at low wavenumbers while the converse is true at high wavenumbers for small particles. A spectral analysis shows that the increase of the high wavenumber portion of the fluid energy spectrum can be attributed to transfer of the fluid-particle covariance by the fluid turbulence. This in turn explains the relative increase of small-scale energy caused by small particles observed in the present simulations as well as those of others.

Squires, Kyle D.↗

A De-centralized Scheduling and Load Balancing Algorithm for Heterogeneous Grid Environments

In the past two decades, numerous scheduling and load balancing techniques have been proposed for locally distributed multiprocessor systems. However, they all suffer from significant deficiencies when extended to a Grid environment: some use a centralized approach that renders the algorithm unscalable, while others assume the overhead involved in searching for appropriate resources to be negligible. Furthermore, classical scheduling algorithms do not consider a Grid node to be N-resource rich and merely work towards maximizing the utilization of one of the resources. In this paper, we propose a new scheduling and load balancing algorithm for a generalized Grid model of N-resource nodes that not only takes into account the node and network heterogeneity, but also considers the overhead involved in coordinating among the nodes. Our algorithm is decentralized, scalable, and overlaps the node coordination time with that of the actual processing of ready jobs, thus saving valuable clock cycles needed for making decisions. The proposed algorithm is studied by conducting simulations using the Message Passing Interface (MPI) paradigm.

Arora, Manish↗

A De-Centralized Scheduling and Load Balancing Algorithm for Heterogeneous Grid Environments

In the past two decades, numerous scheduling and load balancing techniques have been proposed for locally distributed multiprocessor systems. However, they all suffer from significant deficiencies when extended to a Grid environment: some use a centralized approach that renders the algorithm unscalable, while others assume the overhead involved in searching for appropriate resources to be negligible. Furthermore, classical scheduling algorithms do not consider a Grid node to be N-resource rich and merely work towards maximizing the utilization of one of the resources. In this paper we propose a new scheduling and load balancing algorithm for a generalized Grid model of N-resource nodes that not only takes into account the node and network heterogeneity, but also considers the overhead involved in coordinating among the nodes. Our algorithm is de-centralized, scalable, and overlaps the node coordination time of the actual processing of ready jobs, thus saving valuable clock cycles needed for making decisions. The proposed algorithm is studied by conducting simulations using the Message Passing Interface (MPI) paradigm.

Arora, Manish↗

Synthesis of Optimal Constant-Gain Positive-Real Controllers for Passive Systems

This paper presents synthesis methods for the design of constant-gain positive real controllers for passive systems. The results presented in this paper, in conjunction with the previous work by the authors on passification of non-passive systems, offer a useful synthesis tool for the design of passivity-based robust controllers for non-passive systems as well. Two synthesis approaches are given for minimizing an LQ-type performance index, resulting in optimal controller gains. Two separate algorithms, one for each of these approaches, are given. The synthesis techniques are demonstrated using two numerical examples: control of a flexible structure and longitudinal control of a fighter aircraft.

Mao, Y.↗