Search NASA⌕ Search

SEARCH · Search NASA

Results for “mixed 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 865 records · Page 48

A new family of stable elements for the Stokes problem based on a mixed Galerkin/least-squares finite element formulation

Adding to the classical Hellinger-Reissner formulation, a residual form of the equilibrium equation, a new Galerkin/least-squares finite element method is derived. It fits within the framework of a mixed finite element method and is stable for rather general combinations of stress and velocity interpolations, including equal-order discontinuous stress and continuous velocity interpolations which are unstable within the Galerkin approach. Error estimates are presented based on a generalization of the Babuska-Brezzi theory. Numerical results (not presented herein) have confirmed these estimates as well as the good accuracy and stability of the method.

Franca, Leopoldo P.↗

Detection and evaluation of mixed pixels in Landsat agricultural scenes

A major problem area encountered in the identification and estimation of agricultural crop proportions in Landsat imagery involves the large proportion of the pixels which are mixed pixels, whose spectral response is influenced by more than one ground cover type. The development of methods for the detection and estimation of crop proportions in mixed pixels is presently reported. The procedure designated CASCADE, based on the estimation of the gradient image for the detection of mixed pixels, considers the consequences of a linear mixing model and is found to provide a method for the allocation of mixed pixels to the surrounding homogeneous region.

Merickel, M. B.↗

Supersonic free turbulent mixing layers

A mean turbulent closure method with a length scale function is developed for predicting the reduced spreading rate in high Mach number, high Reynolds number free shear layers. The new feature of the method is the inclusion of effects of the pressure-dilation term which appears in the compressible form of the turbulence energy equation. The model indicates that the pressure-dilation correlation may act as a turbulence energy sink or source depending on whether the transverse gradients of mean Mach number and mean density have the same or opposite signs.

Oh, Y. H.↗

Calculation of Compressible Flows past Aerodynamic Shapes by Use of the Streamline Curvature

A simple approximate method is given for the calculation of isentropic irrotational flows past symmetrical airfoils, including mixed subsonic-supersonic flows. The method is based on the choice of suitable values for the streamline curvature in the flow field and the subsequent integration of the equations of motion. The method yields limiting solutions for potential flow. The effect of circulation is considered. A comparison of derived velocity distributions with existing results that are based on calculation to the third order in the thickness ratio indicated satisfactory agreement. The results are also presented in the form of a set of compressibility correction rules that lie between the Prandtl-Glauert rule and the von Karman-Tsien rule (approximately). The different rules correspond to different values of the local shape parameter square root sign YC sub a, in which Y is the ordinate and C sub a is the curvature at a point on an airfoil. Bodies of revolution, completely supersonic flows, and the significance of the limiting solutions for potential flow are also briefly discussed.

Perl, W↗

Scale–Resolving Simulations of a Fundamental Trailing–Edge Cooling Slot Using a Discontinuous–Galerkin Spectral–Element Method

The accurate prediction of turbulent mixing in high-pressure turbines that incorporate various airfoil surface-cooling strategies is becoming increasing critical to the design of modern gas turbine engines where the quest for improved efficiency is driving compressor overall pressure ratios and turbine inlet temperatures to much higher levels than ever before. In the present paper, a recently developed computational capability for accurate and efficient scale–resolving simulations of turbomachinery is extended to study the turbulent mixing mechanism of a simplified abstraction of an airfoil trailing-edge cooling slot - a plane wall jet with finite lip thickness discharging into an ambient flow. The computational capability is based on an entropy– stable, discontinuous–Galerkin approach that extends to arbitrarily high orders of spatial and temporal accuracy. The numerical results show that the present simulations capture the trends observed in the experiments. Discrepancies between the simulations and experiments are believed to be due to differences in the inflow profiles and tunnel side–wall effects. The thick lip configuration leads to a thicker wake and higher unsteadiness in the wall jet compared to the thin lip. A detailed comparison of the turbulent flowfields is presented to highlight differences arising due to lip thickness variations.

Scale Resolving↗

Quantifying uncertainty in stable isotope mixing models

