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 1,135 records · Page 63

Aerodynamics of loaded cascades in subsonic flows subject to unsteady three-dimensional vortical disturbances

A highly efficient numerical method is developed for three-dimensional, periodic, vortical flows around a cascade of loaded airfoils. The method linearizes Euler's equations about the mean flow of the cascade and thus fully accounts for the effects of distortion of the vortical disturbances as they propogate and interact with the cascade mean flow. The numerical scheme is based on splitting the unsteady velocity into vortical and potential parts. The latter is governed by a non-constant coefficient inhomogeneous convective wave equation. A new and computationally suitable out-flow conditions are derived and avoid the difficulties associated with the singular velocity downstream. Solutions were obtained in the frequency domain by using a body-fitted coordinate system. Results are presented to demonstrate the effects of the out-flow boundary conditions, cascade spacing, mean blade loading and gust upstream conditions on the aerodynamics response and unsteady pressure field of a cascade.

Fang, J.↗

Fluid control structures in microfluidic devices

Methods and apparatus for implementing microfluidic analysis devices are provided. A monolithic elastomer membrane associated with an integrated pneumatic manifold allows the placement and actuation of a variety of fluid control structures, such as structures for pumping, isolating, mixing, routing, merging, splitting, preparing, and storing volumes of fluid. The fluid control structures can be used to implement a variety of sample introduction, preparation, processing, and storage techniques.

Mathies, Richard A.↗

Testing of ACC-6 Carbon-Carbon at High Strain Rates and High Temperatures

Some materials, such as steel, titanium, and tungsten, have a yield stress that depends on the strain rate, and actually are stronger at high rates of load or deformation. The increase in strength is relevant; and if not taken into account, the computer models will not match the experimental behavior observed in typical shock loading, ordnance, or ballistic tests. For higher velocity (hypervelocity impacts), an equation of state that accounts for non-linearities and temperature in the pressure-volume relation, may also be needed. This rate-dependent behavior of the strength has also been observed in some composite materials. Southwest Research Institute® (SwRI®) has started an internal research effort to measure and better model the behavior of materials, in particular carbon/carbon composites like ACC-6, at ultra-high temperatures for both static and high rates of deformation. The ultimate objective is to obtain reliable computer models of carbon/carbon material under hypersonic flight conditions, hypervelocity impacts, shock-loading situations such as booster-vehicle separation, and munition operation events. SwRI® has successfully modeled in the past ballistic impact on other composites, namely S-2 glass/epoxy, Kevlar/resin, polyethylene fibers/resin, and carbon fiber composites. This presentation will show test results relevant to the strength properties of the ACC-6 C/C material which were measured at quasistatic, medium, and the moderately high strain rates (~103 s—1) achievable in a conventional Split Hopkinson bar apparatus. Some of the tests were performed at both high-rates and high temperatures typical of hypersonic flight, something that, to the authors’ knowledge, has never been published in the past for this particular material. The test method will be presented, as well as the comparison of the stress-strain curves for the different temperatures and strain-rates of interest. A discussion of how the results can be used in the typical hydrocodes of interest will be included.

carbon/carbon↗

Fluid control structures in microfluidic devices

Methods and apparatus for implementing microfluidic analysis devices are provided. A monolithic elastomer membrane associated with an integrated pneumatic manifold allows the placement and actuation of a variety of fluid control structures, such as structures for pumping, isolating, mixing, routing, merging, splitting, preparing, and storing volumes of fluid. The fluid control structures can be used to implement a variety of sample introduction, preparation, processing, and storage techniques.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fluid control structures in microfluidic devices

Methods and apparatus for implementing microfluidic analysis devices are provided. A monolithic elastomer membrane associated with an integrated pneumatic manifold allows the placement and actuation of a variety of fluid control structures, such as structures for pumping, isolating, mixing, routing, merging, splitting, preparing, and storing volumes of fluid. The fluid control structures can be used to implement a variety of sample introduction, preparation, processing, and storage techniques.

