Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high-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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Comparison of Implicit Collocation Methods for the Heat Equation

We combine a high-order compact finite difference scheme to approximate spatial derivatives arid collocation techniques for the time component to numerically solve the two dimensional heat equation. We use two approaches to implement the collocation methods. The first one is based on an explicit computation of the coefficients of polynomials and the second one relies on differential quadrature. We compare them by studying their merits and analyzing their numerical performance. All our computations, based on parallel algorithms, are carried out on the CRAY SV1.

Kouatchou, Jules↗
NASA NTRS1995

Chaotic Rotation of Nereid

The shape and spin of Neptune's outermost satellite Nereid are still unknown. Ground-based photometry indicates large brightness variations, but different observers report very different lightcurve amplitudes and periods. On the contrary, Voyager 2 images spanning 12 days show no evidence of variations greater than 0.1 mag. The latter suggest either that Nereid is nearly spherical, or that it is rotating slowly. We propose that tides have already despun Nereid's rotation to a period of a few weeks, during the time before the capture of Triton when Nereid was closer to Neptune. Since Nereid reached its present orbit, tides have further despun Nereid to a period on the order of a month. For Nereid's orbital eccentricity of 0.75, tidal evolution ceases when the spin period is still approximately 1/8 of the orbital period. Furthermore, the synchronous resonance becomes quite weak for such high eccentricities, along with other low-order spin orbit commensurabilities. In contrast, high-order resonances become very strong particularly the 6:1, 6.5:1, 7:1, 7.5:1, and 8:1 spin states. If Nereid departs by more than approximately 1% from a sphere, however, these resonances overlap, generating chaos. Our simulations show that Nereid is likely to be in chaotic rotation for any spin period longer than about 2 weeks.

Dobrovolskis, Anthony R.↗
NASA NTRS1995

State Estimation of Main Rotor Flap and Lead-Lag Using Accelerometers and Laser Transducers on the RASCAL UH-60 Helicopter

Modern rotorcraft flight control system designs which promise to yield high vehicle response bandwidth and good gust rejection can benefit from the use of rotor-state feedbacks. The measurement of main rotor blade motions is also desirable to validate and improve rotorcraft simulation models, to identify high-order linear flight dynamics models, to provide rotor system health monitoring; during flight test, and to provide for correlation with acoustic measurements from wind tunnel and flight tests. However, few attempts have been made to instrument a flight vehicle in this manner, and no previous system has had the robustness and accuracy required for these diverse applications. A rotor blade motion measurement and estimation system has been developed by NASA and the U.S. Army for use on the Rotorcraft Aircrew Systems Concepts Airborne Laboratory (RASCAL) helicopter. RASCAL is a UH-60 Blackhawk which is being modified at Ames Research Center in a phased development program for use in flight dynamics and controls, navigation, airspace management, and rotorcraft human factors research. The aircraft will feature a full-authority, digital, fly-by-wire research flight control system; a coupled ring laser gyro, differential GPS based navigation system; a stereoscopic color wide field of view helmet, mounted display; programmable panel mounted displays; and advanced navigation sensors. The rotor blade motion system is currently installed for data acquisition only, but will be integrated with the research flight control system when it is installed later this year.

Fletcher, Jay W.↗
NASA NTRS2000

Two-Species, 3D, MHD Simulation of Europa's Interaction with Jupiter's Magnetosphere

The interaction of Europa with the Jovian a magnetosphere has been studied by using a two species in ideal magnetohydrodynamic (MHD) numerical model. This model considers the upstream plasma in the Jovian magnetosphere and the molecular oxygen ions in the ionosphere of Europa, separately. We present results a from simulation studies, which take into account impact ionization, recombination, and the effect of a possible induced dipole magnetic field of Europa. The total mass loading of the magnetospheric flow and the ionization frequency used in the model are consistent with the estimates of Europa's ionosphere and atmosphere. The multi-species MHD equations are solved by using a finite volume, high-order, Godunov-type method on an adoptively refined unstructured grid, which allows a detailed modeling of the region near Europa's surface, while still resolving both the upstream region and the satellite's wake. We have paid special attention to the wake of Europa, in order to be able to make comparisons with the Galileo's E4 flyby observations, as well as other model calculations. The calculated escape flux of a O2+ down the tail was found to be about 5.6 x 10(exp 25) s(sup -1).