Mixing models are powerful tools for identifying biogeochemical sources and determining mixing fractions in a sample. However, identification of actual source contributors is often not simple, and source compositions typically vary or even overlap, significantly increasing model uncertainty in calculated mixing fractions. This study compares three probabilistic methods, SIAR [ Parnell et al ., 2010] a pure Monte Carlo technique (PMC), and Stable Isotope Reference Source (SIRS) mixing model, a new technique that estimates mixing in systems with more than three sources and/or uncertain source compositions. In this paper, we use nitrate stable isotope examples (δ 15 N and δ 18 O) but all methods tested are applicable to other tracers. In Phase I of a three-phase blind test, we compared methods for a set of six-source nitrate problems. PMC was unable to find solutions for two of the target water samples. The Bayesian method, SIAR, experienced anchoring problems, and SIRS calculated mixing fractions that most closely approximated the known mixing fractions. For that reason, SIRS was the only approach used in the next phase of testing. In Phase II, the problem was broadened where any subset of the six sources could be a possible solution to the mixing problem. Results showed a high rate of Type I errors where solutions included sources that were not contributing to the sample. In Phase III some sources were eliminated based on assumed site knowledge and assumed nitrate concentrations, substantially reduced mixing fraction uncertainties and lowered the Type I error rate. These results demonstrate that valuable insights into stable isotope mixing problems result from probabilistic mixing model approaches like SIRS. The results also emphasize the importance of identifying a minimal set of potential sources and quantifying uncertainties in source isotopic composition as well as demonstrating the value of additional information in reducing the uncertainty in calculated mixing fractions.

54 ENVIRONMENTAL SCIENCES↗

Decontamination of Mixed Paper and Plastic Municipal Solid Waste Increases Low and High Temperature Conversion Yields

With the implementation of China’s Green Fence Policy and the following National Sword Policy, there is a need to divert previously accepted waste materials away from the landfill. Mixed plastic and paper wastes that are too contaminated to be economically recycled can be used for conversion into fuels and products. The effect of common contaminants present in waste streams must be determined to make municipal solid waste (MSW) a viable alternative to agricultural feedstocks for conversion. In this study, MSW was sourced from a dual stream materials recovery facility (MRF) in Emmet County, Michigan and characterized via mass balance. The most common contaminants in mixed paper and plastic were quantified and targeted for decontamination when mixed paper was explored for low temperature conversion pathways and plastics for high temperature conversion pathways. Ink, stickies, and plastic contamination were sorted out of the mixed paper stream and did not have an impact on sugar yields during enzymatic hydrolysis (EH). Although no improvement in yield was found, it was noted that copy paper within the mixed paper stream had an elevated pH and higher concentration of fermentation inhibitors, suggesting that intrinsic properties of some paper types may not be ideal for EH and fermentation conversion. Dilute alkaline pretreatment and mechanical refining of mixed paper increased the EH yield of mixed paper, suggesting this previously established process could be used as a decontamination method. “Dirt” and particulates were washed from mixed plastic using either an aqueous detergent or a non-aqueous chemical washing method. Liquid oil yields were increased during microwave assisted fast pyrolysis. All decontamination techniques met the cost target of less than $30/dry ton with plastic decontamination ranging from $18.16 to $24.81/dry ton. Because dilute alkaline pre-treatment and mechanical refining is considered part of the conversion process, the decontamination cost during feedstock preprocessing was considered $0/dry ton.

09 BIOMASS FUELS↗

In situ analysis of free radicals from the photodecomposition of hydrogen peroxide using a frequency-mixing magnetic detector

We present an analytical method for the real-time detection of free radicals from the photodecomposition (ultraviolet radiation, {lambda} = 254 nm) of hydrogen peroxide (H{sub 2}O{sub 2}) using frequency-mixing magnetic detection. We monitored the free radicals in situ without catalysts or probe molecules. Both water and H{sub 2}O{sub 2} produced frequency-mixing signals under UV radiation, but the water signal was much weaker. The root mean square amplitude of the frequency-mixing signal was found to depend on the initial H{sub 2}O{sub 2} concentration. Considering the physical properties of the reactants, the frequency-mixing signal is attributed to the generation of paramagnetic free radicals by the photodecomposition of H{sub 2}O{sub 2}.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Determining the slag fraction, water/binder ratio and degree of hydration in hardened cement pastes

A method for determining the original mix composition of hardened slag-blended cement-based materials based on analysis of backscattered electron images combined with loss on ignition measurements is presented. The method does not require comparison to reference standards or prior knowledge of the composition of the binders used. Therefore, it is well-suited for application to real structures. The method is also able to calculate the degrees of reaction of slag and cement. Results obtained from an experimental study involving sixty samples with a wide range of water/binder (w/b) ratios (0.30 to 0.50), slag/binder ratios (0 to 0.6) and curing ages (3 days to 1 year) show that the method is very promising. The mean absolute errors for the estimated slag, water and cement contents (kg/m{sup 3}), w/b and s/b ratios were 9.1%, 1.5%, 2.5%, 4.7% and 8.7%, respectively. 91% of the estimated w/b ratios were within 0.036 of the actual values. -- Highlights: •A new method for estimating w/b ratio and slag content in cement pastes is proposed. •The method is also able to calculate the degrees of reaction of slag and cement. •Reference standards or prior knowledge of the binder composition are not required. •The method was tested on samples with varying w/b ratios and slag content.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Optimal Network Topology for Node-Breaker Representations With AC Power Flow Constraints

