Search NASA⌕ Search

SEARCH · Search NASA

Results for “fully linear model”

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 109 records · Page 6

Nonlinear Dynamic Models in Advanced Life Support

To facilitate analysis, ALS systems are often assumed to be linear and time invariant, but they usually have important nonlinear and dynamic aspects. Nonlinear dynamic behavior can be caused by time varying inputs, changes in system parameters, nonlinear system functions, closed loop feedback delays, and limits on buffer storage or processing rates. Dynamic models are usually cataloged according to the number of state variables. The simplest dynamic models are linear, using only integration, multiplication, addition, and subtraction of the state variables. A general linear model with only two state variables can produce all the possible dynamic behavior of linear systems with many state variables, including stability, oscillation, or exponential growth and decay. Linear systems can be described using mathematical analysis. Nonlinear dynamics can be fully explored only by computer simulations of models. Unexpected behavior is produced by simple models having only two or three state variables with simple mathematical relations between them. Closed loop feedback delays are a major source of system instability. Exceeding limits on buffer storage or processing rates forces systems to change operating mode. Different equilibrium points may be reached from different initial conditions. Instead of one stable equilibrium point, the system may have several equilibrium points, oscillate at different frequencies, or even behave chaotically, depending on the system inputs and initial conditions. The frequency spectrum of an output oscillation may contain harmonics and the sums and differences of input frequencies, but it may also contain a stable limit cycle oscillation not related to input frequencies. We must investigate the nonlinear dynamic aspects of advanced life support systems to understand and counter undesirable behavior.

Jones, Harry↗

Theory and measurement of plasmonic terahertz detector response to large signals

Electron gas in the conduction channel of a Field Effect Transistor (FET) can support collective plasma oscillations tunable by the gate voltage. In the Dyakonov-Shur terahertz (THz) detector, nonlinearities in the plasma wave propagation in the gated channel of a FET lead to a constant source-to-drain voltage providing the detector output. We present the detector theory in the frame of the hydrodynamic model using the electron plasma Navier-Stokes and thermal transport equations, thus fully accounting for the hydrodynamic non-linearity, the viscosity, and pressure gradients in the detector response. Both resonant and broadband operations of the high electron mobility transistor (HEMT) based plasmonic detectors are described by this model. The relation between the electron channel density and gate voltage was modeled by the unified charge control model applicable both above and below the threshold voltage. The theoretical results are compared with the response measured in the short channel InGaAs HEMT and the analytical approximation. The THz source was operating at 1.63 THz, and the response was measured at varying signal intensities. The response of the detector operated in the open drain mode was measured above and below the threshold, and the theoretical and experimental results are shown to be in good agreement.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A Fast Scalable Quasi-Static Time Series Analysis Method for PV Impact Studies Using Linear Sensitivity Model

Understanding the impact of distributed photovoltaic (PV) resources on various elements of the distribution feeder is imperative for their cost effective integration. A year-long quasi-static time series (QSTS) simulation at 1-second granularity is often necessary to fully study these impacts. However, the significant computational burden associated with running QSTS simulations is a major challenge to their adoption. In this paper, we propose a fast scalable QSTS simulation algorithm that is based on a linear sensitivity model for estimating voltage-related PV impact metrics of a three-phase unbalanced, nonradial distribution system with various discrete step control elements including tap changing transformers and capacitor banks. Here, the algorithm relies on computing voltage sensitivities while taking into account all the effects of discrete controllable elements in the circuit. Consequently, the proposed sensitivity model can accurately estimate the state of controllers at each time step and the number of control actions throughout the year. For the test case of a real distribution feeder with 2969 buses (5469 nodes), 6 load/PV time series power profiles, and 9 voltage regulating elements including controller delays, the proposed algorithm demonstrates a dramatic time reduction, more than 180 times faster than traditional QSTS techniques.

14 SOLAR ENERGY↗

Numerical Study of the Detonation Wave Structure in a Linear Model Detonation Engine

Rotating detonation engines (RDE) are increasingly favored as a viable approach to gener- ating pressure gain combustion for both propulsion and power generation purposes. Practical designs involve discrete injection of fuel and air, which then partially mix to produce the reac- tive mixture that is processed by a continuously moving detonation wave within the detonation chamber. In this study, the wave structure of the detonation wave within a linearized RDE with fully premixed stoichiometric fuel and air injection is analyzed. The direct numerical simula- tion of the linearized model detonation engine (LMDE) is compared to experimental results, and wave behavior is correlated to the wave structure and flow turbulence. The flow properties and chemical composition of the gases across the detonation wave is studied to examine the reaction zone characteristics. This analysis is extended to equivalence ratios of φ = 0.8 and 1.2.

Prakash, Supraj↗

