Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion problems”

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,081 records · Page 60

Determining pitch-angle diffusion coefficients from test particle simulations

The transport and acceleration of charged particles in turbulent media are topics of great interest in space physics and interstellar astrophysics. These processes are dominated by the scattering of particles off magnetic irregularities. The scattering process itself is usually described by small-angle scattering, with the pitch-angle coefficient D{sub μμ} playing a major role. Since the diffusion coefficient D{sub μμ} can be determined analytically only for the approximation of quasilinear theory, the determination of this coefficient from numerical simulations has become more important. So far these simulations have yielded particle tracks for small-scale scattering, which can then be interpreted using the running diffusion coefficients. This method has a limited range of validity. This paper presents two new methods that allow for the calculation of the pitch-angle diffusion coefficient from numerical simulations. These methods no longer analyze particle trajectories and instead examine the change of particle distribution functions. It is shown that these methods provide better resolved results and allow for the analysis of strong turbulence. The application of these methods to Monte Carlo simulations of particle scattering and hybrid MHD-particle simulations is presented. Both analysis methods are able to recover the diffusion coefficients used as input for the Monte Carlo simulations and provide better results in MHD simulations, especially for stronger turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗

HEATING OF THE PARTIALLY IONIZED SOLAR CHROMOSPHERE BY WAVES IN MAGNETIC STRUCTURES

In this paper, we show a “proof of concept” of the heating mechanism of the solar chromosphere due to wave dissipation caused by the effects of partial ionization. Numerical modeling of non-linear wave propagation in a magnetic flux tube, embedded in the solar atmosphere, is performed by solving a system of single-fluid quasi-MHD equations, which take into account the ambipolar term from the generalized Ohm’s law. It is shown that perturbations caused by magnetic waves can be effectively dissipated due to ambipolar diffusion. The energy input by this mechanism is continuous and shown to be more efficient than dissipation of static currents, ultimately leading to chromospheric temperature increase in magnetic structures.

79 ASTRONOMY AND ASTROPHYSICS↗

On optimal nonlinear estimation. I - Continuous observation.

A generalization of Bucy's (1965) representation theorem is obtained under very weak hypotheses. The generalized theorem is shown to play the same role in the case of general optimal estimation for an arbitrary random process as does the Bucy theorem in the case of optimal filtering for a diffusion process. At least for the models considered, the possibility is pointed out to reduce all sequential estimation problems to the problem of filtering. Hence, filtering theory is seen to represent the core of estimation theory, and is believed to define the direction in which future research should be focused.

Lo, J. T.↗

Unsteady, combined radiation and conduction in an absorbing, scattering, and emitting medium.

The transient cooldown of a gray, absorbing, isotropic scattering, emitting, and conducting medium bounded by gray, diffusely emitting and reflecting parallel plates is considered. Numerical solutions are obtained for the initial boundary-value problem with a discontinuous decrease in temperature at one boundary. The quasi-steady equation of radiative transfer is solved using Gaussian quadrature and a matrix eigenvector technique together with explicit numerical solution of the unsteady energy equation. Temperature and energy flux distributions are presented for variations of optical thickness, boundary emissivity, albedo, and conduction-radiation parameter.

Weston, K. C.↗

Interpretation of recent positron-electron measurements between 20 and 800 MeV

Recently measured positron and negatron spectra are discussed with regard to the problem of solar modulation. At energies above 180 MeV, the spherically symmetric Fokker-Planck equation with a diffusion coefficient proportional to particle rigidity provides reasonable fits to both the positron and total electron data. At energies below 180 MeV, the data are consistent with a continuation of the same diffusion coefficient and a local source of negatrons or with a change in the diffusion coefficient to a constant value.

Pellerin, C. J.↗

The prospects for composites based on boron fibers

The fabrication of boron filaments and the production of composite materials consisting of boron filaments and organic or metallic matrices are discussed. Problem involving the use of tungsten substrates in the filament fabrication process, the protection of boron fibers with diffusion barrier cladings, and the application of alloy additives in the matrix to lessen the effects of diffusion are considered. Data on the kinetics of the boron fiber/matrix interaction at high temperatures, and the influence of the fiber/matrix interaction on the mechanical properties of the composite are presented.

Naslain, R.↗

Combustor modelling for scramjet engines

A system of computer programs is being developed to analyse and predict the complex flow fields found in hydrogen-fueled scramjet combustors. Each program is designed to solve the governing equation system for the type of flow present in a particular combustor region. A two-dimensional parabolic program has been found to be valuable in the development and experimental evaluation of turbulence and chemistry models for supersonic flow, and in the development of a program to model supersonic flow downstream of the fuel injection struts by means of solutions to the three-dimensional parabolic Navier-Stokes equations and species equations. A partially elliptic code has been derived to account for local subsonic flow regions, and fully elliptic programs have been developed by the consideration of streamwise diffusion effects for the recirculating flow fields near transverse fuel injectors. The programs are currently being applied to problems of scramjet engine development.

