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 937 records · Page 52

Split representation of adaptively compressed polarizability operator

The polarizability operator plays a central role in density functional perturbation theory and other perturbative treatment of first principle electronic structure theories. The cost of computing the polarizability operator generally scales as O(Ne4) where Ne is the number of electrons in the system. The recently developed adaptively compressed polarizability operator (ACP) formulation [L. Lin, Z. Xu and L. Ying, Multiscale Model. Simul. 2017] reduces such complexity to O(Ne 3 ) in the context of phonon calculations with a large basis set for the first time, and demonstrates its effectiveness for model problems. In this paper, we improve the performance of the ACP formulation by splitting the polarizability into a near singular component that is statically compressed, and a smooth component that is adaptively compressed. The new split representation maintains the O(Ne 3 ) complexity, and accelerates nearly all components of the ACP formulation, including Chebyshev interpolation of energy levels, iterative solution of Sternheimer equations, and convergence of the Dyson equations. For simulation of real materials, we discuss how to incorporate nonlocal pseudopotentials and finite temperature effects. In this work, we demonstrate the effectiveness of our method using one-dimensional model problem in insulating and metallic regimes, as well as its accuracy for real molecules and solids.

97 MATHEMATICS AND COMPUTING↗

Compound Specific Hydrogen Isotope Composition of Type II and III Kerogen Extracted by Pyrolysis-GC-MS-IRMS

The use of Hydrogen (H) isotopes in understanding oil and gas resource plays is in its infancy. Described here is a technique for H isotope analysis of organic compounds pyrolyzed from oil and gas shale-derived kerogen. Application of this technique will progress our understanding. This work complements that of Pernia et al. (2013, this meeting) by providing a novel method for the H isotope analysis of specific compounds in the characterization of kerogen extracted by analytically diverse techniques. Hydrogen isotope analyses were carried out entirely "on-line" utilizing a CDS 5000 Pyroprobe connected to a Thermo Trace GC Ultra interfaced with a Thermo MAT 253 IRMS. Also, a split of GC-separated products was sent to a DSQ II quadrupole MS to make semi-quantitative compositional measurements of the extracted compounds. Kerogen samples from five different basins (type II and III) were dehydrated (heated to 80 C overnight in vacuum) and analyzed for their H isotope compositions by Pyrolysis-GC-MS-TC-IRMS. This technique takes pyrolysis products separated via GC and reacts them in a high temperature conversion furnace (1450 C) which quantitatively forms H2, following a modified method of Burgoyne and Hayes, (1998, Anal. Chem., 70, 5136-5141). Samples ranging from approximately 0.5 to 1.0mg in size, were pyrolyzed at 800 C for 30s. Compounds were separated on a Poraplot Q GC column. Hydrogen isotope data from all kerogen samples typically show enrichment in D from low to high molecular weight compounds. Water (H2O) average deltaD = -215.2 (V-SMOW), ranging from -271.8 for the Marcellus Shale to -51.9 for the Polish Shale. Higher molecular weight compounds like toluene (C7H8) have an average deltaD of -89.7 0/00, ranging from -156.0 for the Barnett Shale to -50.0 for the Monterey Shale. We interpret these data as representative of potential H isotope exchange between hydrocarbons and sediment pore water during formation within each basin. Since hydrocarbon H isotopes readily exchange with water, these data may provide some useful information on gas-water or oil-water interaction in resource plays, and further as a possible indicator of paleo-environmental conditions. Alternatively, our data may be an indication of H isotope exchange with water and/or acid during the kerogen isolation process. Either of these interpretations will prove useful when deciphering H isotope data derived from kerogen analysis. More experiments are planned to discern these two or other possible scenarios.

Socki, Richard A.↗

Formulation of an explicit-multiple-time-step time integration method for use in a global primitive equation grid model

