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 343 records · Page 19

Accuracy and stability of time-split finite-difference schemes

In a recently published work by Abarbanel and Gottlieb (1980), a new class of explicit time-split algorithms designed for application to the compressible Navier-Stokes equations was developed. These algorithms, which utilize locally-one-dimensional (LOD) spatial steps, were shown to possess stability characteristics superior to those of other time-split schemes. In the present work, the properties of an implicit LOD method, analogous to the Abarbanel-Gottlieb algorithm, are examined using the two-dimensional heat conduction equation as the test problem. Both temporal and spatial inconsistencies inherent in the scheme are identified, and a new consistent, implicit splitting approach is developed and applied to the linear Burgers' equation. The relationship between this new method and other time-split implicit schemes is explained and stability problems encountered with the method in three dimensions are discussed.

Dwoyer, D. L.↗

Modeling and Feedback Compensator Design for Power-Hardware-in-the-Loop System for Medium-Voltage Grid-Connected Power Converters

Hardware-in-the-loop (HIL) evaluation is a method in which the test subject is split into a physical part and a simulated part, and these parts are connected with interfaces to form a combined physical-numerical system. Power-hardware-in-the-loop (PHIL) systems, where actual hardware is connected to a real-time model, can exhibit issues of instability, inaccuracy and errors when operating in a closed-loop. The challenges of generating a PHIL setup are that, first, because of the limited dynamic response of the different parts of the system, the test results might be inaccurate, and, second, because of the high frequency noise introduced by the sensors to the closed-loop system, it can be difficult to design a compensator for the real-time emulator response, while stabilizing the closed-loop system at the same time. In this paper, different parts of a PHIL system are characterized, and the feedback compensator system design is proposed for the stable operation of the closed-loop PHIL system. The issues as observed in a PHIL system without any compensator are demonstrated using experimental results, and the effectiveness of a first order phase lead compensator is validated.

compensator design↗

Unbalanced and Minimal Point Equivalent Estimation Second-Order Split-Plot Designs

Restricting the randomization of hard-to-change factors in industrial experiments is often performed by employing a split-plot design structure. From an economic perspective, these designs minimize the experimental cost by reducing the number of resets of the hard-to- change factors. In this paper, unbalanced designs are considered for cases where the subplots are relatively expensive and the experimental apparatus accommodates an unequal number of runs per whole-plot. We provide construction methods for unbalanced second-order split- plot designs that possess the equivalence estimation optimality property, providing best linear unbiased estimates of the parameters; independent of the variance components. Unbalanced versions of the central composite and Box-Behnken designs are developed. For cases where the subplot cost approaches the whole-plot cost, minimal point designs are proposed and illustrated with a split-plot Notz design.

Parker, Peter A.↗

Split-Step Simulations of Sonic Boom Propagation Beyond the Lateral Cutoff in a Turbulent Atmosphere

Recent flight tests during the Quiet Supersonic Flights 2018 (QSF18) study reported sonic booms heard outside of the primary carpet region. In the absence of turbulence, the lateral cutoff region separates the primary sonic boom carpet from the shadow zone, where the sonic boom signal experiences significant attenuation. However, when turbulence is present in the atmospheric boundary layer (ABL), additional scattering of the sonic boom to the shadow zone region occurs. A method is presented for simulating sonic boom propagation in a turbulent atmospheric boundary layer beyond the lateral cutoff region into the shadow zone. A split-step method is used to integrate a partially one-way equation for the acoustic pressure. Inhomogeneous turbulence, representative of the ABL, is generated in the computational domain with a Fourier synthesis approach. Distributions of several loudness metrics in the shadow zone region for a sonic boom N-wave and a shaped boom are examined. Increasing both turbulence root-mean-square velocity and integral length scale are found to increase the average loudness of booms in the shadow zone. (This research is supported by the Commercial Supersonic Technology Project of the National Aeronautics and Space Administration under Grant No. 80NSSC19K1685.)

sonic boom↗

A quasi-static method for determining the characteristics of a motion capture camera system in a "split-volume" configuration

