Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase field simulation”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Double layers on auroral field lines

Time-stationary solutions to the Vlasov-Poisson equation for ion holes and double layers were examined along with particle simulations which pertain to recent observations of small amplitude (e phi)/t sub e approx. 1 electric field structures on auroral field lines. Both the time-stationary analysis and the simulations suggest that double layers evolve from holes in ion phase space when their amplitude reaches (e phi)/t sub e approx. 1. Multiple small amplitude double layers which are seen in long simulation systems and are seen to propagate past spacecraft may account for the acceleration of plasma sheet electrons to produce the discrete aurora.

Hudson, M. K.↗
NASA NTRS1984

Numerical simulation of reconnection in an emerging magnetic flux region

Numerical solutions in two dimensions are presented for the resistive MHD equations of an initial boundary value problem, simulating reconnection between an emerging magnetic flux region and an overlying coronal magnetic field. This numerical model displays four main phases, which are interpreted in terms of (1) a slowly evolving quasi-steady phase during which most of the magnetic flux emerges, with approximate equilibrium between magnetic and pressure forces; (2) an impulsive phase in which either the reconnection or continuing emergence of the first phase comes to disrupt the equilibrium, and extensive acceleration occurs as the high pressure region within the emerging region drives the fluid upwards and outwards; (3) a second quasi-steady phase; and (4) a potential-static phase in which continuing reconnection, ohmic dissipation, and fluid transport through boundaries depletes the system of all currents and flows.

Forbes, T. G.↗
NASA NTRS1995

A Large-Particle Monte Carlo Code for Simulating Non-Linear High-Energy Processes Near Compact Objects

High-energy radiation processes in compact cosmic objects are often expected to have a strongly non-linear behavior. Such behavior is shown, for example, by electron-positron pair cascades and the time evolution of relativistic proton distributions in dense radiation fields. Three independent techniques have been developed to simulate these non-linear problems: the kinetic equation approach; the phase-space density (PSD) Monte Carlo method; and the large-particle (LP) Monte Carlo method. In this paper, we present the latest version of the LP method and compare it with the other methods. The efficiency of the method in treating geometrically complex problems is illustrated by showing results of simulations of 1D, 2D and 3D systems. The method is shown to be powerful enough to treat non-spherical geometries, including such effects as bulk motion of the background plasma, reflection of radiation from cold matter, and anisotropic distributions of radiating particles. It can therefore be applied to simulate high-energy processes in such astrophysical systems as accretion discs with coronae, relativistic jets, pulsar magnetospheres and gamma-ray bursts.

Stern, Boris E.↗
NASA NTRS2007

The NIRCam Optical Telescope Simulator (NOTES)

The Near Infra-Red Camera (NIRCam), the 0.6-5.0 micron imager and wavefront sensing instrument for the James Webb Space Telescope (JWST), will be used on orbit both as a science instrument, and to tune the alignment of the telescope. The NIRCam Optical Telescope Element Simulator (NOTES) will be used during ground testing to provide an external stimulus to verify wavefront error, imaging characteristics, and wavefront sensing performance of this crucial instrument. NOTES is being designed and built by NASA Goddard Space Flight Center with the help of Swales Aerospace and Orbital Sciences Corporation. It is a single-point imaging system that uses an elliptical mirror to form an U20 image of a point source. The point source will be fed via optical fibers from outside the vacuum chamber. A tip/tilt mirror is used to change the chief ray angle of the beam as it passes through the aperture stop and thus steer the image over NIRCam's field of view without moving the pupil or introducing field aberrations. Interchangeable aperture stop elements allow us to simulate perfect JWST wavefronts for wavefront error testing, or introduce transmissive phase plates to simulate a misaligned JWST segmented mirror for wavefront sensing verification. NOTES will be maintained at an operating temperature of 80K during testing using thermal switches, allowing it to operate within the same test chamber as the NIRCam instrument. We discuss NOTES' current design status and on-going development activities.

Kubalak, David↗
DOE OSTI2022

