Search NASA⌕ Search

SEARCH · Search NASA

Results for “Phase field simulations”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Polarization‐Engineered Near‐Field Generation Using a Hybrid Tip–Antenna System

Precise control of light polarization at the nanoscale is critical for accessing chiral optical responses and manipulating spin–photon interactions in advanced materials. Yet, conventional scattering-type near-field probes predominantly generate out-of-plane linear polarization and offer little control over phase or polarization state. Here, in this study, we introduce a polarization-engineered near-field methodology based on a combined metallic tip and planar dipole nanoantenna system. Using full-wave electromagnetic simulations, we show that the tip acts as a vertically oriented plasmonic resonator, while the antenna supports an in-plane dipolar mode. By tuning the tip–antenna geometry and tip height, the two orthogonal field components attain comparable amplitudes and a controllable ∼90° phase offset, producing circularly polarized nano-light in the antenna gap. The proposed system effectively functions as a nanoscale quarter-wave plate, converting linearly polarized illumination into circularly polarized hotspots without external polarization optics. This method establishes an experimentally accessible route toward polarization-programmable near-field nanoscopy, enabling chiral spectroscopy, selective excitation of spin/valley degrees of freedom, and quantum optical investigations at the nanoscale.

36 MATERIALS SCIENCE↗
NASA NTRS1998

Modeling Macro- and Micro-Scale Turbulent Mixing and Chemistry in Engine Exhaust Plumes

Simulation of turbulent mixing and chemical processes in the near-field plume and plume-vortex regimes has been successfully carried out recently using a reduced gas phase kinetics mechanism which substantially decreased the computational cost. A detailed mechanism including gas phase HOx, NOx, and SOx chemistry between the aircraft exhaust and the ambient air in near-field aircraft plumes is compiled. A reduced mechanism capturing the major chemical pathways is developed. Predictions by the reduced mechanism are found to be in good agreement with those by the detailed mechanism. With the reduced chemistry, the computer CPU time is saved by a factor of more than 3.5 for the near-field plume modeling. Distributions of major chemical species are obtained and analyzed. The computed sensitivities of major species with respect to reaction step are deduced for identification of the dominant gas phase kinetic reaction pathways in the jet plume. Both the near field plume and the plume-vortex regimes were investigated using advanced mixing models. In the near field, a stand-alone mixing model was used to investigate the impact of turbulent mixing on the micro- and macro-scale mixing processes using a reduced reaction kinetics model. The plume-vortex regime was simulated using a large-eddy simulation model. Vortex plume behind Boeing 737 and 747 aircraft was simulated along with relevant kinetics. Many features of the computed flow field show reasonable agreement with data. The entrainment of the engine plumes into the wing tip vortices and also the partial detrainment of the plume were numerically captured. The impact of fluid mechanics on the chemical processes was also studied. Results show that there are significant differences between spatial and temporal simulations especially in the predicted SO3 concentrations. This has important implications for the prediction of sulfuric acid aerosols in the wake and may partly explain the discrepancy between past numerical studies (that employed parabolic or temporal approximations) and the measured data. Finally to address the major uncertainty in the near-field plume modeling related to the plume processing of sulfur compounds and advanced model was developed to evaluate its impact on the chemical processes in the near wake. A comprehensive aerosol model is developed and it is coupled with chemical kinetics and the axisymmetric turbulent jet flow models. The integrated model is used to simulate microphysical processes in the near-field jet plume, including sulfuric acid and water binary homogeneous nucleation, coagulation, non-equilibrium heteromolecular condensation, and sulfur-induced soot activation. The formation and evolution of aerosols are computed and analyzed. The computed results show that a large number of ultra-fine (0.3--0.6 nm in radius) volatile HSO4 - HO embryos are generated in the near-field plume. These embryos further grow in size by self coagulation and condensation. Soot particles can be activated by both heterogeneous nucleation and scavenging of H2SO4-H2O aerosols. These activated soot particles can serve as water condensation nuclei for contrail formation. Conditions under which ice contrails can form behind aircrafts are studied. The sensitivities of the threshold temperature for contrail formation with respect to aircraft propulsion efficiency, relative humidity, and ambient pressure are evaluated. The computed aerosol properties for different extent of fuel sulfur conversion to S(VI) (SO3 and H2SO4) in engine are examined and the results are found to be sensitive to this conversion fraction.

Menon, Suresh↗
NASA NTRS2013

North Atlantic Simulations in Coordinated Ocean-Ice Reference Experiments Phase II (CORE-II). Part I: Mean States

Simulation characteristics from eighteen global ocean-sea-ice coupled models are presented with a focus on the mean Atlantic meridional overturning circulation (AMOC) and other related fields in the North Atlantic. These experiments use inter-annually varying atmospheric forcing data sets for the 60-year period from 1948 to 2007 and are performed as contributions to the second phase of the Coordinated Oceanice Reference Experiments (CORE-II). The protocol for conducting such CORE-II experiments is summarized. Despite using the same atmospheric forcing, the solutions show significant differences. As most models also differ from available observations, biases in the Labrador Sea region in upper-ocean potential temperature and salinity distributions, mixed layer depths, and sea-ice cover are identified as contributors to differences in AMOC. These differences in the solutions do not suggest an obvious grouping of the models based on their ocean model lineage, their vertical coordinate representations, or surface salinity restoring strengths. Thus, the solution differences among the models are attributed primarily to use of different subgrid scale parameterizations and parameter choices as well as to differences in vertical and horizontal grid resolutions in the ocean models. Use of a wide variety of sea-ice models with diverse snow and sea-ice albedo treatments also contributes to these differences. Based on the diagnostics considered, the majority of the models appear suitable for use in studies involving the North Atlantic, but some models require dedicated development effort. atmospheric forcing atmospheric temperatures