The purpose of this study was to determine the accuracy, repeatability and resolution of a six-camera Motion Analysis system in a vertical split-volume configuration using a unique quasi-static methodology. The position of a reflective marker was recorded while it was moved quasi-statically over a range of 2.54 mm (0.100 inches) via a linearly-translating table. The table was placed at five different heights to cover both sub-volumes and the overlapping region. Data analysis showed that accuracy, repeatability and resolution values were consistent across all regions of the split-volume, including the overlapping section.

NASA Center JSC↗

A Quasi-Static Method for Determining the Characteristics of a Motion Capture Camera System in a "Split-Volume" Configuration

To confidently report any data collected from a video-based motion capture system, its functional characteristics must be determined, namely accuracy, repeatability and resolution. Many researchers have examined these characteristics with motion capture systems, but they used only two cameras, positioned 90 degrees to each other. Everaert used 4 cameras, but all were aligned along major axes (two in x, one in y and z). Richards compared the characteristics of different commercially available systems set-up in practical configurations, but all cameras viewed a single calibration volume. The purpose of this study was to determine the accuracy, repeatability and resolution of a 6-camera Motion Analysis system in a split-volume configuration using a quasistatic methodology.

Miller, Chris↗

System and method for sub-wavelength detection for jetting-based additive manufacturing using a split ring resonator probe

The present disclosure relates to a system for detecting and analyzing droplets of feedstock material being ejected from an additive manufacturing device. The system makes use of a split ring resonator (SRR) probe including a ring element having a gap, with the gap being positioned adjacent a path of travel of the droplets of feedstock material. An excitation signal source is used for supplying an excitation signal to the SRR probe. An analyzer analyzes signals generated by the SRR probe in response to perturbations in an electric field generated by the SRR probe as the droplets of feedstock material pass the ring element. The signals are indicative of dimensions of the droplets of feedstock material.

Mukherjee, Saptarshi↗

In-orbit performance of the ITOS improved attitude control system with Hall generator brushless motor and earth-splitting technique

The National Aeronautics and Space Administration (NASA), launched ITOS-D with an improved attitude control system. A Hall generator brushless dc torque motor replaced the brush dc torque motor on Tiros-M and ITOS-A. Two CO2 attitude horizon sensors and one mirror replaced the four wideband horizon sensors and two mirrors on ITOS-1 and NOAA-1. Redundant pitch-control electronic boxes containing additional electronic circuitry for earth-splitting and brushless motor electronics were used. A method of generating a spacecraft earth-facing side reference for comparison to the time occurrence of the earth-splitting pulse was used to automatically correct pitch-attitude error. A single rotating flywheel, supported by a single bearing, provided gyroscopic stability and the required momentum interchange to keep one side of the satellite facing the earth. Magnetic torquing against the earth's magnetic field eliminated the requirement for expendable propellants which would limit satellite life in orbit.

Peacock, W. M.↗

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition

In a previous paper Christlieb et al. (A particle-in-cell method for plasmas with a generalized momentum formulation, part I: Model formulation, 2024), we developed a new particle-in-cell (PIC) method for the relativistic Vlasov–Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This new method evolved the potentials under the Lorenz gauge using integral equation methods. New methods to construct spatial derivatives of the potentials that converge at the same rates as the fields were also presented. The new particle method was compared against standard explicit discretizations, including the well-known FDTD-PIC method, for a range of applications involving sheaths and particle beams. Here, this paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell’s equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. Both the finite difference and spectral implementations satisfy this discrete gauge condition to machine precision. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The potential benefit of this approach is that it eliminates spatial derivatives that appear on the particle data, namely the current density, which is often calculated by linear combinations of low-order spline basis functions. This method is ideally suited to boundary integral equation methods that invert multi-dimensional operators without dimensional splitting techniques and will be the subject of future work. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. This results in a gauge error, as the maps do not enforce the continuity equation. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. This method also enforces the gauge condition to machine precision. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in a non-staggered PIC formulation.

97 MATHEMATICS AND COMPUTING↗

Mersen Grade 2114: A Comparison of Tensile Strength Data

