Search NASA⌕ Search

SEARCH · Search NASA

Results for “low order model”

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 613 records · Page 34

Turbulence Model Validation through Joint Experimental / Computational Studies of Separated Flow Over a Three-Dimensional Tapered Bump: Part II - Computational Investigation

This report serves as documentation of activities related to Computational Fluid Dynamics (CFD) analysis applied to the Smooth Body Separation Experiment designed to produce validation-quality experimental data from wind-tunnel tests performed at the University of Notre Dame. All CFD analysis utilized the Boeing CFD (BCFD) flow solver. BCFD was used to analyze flow over the Boeing Speed Bump with second-order accuracy in space and time. Inviscid fluxes were computed using the HLLE scheme for steady-state simulations, and the Low-Dissipation Roe scheme for time-accurate simulations. The working fluid was modeled as calorically perfect air using Sutherland's law for viscosity. Thermal diffusivity was determined by assuming a constant molecular Prandtl number of 0.72, and a turbulent Prandtl number of 0.9. Reynolds Averaged Navier-Stokes (RANS) simulations were conducted using the one-equation Spalart-Allmaras turbulence model with rotation correction3 and the quadratic constitutive relation (SA-RC-QCR), as well as Menter's two-equation, Shear-Stress Transport (SST) model. Time-accurate computations were conducted using Delayed Detached eddy Simulation (DDES) with the one-equation RANS model, and Large Eddy Stress balance (LESb) with the two-equation RANS model. Implicit dual-time stepping was used to achieve second-order temporal accuracy. The experimental investigation conducted during the course of this study is described in detail in an accompanying report.

Flow separation↗

Reduced size first-order subsonic and supersonic aeroelastic modeling

Various aeroelastic, aeroservoelastic, dynamic-response, and sensitivity analyses are based on a time-domain first-order (state-space) formulation of the equations of motion. The formulation of this paper is based on the minimum-state (MS) aerodynamic approximation method, which yields a low number of aerodynamic augmenting states. Modifications of the MS and the physical weighting procedures make the modeling method even more attractive. The flexibility of constraint selection is increased without increasing the approximation problem size; the accuracy of dynamic residualization of high-frequency modes is improved; and the resulting model is less sensitive to parametric changes in subsequent analyses. Applications to subsonic and supersonic cases demonstrate the generality, flexibility, accuracy, and efficiency of the method.

Karpel, Mordechay↗

Comparing numerical accuracy and stability for different horizontal discretizations in MPAS-Ocean

This manuscript investigates the effectiveness of two possible horizontal discretizations for the global ocean model MPAS-Ocean, both applied to Spherical Centroidal Voronoi Tessellations (SCVTs). The first discretization is TRiSK, a C-grid, finite-volume method, that possesses many desirable mimetic properties, but has a low order accuracy. The second discretization was introduced for the first time by Peixoto (2016), and consists of modifications to the TRiSK scheme designed to achieve at least first-order accuracy in the L ∞ norm, with the loss of some of the mimetic properties. Tests on shallow-water and primitive-equation models show that the scheme due to Peixoto is indeed more accurate, but presents stability issues with respect to TRiSK. Here, TRiSK is indeed found to be often more stable in time and more robust with respect to errors in the geometric properties of the grid.

97 MATHEMATICS AND COMPUTING↗

Research highlights: June 1990 - May 1991

Linear instability calculations at MSFC have suggested that the Geophysical Fluid Flow Cell (GFFC) should exhibit classic baroclinic instability at accessible parameter settings. Interest was in the mechanisms of transition to temporal chaos and the evolution of spatio-temporal chaos. In order to understand more about such transitions, high resolution numerical experiments for the physically simplest model of two layer baroclinic instability were conducted. This model has the advantage that the numerical code is exponentially convergent and can be efficiently run for very long times, enabling the study of chaotic attractors without the often devastating effects of low-order trunction found in many previous studies. Numerical algorithms for implementing an empirical orthogonal function (EOF) analysis of the high resolution numerical results were completed. Under conditions of rapid rotation and relatively low differential heating, convection in a spherical shell takes place as columnar banana cells wrapped around the annular gap, but with axes oriented along the axis of rotation; these were clearly evident in the GFFC experiments. The results of recent numerical simulations of columnar convection and future research plans are presented.

