Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple solutions”

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 325 records · Page 18

A Human-in-the Loop Evaluation of a Coordinated Arrival Departure Scheduling Operations for Managing Departure Delays at LaGuardia Airport

LaGuardia (LGA) departure delay was identified by the stakeholders and subject matter experts as a significant bottleneck in the New York metropolitan area. Departure delay at LGA is primarily due to dependency between LGA's arrival and departure runways: LGA departures cannot begin takeoff until arrivals have cleared the runway intersection. If one-in one-out operations are not maintained and a significant arrival-to-departure imbalance occurs, the departure backup can persist through the rest of the day. At NASA Ames Research Center, a solution called "Departure-sensitive Arrival Spacing" (DSAS) was developed to maximize the departure throughput without creating significant delays in the arrival traffic. The concept leverages a Terminal Sequencing and Spacing (TSS) operations that create and manage the arrival schedule to the runway threshold and added an interface enhancement to the traffic manager's timeline to provide the ability to manually adjust inter-arrival spacing to build precise gaps for multiple departures between arrivals. A more complete solution would include a TSS algorithm enhancement that could automatically build these multi-departure gaps. With this set of capabilities, inter-arrival spacing could be controlled for optimal departure throughput. The concept was prototyped in a human-in-the- loop (HITL) simulation environment so that operational requirements such as coordination procedures, timing and magnitude of TSS schedule adjustments, and display features for Tower, TRACON and Traffic Management Unit could be determined. A HITL simulation was conducted in August 2014 to evaluate the concept in terms of feasibility, controller workload impact, and potential benefits. Three conditions were tested, namely a Baseline condition without scheduling, TSS condition that schedules the arrivals to the runway threshold, and TSS+DSAS condition that adjusts the arrival schedule to maximize the departure throughput. The results showed that during high arrival demand period, departure throughput could be incrementally increased under TSS and TSS+DSAS conditions without compromising the arrival throughput. The concept, operational procedures, and summary results were originally published in ATM20151 but detailed results were omitted. This paper expands on the earlier paper to provide the detailed results on throughput, conformance, safety, flight time/distance, etc. that provide extra insights into the feasibility and the potential benefits on the concept.

trajectory-based operations↗

Multiple Equilibria of the Barotropic Vorticity Equation on a Sphere

Multiple states of the barotropic vorticity equation in which the balance is between the first and third terms on the r.h.s. of (1) are given. Solutions of this type were also considered by Charney and Devore (1979) in their discussion of thermally, rather than topographically, forced waves. Whereas in the case of topographic forcing the multiplicity arises from what they called form-drag instability in the wave-zonal flow interaction, in the thermal forcing case the associated instability appears to be the Rossby-wave instability discussed by Lorenz (1972), and the multiple states of the highly truncated model proved to be unstable when more degrees of freedom were added. The multiple statistically steady solutions described are thus novel in that they do not involve form-drag instability, they have stable statistics in calculations with a large number of degrees of freedoms, and they occur on the sphere, with no artificial confinement in a resonant cavity. Furthermore, the solutions are obtained with no external forcing of the zonal flow. The full non-linear equations for two-dimensional non-divergent motion between smooth, rigid boundaries on a sphere were used.

Suarez, M. J.↗

System Simulation by Recursive Feedback: Coupling A Set of Stand-Alone Subsystem Simulations

Recursive feedback is defined and discussed as a framework for development of specific algorithms and procedures that propagate the time-domain solution for a dynamical system simulation consisting of multiple numerically coupled self-contained stand-alone subsystem simulations. A satellite motion example containing three subsystems (other dynamics, attitude dynamics, and aerodynamics) has been defined and constructed using this approach. Conventional solution methods are used in the subsystem simulations. Centralized and distributed versions of coupling structure have been addressed. Numerical results are evaluated by direct comparison with a standard total-system simultaneous-solution approach.

Nixon, Douglas D.↗

System Simulation by Recursive Feedback: Coupling a Set of Stand-Alone Subsystem Simulations

