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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Low-Frequency Oscillations and Transport Processes Induced by Multiscale Transverse Structures in the Polar Wind Outflow: A Three-Dimensional Simulation

Multiscale transverse structures in the magnetic-field-aligned flows have been frequently observed in the auroral region by FAST and Freja satellites. A number of multiscale processes, such as broadband low-frequency oscillations and various cross-field transport effects are well correlated with these structures. To study these effects, we have used our three-dimensional multifluid model with multiscale transverse inhomogeneities in the initial velocity profile. Self-consistent-frequency mode driven by local transverse gradients in the generation of the low field-aligned ion flow and associated transport processes were simulated. Effects of particle interaction with the self-consistent time-dependent three-dimensional wave potential have been modeled using a distribution of test particles. For typical polar wind conditions it has been found that even large-scale (approximately 50 - 100 km) transverse inhomogeneities in the flow can generate low-frequency oscillations that lead to significant flow modifications, cross-field particle diffusion, and other transport effects. It has also been shown that even small-amplitude (approximately 10 - 20%) short-scale (approximately 10 km) modulations of the original large-scale flow profile significantly increases low-frequency mode generation and associated cross-field transport, not only at the local spatial scales imposed by the modulations but also on global scales. Note that this wave-induced cross-field transport is not included in any of the global numerical models of the ionosphere, ionosphere-thermosphere, or ionosphere-polar wind. The simulation results indicate that the wave-induced cross-field transport not only affects the ion outflow rates but also leads to a significant broadening of particle phase-space distribution and transverse particle diffusion.

Ganguli, Supriya B.↗
NASA NTRS2012

Numerical Simulation of Rocket Exhaust Interaction with Lunar Soil

This technology development originated from the need to assess the debris threat resulting from soil material erosion induced by landing spacecraft rocket plume impingement on extraterrestrial planetary surfaces. The impact of soil debris was observed to be highly detrimental during NASA s Apollo lunar missions and will pose a threat for any future landings on the Moon, Mars, and other exploration targets. The innovation developed under this program provides a simulation tool that combines modeling of the diverse disciplines of rocket plume impingement gas dynamics, granular soil material liberation, and soil debris particle kinetics into one unified simulation system. The Unified Flow Solver (UFS) developed by CFDRC enabled the efficient, seamless simulation of mixed continuum and rarefied rocket plume flow utilizing a novel direct numerical simulation technique of the Boltzmann gas dynamics equation. The characteristics of the soil granular material response and modeling of the erosion and liberation processes were enabled through novel first principle-based granular mechanics models developed by the University of Florida specifically for the highly irregularly shaped and cohesive lunar regolith material. These tools were integrated into a unique simulation system that accounts for all relevant physics aspects: (1) Modeling of spacecraft rocket plume impingement flow under lunar vacuum environment resulting in a mixed continuum and rarefied flow; (2) Modeling of lunar soil characteristics to capture soil-specific effects of particle size and shape composition, soil layer cohesion and granular flow physics; and (3) Accurate tracking of soil-borne debris particles beginning with aerodynamically driven motion inside the plume to purely ballistic motion in lunar far field conditions. In the earlier project phase of this innovation, the capabilities of the UFS for mixed continuum and rarefied flow situations were validated and demonstrated for lunar lander rocket plume flow impingement under lunar vacuum conditions. Applications and improvements to the granular flow simulation tools contributed by the University of Florida were tested against Earth environment experimental results. Requirements for developing, validating, and demonstrating this solution environment were clearly identified, and an effective second phase execution plan was devised. In this phase, the physics models were refined and fully integrated into a production-oriented simulation tool set. Three-dimensional simulations of Apollo Lunar Excursion Module (LEM) and Altair landers (including full-scale lander geometry) established the practical applicability of the UFS simulation approach and its advanced performance level for large-scale realistic problems.

Liever, Peter↗
NASA NTRS2000

Adaptive Detector Arrays for Optical Communications Receivers

