Search NASA⌕ Search

SEARCH · Search NASA

Results for “initial conditions”

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 487 records · Page 27

Local-thermal-gradient and large-scale-circulation impacts on turbine-height wind speed forecasting over the Columbia River Basin

Abstract. We investigate the sensitivity of turbine-height wind speed forecast to initial condition (IC) uncertainties over the Columbia River Gorge (CRG) and Columbia River Basin (CRB) for two typical weather phenomena, i.e., local-thermal-gradient-induced marine air intrusion and a cold frontal passage. Four types of turbine-height wind forecast anomalies and their associated IC uncertainties related to local thermal gradients and large-scale circulations are identified using the self-organizing map (SOM) technique. The four SOM types are categorized into two patterns, each accounting for half of the ensemble members. The first pattern corresponds to IC uncertainties that alter the wind forecast through a modulating weather system, which produces the strongest wind anomalies in the CRG and CRB. In the second pattern, the moderate uncertainties in local thermal gradient and large-scale circulation jointly contribute to wind forecast anomaly. We analyze the cross section of wind and temperature anomalies through the gorge to explore the evolution of vertical features of each SOM type. The turbine-height wind anomalies induced by large-scale IC uncertainties are more concentrated near the front. In contrast, turbine-height wind anomalies induced by the local IC thermal uncertainties are found above the surface thermal anomalies. Moreover, the wind forecast accuracy in the CRG and CRB is limited by IC uncertainties in a few specific regions, e.g., the 2 m temperature within the basin and large-scale circulation over the northeast Pacific around 140∘ W.

17 WIND ENERGY↗

The structure of particle-laden jets and nonevaporating sprays

Mean and fluctuating gas velocities, liquid mass fluxes and drop sizes were in nonevaporating sprays. These results, as well as existing measurements in solid particle-laden jets, were used to evaluate models of these processes. The following models were considered: (1) a locally homogeneous flow (LHF) model, where slip between the phases was neglected; (2) a deterministic separated flow (DSF) model, where slip was considered but effects of particle dispersion by turbulence were ignored; and (3) a stochastic separated flow (SSF) model, where effects of interphase slip and turbulent dispersion were considered using random-walk computations for particle motion. The LHF and DSF models did not provide very satisfactory predictions over the present data base. In contrast, the SSF model performed reasonably well - including conditions in nonevaporating sprays where enhanced dispersion of particles by turbulence caused the spray to spread more rapidly than single-phase jets for comparable conditions. While these results are encouraging, uncertainties in initial conditions limit the reliability of the evaluation. Current work is seeking to eliminate this deficiency.

Shuen, J. S.↗

Numerical techniques for large cosmological N-body simulations

Techniques for carrying out large N-body simulations of the gravitational evolution of clustering in the fundamental cube of an infinite periodic universe are described and compared. The accuracy of the forces derived from several commonly used particle mesh schemes is examined, showing how submesh resolution can be achieved by including short-range forces between particles by direct summation techniques. The time integration of the equations of motion is discussed, and the accuracy of the codes for various choices of 'time' variable and time step is tested by considering energy conservation as well as by direct analysis of particle trajectories. Methods for generating initial particle positions and velocities corresponding to a growing mode representation of a specified power spectrum of linear density fluctuations are described. The effects of force resolution are studied and different simulation schemes are compared. An algorithm is implemented for generating initial conditions by varying the number of particles, the initial amplitude of density fluctuations, and the initial peculiar velocity field.

Efstathiou, G.↗

Numerical simulation of the compressible Orszag-Tang vortex 2. Supersonic flow

The numerical investigation of the Orszag-Tang vortex system in compressible magnetofluids will consider initial conditions with embedded supersonic regions. The simulations have initial average Mach numbers 1.0 and 1.5 and beta 10/3 with Lundquist numbers 50, 100, or 200. The behavior of the system differs significantly from that found previously for the incompressible and subsonic analogs. Shocks form at the downstream boundaries of the embedded supersonic regions outside the central magnetic X-point and produce strong local current sheets which dissipate appreciable magnetic energy. Reconnection at the central X-point, which dominates the incompressible and subsonic systems, peaks later and has a smaller impact as M increases from 0.6 to 1.5. Similarly, correlation between the momentum and magnetic field begins significant growth later than in subsonic and incompressible flows. The shocks bound large compression regions, which dominate the wavenumber spectra of autocorrelations in mass density, velocity, and magnetic field.

