Search NASA⌕ Search

SEARCH · Search NASA

Results for “splitting methods”

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 361 records · Page 20

On Variable Geometric Factor Systems for Top-Hat Electrostatic Space Plasma Analyzers

Even in the relatively small region of space that is the Earth's magnetosphere, ion and electron fluxes can vary by several orders of magnitude. Top-hat electrostatic analyzers currently do not possess the dynamic range required to sample plasma under all conditions. The purpose of this study was to compare, through computer simulation, three new electrostatic methods that would allow the sensitivity of a sensor to be varied through control of its geometric factor (GF) (much like an aperture on a camera). The methods studied were inner filter plates, split hemispherical analyzer (SHA) and top-cap electrode. This is the first discussion of the filter plate concept and also the first study where all three systems are studied within a common analyzer design, so that their relative merits could be fairly compared. Filter plates were found to have the important advantage that they facilitate the reduction in instrument sensitivity whilst keeping all other instrument parameters constant. However, it was discovered that filter plates have numerous disadvantages that make such a system impracticable for a top-hat electrostatic analyzer. It was found that both the top-cap electrode and SHA are promising variable geometric factor system (VGFS) concepts for implementation into a top-hat electrostatic analyzer, each with distinct advantages over the other.

Kataria, Dhiren O.↗

Resistive Magnetohydrodynamic Simulations of Relativistic Magnetic Reconnection

Resistive relativistic magnetohydrodynamic (RRMHD) simulations are applied to investigate the system evolution of relativistic magnetic reconnection. A time-split Harten-Lan-van Leer method is employed. Under a localized resistivity, the system exhibits a fast reconnection jet with an Alfv enic Lorentz factor inside a narrow Petschek-type exhaust. Various shock structures are resolved in and around the plasmoid such as the post-plasmoid vertical shocks and the "diamond-chain" structure due to multiple shock reflections. Under a uniform resistivity, Sweet-Parker-type reconnection slowly evolves. Under a current-dependent resistivity, plasmoids are repeatedly formed in an elongated current sheet. It is concluded that the resistivity model is of critical importance for RRMHD modeling of relativistic magnetic reconnection.

Zenitani, Seiji↗

Probabilistic Forecasting of Ground Magnetic Perturbation Spikes at Mid-Latitude Stations

The prediction of large fluctuations in the ground magnetic field (dB/dt) is essential for preventing damage from Geomagnetically Induced Currents. Directly forecasting these fluctuations has proven difficult, but accurately determining the risk of extreme events can allow for the worst of the damage to be prevented. Here we trained Convolutional Neural Network models for eight mid-latitude magnetometers to predict the probability that dB/dt will exceed the 99th percentile threshold 30–60 min in the future. Two model frameworks were compared, a model trained using solar wind data from the Advanced Composition Explorer (ACE) satellite, and another model trained on both ACE and SuperMAG ground magnetometer data. The models were compared to examine if the addition of current ground magnetometer data significantly improved the forecasts of dB/dt in the future prediction window. A bootstrapping method was employed using a random split of the training and validation data to provide a measure of uncertainty in model predictions. The models were evaluated on the ground truth data during eight geomagnetic storms and a suite of evaluation metrics are presented. The models were also compared to a persistence model to ensure that the model using both datasets did not over-rely on dB/dt values in making its predictions. Overall, we find that the models using both the solar wind and ground magnetometer data had better metric scores than the solar wind only and persistence models, and was able to capture more spatially localized variations in the dB/dt threshold crossings.

Michael Coughlan↗

Methods and systems for coordinating predictive cruise control, engine-off coasting, and hybrid power split

The enclosed disclosure relates to hybrid vehicles and systems with an engine, a drivetrain with a clutch and a transmission, an electric machine, and a controller. The controller receives lookahead information within a lookahead window and present state information of the hybrid vehicle. The controller determines a predicted coasting opportunity exceeding a predetermined threshold within the lookahead window and determines a cruise control reference speed, a power split between the engine and the electric machine, and a timing of enabling engine-off coasting during the coasting opportunity. The controller deactivates the engine and disengages the clutch at a start of the coasting opportunity when the engine-off coasting is enabled.

Gupta, Rohinish↗

Sensitivity of Simulations of Double-detonation Type Ia Supernovae to Integration Methodology