With appropriate modifications, a recently proposed explicit-multiple-time-step scheme (EMTSS) is incorporated into the UCLA model. In this scheme, the linearized terms in the governing equations that generate the gravity waves are split into different vertical modes. Each mode is integrated with an optimal time step, and at periodic intervals these modes are recombined. The other terms are integrated with a time step dictated by the CFL condition for low-frequency waves. This large time step requires a special modification of the advective terms in the polar region to maintain stability. Test runs for 72 h show that EMTSS is a stable, efficient and accurate scheme.

Chao, W. C.↗

Systems, method and computer-accessible medium for providing balanced asymmetric interferometry for vibrationally isolated optical scanning probe(s)

An exemplary apparatus can provide radiation to a sample(s), which can include, for example, a radiation source arrangement configured to provide radiation, a beam splitter configured to split the radiation into (i) a first radiation, and (ii) a second radiation. An optical element can also be provided which, in operation, can, e.g., (a) receive the first radiation and the second radiation, (b) reflect the first radiation as a reference radiation, (c) provide the second radiation as illumination for the sample(s), (d) receive a resultant radiation from the sample(s) that can be based on the illumination from the second radiation, and (e) provide the reference radiation and the resultant radiation to be detected and used for interferometric imaging or spectroscopy.

McLeod, Alexander Swinton↗

An efficient method for solving the steady Euler equations

An efficient numerical procedure for solving a set of nonlinear partial differential equations is given, specifically for the steady Euler equations. Solutions of the equations were obtained by Newton's linearization procedure, commonly used to solve the roots of nonlinear algebraic equations. In application of the same procedure for solving a set of differential equations we give a theorem showing that a quadratic convergence rate can be achieved. While the domain of quadratic convergence depends on the problems studied and is unknown a priori, we show that firstand second-order derivatives of flux vectors determine whether the condition for quadratic convergence is satisfied. The first derivatives enter as an implicit operator for yielding new iterates and the second derivatives indicates smoothness of the flows considered. Consequently flows involving shocks are expected to require larger number of iterations. First-order upwind discretization in conjunction with the Steger-Warming flux-vector splitting is employed on the implicit operator and a diagonal dominant matrix results. However the explicit operator is represented by first- and seond-order upwind differencings, using both Steger-Warming's and van Leer's splittings. We discuss treatment of boundary conditions and solution procedures for solving the resulting block matrix system. With a set of test problems for one- and two-dimensional flows, we show detailed study as to the efficiency, accuracy, and convergence of the present method.

Liou, M. S.↗

Structure and tunable temperature coefficient of magnetization of Mn 4-x Ga x C alloys prepared by induction melting method

The magnetization of most magnetic materials decreases monotonically with increasing temperature. In this work, we found that the temperature coefficient of magnetization of Mn 4-x Ga x C alloys can be tuned from negative values to positive values by controlling the composition x. The antiperovskite type Mn 4-x Ga x C (0.05 ≤ x ≤ 0.75) alloys were prepared by using induction melting method, which is more efficient in large-scale production and obtaining full-density alloys in comparison with the traditional solid-state-reaction method. The values of the temperature coefficient of magnetization of Mn 4-x Ga x C change continuously from negative to positive with decreasing x. The Mn 4-x Ga x C alloys with highly thermal-stable magnetization is expected to present in the composition range of 0.15 < x < 0.25. The saturation magnetization of Mn 4-x Ga x C increases with increasing x, owing to the reduced number of antiferromagnetically coupled Mn atoms at the cubic corner with the face centered Mn atoms. Most Mn 4-x Ga x C alloys with varying x display near-zero remanent magnetization and coercivity at room temperature. The Currie temperature of Mn 4-x Ga x C decreases with increasing x. The x-ray photoelectron spectra of Mn 2p, Ga 2p, and C 1 s reveal distinct splitting due to the diverse chemical states of these atoms at different lattice positions and/or phases. Our work has developed a class of alloys capable of offering a desired temperature coefficient of magnetization across a broad temperature range, thereby offering a method to manipulate the thermodynamics of magnetization.

Magnetic properties↗

Finite element analyses of single particle crushing tests incorporating computed tomography imaging and damage mechanics

