Search NASA⌕ Search

SEARCH · Search NASA

Results for “time-dependent problem”

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 415 records · Page 23

The photospheric magnetic flux budget

The ensemble of bipolar regions and the magnetic network both contain a substantial and strongly variable part of the photospheric magnetic flux at any phase in the solar cycle. The time-dependent distribution of the magnetic flux over and within these components reflects the action of the dynamo operating in the solar interior. We perform a quantitative comparison of the flux emerging in the ensemble of magnetic bipoles with the observed flux content of the solar photosphere. We discuss the photospheric flux budget in terms of flux appearance and disappearance, and argue that a nonlinear dependence exists between the flux present in the photosphere and the rate of flux appearance and disappearance. In this context, we discuss the problem of making quantitative statements about dynamos in cool stars other than the Sun.

Schrijver, C. J.↗

Magnetohydrodynamic Riemann problem and the structure of the magnetic reconnection layer

We present a complete solution for a set of magnetohydrodynamic (MHD) Riemann problems in which the upstream and downstream states have the same total pressure, and in which the normal component of the magnetic field is very small. These solutions are pertinent to subfast flows in the earth's magnetic tail and near the magnetopause. In a coplanar situation a family of solutions exists that depend on two parameters as well as on dissapation mechanisms. In the parallel case the transverse magnetic field either does not change direction or changes the direction twice by involving two intermediate shocks. In the antiparallel case an intermediate shock is always required, except when the solution consists of two switch-off shocks. In a noncoplanar case the solution is not self-similar as a function of x/t, but continues to evolve. At early times the evolution is similar to the coplanar case. In general two time-dependent intermediate shocks are required to rotate the magnetic fields. The velocity shear has a strong effect on the Riemann solution. In some cases no Riemann solution can exist because of the cavitation caused by the slow refraction waves. The calculated magnetopause structure resembles the observed structure for northward interplanetary magnetic field (IMF). However, for southward IMF, the MHD result shows the existence of a depletion layer, which is not supported by observations. We also show that on the magnetosheath side, the Walen relation, which is exact for a rotational discontinuity, can also be well satisfied by a slow shock, an intermediate shock, or the head of a slow rarefaction wave.

Wu, C. C.↗

Analytical and phenomenological studies of rotating turbulence

A framework, which combines mathematical analysis, closure theory, and phenomenological treatment, is developed to study the spectral transfer process and reduction of dimensionality in turbulent flows that are subject to rotation. First, we outline a mathematical procedure that is particularly appropriate for problems with two disparate time scales. The approach which is based on the Green's method leads to the Poincare velocity variables and the Poincare transformation when applied to rotating turbulence. The effects of the rotation are now reflected in the modifications to the convolution of a nonlinear term. The Poincare transformed equations are used to obtain a time-dependent analog of the Taylor-Proudman theorem valid in the asymptotic limit when the non-dimensional parameter mu is identical to Omega(t) approaches infinity (Omega is the rotation rate and t is the time). The 'split' of the energy transfer in both direct and inverse directions is established. Secondly, we apply the Eddy-Damped-Quasinormal-Markovian (EDQNM) closure to the Poincare transformed Euler/Navier-Stokes equations. This closure leads to expressions for the spectral energy transfer. In particular, an unique triple velocity decorrelation time is derived with an explicit dependence on the rotation rate. This provides an important input for applying the phenomenological treatment of Zhou. In order to characterize the relative strength of rotation, another non-dimensional number, a spectral Rossby number, which is defined as the ratio of rotation and turbulence time scales, is introduced. Finally, the energy spectrum and the spectral eddy viscosity are deduced.

Mahalov, Alex↗

Two-Dimensional Dynamic Simulation of a Continuous Foil Bearing

In this paper, the two dimensional(radial and circumferential) transient Navier-Stokes equations are used to solve the hydrodynamic problem in conjunction with the time dependent motion of the journal, and the deformable, spring supported foil. The elastic deformation of the foil and its supports are simulated by a finite element model. The time-dependent Navier-Stokes formulation is used to solve for the interaction between the fluid lubricant, the motion of the journal and the deformable foil boundary. The steady state, the quasi-transient and the full transient dynamic simulation of the foil-fluid journal interaction are examined on a comparative basis. For the steady state simulation, the fluid lubricant pressures are evaluated for a particular journal position, by means of an iterative scheme until convergence is achieved in both the fluid pressures and the corresponding foil deformation. For the quasi-transient case, the transient motion of the journal is calculated using a numerical integration scheme for the velocity and displacement of the journal. The deformation of the foil is evaluated through numerical iteration in feedback mode with the fluid film pressure generated by the journal motion until convergence at every time step is achieved. For the full transient simulation, a parallel real-time integration scheme is used to evaluate simultaneously the new journal position and the new deformed shape of the foil at each time step. The pressure of the fluid lubricant is iterated jointly with the corresponding journal position and the deformed foil geometry until convergence is achieved. A variable time-stepping Newmark-Beta integration procedure is used to evaluate the transient dynamics at each time step of the bearing.

