A modified quasi-linearization method for solving trajectory optimization problems.
Modified quasi-linearization method for numerically solving trajectory optimization problems with undetermined terminal time
LINEARIZATION↗
SEARCH · Engineering Papers
Results for “Boundary value problem”
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.
A method for determining acoustic-liner admittance in ducts with sheared flow in two-cross-sectional directions
A method is developed for determining the acoustic admittance of a test liner installed in the wall of a grazing flow impedance tube. The mean flow is permitted flow gradients in both cross-sectional directions of the tube. The unknown admittance value is obtained by solving an eigenvalue problem. This eigenvalue problem results from the application of the finite-element method to the partial differential equation and boundary conditions governing the acoustic field. The credibility of the method is established by comparing results with exact solutions obtained for a constant mean-flow profile and with previous results for cases involving shear in only one cross-sectional direction. Excellent comparisons are obtained in both cases. The analysis is used in conjunction with a limited amount of experimental data and shows that the flow must be accurately modeled in order to determine the acoustic-liner properties.
Watson, W. R.↗
How good is the impedance boundary condition?
The impedance boundary condition (IBC) is often used in scattering problems involving material-coated conducting bodies. It is shown that for some commonly encountered coating configurations, the value of the impedance varies significantly as functions of the incident angle and polarization. Hence, the use of IBC in a rigorously formulated problem may affect the accuracy of the final solution.
Lee, Shung-Wu↗
Receptivity and Bypass Dynamics
Problems concerning laminar-turbulent transition are addressed by solving a series of initial value problems. Solutions to the temporal, initial-value problem .with an inhomogeneous forcing term imposed upon the flow are sought. It is shown that: (1) A transient disturbance lying located outside of the boundary layer can lead to the growth of an unstable Tollmein-Schlicting wave; (2) A resonance with the continuous spectrum may provide a mechanism for bypass transition; and (3) The continuum modes of a disturbance feed directly into the Tollmein-Schlicting wave downstream through non-parallel effects.
Lasseigne, D. G.↗
Correction of Excessive Precipitation over Steep Mountains in a General Circulation Model (GCM)
Excessive precipitation over steep and high mountains (EPSM) is a well-known problem in GCMs and regional climate models even at a resolution as high as 19km. The affected regions include the Andes, the Himalayas, Sierra Madre, New Guinea and others. This problem also shows up in some data assimilation products. Among the possible causes investigated in this study, we found that the most important one, by far, is a missing upward transport of heat out of the boundary layer due to the vertical circulations forced by the daytime subgrid-scale upslope winds, which in turn is forced by heated boundary layer on the slopes. These upslope winds are associated with large subgrid-scale topographic variance, which is found over steep mountains. Without such subgrid-scale heat ventilation, the resolvable-scale upslope flow in the boundary layer generated by surface sensible heat flux along the mountain slopes is excessive. Such an excessive resolvable-scale upslope flow in the boundary layer combined with the high moisture content in the boundary layer results in excessive moisture transport toward mountaintops, which in turn gives rise to excessive precipitation over the affected regions. We have parameterized the effects of subgrid-scale heated-slope-induced vertical circulation (SHVC) by removing heat from the boundary layer and depositing it in the layers higher up when topographic variance exceeds a critical value. Test results using NASA/Goddard's GEOS-5 GCM have shown that the EPSM problem is largely solved.
Chao, Winston C.↗
A study of three dimensional turbulent boundary layer separation and vortex flow control using the reduced Navier Stokes equations
A reduced Navier-Stokes (RNS) initial value space marching solution technique was applied to vortex generator and separated flow problems and demonstrated good predictions of the engine face flow field. This RNS solution technique using FLARE approximations can adequately describe the topological and topographical structure flow separation associated with vortex liftoff, and this conclusion led to the concept of a subclass of separations which can be called vorticity separations: separations dominated by the transport of vorticity. Adequate near wall resolution of vorticity separations appears necessary for good predictions of these flows.
Anderson, Bernhard H.↗
Numerical Studies of the Non-linear Vlasov Equation
The subject of this dissertation is the numerical integration of the initial-value problem for the non-linear Vlasov equation. The Vlasov equation is used to describe the dynamics of a "collisionless", one-dimensional, classical electron gas confined between two perfectly reflecting boundaries. Only the long-range Coulomb interactions of the electrons are considered; effects associated with the discrete structure are neglected. The numerical results obtained for non-linear Landau damping compare well with similar results obtained by Knorr. A general statement of the results on stable initial conditions is: As the degree of non-linearity of the initial conditions is increased, the deviation from linear Landau damping appears sooner and is more severe. In some cases damping was observed to cease. Curves showing the time dependence of the damping decrement are derived and compared with predictions of non-linear theories. New results obtained in this study include the observation that for strongly non-linear cases, the damping of the electric field causes an initially Maxwellian fo (v, O) to develop a peak in the neighborhood of the phase velocity; strong growth of the second harmonic is seen after fo (v, t) develops such a peak. Also new in this study is the interpretation of the development of a certain class of strongly unstable initial conditions as approaching an inhomogeneous equilibrium.
NUMERICAL INTEGRATION↗
Analysis of High-Speed Aerodynamics of a Swept Wing with Seamless Flaps
Computational fluid dynamics (CFD) analysis was conducted to evaluate the high-speed aerodynamics of a Gulfstream GIII airplane (Gulfstream Aerospace Corporation, Savannah, Georgia) swept wing modified with an experimental seamless, compliant flap called the Adaptive Compliant Trailing Edge (ACTE) flap for airworthiness. The high-speed aerodynamics of the modified ACTE wing were analyzed at a Mach number of 0.85 and an altitude of 40,000 ft. A polyhedral finite-volume unstructured Reynolds-averaged Navier-Stokes CFD code was used with a two-equation turbulence model. Wing aerodynamics sensitivities to lift coefficients, airplane engine operation, vortex generators (VGs), and ACTE flap deflections were examined. High-speed flow over the clean GIII wing without VGs was found to be complex with outboard shock and boundary-layer interactions causing flow separation over the upper wing surface upstream of the aileron, which could result in buffet and controllability problems at higher lift coefficient values. The VGs were found to improve the high-speed aerodynamics of the clean wing, as expected; however, ACTE flap deflections of as little as 1 to 2 deg were found to negate all of the aerodynamic benefits provided by the VGs at high speeds.
wings↗
Iterative solution of nonlinear optimal control problems.
Mathematical programming methods for solving nonlinear state-constrained discrete optimal control problems
DYNAMIC PROGRAMMING↗
Supersonic solutions of the solar-wind equations
We re-examine the inviscid solar-wind equations with heat conduction from below, and establish a fundamentally new approach for finding stellar solutions. Although the problem is fourth order, only two independent integration constants can be assigned, since two boundary conditions that are required to specify well-behaved supersonic solutions determine the values of the other two constants. The solar-wind-model of Noble and Scarf are essentially accurate for practical purposes, but in a fundamental sense they are not self-consistent. At the supersonic level, the ratio of thermal energy to gravitational potential, (kTr/GMm)(sub s), must lie in a narrow range, between 0.4375 and about 0.40, to permit well-behaved solutions for ionized hydrogen. For a supersonic solution this ratio uniquely determines the constants A and epsilon, which in turn govern not only the escape flux and the energy flux per particle but also T'(lambda approaches 0), the temperature gradient at infinity. It is very difficult to find solutions by trial and error, without employing T'(0) as a guide to the iteration of A and epsilon, because not only must (kTr/GMm)(sub s) be chosen within the acceptable range but then either A or epsilon must be guessed with phenomenal accuracy. Otherwise the numerical integration of temperature and expansion velocity do not yield acceptable values at infinity.
Chamberlain, Joseph W.↗
Ocean Energy Systems Wave Energy Modeling Task 10.4: Numerical Modeling of a Fixed Oscillating Water Column
This paper reports on an ongoing international effort to establish guidelines for numerical modeling of wave energy converters, initiated by the International Energy Agency Technology Collaboration Program for Ocean Energy Systems. Initial results for point absorbers were presented in previous work, and here we present results for a breakwater-mounted Oscillating Water Column (OWC) device. The experimental model is at scale 1:4 relative to a full-scale installation in a water depth of 12.8 m. The power-extracting air turbine is modeled by an orifice plate of 1–2% of the internal chamber surface area. Measurements of chamber surface elevation, air flow through the orifice, and pressure difference across the orifice are compared with numerical calculations using both weakly-nonlinear potential flow theory and computational fluid dynamics. Both compressible- and incompressible-flow models are considered, and the effects of air compressibility are found to have a significant influence on the motion of the internal chamber surface. Recommendations are made for reducing uncertainties in future experimental campaigns, which are critical to enable firm conclusions to be drawn about the relative accuracy of the numerical models. It is well-known that boundary element method solutions of the linear potential flow problem (e.g., WAMIT) are singular at infinite frequency when panels are placed directly on the free surface. This is problematic for time-domain solutions where the value of the added mass matrix at infinite frequency is critical, especially for OWC chambers, which are modeled by zero-mass elements on the free surface. A straightforward rational procedure is described to replace ad-hoc solutions to this problem that have been proposed in the literature.
16 TIDAL AND WAVE POWER↗
Spurious frequencies as a result of numerical boundary treatments
The stability theory for finite difference Initial Boundary-Value approximations to systems of hyperbolic partial differential equations states that the exclusion of eigenvalues and generalized eigenvalues is a sufficient condition for stability. The theory, however, does not discuss the nature of numerical approximations in the presence of such eigenvalues. In fact, as was shown previously, for the problem of vortex shedding by a 2-D cylinder in subsonic flow, stating boundary conditions in terms of the primitive (non-characteristic) variables may lead to such eigenvalues, causing perturbations that decay slowly in space and remain periodic time. Characteristic formulation of the boundary conditions avoided this problem. A more systematic study of the behavior of the (linearized) one-dimensional gas dynamic equations under various sets of oscillation-inducing legal boundary conditions is reported.
Abarbanel, Saul↗
Higher-order particle representation for particle-in-cell simulations
In this paper we present an alternative approach to the representation of simulation particles for unstructured electrostatic and electromagnetic PIC simulations. In our modified PIC algorithm we represent particles as having a smooth shape function limited by some specified finite radius, r 0 . A unique feature of our approach is the representation of this shape by surrounding simulation particles with a set of virtual particles with delta shape, with fixed offsets and weights derived from Gaussian quadrature rules and the value of r 0 . As the virtual particles are purely computational, they provide the additional benefit of increasing the arithmetic intensity of traditionally memory bound particle kernels. The modified algorithm is implemented within Sandia National Laboratories' unstructured EMPIRE-PIC code, for electrostatic and electromagnetic simulations, using periodic boundary conditions. We show results for a representative set of benchmark problems, including electron orbit, a transverse electromagnetic wave propagating through a plasma, numerical heating, and a plasma slab expansion. In this work, good error reduction across all of the chosen problems is achieved as the particles are made progressively smoother, with the optimal particle radius appearing to be problem-dependent.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Numerical spatial marching techniques in duct acoustics
Direct calculation of the internal structure of a ducted noise source from farfield pressure measurements is regarded as an initial value problem, where the pressure and pressure gradient (farfield impedance) are assumed to be known along a line in the farfield. If pressure and impedance are known at the boundary of the farfield, the pressure can be uniquely determined in the vicinity of the inlet and inside the inlet ducting. A marching procedure is developed which, with this information obtained from measurements, enables a description of a ducted noise source. The technique uses a finite difference representation of the homogeneous Helmholtz equation.
Baumeister, K. J.↗
Eigenvalues of a baroclinic stability problem with Ekman damping
An analytical solution is presented for the baroclinic stability problem of a Boussinesq fluid in a beta-plane channel with Ekman suction boundary conditions. All of the modes, stable and unstable, belonging to this problem are identified. It is found that an unstable mode exists for only a certain range of values of the Burger number. The value of the Burger number at the upper limit of this range increases as the Ekman number decreases. Beyond this upper limit only a damped mode exists. It is also found that this transition in parameter space from the unstable to the stable mode occurs in a discontinuous manner.
Antar, B. N.↗
The 3DGRAPE book: Theory, users' manual, examples
A users' manual for a new three-dimensional grid generator called 3DGRAPE is presented. The program, written in FORTRAN, is capable of making zonal (blocked) computational grids in or about almost any shape. Grids are generated by the solution of Poisson's differential equations in three dimensions. The program automatically finds its own values for inhomogeneous terms which give near-orthogonality and controlled grid cell height at boundaries. Grids generated by 3DGRAPE have been applied to both viscous and inviscid aerodynamic problems, and to problems in other fluid-dynamic areas. The smoothness for which elliptic methods are known is seen here, including smoothness across zonal boundaries. An introduction giving the history, motivation, capabilities, and philosophy of 3DGRAPE is presented first. Then follows a chapter on the program itself. The input is then described in detail. A chapter on reading the output and debugging follows. Three examples are then described, including sample input data and plots of output. Last is a chapter on the theoretical development of the method.
Sorenson, Reese L.↗
Correction of Excessive Precipitation Over Steep and High Mountains in a General Circulation Model
Excessive precipitation over steep and high mountains (EPSM) is a well-known problem in GCMs and meso-scale models. This problem impairs simulation and data assimilation products. Among the possible causes investigated in this study, we found that the most important one, by far, is a missing upward transport of heat out of the boundary layer due to the vertical circulations forced by the daytime upslope winds, which are forced by the heated boundary layer on subgrid-scale slopes. These upslope winds are associated with large subgrid-scale topographic variation, which is found over steep and high mountains. Without such subgridscale heat ventilation, the resolvable-scale upslope flow in the boundary layer generated by surface sensible heat flux along the mountain slopes is excessive. Such an excessive resolvablescale upslope flow combined with the high moisture content in the boundary layer results in excessive moisture transport toward mountaintops, which in turn gives rise to EPSM. Other possible causes of EPSM that we have investigated include 1) a poorly-designed horizontal moisture flux in the terrain-following coordinates, 2) the condition for cumulus convection being too easily satisfied at mountaintops, 3) the presence of conditional instability of the computational kind, and 4) the absence of blocked flow drag. These are all minor or inconsequential. We have parameterized the ventilation effects of the subgrid-scale heated-slope-induced vertical circulation (SHVC) by removing heat from the boundary layer and depositing it in layers higher up when the topographic variance exceeds a critical value. Test results using NASA/Goddard's GEOS-S GCM have shown that this largely solved the EPSM problem.
Chao, Winston C.↗
Boundary effects on Marangoni-Rayleigh and morphological instabilities
Effects of boundaries (vertical sides) on the Marangoni-Rayleigh problem are reviewed. Spectral crowding and a dimensional singularity when Rayleigh effects are neglected are discussed. In the unbounded case, there is a unique Ma with unique eigenfunction for fixed Ra but infinite Ra for fixed Ma. These infinite values of Ra correspond to vertical cellular modes. In the bounded case, there is infinite Ma for a fixed Ra. The stress free side wall case gives double infinite of Ra for fixed Ma and wave number. Clearly, a dimensional singularity occurs for the Marangoni problem when the side walls go to infinity and a fixed Ra is considered. There is no direct analogy regarding this singularity in the morphological instability problem. Results on the behavior of Mac, Rac with boundary conditions also show only marginal connection to the morphological stability problem.
Narayanan, R.↗