Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 217 records · Page 12

Design, construction and testing of annular diffusers for high speed civil transportation combustor applications

A theoretical and an experimental design study of subsonic flow through curved-wall annular diffusers is being carried out in order to establish the most pertinent design parameters for such devices and the implications of their application in the design of engine components in the aerospace industries. This investigation consists of solving numerically the full Navier Stokes and Continuity equations for the time-mean flow. Various models of turbulence are being evaluated for adoption throughout the study and comparisons would be made with experimental data where they exist. Assessment of diffuser performance based on the dissipated mechanical energy would also be made. The experimental work involves the application of Computer Aided Design software tool to the development of a suitable annular diffuser geometry and the subsequent downloading of such data to a CNC machine at Central State University. The results of the investigations are expected to indicate that more cost effective component design of such devices as effective component design of such devices as diffusers which normally contain complex flows can still be achieved. In this regard a review paper was accepted and presented at the First International Conference on High Speed Civil Transportation Research held at North Carolina A&T in December of 1994.

Okhio, Cyril B.↗

Effects of a Parallel Electric Field and the Geomagnetic Field in the Topside Ionosphere on Auroral and Photoelectron Energy Distributions

The consequences of electric field acceleration and an inhomogencous magnetic field on auroral electron energy distributions in the topside ionosphere are investigated. The one- dimensional, steady state electron transport equation includes elastic and inelastic collisions, an inhomogencous magnetic field, and a field-aligned electric field. The case of a self-consistent polarization electric field is considered first. The self-consistent field is derived by solving the continuity equation for all ions of importance, including diffusion of 0(+) and H(+), and the electron and ion energy equations to derive the electron and ion temperatures. The system of coupled electron transport, continuity, and energy equations is solved numerically. Recognizing observations of parallel electric fields of larger magnitude than the baseline case of the polarization field, the effect of two model fields on the electron distribution function in investigated. In one case the field is increased from the polarization field magnitude at 300 km to a maximum at the upper boundary of 800 km, and in another case a uniform field is added to the polarization field. Substantial perturbations of the low energy portion of the electron flux are produced: an upward directed electric field accelerates the downward directed flux of low-energy secondary electrons and decelerates the upward directed component. Above about 400 km the inhomogencous magnetic field produces anisotropies in the angular distribution of the electron flux. The effects of the perturbed energy distributions on auroral spectral emission features are noted.

Min, Q.-L.↗

A Survey of Symplectic and Collocation Integration Methods for Orbit Propagation

Demands on numerical integration algorithms for astrodynamics applications continue to increase. Common methods, like explicit Runge-Kutta, meet the orbit propagation needs of most scenarios, but more specialized scenarios require new techniques to meet both computational efficiency and accuracy needs. This paper provides an extensive survey on the application of symplectic and collocation methods to astrodynamics. Both of these methods benefit from relatively recent theoretical developments, which improve their applicability to artificial satellite orbit propagation. This paper also details their implementation, with several tests demonstrating their advantages and disadvantages.

integration↗

Validated Numerical Models for the Convective Extinction of Fuel Droplets (CEFD)

The NASA Nebraska Space Grant (NSGC) & EPSCoR programs have continued their effort to support outstanding research endeavors by funding the Numerical Simulation of the Combustion of Fuel Droplets study at the University of Nebraska at Lincoln (UNL). This team of researchers has developed a transient numerical model to study the combustion of suspended and moving droplets. The engines that propel missiles, jets, and many other devices are dependent upon combustion. Therefore, data concerning the combustion of fuel droplets is of immediate relevance to aviation and aeronautical personnel, especially those involved in flight operations. The experiments being conducted by Dr. Gogos and Dr. Nayagam s research teams, allow investigators to gather data for comparison with theoretical predictions of burning rates, flame structures, and extinction conditions. The consequent improved fundamental understanding of droplet combustion may contribute to the clean and safe utilization of fossil fuels (Williams, Dryer, Haggard & Nayagam, 1997, f 2). The present state of knowledge on convective extinction of fuel droplets derives from experiments conducted under normal gravity conditions. However, any data obtained with suspended droplets under normal gravity are grossly affected by gravity. The need to obtain experimental data under microgravity conditions is therefore well justified and addresses one of the goals of NASA's Human Exploration and Development of Space (HEDS) microgravity combustion experiment.