Drummond, J. P.↗

A numerical study of the pitch-angle scattering of cosmic rays

The results are presented of a careful study of finite-difference solutions to the problem of cosmic-ray transport, including pitch-angle scattering. In contrast to some recent studies, the diffusion approximation is confirmed in cases where the scattering mean free path is small compared with other length scales. The reasons for the discrepancy with the conclusions of Gombosi and Owens are discussed.

Kota, J.↗

The thermal stability of coronal loops by nonlinear diffusion asymptotics

A nonlinear reaction-diffusion equation and some additional constraints are derived which describe the time-dependent behavior of the temperature structure of the plasma in coronal loops. The equation is analyzed using nonlinear diffusion asymptotics, in particular singular perturbation techniques, and the results are interpreted in the context of the physical problem of the thermal stability and temporal behavior of the plasma. The results are consistent with the possibility of cyclic thermal behavior of the plasma, as suggested by Kuin and Martens (1982).

Pakkert, J. W.↗

A hyperbolic-equation system approach for magnetized electron fluids in quasi-neutral plasmas

A new approach using a hyperbolic-equation system (HES) is proposed to solve for the electron fluids in quasi-neutral plasmas. The HES approach avoids treatments of cross-diffusion terms which cause numerical instabilities in conventional approaches using an elliptic equation (EE). A test calculation reveals that the HES approach can robustly solve problems of strong magnetic confinement by using an upwind method. The computation time of the HES approach is compared with that of the EE approach in terms of the size of the problem and the strength of magnetic confinement. The results indicate that the HES approach can be used to solve problems in a simple structured mesh without increasing computational time compared to the EE approach and that it features fast convergence in conditions of strong magnetic confinement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modification of the flow pass method as applied to problems of chemistry of planet atmospheres

It was shown that the modified flow pass method possesses considerable effectiveness, both in the case when the coefficient of diffusion changes severely in the examined region and in the case when diffusion is the prevalent process, as compared with chemical reactions. The case when a regular pass proves inapplicable, or applicable in a limited interval of the decisive parameters, was examined.

Parshev, V. A.↗

Modeling of Nonlocal Electron Conduction for Inertial Confinement Fusion

In this report Chenhall presents a new Monte Carlo (MC) transport method developed as part of his PhD thesis work funded via fellowship by Sandia National Laboratories. The model is used to simulate the electron thermal transport within inertial confinement fusion (ICF) type problems. The new model aims to improve upon the currently used implicit Schurtz, Nicolai, and Busquet (iSNB) method, in particular by using finite particle ranges in comparison to the exponential solution of a diffusion method and by improved higher order angular modeling. The new method is being developed using the 1D LILAC and 2D DRACO codes developed by the Rochester Laboratory for Laser Energetics. The model is compared to iSNB for several ICF type problems: Omega shot 60303 a shock timing experiment, Omega shot 68951 a cryo target implosion and a NIF phase plate polar direct drive design. Overall, the MC method performs at least as well as the iSNB method and appears to indicate a lower predicted preheat ahead of the shock fronts. However, the difference in overall results between iSNB and MC transport is small for the choice of test problems and further simulations will need to be run to fully understand the difference between the models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Recent Progress on Numerical Modeling for Microgravity Electric Field Flames: Preprint

This paper presents the recent progress on the NASA Physical Science Informatics (PSI) project PeleLM CFD of Ion Driven Winds from Diffusion Flames in simulating the E-FIELD Flames microgravity results. The focus of this project is to comprehensively simulate the behavior of a small diffusion flame under the influence of an externally applied electric field in zero-gravity. To date, the capability of accurately simulating electric field flames has eluded researchers because the system exhibits dramatic ranges of coupled temporal and spatial scales. Moreover, in earth gravity the hot combustion products are subject to buoyancy effects that are difficult to isolate from those generated by the electric field. This work is an implementation in an existing powerful simulation framework (PeleLM) for this problem using PeleLMeX in order to validate the model and investigate the complex coupled system. The current effort includes establishing the domain of a diffusion coflow burner, examining the boundary conditions, flame geometry and ignition with gravitational forces and also with an electric field applied. A detailed model is used that includes the chemistry of charged ions and chemiluminescent flame intermediates to capture any feedback between ion-driven convection and combustion behavior, and to allow quantitative comparisons with experimental measurements.

coflow flames↗

Diffusion of Innovations: Interplay of Social, Economic, Technological, and Policy Drivers in the Solar Industry. Summary of UT Austin Student Capstone Research Projects