Modeling impacts of ice-nucleating particles from marine aerosols on mixed-phase orographic clouds during 2015 ACAPEX field campaign

Abstract. A large fraction of annual precipitation over the western United States comes from wintertime orographic clouds associated with atmospheric rivers (ARs). Transported African and Asian dust and marine aerosols from the Pacific Ocean may act as ice-nucleating particles (INPs) to affect cloud and precipitation properties over the region. Here we explored the effects of INPs from marine aerosols on orographic mixed-phase clouds and precipitation at different AR stages for an AR event observed during the 2015 ACAPEX field campaign under low dust (<0.02 cm−3) conditions. Simulations were conducted using the chemistry version of the Weather Research and Forecasting Model coupled with the spectral-bin microphysics at 1 km grid spacing, with ice nucleation connected with dust and marine aerosols. By comparing against airborne and ground-based observations, accounting for marine INP effects improves the simulation of AR-precipitation. The marine INPs enhance the formation of ice and snow, leading to less shallow warm clouds but more mixed-phase and deep clouds, as well as to a large spillover effect of precipitation after AR landfall. The responses of cloud and precipitation to marine INPs vary with the AR stages, with more significant effects before AR landfall and post-AR than after AR landfall, mainly because the moisture and temperature conditions change with the AR evolution. This work suggests weather and climate models need to consider the impacts of marine INPs since their contribution is notable under low dust conditions despite the much lower relative ice nucleation efficiency of marine INPs.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1996

Numerical Investigation of Slag Behavior for RSRM

It is known that the flow field of the redesigned solid rocket motor (RSRM) is very complicated due to the complex characteristics of turbulent multi-phase flow, chemical reaction, particle combustion, evaporation, breakup and agglomeration etc. It requires multi-phase calculations, chemical reaction simulation, and particle combustion, evaporation, and breakup models to obtain a better understanding of thermophysics for the RSRM design using numerical methods. Also, the slag buildup due to the molten particles is another factor affecting the performance of the RSRM. To achieve this goal, the volume of fluid (VOF) method is used to capture the free surface motion so as to simulate the accumulation of the molten particles (slag) of the RSRM. A finite rate chemistry model is used to simulate the chemical reaction effects. For multi-phase calculations, the Hermsen combustion model is used for the aluminum particle combustion analysis and the Taylor Analogy Breakup (TAB) model is used for the particle breakup analysis. An interphase mas-exchange model introduced by Spalding is used for the evaporation calculation. The particle trajectories are calculated using a one-step implicit method for several groups of particle sizes by which the drag forces and heat fluxes are then coupled with the gas phase equations. The preliminary results predicted a reasonable physical simulation of the particle effects using a simple two dimensional solid rocket motor configuration. It shows that the AL/AL2O3 particle sizes are reduced due to the combustion, evaporation, and breakup. The flow field is disturbed by the particles. Mach number distributions in the nozzle are deformed due to the effect of particle concentrations away from the center line.

Liaw, P.↗
DOE OSTI2020

Dehydration entropy drives liquid-liquid phase separation by molecular crowding

Complex coacervation driven liquid-liquid phase separation (LLPS) of biopolymers has been attracting attention as a novel phase in living cells. Studies of LLPS in this context are typically of proteins harboring chemical and structural complexity, leaving unclear which properties are fundamental to complex coacervation versus protein-specific. This study focuses on the role of polyethylene glycol (PEG)—a widely used molecular crowder—in LLPS. Significantly, entropy-driven LLPS is recapitulated with charged polymers lacking hydrophobicity and sequence complexity, and its propensity dramatically enhanced by PEG. Experimental and field-theoretic simulation results are consistent with PEG driving LLPS by dehydration of polymers, and show that PEG exerts its effect without partitioning into the dense coacervate phase. Furthermore, it is then up to biology to impose additional variations of functional significance to the LLPS of biological systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

Patterns in simulated turbulent channel flow