Picone, J. M.↗

Lifetimes of lunar satellite orbits

The Space Exploration Initiative has generated a renewed interest in lunar mission planning. The lunar missions currently under study, unlike the Apollo missions, involve long stay times. Several lunar gravity models have been formulated, but mission planners do not have enough confidence in the proposed models to conduct detailed studies of missions with long stay times. In this report, a particular lunar gravitational model, the Ferrari 5 x 5 model, was chosen to determine the lifetimes for 100-km and 300-km perilune altitude, near-circular parking orbits. The need to analyze orbital lifetimes for a large number of initial orbital parameters was the motivation for the formulation of a simplified gravitational model from the original model. Using this model, orbital lifetimes were found to be heavily dependent on the initial conditions of the nearly circular orbits, particularly the initial inclination and argument of perilune. This selected model yielded lifetime predictions of less than 40 days for some orbits, and other orbits had lifetimes exceeding a year. Although inconsistencies and limitations are inherent in all existing lunar gravity models, primarily because of a lack of information about the far side of the moon, the methods presented in this analysis are suitable for incorporating the moon's nonspherical gravitational effects on the preliminary design level for future lunar mission planning.

Meyer, Kurt W.↗

Quantification of Aerosol Transmission through Stress Corrosion Crack-Like Geometries

The formation of a stress corrosion crack (SCC) in the canister wall of a dry cask storage system (DCSS) has been identified as a potential issue for the long-term storage of spent nuclear fuel. The presence of an SCC in a storage system could represent a through-wall flow path from the canister interior to the environment. Modern, vertical DCSSs are of particular interest due to the commercial practice of using more significant backfill pressures in the canister, up to approximately 800 kPa. This pressure differential offers a relatively high driving potential for blowdown of any particulates that might be present in the canister. In this study, the rates of gas flow and aerosol transmission of a spent fuel surrogate through an engineered microchannel with dimensions representative of an SCC were evaluated experimentally using coupled mass flow and aerosol analyzers. The microchannel was formed by mating two gage blocks with a linearly tapering slot orifice nominally 13 μm (0.005 in.) tall on the upstream side and 25 μm (0.0010 in.) tall on the downstream side. The orifice is 12.7 mm (0.500 in.) wide by 8.86 mm (0.349 in.) long (flow length). Surrogate aerosols of cerium oxide, CeO 2 , were seeded and mixed with either helium or air inside a pressurized tank. The aerosol characteristics were measured immediately upstream and downstream of the simulated SCC at elevated and ambient pressures, respectively. These data sets are intended to add to previous testing that characterized SCCs under well-controlled boundary conditions through the inclusion of testing improvements that establish initial conditions in a more consistent way. These ongoing testing efforts are focused on understanding the evolution in both size and quantity of a hypothetical release of aerosolized spent fuel particles from failed fuel to the canister interior and ultimately through an SCC.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Preshot Report for Carbonara: Dynamics of a Boron Carbide Sandwich

The main purpose of the Carbonara experiment is to provide much-needed data to inform physics modeling of boron carbide (B 4 C) subjected to a sweeping high explosive (HE) detonation. The proposed experiment will include three separate shots, each employing a rectangular slab of HE (4” × 3” × 1”) sandwiched between two wedges of B 4 C, see Fig.1. In each case the top B 4 C wedge tapers from a height of 2.5 cm to 0.5 cm, while the bottom wedge tapers from 1.5 cm to 0.5 cm. The wedges are 4” in the direction of the proton beam, and 3” in the direction transverse to the beam. Three shots are currently planned, each one using a different HE: PBX 9501, PBX 9502, and PBX 9701. The HE will be initiated using a line wave generator (LWG) to preserve planar symmetry, in conjunction with a Comp. B booster to smooth out the profile for the initial condition on the main charge and to aid in initiation of the PBX 9502 shot. In each case the booster spans the entire length and height of the main HE slab and extends 0.5” in toward the LWG.

36 MATERIALS SCIENCE↗

Modeling the Effects of Microwave Heating on Densification in Chemical Vapor Infiltration