Braun, M. Jack↗

Computational Issues Associated with Temporally Deforming Geometries Such as Thrust Vectoring Nozzles

During the past decade, computational simulation of fluid flow around complex configurations has progressed significantly and many notable successes have been reported, however, unsteady time-dependent solutions are not easily obtainable. The present effort involves unsteady time dependent simulation of temporally deforming geometries. Grid generation for a complex configuration can be a time consuming process and temporally varying geometries necessitate the regeneration of such grids for every time step. Traditional grid generation techniques have been tried and demonstrated to be inadequate to such simulations. Non-Uniform Rational B-splines (NURBS) based techniques provide a compact and accurate representation of the geometry. This definition can be coupled with a distribution mesh for a user defined spacing. The present method greatly reduces cpu requirements for time dependent remeshing, facilitating the simulation of more complex unsteady problems. A thrust vectoring nozzle has been chosen to demonstrate the capability as it is of current interest in the aerospace industry for better maneuverability of fighter aircraft in close combat and in post stall regimes. This current effort is the first step towards multidisciplinary design optimization which involves coupling the aerodynamic heat transfer and structural analysis techniques. Applications include simulation of temporally deforming bodies and aeroelastic problems.

Boyalakuntla, Kishore↗

Time-Dependent Parabolic Finite Difference Formulation for Harmonic Sound Propagation in a Two-Dimensional Duct with Flow

An explicit finite difference real time iteration scheme is developed to study harmonic sound propagation in aircraft engine nacelles. To reduce storage requirements for future large 3D problems, the time dependent potential form of the acoustic wave equation is used. To insure that the finite difference scheme is both explicit and stable for a harmonic monochromatic sound field, a parabolic (in time) approximation is introduced to reduce the order of the governing equation. The analysis begins with a harmonic sound source radiating into a quiescent duct. This fully explicit iteration method then calculates stepwise in time to obtain the 'steady state' harmonic solutions of the acoustic field. For stability, applications of conventional impedance boundary conditions requires coupling to explicit hyperbolic difference equations at the boundary. The introduction of the time parameter eliminates the large matrix storage requirements normally associated with frequency domain solutions, and time marching attains the steady-state quickly enough to make the method favorable when compared to frequency domain methods. For validation, this transient-frequency domain method is applied to sound propagation in a 2D hard wall duct with plug flow.

Kreider, Kevin L.↗

Gravity Induced Formation of Concentration Gradients in Supersaturated Binary Solutions

Experimental and theoretical studies of the formation of solute concentration gradient in supersaturated binary solutions in a gravitational field were carried out. The formation of solute concentration gradient was associated with the gravity induced redistribution of subcritical solute clusters. The birth-death process of the new solute-rich phase domains (subcritical solute clusters) was described in terms of the time-dependent Ginzburg Landau model developed for metastable state relaxation in binary (solute + solvent) non-critical solutions in the presence of a gravitational field. A new mathematical Ansatz was developed for solution of the model equations. This Ansatz has allowed to approach for the first time the following important problems: (1) Microstructure of solute distribution inside of the subcritical solute clusters. The analytical results obtained demonstrate that solute inside of the subcritical solute clusters is heterogeneously distributed with a spatially periodic structure. (2) Macrostructure of the solute subcritical clusters distribution in a gravitational field. The subcritical solute clusters are found to be distributed heterogeneously in a gravitational field. This heterogeneity, which is due to the heterogeneous birth-death process of the subcritical solute clusters in a gravitational field, initiates a noticeable solute concentration gradient in vertical columns of supersaturated binary solutions. An analysis and comparison of theoretical results and experimental data related to the solute concentration gradient formation in a gravitational field are presented. It is also demonstrated that the critical radius of solute clusters (radius of nucleation) and induction time are gravity-dependent.

Izmailov, Alexander F.↗

Statistical Perspectives on Stratospheric Transport