The surface stress is a key diagnostic in wall bounded turbulent flows. Large fluctuations in the stress are believed to be associated with intermittent 'bursting' events during which a large proportion of the turbulence production takes place. If this is so, then a detailed investigation of the structure of the surface stress and its spatial relationship to events within the flow could have wide application in drag reduction and other aspects of flow control. The initial phase of this project, therefore, concentrated on the surface stress field. The first objective is to carry out a statistical analysis of the instantaneous surface stress in a simulated turbulent channel flow, including comparison with multipoint experimental data from a zero pressure gradient turbulent boundary layer. The second objective is to apply a simple pattern recognition procedure to educe the characteristic spatial structure of various flow diagnostics. The final objective is to extend the pattern recognition analysis to examine the whole three dimensional structure of the flow.

Stretch, D.↗
NASA NTRS1992

Parallel arrays of Josephson junctions for submillimeter local oscillators

In this paper we discuss the influence of the DC biasing circuit on operation of parallel biased quasioptical Josephson junction oscillator arrays. Because of nonuniform distribution of the DC biasing current along the length of the bias lines, there is a nonuniform distribution of magnetic flux in superconducting loops connecting every two junctions of the array. These DC self-field effects determine the state of the array. We present analysis and time-domain numerical simulations of these states for four biasing configurations. We find conditions for the in-phase states with maximum power output. We compare arrays with small and large inductances and determine the low inductance limit for nearly-in-phase array operation. We show how arrays can be steered in H-plane using the externally applied DC magnetic field.

Pance, Aleksandar↗
DOE OSTI2023

Microstructural Pattern Formation during Far-from-Equilibrium Alloy Solidification

We introduce a new phase-field formulation of rapid alloy solidification that quantitatively in- corporates non-equilibrium effects at the solid-liquid interface over a very wide range of interface velocities. Simulations identify a new dynamical instability of dendrite tip growth driven by solute trapping at velocities approaching the absolute stability limit. They also reproduce the formation of the widely observed banded microstructures, revealing how this instability triggers transitions between dendritic and microsegregation-free solidification. Predicted band spacings agree quantita- tively with observations in rapidly solidified Al-Cu thin films.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Unitary quantum lattice simulations for Maxwell equations in vacuum and in dielectric media

Utilizing the similarity between the spinor representation of the Dirac and the Maxwell equations that has been recognized since the early days of relativistic quantum mechanics, a quantum lattice algorithm (QLA) representation of unitary collision-stream operators of Maxwell's equations is derived for both homogeneous and inhomogeneous media. A second-order accurate 4-spinor scheme is developed and tested successfully for two-dimensional (2-D) propagation of a Gaussian pulse in a uniform medium whereas for normal (1-D) incidence of an electromagnetic Gaussian wave packet onto a dielectric interface requires 8-component spinors because of the coupling between the two electromagnetic polarizations. In particular, the well-known phase change, field amplitudes and profile widths are recovered by the QLA asymptotic profiles without the imposition of electromagnetic boundary conditions at the interface. The QLA simulations yield the time-dependent electromagnetic fields as the wave packet enters and straddles the dielectric boundary. QLA involves unitary interleaved non-commuting collision and streaming operators that can be coded onto a quantum computer: the non-commutation being the very reason why one perturbatively recovers the Maxwell equations.

Physics↗
DOE OSTI2021

A Review of Lattice-Boltzmann Models Coupled with Geochemical Modeling Applied for Simulation of Advanced Waterflooding and Enhanced Oil Recovery Processes

To maintain economic profit and improve the oil production efficiency after the primary and secondary production phase, advanced waterflooding techniques such as low salinity waterflooding in carbonate reservoirs have been investigated in numerical simulations, laboratory experiments, and field pilot tests. Multiple underlying mechanisms have been proposed based on these studies, and they are still under debate. Various numerical modeling approaches are introduced, but there exists a lack of a pore-scale comprehensive modeling scheme to fully understand the processes. Lattice-Boltzmann method (LBM) is a type of numerical fluid flow modeling technique that shows capabilities and flexibilities in modeling pore-scale fluid flow to integrate physical–chemical processes within complex structures. The intrinsic feature of LBM makes it a promising framework for simulating advanced waterflooding due to its flexibility, accuracy, and parallel efficiency. LBM works either by itself for solving reactive transport problems or by coupling with a third-party reaction solver. This review mainly introduces the LBM fluid flow and reactive transport capabilities and the concept and modeling approaches to simulate advanced waterflooding techniques. Meanwhile, an evaluation of the coupled LBM models for enhanced oil recovery (EOR) simulations is discussed with future research challenges and directions concluded.