Source record↗

Mass accommodation and gas–particle partitioning in secondary organic aerosols: dependence on diffusivity, volatility, particle-phase reactions, and penetration depth

Abstract. Mass accommodation is an essential process for gas–particle partitioning of organic compounds in secondary organic aerosols (SOA). The mass accommodation coefficient is commonly described as the probability of a gas molecule colliding with the surface to enter the particle phase. It is often applied, however, without specifying if and how deep a molecule has to penetrate beneath the surface to be regarded as being incorporated into the condensed phase (adsorption vs. absorption). While this aspect is usually not critical for liquid particles with rapid surface–bulk exchange, it can be important for viscous semi-solid or glassy solid particles to distinguish and resolve the kinetics of accommodation at the surface, transfer across the gas–particle interface, and further transport into the particle bulk. For this purpose, we introduce a novel parameter: an effective mass accommodation coefficient αeff that depends on penetration depth and is a function of surface accommodation coefficient, volatility, bulk diffusivity, and particle-phase reaction rate coefficient. Application of αeff in the traditional Fuchs–Sutugin approximation of mass-transport kinetics at the gas–particle interface yields SOA partitioning results that are consistent with a detailed kinetic multilayer model (kinetic multilayer model of gas–particle interactions in aerosols and clouds, KM-GAP; Shiraiwa et al., 2012) and two-film model solutions (Model for Simulating Aerosol Interactions and Chemistry, MOSAIC; Zaveri et al., 2014) but deviate substantially from earlier modeling approaches not considering the influence of penetration depth and related parameters. For highly viscous or semi-solid particles, we show that the effective mass accommodation coefficient remains similar to the surface accommodation coefficient in the case of low-volatility compounds, whereas it can decrease by several orders of magnitude in the case of semi-volatile compounds. Such effects can explain apparent inconsistencies between earlier studies deriving mass accommodation coefficients from experimental data or from molecular dynamics simulations. Our findings challenge the approach of traditional SOA models using the Fuchs–Sutugin approximation of mass transfer kinetics with a fixed mass accommodation coefficient, regardless of particle phase state and penetration depth. The effective mass accommodation coefficient introduced in this study provides an efficient new way of accounting for the influence of volatility, diffusivity, and particle-phase reactions on SOA partitioning in process models as well as in regional and global air quality models. While kinetic limitations may not be critical for partitioning into liquid SOA particles in the planetary boundary layer (PBL), the effects are likely important for amorphous semi-solid or glassy SOA in the free and upper troposphere (FT–UT) as well as in the PBL at low relative humidity and low temperature.

54 ENVIRONMENTAL SCIENCES↗

Confinement, reduced entanglement, and spin-glass order in a random quantum spin-ice model

In this work, we study an effective spin model derived perturbatively from a random transverse-field Ising model on the pyrochlore lattice. The model consists of spin configurations on the pyrochlore lattice, restricted to the spin-ice subspace, with spins interacting with random Ising exchange couplings as well as ring exchanges along the hexagons of the lattice. This model is studied by exact diagonalization up to N=64 site systems. Furthermore, we calculate spin-glass correlation functions and local entanglement entropy S T between spins in a single tetrahedron and the rest of the system. We find that the model undergoes two phase transitions. At weak randomness the model is in a quantum spin-ice phase where S T =ln 6. Increasing randomness, at low transverse fields, first leads to a frozen phase, with long-range spin-glass order and S T =ln 2 corresponding to the Cat states associated with Ising order. Further increase in randomness leads to a random resonating-hexagon phase with a frozen backbone of spins and a broad distribution of entanglement entropies. The implications of these studies for non-Kramers rare-earth pyrochlores are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Three-dimensional computer model for the atmospheric general circulation experiment