Liu, Yifan↗
NASA NTRS2001

A Highly Accurate Technique for the Treatment of Flow Equations at the Polar Axis in Cylindrical Coordinates using Series Expansions

Numerical methods for solving the flow equations in cylindrical or spherical coordinates should be able to capture the behavior of the exact solution near the regions where the particular form of the governing equations is singular. In this work we focus on the treatment of these numerical singularities for finite-differences methods by reinterpreting the regularity conditions developed in the context of pseudo-spectral methods. A generally applicable numerical method for treating the singularities present at the polar axis, when nonaxisymmetric flows are solved in cylindrical, coordinates using highly accurate finite differences schemes (e.g., Pade schemes) on non-staggered grids, is presented. Governing equations for the flow at the polar axis are derived using series expansions near r=0. The only information needed to calculate the coefficients in these equations are the values of the flow variables and their radial derivatives at the previous iteration (or time) level. These derivatives, which are multi-valued at the polar axis, are calculated without dropping the accuracy of the numerical method using a mapping of the flow domain from (0,R)*(0,2pi) to (-R,R)*(0,pi), where R is the radius of the computational domain. This allows the radial derivatives to be evaluated using high-order differencing schemes (e.g., compact schemes) at points located on the polar axis. The proposed technique is illustrated by results from simulations of laminar-forced jets and turbulent compressible jets using large eddy simulation (LES) methods. In term of the general robustness of the numerical method and smoothness of the solution close to the polar axis, the present results compare very favorably to similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r=0, following the method proposed recently by Mohseni and Colonius (1). Extension of the method described here for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical or spherical coordinates with finite-differences schemes of various level of accuracy is immediate.

Constantinescu, George S.↗
NASA NTRS2002

Efficient High Order Central Schemes for Multi-Dimensional Hamilton-Jacobi Equations: Talk Slides

This viewgraph presentation presents information on the attempt to produce high-order, efficient, central methods that scale well to high dimension. The central philosophy is that the equations should evolve to the point where the data is smooth. This is accomplished by a cyclic pattern of reconstruction, evolution, and re-projection. One dimensional and two dimensional representational methods are detailed, as well.

Bryson, Steve↗
NASA NTRS2001

Time Series Analysis of the UV Flickering in AGN

Many active galactic nuclei (AGN) exhibit large-amplitude luminosity fluctuations on short timescales. The fluctuations lead to a profound conclusion: The size of the emitting region is remarkably small. This observational fact is one of the pillars supporting the AGN paradigm: Prodigious amounts of gravitational potential energy are liberated in an accretion disk around a supermassive black hole. The goals of the research were to extract from the IUE Archive the very best observational characterizations of AGN flickering, and to use these to test and develop models for AGN variability. We hoped to answer these specific questions: 1) What does the intrinsic flickering continuum spectrum look like? 2) What do the intrinsic flickering emission-line profiles look like? 3) What is the power spectrum of the flickering? 4) What is the wavelength dependence of the power spectrum? 5) Is the flickering spectrum timescale dependent? and 6) What do the high-order cross correlation functions look like? A short summary of the papers produced by this research is presented.

Robinson, Edward L.↗
NASA NTRS1999

Numerical Solution of Incompressible Navier-Stokes Equations Using a Fractional-Step Approach

A fractional step method for the solution of steady and unsteady incompressible Navier-Stokes equations is outlined. The method is based on a finite volume formulation and uses the pressure in the cell center and the mass fluxes across the faces of each cell as dependent variables. Implicit treatment of convective and viscous terms in the momentum equations enables the numerical stability restrictions to be relaxed. The linearization error in the implicit solution of momentum equations is reduced by using three subiterations in order to achieve second order temporal accuracy for time-accurate calculations. In spatial discretizations of the momentum equations, a high-order (3rd and 5th) flux-difference splitting for the convective terms and a second-order central difference for the viscous terms are used. The resulting algebraic equations are solved with a line-relaxation scheme which allows the use of large time step. A four color ZEBRA scheme is employed after the line-relaxation procedure in the solution of the Poisson equation for pressure. This procedure is applied to a Couette flow problem using a distorted computational grid to show that the method minimizes grid effects. Additional benchmark cases include the unsteady laminar flow over a circular cylinder for Reynolds Numbers of 200, and a 3-D, steady, turbulent wingtip vortex wake propagation study. The solution algorithm does a very good job in resolving the vortex core when 5th-order upwind differencing and a modified production term in the Baldwin-Barth one-equation turbulence model are used with adequate grid resolution.