02 PETROLEUM↗
DOE OSTI2021

A Hotspot’s Better Half: Non-Equilibrium Intra-Molecular Strain in Shock Physics

Shockwave interactions with a material’s microstructure localizes energy into hotspots, which act as nucleation sites for complex processes such as phase transformations and chemical reactions. To date, hotspots have been described via their temperature fields. Nonreactive, all-atom molecular dynamics simulations of shock-induced pore collapse in a molecular crystal show that more energy is localized as potential energy (PE) than can be inferred from the temperature field and that PE localization persists beyond thermal diffusion. Furthermore, the origin of the PE hotspot is traced to large intramolecular strains, storing energy in modes readily available for chemical decomposition.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Generation of large-volume high-pressure plasma by spatiotemporal control of space charge

Any attempt to scale pressure and volume of nonthermal plasma usually leads to instabilities due to the formation of localized space charge. The control of the plasma is limited by the discharge geometry, type of excitation, and gas composition. This article explores the possibility of controlling the space charge in a discharge with a spatially and temporally varying electric field. It is shown that a phase-staggered sinusoidal excitation to a set of conformal azimuthal electrodes in a cylindrical geometry leads to a traveling electric field. Simulations show that in space charge dominated transport, the charged species are dispersed both in the radial and azimuthal directions. This will lead to better control of the space charge and stable discharges near atmospheric pressures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Solitary waves and double layers on auroral field lines

Time stationary solutions to the Vlasov-Poisson equations for ion holes and double layers are examined along with particle simulations that pertain to recent observations of small amplitude electric field structures on auroral field lines. Both the time stationary analysis and the simulations suggest that the observed double layers evolve from holes in ion phase space. Multiple small amplitude double layers, as seen in long simulation systems, are observed to propagate past the spacecraft and may account for the acceleration of plasma sheet electrons to produce inverted-V precipitation.

Hudson, M. K.↗
NASA NTRS1997

Comparison of Orbiter PRCS Plume Flow Fields Using CFD and Modified Source Flow Codes

The Space Shuttle Orbiter will use Reaction Control System (RCS) jets for docking with the planned International Space Station (ISS). During approach and backout maneuvers, plumes from these jets could cause high pressure, heating, and thermal loads on ISS components. The object of this paper is to present comparisons of RCS plume flow fields used to calculate these ISS environments. Because of the complexities of 3-D plumes with variable scarf-angle and multi-jet combinations, NASA/JSC developed a plume flow-field methodology for all of these Orbiter jets. The RCS Plume Model (RPM), which includes effects of scarfed nozzles and dual jets, was developed as a modified source-flow engineering tool to rapidly generate plume properties and impingement environments on ISS components. This paper presents flow-field properties from four PRCS jets: F3U low scarf-angle single jet, F3F high scarf-angle single jet, DTU zero scarf-angle dual jet, and F1F/F2F high scarf-angle dual jet. The RPM results compared well with plume flow fields using four CFD programs: General Aerodynamic Simulation Program (GASP), Cartesian (CART), Unified Solution Algorithm (USA), and Reacting and Multi-phase Program (RAMP). Good comparisons of predicted pressures are shown with STS 64 Shuttle Plume Impingement Flight Experiment (SPIFEX) data.

Rochelle, Wm. C.↗
NASA NTRS2000

Electric Plasma Arc-Lamp Combustor Liner Durability Test System Developed

