Search NASA⌕ Search

SEARCH · Search NASA

Results for “explicit solution”

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 559 records · Page 31

Self-similar Reynolds-averaged mechanical–scalar turbulence models for reshocked Richtmyer–Meshkov instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Richtmyer–Meshkov instability-induced mixing in planar geometry derived in the small Atwood number limit are extended to construct models for reshocked Richtmyer–Meshkov mixing. In this study, the models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ modeled either differentially or algebraically. The three- and four-equation models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing. Mixing layer growth parameters and other physical observables were obtained explicitly as functions of the model coefficients and were used to calibrate the model coefficients. Here, the solutions for the singly shocked Richtmyer–Meshkov case for the mixing layer width and the turbulent fields are used to construct piecewise-continuous generalizations of these quantities for times after reshock. For generality, the post-reshock mixing layer width is not assumed to grow with the same power-law as the pre-reshock width, and an impulsive approximation applied to Rayleigh–Taylor instability growth is used to establish the expression for the post-reshock width. A four-equation model is then used to illustrate the spatiotemporal behavior of the mean and turbulent fields and late-time turbulent equation budgets across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for reshocked Richtmyer–Meshkov instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗

A time dependent difference theory for sound propagation in ducts with flow

A time dependent numerical solution of the linearized continuity and momentum equation is developed for sound propagation in a two-dimensional straight hard or soft wall duct with a sheared mean flow. The time dependent governing acoustic-difference equations and boundary conditions are developed along with a numerical determination of the maximum stable time increments. The analysis begins with a harmonic noise source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the transient as well as the 'steady' state solution of the acoustic field. Example calculations are presented for sound propagation in hard and soft wall ducts, with no flow and with plug flow. Although the problem with sheared flow has been formulated and programmed, sample calculations have not yet been examined. So far, the time dependent finite difference analysis has been found to be superior to the steady state finite difference and finite element techniques because of shorter solution times and the elimination of large matrix storage requirements.

Baumeister, K. J.↗

Interactions of Wind Energy Project Siting, Wind Resource Potential, and the Evolution of the U.S. Power System

The interactions of wind energy project siting, resource potential, and power system evolution are critical to understand given interest in high renewable energy systems, particularly as wind costs have fallen and grid operators implement solutions to manage variable generation. This study applies an integrative approach that combines spatially-explicit resource potential analysis with high spatial resolution U.S. electricity system modeling. Multiple wind supply curves, representing variations in siting regimes that account for interactions with built infrastructure, regulatory, physical, and social land use factors, are evaluated to determine how siting could influence prospective wind development. Different siting regimes lead to variations in future wind capacity, with the greatest impacts observed under scenarios with high demand for clean energy. With the tightest emissions limit modeled, 2050 onshore wind capacity varied by +7% (49 GW) in the least-restrictive siting regime to a decrease of 37% (270 GW) under the most-constrained case. More-stringent siting restrictions lead to higher electricity prices and emissions that should be weighed with local impacts. Under restrictive siting regimes, wind deployment is also sensitive to transmission availability and wind plant design. Overall, the findings highlight the importance of local land use considerations in regional and national power system planning.

17 WIND ENERGY↗

Fast estimation of ion-pairing for screening electrolytes: A cluster can approximate a bulk liquid

The propensity for ion-pairing can often dictate the thermodynamic and kinetic properties of electrolyte solutions. Fast and accurate estimates of ion-pairing can thus be extremely valuable for supplementing design and screening efforts for novel electrolytes. We introduce an efficient cluster model to estimate the local ion-pair potential-of-mean-force between ionic solutes in electrolytes. The model incorporates an ion-pair and a few layers of explicit solvent in a gas-phase cluster and leverages an enhanced sampling approach to achieve high efficiency and accuracy. We employ harmonic restraints to prevent solvent escape from the cluster and restrict sampling of large inter-ion distances. We develop a cluster ion-pair sampling tool that implements our cluster model and demonstrate its potential utility for screening simple and poly-electrolyte systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multirate Exponential Rosenbrock Methods

In this paper we propose a novel class of methods for high-order accurate integration of multirate systems of ordinary differential equation initial-value problems. The proposed methods construct multirate schemes by approximating the action of matrix φ functions within explicit exponential Rosenbrock (ExpRB) methods, thereby called multirate ExpRB (MERB) methods. They consist of the solution to a sequence of modified “fast” initial-value problems, which may themselves be approximated through subcycling any desired initial-value problem solver. In addition to proving how to construct MERB methods from certain classes of ExpRB methods, we provide rigorous convergence analysis of these methods and derive efficient MERB schemes of orders 2 through 6 (the highest-order infinitesimal multirate methods to date). Lastly, we then present numerical simulations to confirm these theoretical convergence rates and to compare the efficiency of MERB methods against other recently introduced high-order multirate methods.

97 MATHEMATICS AND COMPUTING↗

Designing an Optimal Sensor Network via Minimizing Information Loss