Long-lived tropospheric source gases, such as nitrous oxide, enter the stratosphere through the tropical tropopause, are transported throughout the stratosphere by the Brewer-Dobson circulation, and are photochemically destroyed in the upper stratosphere. These chemical constituents, or "tracers" can be used to track mixing and transport by the stratospheric winds. Much of our understanding about the stratospheric circulation is based on large scale gradients and other spatial features in tracer fields constructed from satellite measurements. The point of view presented in this paper is different, but complementary, in that transport is described in terms of tracer probability distribution functions (PDFs). The PDF is computed from the measurements, and is proportional to the area occupied by tracer values in a given range. The flavor of this paper is tutorial, and the ideas are illustrated with several examples of transport-related phenomena, annotated with remarks that summarize the main point or suggest new directions. One example shows how the multimodal shape of the PDF gives information about the different branches of the circulation. Another example shows how the statistics of fluctuations from the most probable tracer value give insight into mixing between different regions of the atmosphere. Also included is an analysis of the time-dependence of the PDF during the onset and decline of the winter circulation, and a study of how "bursts" in the circulation are reflected in transient periods of rapid evolution of the PDF. The dependence of the statistics on location and time are also shown to be important for practical problems related to statistical robustness and satellite sampling. The examples illustrate how physically-based statistical analysis can shed some light on aspects of stratospheric transport that may not be obvious or quantifiable with other types of analyses. An important motivation for the work presented here is the need for synthesis of the large and growing database of observations of the atmosphere and the vast quantities of output generated by atmospheric models.

Sparling, L. C.↗

An Implicit Characteristic Based Method for Electromagnetics

An implicit characteristic-based approach for numerical solution of Maxwell's time-dependent curl equations in flux conservative form is introduced. This method combines a characteristic based finite difference spatial approximation with an implicit lower-upper approximate factorization (LU/AF) time integration scheme. This approach is advantageous for three-dimensional applications because the characteristic differencing enables a two-factor approximate factorization that retains its unconditional stability in three space dimensions, and it does not require solution of tridiagonal systems. Results are given both for a Fourier analysis of stability, damping and dispersion properties, and for one-dimensional model problems involving propagation and scattering for free space and dielectric materials using both uniform and nonuniform grids. The explicit Finite Difference Time Domain Method (FDTD) algorithm is used as a convenient reference algorithm for comparison. The one-dimensional results indicate that for low frequency problems on a highly resolved uniform or nonuniform grid, this LU/AF algorithm can produce accurate solutions at Courant numbers significantly greater than one, with a corresponding improvement in efficiency for simulating a given period of time. This approach appears promising for development of dispersion optimized LU/AF schemes for three dimensional applications.

Beggs, John H.↗

Exact Solutions to Time-dependent Mdps

We describe an extension of the Markov decision process model in which a continuous time dimension is included in the state space. This allows for the representation and exact solution of a wide range of problems in which transitions or rewards vary over time. We examine problems based on route planning with public transportation and telescope observation scheduling.

Boyan, Justin A.↗

Scientific Visualization of Large Scale Unsteady Fluid Flow

In a numerical flow simulation. it is common to generate several thousand time steps of unsteady (time-dependent) flow data. Each time step may require tens to hundreds of megabytes for disk storage, and the total disk requirement for storing the unsteady flow data may be hundreds of gigabytes. Interactive visualization of unsteady flow data of this magnitude is presently impossible with the current hardware technology. This chapter describes the current approaches for unsteady flow visualization. An effective particle tracing technique for unsteady flow is also described. First, the life cycle of a typical numerical flow simulation is outlined. Several unsteady flow data sets from real-world problems are then given. The current approaches for visualizing unsteady flow are then described. There are many existing systems for flow visualization, and some of them are discussed. Streaklines depict time-varying phenomena that are sometimes difficult or impossible to see with other visualization techniques. The algorithms for computing streaklines are described. Several unsteady flow data sets have been visualized using streaklines, and the results are presented. Finally, some current issues in unsteady flow visualization are discussed.

Lane, David A.↗

Progress in the Simulation of Steady and Time-Dependent Flows with 3D Parallel Unstructured Cartesian Methods

The proposed paper will present recent extensions in the development of an efficient Euler solver for adaptively-refined Cartesian meshes with embedded boundaries. The paper will focus on extensions of the basic method to include solution adaptation, time-dependent flow simulation, and arbitrary rigid domain motion. The parallel multilevel method makes use of on-the-fly parallel domain decomposition to achieve extremely good scalability on large numbers of processors, and is coupled with an automatic coarse mesh generation algorithm for efficient processing by a multigrid smoother. Numerical results are presented demonstrating parallel speed-ups of up to 435 on 512 processors. Solution-based adaptation may be keyed off truncation error estimates using tau-extrapolation or a variety of feature detection based refinement parameters. The multigrid method is extended to for time-dependent flows through the use of a dual-time approach. The extension to rigid domain motion uses an Arbitrary Lagrangian-Eulerlarian (ALE) formulation, and results will be presented for a variety of two- and three-dimensional example problems with both simple and complex geometry.

