SEARCH · Search NASA

Results for “Nonlinearly Consistent Implementation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Best Practices for Wake Model and Optimization Algorithm Selection in Wind Farm Layout Optimization

This paper presents the results of two case studies regarding the wind farm layout optimization problem. We asked a general audience to take part in the studies that we designed, and nine individuals participated. Case study 1 considered variations in optimization strategies for a given simple Gaussian wake model. Participants were provided with a wake model that outputs annual energy production (AEP) for an input set of wind turbine locations. Participants used an optimization method of their choosing to find an optimal wind farm layout. Case study 2 looked at trade-offs in performance resulting from variation in both physics model and optimization strategy. For case study 2, participants calculated AEP using a wake model of their choice while also implementing their chosen optimization method. Participants then used their wake model to calculate the AEP all other participant submitted turbine configurations produce for cross-comparison results. for optimized turbine locations were then cross-compared by recalculating the AEP using every other participant's wake model. Results for case study 1 show that the best optimal wind farm layouts in this study were achieved by participants who used gradient-based optimization methods. A front-runner emerged with the Sparse Nonlinear OPTimizer plus Wake Expansion Continuation (SNOPTplusWEC) optimization method, which consistently discovered a higher AEP for each scenario. Results for case study 2 show that for small wind farms with few turbines, turbine placement on the wind farm boundary is superior. Conclusions for case study 2 were drawn from participant cross-comparison of results.

17 WIND ENERGY↗

DOE OSTI2019

Best Practices for Wake Model and Optimization Algorithm Selection in Wind Farm Layout Optimization: Preprint

17 WIND ENERGY↗

DOE OSTI2012

Development of a fully-consistent reduced order model to study instabilities in boiling water reactors

A simple nonlinear Reduced Order Model to study global, regional and local instabilities in Boiling Water Reactors is described. The ROM consists of three submodels: neutron-kinetic, thermal-hydraulic and heat-transfer models. The neutron-kinetic model allows representing the time evolution of the three first neutron kinetic modes: the fundamental, the first and the second azimuthal modes. The thermal-hydraulic model describes four heated channels in order to correctly simulate out-of-phase behavior. The coupling between the different submodels is performed via both void and Doppler feedback mechanisms. After proper spatial homogenization, the governing equations are discretized in the time-domain. Several modifications, compared to other existing ROMs, have been implemented, and are reported in this paper. One novelty of the ROM is the inclusion of both azimuthal modes, which allows to study combined instabilities (in-phase and out-of-phase), as well as to investigate the corresponding interference effects between them. The second modification concerns the precise estimation of so-called reactivity coefficients or C{sub mn}{sup *V,D} - coefficients by using direct cross-section data from SIMULATE-3 combined with the CORE SIM core simulator in order to calculate Eigenmodes. Furthermore, a non-uniform two-step axial power profile is introduced to simulate the separate heat production in the single and two-phase regions, respectively. An iterative procedure was developed to calculate the solution to the coupled neutron-kinetic/thermal-hydraulic static problem prior to solving the time-dependent problem. Besides, the possibility of taking into account the effect of local instabilities is demonstrated in a simplified manner. The present ROM is applied to the investigation of an actual instability that occurred at the Swedish Forsmark-1 BWR in 1996/1997. The results generated by the ROM are compared with real power plant measurements performed during stability tests and show a good qualitative agreement. The present study provides some insight in a deeper understanding of the physical principles which drive both core-wide and local instabilities. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRSDate not supplied

Control Effector Unsaturation Modification to the Cascading Generalized Inverse Control Allocation Algorithm

