Results for “time integrator” — Page 21

NASA NTRS1979

Observability of the neutrino flux from the inner region of the galactic disk

The observability of galactic neutrinos in a detector of 10 to the 10th tons of water with an observing time of a few years is explored. Although the atmospheric flux exceeds the galactic flux considerably at energies greater than or equal to 1 TeV, the latter may still provide a marginally observable signal owing to its directionality. Galactic muon neutrinos with energy greater than approximately 1 TeV will produce a signal approximately two standard deviations above the atmospheric background over a four year period. If electron neutrinos can also be studied with DUMAND, then galactic electron neutrinos above 1 TeV would give an approximately four to five standard deviations signal above the electron neutrino background over a four-year integration time.

Silberberg, R.↗

NASA NTRS2002

RHESSI Investigations of the Neupert Effect in Solar Flares

In many solar flares the time integral of the hard X-ray emission closely matches the temporal variation of the soft X-ray emission. This is known as the Neupert Effect and is believed to result when the same accelerated electrons that produce the bremsstrahlung hard X-rays also heat the plasma that produces the thermal bremsstrahlung soft X-rays. The Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) is uniquely capable of investigating this effect since is covers the soft (down to 3 keV) and hard X-ray energy ranges with keV energy resolution, arcsecond-class angular resolution, and sub-second time resolution. Several GOES M-class flares have already been detected that show this general effect, but the detailed time histories suggest that there must be other heating mechanisms besides the thermalization of accelerated electrons that are heating the soft X-ray emitting plasma, even during the impulsive phase. Spatially-resolved hard and soft X-ray observations of these flares are also being investigated to determine what fraction of the heating is produced by the accelerated electrons.

Dennis, Brian R.↗

NASA NTRS2012

The Search for Habitable Worlds. 1. The Viability of a Starshade Mission

As part of NASA's mission to explore habitable planets orbiting nearby stars, this article explores the detection and characterization capabilities of a 4 m space telescope plus 50 m starshade located at the Earth-Sun L2 point, known as the New Worlds Observer (NWO). Our calculations include the true spectral types and distribution of stars on the sky, an iterative target selection protocol designed to maximize efficiency based on prior detections, and realistic mission constraints. We conduct simulated observing runs for a wide range in exozodiacal background levels (epsilon = 1-100 times the local zodi brightness) and overall prevalence of Earth-like terrestrial planets (eta(sub solar halo))0.1-1). We find that even without any return visits, the NWO baseline architecture (IWA = 65 mas, limiting FPB = 4 x 10(exp -11) can achieve a 95% probability of detecting and spectrally characterizing at least one habitable Earth-like planet and an expectation value of approximately 3 planets found, within the mission lifetime and delta V budgets, even in the worst-case scenario (eta(sub solar halo) = 0.1 and = epsilon = 100 zodis for every target). This achievement requires about 1 yr of integration time spread over the 5 yr mission, leaving the remainder of the telescope time for UV-NIR general astrophysics. Cost and technical feasibility considerations point to a "sweet spot" in starshade design near a 50 m starshade effective diameter. with 12 or 16 petals, at a distance of 70,000-100,000 km from the telescope.

Turnbull, Margaret C.↗

NASA NTRS1991

A fast implicit upwind solution algorithm for three-dimensional unstructured dynamic meshes

A fast implicit upwind algorithm for the solution of the time-dependent Euler equations is presented for aerodynamic analysis involving unstructured dynamic meshes. The spatial discretization of the scheme is based on the upwind approach of Roe, referred to as flux-difference splitting (FDS). The FDS approach is naturally dissipative and captures shock waves and contact discontinuities sharply. The temporal discretization of the scheme involves an implicit time-integration using a two-sweep Gauss-Seidel relaxation procedure. The procedure is computationally efficient for either steady or unsteady flow problems. A detailed description is given of the implicit upwind solution algorithm along with results which assess the capability. The results are presented for the NACA 0012 airfoil and for the Boeing 747 aircraft. The 747 geometry includes the fuselage, wing, horizontal and vertical tails, under-wing pylons, and flow-through engine nacelles. Euler solutions for the 747 aircraft on an unstructured tetrahedral mesh containing approximately 100,000 cells were obtained to engineering accuracy in less than one hour CPU time on a Cray-2 computer.

Batina, John T.↗

NASA NTRS1992

