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 541 records · Page 30

Fast methods incorporating direct elliptic solvers for nonlinear applications in fluid dynamics

Semidirect methods are discussed, their present role, as well as some developments for their application in computational fluid dynamics. A semidirect method is a computational scheme that uses a fast, direct, elliptic solver as the driving algorithm for the iterative solution of finite difference equations. Specific subtopics include: (1) direct Cauchy Riemann solvers for first order elliptic equations; (2) application of the semidirect method to the mixed elliptic hyperbolic problem of steady, inviscid transonic flow; and (3) the treatment of interior conditions, such as those on an airfoil or wing, in semidirect methods.

Martin, E. D.↗

A crucible salt saturation method for determining sulfur solubility in glass melt

Experiments were conducted to determine sulfur solubility in Hanford low-activity waste (LAW) glass melts by a sulfur saturation method. Sulfur-incorporated glass melts were prepared by salt saturation and bubbling methods. The salt saturation method was performed by mixing crushed pre-melted baseline glasses with an excess amount of Na2SO4 prior to melting the mixture at 1150°C for one hour. The bubbling method involved bubbling the glass melt at 1150°C in a Pt crucible with an SO2/O2/N2 gas mix to equilibrate the melt at a known pressure of SO3. Preliminary results suggested that performing one cycle of mixing and melting was not sufficient to saturate the glass. The bubbling method successfully incorporated sulfur into the glass but caused significant losses of sodium from the melt. In order to saturate the glass melt with sulfate without causing noticeable sodium loss, a modified crucible salt saturation method was developed by repeating the mixing and melting of the glass and salt mixture. For all three representative LAW glasses tested in this study, it was found that after three mixing and melting cycles, the sulfur concentration reaches a plateau, indicating reasonable sulfur saturation.

Jin, Tongan↗

Subgrid Combustion Modeling for the Next Generation National Combustion Code

In the first year of this research, a subgrid turbulent mixing and combustion methodology developed earlier at Georgia Tech has been provided to researchers at NASA/GRC for incorporation into the next generation National Combustion Code (called NCCLES hereafter). A key feature of this approach is that scalar mixing and combustion processes are simulated within the LES grid using a stochastic 1D model. The subgrid simulation approach recovers locally molecular diffusion and reaction kinetics exactly without requiring closure and thus, provides an attractive feature to simulate complex, highly turbulent reacting flows of interest. Data acquisition algorithms and statistical analysis strategies and routines to analyze NCCLES results have also been provided to NASA/GRC. The overall goal of this research is to systematically develop and implement LES capability into the current NCC. For this purpose, issues regarding initialization and running LES are also addressed in the collaborative effort. In parallel to this technology transfer effort (that is continuously on going), research has also been underway at Georgia Tech to enhance the LES capability to tackle more complex flows. In particular, subgrid scalar mixing and combustion method has been evaluated in three distinctly different flow field in order to demonstrate its generality: (a) Flame-Turbulence Interactions using premixed combustion, (b) Spatially evolving supersonic mixing layers, and (c) Temporal single and two-phase mixing layers. The configurations chosen are such that they can be implemented in NCCLES and used to evaluate the ability of the new code. Future development and validation will be in spray combustion in gas turbine engine and supersonic scalar mixing.

Menon, Suresh↗

Mixed-metal oxide based catalysts

One or more embodiments relates to a method of catalytically converting a reactant gas mixture for pollution abatement of products of hydrocarbon fuel combustion. The method provides substituted mixed-metal oxides where catalytically active metals are substituted within the crystal lattice to create an active and well dispersed metal catalyst available to convert the reactant gas mixture. Embodiments may be used with gasoline and diesel fueled internal combustion engine exhaust, although specific embodiments may differ somewhat for each.

Berry, David A.↗

An eddy viscosity calculation method for a turbulent duct flow