Statistical simulation of the magnetorotational dynamo

We analyze turbulence and dynamo induced by the magnetorotational instability (MRI) using quasi-linear statistical simulation methods. We find that homogenous turbulence is unstable to a large scale dynamo instability, which saturates to an inhomogenous equilibrium with a very strong dependence on the magnetic Prandtl number (Pm). Despite its enormously reduced nonlinearity, the quasi-linear model exhibits the same qualitative scaling of angular momentum transport with Pm as fully nonlinear turbulence. This demonstrates the relationship of recent convergence problems to the large scale dynamo and suggests possible methods for studying astrophysically relevant regimes at very low or high Pm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A fuzzy set preference model for market share analysis

Consumer preference models are widely used in new product design, marketing management, pricing, and market segmentation. The success of new products depends on accurate market share prediction and design decisions based on consumer preferences. The vague linguistic nature of consumer preferences and product attributes, combined with the substantial differences between individuals, creates a formidable challenge to marketing models. The most widely used methodology is conjoint analysis. Conjoint models, as currently implemented, represent linguistic preferences as ratio or interval-scaled numbers, use only numeric product attributes, and require aggregation of individuals for estimation purposes. It is not surprising that these models are costly to implement, are inflexible, and have a predictive validity that is not substantially better than chance. This affects the accuracy of market share estimates. A fuzzy set preference model can easily represent linguistic variables either in consumer preferences or product attributes with minimal measurement requirements (ordinal scales), while still estimating overall preferences suitable for market share prediction. This approach results in flexible individual-level conjoint models which can provide more accurate market share estimates from a smaller number of more meaningful consumer ratings. Fuzzy sets can be incorporated within existing preference model structures, such as a linear combination, using the techniques developed for conjoint analysis and market share estimation. The purpose of this article is to develop and fully test a fuzzy set preference model which can represent linguistic variables in individual-level models implemented in parallel with existing conjoint models. The potential improvements in market share prediction and predictive validity can substantially improve management decisions about what to make (product design), for whom to make it (market segmentation), and how much to make (market share prediction).

Turksen, I. B.↗

Stability analysis of WEST L-mode discharges with improved confinement from boron powder injection

WEST L-mode plasmas with dominant electron heating and no core torque source have observed improvements in confinement during boron (B) powder injection. These results are reminiscent of previous powder injection experiments on other devices and gaseous impurity seeding experiments on WEST. During powder injection, the stored energy increased up to 25% due to enhanced ion and electron heat and particle confinement. The improvements in confinement were not indicative of an L-H transition. To identify the dominant mechanisms and the causality chain behind these improvements in confinement, we employ interpretative modeling using METIS, predictive integrated modeling using a high-fidelity plasma simulator (HFPS), and stand-alone gyrokinetic simulations using quasi-linear gyrokinetic code. Interpretative modeling with METIS allowed for the estimation of missing data while maintaining good overall consistency with experiment. These results provided the initial conditions for fully predictive flux driven simulations using the HFPS. From these simulations, quasi-linear gyrokinetic analysis was performed at ρ = 0.5 and ρ = 0.65. Collisionality was found to be a strong candidate for the turbulence suppression mechanism at ρ = 0.5, while a combination of collisionality and the T e /T i ratio was found to be the likely mechanism at ρ = 0.65. The results additionally suggested that increased Z eff (through main ion dilution) could play a role in the improved confinement, but this could not be confirmed due to a lack of experimental measurements. The modeling framework established here can now be used to evaluate and exploit a variety of future powder injection experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimal Location of Radiation Therapy Centers With Respect to Geographic Access

Purpose: To develop a framework with which to evaluate locations of radiation therapy (RT) centers in a region based on geographic access. Methods and Materials: Patient records were obtained for all external beam radiation therapy started in 2011 for the province of British Columbia, Canada. Two metrics of geographic access were defined. The primary analysis was percentage of patients (coverage) within a 90-minute drive from an RT center (C90), and the secondary analysis was the average drive time (ADT) to an RT center. An integer programming model was developed to determine optimal center locations, catchment areas, and capacity required under different scenarios. Results: Records consisted of 11,096 courses of radiation corresponding to 161,616 fractions. Baseline geographic access was estimated at 102.5 minutes ADT (each way, per fraction) and 75.9% C90. Adding 2 and 3 new centers increased C90 to 88% and 92%, respectively, and decreased ADT by between 43% and 61%, respectively. A scenario in which RT was provided in every potential location that could support at least 1 fully utilized linear accelerator resulted in 35.3 minutes' ADT and 93.6% C90. Conclusions: The proposed framework and model provide a data-driven means to quantitatively evaluate alternative configurations of a regional RT system. Results suggest that the choice of location for future centers can significantly improve geographic access to RT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Newton-Krylov method in solving drift-flux two-phase flow problems for both pre-CHF and post-CHF conditions at steady state