Mathies, Richard A.↗

Method and Apparatus for Accurately Calibrating a Spectrometer

A calibration assembly for a spectrometer is provided. The assembly includes a spectrometer having n detector elements, where each detector element is assigned a predetermined wavelength value. A first source emitting first radiation is used to calibrate the spectrometer. A device is placed in the path of the first radiation to split the first radiation into a first beam and a second beam. The assembly is configured so that one of the first and second beams travels a path-difference distance longer than the other of the first and second beams. An output signal is generated by the spectrometer when the first and second beams enter the spectrometer. The assembly includes a controller operable for processing the output signal and adapted to calculate correction factors for the respective predetermined wavelength values assigned to each detector element.

Youngquist, Robert C.↗

The Uniform Methods Project: Methods for Determining Energy Efficiency Savings for Specific Measures: January 2012 - September 2016

This is the complete set of Uniform Methods Project protocols for determining energy efficiency savings for specific measures, all within one report. The report includes the following chapters: Chapter 1: Introduction. Chapter 2: Commercial and Industrial Lighting Evaluation Protocol. Chapter 3: Commercial and Industrial Lighting Controls Evaluation Protocol. Chapter 4: Small Commercial and Residential Unitary and Split System HVAC Cooling Equipment-Efficiency Upgrade Evaluation Protocol. Chapter 5: Residential Furnaces and Boilers Evaluation Protocol. Chapter 6: Residential Lighting Evaluation Protocol. Chapter 7: Refrigerator Recycling Evaluation Protocol. Chapter 8: Whole-Building Retrofit with Consumption Data Analysis Evaluation Protocol. Chapter 9: Metering Cross-Cutting Protocol. Chapter 10: Peak Demand and Time-Differentiated Energy Savings Cross-Cutting Protocol. Chapter 11: Sample Design Cross-Cutting Protocol. Chapter 12: Survey Design and Implementation for Estimating Gross Savings Cross-Cutting Protocol. Chapter 13: Assessing Persistence and Other Evaluation Issues Cross-Cutting Protocol. Chapter 14: Chiller Evaluation Protocol. Chapter 15: Commercial New Construction Evaluation Protocol. Chapter 16: Retrocommissioning Evaluation Protocol. Chapter 17: Residential Behavior Evaluation Protocol. Chapter 18: Variable Frequency Drive Evaluation Protocol. Chapter 19: HVAC Controls (DDC/EMS/BAS) Evaluation Protocol. Chapter 20: Data. Center IT Efficiency Measures Evaluation Protocol. Chapter 21: Estimating Net Savings - Common Practices. Chapter 22: Compressed Air Evaluation Protocol. Chapter 23: Combined Heat and Power Evaluation Protocol. Chapter 24: Strategic Energy Management (SEM) Evaluation Protocol

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Stochastic projective splitting

Here, we present a new, stochastic variant of the projective splitting (PS) family of algorithms for inclusion problems involving the sum of any finite number of maximal monotone operators. This new variant uses a stochastic oracle to evaluate one of the operators, which is assumed to be Lipschitz continuous, and (deterministic) resolvents to process the remaining operators. Our proposal is the first version of PS with such stochastic capabilities. We envision the primary application being machine learning (ML) problems, with the method’s stochastic features facilitating “mini-batch” sampling of datasets. Since it uses a monotone operator formulation, the method can handle not only Lipschitz-smooth loss minimization, but also min–max and noncooperative game formulations, with better convergence properties than the gradient descent-ascent methods commonly applied in such settings. The proposed method can handle any number of constraints and nonsmooth regularizers via projection and proximal operators. We prove almost-sure convergence of the iterates to a solution and a convergence rate result for the expected residual, and close with numerical experiments on a distributionally robust sparse logistic regression problem.

97 MATHEMATICS AND COMPUTING↗

ν-point energy correletors with F AST EEC: Small-x physics from LHC jets