The mean velocity profile across a fully developed turbulent duct flow is obtained from an eddy viscosity relation combined with an empirical outer region wake function. Results are in good agreement with experiments and with direct numerical simulations in the same flow at two Reynolds numbers. In particular, the near-wall trend of the Reynolds shear stress and its variation with Reynolds number are similar to those of the simulations. The eddy viscosity method is more accurate than previous mixing length or implicit function methods.

Antonia, R. A.↗

Mechanical and Combustion Performance of Multi-Walled Carbon Nanotubes as an Additive to Paraffin-Based Solid Fuels for Hybrid Rockets

Paraffin-based solid fuels for hybrid rocket motor applications are recognized as a fastburning alternative to other fuel binders such as HTPB, but efforts to further improve the burning rate and mechanical properties of paraffin are still necessary. One approach that is considered in this study is to use multi-walled carbon nanotubes (MWNT) as an additive to paraffin wax. Carbon nanotubes provide increased electrical and thermal conductivity to the solid-fuel grains to which they are added, which can improve the mass burning rate. Furthermore, the addition of ultra-fine aluminum particles to the paraffin/MWNT fuel grains can enhance regression rate of the solid fuel and the density impulse of the hybrid rocket. The multi-walled carbon nanotubes also present the possibility of greatly improving the mechanical properties (e.g., tensile strength) of the paraffin-based solid-fuel grains. For casting these solid-fuel grains, various percentages of MWNT and aluminum particles will be added to the paraffin wax. Previous work has been published about the dispersion and mixing of carbon nanotubes.1 Another manufacturing method has been used for mixing the MWNT with a phenolic resin for ablative applications, and the manufacturing and mixing processes are well-documented in the literature.2 The cost of MWNT is a small fraction of single-walled nanotubes. This is a scale-up advantage as future applications and projects will require low cost additives to maintain cost effectiveness. Testing of the solid-fuel grains will be conducted in several steps. Dog bone samples will be cast and prepared for tensile testing. The fuel samples will also be analyzed using thermogravimetric analysis and a high-resolution scanning electron microscope (SEM). The SEM will allow for examination of the solid fuel grain for uniformity and consistency. The paraffin-based fuel grains will also be tested using two hybrid rocket test motors located at the Pennsylvania State University s High Pressure Combustion Lab.

Larson, Daniel B.↗

Composite interlaminar fracture toughness: Three-dimensional finite element modeling for mixed mode 1, 2 and 3 fracture

A computational method/procedure is described which can be used to simulate individual and mixed mode interlaminar fracture progression in fiber composite laminates. Different combinations of Modes 1, 2, and 3 fracture are simulated by varying the crack location through the specimen thickness and by selecting appropriate unsymmetric laminate configurations. The contribution of each fracture mode to strain energy release rate is determined by the local crack closure methods while the mixed mode is determined by global variables. The strain energy release rates are plotted versus extending crack length, where slow crack growth, stable crack growth, and rapid crack growth regions are easily identified. Graphical results are presented to illustrate the effectiveness and versatility of the computational simulation for: (1) evaluating mixed-mode interlaminar fracture, (2) for identifying respective dominant parameters, and (3) for selecting possible simple test methods.

Murthy, P. L. N.↗

Composite interlaminar fracture toughness - Three-dimensional finite-element modeling for mixed mode I, II, and fracture

A computational method/procedure is described which can be used to simulate individual and mixed mode interlaminar fracture progression in fiber composite laminates. Different combinations of Modes 1, 2, and 3 fracture are simulated by varying the crack location through the specimen thickness and by selecting appropriate unsymmetric laminate configurations. The contribution of each fracture mode to strain energy release rate is determined by the local crack closure methods while the mixed mode is determined by global variables. The strain energy release rates are plotted versus extending crack length, where slow crack growth, stable crack growth, and rapid crack growth regions are easily identified. Graphical results are presented to illustrate the effectiveness and versatility of the computational simulation for: (1) evaluating mixed-mode interlaminar fracture, (2) for identifying respective dominant parameters, and (3) for selecting possible simple test methods.

Murthy, Pappu L. N.↗

Low external pH replaces 2,4-D in maintaining and multiplying 2,4-D-initiated embryogenic cells of carrot