Silicon carbide matrix composites are candidate materials for high-temperature combustor liners. Because through-the-thickness thermal gradients are the primary cause of stress on combustor liners, a unique test facility was developed at the NASA Glenn Research Center at Lewis Field to simulate in-service pure thermal stress distributions in fiber reinforced silicon carbide cylinders. It was developed initially under Phase II of the High- Speed Research Program. This test stand can accommodate 8-in.-long test cylinders that have outer diameters of 4 in. and a wall thickness of about 0.08 to 0.12 in. One cylinder at a time is loaded vertically into the test stand. Water-cooled plates enclose the open ends of the cylinder and provide cooling. Load plates on the exterior side of the water-cooled plates provide support and compression loads. To evaluate a combustor liner material s potential performance, researchers induce thermal gradients with an axisymmetric, direct-current, electric arc within the cylinder while refrigerated air at a rate of 1.5 lb/sec impinges on the outside surface of the cylinder. The achievable through-the-thickness thermal gradient is predicted to be in excess of 200 C. The 8-in. long, 0.5-in.-diameter plasma arc emits full spectrum visible light; radiant intensity exceeds 300 Watts per square centimeters to produce temperatures in excess of 1500 C on materials with emissivity near unity. Because the system does not rely upon the combustion of fuels to achieve the related thermal conditions, ancillary environmental reactions with the sample are eliminated. The system incorporates a standard mechanical test frame, which can impose constant as well cyclical axial stresses up to 2200 lb upon the test piece. Silicon-carbide-fiber reinforced silicon carbide matrix composite cylinders were instrumented with thin-film thermocouples to obtain through-the-thickness thermal flux measurements. Inside wall temperatures reached 1200 C with only 250 A of current. One of the special features of this configuration is the creation of hoop stress states within the cylinder, which up this point have not been obtainable in planar coupon tests. This facility will allow various operational modes, including accelerated tests of thermal transients simulating the effects of repeated engine ignition as well as prescribed thermal and mechanical histories to simulate various duty cycle profiles. Tests can now be performed on thermal-barrier-coated metallic liners and ceramic composite liners that require a combination of high heat flux and controlled mechanical stresses.

Goldsby, Jon C.↗
NASA NTRS2016

Simulation of Satellite, Airborne and Terrestrial LiDAR with DART (I):Waveform Simulation with Quasi-Monte Carlo Ray Tracing

Light Detection And Ranging (LiDAR) provides unique data on the 3-D structure of atmosphere constituents and the Earth's surface. Simulating LiDAR returns for different laser technologies and Earth scenes is fundamental for evaluating and interpreting signal and noise in LiDAR data. Different types of models are capable of simulating LiDAR waveforms of Earth surfaces. Semi-empirical and geometric models can be imprecise because they rely on simplified simulations of Earth surfaces and light interaction mechanisms. On the other hand, Monte Carlo ray tracing (MCRT) models are potentially accurate but require long computational time. Here, we present a new LiDAR waveform simulation tool that is based on the introduction of a quasi-Monte Carlo ray tracing approach in the Discrete Anisotropic Radiative Transfer (DART) model. Two new approaches, the so-called "box method" and "Ray Carlo method", are implemented to provide robust and accurate simulations of LiDAR waveforms for any landscape, atmosphere and LiDAR sensor configuration (view direction, footprint size, pulse characteristics, etc.). The box method accelerates the selection of the scattering direction of a photon in the presence of scatterers with non-invertible phase function. The Ray Carlo method brings traditional ray-tracking into MCRT simulation, which makes computational time independent of LiDAR field of view (FOV) and reception solid angle. Both methods are fast enough for simulating multi-pulse acquisition. Sensitivity studies with various landscapes and atmosphere constituents are presented, and the simulated LiDAR signals compare favorably with their associated reflectance images and Laser Vegetation Imaging Sensor (LVIS) waveforms. The LiDAR module is fully integrated into DART, enabling more detailed simulations of LiDAR sensitivity to specific scene elements (e.g., atmospheric aerosols, leaf area, branches, or topography) and sensor configuration for airborne or satellite LiDAR sensors.

Gastellu-Etchegorry, Jean-Philippe↗