Microwave heating has great potential to accelerate the synthesis of ceramic matrix composites (CMCs) by a process called Chemical Vapor Infiltration (CVI). In CVI, reactive gases ingress a porous preform and undergo chemical transformation to deposit solid ceramic phase within the pores at high temperature thus, densifying the preform. However, the competing effects of chemical kinetics and gas transport are known to result in non-uniform densification as the outer surfaces of the preform experience faster depositions compared to the core. Achieving spatial temperature control plays a key role improving the densification quality. Microwave heating can potentially create temperature inversion such that the core of the preform is hotter than the outer surface and subsequently, lead to improved densification while keeping the manufacturing times and costs low. In the present work, a computational modeling strategy has been developed that accounts for the key physical phenomena responsible for densification of porous preforms using microwave heating. Specifically, a chemical kinetics model has been formulated for Silicon Carbide (SiC) deposition from MTS/H2 precursor. The model is implemented in a pore-resolved reactive transport solver, Quilt. The CVI simulations are performed for several conditions. Initially, parameterized temperature control is imposed to identify optimum conditions for good densification quality at a fraction of processing time. Further, simulations of microwave heating of porous SiC preforms are performed using OpenFOAM. A strategy to achieve and enhance temperature inversion is identified. The resulting temperature profiles are used in the pore-resolved densification simulations to analyze the densification behavior. It is observed that the temperature inversion achieved by microwave heating leads to significant improvements in densification quality and at the same time, keeps the manufacturing time low.

36 MATERIALS SCIENCE↗

Temperature distribution of a hot water storage tank in a simulated solar heating and cooling system

A 2,300-liter hot water storage tank was studied under conditions simulating a solar heating and cooling system. The initial condition of the tank, ranging from 37 C at the bottom to 94 C at the top, represented a condition midway through the start-up period of the system. During the five-day test period, the water in the tank gradually rose in temperature but in a manner that diminished its temperature stratification. Stratification was found not to be an important factor in the operation of the particular solar system studied.

Namkoong, D.↗

The Wiener-Hermite expansion applied to decaying isotropic turbulence using a renormalized time-dependent base

The problem of decaying isotropic turbulence has been studied using a Wiener-Hermite expansion with a renormalized time-dependent base. The theory is largely deductive and uses no modeling approximations. It has been found that many properties of large-Reynolds-number turbulence can be calculated (at least for moderate time) using the moving-base expansion alone. Such properties found are the spectrum shape in the dissipation range, the Kolmogorov constant, and the energy cascade in the inertial subrange. Furthermore, by using a renormalization scheme, it is possible to extend the calculation to larger times and to initial conditions significantly different from the equilibrium form. If the initial spectrum is the Kolmogorov spectrum perturbed with a spike or dip in the inertial subrange, the process proceeds to eliminate the perturbation and relax to the preferred spectrum shape. The turbulence decays with the proper dissipation rate, and several other properties are found to agree with measured data. The theory is also used to calculate the energy transfer and the flatness factor of turbulence.

Hogge, H. D.↗

Bridgman growth of semiconductors

The purpose of this study was to improve the understanding of the transport phenomena which occurs in the directional solidification of alloy semiconductors. In particular, emphasis was placed on the strong role of convection in the melt. Analytical solutions were not deemed possible for such an involved problem. Accordingly, a numerical model of the process was developed which simulated the transport. This translates into solving the partial differential equations of energy, mass, species, and momentum transfer subject to various boundary and initial conditions. A finite element method with simple elements was initially chosen. This simulation tool will enable the crystal grower to systematically identify and modify the important design factors within her control to produce better crystals.

Carlson, F. M.↗

Numerical models for the collapse and fragmentation of centrally condensed molecular cloud cores

The gravitational collapse and fragmentation of centrally condensed molecular cloud cores are investigated using a new hydrodynamical code. The numerical scheme is second-order accurate and uses explicit finite difference methods to advance the fluid variables on a 3D Cartesian grid. Two initial power-law density profiles, rho varies as r exp -1 and rho varies as r exp -2, are considered, as well as two initial density perturbations in the azimuthal coordinate theta, rho-i - rho(1 + a cos 2theta) where a = 0.1 and 0.5. Fragmentation is found to be possible in these centrally condensed cores if the initial conditions also include differential rotation. Models which collapse with initial uniform rotation do not produce fragments. If molecular cloud cores are indeed centrally condensed, as suggested by observations of star-forming regions and by studies of ambipolar diffusion, then differential rotation may be a mechanism for producing binary protostars during gravitational collapse.