Recursive feedback is defined and discussed as a framework for development of specific algorithms and procedures that propagate the time-domain solution for a dynamical system simulation consisting of multiple numerically coupled, self-contained, stand-alone subsystem simulations. A satellite motion example containing three subsystems (orbit dynamics, attitude dynamics, and aerodynamics) has been defined and constructed using this approach. Conventional solution methods are used in the subsystem simulations. Centralized and distributed versions of coupling structure have been addressed. Numerical results are evaluated by direct comparison with a standard total-system, simultaneous-solution approach.

Nixon, Douglas D.↗

Progress Toward Overset-Grid Moving Body Capability for USM3D Unstructured Flow Solver

A static and dynamic Chimera overset-grid capability is added to an established NASA tetrahedral unstructured parallel Navier-Stokes flow solver, USM3D. Modifications to the solver primarily consist of a few strategic calls to the Donor interpolation Receptor Transaction library (DiRTlib) to facilitate communication of solution information between various grids. The assembly of multiple overlapping grids into a single-zone composite grid is performed by the Structured, Unstructured and Generalized Grid AssembleR (SUGGAR) code. Several test cases are presented to verify the implementation, assess overset-grid solution accuracy and convergence relative to single-grid solutions, and demonstrate the prescribed relative grid motion capability.

Pandyna, Mohagna J.↗

NASA’s Progress Toward Commercial Space Communications—SATCOM Demonstrations and Wideband Multilingual Terminal Development

The National Aeronautics and Space Administration (NASA) Space Communications and Navigation (SCaN) Program is continuing to forge a path toward programmatic and operational transformation enabling commercialization of communications and navigation services to the greatest extent possible for near Earth users. National Space Policy is driving government agencies to incorporate commercial offerings and services to the maximum extent possible. As such, NASA has no plans to build or deploy Tracking and Data Relay Satellites, but will seek to fill NASA mission space-relay capability needs using commercial satellite communications (SATCOM) providers and services. The goal is to support an approach that is both flexible to commercial service constructs and provides continuity of support with current assets as long as required. Progress toward commercialization is advancing in two key areas: (1) awards have been made for demonstrations of SATCOM services, and (2) development of wideband and multilingual user terminals has advanced to the flight demonstration phase. The Communications Services Project (CSP) at Glenn Research Center (GRC), has awarded Funded Space Act Agreements (FSAA) to six industry vendors. Inmarsat Government Inc., Kuiper Government Solutions (KGS) LLC, SES Government Solutions, Space Exploration Technologies, Telesat U.S. Services LLC, and Viasat Incorporated received a combined $278.5 million to complete technology development and in-space demonstrations over the next several years. The end-to-end service capabilities being targeted are based on existing NASA mission operational needs. Accordingly, each company has proposed a technical approach to lower costs, increase flexibility, and improve performance for a broad range of missions. In the current market, industry vendors operate at a range of different frequencies and use variable, sometime proprietary, coding and modulations schemes. As NASA strives for a network comprised of services provided by multiple distinct suppliers, interoperability is a fitting solution to unify the network. Wideband and multi-lingual user terminals are being developed to bridge differences in industry services. Building on ground demonstrations completed in 2021, the Johns Hopkins Applied Physics Lab (APL) will be flight testing a multi-lingual wideband terminal (payload) and demonstrating connectivity to both government and commercial relay services. This paper provides a discussion of the progress toward NASA’s commercial space communications goals and consideration of key challenges and next steps.

space communications↗

NASA's Progress Toward Commercial Space Communications — SATCOM Demonstrations and Wideband Multilingual Terminal Development