A mixed culture comprised of both embryonic globules and nonembryogenic callus was derived from seedling hypocotyls of Daucus carota cv. Scarlet Nantes on 2,4-D- containing medium using well-established methods. Then the mixed cultures were transferred to, and serially subcultured on, a hormone-free medium near pH 4. The medium contained 1 mM NH4+ as the sole nitrogen source. When cultured in this way, embryonic globules were able to multiply without development into later embryo stages. Nonembryogenic callus did not survive. Continuous culture of embryonic globules on this low pH hormone-free medium yielded cultures consisting entirely of preglobular stage proembryos (PGSPs). PGSP cultures have been maintained as such with continuous multiplication for nearly 2 years without loss of embryogenic potential. These hormone-free-maintained PGSPs continue their development to later embryo stages when cultured on the same hormone-free medium buffered at pH 5.8. We show that hormone-free medium near pH 4 can replace 2,4-D in its ability to sustain multiplication of 2,4-D-initiated embryogenic cells of carrot at an acceptable growth rate without their development into later embryo stages. This procedure provides selective conditions that do not permit the growth of non-embryogenic cells while providing an adequate environment for embryogenic cell proliferation and should prove invaluable in studying habituation.

Non-NASA Center↗

A stabilized Runge–Kutta–Legendre method for explicit super-time-stepping of parabolic and mixed equations

Parabolic partial differential equations appear in several physical problems, including problems that have a dominant hyperbolic part coupled to a sub-dominant parabolic component. Explicit methods for their solution are easy to implement but have very restrictive time step constraints. Implicit solution methods can be unconditionally stable but have the disadvantage of being computationally costly or difficult to implement. Super-time-stepping methods for treating parabolic terms in mixed type partial differential equations occupy an intermediate position. In such methods each superstep takes “s” explicit Runge–Kutta-like time-steps to advance the parabolic terms by a time-step that is s{sup 2} times larger than a single explicit time-step. The expanded stability is usually obtained by mapping the short recursion relation of the explicit Runge–Kutta scheme to the recursion relation of some well-known, stable polynomial. Prior work has built temporally first- and second-order accurate super-time-stepping methods around the recursion relation associated with Chebyshev polynomials. Since their stability is based on the boundedness of the Chebyshev polynomials, these methods have been called RKC1 and RKC2. In this work we build temporally first- and second-order accurate super-time-stepping methods around the recursion relation associated with Legendre polynomials. We call these methods RKL1 and RKL2. The RKL1 method is first-order accurate in time; the RKL2 method is second-order accurate in time. We verify that the newly-designed RKL1 and RKL2 schemes have a very desirable monotonicity preserving property for one-dimensional problems – a solution that is monotone at the beginning of a time step retains that property at the end of that time step. It is shown that RKL1 and RKL2 methods are stable for all values of the diffusion coefficient up to the maximum value. We call this a convex monotonicity preserving property and show by examples that it is very useful in parabolic problems with variable diffusion coefficients. This includes variable coefficient parabolic equations that might give rise to skew symmetric terms. The RKC1 and RKC2 schemes do not share this convex monotonicity preserving property. One-dimensional and two-dimensional von Neumann stability analyses of RKC1, RKC2, RKL1 and RKL2 are also presented, showing that the latter two have some advantages. The paper includes several details to facilitate implementation. A detailed accuracy analysis is presented to show that the methods reach their design accuracies. A stringent set of test problems is also presented. To demonstrate the robustness and versatility of our methods, we show their successful operation on problems involving linear and non-linear heat conduction and viscosity, resistive magnetohydrodynamics, ambipolar diffusion dominated magnetohydrodynamics, level set methods and flux limited radiation diffusion. In a prior paper (Meyer, Balsara and Aslam 2012 [36]) we have also presented an extensive test-suite showing that the RKL2 method works robustly in the presence of shocks in an anisotropically conducting, magnetized plasma.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Domain-Decomposed A-ϕ Formulation Based on Lagrange Multipliers for Low-Frequency Problems