Fracture of individual grains within a granular assembly affect the strength and deformation behavior of granular materials under large stress states. Interparticle forces, leading to contact stresses and ultimately fracture initiation, are influenced by particle morphology. A numerical method that incorporates both particle fracture and morphology can provide a more accurate model of a granular material’s macro-scale response. This research investigates an approach to modeling fracture initiation and propagation within individual grains, utilizing high resolution X-ray computed tomography (CT) imaging to consider grain morphology and an explicit finite element code optimized for high performance computing which incorporates damage mechanics. Using precisely measured force-displacement curves from single particle crushing tests of manufactured quartz spheres, the method is validated, and the effects of fracture properties on grain fragmentation are determined. This approach is then used to simulate single particle crushing of Ottawa sand, with grain morphology incorporated through discretized CT images. Here the effect of grain orientation on the force and applied displacement at catastrophic splitting is also explored. The proposed discrete particle based finite element framework can be applied to modeling granular assemblies, considering the effect of individual particle shape and fracture on the assembly’s deformation response through well calibrated numerical simulations.

42 ENGINEERING↗

OmicsMLMentor: A Web Application for Guided Machine Learning Analysis of Omics Data

Expression-based omics technologies (e.g. proteomics, metabolomics, transcriptomics, etc.) increasingly rely on supervised and unsupervised machine learning (ML) models to find key biomolecules distinguishing conditions, identify natural groupings in biological data, or generate predictions for outcomes of interest. Fitting ML models to omics data presents several challenges, including handling missing data, selecting a normalization method, choosing a valid model, and optimizing hyperparameters, all requiring statistical programming skills to address these challenges. Thus, the open-source web application SLOPE was designed to lower the barrier to ML modeling for omics data. SLOPE supports the fitting of 15 ML models (10 supervised and 5 unsupervised) tailored to omics datasets, such as proteomics, metabolomics, lipidomics, and transcriptomics. SLOPE offers several omics-specific features, including methods for handling missingness (imputation, conversion, removal), normalization tests, ranking of models based on the structure of a user’s data and user input, and optimal hyperparameter selections using cross-validation splits. By streamlining ML workflows for omics analysis, SLOPE address critical gaps in existing online web tools, facilitating a broader adoption of these models for omics research. Here, SLOPE is applied to data from a lignin exposure study to highlight the workflow for fitting both supervised and unsupervised models to data.

lipidomics↗

Zero-field spin resonance in graphene with proximity-induced spin-orbit coupling

In this work we investigate collective spin excitations in graphene with proximity-induced spin-orbit coupling (SOC) of the Rashba and valley-Zeeman types, as it is the case, e.g., for graphene on transition-metal-dichalcogenide substrates. It is shown that, even in the absence of an external magnetic field, such a system supports collective spin modes, which correspond to coupled oscillations of the uniform and valley-staggered magnetizations. These modes can be detected via both zero-field electron spin resonance (ESR) and zero-field electric-dipole spin resonance (EDSR), with EDSR response arising solely due to Rashba SOC. We analyze the effect of electron-electron interaction within the Fermi-liquid kinetic equation and show that the interaction splits both the ESR and EDSR peaks into two. The magnitude of splitting and the relative weights of the resonances can be used to extract the spin-orbit coupling constants and many-body interaction parameters that may not be accessible by other methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Development of secondary flow and vorticity in curved ducts, cascades, and rotors, including effects of viscosity and rotation

This paper is concerned with the numerical solution of the secondary vorticity equations in curved ducts, cascades, and rotors. The classical approach of splitting the flow into primary and secondary flow fields is employed and extended to include effects of viscosity, rotation, and density stratification. All elliptical effects are neglected and the Crank-Nicolson method is used to solve the secondary vorticity equation. The secondary flow field is obtained by solving a Poisson equation using a successive over-relaxation method. The results are compared with theoretical and experimental data from stationary ducts, compressor and turbine cascades. The agreement is good for most of the cases.

Pouagare, M.↗