Optimal experimental design is a classic topic in statistics, with many well-studied problems, applications, and solutions. The design problem we study is the placement of sensors to monitor spatiotemporal processes, explicitly accounting for the temporal dimension in our modeling and optimization. We observe that recent advancements in computational sciences often yield large datasets based on physics-based simulations, which are rarely leveraged in experimental design. We introduce a novel model-based sensor placement criterion, along with a highly-efficient optimization algorithm, which integrates physics-based simulations and Bayesian experimental design principles to identify sensor networks that “minimize information loss” from simulated data. Our technique relies on sparse variational inference and (separable) Gauss-Markov priors, and thus may adapt many techniques from Bayesian experimental design. We validate our method through a case study monitoring air temperature in Phoenix, Arizona, using state-of-the-art physics-based simulations. Our results show our framework to be superior to random or quasi-random sampling, particularly with a limited number of sensors. We conclude by discussing practical considerations and implications of our framework, including more complex modeling tools and real-world deployments.

54 ENVIRONMENTAL SCIENCES↗

Upwind second-order difference schemes and applications in unsteady aerodynamic flows

Explicit second-order upwind difference schemes in combination with spatially symmetric schemes can produce larger stability bounds and better numerical resolution than symmetric schemes alone. However, if conservation form is essential, a special operator is required for transition between schemes. An operational approach has been devised for deriving transition operators so that strict conservation and local consistency are maintained. Various aspects of hybrid schemes are studied numerically for model linear and nonlinear equations. To demonstrate the utility of combining two different algorithms, MacCormack's explicit, noncentered, second-order method is combined with a completely upwind version, and numerical solutions of the Euler equations are obtained for two-dimensional, transonic flows with embedded supersonic regions and shock waves. The general utility of the operational approach for combining schemes is emphasized by deriving a second-order conservative scheme for the steady transonic small-disturbance potential equation.

Warming, R. F.↗

High frequency sound emission from moving point multipole sources embedded in arbitrary transversely sheared mean flows

Formulas are derived for the high frequency sound emission from moving point multipole sources embedded in an arbitrary unidirectional transversely sheared mean flow. The results are used to study the sound generated by non-axisymmetric turbulent jets. The effect of the asymmetry in both the mean flow and the source distribution is accounted for by a 'circumferential directivity factor', which is easily calculated from the solution of a second order ordinary differential equation in the general case and from an explicit formula when the mean flow is symmetric but the source location is not. This factor is used to assess the potential of employing asymmetric velocity profiles that redirect the sound upward to reduce the noise radiation below the flight path of a jet aircraft.

Goldstein, M. E.↗

Convergence acceleration of viscous flow computations

A multiple-grid convergence acceleration technique introduced for application to the solution of the Euler equations by means of Lax-Wendroff algorithms is extended to treat compressible viscous flow. Computational results are presented for the solution of the thin-layer version of the Navier-Stokes equations using the explicit MacCormack algorithm, accelerated by a convective coarse-grid scheme. Extensions and generalizations are mentioned.

Johnson, G. M.↗

Implicit Total Variation Diminishing (TVD) schemes for steady-state calculations

The application of a new implicit unconditionally stable high resolution total variation diminishing (TVD) scheme to steady state calculations. It is a member of a one parameter family of explicit and implicit second order accurate schemes developed by Harten for the computation of weak solutions of hyperbolic conservation laws. This scheme is guaranteed not to generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments show that this scheme not only has a rapid convergence rate, but also generates a highly resolved approximation to the steady state solution. A detailed implementation of the implicit scheme for the one and two dimensional compressible inviscid equations of gas dynamics is presented. Some numerical computations of one and two dimensional fluid flows containing shocks demonstrate the efficiency and accuracy of this new scheme.

Yee, H. C.↗

Implicit total variation diminishing (TVD) schemes for steady-state calculations

The application of a new implicit unconditionally stable high resolution total variation diminishing (TVD) scheme to steady state calculations. It is a member of a one parameter family of explicit and implicit second order accurate schemes developed by Harten for the computation of weak solutions of hyperbolic conservation laws. This scheme is guaranteed not to generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments show that this scheme not only has a rapid convergence rate, but also generates a highly resolved approximation to the steady state solution. A detailed implementation of the implicit scheme for the one and two dimensional compressible inviscid equations of gas dynamics is presented. Some numerical computations of one and two dimensional fluid flows containing shocks demonstrate the efficiency and accuracy of this new scheme. Previously announced in STAR as N83-23085

Yee, H. C.↗

Investigation of the influence of the zeta-potential on the filtration rate in the presence of collectors

The value of the zeta-potential does not have an explicit effect, which is expressed by a simple math correlation, on filtration rate when a solution of the tested collector is filtered through a cake prepared under standard conditions from the examined particulate material. The zeta-potential measurements and filtration tests were carried out on silica and galena with solutions contg. a cationic container ANP and Et xanthane, resp. at PH = 6.5, varying concentration of the agent (0-2500 g/ton), and under a vacuum of 100 to 600 mm Hg.

Provirnina, E. V.↗

