Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable order”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Three-Dimensional Simulations of Marangoni-Benard Convection in Small Containers by the Least-Squares Finite Element Method

This paper reports a numerical study of the Marangoni-Benard (MB) convection in a planar fluid layer. The least-squares finite element method (LSFEM) is employed to solve the three-dimensional Stokes equations and the energy equation. First, the governing equations are reduced to be first-order by introducing variables such as vorticity and heat fluxes. The resultant first-order system is then cast into a div-curl-grad formulation, and its ellipticity and permissible boundary conditions are readily proved. This numerical approach provides an equal-order discretization for velocity, pressure, vorticity, temperature, and heat conduction fluxes, and therefore can provide high fidelity solutions for the complex flow physics of the MB convection. Numerical results reported include the critical Marangoni numbers (M(sub ac)) for the onset of the convection in containers with various aspect ratios, and the planforms of supercritical MB flows. The numerical solutions compared favorably with the experimental results reported by Koschmieder et al..

Yu, Sheng-Tao↗
DOE OSTI2024

Analysis of single-mode Richtmyer–Meshkov instability using high-order incompressible vorticity-streamfunction and shock-capturing simulations

In this study, two- and three-dimensional simulation results obtained using a new high-order incompressible, variable-density vorticity–streamfunction (VS) method and data from previous ninth-order weighted essentially nonoscillatory (WENO) shock-capturing simulations are used to investigate the nonlinear dynamics of single-mode Richtmyer–Meshkov instability using a model of a Mach 1.3 air(acetone)/SF6 shock tube experiment. A comparison of the density fields from both simulations with the experimental images demonstrates very good agreement in the large-scale structure with both methods but differences in the small-scale structure. The WENO method captures the small-scale disordered structure observed in the experiment, while the VS method partially captures such structure and yields a strong rotating core. The perturbation amplitude growth from the simulations generally agrees well with the experiment. The simulation bubble and spike amplitudes agree well at early times. At later times, the WENO bubble amplitude is smaller than the VS amplitude and vice versa for the spike amplitude. The predictions of nonlinear single-mode instability growth models are shown to agree with the simulation amplitudes at early-to-intermediate times but underpredict the amplitudes at later times in the nonlinear regime. Visualizations of the mass fraction and enstrophy isosurfaces, velocity and vorticity fields, and baroclinic vorticity production and vortex stretching terms from the three-dimensional simulations indicate that, with the exception of the small-scale structure within the rollups, the VS and WENO results are in good agreement.

42 ENGINEERING↗
NASA NTRS1993

Chondrule formation by clumpy accretion onto the solar nebula

Chondrule textures and compositions appear to require rapid heating of precursor grain aggregates to temperatures in the range 1500 K to 2100 K, cooling times on the order of hours, and episodic and variable intensity events in order to produce chondrule rims and chemically distinct groups. Nebula shock waves have been proposed by Hood and Horanyi as a physical mechanism that may be capable of meeting the meteoritical constraints. Motivated by astronomical observations of the close environments of young stars, we suggest that the source of the nebula shock waves may be clumpy accretion onto the solar nebula - that is, episodic impacts onto the nebula by discrete cloud clumps with masses of at least 10(exp 22) g. If the cloud clumps are massive enough (10(exp 26) g), the resulting shockwave may be able to propagate to the midplane and process precursor aggregates residing in a dust sub-disk.

Boss, A. P.↗
NASA NTRS2011

Method for Measuring Collimator-Pointing Sensitivity to Temperature Changes

For a variety of applications, it is important to measure the sensitivity of the pointing of a beam emerging from a collimator, as a function of temperature changes. A straightforward method for carrying out this measurement is based on using interferometry for monitoring the changes in beam pointing, which presents its own problems. The added temperature dependence and complexity issues relating to using an interferometer are addressed by not using an interferometer in the first place. Instead, the collimator is made part of an arrangement that uses a minimum number of low-cost, off-the-shelf materials and by using a quad diode to measure changes in beam pointing. In order to minimize the influence of the test arrangement on the outcome of the measurement, several steps are taken. The collimator assembly is placed on top of a vertical, 1-m-long, fused silica tube. The quad diode is bonded to a fused silica bar, which, in turn, is bonded to the lower end of the fused silica tube. The lower end of the tube rests on a self-aligning support piece, while the upper end of the tube is kept against two rounded setscrew tips, using a soft rubber string. This ensures that very little stress is applied to the tube as the support structure changes dimensions due to thermal expansion. Light is delivered to the collimator through a bare fiber in order to minimize variable bending torque caused by a randomly relaxing, rigid fiber jacket. In order to separate the effect of temperature on the collimator assembly from the effect temperature has on the rest of the setup, multiple measurements are taken with the collimator assembly rotated from measurement to measurement. Laboratory testing, with 1-m spacing between the collimator and the quad diode, has shown that the sensitivity of the arrangement is better than 100 nm rms, over time spans of at least one hour, if the beam path is protected from atmospheric turbulence by a tube. The equivalent sensitivity to detecting changes in pointing angle is 100 nanoradians.