The structure of an optimal adaptive array receiver for ground-based optical communications is described and its performance investigated. Kolmogorov phase screen simulations are used to model the sample functions of the focal-plane signal distribution due to turbulence and to generate realistic spatial distributions of the received optical field. This novel array detector concept reduces interference from background radiation by effectively assigning higher confidence levels at each instant of time to those detector elements that contain significant signal energy and suppressing those that do not. A simpler suboptimum structure that replaces the continuous weighting function of the optimal receiver by a hard decision on the selection of the signal detector elements also is described and evaluated. Approximations and bounds to the error probability are derived and compared with the exact calculations and receiver simulation results. It is shown that, for photon-counting receivers observing Poisson-distributed signals, performance improvements of approximately 5 dB can be obtained over conventional single-detector photon-counting receivers, when operating in high background environments.

Vilnrotter, V.↗
NASA NTRS2004

Limitation of Liquid Crystal on Silicon Spatial Light Modular for Holographic Three-dimensional Displays

In a 3-D display system based on an opto-electronic reconstruction of a digitally recorded hologram, the field of view of such a system is limited by the spatial resolution of the liquid crystal on silicon (LCOS) spatial light modular (SLM) used to perform the opto-electronic reconstruction. In this article, the special resolution limitation of LCOS SLM associated with the fringe field effect and interpixel coupling is determined by the liquid crystal detector simulation and the Finite Difference Time Domain (FDTD) simulation. The diffraction efficiency loss associated with the imperfection in the phase profile is studied with an example of opto-electronic reconstruction of an amplitude object. A high spatial resolution LCOS SLM with a wide reconstruction angle is proposed.

Wang, Xinghua↗
NASA NTRS2023

Understanding the Deflection of the “Cartwheel CME”: Data Analysis and Modeling

We study the low corona evolution of the “Cartwheel” coronal mass ejection (CME; 2008 April 9) by- reconstructing its three-dimensional path and modeling it with magnetohydrodynamic simulations. This event exhibited a double deflection that has been reported and analyzed in previous works but whose underlying cause remained unclear. The Cartwheel CME traveled toward a coronal hole (CH) and against the magnetic gradients. Using a high-cadence, full-trajectory reconstruction, we accurately determine the location of the magnetic flux rope (MFR) and, consequently, the magnetic environment in which it is immersed. We find a pseudostreamer (PS) structure whose null point may be responsible for the complex evolution of the MFR at the initial phase. From the preeruptive magnetic field reconstruction, we estimate the dynamic forces acting on the MFR and provide a new physical insight into the motion exhibited by the 2008 April 9 event. By setting up a similar magnetic configuration in a 2.5D numerical simulation we are able to reproduce the observed behavior, confirming the importance of the PS null point. We find that the magnetic forces directed toward the null point cause the first deflection, directing the MFR toward the CH. Later, the magnetic pressure gradient of the CH produces the reversal motion of the MFR.

Solar prominences↗
NASA NTRS1998

Can Kelvin-Helmholtz Instabilities of Jet-Like Structures and Plumes Cause Solar Wind Fluctuations at 1AU