A fast implicit upwind solution algorithm for three-dimensional unstructured dynamic meshes

A fast implicit upwind algorithm for the solution of the time-dependent Euler equations is presented for aerodynamic analysis involving unstructured dynamic meshes. The spatial discretization of the scheme is based on the upwind approach of Roe, referred to as flux-difference splitting (FDS). The FDS approach is naturally dissipative and captures shock waves and contact discontinuities sharply. The temporal discretization of the scheme involves an implicit time-integration using a two-sweep Gauss-Seidel relaxation procedure. The procedure is computationally efficient for either steady or unsteady flow problems. A detailed description is given of the implicit upwind solution algorithm along with results which assess the capability. The results are presented for the NACA 0012 airfoil and for the Boeing 747 aircraft. The 747 geometry includes the fuselage, wing, horizontal and vertical tails, under-wing pylons, and flow-through engine nacelles. Euler solutions for the 747 aircraft on an unstructured tetrahedral mesh containing approximately 100,000 cells were obtained to engineering accuracy in less than one hour CPU time on a Cray-2 computer.

Batina, John T.↗

NASA NTRS1978

Observability of the neutrino flux from the inner region of the galactic disk

The observability of galactic neutrinos in a detector of 10 billion tons of water with an observing time of a few years is explored. Although the atmospheric flux exceeds the galactic flux considerably at energies greater than or equal to 1 TeV, the latter may still provide a marginally observable signal owing to its directionality. Galactic muon neutrinos with energy greater than or equal to 1 TeV will produce a signal approximately 2 sigma above the atmospheric background over a four year period. If electron neutrinos can also be studied with the deep underwater muon and neutrino detector, then galactic electron neutrinos above 1 TeV would give an approximate 4 to 5 sigma signal above the electron neutrino background over a four year integration time.

Silberberg, R.↗

NASA NTRS1987

Dynamic delamination buckling in composite laminates under impact loading: Computational simulation

A unique dynamic delamination buckling and delamination propagation analysis capability has been developed and incorporated into a finite element computer program. This capability consists of the following: (1) a modification of the direct time integration solution sequence which provides a new analysis algorithm that can be used to predict delamination buckling in a laminate subjected to dynamic loading, and (2) a new method of modeling the composite laminate using plate bending elements and multipoint constraints. This computer program is used to predict both impact induced buckling in composite laminates with initial delaminations and the strain energy release rate due to extension of the delamination. It is shown that delaminations near the outer surface of a laminate are susceptible to local buckling and buckling-induced delamination propagation when the laminate is subjected to transverse impact loading. The capability now exists to predict the time at which the onset of dynamic delamination buckling occurs, the dynamic buckling mode shape, and the dynamic delamination strain energy release rate.

Grady, Joseph E.↗

NASA NTRS1989

Dynamic delamination buckling in composite laminates under impact loading - Computational simulation

A unique dynamic delamination buckling and delamination propagation analysis capability has been developed and incorporated into a finite element computer program. This capability consists of the following: (1) a modification of the direct time integration solution sequence which provides a new analysis algorithm that can be used to predict delamination buckling in a laminate subjected to dynamic loading, and (2) a new method of modeling the composite laminate using plate bending elements and multipoint constraints. This computer program is used to predict both impact induced buckling in composite laminates with initial delaminations and the strain energy release rate due to extension of the delamination. It is shown that delaminations near the outer surface of a laminate are susceptible to local buckling and buckling-induced delamination propagation when the laminate is subjected to transverse impact loading. The capability now exists to predict the time at which the onset of dynamic delamination buckling occurs, the dynamic buckling mode shape, and the dynamic delamination strain energy release rate.

Grady, Joseph E.↗

NASA NTRS2007

Synthetic Vision Systems - Operational Considerations Simulation Experiment