Recent developments in FEM-CFD

The current status of CFD with regard to unstructured grids employing finite element methods and Eulerian frames is reviewed. Algorithms suitable for the computation of large three-dimensional problems involving flow past arbitrary geometries are developed. Adaptive mesh refinement strategy is reviewed, and domain splitting or local time-stepping are briefly addressed. The development of search algorithms of optimal order, variable time-stepping Jacobi smoothers for elliptic problems, and transport concepts for hyperbolics to help achieve good performance for unstructured multigrid processes is discussed. As examples, transient supersonic flow in a channel, regular shock reflection of a wall, viscous flow past a protruberance, potential flow past a cylinder, and Burgers equation are considered.

Loehner, R.↗

Conditionally-averaged structures in wall-bounded turbulent flows

The quadrant-splitting and the wall-shear detection techniques were used to obtain ensemble-averaged wall layer structures. The two techniques give similar results for Q4 events, but the wall-shear method leads to smearing of Q2 events. Events were found to maintain their identity for very long times. The ensemble-averaged structures scale with outer variables. Turbulence producing events were associated with one dominant vortical structure rather than a pair of counter-rotating structures. An asymmetry-preserving averaging scheme was devised that allowed a picture of the average structure which more closely resembles the instantaneous one, to be obtained.

Guezennec, Yann G.↗

Turbomachine Sealing and Secondary Flows: Review of Rotordynamics Issues in Inherently Unsteady Flow Systems With Small Clearances - Part 2

Today's computational methods enable the determination of forces in complex systems, but without field validation data, or feedback, there is a high risk of failure when the design envelope is challenged. The data of Childs and Bently and field data reported in NASA Conference Proceedings serve as sources of design information for the development of these computational codes. Over time all turbomachines degrade and instabilities often develop, requiring responsible, accurate, turbomachine diagnostics with proper decisions to prevent failures. Tam et al. (numerical) and Bently and Muszynska (analytical) models corroborate and implicate that destabilizing factors are related through increases in the fluid-force average circumferential velocity. The stability threshold can be controlled by external swirl and swirl brakes and increases in radial fluid film stiffness (e.g., hydrostatic and ambient pressures) to enhance rotor stability. Also cited are drum rotor self-excited oscillations, where the classic fix is to add a split or severed damper ring or cylindrical damper drum, and the Benkert-Wachter work that engendered swirl brake concepts. For a smooth-operating, reliable, long-lived machine, designers must pay very close attention to sealing dynamics and diagnostic methods. Correcting the seals enabled the space shuttle main engine high-pressure fuel turbopump (SSME HPFTP) to operate successfully.

Hendricks, R. C.↗

A method of obtaining signal components of residual carrier signal with their power content and computer simulation

A novel algorithm to obtain all signal components of a residual carrier signal with any number of channels is presented. The phase modulation type may be NRZ-L or split phase (Manchester). The algorithm also provides a simple way to obtain the power contents of the signal components. Steps to recognize the signal components that influence the carrier tracking loop and the data tracking loop at the receiver are given. A computer program for numerical computation is also provided.

Kantak, Anil V.↗

GC/MS Method Development for Separating Lunar Volatile Ice Simulant Headspace Gases