ocean model comparisons↗
DOE OSTI2026

Interface and Thermophysical Properties of R 32 Refrigerant

Driven by the urgent demand for efficient cooling in microelectronics and advanced thermal management systems, difluoromethane (R32/CH 2 F 2 ) has emerged as a promising candidate owing to its favorable thermophysical properties, including high heat transfer efficiency and low viscosity. While bulk properties such as density, viscosity, and thermal conductivity have been widely studied, interfacial properties, including surface tension and interfacial thickness, remain comparatively underexplored, despite their importance in phase-transition dynamics. Here, we perform molecular dynamics (MD) simulations from 180 to 300 K using an optimized transferable force field for fluoropropenes with enhanced electrostatics to assess both bulk and interfacial behavior of R32. Simulations reproduced density within ±2.1%, viscosity within 3.05%, and thermal conductivity within 7.41% of NIST reference data. Heat capacities (C p and C v ) were predicted within 5%. For interfacial properties, surface tension trends were reproduced within 13.58% deviation, and the vapor–liquid coexistence curve closely matched reference data, yielding a critical temperature of 345.7 K (1.6% deviation) and a critical density of 0.397 g/cm 3 (6.4% deviation). Importantly, the vapor–liquid interface exhibited pronounced temperature-dependent broadening across the 180–290 K range. This behavior correlates with increasing molecular kinetic energy, reduction in intermolecular cohesive interactions, and a progressive loss of preferential dipole alignment, which collectively enhance thermal fluctuation amplitudes at elevated temperatures. These validated results provide predictive molecular-level insights, particularly for interfacial properties that remain less characterized. By reducing property prediction errors in key parameters such as critical temperature, this work provides reliable inputs for heat-exchanger and system models. Such correlations can support optimized component sizing, improved performance, and reduced refrigerant charge. Beyond R32, the methodology offers a transferable framework for blended and next-generation low-GWP refrigerants, contributing to sustainable thermal management aligned with the 2027 EU F-Gas regulation and 2030 Kigali Amendment.

Fluids↗
DOE OSTI2026

Feynman diagrams for matter wave interferometry

We introduce a new theoretical framework based on Feynman diagrams to compute phase shifts in matter wave interferometry. The method allows for analytic computation of higher order quantum corrections, beyond the traditional semi-classical approximation. These additional terms depend on the finite size of the initial matter wavefunction and/or have higher order dependence on ℏ. We apply the method to compute the response of matter wave interferometers to power law potentials and potentials with an arbitrary spatial dependence. The analytic expressions are validated by comparing to numerical simulations, and estimates are provided for the scale of the quantum corrections to the phase shift response to the gravitational field of the earth, anharmonic trapping potentials, and gravitational fields from local proof masses. We also find that for certain experimentally feasible parameters, these corrections are large enough to be measured and could lead to systematic errors if they are not mitigated. We find that to first order in a spatially dependent potential, quantum corrections vanish when the initial matter wavepacket has spherical symmetry and the potential satisfies Laplace's equation. We anticipate these quantum corrections will be especially important for trapped matter wave interferometers and for free-space matter wave interferometers in the presence of proof masses. These interferometers are becoming increasingly sensitive tools for mobile inertial sensing, gravity surveying, tests of gravity and its interplay with quantum mechanics, and searches for dark energy.

Glick, Jonah [Northwestern U.; Fermilab] (ORCID:00↗
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↗
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.↗
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↗
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↗
NASA NTRS2017

Data Assimilation and Uncertainties in Early Solar Cycle Predictions

Stochastic nature of solar activity variations together with our limited knowledge of the dynamo mechanism and subsurface dynamics causes uncertainty in predictions of the solar cycle. For improving the physics-­‐based predictions we can take advantage of the mathematical data assimilation approach that allows us to take into account both, observational errors and model uncertainties, and provide estimates of the next solar cycle along with prediction uncertainties. In this study we use the Parker's migratory dynamo model together with the equation of magnetic helicity balance, which reproduces main properties of the sunspot cycles and allow us to minimize discrepancies between the observed global activity variations and the model solution. The test simulation runs show that a reliable prediction can be obtained for two phases of preceding solar cycle: 1) if the polar field reversals shortly after the solar maxima (strong toroidal field and weak poloidal field), and 2) during the solar minima (strongest poloidal and weak toroidal fields). The early estimate of Cycle 25 obtained by this method shows that this cycle will start in 2019 – 2020, reach the maximum in 2023 – 2024, and that the mean sunspot number at the maximum will be about 90 (for the v2.0 sunspot number series).

Data Assimilation↗
NASA NTRS2021

A Field-Particle Correlation Analysis of a Perpendicular Magnetized Collisionless Shock

Using the field-particle correlation technique, we examine the particle energization in a 1D-2V continuum Vlasov-Maxwell simulation of a perpendicular magnetized collisionless shock. The combination of the field-particle correlation technique with the high fidelity representation of the particle distribution function provided by a direct discretization of the Vlasov equation allows us to ascertain the details of the exchange of energy between the electromagnetic fields and the particles in phase space. We identify the velocity-space signatures of shock-drift acceleration of the ions and adiabatic heating of the electrons due to the perpendicular collisionless shock by constructing a simplified model with the minimum ingredients necessary to produce the observed energization signatures in the self-consistent Vlasov-Maxwell simulation. We are thus able to completely characterize the energy transfer in the perpendicular collisionless shock considered here and provide predictions for the application of the field-particle correlation technique to spacecraft measurements of collisionless shocks.

James Juno↗