Abstract We study the coupling of hydrodynamics and reactions in simulations of the double-detonation model for Type Ia supernovae. When assessing the convergence of simulations, the focus is usually on spatial resolution; however, the method of coupling the physics together as well as the tolerances used in integrating a reaction network also play an important role. In this paper, we explore how the choices made in both coupling and integrating the reaction portion of a simulation (operator/Strang splitting versus the simplified spectral deferred corrections method we introduced previously) influences the accuracy, efficiency, and nucleosynthesis of simulations of double detonations. We find no need to limit reaction rates or reduce the simulation time step to the reaction timescale. The entire simulation methodology used here is GPU-accelerated and made freely available as part of the Castro simulation code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Toward Accurate Spin–Orbit Splittings from Relativistic Multireference Electronic Structure Theory

Most nonrelativistic electron correlation methods can be adapted to account for relativistic effects, as long as the relativistic molecular spinor integrals are available, from either a four-, two-, or one-component mean-field calculation. Furthermore, relativistic multireference correlation methods remain a relatively unexplored area, with mixed evidence regarding the improvements brought by perturbative treatments. We report, for the first time, the implementation of state-averaged four-component relativistic multireference perturbation theories to second and third order based on the driven similarity renormalization group (DSRG). With our methods, named 4c-SA-DSRG-MRPT2 and 3, we find that the dynamical correlation included on top of 4c-CASSCF references can significantly improve the spin-orbit splittings in p-block elements and potential energy surfaces when compared to 4c-CASSCF and 4c-CASPT2 results. We further show that 4c-DSRG-MRPT2 and 3 are applicable to these systems over a wide range of the flow parameter, with systematic improvement from second to third order in terms of both improved error statistics and reduced sensitivity with respect to the flow parameter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Survey of computational methods for three-dimensional supersonic inviscid flows with shocks

Hyperbolic systems of partial differential equations governing supersonic inviscid flows are discussed and analyzed. Finite-difference analogues for integrating these systems in the interior of fluid domains are described from two points of view: a differential form approach and an integral form approach. The algorithms presented are analyzed for stability and accuracy. The concept of time splitting is discussed and supplied to these methods to achieve increased numerical efficiency. Techniques for treating conditions at the boundaries of the fluid domain and shock-wave discontinuities at surfaces within the domain are described.

Maccormack, R. W.↗

Flux-difference split scheme for turbulent transport equations

This paper describes a high accuracy upwind method to solve the mean compressible flow equations closed with two-equation turbulence models. A flux-difference splitting scheme is developed to compute the flowfield without introducing oscillations into the solution. A high-order accuracy scheme is used to predict both the mean flow variables and the turbulence variables. Results comparing two different two-equation turbulence models are presented for supersonic flow over a flat plate and a compressible, free shear-layer. The non-oscillatory behavior of the upwind scheme is demonstrated on the free shear-layer.

Morrison, J. H.↗

New Time Integrators and Capabilities in SUNDIALS Versions 6.2.0-7.4.0

SUNDIALS is a well-established numerical library that provides robust and efficient time integrators and nonlinear solvers. This article overviews several significant improvements and new features added over the last 3 years to support scientific simulations run on high-performance computing systems. Notably, three new classes of one-step methods have been implemented: low storage Runge–Kutta, symplectic partitioned Runge–Kutta, and operator splitting. In addition, we describe new timestep adaptivity support for multirate methods, adjoint sensitivity analysis capabilities for explicit Runge–Kutta methods, additional options for Anderson acceleration in nonlinear solvers, and improved error handling and logging.

Computer science↗

Evaluating Corn (Zea Mays L.) N Variability Via Remote Sensed Data

Transformations and losses of nitrogen (N) throughout the growing season can be costly. Methods in place to improve N management and facilitate split N applications during the growing season can be time consuming and logistically difficult. Remote sensing (RS) may be a method to rapidly assess temporal changes in crop N status and promote more efficient N management. This study was designed to evaluate the ability of three different RS platforms to predict N variability in corn (Zea mays L.) leaves during vegetative and early reproductive growth stages. Plots (15 x 15m) were established in the Coastal Plain (CP) and Appalachian Plateau (AP) physiographic regions each spring from 2000 to 2002 in a completely randomized design. Treatments consisted of four N rates (0, 56, 112, and 168 kg N/ha) applied as ammonium nitrate (NH4N03) replicated four time. Spectral measurements were acquired via spectroradiometer (lambda = 350 - 1050 nm), Airborne Terrestrial Applications Sensor (ATLAS) (lambda = 400 - 12,500 nm), and the IKONOS satellite (lambda = 450 - 900 nm). Spectroradiometer data were collected on a biweekly basis from V4 through R1. Due to the nature of - satellite and aircraft acquisitions, these data were acquired per availability. Chlorophyll meter (SPAD) and tissue N were collected as ancillary data along with each RS acquisition. Results showed vegetation indices derived from hand-held spectroradiometer measurements as early as V6-V8 were linearly related to yield and tissue N content. ATLAS data was correlated with tissue N at the AP site during the V6 stage (r2 = 0.66), but no significant relationships were observed at the CP site. No significant relationships were observed between plant N and IKONOS imagery. Using a combination of the greenness vegetation index (GNDVI) and the normalized difference vegetation index (NDVI), RS data acquired via ATLAS and the spectroradiometer could be used to evaluate tissue N variability and estimate corn yield variability under ideal growing conditions.