Control allocation has sufficiently progressed such that it is used in front-line fighter aircraft such as the F-18Superhornet and the F-35 Joint Strike Fighter. Published literature shows the F-35 utilizes Nonlinear Dynamic Inversion in conjunction with an Effector Blender that incorporates the Cascading Generalized Inverse control allocation algorithm. While the Cascading Generalized Inverse algorithm is one of the premier generalized inverse methods, it does suffer from three deficiencies. In particular, it suffers from an inability to achieve some desired outcomes, it intermittently provides non-optimal solutions and generally fails to preserve moment direction near maximal achievable moments. An effector unsaturation method based on a Scalar Difference Quadratic was first introduced and implemented on the iterative Prediction Method control allocation algorithm which was shown to consistently achieve optimal (weighted) control allocation solutions throughout the entire Attainable Moment Set while preserving desired moment direction. In this paper, the shortcomings of the Cascading Generalized Inverse algorithm are addressed by augmenting the baseline algorithm with Scalar Difference Quadratic unsaturation identification and location at each iteration. Numerical case studies demonstrate that the Modified Cascading Generalized Inverse algorithm resolves the aforementioned deficiencies.

Michael J Acheson↗

DOE OSTI2016

A DEEP SEARCH FOR ADDITIONAL SATELLITES AROUND THE DWARF PLANET HAUMEA

Haumea is a dwarf planet with two known satellites, an unusually high spin rate, and a large collisional family, making it one of the most interesting objects in the outer solar system. A fully self-consistent formation scenario responsible for the satellite and family formation is still elusive, but some processes predict the initial formation of many small moons, similar to the small moons recently discovered around Pluto. Deep searches for regular satellites around Kuiper belt objects are difficult due to observational limitations, but Haumea is one of the few for which sufficient data exist. We analyze Hubble Space Telescope (HST) observations, focusing on a 10-consecutive-orbit sequence obtained in 2010 July, to search for new very small satellites. To maximize the search depth, we implement and validate a nonlinear shift-and-stack method. No additional satellites of Haumea are found, but by implanting and recovering artificial sources, we characterize our sensitivity. At distances between ∼10,000 and ∼350,000 km from Haumea, satellites with radii as small as ∼10 km are ruled out, assuming an albedo (p≃0.7) similar to Haumea. We also rule out satellites larger than ≳40 km in most of the Hill sphere using other HST data. This search method rules out objects similar in size to the small moons of Pluto. By developing clear criteria for determining the number of nonlinear rates to use, we find that far fewer shift rates are required (∼35) than might be expected. The nonlinear shift-and-stack method to discover satellites (and other moving transients) is tractable, particularly in the regime where nonlinear motion begins to manifest itself.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

A Fast Solver for Implicit Integration of the Vlasov--Poisson System in the Eulerian Framework

In this paper, we present a domain decomposition algorithm to accelerate the solution of Eulerian-type discretizations of the linear, steady-state Vlasov equation. The steady-state solver then forms a key component in the implementation of fully implicit or nearly fully implicit temporal integrators for the nonlinear Vlasov--Poisson system. The solver relies on a particular decomposition of phase space that enables the use of sweeping techniques commonly used in radiation transport applications. The original linear system for the phase space unknowns is then replaced by a smaller linear system involving only unknowns on the boundary between subdomains, which can then be solved efficiently with Krylov methods such as GMRES. Steady-state solves are combined to form an implicit Runge--Kutta time integrator, and the Vlasov equation is coupled self-consistently to the Poisson equation via a linearized procedure or a nonlinear fixed-point method for the electric field. Finally, numerical results for standard test problems demonstrate the efficiency of the domain decomposition approach when compared to the direct application of an iterative solver to the original linear system.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1998

An approach to the design and implementation of spacecraft attitude control systems

Over 39 years and a long list of missions, the guidance, navigation, and control (GN&C) groups at the Goddard Space Flight Center have gradually developed approaches to the design and implementation of successful spacecraft attitude control systems. With the recent creation of the Guidance, Navigation, and Control Center at Goddard, there is a desire to document some of these design practices to help to ensure their consistent application in the future. In this paper, we will discuss the beginnings of this effort, drawing primarily on the experience of one of the past attitude control system (ACS) groups at Goddard (what was formerly known as Code 712, the Guidance, Navigation, and Control Branch). We will discuss the analysis and design methods and criteria used, including guidelines for linear and nonlinear analysis, as well as the use of low- and high-fidelity simulation for system design and verification of performance. Descriptions of typical ACS sensor and actuator hardware will be shown, and typical sensor/actuator suites for a variety of mission types detailed. A description of the software and hardware test effort will be given, along with an attempt to make some qualitative estimates on how much effort is involved. The spacecraft and GN&C subsystem review cycles will be discussed, giving an outline of what design reviews are typically held and what information should be presented at each stage. Finally, we will point out some of the lessons learned at Goddard.