Synthetic vision is a computer-generated image of the external scene topography that is generated from aircraft attitude, high-precision navigation information, and data of the terrain, obstacles, cultural features, and other required flight information. A synthetic vision system (SVS) enhances this basic functionality with real-time integrity to ensure the validity of the databases, perform obstacle detection and independent navigation accuracy verification, and provide traffic surveillance. Over the last five years, NASA and its industry partners have developed and deployed SVS technologies for commercial, business, and general aviation aircraft which have been shown to provide significant improvements in terrain awareness and reductions in the potential for Controlled-Flight-Into-Terrain incidents/accidents compared to current generation cockpit technologies. It has been hypothesized that SVS displays can greatly improve the safety and operational flexibility of flight in Instrument Meteorological Conditions (IMC) to a level comparable to clear-day Visual Meteorological Conditions (VMC), regardless of actual weather conditions or time of day. An experiment was conducted to evaluate SVS and SVS-related technologies as well as the influence of where the information is provided to the pilot (e.g., on a Head-Up or Head-Down Display) for consideration in defining landing minima based upon aircraft and airport equipage. The "operational considerations" evaluated under this effort included reduced visibility, decision altitudes, and airport equipage requirements, such as approach lighting systems, for SVS-equipped aircraft. Subjective results from the present study suggest that synthetic vision imagery on both head-up and head-down displays may offer benefits in situation awareness; workload; and approach and landing performance in the visibility levels, approach lighting systems, and decision altitudes tested.

Kramer, Lynda J.↗

NASA NTRS2024

Discounting Water for Optimal Carbon Gain as a Basis of Stomatal Closure