Furthermore, this paper presents numerical simulations of two-phase flow heat transfer for both pre-CHF and post-CHF conditions typical in light water nuclear reactors. The conjugate heat transfer problem is modeled with one-dimensional drift-flux two-phase flow model for coolant flow and two-dimensional heat conduction equation in solid walls. Fully implicit time integration schemes were employed, and the resulted non-linear equations were solved with a Newton-Krylov method. For the drift-flux two-phase flow model, the EPRI drift-flux closure correlation is used to model the relative motion between the two phases. Additional closure correlations were adopted from REALP5-3D to determine two-phase flow regime, wall boiling heat and mass transfer, interfacial heat and mass transfer, and two-phase flow frictional pressure drop. Three sets of experiments were used for code validation, including Bartolomei subcooled flow boiling in vertical pipes, FRIGG boiling tests in vertical bundles, and Bennett heated tube tests. RELAP5-3D simulation results were also provided to perform code-to-code benchmark. Overall, numerical results of this work showed good agreements with both experimental data and RELAP5-3D simulation results. Discrepancy between them were also identified in post-CHF region, which suggests necessary further improvement.

42 ENGINEERING↗

Dynamical attractors in contracting spacetimes dominated by kinetically coupled scalar fields

We present non-perturbative numerical relativity simulations of slowly contracting spacetimes in which the scalar field driving slow contraction is coupled to a second scalar field through an exponential non-linear σ model-type kinetic interaction. These models are important because they can generate a nearly scale-invariant spectrum of super-Hubble density fluctuations fully consistent with cosmic microwave background observations. We show that the non-linear evolution rapidly approaches a homogeneous, isotropic and flat Friedmann-Robertson-Walker (FRW) geometry for a wide range of inhomogeneous and anisotropic initial conditions. Ultimately, we find, the kinetic coupling causes the evolution to deflect away from flat FRW and towards a novel Kasner-like stationary point, but in general this occurs on time scales that are too long to be observationally relevant.

79 ASTRONOMY AND ASTROPHYSICS↗

Analytical and numerical studies of the thermocapillary flow in a uniformly floating zone

The microgravity environment of an orbiting vehicle permits crystal growth experiments in the presence of greatly reduced buoyant convection in the liquid melt. Crystals grown in ground-based laboratories do not achieve their potential properties because of dopant variations caused by flow in the melt. The floating zone crystal growing system is widely used to produce crystals of silicon and other materials. However, in this system the temperature gradient on the free sidewall surface of the melt is the source of a thermocapillary flow which does not disappear in the low-gravity environment. The idea of using a uniform rotation of the floating zone system to confine the thermocapillary flow to the melt sidewall leaving the interior of the melt passive is examined. A cylinder of fluid with an axial temperature gradient imposed on the cylindrical sidewall is considered. A half zone and the linearized, axisymmetric flow in the absence of crystal growth is examined. Rotation is found to confine the linear thermocapillary flow. A simplified model is extended to a full zone and both linear and nonlinear thermocapillary flows are studied theoretically. Analytical and numerical methods are used for the linear flows and numerical methods for the nonlinear flows. It was found that the linear flows in the full zone have more complicated and thicker boundary layer structures than in the half zone, and that these flows are also confined by the rotation. However, for the simplified model considered and for realistic values for silicon, the thermocapillary flow is not linear. The fully nonlinear flow is strong and unsteady (a weak oscillation is present) and it penetrates the interior. Some non-rotating flow results are also presented. Since silicon as a large value of thermal conductivity, one would expect the temperature fields to be determined by conduction alone. This is true for the linear and weakly nonlinear flows, but for the stronger nonlinear flow the results show that temperature advection is also important. Uniform rotation may still be a means of confining the flow and the results obtained define the procedure to be used to examine this hypothesis.

Fowlis, W. W.↗

Galileons coupled to massive gravity: general analysis and cosmological solutions

We further develop the framework for coupling galileons and Dirac-Born-Infeld (DBI) scalar fields to a massive graviton while retaining both the non-linear symmetries of the scalars and ghost-freedom of the theory. The general construction is recast in terms of vielbeins which simplifies calculations and allows for compact expressions. Expressions for the general form of the action are derived, with special emphasis on those models which descend from maximally symmetric spaces. We demonstrate the existence of maximally symmetric solutions to the fully non-linear theory and analyze their spectrum of quadratic fluctuations. Finally, we consider self-accelerating cosmological solutions and study their perturbations, showing that the vector and scalar modes have vanishing kinetic terms.

79 ASTRONOMY AND ASTROPHYSICS↗