Myhill, Elizabeth A.↗

Numerical simulation of torus-driven plasma transport in the Jovian magnetosphere

The Rice convection model has been modified for application to the transport of Io-generated plasma through the Jovian magnetosphere. The new code, called the RCM-J, has been used for several ideal-magnetohydrodynamic (MHD) numerical simulations to study how interchange instability causes an initially assumed torus configuration to break up. In simulations that start from a realistic torus configuration but include no energetic particles, the torus disintegrates too quickly (approximately 50 hours). By adding an impounding distribution of energetic particles to suppress the interchange instability, resonable lifetimes were obtained. For cases in which impoundment is insufficient to produce ideal-MHD stability, the torus breaks up predominantly into long fingers, unless the initial condition strongly favors some other geometrical form. If the initial torus has more mass on one side of the planet than the other, fingers form predominatly on the heavy side (which we associate with the active sector). Coriolis force bends the fingers to lag corotation. The simulation results are consistent with the idea that the fingers are formed with a longitudinal thickness that is roughly equal to the latitudinal distance over which the invariant density declines at the outer edges of the initial torus. Our calculations give an average longitudinal distance between plasma fingers of about 15 deg which corresponds to 20 to 30 minutes of rotation of the torus. We point to some Voyager and Ulysses data that are consistent with this scale of torus longitudinal irregularity.

Yang, Y. S.↗

Leading Your Leaders

Even though working on a problem has been your primary effort for the past year, your leadership may have heard about this once in a briefing a decade ago. Now they are basically clueless. Pretend that you are talking to your daughter's fifth-grade class. Explain how your complicated gizmo works. If possible, do not use acronyms. Define your terms. Put your work in context. Assume your leader has no idea what you do, who you work for, or what your gizmo does. That is a good place to start. Remember, taking the next century to study the problem or spending the Gross National Product to invent a new solution are probably not going to be acceptable solutions. Real engineers and technicians build real hardware that works in the real world in a reasonable manner within a reasonable time at a reasonable cost. True, skimping on time or money can cause mistakes, but folks whose gizmos are delayed unreasonably or cost more than is practical get their programs canceled, force the business into bankruptcy, or give the market over to the competition. Real engineers and technicians always consider cost and schedule in their work. Raising questions is important. However, we are in the business of doing things. Engineers and technicians are paid to get things done. Yes, you have to identify the problem, frame the design, identify the tests, perform the analysis, and assemble the hardware. But the goal is to solve the problem. Nobody ever said flying in space was easy. We make it look easy the same way that an Olympic champion makes her sport look easy: by working hard at improving performance every day. Better are the results of a well-defined test. Remember that a test on a laboratory bench is always an approximation of reality, and rules similar to those for good analysis also apply. One should always be mindful of Mechelay's rule: "It is better to be stupid than to run a stupid test." Often we try to overtest. If a piece of hardware passes an unbelievably difficult test, then life is good. More often when an unbelievably difficult test fails, we are left with a very long discussion of why and what was wrong in the design or execution of the test. Make sure that the test is well defined. Even then, it is important to explain to your leaders what inherent accuracy (or error) the test conditions or equipment have and what the assumptions or initial conditions were for the test. Test results without a good understanding of the test's accuracy or the pedigree of the test assumptions are worth very little. Finally, there is flight test data. Always limited, never at the edge of the envelope, it still shows how the real hardware works in a combined environment. Flight experience is dangerous because it typically doesn't show how close to the edge of the cliff the equipment is operating, but it does demonstrate how the hardware really works. A flight test is the ultimate test, again taken with the knowledge that it is probably not the extreme but something more like the middle of the environmental and systems performance. Good understanding of a problem and its solution always relies on a combination of all these methods. Be sure to lead your leaders by using all the tools you have at your disposal. At the end of the day, decisions in space flight always come down to a risk trade. Our business is not remotely safe, not in the sense that the public, the media, or our legislators use the term. Everything we do has a risk, cost, schedule, or performance trade-off. For your leaders to make an appropriate decision, you need to educate them, lead them, talk with them, and engage them in the discussion until full understanding takes place. It's your job. *

Hale, Wayne N.↗

Peak Seeking Control for Reduced Fuel Consumption with Preliminary Flight Test Results