Abramovici, Alex↗
NASA NTRS2002

Coastal Atmosphere and Sea Time Series (CoASTS)

In this document, the first three years of a time series of bio-optical marine and atmospheric measurements are presented and analyzed. These measurements were performed from an oceanographic tower in the northern Adriatic Sea within the framework of the Coastal Atmosphere and Sea Time Series (CoASTS) project, an ocean color calibration and validation activity. The data set collected includes spectral measurements of the in-water apparent (diffuse attenuation coefficient, reflectance, Q-factor, etc.) and inherent (absorption and scattering coefficients) optical properties, as well as the concentrations of the main optical components (pigment and suspended matter concentrations). Clear seasonal patterns are exhibited by the marine quantities on which an appreciable short-term variability (on the order of a half day to one day) is superimposed. This short-term variability is well correlated with the changes in salinity at the surface resulting from the southward transport of freshwater coming from the northern rivers. Concentrations of chlorophyll alpha and total suspended matter span more than two orders of magnitude. The bio-optical characteristics of the measurement site pertain to both Case-I (about 64%) and Case-II (about 36%) waters, based on a relationship between the beam attenuation coefficient at 660nm and the chlorophyll alpha concentration. Empirical algorithms relating in-water remote sensing reflectance ratios and optical components or properties of interest (chlorophyll alpha, total suspended matter, and the diffuse attenuation coefficient) are presented.

Hooker, Stanford B.↗
NASA NTRS1976

On the optimal minimum order observer-based compensator and the limited state variable feedback controller

Four design problems are considered: two on the optimal minimum order observer-based compensator design and two on the optimal limited state variable feedback controller. The problem of designing an optimal discrete time linear time-invariant observer-based compensator for the regulation of an n dimensional linear discrete time time-invariant plant with m independent outputs is considered. This is a stochastic design problem to the extent that the initial plant state is assumed to be a random vector with known first and second order statistics. The compensator parameters are obtained by minimizing the expectation, with respect to the initial conditions, of the standard cost, quadratic in the state and control vectors with the inclusion of cross terms.

Llorens-Ortiz, B.↗
NASA NTRS1990

Compensating for pneumatic distortion in pressure sensing devices

A technique of compensating for pneumatic distortion in pressure sensing devices was developed and verified. This compensation allows conventional pressure sensing technology to obtain improved unsteady pressure measurements. Pressure distortion caused by frictional attenuation and pneumatic resonance within the sensing system makes obtaining unsteady pressure measurements by conventional sensors difficult. Most distortion occurs within the pneumatic tubing which transmits pressure impulses from the aircraft's surface to the measurement transducer. To avoid pneumatic distortion, experiment designers mount the pressure sensor at the surface of the aircraft, (called in-situ mounting). In-situ transducers cannot always fit in the available space and sometimes pneumatic tubing must be run from the aircraft's surface to the pressure transducer. A technique to measure unsteady pressure data using conventional pressure sensing technology was developed. A pneumatic distortion model is reduced to a low-order, state-variable model retaining most of the dynamic characteristics of the full model. The reduced-order model is coupled with results from minimum variance estimation theory to develop an algorithm to compensate for the effects of pneumatic distortion. Both postflight and real-time algorithms are developed and evaluated using simulated and flight data.

Whitmore, Stephen A.↗
NASA NTRS1990

Compensating for pneumatic distortion in pressure sensing devices

A general numerical technique for obtaining unsteady pressure measurements using conventional pressure sensing technology has been developed. A pneumatic distortion model, based on the Navier-Stokes equations of momentum and continuity, was reduced to a low-order, state-variable model retaining most of the dynamic characteristics of the full model. The reduced-order model is coupled with results from minimum variance estimation theory to develop an algorithm to compensate the effects of pneumatic distortion. Both postflight and real-time algorithms were developed and evaluated using simulated and flight data.

Whitmore, Stephen A.↗
NASA NTRS1990

Formulation of a minimum variance deconvolution technique for compensation of pneumatic distortion in pressure sensing devices

Increasingly, aircraft system designs require that aerodynamic parameters derived from pneumatic measurements be employed as control-system feedbacks. Such high frequency pressure measurements' accuracy is compromised by pressure distortion due to frictional attenuation and pneumatic resonance within the sensing system. A pneumatic distortion model is here formulated and reduced to a low-order state-variable model which retains most of the full model's dynamic characteristics. This reduced-order model is coupled with standard results from minimum variance estimation theory to develop an algorithm to compensate for pneumatic-distortion effects.