In recent years, energy correlators have emerged as a powerful tool for studying jet substructure, with promising applications such as probing the hadronization transition, analyzing the quark-gluon plasma, and improving the precision of top quark mass measurements. The projected N-point correlator measures correlations between N final-state particles by tracking the largest separation between them, showing a scaling behavior related to DGLAP splitting functions. These correlators can be analytically continued in N, commonly referred to as ν-correlators, allowing access to non-integer moments of the splitting functions. Of particular interest is the ν → 0 limit, where the small momentum fraction behavior of the splitting functions requires resummation. Originally, the computational complexity of evaluating ν-correlators for M particles scaled as 2 2M , making it impractical for real-world analyses. However, by using recursion, we reduce this to M 2M , and through the FastEEC method of dynamically resolving subjets, M is replaced by the number of subjets. This breakthrough enables, for the first time, the computation of ν-correlators for LHC data. In practice, limiting the number of subjets to 16 is sufficient to achieve percent-level precision, which we validate using known integer-ν results and convergence tests for non-integer ν. We have implemented this in an update to FastEEC and conducted an initial study of power-law scaling in the perturbative regime as a function of ν, using CMS Open Data on jets. The results agree with DGLAP evolution, except at small ν, where the anomalous dimension saturates to a value that matches the BFKL anomalous dimension.

Energy correlators↗

Developing Methods to Assess Changes in Mechanical Properties of Shale Modified by Engineered Mineral Precipitation

Fractures in subsurface shale formations serve multiple purposes, for example, in the recovery of resources in hydraulic fracturing or as potential harmful leakage passages through caprocks that may contribute undesired fluids to the atmosphere or functional groundwater aquifers. A proposed method to seal or influence fracture properties is Ureolysis-Induced Calcium Carbonate Precipitation (UICP), a bio-mineralization technology driven by the enzymatic hydrolysis of urea, resulting in the formation of calcium carbonate. Sporosarcina pasteurii is a common microbe used as the source of the urease enzyme that catalyzes the chemical reaction. The resulting calcium carbonate can bridge the gaps in fractured shale and reduce fluid flow through fractures. However, there is little information on how this process affects the mechanical properties of the resulting biomineralized shale. This study represents the first step toward determining the influence of UICP treatment on shale material and its subsequent mechanical strength properties. This methods development study aims to determine the effect that temperature has on the tensile strength of intact, unfractured shale cores (2.54 cm (1 in) diameter, 5.08 cm (2 in) long). Tensile strength was determined indirectly using a modified Brazilian test where the splitting tensile strength is attained by applying a compressive load onto the core. Shale cores from Eagle Ford and Wolfcamp formations were tested at both room temperature and 60°C to determine if increased temperature influences the tensile strength of the rock. This data will help to assess the necessity of testing biomineralized cores at temperature. Though 60°C may not mimic subsurface temperatures of the shales used in this study, it was chosen due to limitations of the UICP process while still approaching temperatures of shale formations. This project aims to evaluate what effect temperature has on the mechanical properties of intact shale cores so that engineered or natural rock fractures that are sealed by biomineralization can be better understood.

Bedey, Kayla↗

Unraveling Electronic and Vibrational Coherences Following a Charge Transfer Process in a Photosystem II Reaction Center