The National Aeronautics and Space Administration (NASA) Space Communications and Navigation (SCaN) Program is continuing to forge a path toward programmatic and operational transformation enabling commercialization of communications and navigation services to the greatest extent possible for near Earth users. U.S. National Space Policy is driving government agencies to incorporate commercial offerings and services to the maximum extent possible. As such, NASA has no plans to build or deploy additional Tracking and Data Relay Satellites (TDRS), but will seek to fill NASA mission space-relay capability needs using commercial satellite communications (SATCOM) providers and services. The goal is to support an approach that is both flexible to commercial service constructs and provides continuity of support with current assets as long as required. Progress toward commercial services is advancing in two key areas: (1) awards have been made for demonstrations of SATCOM services, and (2) development of wideband and multilingual user terminals has advanced to the flight demonstration phase. The Communications Services Project (CSP) at Glenn Research Center (GRC), has awarded Funded Space Act Agreements (FSAA) to six industry vendors. Inmarsat Government Inc., Kuiper Government Solutions (KGS) LLC, SES Government Solutions, Space Exploration Technologies, Telesat U.S. Services LLC, and Viasat Incorporated will receive a combined $278.5 million to complete technology development and in-space demonstrations over the next several years. The end-to-end service capabilities being targeted are based on existing NASA mission operational needs. Accordingly, each company has proposed a technical approach to lower costs, increase flexibility, and improve performance for a broad range of missions. In the current market, industry vendors operate at a range of different frequencies and use variable, sometimes proprietary, coding and modulations schemes. As NASA strives for a network comprised of services provided by multiple distinct suppliers, interoperability is a fitting solution to unify the network. Wideband and multi-lingual user terminals are being developed to bridge differences in industry services. Building on ground demonstrations completed in 2021, the Johns Hopkins Applied Physics Lab (APL) will be flight testing a multi-lingual wideband terminal (payload) and demonstrating connectivity to both government and commercial relay services. This paper provides a discussion of the progress toward NASA’s commercial space communications goals and consideration of key challenges and next steps.

commercialization↗

LRAT: Lightning Radiative Transfer

In this report, we extend to cloud physics the work done for single and multiple scattering of electromagnetic waves. We consider the scattering of light, visible or infrared, by a spherical cloud represented by a statistically homogeneous ensemble of configurations of N identical spherical water droplets whose centers are uniformly distributed in its volume V. The ensemble is specified by the average number rho of scatterers in unit volume and by rho f(R) with f(R) as the distribution function for separations R of pairs. The incident light, vector-phi(sub 0) a plane electromagnetic wave with harmonic time dependence, is from outside the cloud. The propagation parameter kappa(sub 0) and the index of refraction eta(sub 0) determine physically the medium outside the distribution of scatterers. We solve the interior problem separately to obtain the bulk parameters for the scatterer equivalent to the ensemble of spherical droplets. With the interior solution or the equivalent medium approach, the multiple scattering problem is reduced to that of an equivalent single scatterer excited from outside illumination. A dispersion relation which determines the bulk propagation parameter K and the bulk index of refraction eta of the cloud is given in terms of the vector equivalent scattering amplitude vector-G and the dyadic scattering amplitude tilde-g of the single object in isolation. Based on this transfer model we will have the ability to consider clouds composed of inhomogeneous distribution of water and/or ice particles and we will be able to take into account particle size distributions within the cloud. We will also be able to study the effects of cloud composition (i.e., particle shape, size, composition, orientation, location) on the polarization of the single or the multiple scattered waves. Finally, this study will provide a new starting point for studying the problem of lightning radiative transfer.

Phanord, Dieudonne D.↗

Compensator improvement for multivariable control systems

A theory and the associated numerical technique are developed for an iterative design improvement of the compensation for linear, time-invariant control systems with multiple inputs and multiple outputs. A strict constraint algorithm is used in obtaining a solution of the specified constraints of the control design. The result of the research effort is the multiple input, multiple output Compensator Improvement Program (CIP). The objective of the Compensator Improvement Program is to modify in an iterative manner the free parameters of the dynamic compensation matrix so that the system satisfies frequency domain specifications. In this exposition, the underlying principles of the multivariable CIP algorithm are presented and the practical utility of the program is illustrated with space vehicle related examples.

Mitchell, J. R.↗

Attitude Estimation Using Modified Rodrigues Parameters