Various investigators propose the lunar surface contains widely distributed volatiles, especially water- like species, i.e. OH and H2O. Surface volatiles are theorized to exist as a hydrated regolith layer, concentrated in extremely cold polar permanently shadowed regions (PSR), and/or solar wind implantation reservoirs in lunar glasses. The proposed sources of lunar surface volatiles range from cometary impacts, solar wind, or a supply present during moon formation. Future Artemis missions aim to collect and return the samples containing volatiles collected near lunar polar craters or PSRs. We, as advanced curation scientists, are responsible for developing techniques and methodologies for preserving returned sample integrity as much as possible. Pristine volatile-bearing samples are invaluable to the scientific community seeking to unravel the history of the solar system. Realistically, a sample will experience alteration during collection, transportation back to earth, and storage. The Planetary Exploration and Astromaterials Research Lab (PEARL) seeks to understand temperature and pressure effects on high-fidelity volatile-containing regolith simulants, the foundation for the future of cold curation. This abstract outlines the separation, identification, and quantification of headspace gases over volatile ice feed stock material using gas chromatography/mass spectrometry (GC/MS). Preliminary objectives concentrated on sample handling, reproducibility, and understanding the elution characteristics for each analyte. Initial GC/MS method development experiments utilized diluted static headspace sample preparation. Diluted samples were used because sampling headspace gases directly from a vial containing liquid analyte resulted in overloading of the column and detector. Overloading is evident based on chromatogram peak shapes and instrument contamination, or carry over, between experiments. A mixture of three alcohols were used for a majority of the sample handling and reproducibility studies. Reproducibility was tested via multiple users, calibration curves, and check standards. Stock solutions of condensed lunar volatile analytes included methanol, ammonia in methanol, hydrogen sulfide in water, and an equal volume mixture of methanol, ethanol, and isopropanol. Current samples use room air as the headspace sample matrix, however future experiments will incorporate an inert purge gas, such as argon or nitrogen. Three mL of each analyte solution were capped in separate 20 mL crimp top GC vials. Dilutions were carried out by removing an aliquot of headspace gases with a calibrated 1 mL gastight syringe and immediately transferring to a 20 mL capped crimp top vial. The GC/MS is a Thermo Fisher Trace 1310/ISQ 7000 with a TriPlus RSH autosampler and split/splitless injector module. The experiments outlined in this abstract use the following hardware: a 2.5 mL gastight headspace syringe tool, 1 mm ID x 78.5 mm length ultra-inert straight injection liner, and a TG-BondQ 30 m × 0.32 mm × 10 μm column. Various parameters, such as hardware selection and the temperature, pressure, and split ratio set points, continue to evolve as the overall experiment is refined. Diluted headspace chromatograms were collected for the individual stock solutions. Retention times, peak shapes, and mass spectra were evaluated and added to the data processing method for each molecule of interest. Figure 1 shows the total ion chromatograms for the three major lunar volatile simulant stock solutions: methanol, 7 N ammonia in methanol, and 0.4% hydrogen sulfide in water. Tailing peak shapes for ammonia (2.98 min rt) and water (4.06 min rt) indicate the molecules are not properly eluting from the selected column with the current separation method. Additionally, hydrogen sulfide and ammonia have overlapping peak windows, which could impact quantification. Ongoing experiments aim to address the peak shape and overlapping via the separation method and hardware selection. Sample preparation reproducibility experiments used stock solution containing equal volumes of a non- interactive mixture of methanol, ethanol, and isopropanol. Mass spectrum ion traces were used to identify and quantify all three alcohols. Peaks were automatically detected, identified, and integrated through the mass spectra detection and processing parameters. Calibration response curves and check standards were used to evaluate the validity of the sample preparation procedure. Figure 2 shows the methanol chromatogram peak area versus total headspace dilution volume transferred from the alcohol mixture vial. The calibration response curves and check standards validate sample preparation procedure. Continuing data analysis efforts are working towards correlating the peak area and instrument response factor to the headspace analyte concentration and condensed phase composition. Static headspace gas chromatography theory relies on Dalton’s law, Raoult’s law, Henry’s Law, and the Kolb and Ettre equation to associate peak area to the analyte composition in a non-ideal solution. Equation 1 is a simplified expression derived from the aforementioned theories. Future experiments involve liquid injections of the individual stock solutions, liquid and headspace analysis of various stock solution combinations, and the addition of regolith simulants to the mixtures. Temperature is another variable expected to affect reaction rates and will be explored.

Cecilia L. Amick↗

A method for the measurement and the statistical analysis of atmospheric turbulence