During its south polar passage in 1994, Ulysses sampled the solar wind emerging from the south polar coronal hole. In analysing these data two types of density fluctuations lasting a few hours have been reported, one characterized by fluctuations in velocity gradients ("microstreams") and the other by magnetic plus thermal-pressure balance structures ("PBS"). Microstreams were both temporal and spatial in nature. At higher frequencies, MHD turbulence was observed and found to be less evolved than in the ecliptic, but essentially independent of heliographic latitude. It is argued here that microstreams, PBS, and MHD turbulence could all be the remnant of shears associated with plumes and other filamentary structures ("jets") which have been reported to exist in coronal holes. The shear between a jet and its ambient can become unstable to the MHD Kelvin-Helmholtz ("KH') instability at 5-10 solar radius and the propagating instability can cause fluctuations like those seen by Ulysses. This motivates us to simulate coronal jets using a 3D MHD ZEUS code. The first 2D results have just started to come and are promising. To study the KH instability the jet is perturbed at the boundary with a linear amplitude and fixed frequency. The jet seems to pass through various distinct phases, one of which Is apparently dominated by KH instabilities. These instabilities drive oblique shocks into the jet as the turbulent eddies contact the jet surface. It is known that KH instabilities and internal shock waves are partially suppressed by magnetic field tension. Hence, in simulating far along the jet the Instability is expected to produce Alfvenic fluctuations like those seen near 1 AU.

Parhi, Shyamsundar↗
NASA NTRS2013

Silent and Efficient Supersonic Bi-Directional Flying Wing

The supersonic bi-directional (SBiDir) flying wing (FW) concept has a great potential to achieve low sonic boom with high supersonic aerodynamic performance due to removal of performance conflict between high speed and low speed by rotating goo in flight. This NIAC Phase 1 research has achieved three objectives: 2) prove the concept based on simulation that it can achieve very low boom with smooth Sine wave ground over-pressure signature and excellent aerodynamic efficiency; 3) conduct trade study to correlate the geometric parameters with sonic boom and aerodynamic performance for further automated design optimization in Phase II. The design methodology developed in Phase I includes three parts: 1) an advanced geometry model, which can vary airfoil meanline angle distribution to control the expansion and shock waves on the airplane surface to mitigate sonic boom and improve aerodynamic efficiency. 2) a validated CFD procedure to resolve near field flow with accurate shock strength. The sonic boom propagation from near field to far field ground is simulated by NASA NF Boom code. The surface friction drag prediction is based on fiat plate correlation adopted by Seebass and supported by the experimental study of Winter and Smith, which is on the conservative side and is more reliable than CFD RANS simulation. 3) a mission analysis tool based on Corke's model that provides design requirements and constraints of supersonic airplanes for range, payload, volume, size, weight, etc. The design mission target is a supersonic transport with cruise Mach number 1.6, 100 passengers, and 4000nm range. The trade study has several very important findings: 1) The far field ground sonic boom signature is directly related to the smoothness of the flow on the airplane surface. The meanline angle distribution is a very effective control methodology to mitigate surface shock and expansion wave strength, and mitigating compression wave coalescing by achieving smooth loading distribution chord-wise. Compared with a linear meanline angle distribution, a design using nonlinear and non-monotonic meanline angle distribution is able to reduce the sonic boom ground loudness by over 20dBP1. The design achieves sonic boom ground loudness less than 70dBP1 and aerodynamic dynamic efficiency 1/D of 8.4. 2) Decreasing sweep angle within the Mach cone will increase 1/D as well as sonic boom. A design with variable sweep from 84° at the very leading edge to 68° at the tip achieves an extraordinarily high 1/D of 10.4 at Mach number 1.6 due to the low wave drag. If no sonic boom constraint is attached, SBiDir-FW concept still has a lot of room to increase the 1/D. 3) The round leading edge and trailing edge under high sweep angle are beneficial to improve aerodynamic performance, sonic boom, and to increase volume of the airplane. 4) Subsonic performance is benefited greatly from the high slenderness of supersonic configuration after rotating goo. A design with excellent supersonic aspect ratio of 0.44, 1/D of 8.g, gives an extraordinary subsonic aspect ration of 10 and 1/D of 1g.7. Two configurations are designed in details to install internal seats, landing gears, and engine installation to demonstrate the feasibility of SBiDir-FW configuration to accommodate all the required volume for realistic airplane. Here we emphasize that the qualitative findings in Phase I are very encouraging, more important than the quantitative results. Qualitative findings give the understanding of physics and provide the path to achieve the ultimate high performance design. The promising quantitative results achieved in Phase I need to be confirmed by wind tunnel testing in Phase II and ultimately proved by flight test. The other important step forward will be made to study the rotation transition from both CFD unsteady simulation and wind tunnel testing.

Zha, Gecheng↗
NASA NTRS2001

Two-Dimensional Model Simulations of Interannual Variability in the Tropical Stratosphere