The exchange of carbon dioxide and water vapor between terrestrial ecosystems and the atmosphere is regulated by stomata (small pores in the leaves of plants). Unsurprisingly, environmental factors controlling the opening and closure of stomata has been sought as early as 1800. One approach, popularized in the early 1970s, is a stomatal optimization framework. This framework is based on the hypothesis that plants optimize carbon gain subject to water loss or water availability constraints. This constraint optimization problem was solved in various forms assuming instantaneous adjustments of stomatal aperture to maximize a reward function with no future foresight or legacy effects. Holtzman et al. (2024, https://doi.org/10.1029/2023av001113) offers a novel approach that can diagnose the effective timescale over which the reward function maximization must be time-integrated. The developed method thus optimizes an integrated carbon gain function but adjusted by a discount factor subject to water availability in the root zone. The discount factor considers how the plant values carbon gain to save water and its timescale can be inferred from observations because the model is analytically tractable. The results suggest that the most important climate factor that determines this discount timescale is multi-annual mean of the longest dry period during the growing season. The findings highlight how local climate traits influence the spatial variation in ecosystem-level water use strategies. This sets the stage for expanding such a framework to cases where multiple constraints act in concert while operating at distinct time scales.

stomata↗

NASA NTRS1990

Conical Euler simulation and active suppression of delta wing rocking motion

A conical Euler code was developed to study unsteady vortex-dominated flows about rolling highly-swept delta wings, undergoing either forced or free-to-roll motions including active roll suppression. The flow solver of the code involves a multistage Runge-Kutta time-stepping scheme which uses a finite volume spatial discretization of the Euler equations on an unstructured grid of triangles. The code allows for the additional analysis of the free-to-roll case, by including the rigid-body equation of motion for its simultaneous time integration with the governing flow equations. Results are presented for a 75 deg swept sharp leading edge delta wing at a freestream Mach number of 1.2 and at alpha equal to 10 and 30 deg angle of attack. A forced harmonic analysis indicates that the rolling moment coefficient provides: (1) a positive damping at the lower angle of attack equal to 10 deg, which is verified in a free-to-roll calculation; (2) a negative damping at the higher angle of attack equal to 30 deg at the small roll amplitudes. A free-to-roll calculation for the latter case produces an initially divergent response, but as the amplitude of motion grows with time, the response transitions to a wing-rock type of limit cycle oscillation. The wing rocking motion may be actively suppressed, however, through the use of a rate-feedback control law and antisymmetrically deflected leading edge flaps. The descriptions of the conical Euler flow solver and the free-to-roll analysis are presented. Results are also presented which give insight into the flow physics associated with unsteady vortical flows about forced and free-to-roll delta wings, including the active roll suppression of this wing-rock phenomenon.

Lee, Elizabeth M.↗

NASA NTRS1996

Parallel Simulation of Unsteady Turbulent Flames

Time-accurate simulation of turbulent flames in high Reynolds number flows is a challenging task since both fluid dynamics and combustion must be modeled accurately. To numerically simulate this phenomenon, very large computer resources (both time and memory) are required. Although current vector supercomputers are capable of providing adequate resources for simulations of this nature, the high cost and their limited availability, makes practical use of such machines less than satisfactory. At the same time, the explicit time integration algorithms used in unsteady flow simulations often possess a very high degree of parallelism, making them very amenable to efficient implementation on large-scale parallel computers. Under these circumstances, distributed memory parallel computers offer an excellent near-term solution for greatly increased computational speed and memory, at a cost that may render the unsteady simulations of the type discussed above more feasible and affordable.This paper discusses the study of unsteady turbulent flames using a simulation algorithm that is capable of retaining high parallel efficiency on distributed memory parallel architectures. Numerical studies are carried out using large-eddy simulation (LES). In LES, the scales larger than the grid are computed using a time- and space-accurate scheme, while the unresolved small scales are modeled using eddy viscosity based subgrid models. This is acceptable for the moment/energy closure since the small scales primarily provide a dissipative mechanism for the energy transferred from the large scales. However, for combustion to occur, the species must first undergo mixing at the small scales and then come into molecular contact. Therefore, global models cannot be used. Recently, a new model for turbulent combustion was developed, in which the combustion is modeled, within the subgrid (small-scales) using a methodology that simulates the mixing and the molecular transport and the chemical kinetics within each LES grid cell. Finite-rate kinetics can be included without any closure and this approach actually provides a means to predict the turbulent rates and the turbulent flame speed. The subgrid combustion model requires resolution of the local time scales associated with small-scale mixing, molecular diffusion and chemical kinetics and, therefore, within each grid cell, a significant amount of computations must be carried out before the large-scale (LES resolved) effects are incorporated. Therefore, this approach is uniquely suited for parallel processing and has been implemented on various systems such as: Intel Paragon, IBM SP-2, Cray T3D and SGI Power Challenge (PC) using the system independent Message Passing Interface (MPI) compiler. In this paper, timing data on these machines is reported along with some characteristic results.

Menon, Suresh↗

NASA NTRS2016

Parallel Aircraft Trajectory Optimization with Analytic Derivatives

Trajectory optimization is an integral component for the design of aerospace vehicles, but emerging aircraft technologies have introduced new demands on trajectory analysis that current tools are not well suited to address. Designing aircraft with technologies such as hybrid electric propulsion and morphing wings requires consideration of the operational behavior as well as the physical design characteristics of the aircraft. The addition of operational variables can dramatically increase the number of design variables which motivates the use of gradient based optimization with analytic derivatives to solve the larger optimization problems. In this work we develop an aircraft trajectory analysis tool using a Legendre-Gauss-Lobatto based collocation scheme, providing analytic derivatives via the OpenMDAO multidisciplinary optimization framework. This collocation method uses an implicit time integration scheme that provides a high degree of sparsity and thus several potential options for parallelization. The performance of the new implementation was investigated via a series of single and multi-trajectory optimizations using a combination of parallel computing and constraint aggregation. The computational performance results show that in order to take full advantage of the sparsity in the problem it is vital to parallelize both the non-linear analysis evaluations and the derivative computations themselves. The constraint aggregation results showed a significant numerical challenge due to difficulty in achieving tight convergence tolerances. Overall, the results demonstrate the value of applying analytic derivatives to trajectory optimization problems and lay the foundation for future application of this collocation based method to the design of aircraft with where operational scheduling of technologies is key to achieving good performance.

aircraft↗

NASA NTRS1993

Wing flutter boundary prediction using an unsteady Euler aerodynamic method

Modifications to an existing three-dimensional, implicit, upwind Euler/Navier-Stokes code (CFL3D Version 2.1) for the aeroelastic analysis of wings are described. These modifications, which were previously added to CFL3D Version 1.0, include the incorporation of a deforming mesh algorithm and the addition of the structural equations of motion for their simultaneous time-integration with the government flow equations. The paper gives a brief description of these modifications and presents unsteady calculations which check the modifications to the code. Euler flutter results for an isolated 45 degree swept-back wing are compared with experimental data for seven freestream Mach numbers which define the flutter boundary over a range of Mach number from 0.499 to 1.14. These comparisons show good agreement in flutter characteristics for freestream Mach numbers below unity. For freestream Mach numbers above unity, the computed aeroelastic results predict a premature rise in the flutter boundary as compared with the experimental boundary. Steady and unsteady contours of surface Mach number and pressure are included to illustrate the basic flow characteristics of the time-marching flutter calculations and to aid in identifying possible causes for the premature rise in the computational flutter boundary.

Lee-Rausch, Elizabeth M.↗

NASA NTRS1993

Wing flutter boundary prediction using unsteady Euler aerodynamic method

Modifications to an existing 3D implicit upwind Euler/Navier-Stokes code for the aeroelastic analysis of wings are described. These modifications include the incorporation of a deforming mesh algorithm and the addition of the structural equations of motion for their simultaneous time-integration with the governing flow equations. The paper gives a brief description of these modifications and presents unsteady calculations which check the modifications to the code. Euler flutter results for an isolated 45 deg swept-back wing are compared with experimental data for seven freestream Mach numbers which define the flutter boundary over a range of Mach number from 0.499 to 1.14. These comparisons show good agreement in flutter characteristics for freestream Mach numbers below unity. For freestream Mach numbers above unity, the computed aeroelastic results predict a premature rise in the flutter boundary as compared with the experimental boundary. Steady and unsteady contours of surface Mach number and pressure are included to illustrate the basic flow characteristics of the time-marching flutter calculations and to aid in identifying possible causes for the premature rise in the computational flutter boundary.

Lee-Rausch, Elizabeth M.↗

NASA NTRS2006

Fundamental Limits to Wavefront Sensing in the Submillimeter

With the advent of large-format submillimeter wavelength detector arrays, and a new 25 m diameter submillimeter telescope under consideration, the question of optimal wavefront sensing methods is timely. Indeed, not only should bolometric array detectors allow the use of a variety of wavefront sensing techniques already in use in the optical/infrared, but in some cases it should actually be easier to apply these techniques because of the more benign temporal properties of the atmosphere at long wavelengths. This paper thus addresses the fundamental limits to wavefront sensing at submillimeter wavelengths, in order to determine how well a telescope surface can be measured in the submillimeter band. First several potential measurement approaches are discussed and compared. Next the theoretical accuracy of a fringe phase measurement in the submillimeter is discussed. It is concluded that with Mars as the source, wavefront sensing at the micron level should be achievable at submillimeter wavelengths in quite reasonable integration times.

wavefront sensing↗

NASA NTRS1991

Fourier functional analysis for unsteady aerodynamic modeling

A method based on Fourier analysis is developed to analyze the force and moment data obtained in large amplitude forced oscillation tests at high angles of attack. The aerodynamic models for normal force, lift, drag, and pitching moment coefficients are built up from a set of aerodynamic responses to harmonic motions at different frequencies. Based on the aerodynamic models of harmonic data, the indicial responses are formed. The final expressions for the models involve time integrals of the indicial type advocated by Tobak and Schiff. Results from linear two- and three-dimensional unsteady aerodynamic theories as well as test data for a 70-degree delta wing are used to verify the models. It is shown that the present modeling method is accurate in producing the aerodynamic responses to harmonic motions and the ramp type motions. The model also produces correct trend for a 70-degree delta wing in harmonic motion with different mean angles-of-attack. However, the current model cannot be used to extrapolate data to higher angles-of-attack than that of the harmonic motions which form the aerodynamic model. For linear ramp motions, a special method is used to calculate the corresponding frequency and phase angle at a given time. The calculated results from modeling show a higher lift peak for linear ramp motion than for harmonic ramp motion. The current model also shows reasonably good results for the lift responses at different angles of attack.

Lan, C. Edward↗

NASA NTRS1991

Fourier functional analysis for unsteady aerodynamic modeling

A method based on Fourier analysis is developed to analyze the force and moment data obtained in large amplitude forced oscillation tests at high angles of attack. The aerodynamic models for normal force, lift, drag, and pitching moment coefficients are built up from a set of aerodynamic responses to harmonic motions at different frequencies. Based on the aerodynamic models of harmonic data, the indicial responses are formed. The final expressions for the models involve time integrals of the indicial type advocated by Tobak and Schiff. Results from linear two- and three-dimensional unsteady aerodynamic theories as well as test data for a 70-degree delta wing are used to verify the models. It is shown that the present modeling method is accurate in producing the aerodynamic responses to harmonic motions and the ramp type motions. The model also produces correct trend for a 70-degree delta wing in harmonic motion with different mean angles-of-attack. However, the current model cannot be used to extrapolate data to higher angles-of-attack than that of the harmonic motions which form the aerodynamic model. For linear ramp motions, a special method is used to calculate the corresponding frequency and phase angle at a given time. The calculated results from modeling show a higher lift peak for linear ramp motion than for harmonic ramp motion. The current model also shows reasonably good results for the lift responses at different angles of attack.

Chin, Suei↗