SEARCH · Search NASA

Results for “numerical dispersion”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Axially symmetric pulse propagation in semi-infinite hollow cylinders

The effects of dispersion and end conditions on axially symmetric pulse propagation in semiinfinite tubes are investigated. Numerical modeling of the dynamic response is accomplished using a stable and efficient finite element and wave propagation based method. This method yields the dispersive and modal data required to express displacements in a frequency domain modal expansion. Boundary conditions are applied to the tube end to obtain mode amplitudes. A fast Fourier transform is then used to get the time response. The specific end conditions considered are intended to give some insight into how the tube would behave as a member of a jointed truss structure.

Kohl, T.↗

DOE OSTI2025

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1977

Guidance and trajectory considerations in lunar mass transportation

Flight-mechanics problems associated with large-scale transport of lunar mass to a space colony or manufacturing facility are discussed. The proposed transport method involves launch of payloads from a mass-driver on the lunar surface, onto ballistic trajectories to a passive mass-catcher located near the L2 libration point, with the caught mass subsequently being transported to the colony. Arrival velocities at L2, sensitivities in arrival dispersion due to launch errors, and effects of launch site location are treated, via numerically integrated orbits in the restricted three-body problem. From any launch site it is possible to define a target point reached with zero dispersion due to errors in a selected component of launch velocity. Effects of lunar geometrical librations and of obliquity, as well as the conditions for biasing a trajectory away from L2 so as to reduce stationkeeping costs, are dealt with along with transfer orbits from L2 to the colony. The theory of capture and the theory of resonance lead to a colony orbit, with period approximately two weeks, reached from L2 with velocity increment as low as 9.02 m/sec.

Heppenheimer, T. A.↗

DOE OSTI2024

Vibrational Fermi Resonance in Atomically Thin Black Phosphorus

Fermi resonance is a phenomenon involving the hybridization of two coincidentally quasi-degenerate states that is observed in the vibrational or electronic spectra of molecules. Despite numerous examples in molecular systems, vibrational Fermi resonances in dispersive semiconducting systems remain largely unexplored due to the rarity of occurrence. Here we report a vibrational Fermi resonance in atomically thin black phosphorus. The Fermi resonance arises via anharmonic mixing of a fundamental Raman mode and a Davydov component of an infrared mode, leading to a doublet with mixed character. The extent of Fermi coupling can be modulated by the application of external biaxial strain. The consequences of Fermi hybridization are revealed by electronic resonance effects in the thickness-dependent and excitation-wavelength-dependent Raman spectrum, which is predicted by ab initio hybrid functional simulations including excitonic interactions. This work reveals new insight into electron–phonon coupling in black phosphorus and demonstrates a novel method for modulating Fermi resonances in 2D semiconductors.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Shock formation in flowing plasmas by temporally and spatially smoothed laser beams

The cumulative impact of multiple laser speckles on a supersonic plasma flow across optically smoothed laser beams is investigated. The bending of laser beams caused by ponderomotive laser–plasma coupling, together with flow, leads to plasma a momentum-conserving response that results in a deceleration of the flow. Once the flow velocity decreases to a subsonic level, the action of the laser beams can generate a shock within the plasma. This scenario has been predicted theoretically and confirmed by hydrodynamic simulations. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size, and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging only in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser–plasma coupling in the inertial confinement fusion experiments are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Modeling and simulations of hydrodynamic shocks in a plasma flowing across randomized ICF scale laser beams

High-energy laser beams interacting with flowing plasmas can produce a plasma response that leads to deflection of the beam, beam bending. Such beams have usually a speckle structure generated by optical smoothing techniques that reduce the spatial and temporal coherence in the laser field pattern. The cumulative plasma response from laser speckles slows down the velocity of the incoming flow by momentum conservation. For slightly super-sonic flow the cumulative plasma response to the ponderomotive force exerted by the beam speckle ensemble is the strongest, such that slowing down the flow to subsonic velocities leads eventually to the generation of a shock around the cross section of the beam. This scenario has been predicted theoretically and is confirmed here by our hydrodynamic simulations in two dimensions with speckled beams and in one dimension with a reduced model. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed laser beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging from those produced by beams with only a random phase plates in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser plasma coupling in inertial confinement fusion (ICF) experiments are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1985

Numerical methods for analyzing electromagnetic scattering