Meteorological data from the United Kingdom Meteorological Office (UKMO) and constituent data from the Upper Atmospheric Research Satellite (UARS) are used to construct yearly zonal mean dynamical fields for the 1990s for use in the GSFC 2-D chemistry and transport model. This allows for interannual dynamical variability to be included in the model constituent simulations. In this study, we focus on the tropical stratosphere. We find that the phase of quasi-biennial oscillation (QBO) signals in equatorial CH4, and profile and total column 03 data is resolved quite well using this empirically- based 2-D model transport framework. However. the QBO amplitudes in the model constituents are systematically underestimated relative to the observations at most levels. This deficiency is probably due in part to the limited vertical resolutions of the 2-D model and the UKMO and UARS input data sets. We find that using different heating rate calculations in the model affects the interannual and QBO amplitudes in the constituent fields, but has little impact on the phase. Sensitivity tests reveal that the QBO in transport dominates the ozone interannual variability in the lower stratosphere. with the effect of the temperature QBO being dominant in the tipper stratosphere via the strong temperature dependence of the ozone loss reaction rates. We also find that the QBO in odd nitrogen radicals, which is caused by the QBO modulated transport of NOy, plays a significant but not dominant role in determining the ozone QBO variability in the middle stratosphere. The model mean age of air is in good overall agreement with that determined from tropical lower,middle stratospheric OMS balloon observations of SF6 and CO2. The interannual variability of tile equatorial mean age in the model increases with altitude and maximizes near 40 km, with a range, of 4-5 years over the 1993-2000 time period.

Fleming, Eric L.↗
NASA NTRS1982

Initial stages of cavitation damage and erosion on copper and brass tested in a rotating disk device

In view of the differences in flow and experimental conditions, there has been a continuing debate as to whether or not the ultrasonic method of producing cavitation damage is similar to the damage occurring in cavitating flow systems, namely, venturi and rotating disk devices. In this paper, the progress of cavitation damage during incubation periods on polycrystalline copper and brass tested in a rotating disk device is presented. The results indicate several similarities and differences in the damage mechanism encountered in a rotating disk device (which simulates field rotary devices) and a magnetostriction apparatus. The macroscopic erosion appears similar to that in the vibratory device except for nonuniform erosion and apparent plastic flow during the initial damage phase.

Rao, P. V.↗
NASA NTRS2017

Prominence Eruption Initiated by Helical Kink Instability of an Embedded Flux Rope

We study the triggering mechanism of a limb-prominence eruption and the associated coronal mass ejection (CME) near AR 12342 using Solar Dynamics Observatory and Large Angle and Spectrometric Coronagraph/Solar Heliospheric Observatory observations. The prominence is seen with an embedded flux thread (FT) at one end and bifurcates from the middle to a different footpoint location. The morphological evolution of the FT is similar to that of an unstable flux rope (FR), which we regard as a prominence-embedded FR. The FR twist exceeds the critical value. In addition, the morphology of the prominence plasma in 304 Å images marks the helical nature of the magnetic skeleton, with a total of 2.96 turns along arc length. The potential field extrapolation model indicates that the critical height of the background magnetic field gradient falls within the inner corona (105 Mm), which is consistent with the extent of coronal plasma loops. These results suggest that the helical kink instability in the embedded FR caused the slow rise of the prominence to the height of the torus instability domain. Moreover, the differential emission measure analysis unveils heating of the prominence plasma to coronal temperatures during an eruption, suggesting reconnection-related heating underneath the upward rising embedded FR. The prominence starts with a slow rise motion of 10 km s(exp. -1), which is followed by fast and slow acceleration phases that have an average acceleration of 28.9 m s(exp. -2) and 2.4 m s(exp. -2) in C2 and C3 field of view, respectively. As predicted by previous numerical simulations, the observed synchronous kinematic profiles of the CME leading edge and the core support the involved FR instability in the prominence initiation.

Vemareddy, P.↗
NASA NTRS2022

Assessment of Two Stochastic Cloud Subcolumn Generators Using Observed Fields of Vertically Resolved Cloud Extinction

We evaluate two stochastic subcolumn generators used in GCMs to emulate subgrid cloud variability enabling comparisons with satellite observations and simulations of certain physical processes. Our evaluation necessitated the creation of a reference observational dataset that resolves horizontal and vertical cloud variability. The dataset combines two CloudSat cloud products that resolve two-dimensional cloud optical depth variability of liquid, ice, and mixed phase clouds when blended at ~200 m vertical and ~ 2 km horizontal scales. Upon segmenting the dataset to individual “scenes”, mean profiles of the cloud fields are passed as input to generators that produce scene-level cloud subgrid variability. The assessment of generator performance at the scale of individual scenes and in a mean sense is largely based on inferred joint histograms that partition cloud fraction within predetermined combinations of cloud top pressure –cloud optical thickness ranges. Our main finding is that both generators tend to underestimate optically thin clouds, while one of them also tends to overestimate some cloud types of moderate and high optical thickness. Associated radiative flux errors are also calculated by applying a simple transformation to the cloud fraction histogram errors, and are found to approach values almost as high as 3 Wm-2 for the cloud radiative effect in the shortwave part of the spectrum.