The University of Texas at Austin’s Policy Research Project (PRP), a nine-month (two semesters) capstone, is a keystone of the core curriculum at the LBJ School of Public Affairs. In PRPs, small groups of students, under the mentorship of a faculty director, take on real-world problems that require special knowledge and skill sets. PRPs expose students to challenges in formulating and executing research, and in communicating academic research and related complex data to broader stakeholder communities and decision makers. The PRP structure is an innovative and effective approach for integrating research within the teaching and training of graduate students who are preparing themselves to address important real-world problems at the intersection of society, economics, technology, and policy. The project summaries below describe seven papers developed during September 2017 – May 2018 as part of a PRP on “Diffusion of Innovations: Interplay of Social, Economic, Technological, and Policy Drivers in the Solar Industry.” Twenty graduate students, drawn from the LBJ School’s Masters in Public Affairs and Masters in Global Policy Studies programs and the Jackson School Geoscience’s Energy and Earth Resources program, participated in this PRP. Dr. Varun Rai, Associate Professor and Associate Dean for Research at the LBJ School, directed the PRP, with support from his research team including: Dr. Ariane Beck, Dr. Ashok Sekar, D. Cale Reeves, and Erik Funkhouser. Clients for the project included the U.S. Department of Energy (Casey Canfield), Lawrence Berkeley National Laboratory (Ben Hoen, Galen Barbose Joachim Seel, Naïm Darghouth, Ryan Wiser), and National Renewable Energy Laboratory (Benjamin Sigrin, Eric O’Shaughnessy). The seven projects separately addressed one of the following topics: (1) low- and middle-income PV adoption, (2) modeling economic and information intervention design, (3) evaluation of DOE’s Solar in Your Community Challenge, (4) property value impacts near large-scale solar facilities, (5) solar market maturity and evolution of business models, (6) social media data for predicting PV adoption, and (7) individual-level variation in adoption of innovations. Many of the papers relied on data collected and curated by Lawrence Berkeley National Laboratory, including data embedded within the annual Tracking the Sun and Utility-Scale Solar reports. Each of the seven teams in the PRP prepared a research paper. The PRP culminated with a full-day conference at UT Austin in May 2018 to present findings from the seven projects in this PRP to a broad audience of about 75 experts from academia, national labs, industry, and government from across the country.

14 SOLAR ENERGY↗

Kinetic modelling for zinc (II) ions biosorption onto Luffa cylindrica

The biosorption of Zinc (II) ions onto a biomaterial - Luffa cylindrica has been studied. This biomaterial was characterized by elemental analysis, surface area, pore size distribution, scanning electron microscopy, and the biomaterial before and after sorption, was characterized by Fourier Transform Infra Red (FTIR) spectrometer. The kinetic nonlinear models fitted were Pseudo-first order, Pseudo-second order and Intra-particle diffusion. A comparison of non-linear regression method in selecting the kinetic model was made. Four error functions, namely coefficient of determination (R{sup 2}), hybrid fractional error function (HYBRID), average relative error (ARE), and sum of the errors squared (ERRSQ), were used to predict the parameters of the kinetic models. The strength of this study is that a biomaterial with wide distribution particularly in the tropical world and which occurs as waste material could be put into effective utilization as a biosorbent to address a crucial environmental problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Competitive autocatalytic reactions in chaotic flows with diffusion: Prediction using finite-time Lyapunov exponents

We investigate chaotic advection and diffusion in autocatalytic reactions for time-periodic sine flow computationally using a mapping method with operator splitting. We specifically consider three different autocatalytic reaction schemes: a single autocatalytic reaction, competitive autocatalytic reactions, which can provide insight into problems of chiral symmetry breaking and homochirality, and competitive autocatalytic reactions with recycling. In competitive autocatalytic reactions, species B and C both undergo an autocatalytic reaction with species A such that A+B→2B and A+C→2C. Small amounts of initially spatially localized B and C and a large amount of spatially homogeneous A are advected by the velocity field, diffuse, and react until A is completely consumed and only B and C remain. We find that local finite-time Lyapunov exponents (FTLEs) can accurately predict the final average concentrations of B and C after the reaction completes. The species that starts in the region with the larger FTLE has, with high probability, the larger average concentration at the end of the reaction. If B and C start in regions with similar FTLEs, their average concentrations at the end of the reaction will also be similar. When a recycling reaction is added, the system evolves towards a single species state, with the FTLE often being useful in predicting which species fills the entire domain and which is depleted. The FTLE approach is also demonstrated for competitive autocatalytic reactions in journal bearing flow, an experimentally realizable flow that generates chaotic dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Numerical experiments on a new class of nonoscillatory schemes

Numerical experiments for the SONIC schemes on 2D inviscid, compressible, steady, and unsteady problems are presented. These schemes belong to a new class of uniformly second-order accurate nonoscillatory schemes introduced by Huynh, with the well known UNO2 scheme of Harten and Osher being the most 'diffusive' in this class. The SONIC schemes can also be considered as uniformly second order accurate extensions of the popular TVD schemes. For simplicity, a MUSCL approach for spatial discretization and a Runge-Kutta method for time integration are used. Test problems include steady oblique shock reflection and the well known unsteady double Mach reflection problem. Results confirm that the SONIC schemes are more accurate than their TVD counterparts.

Suresh, Ambady↗