The ASTM Brazilian disc graphite strength (splitting tensile strength, σ sts ) test method (ASTM D8289) is of interest because the small-specimen geometry is compatible with that of environmental effects specimens, thus allowing for environmental effects such as irradiation, irradiation creep, or oxidation on tensile strength to be investigated. The Brazilian disc strength of Mersen grade 2114 graphite is reported and compared with strength data previously obtained using larger cylindrical ASTM dog-bone specimens.

36 MATERIALS SCIENCE↗

Refining Principal Stress Measurements in Reservoir Underburden in Regions of Induced Seismicity through Seismological Tools, Laboratory Experiments - Final Technical Report

This project developed methodologies to measure the in-situ principal stress in the deep subsurface through use of multiple independent, but complementary, seismic methods, laboratory verification, and development of theoretical frameworks. By leveraging existing regional and local datasets we developed, tested, and refined a set of diagnostic tools for determining the in-situ stress state with reduced uncertainty at and below reservoir depths (1.5-6 km). A set of novel tools was produced that are scale independent, such that their utility is equivalent on regional, field scale, and near borehole monitoring of principal stresses in reservoir underburden for carbon storage projects. During a 4-year Department of Energy (DOE) and Southern Company funded project, carried out by the Electric Power Research Institute (EPRI), Lawrence Livermore National Laboratory (LLNL), the University of Oklahoma (OU), and the U.S. Geological Survey (USGS), the project team developed methodologies to measure the far-field in-situ principal stress in the deep subsurface, leveraging induced seismicity data from waste-water disposal projects. These methodologies consisted in the use of well-established and technically advanced seismic processing methods, such as virtual seismometer method-moment tensor (VSM-MT) and shear wave splitting (SWS), that are adept at recovering the stress orientation and certain components of the stress tensor. These methods were applied to robust seismicity catalogs created with matched filter techniques near sites of active fluid disposal—a proxy for carbon storage sites where such datasets are more limited. Estimates of the stress orientation made with seismic processing tools were considered along with laboratory acoustic emission experiments conducted on rock samples from the region of interest. Stress orientations in the studied region do not vary significantly across distances of ~100 km, nor are they found to rotate through time as a consequence of local wastewater disposal, as previously speculated. Finally, the project team investigated the trade-offs among the different seismic methods and evaluated the range of uncertainty that is generated with these methodologies, which led to a practical use and refinement of the VSM-MT technique when it is applied to field datasets. Understanding the trade-offs between these different methods highlighted the potential benefits of improved quantification of uncertainties on stress field estimations.

58 GEOSCIENCES↗

Mersen Grade 2114: Tensile Strength from Brazilian Disc Testing

The ASTM Brazilian disc graphite strength (splitting tensile strength, σ sts ) test method (ASTM D8289) is of interest because the small-specimen geometry is compatible with that of environmental effects specimens, thus allowing for environmental effects such as irradiation, irradiation creep, or oxidation on tensile strength to be investigated. The Brazilian disc strength of Mersen grade 2114 graphite (billet 116310 slab 1) is reported and compared with strength data previously obtained using larger cylindrical ASTM dog-bone specimens and with Brazilian disc strength data obtained from billet 116310 slab 5.

36 MATERIALS SCIENCE↗

Applications of partial orderings to the study of positive definiteness, monotonicity, and convergence of iterative methods for linear systems.

Most of the known results concerning convergence of iterative methods for solving linear systems involve either positive definiteness or monotonicity. In this paper a new concept, called K-semipositivity, is introduced, which provides a link between convergence theory, monotonicity, and positive definiteness. By using this concept, together with partial orderings in Euclidean n-space, several new convergence theorems are proved. Application to Jacobi's methods and the theory of regular splittings shows the usefulness of these new results.

Vandergraft, J. S.↗

Adaptive grid embedding for the two-dimensional flux-split Euler equations