cloud variability↗
NASA NTRS1993

Propagation characteristics of waves upstream and downstream of quasi-parallel shocks

The propagation characteristics of waves upstream and downstream of quasi-parallel shocks are investigated by using 2D hybrid simulations. At low Alfven Mach numbers, M(A) below about 2, the shock is initially associated with upstream phase-standing whistlers. At later times, backstreaming ions excite longer-wavelength whistlers via the right-hand resonant ion/ion instability. These waves propagate along the magnetic field at a group velocity no smaller than the upstream flow speed, so that the waves remain in the upstream region. At higher MA (above about 3), these waves are convected back into the shock, causing its reformation and downstream perturbations. Shock transmitted waves mode-convert into Alfven/ion-cyclotron waves which have a wave vector along the shock normal (pointing upstream) and convect downstream. The 2D simulation results confirm our earlier suggestion that the upstream waves should be field aligned, and that their convection into the downstream is associated with linear mode conversion into the Alfven/ion-cyclotron branch.

Krauss-Varban, D.↗
NASA NTRS2020

Polarimetric Radar Characteristics of Simulated and Observed Intense Convective Cores for a Midlatitude Continental and Tropical Maritime Environment

This study contrasts midlatitude continental and tropical maritime deep convective cores using polarimetric radar observables and retrievals from selected deep convection episodes during the MC3E and TWPICE field campaigns. The continental convective cores produce stronger radar reflectivities throughout the profiles, while maritime convective cores produce more positive differential reflectivity Z(dr) and larger specific differential phase K(dp) above the melting level. Hydrometeor identification retrievals revealed the presence of large fractions of rimed ice particles (snow aggregates) in the continental (maritime) convective cores, consistent with the Z(dr) and K(dp) observations. The regional cloud-resolving model simulations with bulk and size-resolved bin microphysics are conducted for the selected cases, and the simulation outputs are converted into polarimetric radar signals and retrievals identical to the observational composites. Both the bulk and the bin microphysics reproduce realistic land and ocean (L-O) contrasts in reflectivity, polarimetric variables of rain drops, and hydrometeor profiles, but there are still large uncertainties in describing Zdr and Kdp of ice crystals associated with the ice particle shapes/orientation assumptions. Sensitivity experiments are conducted by swapping background aerosols between the continental and maritime environments, revealing that background aerosols play a role in shaping the distinct L-O contrasts in radar reflectivity associated with raindrop sizes, in addition to the dominant role of background thermodynamics.

Toshihisa Matsui↗
NASA NTRS1995

Nonequilibrium, large-amplitude MHD fluctuations in the solar wind

Compressible MHD simulations in one dimension with three-dimensional vectors are used to investigate a number of processes relevant to problems in interplanetary physics. The simulations indicate that a large-amplitude nonequilibrium (e.g., linearly polarized) Alfvenic wave, which always starts with small relative fluctuations in the magnitude B of the magnetic field, typically evolves to flatten the magnetic profile in most regions. Under a wide variety of conditions B and the density rho become anticorrelated on average. If the mean magnetic field is allowed to decrease in time, the point where the transverse magnetic fluctuation amplitude delta B(sub T) is greater than the mean field B(sub 0) is not special, and large values of delta B(sub T)/B(sub 0) do not cause the compressive thermal energy to increase remarkably or the wave energy to dissipate at an unusually high rate. Nor does the 'backscatter' of the waves that occurs when the sound speed is less than the Alfven speed result, in itself, in substantial energy dissipation, but rather primarily in a phase change between the magnetic and velocity fields. For isolated wave packets the backscatter does not occur for any of the parameters examined; an initial radiation of acoustic waves away from the packet establishes a stable traveling structure. Thus these simulations, although greatly idealized compared to reality, suggest a picture in which the interplanetary fluctuations should have small deltaB and increasingly quasi-pressure balanced compressive fluctuations, as observed, and in which the dissipation and 'saturation' at delta B(sub T)/B(sub 0) approximately = 1 required by some theories of wave acceleration of the solar wind do not occur. The simulations also provide simple ways to understand the processes of nonlinear steepening and backscattering of Alfven waves and demonstrate the existence of previously unreported types of quasi-steady MHD states.