Sullivan, D. G.↗

An analysis of flux-split algorithms for Euler's equations with real gases

An analysis of flux-splitting procedures for the solution of Euler's equations with real gas effects is presented. An alternative real-gas flux-splitting is derived which can easily be implemented into existing codes. This approach, which takes the form of an 'equivalent' gamma representation is not an ad hoc model, but is based on theoretical considerations. Details of this method with the Steger-Warming and Van Leer flux vector splittings and the Roe flux-difference splitting are given. Applications of the method to several high Mach number, high temperature flows are presented for one and two space dimensions.

Grossman, B.↗

Accurate upwind methods for the Euler equations

A new class of piecewise linear methods for the numerical solution of the one-dimensional Euler equations of gas dynamics is presented. These methods are uniformly second-order accurate, and can be considered as extensions of Godunov's scheme. With an appropriate definition of monotonicity preservation for the case of linear convection, it can be shown that they preserve monotonicity. Similar to Van Leer's MUSCL scheme, they consist of two key steps: a reconstruction step followed by an upwind step. For the reconstruction step, a monotonicity constraint that preserves uniform second-order accuracy is introduced. Computational efficiency is enhanced by devising a criterion that detects the 'smooth' part of the data where the constraint is redundant. The concept and coding of the constraint are simplified by the use of the median function. A slope steepening technique, which has no effect at smooth regions and can resolve a contact discontinuity in four cells, is described. As for the upwind step, existing and new methods are applied in a manner slightly different from those in the literature. These methods are derived by approximating the Euler equations via linearization and diagonalization. At a 'smooth' interface, Harten, Lax, and Van Leer's one intermediate state model is employed. A modification for this model that can resolve contact discontinuities is presented. Near a discontinuity, either this modified model or a more accurate one, namely, Roe's flux-difference splitting. is used. The current presentation of Roe's method, via the conceptually simple flux-vector splitting, not only establishes a connection between the two splittings, but also leads to an admissibility correction with no conditional statement, and an efficient approximation to Osher's approximate Riemann solver. These reconstruction and upwind steps result in schemes that are uniformly second-order accurate and economical at smooth regions, and yield high resolution at discontinuities.

Huynh, Hung T.↗

Manifestation of the Fermi resonance in surface polariton spectra

The method of disturbed full internal reflection (DFIR) is used to detect and interpret the resonance splitting of the surface polariton. The effect in the spectra of oscillatory SP which are reflected by the DFIR method in Otto geometry was experimentally recorded. It is concluded that the resonance splitting of the dispersion branch of SP may serve as an effective method for detecting weak oscillations and for measuring their parameters.

Puchkovskaya, G. A.↗

Simulation of a Periodic Jet in a Crossflow with a RANS Solver Using an Unstructured Grid

A second-order unstructured-grid code, developed and used primarily for steady aerodynamic simulations, is applied to the synthetic jet in a cross flow. The code, FUN3D, is a vertex-centered finite-volume method originally developed by Anderson[1, 2], and is currently supported by members of the Fast Adaptive Aerospace Tools team at NASA Langley. Used primarily for design[3] and analysis[4] of steady aerodynamic configurations, FUN3D incorporates a discrete adjoint capability, and supports parallel computations using MPI. A detailed description of the FUN3D code can be found in the references given above. The code is under continuous development and contains a variety of flux splitting algorithms for the inviscid terms, two methods for computing gradients, several turbulence models, and several solution methodologies; all in varying states of development. Only the most robust and reliable components, based on experiences with steady aerodynamic simulations, were employed in this work. As applied in this work, FUN3D solves the Reynolds averaged Navier-Stokes equations using the one equation turbulence model of Spalart and Allmaras[5]. The spatial discretization is formed on unstructured meshes using a vertex-centered approach. The inviscid terms are evaluated by a flux-difference splitting formulation using least-squares reconstruction and Roe-type approximate Riemann fluxes. Green-Gauss gradient evaluations are used for viscous and turbulence modeling terms. The discrete spatial operator is combined with a backward time operator which is then solved iteratively using point or line Gauss-Seidel and local time stepping in a pseudo time. For steady flows, the physical time step is set to infinity and the pseudo time step is ramped up with the iteration count. A second-order backward in time operator is used for time accurate flows with 20 to 50 steps in the pseudo time applied at each physical time step. For this effort, FUN3D was modified to support spatially varying boundary and initial conditions, and unsteady boundary conditions. Also, a specialized in/out flow boundary condition was implemented to model the action of the diaphragm. This boundary condition is described below in more detail. The grids were generated using the internally developed codes GridEX[6] for meshing the surfaces and inviscid regions of the domain, and for CAD access; and MesherX[7] for meshing the viscous regions. Grid spacing in on the surfaces and in the inviscid regions are indirectly controlled by specifying sources. The viscous layers are generated using an advancing layer technique. MeshersX allows the user to control the spatial variation of the first step off the surface, growth rates, and the termination criterion by providing small problem dependent subroutines.