An efficient, flexible, three-dimensional, hydrodynamic, computer code has been developed for a spherical cap geometry. The code will be used to simulate NASA's Atmospheric General Circulation Experiment (AGCE). The AGCE is a spherical, baroclinic experiment which will model the large-scale dynamics of our atmosphere; it has been proposed to NASA for future Spacelab flights. In the AGCE a radial dielectric body force will simulate gravity, with hot fluid tending to move outwards. In order that this force be dominant, the AGCE must be operated in a low gravity environment such as Spacelab. The full potential of the AGCE will only be realized by working in conjunction with an accurate computer model. Proposed experimental parameter settings will be checked first using model runs. Then actual experimental results will be compared with the model predictions. This interaction between experiment and theory will be very valuable in determining the nature of the AGCE flows and hence their relationship to analytical theories and actual atmospheric dynamics.

Roberts, G. O.↗

Regulatory Decision-Making Process Using 'Weight of Evidence: An Evidence Integration Approach' by the State of Washington at the Hanford Site - 20163

Performance Assessment and Environmental Impact Statement Modeling, Human and Ecological Risk Assessments are the tools used by the Washington State Department of Ecology (Ecology) to make regulatory decisions for permitting and cleanup. These efforts all inherently have large amounts of uncertainty no matter how many alternative approaches are attempted. Professional feedback from tribal governments, stakeholders, federal and state regulators, the public and independent peer review groups are taken into account in this process of informed decision making. The process also includes full evaluation of the comments and the associated reviews made by the various peer review groups such as the Low -Level Waste Disposal Facility Federal Review Group (DOE order 435.1) and the Nuclear Regulatory Commission (NRC). The evolution of technology has produced quicker and more efficient fate and transport modeling tools with related quantification of uncertainties and their significance. Groundwater and vadose zone models have been expanded across the entire Hanford site. Many of the technical approaches identified in the Washington State Regulations are updated using current knowledge to better understand the complexities of the Hanford Site using the latest state of the art tools for risk assessment and modeling. More emphasis and time are directed on defining uncertainty and then managing the uncertainty. One added tool is to run sensitivity cases based on informed decision making. The weight of evidence approach considers all relevant information in an integrative assessment that takes into account the kinds of evidence available, the quality and quantity, the uncertainty and error associated with each type, and how they fit together. Evidence of integration involves looking at systematic reviews and evaluating their strengths and weaknesses. In conclusion, Ecology, Washington State's lead regulatory agency at the Hanford site, has had to modify its site closure approach compared to the rest of the state. A more holistic approach dealing with the Hanford site as an integrated system has evolved leading to a defensible informed approach for regulatory closure decisions. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Plasma Reactor Modeling and Validation Experiments

Plasma processing is a key processing stop in integrated circuit manufacturing. Low pressure, high density plum reactors are widely used for etching and deposition. Inductively coupled plasma (ICP) source has become popular recently in many processing applications. In order to accelerate equipment and process design, an understanding of the physics and chemistry, particularly, plasma power coupling, plasma and processing uniformity and mechanism is important. This understanding is facilitated by comprehensive modeling and simulation as well as plasma diagnostics to provide the necessary data for model validation which are addressed in this presentation. We have developed a complete code for simulating an ICP reactor and the model consists of transport of electrons, ions, and neutrals, Poisson's equation, and Maxwell's equation along with gas flow and energy equations. Results will be presented for chlorine and fluorocarbon plasmas and compared with data from Langmuir probe, mass spectrometry and FTIR.

Meyyappan, M.↗

Breakdown of the Static Dielectric Screening Approximation of Coulomb Interactions in Atomically Thin Semiconductors