Prediction of Transonic Vortex Flows Using Linear and Nonlinear Turbulent Eddy Viscosity Models

Three-dimensional transonic flow over a delta wing is investigated with a focus on the effect of transition and influence of turbulence stress anisotropies. The performance of linear eddy viscosity models and an explicit algebraic stress model is assessed at the start of vortex flow, and the results compared with experimental data. To assess the effect of transition location, computations that either fix transition or are fully turbulent are performed. To assess the effect of the turbulent stress anisotropy, comparisons are made between predictions from the algebraic stress model and the linear eddy viscosity models. Both transition location and turbulent stress anisotropy significantly affect the 3D flow field. The most significant effect is found to be the modeling of transition location. At a Mach number of 0.90, the computed solution changes character from steady to unsteady depending on transition onset. Accounting for the anisotropies in the turbulent stresses also considerably impacts the flow, most notably in the outboard region of flow separation.

Bartels, Robert E.↗

Active Flutter Suppression Controllers Derived from Linear and Nonlinear Aerodynamics: Application to a Transport Aircraft Model

Active flutter suppression has been demonstrated in simulation by many researchers, generally using methods based on linear aerodynamics and often with simplistic geometries. In this paper, active flutter suppression is demonstrated in a simulation using a Navier-Stokes aerodynamics code, FUN3D (Fully Unstructured Navier-Stokes Three-Dimensional), and a realistic transport aircraft configuration. This is accomplished using simple observer-feedback controllers derived from linear aeroelastic models, including reduced order models built via FUN3D data. The development of these reduced order models is described here. It is shown that controllers derived from reduced order models of the nonlinear aerodynamics outperform controllers based on linear aerodynamics.

Waite, Josiah M.↗

Fully compressible solutions for early stage Richtmyer–Meshkov instability

Here, we will consider the effects of compressibility and viscosity on the early dynamics of the Richtmyer–Meshkov instability (RMI). In particular, we will combine theory, scaling, and high resolution simulation of RMI to probe the details of the initial compression and the subsequent viscous damping as a shock interacts with a density discontinuity. We will propose a refinement of the classic 1D model for the linear regime of RMI that, for small initial perturbation wavelengths, more accurately reproduces the 2D dynamics of a fully resolved numerical simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Artificial dissipation models for the Euler equations

Various artificial dissipation models which are used with central difference algorithms for the Euler equations are analyzed for their effect on accuracy, stability and convergence rates. In particular, linear and nonlinear models are investigated using an implicit approximate factorization code (ARC2D) for transonic airfoils. Fully implicit application of the dissipation models is shown to improve robustness and convergence rates. The treatment of dissipation models at boundaries will be examined. It will be shown that accurate, error free solutions with sharp shocks can be obtained using a central difference algorithm coupled with an appropriate nonlinear artificial dissipation model.

Pulliam, T. H.↗

Preconditioned conjugate gradient methods for the compressible Navier-Stokes equations

The compressible Navier-Stokes equations are solved for a variety of two-dimensional inviscid and viscous problems by preconditioned conjugate gradient-like algorithms. Roe's flux difference splitting technique is used to discretize the inviscid fluxes. The viscous terms are discretized by using central differences. An algebraic turbulence model is also incorporated. The system of linear equations which arises out of the linearization of a fully implicit scheme is solved iteratively by the well known methods of GMRES (Generalized Minimum Residual technique) and Chebyschev iteration. Incomplete LU factorization and block diagonal factorization are used as preconditioners. The resulting algorithm is competitive with the best current schemes, but has wide applications in parallel computing and unstructured mesh computations.

Venkatakrishnan, V.↗

(abstract) Observation by Clementine of a Hugh Opposition Surge on the Moon at Very Small Solar Phase Angles

The Clementine mission enabled the first quantative observations of the Moon at very small solar phase angles. It is well established from Earth-based observations that the Moon exhibits a non-linear increase in brightness as its face becomes fully illuminated to a terrestrial observer. Because the models of the opposition effect are sensitive indicators of surfical compaction state and particle size, observations at small solar phase angles are important to obtain. A recent model for optical coherent backscatter seeks to explain a narrow opposition spike at very small phase angles (< 1 degree) seen on several icy satellites , including Europa, Icarus 90, and Oberon . Over 90 images of the Moon's surface at phase angles less than 0.5(deg) (the minimum phase angle observable from Earth due to the Moon's angular size) were obtained by the Clementine spacecraft. Our analysis of these images shows the moon exhibits a surge in brightness of approximately 20% below 0.25(deg), comparable to the values observed on icy satellites. No color dependence was detected in the lunar phase curve below 2 degrees.

partical size icy satellites Europa Icarus 90 Ober↗