Roberts, D. Aaron↗
NASA NTRS1992

Images constructed from computed flow fields

A method for constructing interferograms, schlieren, and shadowgraphs from ideal- and real-gas, two- and three-dimensional computed flow fields is described. The computational grids can be structured or unstructured, and multiple grids are an option. The constructed images are compared to experimental images for several types of flow, including a ramp, a blunt-body, a nozzle, and a reacting flow. The constructed images simulate the features observed in the experimental images. They are sensitive to errors in the flow-field solutions and can be used to identify solution errors. In addition, techniques for obtaining phase shifts from experimental finite-fringe interferograms and for removing experimentally induced phase-shift errors are discussed. Both the constructed images and calculated phase shifts can be used for validation of computational fluid dynamics (CFD) codes.

Yates, Leslie A.↗
NASA NTRS1994

Magnetotail dynamics under isobaric constraints

Using linear theory and nonlinear MHD simulations, we investigate the resistive and ideal MHD stability of two-dimensional plasma configurations under the isobaric constraint dP/dt = 0, which in ideal MHD is equivalent to conserving the pressure function P = P(A), where A denotes the magnetic flux. This constraint is satisfied for incompressible modes, such as Alfven waves, and for systems undergoing energy losses. The linear stability analysis leads to a Schroedinger equation, which can be investigated by standard quantum mechanics procedures. We present an application to a typical stretched magnetotail configuration. For a one-dimensional sheet equilibrium characteristic properties of tearing instability are rediscovered. However, the maximum growth rate scales with the 1/7 power of the resistivity, which implies much faster growth than for the standard tearing mode (assuming that the resistivity is small). The same basic eigen-mode is found also for weakly two-dimensional equilibria, even in the ideal MHD limit. In this case the growth rate scales with the 1/4 power of the normal magnetic field. The results of the linear stability analysis are confirmed qualitatively by nonlinear dynamic MHD simulations. These results suggest the interesting possibility that substorm onset, or the thinning in the late growth phase, is caused by the release of a thermodynamic constraint without the (immediate) necessity of releasing the ideal MHD constraint. In the nonlinear regime the resistive and ideal developments differ in that the ideal mode does not lead to neutral line formation without the further release of the ideal MHD constraint; instead a thin current sheet forms. The isobaric constraint is critically discussed. Under perhaps more realistic adiabatic conditions the ideal mode appears to be stable but could be driven by external perturbations and thus generate the thin current sheet in the late growth phase, before a nonideal instability sets in.

Birn, Joachim↗
NASA NTRS2004

A Practical Approach to Obstacle Field Route Planning for Unmanned Rotorcraft

Autonomous helicopter operations in the near-earth environment will require robust and efficient obstacle field route planning. A method for obstacle field route planning is presented, which is composed of a mesh generation phase, a graph search phase, and a route refinement phase. The method mixes optimization and heuristics to obtain a satisfactory solution quickly. Simulations based on an unmanned helicopter model are presented.

Howlet, Jason K.↗
NASA NTRS2023

Magnetosonic ULF Waves With Anomalous Plasma–Magnetic Field Correlations: Standing Waves and Inhomogeneous Plasmas

Ultra-low frequency (ULF) wave observations across the heliosphere often rely on the sign of correlations between plasma (density/pressure) and magnetic field perturbations to distinguish between fast and slow magnetosonic modes. However, the assumptions behind this magnetohydrodynamic result are not always valid, particularly within the magnetosphere which is inhomogeneous and supports standing waves along the geomagnetic field. Through theory and a global simulation, we find both effects can result in anomalous plasma–magnetic field correlations. The interference pattern in standing waves can lead both body and surface magnetosonic waves to have different cross-phases than their constituent propagating waves. Furthermore, if the scale of gradients in the background are shorter than the wavelength or the waves are near-incompressible, then advection by the wave of inhomogeneities can overcome the wave's inherent sense of compression. These effects need to be allowed for and taken into account when applying the typical diagnostic to observations.

MHD Waves↗