Whitmore, Stephen A.↗
NASA NTRS2015

Higher-Order Methods for Compressible Turbulent Flows Using Entropy Variables

Turbulent flows have a large range of spatial and temporal scales which need to be resolved in order to obtain accurate predictions. Higher-order methods can provide greater efficiency for simulations requiring high spatial and temporal resolution, allowing for solutions with fewer degrees of freedom and lower computational cost than traditional second-order computational fluid dynamics (CFD) methods.1 Higher-order methods have been widely used for turbulent flows. However, the reduced numerical stabilization present in higher-order schemes implies that special care needs to be taken in the development of numerical methods to suppress nonlinear instabilities.2–6 In this work we present the development of a higher-order space-time discontinuous Galerkin method with a focus on the aspects of our numerical scheme required for ensuring nonlinear stability for turbulent simulations at high Reynolds numbers.

Diosady, Laslo T.↗
DOE OSTI2021

An embedded variable step IMEX scheme for the incompressible Navier–Stokes equations

This report presents a series of implicit–explicit (IMEX) variable stepsize algorithms for the incompressible Navier–Stokes equations (NSE). Here, IMEX means the nonlinear term is treated fully explicitly while the remaining terms are treated implicitly. With the advent of new computer architectures there has been growing demand for low memory solvers of this type. The addition of time adaptivity improves the accuracy and greatly enhances the efficiency of the algorithm. We prove energy stability of an embedded first–second order IMEX pair. For the first order member of the pair, we prove stability for variable stepsizes, and analyze convergence. We believe this to be the first proof of this type for a variable stepsize IMEX scheme for the incompressible NSE. We then define and test a variable stepsize, variable order IMEX scheme using these methods. Our work contributes several firsts for IMEX NSE schemes, including an energy argument and error analysis of a two-step, variable stepsize method, and embedded error estimation for an IMEX multi-step method.

42 ENGINEERING↗
NASA NTRS2002

Scale Closure in Upper Ocean Optical Properties: From Single Particles to Ocean Color

Predictions of chlorophyll concentration from satellite ocean color are an indicator of primary productivity, with implications for foodwebs, fisheries, and the global carbon cycle. Models describing the relationship between optical properties and chlorophyll do not account for much of the optical variability observed in natural waters, because of the presence of seawater constituents that do not covary with phytoplankton pigments. in order to understand variability in these models, the optical contributions of seawater constituents were investigated. A combination of Mie theory and flow cytometry was used to determine the diameter, complex refractive index, and optical cross-sections of individual particles. In New England continental shelf waters, eukaryotic phytoplankton were the main particle contributors to absorption and scaftering. Minerals were the main contributor to backscattering (bb) in the spring, whereas in the summer both minerals and detritus contributed to bb. Synechococcus and heterotrophic bacteria were relatively unimportant optically. Seasonal differences in the spectral shape of remote sensing reflectance, Rrs, were contributed to approximately equally by eukaryotic phytoplankton absorption, dissolved absorption, and non-phytoplankton bb. Differences between measurements of bb and Prs and modeled values based on chlorophyll concentration were caused by higher dissolved absorption and non-phytoplankton bb than were assumed by the model.

Green, Rebecca E.↗
DOE OSTI2026

Coexistence of Synchronization and Stochasticity in Thermally Coupled Mott Oscillators

Synchronization is conventionally regarded as a mechanism for suppressing variability and enforcing order in coupled systems, from pendula and lasers to neurons and electronic oscillators. Here, we show that synchronization can also embed stochasticity at finer scales. We observe this phenomenon in thermally coupled VO 2 neuristors, where robust in-phase synchronization at the microsecond scale coexists with spike onset fluctuations at the nanosecond scale, with no fixed leader. The coexistence of order and disorder originates from stochastic domain-level physics of the insulator–metal and metal–insulator transitions, where local variations in transition temperature drive cycle-to-cycle randomness in nucleation, percolation, and relaxation. A stochastic domain model reproduces this effect by generating synchronized spike trains with random lead–lag jitter, and experimental interspike interval statistics confirm the persistence of fine-scale variability despite macroscopic phase locking. These findings establish that synchronization and stochasticity can coexist within the same physical platform, revealing hidden disorder within collective order. Furthermore, this insight reframes synchronization as not purely deterministic, but as a universal context where microscopic variability can persist, with implications for electronics, cryptography, and the fundamental physics of order–disorder coexistence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Compatible, energy conserving, bounds preserving remap of hydrodynamic fields for an extended ALE scheme