The Environmentally Responsible Aviation project seeks to accomplish the simultaneous reduction of fuel burn, noise, and emissions. A project at NASA Dryden Flight Research Center is contributing to ERAs goals by exploring the practical application of real-time trim configuration optimization for enhanced performance and reduced fuel consumption. This peak-seeking control approach is based on Newton-Raphson algorithm using a time-varying Kalman filter to estimate the gradient of the performance function. In real-time operation, deflection of symmetric ailerons, trailing-edge flaps, and leading-edge flaps of a modified F-18 are directly optimized, and the horizontal stabilators and angle of attack are indirectly optimized. Preliminary results from three research flights are presented herein. The optimization system found a trim configuration that required approximately 3.5% less fuel flow than the baseline trim at the given flight condition. The algorithm consistently rediscovered the solution from several initial conditions. These preliminary results show the algorithm has good performance and is expected to show similar results at other flight conditions and aircraft configurations.

Brown, Nelson↗

Peak-Seeking Optimization of Trim for Reduced Fuel Consumption: Flight-Test Results

A peak-seeking control algorithm for real-time trim optimization for reduced fuel consumption has been developed by researchers at the National Aeronautics and Space Administration (NASA) Dryden Flight Research Center to address the goals of the NASA Environmentally Responsible Aviation project to reduce fuel burn and emissions. The peak-seeking control algorithm is based on a steepest-descent algorithm using a time-varying Kalman filter to estimate the gradient of a performance function of fuel flow versus control surface positions. In real-time operation, deflections of symmetric ailerons, trailing-edge flaps, and leading-edge flaps of an F/A-18 airplane (McDonnell Douglas, now The Boeing Company, Chicago, Illinois) are used for optimization of fuel flow. Results from six research flights are presented herein. The optimization algorithm found a trim configuration that required approximately 3 percent less fuel flow than the baseline trim at the same flight condition. The algorithm consistently rediscovered the solution from several initial conditions. These results show that the algorithm has good performance in a relevant environment.

flight optimization↗

Peak-Seeking Optimization of Trim for Reduced Fuel Consumption: Flight-test Results

A peak-seeking control algorithm for real-time trim optimization for reduced fuel consumption has been developed by researchers at the National Aeronautics and Space Administration (NASA) Dryden Flight Research Center to address the goals of the NASA Environmentally Responsible Aviation project to reduce fuel burn and emissions. The peak-seeking control algorithm is based on a steepest-descent algorithm using a time-varying Kalman filter to estimate the gradient of a performance function of fuel flow versus control surface positions. In real-time operation, deflections of symmetric ailerons, trailing-edge flaps, and leading-edge flaps of an F/A-18 airplane (McDonnell Douglas, now The Boeing Company, Chicago, Illinois) are used for optimization of fuel flow. Results from six research flights are presented herein. The optimization algorithm found a trim configuration that required approximately 3 percent less fuel flow than the baseline trim at the same flight condition. The algorithm consistently rediscovered the solution from several initial conditions. These results show that the algorithm has good performance in a relevant environment.

flight optimization↗

Hybrid Simulation of Solar-Wind-Like Turbulence

We present 2.5D hybrid simulations of the spectral and thermodynamic evolution of an initial state of magnetic field and plasma variables that in many ways represents solar wind fluctuations. In accordance with Helios near-Sun high-speed stream observations, we start with Alfvénic fluctuations along a mean magnetic field in which the fluctuations in the magnitude of the magnetic field are minimized. Since fluctuations in the radial flow speed are the dominant free energy in the observed fluctuations, we include a field-aligned v∥(k⊥) with an k(exp −1) spectrum of velocity fluctuations to drive the turbulent evolution. The flow rapidly distorts the Alfvénic fluctuations, yielding spectra (determined by spacecraft-like cuts) transverse to the field that become comparable to the k∥ fluctuations, as in spacecraft observations. The initial near constancy of the magnetic field is lost during the evolution; we show this also takes place observationally. We find some evolution in the anisotropy of the thermal fluctuations, consistent with expectations based on Helios data. We present 2D spectra of the fluctuations, showing the evolution of the power spectrum and cross-helicity. Despite simplifying assumptions, many aspects of simulations and observations agree. The greatly faster evolution in the simulations is at least in part due to the small scales being simulated, but also to the non-equilibrium initial conditions and the relatively low overall Alfvénicity of the initial fluctuations.

Roberts, D. Aaron↗