In this paper, a Kalman filter formulation for attitude estimation is derived using the Modified Rodrigues Parameters. The extended Kalman filter uses a gyro-based model for attitude propagation. Two solutions are developed for the sensitivity matrix in the Kalman filter. One is based upon an additive error approach, and the other is based upon a multiplicative error approach. It is shown that the two solutions are in fact equivalent. The Kalman filter is then used to estimate the attitude of a simulated spacecraft. Results indicate that then new algorithm produces accurate attitude estimates by determining actual gyro biases.

Crassidis, John L.↗

Tuning Parameters in Heuristics by Using Design of Experiments Methods

With the growing complexity of today's large scale problems, it has become more difficult to find optimal solutions by using exact mathematical methods. The need to find near-optimal solutions in an acceptable time frame requires heuristic approaches. In many cases, however, most heuristics have several parameters that need to be "tuned" before they can reach good results. The problem then turns into "finding best parameter setting" for the heuristics to solve the problems efficiently and timely. One-Factor-At-a-Time (OFAT) approach for parameter tuning neglects the interactions between parameters. Design of Experiments (DOE) tools can be instead employed to tune the parameters more effectively. In this paper, we seek the best parameter setting for a Genetic Algorithm (GA) to solve the single machine total weighted tardiness problem in which n jobs must be scheduled on a single machine without preemption, and the objective is to minimize the total weighted tardiness. Benchmark instances for the problem are available in the literature. To fine tune the GA parameters in the most efficient way, we compare multiple DOE models including 2-level (2k ) full factorial design, orthogonal array design, central composite design, D-optimal design and signal-to-noise (SIN) ratios. In each DOE method, a mathematical model is created using regression analysis, and solved to obtain the best parameter setting. After verification runs using the tuned parameter setting, the preliminary results for optimal solutions of multiple instances were found efficiently.

Arin, Arif↗

Recursive Branching Simulated Annealing Algorithm

This innovation is a variation of a simulated-annealing optimization algorithm that uses a recursive-branching structure to parallelize the search of a parameter space for the globally optimal solution to an objective. The algorithm has been demonstrated to be more effective at searching a parameter space than traditional simulated-annealing methods for a particular problem of interest, and it can readily be applied to a wide variety of optimization problems, including those with a parameter space having both discrete-value parameters (combinatorial) and continuous-variable parameters. It can take the place of a conventional simulated- annealing, Monte-Carlo, or random- walk algorithm. In a conventional simulated-annealing (SA) algorithm, a starting configuration is randomly selected within the parameter space. The algorithm randomly selects another configuration from the parameter space and evaluates the objective function for that configuration. If the objective function value is better than the previous value, the new configuration is adopted as the new point of interest in the parameter space. If the objective function value is worse than the previous value, the new configuration may be adopted, with a probability determined by a temperature parameter, used in analogy to annealing in metals. As the optimization continues, the region of the parameter space from which new configurations can be selected shrinks, and in conjunction with lowering the annealing temperature (and thus lowering the probability for adopting configurations in parameter space with worse objective functions), the algorithm can converge on the globally optimal configuration. The Recursive Branching Simulated Annealing (RBSA) algorithm shares some features with the SA algorithm, notably including the basic principles that a starting configuration is randomly selected from within the parameter space, the algorithm tests other configurations with the goal of finding the globally optimal solution, and the region from which new configurations can be selected shrinks as the search continues. The key difference between these algorithms is that in the SA algorithm, a single path, or trajectory, is taken in parameter space, from the starting point to the globally optimal solution, while in the RBSA algorithm, many trajectories are taken; by exploring multiple regions of the parameter space simultaneously, the algorithm has been shown to converge on the globally optimal solution about an order of magnitude faster than when using conventional algorithms. Novel features of the RBSA algorithm include: 1. More efficient searching of the parameter space due to the branching structure, in which multiple random configurations are generated and multiple promising regions of the parameter space are explored; 2. The implementation of a trust region for each parameter in the parameter space, which provides a natural way of enforcing upper- and lower-bound constraints on the parameters; and 3. The optional use of a constrained gradient- search optimization, performed on the continuous variables around each branch s configuration in parameter space to improve search efficiency by allowing for fast fine-tuning of the continuous variables within the trust region at that configuration point.