Kiris, Cetin↗
NASA NTRS1997

Large-Eddy Simulation of Transition to Turbulence in Boundary Layers

Large-eddy simulation results for laminar-to-turbulent transition in a spatially developing boundary layer are presented. The disturbances are ingested into a laminar flow through an unsteady suction-and-blowing strip. The filtered, three-dimensional time- dependent Navier-Stokes equations are integrated numerically using spectral, high-order finite-difference, and three-stage low-storage Runge-Kutta methods. The buffer-domain technique is used for the outflow boundary condition. The localized dynamic model used to parameterize the subgrid-scale stresses begins to have a significant impact at the beginning of the nonlinear transition (or intermittency) region. The flow structures commonly found in experiments are also observed in the present simulation; the computed linear instability modes and secondary instability lambda-vortex structures are in agreement with the experiments, and the streak-like-structures and turbulent statistics compare with both the experiments and the theory. The physics captured in the present LES are consistent with the experiments and the full Navier-Stokes simulation (DNS), at a significant fraction of the DNS cost. A comparison of the results obtained with several SGS models shows that the localized model gives accurate results both in a statistical sense and in terms of predicting the dynamics of the energy-carrying eddies, without ad hoc adjustments.

Huai, Xiao-Li↗
NASA NTRS2000

Selected contribution: a three-dimensional model for assessment of in vitro toxicity in balaena mysticetus renal tissue