Atkins, H. L.↗

A pressure flux-split technique for computation of inlet flow behavior

A method for calculating the flow field in aircraft engine inlets is presented. The phenomena of inlet unstart and restart are investigated. Solutions of the reduced Navier-Stokes (RNS) equations are obtained with a time consistent direct sparse matrix solver that computes the transient flow field both internal and external to the inlet. Time varying shocks and time varying recirculation regions can be efficiently analyzed. The code is quite general and is suitable for the computation of flow for a wide variety of geometries and over a wide range of Mach and Reynolds numbers.

Pordal, H. S.↗

Electronic circuit for control or coupling of single charges or spins and methods therefor

A quantum dot structure having a split-gate geometry is provided. The quantum dot is configured for incorporation into a quantum dot array of a quantum processing unit. A gap between a reservoir accumulation gate and a quantum dot accumulation gate provides a tunnel barrier between an electric charge reservoir and a quantum dot well. An electrical potential applied to the gates defines a tunnel barrier height, width and charge tunneling rate between the well and the reservoir without relying on any barrier gate to control the charge tunneling rate.

Pioro-Ladriere, Michel↗

Comparison of CALIPSO-Like, LaRC, and MODIS Retrievals of Ice Cloud Properties over SIRTA in France and Florida during CRYSTAL-FACE

This study compares cirrus particle effective radius retrieved by a CALIPSO-like method with two similar methods using MODIS, MODI Airborne Simulator (MAS), and GOES imagery. The CALIPSO-like method uses lidar measurements coupled with the split-window technique that uses the infrared spectral information contained at the 8.65-micrometer, 11.15-micrometer and 12.05-micrometer bands to infer the microphysical properties of cirrus clouds. The two other methods, sing passive remote sensing at visible and infrared wavelengths, are the operational MODIS cloud products (referred to by its archival product identifier MOD06 for MODIS Terra) and MODIS retrievals performed by the CERES team at LaRC (Langley Research Center) in support of CERES algorithms; the two algorithms will be referred to as MOD06- and LaRC-method, respectively. The three techniques are compared at two different latitudes: (i) the mid-latitude ice clouds study uses 18 days of observations at the Palaiseau ground-based site in France (SIRTA: Site Instrumental de Recherche par Teledetection Atmospherique) including a ground-based 532 nm lidar and the Moderate Resolution Imaging Spectrometer (MODIS) overpasses on the Terra Platform, (ii) the tropical ice clouds study uses 14 different flight legs of observations collected in Florida, during the intensive field experiment CRYSTAL-FACE (Cirrus Regional Study of Tropical Anvils and cirrus Layers-Florida Area Cirrus Experiment), including the airborne Cloud Physics Lidar (CPL) and the MAS. The comparison of the three methods gives consistent results for the particle effective radius and the optical thickness, but discrepancies in cloud detection and altitudes. The study confirms the value of an active remote-sensing method (CALIPSO-like) for the study of sub-visible ice clouds, in both mid-latitudes and tropics. Nevertheless, this method is not reliable in optically very thick tropical ice clouds.

Chiriaco, M.↗

Experimental Studies on Wood Pellets Combustion in a Fixed Bed Combustor Using Taguchi Method

The combustion of wood pellets in a fixed bed combustor of a 20 kW capacity domestic pellet boiler was tested according to several factors including Power, Excess Air (EA), Primary/Secondary air Split Ratio (SR) and Grate Area (GA). The Taguchi method was applied to program the experimental design. Several parameters were measured, including gas emissions (CO), fuel bed temperature (measured at 4 different heights), and efficiency. The experimental results show that the lower CO emission and the higher efficiency were obtained at medium thermal loads and the highest temperature on the fuel bed was obtained at about ¼ of its height (15 mm). The results obtained from the analysis of variance (ANOVA) show that the SR and the Power are the most important factors contributing to the CO reduction and also increase the fuel bed temperature.

Castro, Carlos↗