Gogos, George↗

Power Operations of the Mars Exploration Rovers

The rovers of the National Aeronautics and Space Administration’s (NASA) Mars Exploration Rovers (MER) project, Spirit (MER-A) and Opportunity (MER-B), safely landed on the surface of Mars three weeks apart during January 2004. Spirit and Opportunity were built and operated by the Jet Propulsion Laboratory (JPL), which is managed by the California Institute of Technology (Caltech) for NASA. Spirit landed at Gusev Crater, 14.8 degrees south of the equator, and operated continuously on the surface of Mars from January 4, 2004 until last contact from Spirit on March 22, 2010. Opportunity landed at Meridiani Planum, a location 2.5 degrees south of the Martian equator. Opportunity operated continuously on the Martian surface from January 25, 2004 until the last received transmission from Opportunity on June 10, 2018. The goal of the MER project was to determine if Mars ever had a habitable environment, in particular, if it ever had water. During their missions, both Spirit and Opportunity found evidence that liquid water once flowed on the surface of Mars. Both Spirit and Opportunity used a 1.33 m2 triple-junction solar array as their power sources. Based on observations of the original Mars rover, the Sojourner rover of the Mars Pathfinder mission that landed on Mars on July 4, 1997, the expectation was that the Martian dust would rapidly accumulate on the solar arrays of Spirit and Opportunity, and that the rovers would not have enough energy to continue operations after 90 Martian days (sols). Instead, due in part to lower dust accumulation rates than expected and numerous dust cleaning events, the Spirit and Opportunity rovers continued to operate on the Martian surface for over 2000 sols (MER-A) and 5000 sols (MER-B), respectively. During this time, the rovers experienced multiple Martian winters and several dust storms. Because the sources of solar array energy loss were known, the solar array energy output offered a method to scientifically estimate the loading and aeolian removal of dust from the solar arrays each sol. The MER Power operations team called this value that they calculated the solar array Dust Factor (DF). Dust Factor was defined as the fraction of sunlight that penetrates the accumulated dust on the surface of the solar array. A Dust Factor of 1.0 would indicate that the solar array was perfectly clean. A Dust Factor of 0.6 would indicate that only 60% of the available sunlight was able to penetrate the accumulated dust on the solar arrays. The MER Power subsystem operations team used the Multi-Mission Power Analysis Tool (MMPAT) to perform these Dust Factor calculations. The MMPAT software tool modeled the behavior of the solar arrays and the batteries as they interacted with the spacecraft power loads over the mission timeline. MMPAT also had knowledge (through telemetry and user inputs) of telemetered Power subsystem voltages and currents, atmospheric opacity (Tau), rover surface location, rover attitude, the planetary tilt and distance of Mars from the sun based on day of year, the instantaneous elevation of the sun based on time of day, temperatures (internal and external), terrain masking, and shadowing (due to the camera mast and antennas). Once all of the known sources of array energy loss are accounted, the remaining difference between the expected array energy and the actual array energy determines the solar array Dust Factor. The assumptions made while determining the Dust Factor are 1) that the single measured atmospheric opacity value (Tau) is constant over the course of the entire sol, 2) there is no measurable solar cell degradation, 3) there are no shorted strings and therefore, 4) all unexpected solar array energy losses are due to accumulated dust on the solar array. Although it cannot provide an absolute measure of dust loading, the determination of solar array Dust Factor provides a useful way of tracking dust accumulation and aeolian dust removal trends on Martian spacecraft. Any spacecraft on the Martian surface is vulnerable to dust, especially solar-powered spacecraft. The Spirit and Opportunity rovers were operational on the Martian surface for much longer than expected due in part to aeolian removal of dust from their solar arrays. The first few dust removal events were a pleasant surprise to the MER operations teams; however, over time a pattern began to arise. In over three Mars Years on the surface, the Power operations team tracked the solar array Dust Factor at Gusev Crater (location of MER-A) and observed that there were several significant dust removal and deposition events in Mars Year (MY) 27, even in the absence of a large dust storm. In MY 28 at Gusev Crater, the large atmospheric opacity (Tau) increase lagged significant dust removal events. Overall at Gusev Crater, there was a pattern of steady dust accumulation on the solar arrays, with a small number of significant dust cleaning events. At Meridiani Planum, where Opportunity rover operated for over seven Mars Years (late MY 26 to mid MY 34), a clear and consistent pattern of dust movement emerged. Meridiani Planum had a predictable, seasonally dependent pattern of gradual and continuous dust accumulation and removal. In summary, this paper explains the reasons for the development of the solar array Dust Factor and how it was used in mission operations. In particular, this paper describes the MMPAT software package, how it models array energy, including the important assumptions, model inputs and sources of error. And finally, this paper will show how the calculated solar array Dust Factor was used at Gusev Crater and Meridiani Planum to predict dust accumulation rates and weather patterns, and the importance of this generated data set for current and future solar- powered missions to Mars, such as the InSight lander and the planned Mars Sample Return rover.