From the very origins of numerical hydrodynamics in the Lagrangian work of von Neumann and Richtmyer [83], the issue of total energy conservation as well as entropy production has been problematic. Because of well known problems with mesh deformation, Lagrangian schemes have evolved into Arbitrary Lagrangian–Eulerian (ALE) methods [39] that combine the best properties of Lagrangian and Eulerian methods. Energy issues have persisted for this class of methods. We believe that fundamental issues of energy conservation and entropy production in ALE require further examination. The context of the paper is an ALE scheme that is extended in the sense that it permits cyclic or periodic remap of data between grids of the same or differing connectivity. The principal design goals for a remap method then consist of total energy conservation, bounded internal energy, and compatibility of kinetic energy and momentum. We also have secondary objectives of limiting velocity and stress in a non-directional manner, keeping primitive variables monotone, and providing a higher than second order reconstruction of remapped variables. Particularly, the new contributions fall into three categories associated with: energy conservation and entropy production, reconstruction and bounds preservation of scalar and tensor fields, and conservative remap of nonlinear fields. Our paper presents a derivation of the methods, details of implementation, and numerical results for a number of test problems. The methods requires volume integration of polynomial functions in polytopal cells with planar facets, and the requisite expressions are derived for arbitrary order.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Phase Field Dislocation Dynamics (PFDD) version 2.x

This disclosure is for version 2.x of a mesoscale model called Phase Field Dislocation Dynamics (PFDD). PFDD is used for investigating deformation in nanoscale (grain sizes of ~300 nm and less) materials, such as metals and alloys. This approach models the motion and interaction of individual defects, namely dislocations, in the material using scalar-valued phase field variables, also called order parameters. The system is evolved through energy minimization thus the model calculates the total energy density in terms of the phase field variables. The energy minimization is completed using the Ginzburg-Landau equation, and is implemented with explicit time integration. The total system energy can be comprised of several terms, including the strain energy (which describes dislocation-dislocation interactions), the energy due to an applied stress (dislocation interactions with the applied stress), and a core/lattice (perfect dislocations) or generalized stacking fault (partial dislocations) energy (described the dislocation core structure). The latter term in particular may vary based on the crystal structure being modeled and is typically informed using lower length scale (e.g., atomistic) approaches, although no such (atomistic) calculations are completed within the PFDD algorithm. This basic formulation was previously reviewed by Los Alamos National Laboratory and released under license number C17113. This previously reviewed version we will henceforth refer to as PFDD v1.0. PFDD v1.0 consisted of 2 codes (one parallel and one serial) plus input files, all written in the C language. This new disclosure is addressing the next versions of the PFDD, versions 2.x. There have been several enhancements of PFDD v1.0, which are described here and included in the attached code, which we will refer to as PFDD v2.0. There are also several new features described here that are either planned or already in process and are expected to be subsequent releases, i.e., v2.1, v2.2, ...v2.x.

Hunter, Abigail↗
NASA NTRS1990

A fully Sinc-Galerkin method for Euler-Bernoulli beam models

A fully Sinc-Galerkin method in both space and time is presented for fourth-order time-dependent partial differential equations with fixed and cantilever boundary conditions. The Sinc discretizations for the second-order temporal problem and the fourth-order spatial problems are presented. Alternate formulations for variable parameter fourth-order problems are given which prove to be especially useful when applying the forward techniques to parameter recovery problems. The discrete system which corresponds to the time-dependent partial differential equations of interest are then formulated. Computational issues are discussed and a robust and efficient algorithm for solving the resulting matrix system is outlined. Numerical results which highlight the method are given for problems with both analytic and singular solutions as well as fixed and cantilever boundary conditions.

Smith, R. C.↗
NASA NTRS1999

Constraints on the Martian Plate Tectonic Hypothesis from Gravity and Topography Data

The Mars Global Surveyor Magnetic Fields Experiment/ Electron Reflectometer (MGS MAG/ER) experiment serendipitously discovered unanticipated and unprecedented regions of high amplitude crustal magnetic anomalies, indicating strong sources of remanent crustal magnetism. In one area of the southern hemisphere, the anomalies appear lineated and alternate in direction, resembling the stripes formed at terrestrial oceanic spread-ing regions. However, many significant differences exist. The inferred magnetization are easily an order of magnitude greater in strength than terrestrial counterparts. The width of the anomalies appears to be approximately 200 km, in comparison to a variable width of order 10-1000 km at terrestrial spreading centers. However, the spacecraft altitude of 100-200 km may be such that narrower anomalies are simply unresolved. Although the majority of strong anomalies are found in the southern highlands, there is no clear correlation with landforms at the surface. The lack of a correlation between magnetism and topography hinders the confident interpretation of magnetic sources. Additional information is contained in the original extended abstract.

Smrekar, S.↗