A reaction center is a unique biological system that performs the initial charge separation within a Photosystem II (PSII) multiunit enzyme, which eventually drives the catalytic water-splitting in plants and algae. The possible role of quantum coherences coinciding with the energy and charge transfer processes in PSII reaction center is one of the active areas of research. Here, we study these quantum coherences by using a numerically exact method on an excitonic dimer model, including linear vibronic coupling and employing optimal parameters from experimental two-dimensional coherent spectroscopic measurements. This enables us to precisely capture the excitonic interaction between pigments and the dissipation of the energy from electronic and charge-transfer (CT) states to the protein environment. We employ the time nonlocal (TNL) quantum master equation to calculate the population dynamics, which yields numerically reliable results. The calculated results show that, due to the strong dissipation, the lifetime of electronic coherence is too short to have direct participation in the charge transfer processes. However, there are long-lived vibrational coherences present in the system at frequencies close to the excitionic energy gap. These are strongly coupled with the electronic coherences, which makes the detection of the electronic coherences with conventional techniques very challenging. Additionally, we unravel the strong excitonic interaction of radical pair (PD1 and PD2) in the reaction center, which results in a long-lived electronic coherence of >100 fs, even at room temperature. Our work provide important physical insight to the charge separation process in PSII reaction center, which may be helpful for better understanding of photophysical processes in other natural and artificial light-harvesting systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A split finite element algorithm for the compressible Navier-Stokes equations

An accurate and efficient numerical solution algorithm is established for solution of the high Reynolds number limit of the Navier-Stokes equations governing the multidimensional flow of a compressible essentially inviscid fluid. Finite element interpolation theory is used within a dissipative formulation established using Galerkin criteria within the Method of Weighted Residuals. An implicit iterative solution algorithm is developed, employing tensor product bases within a fractional steps integration procedure, that significantly enhances solution economy concurrent with sharply reduced computer hardware demands. The algorithm is evaluated for resolution of steep field gradients and coarse grid accuracy using both linear and quadratic tensor product interpolation bases. Numerical solutions for linear and nonlinear, one, two and three dimensional examples confirm and extend the linearized theoretical analyses, and results are compared to competitive finite difference derived algorithms.

Baker, A. J.↗

Computing with high-resolution upwind schemes for hyperbolic equations

Computational aspects of modern high-resolution upwind finite-difference schemes for hyperbolic systems of conservation laws are examined. An operational unification is demonstrated for constructing a wide class of flux-difference-split and flux-split schemes based on the design principles underlying total variation diminishing (TVD) schemes. Consideration is also given to TVD scheme design by preprocessing, the extension of preprocessing and postprocessing approaches to general control volumes, the removal of expansion shocks and 'glitches', relaxation methods for implicit TVD schemes, and a new family of high-accuracy TVD schemes.

Chakravarthy, S. R.↗

Development of a linearized unsteady aerodynamic analysis for cascade gust response predictions

A method for predicting the unsteady aerodynamic response of a cascade of airfoils to entropic, vortical, and acoustic gust excitations is being developed. Here, the unsteady flow is regarded as a small perturbation of a nonuniform isentropic and irrotational steady background flow. A splitting technique is used to decompose the linearized unsteady velocity into rotational and irrotational parts leading to equations for the complex amplitudes of the linearized unsteady entropy, rotational velocity, and velocity potential that are coupled only sequentially. The entropic and rotational velocity fluctuations are described by transport equations for which closed-form solutions in terms of the mean-flow drift and stream functions can be determined. The potential fluctuation is described by an inhomogeneous convected wave equation in which the source term depends on the rotational velocity field, and is determined using finite-difference procedures. The analytical and numerical techniques used to determine the linearized unsteady flow are outlined. Results are presented to indicate the status of the solution procedure and to demonstrate the impact of blade geometry and mean blade loading on the aerodynamic response of cascades to vortical gust excitations. The analysis described herein leads to very efficient predictions of cascade unsteady aerodynamic response phenomena making it useful for turbomachinery aeroelastic and aeroacoustic design applications.

Verdon, Joseph M.↗

Three-dimensional unstructured grid Euler method applied to turbine blades

Flow through a turbine annular cascade is calculated using a three-dimensional Euler method based on unstructured tetrahedral meshes. The equations are integrated in time using an explicit Runge-Kutta time-stepping scheme. The inviscid flux terms are discretized using a cell-centered finite-volume formulation with upwind flux-difference splitting. The tetrahedral meshes around the turbine blade are generated using an advancing-front technique with forced geometric periodicity between the blades. Good agreement is obtained between the present calculation and the experiment for both surface pressure distribution and flow behavior in the passage between the blades, demonstrating the capability of the present methodology for turbomachinery flow applications.