ODonnell, James R., Jr.↗

NASA NTRS1998

An Approach to the Design and Implementation of Spacecraft Attitude Control Systems

Over 39 years and a long list of missions, the guidance, navigation, and control (GN&C) groups at the Goddard Space Flight Center have gradually developed approaches to the design and implementation of successful spacecraft attitude control systems. With the recent creation of the Guidance, Navigation, and Control Center at Goddard, there is a desire to document some of these design practices to help to ensure their consistent application in the future. In this paper, we will discuss the beginnings of this effort, drawing primarily on the experience of one of the past attitude control system (ACS) groups at Goddard (what was formerly known as Code 712, the Guidance, Navigation, and Control Branch). We will discuss the analysis and design methods and criteria used, including guidelines for linear and nonlinear analysis, as well as the use of low- and high-fidelity simulation for system design and verification of performance. Descriptions of typical ACS sensor and actuator hardware will be shown, and typical sensor/actuator suites for a variety of mission types detailed. A description of the software and hardware test effort will be given, along with an attempt to make some qualitative estimates on how much effort is involved. The spacecraft and GN&C subsystem review cycles will be discussed, giving an outline of what design reviews are typically held and .what information should be presented at each stage. Finally, we will point out some of the lessons learned at Goddard.

ODonnell, James R., Jr.↗

DOE OSTI2017

Validity of the Born approximation for beyond Gaussian weak lensing observables

Accurate forward modeling of weak lensing (WL) observables from cosmological parameters is necessary for upcoming galaxy surveys. Because WL probes structures in the nonlinear regime, analytical forward modeling is very challenging, if not impossible. Numerical simulations of WL features rely on ray tracing through the outputs of N-body simulations, which requires knowledge of the gravitational potential and accurate solvers for light ray trajectories. A less accurate procedure, based on the Born approximation, only requires knowledge of the density field, and can be implemented more efficiently and at a lower computational cost. In this work, we use simulations to show that deviations of the Born-approximated convergence power spectrum, skewness and kurtosis from their fully ray-traced counterparts are consistent with the smallest nontrivial O(Φ 3 ) post-Born corrections (so-called geodesic and lens-lens terms). Our results imply a cancellation among the larger O(Φ 4 ) (and higher order) terms, consistent with previous analytic work. We also find that cosmological parameter bias induced by the Born-approximated power spectrum is negligible even for a LSST-like survey, once galaxy shape noise is considered. When considering higher order statistics such as the κ skewness and kurtosis, however, we find significant bias of up to 2.5σ. Using the LensTools software suite, we show that the Born approximation saves a factor of 4 in computing time with respect to the full ray tracing in reconstructing the convergence.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2006

Pattern-Recognition Algorithm for Locking Laser Frequency

A computer program serves as part of a feedback control system that locks the frequency of a laser to one of the spectral peaks of cesium atoms in an optical absorption cell. The system analyzes a saturation absorption spectrum to find a target peak and commands a laser-frequency-control circuit to minimize an error signal representing the difference between the laser frequency and the target peak. The program implements an algorithm consisting of the following steps: Acquire a saturation absorption signal while scanning the laser through the frequency range of interest. Condition the signal by use of convolution filtering. Detect peaks. Match the peaks in the signal to a pattern of known spectral peaks by use of a pattern-recognition algorithm. Add missing peaks. Tune the laser to the desired peak and thereafter lock onto this peak. Finding and locking onto the desired peak is a challenging problem, given that the saturation absorption signal includes noise and other spurious signal components; the problem is further complicated by nonlinearity and shifting of the voltage-to-frequency correspondence. The pattern-recognition algorithm, which is based on Hausdorff distance, is what enables the program to meet these challenges.