It has been demonstrated that network topology optimization (NTO) may change the topology of power system networks, and consequently, provide additional flexibility to reduce network congestion and violations. Most NTO problems are formulated based on the bus-branch model in which it is challenging to represent a realistic picture of all substation configurations. In this paper, we explore advantages of substation reconfiguration modeling based on node-breaker representations for NTO problem with full nonlinear alternating current power flow. It also proposes a tailored solution algorithm to solve this nonconvex mixed-integer nonlinear programming through the outer approximation method. The proposed solution approach iterates between a mixed-integer linear programming and a nonlinear subproblem. Additional enhancements to further accelerate the iteration process are illustrated. Numerical case studies demonstrate the relative economic and operational impact of optimal network topology with node-breaker representations.

42 ENGINEERING↗

Higher-order compositional modeling of three-phase flow in 3D fractured porous media based on cross-flow equilibrium

Numerical simulation of multiphase compositional flow in fractured porous media, when all the species can transfer between the phases, is a real challenge. Despite the broad applications in hydrocarbon reservoir engineering and hydrology, a compositional numerical simulator for three-phase flow in fractured media has not appeared in the literature, to the best of our knowledge. In this work, we present a three-phase fully compositional simulator for fractured media, based on higher-order finite element methods. To achieve computational efficiency, we invoke the cross-flow equilibrium (CFE) concept between discrete fractures and a small neighborhood in the matrix blocks. We adopt the mixed hybrid finite element (MHFE) method to approximate convective Darcy fluxes and the pressure equation. This approach is the most natural choice for flow in fractured media. The mass balance equations are discretized by the discontinuous Galerkin (DG) method, which is perhaps the most efficient approach to capture physical discontinuities in phase properties at the matrix-fracture interfaces and at phase boundaries. In this work, we account for gravity and Fickian diffusion. The modeling of capillary effects is discussed in a separate paper. We present the mathematical framework, using the implicit-pressure-explicit-composition (IMPEC) scheme, which facilitates rigorous thermodynamic stability analyses and the computation of phase behavior effects to account for transfer of species between the phases. A deceptively simple CFL condition is implemented to improve numerical stability and accuracy. We provide six numerical examples at both small and larger scales and in two and three dimensions, to demonstrate powerful features of the formulation.

97 MATHEMATICS AND COMPUTING↗

Numerical Diffusion (Mixing) of Material in Numerical Simulations of Hydrodynamics [Slides]

It is often assumed that a material interface between two materials is spread over a few numerical cells in numerical simulations for hydrodynamics. Also, we have the impression that higher order methods introduce less numerical diffusion (mixing) of material. As we know one of the purposes of adaptive mesh refinement (AMR) is to resolve interfaces between materials, but we would like to know how effective AMR is to reduce numerical diffusion of material. We will present our investigation about numerical diffusion (mixing) of material in xRage. The result of the investigation indicates that the assumptions mentioned above are not always valid. In this talk, we will also demonstrate the effectiveness of numerical techniques to reduce numerical diffusion of material, including contact discontinuity steepening, isotropic interface steepening, max interface steepening, material interface reconstruct.

36 MATERIALS SCIENCE↗

Controlling radioisotope proportions when randomly sampling from Dirichlet distributions in PyRIID

As machine learning models for radioisotope quantification become more powerful, likewise the need for high-quality synthetic training data grows as well. For problem spaces that involve estimating the relative isotopic proportions of various sources in gamma spectra it is necessary to generate training data that accurately represents the variance of proportions encountered. In this report, we aim to provide guidance on how to target a desired variance of proportions which are randomly when using the PyRIID Seed Mixer, which samples from a Dirichlet distribution. We provide a method for properly parameterizing the Dirichlet distribution in order to maintain a constant variance across an arbitrary number of dimensions, where each dimension represents a distinct source template being mixed. We demonstrate that our method successfully parameterizes the Dirichlet distribution to target a specific variance of proportions, provided that several conditions are met. This allows us to follow a principled technique for controlling how random mixture proportions are generated which are then used downstream in the synthesis process to produce the final, noisy gamma spectra.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Divergence preserving discrete surface integral methods for Maxwell's curl equations using non-orthogonal unstructured grids

Several new discrete surface integral (DSI) methods for solving Maxwell's equations in the time-domain are presented. These methods, which allow the use of general nonorthogonal mixed-polyhedral unstructured grids, are direct generalizations of the canonical staggered-grid finite difference method. These methods are conservative in that they locally preserve divergence or charge. Employing mixed polyhedral cells, (hexahedral, tetrahedral, etc.) these methods allow more accurate modeling of non-rectangular structures and objects because the traditional stair-stepped boundary approximations associated with the orthogonal grid based finite difference methods can be avoided. Numerical results demonstrating the accuracy of these new methods are presented.

Madsen, Niel K.↗