Kwon, Oh J.↗

Methodology for Sensitivity Analysis, Approximate Analysis, and Design Optimization in CFD for Multidisciplinary Applications

An incremental iterative formulation together with the well-known spatially split approximate-factorization algorithm, is presented for solving the large, sparse systems of linear equations that are associated with aerodynamic sensitivity analysis. This formulation is also known as the 'delta' or 'correction' form. For the smaller two dimensional problems, a direct method can be applied to solve these linear equations in either the standard or the incremental form, in which case the two are equivalent. However, iterative methods are needed for larger two-dimensional and three dimensional applications because direct methods require more computer memory than is currently available. Iterative methods for solving these equations in the standard form are generally unsatisfactory due to an ill-conditioned coefficient matrix; this problem is overcome when these equations are cast in the incremental form. The methodology is successfully implemented and tested using an upwind cell-centered finite-volume formulation applied in two dimensions to the thin-layer Navier-Stokes equations for external flow over an airfoil. In three dimensions this methodology is demonstrated with a marching-solution algorithm for the Euler equations to calculate supersonic flow over the High-Speed Civil Transport configuration (HSCT 24E). The sensitivity derivatives obtained with the incremental iterative method from a marching Euler code are used in a design-improvement study of the HSCT configuration that involves thickness. camber, and planform design variables.

Taylor, Arthur C., III↗

Generation and Radiation of Acoustic Waves from a 2-D Shear Layer

A parallel numerical simulation of the radiation of sound from an acoustic source inside a 2-D jet is presented in this paper. This basic benchmark problem is used as a test case for scattering problems that are presently being solved by using the Impedance Mismatch Method (IMM). In this technique, a solid body in the domain is represented by setting the acoustic impedance of each medium, encountered by a wave, to a different value. This impedance discrepancy results in reflected and scattered waves with appropriate amplitudes. The great advantage of the use of this method is that no modifications to a simple Cartesian grid need to be made for complicated geometry bodies. Thus, high order finite difference schemes may be applied simply to all parts of the domain. In the IMM, the total perturbation field is split into incident and scattered fields. The incident pressure is assumed to be known and the equivalent sources for the scattered field are associated with the presence of the scattering body (through the impedance mismatch) and the propagation of the incident field through a non-uniform flow. An earlier version of the technique could only handle uniform flow in the vicinity of the source and at the outflow boundary. Scattering problems in non-uniform mean flow are of great practical importance (for example, scattering from a high lift device in a non-uniform mean flow or the effects of a fuselage boundary layer). The solution to this benchmark problem, which has an acoustic wave propagating through a non-uniform mean flow, serves as a test case for the extensions of the IMM technique.

Agarwal, Anurag↗

Effects of a Rotating Aerodynamic Probe on the Flow Field of a Compressor Rotor

An investigation of distortions of the rotor exit flow field caused by an aerodynamic probe mounted in the rotor is described in this paper. A rotor total pressure Kiel probe, mounted on the rotor hub and extending up to the mid-span radius of a rotor blade channel, generates a wake that forms additional flow blockage. Three types of high-response aerodynamic probes were used to investigate the distorted flow field behind the rotor. These probes were: a split-fiber thermo-anemometric probe to measure velocity and flow direction, a total pressure probe, and a disk probe for in-flow static pressure measurement. The signals acquired from these high-response probes were reduced using an ensemble averaging method based on a once per rotor revolution signal. The rotor ensemble averages were combined to construct contour plots for each rotor channel of the rotor tested. In order to quantify the rotor probe effects, the contour plots for each individual rotor blade passage were averaged into a single value. The distribution of these average values along the rotor circumference is a measure of changes in the rotor exit flow field due to the presence of a probe in the rotor. These distributions were generated for axial flow velocity and for static pressure.

Lepicovsky, Jan↗