Karayan, Vahag↗

NASA NTRS2016

Model Development for MODIS Thermal Band Electronic Crosstalk

MODerate-resolution Imaging Spectroradiometer (MODIS) has 36 bands. Among them, 16 thermal emissive bands covering a wavelength range from 3.8 to 14.4 m. After 16 years on-orbit operation, the electronic crosstalk of a few Terra MODIS thermal emissive bands developed substantial issues that cause biases in the EV brightness temperature measurements and surface feature contamination. The crosstalk effects on band 27 with center wavelength at 6.7 m and band 29 at 8.5 m increased significantly in recent years, affecting downstream products such as water vapor and cloud mask. The crosstalk effect is evident in the near-monthly scheduled lunar measurements, from which the crosstalk coefficients can be derived. The development of an alternative approach is very helpful for independent verification.In this work, a physical model was developed to assess the crosstalk impact on calibration as well as in Earth view brightness temperature retrieval. This model was applied to Terra MODIS band 29 empirically to correct the Earth brightness temperature measurements. In the model development, the detectors nonlinear response is considered. The impact of the electronic crosstalk is assessed in two steps. The first step consists of determining the impact on calibration using the on-board blackbody (BB). Due to the detectors nonlinear response and large background signal, both linear and nonlinear coefficients are affected by the crosstalk from sending bands. The second step is to calculate the effects on the Earth view brightness temperature retrieval. The effects include those from affected calibration coefficients and the contamination of Earth view measurements. This model links the measurement bias with crosstalk coefficients, detector non-linearity, and the ratio of Earth measurements between the sending and receiving bands. The correction of the electronic cross talk can be implemented empirically from the processed bias at different brightness temperature. The implementation can be done through two approaches. As routine calibration assessment for thermal infrared bands, the trending over select Earth scenes is processed for all the detectors in a band and the band averaged bias is derived at a certain time. In this case, the correction of an affected band can be made using the regression of the model with band averaged bias and then corrections of detector differences are applied. The second approach requires the trending for individual detectors and the bias for each detector is used for regression with the model. A test using the first approach was made for Terra MODIS band 29 with the biases derived from long-term trending of brightness temperature over ocean and Dome-C.

Chang, Tiejun↗

DOE OSTI2017

A scalable nonlinear fluid–structure interaction solver based on a Schwarz preconditioner with isogeometric unstructured coarse spaces in 3D

Nonlinear fluid-structure interaction (FSI) problems on unstructured meshes in 3D appear many applications in science and engineering, such as vibration analysis of aircrafts and patient-specific diagnosis of cardiovascular diseases. In this work, we develop a highly scalable, parallel algorithmic and software framework for FSI problems consisting of a nonlinear fluid system and a nonlinear solid system, that are coupled monolithically. The FSI system is discretized by a stabilized finite element method in space and a fully implicit backward difference scheme in time. To solve the large, sparse system of nonlinear algebraic equations at each time step, we propose an inexact Newton-Krylov method together with a multilevel, smoothed Schwarz preconditioner with isogeometric coarse meshes generated by a geometry preserving coarsening algorithm. Here ''geometry'' includes the boundary of the computational domain and the wet interface between the fluid and the solid. We show numerically that the proposed algorithm and implementation are highly scalable in terms of the number of linear and nonlinear iterations and the total compute time on a supercomputer with more than 10,000 processor cores for several problems with hundreds of millions of unknowns.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Constraining nuclear data via cosmological observations: Neutrino energy transport and big bang nucleosynthesis