Coulomb interactions in atomically thin materials are remarkably sensitive to variations in the dielectric screening of the environment, which can be used to control exotic quantum many-body phases and engineer exciton potential landscapes. For decades, static or frequency-independent approximations of the dielectric response, where increased dielectric screening is predicted to cause an energy redshift of the exciton resonance, have been sufficient. These approximations were first applied to quantum wells and were more recently extended with initial success to layered transition metal dichalcogenides (TMDs). Here, we use charge-tunable exciton resonances to investigate screening effects in TMD monolayers embedded in materials with low-frequency dielectric constants ranging from 4 to more than 1000, a range of 2 orders of magnitude larger than in previous studies. In contrast to the redshift predicted by static models, we observe a blueshift of the exciton resonance exceeding 30 meV in higher dielectric constant environments. We explain our observations by introducing a dynamical screening model based on a solution to the Bethe-Salpeter equation (BSE). When dynamical effects are strong, we find that the exciton binding energy remains mostly controlled by the low-frequency dielectric response, while the exciton self-energy is dominated by the high-frequency one. Our results supplant the understanding of screening in layered materials and their heterostructures, introduce a knob to tune selected many-body effects, and reshape the framework for detecting and controlling correlated quantum many-body states and designing optoelectronic and quantum devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Data-Driven Discovery of Active Nematic Hydrodynamics

Active nematics can be modeled using phenomenological continuum theories that account for the dynamics of the nematic director and fluid velocity through partial differential equations (PDEs). While these models provide a statistical description of the experiments, the relevant terms in the PDEs and their parameters are usually identified indirectly. We adapt a recently developed method to automatically identify optimal continuum models for active nematics directly from spatio-temporal data, via sparse regression of the coarse-grained fields onto generic low order PDEs. After extensive benchmarking, we apply the method to experiments with microtubule-based active nematics, finding a surprisingly minimal description of the system. Furthermore, our approach can be generalized to gain insights into active gels, microswimmers, and diverse other experimental active matter systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dynamic Pressure Drag on Rising Buoyant Thermals in a Neutrally Stable Environment

Abstract This study examines dynamic pressure drag on rising dry buoyant thermals. A theoretical expression for drag coefficient C d as a function of several other nondimensional parameters governing thermal dynamics is derived based on combining the thermal momentum budget with the similarity theory of Scorer. Using values for these nondimensional parameters from previous studies, the theory suggests drag on thermals is small relative to that on solid spheres in laminar or turbulent flow. Two sets of numerical simulations of thermals in an unstratified, neutrally stable environment using an LES configuration of the Cloud Model 1 (CM1) are analyzed. One set has a relatively low effective Reynolds number R e and the other has an order-of-magnitude-higher R e ; these produce laminar and turbulent resolved-scale flows, respectively. Consistent with the theoretical C d , the magnitude of drag is small in all simulations. However, whereas the laminar thermals have C d ≈ 0.01, the turbulent thermals have weakly negative drag ( C d ≈ −0.1). This difference is explained by the laminar thermals having near vertical symmetry but the turbulent thermals exhibiting considerable vertical asymmetry of their azimuthally averaged flows. In the laminar thermals, buoyancy rapidly becomes concentrated around the main centers of rotation located along the horizontal central axis, leading to expansion of thermals via baroclinic vorticity generation but doing little to break vertical symmetry of the flow. Vertical asymmetry of the azimuthally averaged flow of turbulent thermals is attributed mainly to small-scale resolved eddies that are concentrated in the upper part of the thermals.

Meteorology & Atmospheric Sciences↗

Nonequilibrium Crystal Growth During Laser Annealing of Silicon

The Systematic studies of dopant (Group III, IV, or V) incorporation into silicon as a result of pulsed-laser annealing are reviewed. At low concentrations, dopant profiles measured after laser annealing can be compared to model calculations for redistribution by liquid-phase diffusion in order to determine the interfacial distribution coefficient K('). At high concentrations there is a limit (C sub s) to dopant incorporation into substitutional lattice sites. The mechanisms which limit dopant incorporation are discussed as well as the dependence of K(') and (max) on crystal orientation. Finally, measurements of interface temperatures, temperature profiles, and thermal gradients during rapid solidification are summarized.

White, C. W.↗

Modeling and estimation of a low degree geopotential model from terrestrial gravity data