Bolcar, Matthew↗

Random harmonic analysis program, L221 (TEV156). Volume 1: Engineering and usage

A digital computer program capable of calculating steady state solutions for linear second order differential equations due to sinusoidal forcing functions is described. The field of application of the program, the analysis of airplane response and loads due to continuous random air turbulence, is discussed. Optional capabilities including frequency dependent input matrices, feedback damping, gradual gust penetration, multiple excitation forcing functions, and a static elastic solution are described. Program usage and a description of the analysis used are presented.

Miller, R. D.↗

Quantitative molecular characterization of bovine vitreous and lens with non-invasive dynamic light scattering

The non-invasive technique of dynamic light scattering (DLS) was used to quantitatively characterize vitreous and lens structure on a molecular level by measuring the sizes of the predominant particles and mapping the three-dimensional topographic distribution of these structural macromolecules in three spatial dimensions. The results of DLS measurements in five fresh adult bovine eyes were compared to DLS measurements in model solutions of hyaluronan (HA) and collagen (Coll). In the bovine eyes DLS measurements were obtained from excised samples of gel and liquid vitreous and compared to the model solutions. Measurements in whole vitreous were obtained at multiple points posterior to the lens to generate a three-dimensional 'map' of molecular structure. The macromolecule distribution in bovine lens was similarly characterized.In each bovine vitreous (Bo Vit) specimen, DLS predominantly detected two distinct particles, which differed in diffusion properties and hence size. Comparisons with model vitreous solutions demonstrated that these most likely corresponded to the Coll and HA components of vitreous. Three-dimensional mapping of Bo Vit found heterogeneity throughout the vitreous body, with different particle size distributions for Coll and HA at different loci. In contrast, the three-dimensional distribution of lens macromolecules was more homogeneous. Thus, the non-invasive DLS technique can quantitate the average sizes of vitreous and lens macromolecules and map their three-dimensional distribution. This method to assess quantitatively the macromolecular structure of vitreous and lens should be useful for clinical as well as experimental applications in health and disease. Copyright 2001 Academic Press.

Models, Molecular↗

Rao-Blackwellization for Adaptive Gaussian Sum Nonlinear Model Propagation

When dealing with imperfect data and general models of dynamic systems, the best estimate is always sought in the presence of uncertainty or unknown parameters. In many cases, as the first attempt, the Extended Kalman filter (EKF) provides sufficient solutions to handling issues arising from nonlinear and non-Gaussian estimation problems. But these issues may lead unacceptable performance and even divergence. In order to accurately capture the nonlinearities of most real-world dynamic systems, advanced filtering methods have been created to reduce filter divergence while enhancing performance. Approaches, such as Gaussian sum filtering, grid based Bayesian methods and particle filters are well-known examples of advanced methods used to represent and recursively reproduce an approximation to the state probability density function (pdf). Some of these filtering methods were conceptually developed years before their widespread uses were realized. Advanced nonlinear filtering methods currently benefit from the computing advancements in computational speeds, memory, and parallel processing. Grid based methods, multiple-model approaches and Gaussian sum filtering are numerical solutions that take advantage of different state coordinates or multiple-model methods that reduced the amount of approximations used. Choosing an efficient grid is very difficult for multi-dimensional state spaces, and oftentimes expensive computations must be done at each point. For the original Gaussian sum filter, a weighted sum of Gaussian density functions approximates the pdf but suffers at the update step for the individual component weight selections. In order to improve upon the original Gaussian sum filter, Ref. [2] introduces a weight update approach at the filter propagation stage instead of the measurement update stage. This weight update is performed by minimizing the integral square difference between the true forecast pdf and its Gaussian sum approximation. By adaptively updating each component weight during the nonlinear propagation stage an approximation of the true pdf can be successfully reconstructed. Particle filtering (PF) methods have gained popularity recently for solving nonlinear estimation problems due to their straightforward approach and the processing capabilities mentioned above. The basic concept behind PF is to represent any pdf as a set of random samples. As the number of samples increases, they will theoretically converge to the exact, equivalent representation of the desired pdf. When the estimated qth moment is needed, the samples are used for its construction allowing further analysis of the pdf characteristics. However, filter performance deteriorates as the dimension of the state vector increases. To overcome this problem Ref. [5] applies a marginalization technique for PF methods, decreasing complexity of the system to one linear and another nonlinear state estimation problem. The marginalization theory was originally developed by Rao and Blackwell independently. According to Ref. [6] it improves any given estimator under every convex loss function. The improvement comes from calculating a conditional expected value, often involving integrating out a supportive statistic. In other words, Rao-Blackwellization allows for smaller but separate computations to be carried out while reaching the main objective of the estimator. In the case of improving an estimator's variance, any supporting statistic can be removed and its variance determined. Next, any other information that dependents on the supporting statistic is found along with its respective variance. A new approach is developed here by utilizing the strengths of the adaptive Gaussian sum propagation in Ref. [2] and a marginalization approach used for PF methods found in Ref. [7]. In the following sections a modified filtering approach is presented based on a special state-space model within nonlinear systems to reduce the dimensionality of the optimization problem in Ref. [2]. First, the adaptive Gaussian sum propagation is explained and then the new marginalized adaptive Gaussian sum propagation is derived. Finally, an example simulation is presented.