Here, we introduce a new computational capability that moves toward a self-consistent calculation of neutrino transport and nuclear reactions for big bang nucleosynthesis (BBN). Such a self-consistent approach is needed to be able to extract detailed information about nuclear reactions and physics beyond the standard model from precision cosmological observations of primordial nuclides and the cosmic microwave background radiation. We also calculate the evolution of the early universe through the epochs of weak decoupling, weak freeze-out and big bang nucleosynthesis (BBN) by simultaneously coupling a full strong, electromagnetic, and weak nuclear reaction network with a multi-energy group Boltzmann neutrino energy transport scheme. The modular structure of our approach allows the dissection of the relative contributions of each process responsible for evolving the dynamics of the early universe. Such an approach allows a detailed account of the evolution of the active neutrino energy distribution functions alongside and self-consistently with the nuclear reactions and entropy/heat generation and flow between the neutrino and photon/electron/positron/baryon plasma components. Our calculations reveal nonlinear feedback in the time evolution of neutrino distribution functions and plasma thermodynamic conditions. We discuss the time development of neutrino spectral distortions and concomitant entropy production and extraction from the plasma. These effects result in changes in the computed values of the BBN deuterium and helium-4 yields that are on the order of a half-percent relative to a baseline standard BBN calculation with no neutrino transport. This is an order of magnitude larger effect than in previous estimates. For particular implementations of quantum corrections in plasma thermodynamics, our calculations show a 0.4% increase in deuterium and a 0.6% decrease in 4 He over our baseline. The magnitude of these changes are on the order of uncertainties in the nuclear physics for the case of deuterium and are potentially significant for the error budget of helium in upcoming cosmological observations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Proper orthogonal decomposition based reduced-order modeling of flux-Limited gray thermal radiation

Here, in this work, a proper orthogonal decomposition (POD) based reduced-order model (ROM) is developed to solve gray, flux-limited thermal radiation diffusion. We focus on the variable opacity radiation penetration benchmark posed by Olson, Auer, and Hall. The T -3 relationship for opacity in conjunction with high-temperature radiation penetrating an initially cold material produces a strong thermal radiation shock. This class of problems is particularly challenging for standard POD-based reduced-order modeling due to the nonlinearities presented by 1) the T 4 source term, and 2) flux-limited diffusion operator. To address these challenges and develop a cost competitive ROM, we employ a “hyper-reduction” technique through discrete empirical interpolation (DEIM) and allow for adaptive reduced-order projections through principal interval decomposition (PID). Performance of the proposed methodology is quantified by comparing the cost savings and accuracy relative to a full-order computation. Reference solutions and snapshot data are obtained through a full-order calculation performed by the University of Chicago maintained astrophysics code, FLASH. For consistency and potential extensibility, the developed ROM is also implemented in FLASH. We find that in the initialization regime, where the thermal radiation wave is initially created by the warming of the material, this class of problems is highly reducible and suitable for POD-based ROMs. However, the strong convective nature of the wave propagation regime is less reducible and more challenging to create an efficient ROM.

42 ENGINEERING↗

NASA NTRS1993

Precision attitude control for tethered satellites

Tethered spacecraft are particularly well suited to serve as isolation platforms for space-borne observatories. It has previously been shown that, due to the relatively large tether force, conventional means of performing attitude control for tethered satellites are inefficient for any mission with pointing requirements more stringent than about 1 deg. A particularly effective method of implementing attitude control for tethered satellites is to use the tether tension force to generate control moments by moving the tether attach point relative to the subsatellite center of mass. This paper presents the development of a precision pointing control algorithm for tethered satellites and the simulation of the control system with laboratory hardware. The control algorithm consists of a linear quadratic regulator feedback law and a Kalman filter. The control algorithm has been shown to regulate the vehicle orientation to within 0.60 arcsec rms. This level of precision was achieved only after including a mass center estimator and accurately modeling the effects of the nonlinear attach point motion actuator.

Kline-Schoder, Robert J.↗

NASA NTRS1991

Control algorithm implementation for a redundant degree of freedom manipulator