A domain-decomposed A-ϕ formulation based on Lagrange multipliers is proposed to simulate low-frequency elec- tromagnetic problems. This method partitions the computational domain into smaller subdomains, allowing each subdomain to be independently formulated using Lagrange multipliers as Dirichlet boundary conditions, while ensuring continuity of the fields across the interfaces. A mixed finite element method, utilizing both vector and scalar basis functions, is employed to discretize the formulation, resulting in a global system to be solved. The proposed method is validated using TEAM Problem 7 at 50 Hz, demonstrating its effectiveness in handling complex geometries and addressing the low-frequency breakdown issues commonly encountered in traditional finite element methods.

Hossain, Amzad↗

Turbulent mixing in a precessing sphere

By numerically simulating turbulent flows at high Reynolds numbers in a precessing sphere, we propose a method to enhance the mixing of a fluid confined within a smooth cavity by its rotational motion alone. To precisely evaluate the mixing efficiency, we extend the quantification method proposed by Danckwerts [“The definition and measurement of some characteristics of mixtures,” Appl. Sci. Res. A 3, 279–296 (1952)] to the case in which only a finite number of fluid particle trajectories can be known. Our accurate numerical tracking of fluid particles in the flow, which is controlled by the Reynolds number (an indicator of the spin rate) and the Poincaré number (the precession rate), shows the following results. First, the mixing process on the time scale normalized by the spin period is independent of the Reynolds number as long as it is high enough for the flow to be developed turbulence. Second, fastest mixing is achieved under weak precession (Poincaré number ≈0.1); in such cases, perfect mixing requires only 10–15 spins of the container. Third, the power to sustain turbulence is a weakly increasing function of the Poincaré number, and the energy efficiency of the mixing is also maximized when the Poincaré number is about 0.1. Fourth, efficient mixing driven by the weak precession arises from the effective cooperation of complex large-scale flow and small-scale turbulence, which itself is sustained by the large-scale flow.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Chaotic advection at large Péclet number: Electromagnetically driven experiments, numerical simulations, and theoretical predictions

We present a combination of experiment, theory, and modelling on laminar mixing at large Péclet number. The flow is produced by oscillating electromagnetic forces in a thin electrolytic fluid layer, leading to oscillating dipoles, quadrupoles, octopoles, and disordered flows. The numerical simulations are based on the Diffusive Strip Method (DSM) which was recently introduced (P. Meunier and E. Villermaux, “The diffusive strip method for scalar mixing in two-dimensions,” J. Fluid Mech. 662, 134–172 (2010)) to solve the advection-diffusion problem by combining Lagrangian techniques and theoretical modelling of the diffusion. Numerical simulations obtained with the DSM are in reasonable agreement with quantitative dye visualization experiments of the scalar fields. A theoretical model based on log-normal Probability Density Functions (PDFs) of stretching factors, characteristic of homogeneous turbulence in the Batchelor regime, allows to predict the PDFs of scalar in agreement with numerical and experimental results. This model also indicates that the PDFs of scalar are asymptotically close to log-normal at late stages, except for the large concentration levels which correspond to low stretching factors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Determination of normal-shock position in a mixed compression supersonic inlet

Methods of determining inlet normal shock position from wall static pressure profiles were investigated. By using methods investigated in this report, it should be possible to control an inlet with less stability margin. Variation in inlet angle-of-attack caused drastic changes in pressure profile shape and wide variations in the angle of the shock plane with respect to the inlet axis. Thus, four sensors located around the circumference may be required to handle angle-of-attack and sideslip variations. One criterion would determine the shock position by comparing individual statics to a properly selected reference pressure.

Dustin, M. O.↗

A simple method for estimation of coagulation efficiency in mixed aerosols

Aerosols of KBr and AgNO3 were mixed, exposed to light in a glass tube and collected in the dark. About 15% of the collected material was reduced to silver upon development. Thus, two aerosols of particles that react to form a photo-reducible compound can be used to measure coagulation efficiency.

Dimmick, R. L.↗