Chin, Keith B.↗

Numerical simulation of extended corona

A three-dimensional, time-dependent MHD model is presented for the study of coronal dynamics. The model, written in spherical coordinates, extends from the solar surface and was developed with two major issues in mind, namely for interpretation of various steady state and evolutionary dynamical structures in the corona. In order to achieve these objectives, two different numerical techniques are used to seek solutions for these two different, but related, problems; steady state structures and evolutionary structures. These two numerical techniques are: (1) relaxation technique for steady state structures; and (2) FICE (Full-Implicit-Continuous-Eulerian) technique for evolutionary structures. To illustrate this model, numerical results are given for examples of both the steady state and evolutionary structure of the corona which show the additional physical features which cannot be shown by a two-dimensional model.

Wu, S. T.↗

Expert systems and their use in augmenting design optimization

The challenging requirements that are evolving for future aircraft demand that each design be optimally integrated, for the penalties imposed by nonoptimal performance are significant. Classic numerical optimization algorithms have been and will continue to be important tools for aircraft designers. These methods are, however, limited to certain categories of aircraft design variables, leaving the remainder to be determined by the user. A method that makes use of knowledge-based expert systems offers the potential for aiding the conceptual design process in a way that is similar to that of numerical optimization, except that it would address discrete, discontinuous, abstract, or any other unoptimized aspect of vehicle design and integration. Other unique capabilities such as automatic discovery and learning in design may also be achievable in the near term. This paper discusses current practice in conceptual aircraft design and knowledge-based systems, and how knowledge-based systems can be used in conceptual design.

Kidwell, G. H.↗

Hierarchical classification in high dimensional numerous class cases

As progress in new sensor technology continues, increasingly high resolution imaging sensors are being developed. These sensors give more detailed and complex data for each picture element and greatly increase the dimensionality of data over past systems. Three methods for designing a decision tree classifier are discussed: a top down approach, a bottom up approach, and a hybrid approach. Three feature extraction techniques are implemented. Canonical and extended canonical techniques are mainly dependent upon the mean difference between two classes. An autocorrelation technique is dependent upon the correlation differences. The mathematical relationship between sample size, dimensionality, and risk value is derived.

Kim, Byungyong↗

Hierarchical classifier design in high-dimensional, numerous class cases

As progress in new sensor technology continues, increasingly high spectral resolution sensors are being developed. These sensors give more detailed and complex data for each picture element and greatly increase the dimensionality of data over past systems. Three methods for designing a decision tree classifier are discussed; a top down approach, a bottom up approach, and a hybrid approach. Three feature extraction techniques are implemented. Canonical and extended canonical techniques are mainly dependent on the mean difference between two classes. An autocorrelation technique is dependent on the correlation differences. The mathematical relationship among sample size, dimensionality, and risk value is derived.

Kim, Byungyong↗

Analysis of two-dimensional incompressible flow past airfoils using unsteady Navier-Stokes equations