Aftosmis, M. J.↗

A New Time-dependent Model for the Martian Radiation Environment

Manned space activities have been until present time limited to the near-Earth environment, most of them to low Earth orbit (LEO) scenarios, with only some of the Apollo missions targeted to the Moon. In current times most human exploration and development of space (HEDS) activities are related to the development of the International Space Station (ISS), and therefore take place in the LEO environment. A natural extension of HEDS activities will be going beyond LEO, and reach asteroids, Mars, Jupiter, Saturn, the Kuiper belt and the outskirts of the Solar System. Such long journeys onboard spacecraft outside the protective umbrella of the geomagnetic field will require higher levels of protection from the radiation environment found in the deep space for both astronauts and equipment. So, it is important to have available a tool for radiation shielding which takes into account the radiation environments found all along the interplanetary space and at the different bodies encountered in the Solar System. Moreover, the radiation protection is one of the two NASA highest concerns and priorities. A tool integrating different radiation environments with shielding computation techniques especially tailored for deep space mission scenario is instrumental in view of this exigency. In view of manned missions targeted to Mars, for which radiation exposure is one of the greatest problems and challenges to be tackled, it is of fundamental importance to have available a tool which allows to know which are the particle flux and spectra at any time at any point of the Martian surface. With this goal in mind, a new model for the radiation environment to be found on the planet Mars due to Galactic Cosmic Rays (GCR) has been developed. Solar modulated primary particles rescaled for Mars conditions are transported within the Martian atmosphere, with temporal properties modeled with variable timescales, down to the surface, with altitude and backscattering patterns taken into account. The tool allows analysis for manned Mars landing missions, as well as planetary science studies, e.g. subsurface water and volatile inventory studies. This Mars environmental model is available through the SIREST website, a project of NASA Langley Research Center.

G De Angelis↗

Relationship Between Vortex Meander and Ambient Turbulence

Efforts are currently underway to increase the capacity of airports by use of closely-spaced parallel runways. If such an objective is to be achieved safely and efficiently during both visual and instrument flight conditions, it will be necessary to develop more precise methods for the prediction of the motion and spread of the hazard posed by the lift-generated vortex-wakes of aircraft, and their uncertainties. The purpose of the present study is to relate the motion induced in vortex filaments by turbulence in the ambient flow field to the measured turbulence in the flow field. The problem came about when observations made in the two largest NASA wind tunnels indicated that extended exposure of vortex wakes to the turbulence in the wind tunnel air stream causes the centers of the vortices to meander about with time at a given downstream station where wake measurements are being made. Although such a behavior was expected, the turbulence level based on the maximum amplitude of meander was much less than the root-mean-squared value measured in the free-stream of the wind tunnel by use of hot-film anemometers. An analysis of the time-dependent motion of segments of vortex filaments as they interact with an eddy, indicates that the inertia of the filaments retards their motion enough in the early part of their travel to account for a large part of the difference in the two determinations of turbulence level. Migration of vortex filaments from one turbulent eddy to another (probably with a different orientation), is believed to account for the remainder of the difference. Methods that may possibly be developed for use in the measurement of the magnitude of the more intense eddies in turbulent flow fields and how they should be adjusted to predict vortex meander are then discussed.

Rossow, Vernon J.↗

New Observations of the Heliospheric Magnetic Field from the Voyager Spacecraft

We review recent observations of variations of the heliospheric magnetic field B(t) made by Voyager 1 and 2 (V1 and V2), and we discuss the boundary conditions needed for models to explain the observations. Usually, observations from a spacecraft close to the Sun, such as ACE, WIND or Ulysses are used as input to a time-dependent model. Generally, the predicted profile B(t) can be compared directly with the observed profile only when either V1 or V2 is approximately radially aligned with a near-Sun spacecraft; this happens rarely and only for a brief time interval. The Bastille Day events illustrate this situation. In the absence of radial alignment of the spacecraft it is possible to predict the development of a global structure (a GMIR) with data from ACE or WIND, if they obtain a representative sample the flows that merge to form a GMIR. When latitudinal gradients are small and when there is statistical homogeneity in the azimuthal direction, it is possible to predict the statistical properties of the large-scale fluctuations of B(t) observed by V1 or V2 during a year or so. We illustrate this situation with observations from the recent solar maximum and the declining phase of the solar cycle. Predictions of detailed observations made by V1 and V2 under general conditions (e.g., when there is a large latitudinal gradient) require boundary conditions as a function of time on a surface, such as a Sun-centered sphere with a radius of 1 AU. These conditions can only be provided by global solar observations. We suggest the feasibility of such an approach, using V2 observations for 2005 and 2006. The prediction of observations in the heliosheath requires the solution of the 3-D boundary problem for the supersonic solar wind and propagation of solar wind through the termination shock into the heliosphere. The properties of B(t) observed in the heliosheath have not yet been predicted.