This project's purpose is to develop and implement control algorithms for a kinematically redundant robotic manipulator. The manipulator is being developed concurrently by Odetics Inc., under internal research and development funding. This SBIR contract supports algorithm conception, development, and simulation, as well as software implementation and integration with the manipulator hardware. The Odetics Dexterous Manipulator is a lightweight, high strength, modular manipulator being developed for space and commercial applications. It has seven fully active degrees of freedom, is electrically powered, and is fully operational in 1 G. The manipulator consists of five self-contained modules. These modules join via simple quick-disconnect couplings and self-mating connectors which allow rapid assembly/disassembly for reconfiguration, transport, or servicing. Each joint incorporates a unique drive train design which provides zero backlash operation, is insensitive to wear, and is single fault tolerant to motor or servo amplifier failure. The sensing system is also designed to be single fault tolerant. Although the initial prototype is not space qualified, the design is well-suited to meeting space qualification requirements. The control algorithm design approach is to develop a hierarchical system with well defined access and interfaces at each level. The high level endpoint/configuration control algorithm transforms manipulator endpoint position/orientation commands to joint angle commands, providing task space motion. At the same time, the kinematic redundancy is resolved by controlling the configuration (pose) of the manipulator, using several different optimizing criteria. The center level of the hierarchy servos the joints to their commanded trajectories using both linear feedback and model-based nonlinear control techniques. The lowest control level uses sensed joint torque to close torque servo loops, with the goal of improving the manipulator dynamic behavior. The control algorithms are subjected to a dynamic simulation before implementation.

Cohan, Steve↗

DOE OSTI2011

Approaches to cancer assessment in EPA's Integrated Risk Information System

The U.S. Environmental Protection Agency's (EPA) Integrated Risk Information System (IRIS) Program develops assessments of health effects that may result from chronic exposure to chemicals in the environment. The IRIS database contains more than 540 assessments. When supported by available data, IRIS assessments provide quantitative analyses of carcinogenic effects. Since publication of EPA's 2005 Guidelines for Carcinogen Risk Assessment, IRIS cancer assessments have implemented new approaches recommended in these guidelines and expanded the use of complex scientific methods to perform quantitative dose-response assessments. Two case studies of the application of the mode of action framework from the 2005 Cancer Guidelines are presented in this paper. The first is a case study of 1,2,3-trichloropropane, as an example of a chemical with a mutagenic mode of carcinogenic action thus warranting the application of age-dependent adjustment factors for early-life exposure; the second is a case study of ethylene glycol monobutyl ether, as an example of a chemical with a carcinogenic action consistent with a nonlinear extrapolation approach. The use of physiologically based pharmacokinetic (PBPK) modeling to quantify interindividual variability and account for human parameter uncertainty as part of a quantitative cancer assessment is illustrated using a case study involving probabilistic PBPK modeling for dichloromethane. We also discuss statistical issues in assessing trends and model fit for tumor dose-response data, analysis of the combined risk from multiple types of tumors, and application of life-table methods for using human data to derive cancer risk estimates. These issues reflect the complexity and challenges faced in assessing the carcinogenic risks from exposure to environmental chemicals, and provide a view of the current trends in IRIS carcinogenicity risk assessment.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

Neutrino energy transport in weak decoupling and big bang nucleosynthesis

In this study, we calculate the evolution of the early universe through the epochs of weak decoupling, weak freeze-out and big bang nucleosynthesis (BBN) by simultaneously coupling a full strong, electromagnetic, and weak nuclear reaction network with a multienergy group Boltzmann neutrino energy transport scheme. The modular structure of our code provides the ability to dissect the relative contributions of each process responsible for evolving the dynamics of the early universe in the absence of neutrino flavor oscillations. Such an approach allows a detailed accounting of the evolution of the νe, ν¯e, νμ, ν¯μ, ντ, ν¯τ energy distribution functions alongside and self-consistently with the nuclear reactions and entropy/heat generation and flow between the neutrino and photon/electron/positron/baryon plasma components. This calculation reveals nonlinear feedback in the time evolution of neutrino distribution functions and plasma thermodynamic conditions (e.g., electron-positron pair densities), with implications for the phasing between scale factor and plasma temperature; the neutron-to-proton ratio; light-element abundance histories; and the cosmological parameter N eff . We find that our approach of following the time development of neutrino spectral distortions and concomitant entropy production and extraction from the plasma results in changes in the computed value of the BBN deuterium yield. For example, for particular implementations of quantum corrections in plasma thermodynamics, our calculations show a 0.4% increase in deuterium. These changes are potentially significant in the context of anticipated improvements in observational and nuclear physics uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