This study established two- and three-dimensional renal proximal tubular cell cultures of the endangered species bowhead whale (Balaena mysticetus), developed SV40-transfected cultures, and cloned the 61-amino acid open reading frame for the metallothionein protein, the primary binding site for heavy metal contamination in mammals. Microgravity research, modulations in mechanical culture conditions (modeled microgravity), and shear stress have spawned innovative approaches to understanding the dynamics of cellular interactions, gene expression, and differentiation in several cellular systems. These investigations have led to the creation of ex vivo tissue models capable of serving as physiological research analogs for three-dimensional cellular interactions. These models are enabling studies in immune function, tissue modeling for basic research, and neoplasia. Three-dimensional cellular models emulate aspects of in vivo cellular architecture and physiology and may facilitate environmental toxicological studies aimed at elucidating biological functions and responses at the cellular level. Marine mammals occupy a significant ecological niche (72% of the Earth's surface is water) in terms of the potential for information on bioaccumulation and transport of terrestrial and marine environmental toxins in high-order vertebrates. Few ex vivo models of marine mammal physiology exist in vitro to accomplish the aforementioned studies. Techniques developed in this investigation, based on previous tissue modeling successes, may serve to facilitate similar research in other marine mammals.

Whales/genetics/physiology↗
NASA NTRS1999

622-Mbps Orthogonal Frequency Division Multiplexing Modulator Developed

The Communications Technology Division at the NASA Lewis Research Center is developing advanced electronic technologies for the space communications and remote sensing systems of tomorrow. As part of the continuing effort to advance the state-of-the art in satellite communications and remote sensing systems, Lewis is developing a programmable Orthogonal Frequency Division Multiplexing (OFDM) modulator card for high-data-rate communication links. The OFDM modulator is particularly suited to high data-rate downlinks to ground terminals or direct data downlinks from near-Earth science platforms. It can support data rates up to 622 megabits per second (Mbps) and high-order modulation schemes such as 16-ary quadrature amplitude modulation (16-ary QAM) or 8- phase shift keying (8PSK). High order modulations can obtain the bandwidth efficiency over the traditional binary phase shift keying (BPSK) or quadrature phase shift keying (QPSK) modulator schemes. The OFDM modulator architecture can also be precompensated for channel disturbances and alleviate amplitude degradations caused by nonlinear transponder characteristics.

Nguyen, Na T.↗
NASA NTRS1996

Statistical Mechanics and Dynamics of the Outer Solar System.I. The Jupiter/Saturn Zone

We report on numerical simulations designed to understand how the solar system evolved through a winnowing of planetesimals accreeted from the early solar nebula. This sorting process is driven by the energy and angular momentum and continues to the present day. We reconsider the existence and importance of stable niches in the Jupiter/Saturn Zone using greatly improved numerical techniques based on high-order optimized multi-step integration schemes coupled to roundoff error minimizing methods.

numerical simulations planetesimals planet formati↗
NASA NTRS1997

Direct Numerical Simulations of Boundary Layer Transition on a Flat Plate

In recent years the techniques of computational fluid dynamics (CFD) have been used to compute flows associated with geometrically complex configurations. However, success in terms of accuracy and reliability has been limited to cases where the effects of turbulence and transition could be modeled in a straightforward manner. Even in simple flows, the accurate computation of skin friction and heat transfer using existing turbulence models has proved to be a difficult task, one that has required extensive fine-tuning of the turbulence models used. In more complex flows (for example, in turbomachinery flows in which vortices and wakes impinge on airfoil surfaces causing periodic transitions from laminar to turbulent flow) the development of a model that accounts for all scales of turbulence and predicts the onset of transition is an extremely difficult task. Fortunately, current trends in computing suggest that it may be possible to perform direct simulations of turbulence and transition at moderate Reynolds numbers in some complex cases in the near future. This presentation will focus on direct simulations of transition and turbulence using high-order accurate finite-difference methods. The advantage of the finite-difference approach over spectral methods is that complex geometries can be treated in a straightforward manner. Additionally, finite-difference techniques are the prevailing methods in existing application codes. An application of accurate finite-difference methods to direct simulations of transition and turbulence in a spatially evolving boundary layer subjected to high levels of freestream turbulence will be presented.

Rai, Man Mohan↗
NASA NTRS2000

Cellular Responses to Mechanical Stress Selected Contribution: A Three-Dimensional Model for Assessment of in Vitro Toxicity in Balaena Mysticetus Renal Tissue

This study established two- and three-dimensional renal proximal tubular cell cultures of the endangered species bowhead whale (Balaena mysticetus), developed SV40-transfected cultures, and cloned the 61-amino acid open reading frame for the metallothionein protein, the primary binding site for heavy metal contamination in mammals. Microgravity research, modulations in mechanical culture conditions (modeled microgravity), and shear stress have spawned innovative approaches to understanding the dynamics of cellular interactions, gene expression, and differentiation in several cellular systems. These investigations have led to the creation of ex vivo tissue models capable of serving as physiological research analogs for three-dimensional cellular interactions. These models are enabling studies in immune function, tissue modeling for basic research, and neoplasia. Three-dimensional cellular models emulate aspects of in vivo cellular architecture and physiology and may facilitate environmental toxicological studies aimed at elucidating biological functions and responses at the cellular level. Marine mammals occupy a significant ecological niche (72% of the Earth's surface is water) in terms of the potential for information on bioaccumulation and transport of terrestrial and marine environmental toxins in high-order vertebrates. Few ex vivo models of marine mammal physiology exist in vitro to accomplish the aforementioned studies. Techniques developed in this investigation, based on previous tissue modeling successes, may serve to facilitate similar research in other marine mammals.

Goodwin, T. J.↗
DOE OSTI2025

Coupling Approaches with Non-matching Grids for Classical Linear Elasticity and Bond-based Peridynamic Models in 1D

Local-nonlocal coupling approaches provide a means to combine the computational efficiency of local models and the accuracy of nonlocal models. To facilitate the coupling of the two models, non-matching grids are often desirable as nonlocal grids usually require a finer resolution than local grids. In that case, it is often convenient to resort to interpolation operators so that models can exchange information in the overlap regions when nodes from the two grids do not coincide. This paper studies three existing coupling approaches, namely 1) a method that enforces matching displacements in an overlap region, 2) a variant that enforces a constraint on the stresses instead, and 3) a method that considers a variable horizon in the vicinity of the interfaces. Further, the effect of the interpolation order and of the grid ratio on the performance of the three coupling methods with non-matching grids is carefully studied on one-dimensional examples using polynomial manufactured solutions. The numerical results show that the degree of the interpolants should be chosen with care to avoid introducing additional modeling errors, or simply minimize these errors, in the coupling approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