Numerical methods to analyze electromagnetic scattering are presented. The dispersions and attenuations of the normal modes in a circular waveguide coated with lossy material were completely analyzed. The radar cross section (RCS) from a circular waveguide coated with lossy material was calculated. The following is observed: (1) the interior irradiation contributes to the RCS much more than does the rim diffraction; (2) at low frequency, the RCS from the circular waveguide terminated by a perfect electric conductor (PEC) can be reduced more than 13 dB down with a coating thickness less than 1% of the radius using the best lossy material available in a 6 radius-long cylinder; (3) at high frequency, a modal separation between the highly attenuated and the lowly attenuated modes is evident if the coating material is too lossy, however, a large RCS reduction can be achieved for a small incident angle with a thin layer of coating. It is found that the waveguide coated with a lossy magnetic material can be used as a substitute for a corrugated waveguide to produce a circularly polarized radiation yield.

Lee, S. W.↗

NASA NTRS1986

Compositional redistribution during casting of Hg sub 0.8 Cd sub 0.2 Te alloys

A series of Hg(0.8)Cd(0.2)Te ingots was cast both vertically and horizontally under well-defined thermal conditions by using a two-zone furnace with isothermal heat-pipe liners. The main objective of the experiments was to establish correlations between casting parameters and compositional redistribution and to develop ground-based data for a proposed flight experiment of casting of Hg(1-x)Cd(x)Te alloys under reduced gravity conditions. The compositional variations along the axial and radial directions were determined by precision density measurements, infrared transmission spectra, and X-ray energy dispersion spectrometry. Comparison between the experimental results and a numerical simulation of the solidification process of Hg(0.8)Cd(0.2)Te is described.

Su, Ching-Hua↗

NASA NTRS1988

Compositional redistribution during casting of Hg(0.8)Cd(0.2)Te alloys

Su, Ching-Hua↗

NASA NTRS1991

The formation of primordial binaries in globular clusters by star-disk interactions

The formation of primordial binaries in globular clusters is examined using simple numerical models. Clusters of protostars collapse until their velocity dispersion rises sufficiently to reverse the infall and the cluster reaches equilibrium. During the collapse, interactions between stars and protostellar disks lead to stellar capture. It is found that binary fraction of a few percent typically result. Binary formation is terminated when the velocity dispersion rises to a point at which most encounters result in disk destruction rather than capture. As a result, much gas is returned to the cluster ISM, limiting the star formation efficiency to a value significantly below 100 percent.

Murray, Stephen D.↗

NASA NTRS2000

Exchange Forces in Dispersion Relations Investigated Using Circuit Relations

We propose a novel method to compute in an exact manner the left-hand cut discontinuity of the electron-Atom partial wave scattering amplitude in the complex energy plane within the static stage approximation. Zero energy dispersion relations for electron-Hydrogen scattering are computed numerically for illustration.

Vrinceanu, D.↗

DOE OSTI2021

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗

NASA NTRS2004

An Investigation of Wave Propagations in Discontinuous Galerkin Method

Analysis of the discontinuous Galerkin method has been carried out for one- and two-dimensional system of hyperbolic equations. Analytical, as well as numerical, properties of wave propagation in a DGM scheme are derived and verified with direct numerical simulations. In addition to a systematic examination of the dissipation and dispersion errors, behaviours of a DG scheme at an interface of two different grid topologies are also studied. Under the same framework, a quantitative discrete analysis of various artificial boundary conditions is also conducted. Progress has been made in numerical boundary condition treatment that is closely related to the application of DGM in aeroacoustics problems. Finally, Fourier analysis of DGM for the Convective diffusion equation has also be studied in connection with the application of DG schemes for the Navier-Stokes equations. This research has resulted in five(5) publications, plus one additional manuscript in preparation, four(4) conference presentations, and three(3) departmental seminars, as summarized in part II. Abstracts of papers are given in part 111 of this report.

Hu, Fang Q.↗

NASA NTRS1982

A parameterized model for the evolution of isotopic heterogeneities in a convecting system

It is experimentally shown that, although steady convective flows are efficient means to heterogeneity within a single cell, they do not produce a dispersal of heterogeneous material over scales that are large by comparison to their depth, which requires that the flow be time-dependent on a time scale comparable to the overturn time. Convection in an internally heated layer does possess this property, and numerical solutions are presently used to study the way in which it disperses a set of neutrally bouyant particles initially confined to a small space. The derived concept of effective diffusivity is applied to the isotopic evolution of the Sm-Nd and Rb-Sr systems, with spatial variations generated by horizontal variations in degree of melting 1.8 billion years ago.