The development of appropriate modeling and adjustment procedures for the estimation of harmonic coefficients of the geopotential, from surface gravity data was studied, in order to provide an optimum way of utilizing the terrestrial gravity information in combination solutions currently developed at NASA/Goddard Space Flight Center, for use in the TOPEX/POSEIDON mission. The mathematical modeling was based on the fundamental boundary condition of the linearized Molodensky boundary value problem. Atmospheric and ellipsoidal corrections were applied to the surface anomalies. Terrestrial gravity solutions were found to be in good agreement with the satellite ones over areas which are well surveyed (gravimetrically), such as North America or Australia. However, systematic differences between the terrestrial only models and GEMT1, over extended regions in Africa, the Soviet Union, and China were found. In Africa, gravity anomaly differences on the order of 20 mgals and undulation differences on the order of 15 meters, over regions extending 2000 km in diameter, occur. Comparisons of the GEMT1 implied undulations with 32 well distributed Doppler derived undulations gave an RMS difference of 2.6 m, while corresponding comparison with undulations implied by the terrestrial solution gave RMS difference on the order of 15 m, which implies that the terrestrial data in that region are substantially in error.

Pavlis, Nikolaos K.↗

Krylov model reduction algorithm for undamped structural dynamics systems

Krylov vectors furnish an efficient basis for eigenvalue analysis and model reduction of structural dynamics systems. The reduced-order model obtained by the present Krylov model-reduction algorithm for an undamped structural-dynamics system is found to match low-frequency moments. The transformed system equation in Krylov coordinates reflects the structure of a tandem system.

Craig, Roy R., Jr.↗

Modeling a simple coronal streamer during whole sun month

The solar minimum streamer structure observed during the whole sun month was modeled. The Van de Hulst inversion was used in order to determine the coronal electron density profiles and scale-height temperature profiles. The axisymmetric magnetostatic model of Gibson, Bagenal and Low was also used. The density, temperature, and magnetic field distribution were quantified using both coronal white light data and photospheric magnetic field data from the Wilcox Solar Observatory. The densities and temperatures obtained by the Van de Hulst and magnetostatic models are compared to the magnetic field predicted by the magnetostatic model to a potential field extrapolated from the photosphere.

Gibson, S. E.↗

Noise Generation by Fans with Supersonic Tip Speeds

Fan noise continues to be a significant issue for commercial aircraft engines and there still exists a requirement for improved understanding of the fundamental issues associated with fan noise source mechanisms. At the present time, most of the prediction methods identify the dominant acoustic sources to be associated with the stator vanes or blade trailing edges which are downstream of the fan face. However recent studies have shown that acoustic waves are significantly attenuated as they propagate upstream through a rotor, and if the appropriate corrections are applied, sound radiation from the engine inlet is significantly underpredicted. The prediction models can only be applied to fans with subsonic tip speeds. In contrast, most aircraft engines have fan tip speeds which are transonic and this implies an even higher attenuation for upstream propagating acoustic waves. Consequently understanding how sound propagates upstream through the fan is an important, and not well understood phenomena. The objective of this study is to provide improved insight into the upstream propagation effects through a rotor which are relevant to full scale engines. The focus of this study is on broadband fan noise generated by boundary layer turbulence interacting with the trailing edges of the fan blades. If this source mechanism is important upstream of the fan, the sound must propagate upstream through a transonic non uniform flow which includes large gradients and non linearities. Developing acoustic propagation models in this type of flow is challenging and currently limited to low frequency applications, where the frequency is of the same order as the blade passing frequency of the fan. For trailing edge noise, much higher frequencies are relevant and so a suitable approach needs to be developed, which is not limited by an unacceptably large computational effort. In this study we are in the process of developing a computational method which applies for the high frequencies of interest, and allows for any type of flow field associated with the fan. In this progress report the approach to be used and the basic equations will be presented. Some initial results will be given, but these are preliminary and need further verification.

Glegg, Stewart↗

(abstract) Timing and Extent of Wurm-Wisconsinan Deglaciation and Subsequent Rebound in Antarctica

Recent information from glacio-geomorphologists, when coupled to new glacio-sedimentary dating of the grounded ice sheet in the Ross Embayment at the last glacial maximum, leads to the conclusion that global geophysical models might require some revision. We investigate such new models and the possible implications for observable geodetic signatures, such as rebound-driven temporal variations in the Earth's low-order gravity field and in vertical crustal motions.

Antarctica↗