A numerical algorithm is presented for solving the 2-D flux-split Euler equations using a multigrid method with adaptive grid embedding. The method uses an unstructured data set along with a system of pointers for communication on the irregularly shaped grid topologies. An explicit two-stage time advancement scheme is implemented. A multigrid algorithm is used to provide grid level communication and to accelerate the convergence of the solution to steady state. Results are presented for a subcritical airfoil and a transonic airfoil with 3 levels of adaptation. Comparisons are made with a structured upwind Euler code which uses the same flux integration techniques of the present algorithm. Good agreement is obtained with converged surface pressure coefficients. The lift coefficients of the adaptive code are within 2 1/2 percent of the structured code for the sub-critical case and within 4 1/2 percent of the structured code for the transonic case using approximately one-third the number of grid points.

Warren, Gary Patrick↗

Adaptive grid embedding for the two-dimensional Euler equations

A numerical algorithm is presented for solving the two-dimensional flux-split Euler equations using a multigrid method with adaptive grid embedding. The method uses an unstructured data set along with a system of pointers for communication on the irregularly shaped grid topologies. An explicit two-stage time advancement scheme is implemented. A multigrid algorithm is used to provide grid level communication and to accelerate the convergence of the solution to steady state. Results are presented for an NACA 0012 airfoil in a freestream with Mach numbers of 0.95 and 1.054. Excellent resolution of the shock structures is obtained with the adaptive grid embedding method with significantly fewer grid points than the comparable structured grid.

Warren, Gary P.↗

Three-dimensional shock-wave/boundary-layer interactions with bleed through a circular hole

Computations were performed to study three-dimensional (3-D), shock-wave/boundary-layer interactions on a flat plate in which fluid in the boundary layer was bled through a circular hole into a plenum to control shock-wave induced separation. Results are presented which show the details of the 3-D flowfield about the bleed hole and how bleed-hole placement relative to shock-wave impingement affect upstream, spanwise, and downstream influence lengths. This study revealed an underlying mechanisms by which bleed holes can affect shock-wave/boundary-layer interactions. This investigation is based on the ensemble-averaged, "full-compressible" Navier-Stokes equations closed by the Baldwin-Lomax turbulence model. Solutions to these equations were obtained by an implicit finite-volume method based on the partially-split, two-factored algorithm of Steger.

Rimlinger, M. J.↗

Application of multigrid and adaptive grid embedding to the two-dimensional flux-split Euler equations

A numerical algorithm is presented for solving the two-dimensional flux-split Euler equations using a multigrid method with adaptive grid embedding. The method uses an unstructured data set along with a system of pointers for communication on the irregularly shaped grid topologies. An explicit two-stage time-advancement scheme is implemented. A multigrid algorithm is used to provide grid level communication and to accelerate the convergence of the solution to steady state. Results are presented for a NACA 0012 aerofoil in a free stream with a Mach number of 0.85 and an angle of attack of 1.0 degree. Excellent resolution of the shock structures is obtained with the adaptive grid embedding method with significantly fewer grid points than the comparable structured grid.

Warren, Gary P.↗

Three-dimensional fully spectral numerical method for mantle convection with depth-dependent properties

A semi-implicit fully spectral collocation method for the simulation of three-dimensional mantle convection with depth-dependent thermo-dynamic and transport properties is presented. The variable property Navier-Stokes equation expressed in terms of the primitive variable velocity and pressure is solved with the mass continuity and temperature equations. The periodic horizontal boundary conditions allow a Fourier expansion for the two horizontal directions. The stress-free, impermeable isothermal boundary conditions along with the depth dependent coefficients are handled with a Chebyshev expansion in the vertical direction. In the limit of an infinite Prandtl number appropriate to mantle convection, the inertial terms in the momentum equation are unimportant. In this case an explicit solution of a Poisson equation for pressure can be avoided; instead a fourth-order equation for vertical velocity can be solved. Simultaneous imposition of both impermeable and continuity boundary conditions during the vertical velocity evaluation is discussed. The pressure distributions on the top and bottom bounding planes were determined by means of an influence matrix technique. The numerical method employed here avoids time-splitting errors and enforces velocity boundary conditions and continuity over the entire domain, including the boundaries, to machine accuracy. Strongly time-dependent three-dimensional solutions up to a surface Rayleigh number of 1 x 10(exp 7) have been obtained. Strong upwellings, pulsating chaotically, are formed by the collective merging of cylindrical plumes.

Balachandar, S.↗