The instantaneous values of output voltages representing the wind velocity vector and the temperature at different elevations of the 250-foot meteorological tower located at NASA Wallops Flight Center are provided with the three dimensional split-film TSI Model 1080 anemometer system. The output voltages are sampled at a rate of one every 5 milliseconds, digitized and stored on digital magnetic tapes for a time period of approximately 40 minutes, with the use of a specially designed data acqusition system. A new calibration procedure permits the conversion of the digital voltages to the respective values of the temperature and the velocity components in a Cartesian coordinate system connected with the TSI probe with considerable accuracy. Power, cross, coincidence and quadrature spectra of the wind components and the temperature are obtained with the use of the fast Fourier transform. The cosine taper data window and ensemble and frequency smoothing techniques are used to provide smooth estimates of the spectral functions.

Tieleman, H. W.↗

Tension fracture of laminates for transport fuselage. Part 1: Material screening

Transport fuselage structures are designed to contain pressure following a large penetrating damage event. Applications of composites to fuselage structures require a database and supporting analysis on tension damage tolerance. Tests with 430 fracture specimens were used to accomplish the following: (1) identify critical material and laminate variables affecting notch sensitivity; (2) evaluate composite failure criteria; and (3) recommend a screening test method. Variables studied included fiber type, matrix toughness, lamination manufacturing process, and intraply hybridization. The laminates found to have the lowest notch sensitivity were manufactured using automated tow placement. This suggests a possible relationship between the stress distribution and repeatable levels of material inhomogeneity that are larger than found in traditional tape laminates. Laminates with the highest notch sensitivity consisted of toughened matrix materials that were resistant to a splitting phenomena that reduces stress concentrations in major load bearing plies. Parameters for conventional fracture criteria were found to increase with crack length for the smallest notch sizes studied. Most material and laminate combinations followed less than a square root singularity for the largest crack sizes studied. Specimen geometry, notch type, and notch size were evaluated in developing a screening test procedure. Traitional methods of correcting for specimen finite width were found to be lacking. Results indicate that a range of notch sizes must be tested to determine notch sensitivity. Data for a single small notch size (0.25 in. diameter) was found to give no indication of the sensitivity of a particular material and laminate layup to larger notch sizes.

Walker, T. H.↗

Ride-hailing and taxi versus walking: Long term forecasts and implications from large-scale behavioral data

Introduction: Although ride-hailing and taxi trips can potentially reduce single-occupant vehicle trips and auto ownership, they can also replace pedestrian trips. Because physical activity is associated with improved health outcomes, the extent to which ride-hailing and taxi travel captures walking's mode share is of interest to policymakers. Methods: Based on large-scale behavioral data from the 2017 U.S. National Household Travel Survey, this paper reports on the development of a full Bayesian logistic regression model for determining the mode split between (1) ride-hailing and taxi and (2) walk while accounting for unobserved heterogeneity. The results from the stand-alone model inform two longer-term travel forecasting scenarios: a) higher risk of walk trips converting to ride-hailing and taxi, specifically in the future with high prevalence of automated vehicles, b) higher probability of such trips remaining as walking. Results: The results revealed that some of the important characteristics that increase the likelihood of a traveler using the ride-hailing and taxi mode versus walking include having a longer trip, using a smartphone to access the internet, having an interest in technologies, having a medical condition, and living in a metropolitan area with rail access. Further, the results from the first scenario suggest that an overall increase of up to 2.9% in the ride-hailing and taxi mode share may be expected. The second scenario shows that between 68% and 76% of ride-hailing and taxi trips could be diverted to walking if supportive pedestrian infrastructure were provided in the case study locations. The planning process can be adapted to consider not only congestion, crash, and emissions impacts of such shifts but also the effects of a loss of physical activity. Conclusions: The study findings show how the ride-hailing and taxi mode competes with walking. Further, the findings enable planners to update their regional travel forecasting models; policy makers can thus encourage active travel by prioritizing pedestrian infrastructure investments that may divert ride-hailing and taxi trips to walking. However, equity should be a key consideration to ensure that addressing the competition between these two modal choices does not hinder the provision of pedestrian facilities in communities that depend on walking.

99 GENERAL AND MISCELLANEOUS↗