Richter, F. M.↗

NASA NTRS2000

Modeling of Ultrasonically Generated Liquid-Liquid Dispersions During Controlled Directional Solidification

There are innumerable two-component systems in which two very different liquid phases co-exist in equilibrium over a range of temperature and composition, e.g., oil and water, salt fluxes and solders, aluminum and lead. Often it is of practical concern to fabricate a solid component consisting of a uniform dispersion of one phase in the other. Unfortunately, uniform microstructural development during solidification of two immiscible liquids is hampered by inherent, often large, density differences between the phases that lead to severe segregation. Uniformity is also compromised by preferential wetting and coalescence phenomena. It is, however, well known that ultrasonic energy can initiate and maintain a fine liquid-liquid dispersion. The work presented here extends that observation by application of ultrasonic energy to promote uniform phase incorporation during controlled directional solidification. To this end experiments with the transparent organic, immiscible, succinonitrile-glycerol system were conducted and the numerous processing parameters associated with this technique were evaluated in view of optimizing dispersion uniformity. In view of the initial experimental results a model that predicts the dispersed liquid droplet size as a function of material properties, sample geometry, and applied energy has been developed. In the mathematical model we consider the ultrasonic field in an experimental ampoule of length L and diameter D induced by a probe having a vibration frequency of f=2OKhz (circular frequency omega = 2 pi f). The amplitude is adjustable from A=65 to 13Omicrons. The probe tip diameter is d, the liquid has a density of p, in which the speed of sound and surface tension are, respectively, c and sigma. The mathematical model and numerical investigation for the experiments [1] is done using the following assumptions: (i) The droplet size is small in comparison to the sound wave length; (ii) The forces between droplets are neglected (relative concentration is small); (iii) The droplet is stable if the kinetic energy, E(sub K), of the liquid motion due to ultrasonic field influence is less then the binding energy, E(sub S), due to the surface tension (it is easy to show that the surface energy of two droplets resulting from one is larger by about a factor of two.); (iv) The stability limit is characterized by E(sub S) to approx. E(sub K).

Grugel, R. N.↗

NASA NTRS1978

NASA Agricultural Aircraft Research Program in the Langley Vortex Research Facility and the Langley Full Scale Wind Tunnel

The current status of aerial applications technology research at the Langley's Vortex Research Facility and Full-Scale Wind Tunnel is reviewed. Efforts have been directed mainly toward developing and validating the required experimental and theoretical research tools. A capability to simulate aerial dispersal of materials from agricultural airplanes with small-scale airplane models, numerical methods, and dynamically scaled test particles was demonstrated. Tests on wake modification concepts have proved the feasibility of tailoring wake properties aerodynamically to produce favorable changes in deposition and to provide drift control. An aerodynamic evaluation of the Thrush Commander 800 agricultural airplane with various dispersal systems installed is described. A number of modifications intended to provide system improvement to both airplane and dispersal system are examined, and a technique for documenting near-field spray characteristics is evaluated.

Jordan, F. L., Jr.↗

NASA NTRS1991

A detailed numerical investigation of Burke-Schumann gaseous and spray flames

The classical Burke-Schumann gaseous diffusion flame, and the effect of introducing fuel in the form of liquid sprays are investigated numerically. The time-dependent two-dimensional, axisymmetric conservation equations are solved for the gas phase. A Lagrangian approach is used for the dispersed phase. The chemical reactions are modeled through the one-step global reaction scheme. The numerical results show that the computed flame structure is significantly different from that given by the analytical solution. The computed flame is about 15 percent longer and 5 percent narrower than the classical Burke-Schumann flame. The predicted flame shape is, however, in better agreement with experimental observations than the Burke-Schumann flame. The effects of introducing liquid sprays is found to produce thinner and longer flames than its counterpart gaseous flame. A detailed comparison of the structures of gaseous and spray flames is made, and the differences are discussed.

Mawid, M. A.↗

NASA NTRS1971

Waves along a cylindrical 'grid-like' antenna in an anisotropic compressible plasma.

The dispersion equation is obtained analytically of waves along an infinitely long cylindrical ?grid-like' antenna in an anisotropic compressible homogeneous plasma. Numerical solutions of this equation show that an external static magnetic field has great influence upon the anomalously dispersive region of the dispersion relation around the plasma frequency.

Yu, W.↗