Origin of the $ν$ = 1/2 fractional quantum Hall effect in wide quantum wells

The nature of the fractional quantum Hall effect at ν = 1/2, observed in wide quantum wells al- most three decades ago, is still under debate. Previous studies have investigated it by the variational Monte Carlo method, which assumes that the transverse wave function and the gap between the symmetric and antisymmetric subbands obtained in a local density approximation at zero magnetic field remain valid even at high perpendicular magnetic fields; this method also ignores the effect of Landau level mixing. We develop in this work a three-dimensional fixed phase diffusion Monte Carlo method, which gives, in a single framework, the total energies of various candidate states in a finite width quantum well, including Landau level mixing, directly in a large magnetic field. This method can be applied to one-component states, and also to two-component states in the limit where the symmetric and antisymmetric bands are nearly degenerate. Our three-dimensional fixed-phase diffusion Monte Carlo calculations find that the one-component composite-fermion Fermi sea and the one-component Pfaffian states are very close in energy for a range of quantum well widths and densities, suggesting that the observed 1/2 fractional quantum Hall state in wide quantum wells is likely to be the one-component Pfaffian state. Furthermore, we hope that this will motivate further experimental studies of this state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Scalar conservation and boundedness in simulations of compressible flow

With the proper combination of high-order, low-dissipation numerical methods, physics-based subgrid-scale models, and boundary conditions it is becoming possible to simulate many combustion flows at relevant conditions. However, non-premixed flows are a particular challenge because the thickness of the fuel/oxidizer interface scales inversely with Reynolds number. Sharp interfaces can also be present in the initial or boundary conditions. When higher-order numerical methods are used, there are often aphysical undershoots and overshoots in the scalar variables (e.g.passive scalars, species mass fractions or progress variable). These numerical issues are especially prominent when low-dissipation methods are used, since sharp jumps in flow variables are not always coincident with regions of strong variation in the scalar fields: consequently, special detection mechanisms and dissipative fluxes are needed. Most numerical methods diffuse the interface, resulting in artificial mixing and spurious reactions. In this paper, we propose a numerical method that mitigates this issue. As a result, we present methods for passive and active scalars, and demonstrate their effectiveness with several examples.

97 MATHEMATICS AND COMPUTING↗

Uncertainty Visualization of Critical Points of 2D Scalar Fields for Parametric and Nonparametric Probabilistic Models

This paper presents a novel end-to-end framework for closed-form computation and visualization of critical point uncertainty in 2D uncertain scalar fields. Critical points are fundamental topological descriptors used in the visualization and analysis of scalar fields. The uncertainty inherent in data (e.g., observational and experimental data, approximations in simulations, and compression), however, creates uncertainty regarding critical point positions. Uncertainty in critical point positions, therefore, cannot be ignored, given their impact on downstream data analysis tasks. Here, in this work, we study uncertainty in critical points as a function of uncertainty in data modeled with probability distributions. Although Monte Carlo (MC) sampling techniques have been used in prior studies to quantify critical point uncertainty, they are often expensive and are infrequently used in production-quality visualization software. We, therefore, propose a new end-to-end framework to address these challenges that comprises a threefold contribution. First, we derive the critical point uncertainty in closed form, which is more accurate and efficient than the conventional MC sampling methods. Specifically, we provide the closed-form and semianalytical (a mix of closed-form and MC methods) solutions for parametric (e.g., uniform, Epanechnikov) and nonparametric models (e.g., histograms) with finite support. Second, we accelerate critical point probability computations using a parallel implementation with the VTK-m library, which is platform portable. Finally, we demonstrate the integration of our implementation with the ParaView software system to demonstrate near-real-time results for real datasets.

97 MATHEMATICS AND COMPUTING↗

Vortex-scalar element calculations of a diffusion flame stabilized on a plane mixing layer

The vortex-scalar element method, a scheme which utilizes vortex elements to discretize the region of high vorticity and scalar elements to represent species or temperature fields, is utilized in the numerical simulations of a two-dimensional reacting mixing layer. Computations are performed for a diffusion flame at high Reynolds and Peclet numbers without resorting to turbulence models. In the nonreacting flow, the mean and fluctuation profiles of a conserved scalar show good agreement with experimental measurements. Results for the reacting flow indicate that for temperature independent kinetics, the chemical reaction begins immediately downstream of the splitter plate where mixing starts. Results for the reacting flow with Arrhenius kinetics show an ignition delay, which depends on reactant temperature, before significant chemical reaction occurs. Harmonic forcing changes the structure of the layer, and concomitantly the rates of mixing and reaction, in accordance with experimental results. Strong stretch within the braids in the nonequilibrium kinetics case causes local flame quenching due to the temperature drop associated with the large convective fluxes.

Ghoniem, Ahmed F.↗