Burlaga, Leonard F.↗

Progress in Computational Simulation of Earthquakes

GeoFEST(P) is a computer program written for use in the QuakeSim project, which is devoted to development and improvement of means of computational simulation of earthquakes. GeoFEST(P) models interacting earthquake fault systems from the fault-nucleation to the tectonic scale. The development of GeoFEST( P) has involved coupling of two programs: GeoFEST and the Pyramid Adaptive Mesh Refinement Library. GeoFEST is a message-passing-interface-parallel code that utilizes a finite-element technique to simulate evolution of stress, fault slip, and plastic/elastic deformation in realistic materials like those of faulted regions of the crust of the Earth. The products of such simulations are synthetic observable time-dependent surface deformations on time scales from days to decades. Pyramid Adaptive Mesh Refinement Library is a software library that facilitates the generation of computational meshes for solving physical problems. In an application of GeoFEST(P), a computational grid can be dynamically adapted as stress grows on a fault. Simulations on workstations using a few tens of thousands of stress and displacement finite elements can now be expanded to multiple millions of elements with greater than 98-percent scaled efficiency on over many hundreds of parallel processors (see figure).

Donnellan, Andrea↗

Finite-Element Methods for Real-Time Simulation of Surgery

Two finite-element methods have been developed for mathematical modeling of the time-dependent behaviors of deformable objects and, more specifically, the mechanical responses of soft tissues and organs in contact with surgical tools. These methods may afford the computational efficiency needed to satisfy the requirement to obtain computational results in real time for simulating surgical procedures as described in Simulation System for Training in Laparoscopic Surgery (NPO-21192) on page 31 in this issue of NASA Tech Briefs. Simulation of the behavior of soft tissue in real time is a challenging problem because of the complexity of soft-tissue mechanics. The responses of soft tissues are characterized by nonlinearities and by spatial inhomogeneities and rate and time dependences of material properties. Finite-element methods seem promising for integrating these characteristics of tissues into computational models of organs, but they demand much central-processing-unit (CPU) time and memory, and the demand increases with the number of nodes and degrees of freedom in a given finite-element model. Hence, as finite-element models become more realistic, it becomes more difficult to compute solutions in real time. In both of the present methods, one uses approximate mathematical models trading some accuracy for computational efficiency and thereby increasing the feasibility of attaining real-time up36 NASA Tech Briefs, October 2003 date rates. The first of these methods is based on modal analysis. In this method, one reduces the number of differential equations by selecting only the most significant vibration modes of an object (typically, a suitable number of the lowest-frequency modes) for computing deformations of the object in response to applied forces.

Basdogan, Cagatay↗

Electro-Optic Time-to-Space Converter for Optical Detector Jitter Mitigation

A common problem in optical detection is determining the arrival time of a weak optical pulse that may comprise only one to a few photons. Currently, this problem is solved by using a photodetector to convert the optical signal to an electronic signal. The timing of the electrical signal is used to infer the timing of the optical pulse, but error is introduced by random delay between the absorption of the optical pulse and the creation of the electrical one. To eliminate this error, a time-to-space converter separates a sequence of optical pulses and sends them to different photodetectors, depending on their arrival time. The random delay, called jitter, is at least 20 picoseconds for the best detectors capable of detecting the weakest optical pulses, a single photon, and can be as great as 500 picoseconds. This limits the resolution with which the timing of the optical pulse can be measured. The time-to-space converter overcomes this limitation. Generally, the time-to-space converter imparts a time-dependent momentum shift to the incoming optical pulses, followed by an optical system that separates photons of different momenta. As an example, an electro-optic phase modulator can be used to apply longitudinal momentum changes (frequency changes) that vary in time, followed by an optical spectrometer (such as a diffraction grating), which separates photons with different momenta into different paths and directs them to impinge upon an array of photodetectors. The pulse arrival time is then inferred by measuring which photodetector receives the pulse. The use of a time-to-space converter mitigates detector jitter and improves the resolution with which the timing of an optical pulse is determined. Also, the application of the converter enables the demodulation of a pulse position modulated signal (PPM) at higher bandwidths than using previous photodetector technology. This allows the creation of a receiver for a communication system with high bandwidth and high bits/photon efficiency.

Birnbaum, Kevin↗