The conservative form of the unsteady Navier-Stokes equations in terms of vorticity and stream function in generalized curvilinear coordinates are used to analyze the flow structure of steady separation and unsteady flow with massive separation. The numerical method solves the discretized equations using an ADI-BGE method. The method is applied to a symmetric 12 percent thick Joukowski airfoil. A conformal clustered grid is generated; several 1-D stretching transformations are used to obtain a grid that attempts to resolve many of the multiple scales of the unsteady flow with massive separation, while maintaining the transformation metrics to be smooth and continuous in the entire flow field. Detailed numerical results are obtained for three flow configurations (1) Re = 1000, alpha = 5 deg, (2) Re = 1000, alpha = 15 deg, (3) Re = 10,000, alpha = 5 deg. No artificial dissipation was added; however, lack of a fine grid in the normal direction has presently led to results which are considered qualitative, especially for case (3).

Ghia, K. N.↗

Analysis of two-dimensional incompressible flow past airfoils using unsteady Navier-Stokes equations

The conservative form of the unsteady Navier-Stokes equations in terms of vorticity and stream function in generalized curvilinear coordinates are used to analyze the flow structure of steady separation and unsteady flow with massive separation. The numerical method solves the discretized equations using an ADI-BGE method. The method is applied to a symmetric 12 percent thick Joukowski airfoil. A conformal clustered grid is generated; several 1-D stretching transformations are used to obtain a grid that attempts to resolve many of the multiple scales of the unsteady flow with massive separation, while maintaining the transformation metrics to be smooth and continuous in the entire flow field. Detailed numerical results are obtained for three flow configurations (1) Re = 1000, alpha = 5 deg., (2) Re =1000, alpha = 15 deg., (3) Re = 10,000, alpha = 5 deg. No artificial dissipation was added; however, lack of a fine grid in the normal direction has presently led to results which are considered qualitative, especially for case (3).

Ghia, K. N.↗

A Global Increase in Nearshore Tropical Cyclone Intensification

Abstract Tropical Cyclones (TCs) inflict substantial coastal damages, making it pertinent to understand changing storm characteristics in the important nearshore region. Past work examined several aspects of TCs relevant for impacts in coastal regions. However, few studies explored nearshore storm intensification and its response to climate change at the global scale. Here, we address this using a suite of observations and numerical model simulations. Over the historical period 1979–2020, observations reveal a global mean TC intensification rate increase of about 3 kt per 24‐hr in regions close to the coast. Analysis of the observed large‐scale environment shows that stronger decreases in vertical wind shear and larger increases in relative humidity relative to the open oceans are responsible. Further, high‐resolution climate model simulations suggest that nearshore TC intensification will continue to rise under global warming. Idealized numerical experiments with an intermediate complexity model reveal that decreasing shear near coastlines, driven by amplified warming in the upper troposphere and changes in heating patterns, is the major pathway for these projected increases in nearshore TC intensification.

54 ENVIRONMENTAL SCIENCES↗

Induced mass and wave motions in the lower solar atmosphere. I - Effects of shear motion of flux tubes

Observations indicate that various dynamic solar phenomena lead to enhanced emission of electromagnetic waves from radio to X-ray wavelengths which can be traced to magnetic activity in the photospheric level. A number of previous investigations have ignored the dynamic responses in the solar atmosphere. On the other hand, Nakagawa et al. (1978, 1981) have studied the atmospheric responses in the frame of MHD in the supersonic super-Alfvenic region. Studies of the slowly varying dynamic response (subsonic) have been unsuccessful because of the requirements of high accuracy in the numerical scheme in which a rigorous mathematical treatment of the boundary conditions is necessary. Recently, a numerical MHD model was constructed by using the full implicit continuous eulerian method. The present investigation makes use of a method which is written in a more convenient numerical code. A two-dimensional, time-dependent, nonplanar MHD model is used to investigate the induced mass and wave motions in the lower solar atmosphere due to the shear motion of flux tubes.

Wu, S. T.↗

Sensitivity Analysis for Solutions to Heterogeneous Nonlocal Systems. Theoretical and Numerical Studies

The paper presents a collection of results on continuous dependence for solutions to nonlocal problems under perturbations of data and system parameters. The integral operators appearing in the systems capture interactions via heterogeneous kernels that exhibit different types of weak singularities, space dependence, even regions of zero-interaction. Here, the stability results showcase explicit bounds involving the measure of the domain and of the interaction collar size, nonlocal Poincaré constant, and other parameters. In the nonlinear setting, the bounds quantify in different L p norms the sensitivity of solutions under different nonlinearity profiles. The results are validated by numerical simulations showcasing discontinuous solutions, varying horizons of interactions, and symmetric and heterogeneous kernels.