state estimation↗

Enthalpies of formation of CaAl4O7 and CaAl12O19 (hibonite) by high temperature, alkali borate solution calorimetry

Enthalpies of formation were determined for two calcium aluminate phases, CaAl4O7 and CaAl12O19, using high-temperature alkali borate solution calorimetry. The aluminates were synthesized by multiple-cycle heating and grinding stoichiometric mixtures of CaCO3 and Al2O3, and the products were characteized by X-ray diffraction and SEM microbeam analysis. The data on impurities (CaAl4O7 was found to be about 89.00 percent pure by weight and the CaAl12O19 samples about 91.48 percent pure) were used to correct the heat of solution values of the synthetic products. The enthalpies of formation, at 1063 K, from oxides, were found to be equal to -(25.6 + or - 4.7) kJ/g.f.w. for CaAl4O7 and -(33.0 + or - 9.7) kJ/g.f.w. for CaAl12O19; the respective standard enthalpies of formation from elements, at 298 K, were estimated to be -4007 + or - 5.2 kJ/g.f.w. and -10,722 + or - 12 kJ/g.f.w.

Geiger, C. A.↗

A method for determining optimum phasing of a multiphase propulsion system for a single-stage vehicle with linearized inert weight

A general analytical treatment is presented of a single-stage vehicle with multiple propulsion phases. A closed-form solution for the cost and for the performance and a derivation of the optimal phasing of the propulsion are included. Linearized variations in the inert weight elements are included, and the function to be minimized can be selected. The derivation of optimal phasing results in a set of nonlinear algebraic equations for optimal fuel volumes, for which a solution method is outlined. Three specific example cases are analyzed: minimum gross lift-off weight, minimum inert weight, and a minimized general function for a two-phase vehicle. The results for the two-phase vehicle are applied to the dual-fuel rocket. Comparisons with single-fuel vehicles indicate that dual-fuel vehicles can have lower inert weight either by development of a dual-fuel engine or by parallel burning of separate engines from lift-off.

Martin, J. A.↗

COMOC: Three dimensional boundary region variant, programmer's manual

The three-dimensional boundary region variant of the COMOC computer program system solves the partial differential equation system governing certain three-dimensional flows of a viscous, heat conducting, multiple-species, compressible fluid including combustion. The solution is established in physical variables, using a finite element algorithm for the boundary value portion of the problem description in combination with an explicit marching technique for the initial value character. The computational lattice may be arbitrarily nonregular, and boundary condition constraints are readily applied. The theoretical foundation of the algorithm, a detailed description on the construction and operation of the program, and instructions on utilization of the many features of the code are presented.

Orzechowski, J. A.↗