Aerothermal loads analysis for high speed flow over a quilted surface configuration

Attention is given to hypersonic laminar flow over a quilted surface configuration that simulates an array of Space Shuttle Thermal Protection System panels bowed in a spherical shape as a result of thermal gradient through the panel thickness. Pressure and heating loads to the surface are determined. The flow field over the configuration was mathematically modeled by means of time-dependent, three-dimensional conservation of mass, momentum, and energy equations. A boundary mapping technique was then used to obtain a rectangular, parallel piped computational domain, and an explicit MacCormack (1972) explicit time-split predictor corrector finite difference algorithm was used to obtain steady state solutions. Total integrated heating loads vary linearly with bowed height when this value does not exceed the local boundary layer thickness.

Olsen, G. C.↗

A complex T-matrix derivation of a resonance amplitude

Time-dependent perturbation-theory techniques are used to derive a compact expression, valid to arbitrary order and displaying time dependence explicitly, for a quantum-mechanical transition amplitude applicable to the description of resonances. A solution representing a complex-energy generalization of the usual real-energy (nonresonant) amplitude is obtained and generalized to all orders by introducing a complex-energy T matrix. Applications to physical problems such as the extension of the Fermi golden rules to resonances (Norbury and Deutchman, 1984); pion production in relativistic nucleus-nucleus collisions (Norbury et al., 1985); and the determination of resonant cross sections for nuclear, atomic, or molecular processes involving the formation and decay of intermediate discrete or continuum states are indicated.

Norbury, J. W.↗

Probabilistic Structures Analysis Methods (PSAM) for select space propulsion system components

The basic formulation for probabilistic finite element analysis is described and demonstrated on a few sample problems. This formulation is based on iterative perturbation that uses the factorized stiffness on the unperturbed system as the iteration preconditioner for obtaining the solution to the perturbed problem. This approach eliminates the need to compute, store and manipulate explicit partial derivatives of the element matrices and force vector, which not only reduces memory usage considerably, but also greatly simplifies the coding and validation tasks. All aspects for the proposed formulation were combined in a demonstration problem using a simplified model of a curved turbine blade discretized with 48 shell elements, and having random pressure and temperature fields with partial correlation, random uniform thickness, and random stiffness at the root.

Source record↗

Suppression of pressure oscillations in an open cavity by passive pneumatic control

A computational investigation has been conducted to determine the effectiveness of a passive control technique in suppressing the oscillations in an open cavity exposed to a supersonic flow. Time-accurate solutions of the unsteady, Reynolds-averaged, Navier-Stokes equations were obtained with an explicit predictor-corrector algorithm; the passive control was implemented through the use of a simple linear pressure-velocity law along the porous cavity floor. The computational code was validated by comparisons with experimental data for the cavity flow without control. The computational results with control demonstrate that the oscillations are suppressed; the fluid dynamic mechanism of the control is seen to be a stabilization of the free-shear layer. Spectral analysis of the unsteady data shows that the resonant frequencies are essentially unchanged by the control. Beneficial reductions in the cavity resonant drag are also observed.

Chokani, N.↗

Agglomeration Multigrid for an Unstructured-Grid Flow Solver

An agglomeration multigrid scheme has been implemented into the sequential version of the NASA code USM3Dns, tetrahedral cell-centered finite volume Euler/Navier-Stokes flow solver. Efficiency and robustness of the multigrid-enhanced flow solver have been assessed for three configurations assuming an inviscid flow and one configuration assuming a viscous fully turbulent flow. The inviscid studies include a transonic flow over the ONERA M6 wing and a generic business jet with flow-through nacelles and a low subsonic flow over a high-lift trapezoidal wing. The viscous case includes a fully turbulent flow over the RAE 2822 rectangular wing. The multigrid solutions converged with 12%-33% of the Central Processing Unit (CPU) time required by the solutions obtained without multigrid. For all of the inviscid cases, multigrid in conjunction with an explicit time-stepping scheme performed the best with regard to the run time memory and CPU time requirements. However, for the viscous case multigrid had to be used with an implicit backward Euler time-stepping scheme that increased the run time memory requirement by 22% as compared to the run made without multigrid.

Frink, Neal↗

Fusing Quantitative Requirements Analysis with Model-based Systems Engineering

A vision is presented for fusing quantitative requirements analysis with model-based systems engineering. This vision draws upon and combines emergent themes in the engineering milieu. "Requirements engineering" provides means to explicitly represent requirements (both functional and non-functional) as constraints and preferences on acceptable solutions, and emphasizes early-lifecycle review, analysis and verification of design and development plans. "Design by shopping" emphasizes revealing the space of options available from which to choose (without presuming that all selection criteria have previously been elicited), and provides means to make understandable the range of choices and their ramifications. "Model-based engineering" emphasizes the goal of utilizing a formal representation of all aspects of system design, from development through operations, and provides powerful tool suites that support the practical application of these principles. A first step prototype towards this vision is described, embodying the key capabilities. Illustrations, implications, further challenges and opportunities are outlined.

optimization↗