97 MATHEMATICS AND COMPUTING↗

Progress Towards Fuselage Drag Reduction via Active Flow Control: A Combined CFD and Experimental Effort

A combined computational and experimental effort has been undertaken to study fuselage drag reduction on a generic, non-proprietary rotorcraft fuselage by the application of active ow control. Fuselage drag reduction is an area of research interest to both the United States and France and this area is being worked collaboratively as a task under the United States/France Memorandum of Agreement on Helicopter Aeromechanics. In the first half of this task, emphasis is placed on the US generic fuselage, the ROBIN-mod7, with the experimental work being conducted on the US side and complementary US and French CFD analysis of the baseline and controlled cases. Fuselage simulations were made using Reynolds-averaged Navier-Stokes ow solvers and with multiple turbulence models. Comparisons were made to experimental data for numerical simulations of the isolated fuselage and for the fuselage as installed in the tunnel, which includes modeling of the tunnel contraction, walls, and support fairing. The numerical simulations show that comparisons to the experimental data are in good agreement when the tunnel and model support are included. The isolated fuselage simulations compare well to each other, however, there is a positive shift in the centerline pressure when compared to the experiment. The computed flow separation locations on the rear ramp region had only slight differences with and without the tunnel walls and model support. For the simulations, the flow control slots were placed at several locations around the flow separation lines as a series of eight slots that formed a nearly continuous U-shape. Results from the numerical simulations resulted in an estimated 35% fuselage drag reduction from a steady blowing flow control configuration and a 26% drag reduction for unsteady zero-net-mass flow control configuration. Simulations with steady blowing show a delayed flow separation at the rear ramp of the fuselage that increases the surface pressure acting on the ramp, thus decreasing the overall fuselage pressure drag.

Schaeffler, Norman W.↗

Linear tracking systems with applications to aircraft control system design

A class of optimal linear time invariant tracking systems, both in continuous time and discrete time, of which the number of inputs (which are restricted to be step functions) is equal to the number of system outputs, is studied. Along with derivation of equations and design procedures, two discretization schemes are presented, constraining either the control or its time derivative, to be a constant over each sampling period. Descriptions are given for the linearized model of the F-8C aircraft longitudinal dynamics, and the C* handling qualities criterion, which then serve as an illustration of the applications of these linear tracking designs. A suboptimal reduced state design is also presented. Numerical results are given for both the continuous time and discrete time designs.

Lee, W. H.↗

Numerical solution of singular Lyapunov equations

We consider the numerical solution of large scale singular (continuous-time) Lyapunov equations of the form AX + XA T + BB T = 0, where A is semistable, that is, its spectrum is contained in the left half plane, with the exception of a few semisimple eigenvalues at zero. We also consider the case of a few semisimple eigenvalues on the imaginary axis. We assume that we know these few eigenvalues (zero or imaginary), and that we have or can compute the corresponding invariant subspaces. We use this information to build an appropriate newly proposed subspace on which to project the Lyapunov equations, and then compute a low-rank approximation to the least squares solution. Selected illustrative numerical examples are provided.

97 MATHEMATICS AND COMPUTING↗

Insights on the continuous representations of piecewise-smooth nonlinear systems: limits of applicability and effectiveness

Dynamical systems subject to intermittent contact are often modeled with piecewise-smooth contact forces. However, the discontinuous nature of the contact can cause inaccuracies in numerical results or failure in numerical solvers. Representing the piecewise contact force with a continuous and smooth function can mitigate these problems, but not all continuous representations may be appropriate for this use. Here, five representations used by previous researchers (polynomial, rational polynomial, hyperbolic tangent, arctangent, and logarithm-arctangent functions) are studied to determine which ones most accurately capture nonlinear behaviors including super- and subharmonic resonances, multiple solutions, and chaos. The test case is a single-DOF forced Duffing oscillator with freeplay nonlinearity, solved using direct time integration. This work intends to expand on past studies by determining the limits of applicability for each representation and what numerical problems may occur.

42 ENGINEERING↗