Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow conservation”

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 1,621 records · Page 90

Base Flow Investigation of the Apollo AS-202 Command Module

In recent years, both Europe and the US are developing hypersonic research and operational vehicles. These include (re)entry capsules (both ballistic and lifting) and lifting bodies such as ExoMars, EXPERT, ARV, CEV and IXV. The research programs are meant to enable technology and engineering capabilities to support during the next decade the development of affordable (possibly reusable) space transportation systems as well as hypersonic weapons systems for time critical targets. These programs have a broad range of goals, ranging from the qualification of thermal protection systems, the assessment of RCS performances, the development of GNC algorithms, to the full demonstration of the performance and operability of the integrated vehicles. Since the aerothermodynamic characteristics influence nearly all elements of the vehicle design, the accurate prediction of the aerothermal environment is a prerequisite for the design of efficient hypersonic systems. Significant uncertainties in the prediction of the hypersonic aerodynamic and the aerothermal loads can lead to conservative margins in the design of the vehicle including its Outer Mould Line (OML), thermal protection system, structure, and required control system robustness. The current level of aerothermal prediction uncertainties results therefore in reduced vehicle performances (e.g., sub-optimal payload to mass ratio, increased operational constraints). On the other hand, present computational capabilities enable the simulation of three dimensional flow fields with complex thermo-chemical models over complete trajectories and ease the validation of these tools by, e.g., reconstruction of detailed wind tunnel tests performed under identified and controlled conditions (flow properties and vehicle attitude in particular). These controlled conditions are typically difficult to achieve when performing in flight measurements which in turn results in large associated measurement uncertainties. Similar problems arise when attempting to rebuild measurements performed in "hot" ground facilities, where the difficulty level is increased by the addition of the free-flow characterization itself. The implementation of ever more sophisticated thermochemical models is no obvious cure to the aforementioned problems since their effect is often overwhelmed by the large measurement uncertainties incurred in both flight and ground high enthalpy facilities. Concurrent to the previous considerations, a major contributor to the overall vehicle mass of re-entry vehicles is the afterbody thermal protection system. This is due to the large acreage (equal or bigger than that of the forebody) to be protected. The present predictive capabilities for base flows are comparatively lower than those for windward flowfields and offer therefore a substantial potential for improving the design of future re-entry vehicles. To that end, it is essential to address the accuracy of high fidelity CFD tools exercised in the US and EU, which motivates a thorough investigation of the present status of hypersonic flight afterbody heating. This paper addresses the predictive capabilities of after body flow fields of re-entry vehicles investigated in the frame of the NATO/RTO - RTG-043 Task Group and is structured as follows: First, the verification of base flow topologies on the basis of available wind-tunnel results performed under controlled supersonic conditions (i.e., cold flows devoid of reactive effects) is performed. Such tests address the detailed characterization of the base flow with particular emphasis on separation/reattachment and their relation to Mach number effects. The tests have been performed on an Apollo-like re-entry capsule configuration. Second, the tools validated in the frame of the previous effort are exercised and appraised against flight-test data collected during the Apollo AS-202 re-entry.

Walpot, Louis M. G.↗

Towards quantifying mesoscale flows in the troposphere using Raman lidar and Sondes

Water vapor plays an important role in the energetics of the boundary layer processes which in turn play a key role in regulating regional and global climate. It plays a primary role in Earth's hydrological cycle, in radiation balance as a direct absorber of infrared radiation, and in atmospheric circulation as a latent heat energy source, as well as in determining cloud development and atmospheric stability. Water vapor concentration, expressed as a mass mixing ratio (g kg(exp -l)), is conserved in all meteorological processes except condensation and evaporation. This property makes it an ideal choice for studying many of the atmosphere's dynamic features. Raman scattering measurements from lidar also allow retrieval of water vapor mixing ratio profiles at high temporal and vertical resolution. Raman lidars sense water vapor to altitudes not achievable with towers and surface systems, sample the atmosphere at much higher temporal resolution than radiosondes or satellites, and do not require strong vertical gradients or turbulent fluctuations in temperature that is required by acoustic sounders and radars. Analysis of highly-resolved water vapor profiles are used here to characterize two important mesoscale flows: thunderstorm outflows and a cold front passage. The data were obtained at the Atmospheric Radiation Measurement Site (CART) by the groundbased Department of Energy/Sandia National Laboratories lidar (CART Raman lidar or CARL) and Goddard Space Flight Center Scanning Raman Lidar (SRL). A detailed discussion of the SRL and CARL performance during the IOPs is given by others in this meeting.

Demoz, B.↗

Multiscale modelling of evolving foams

We present a set of multi-scale interlinked algorithms to model the dynamics of evolving foams. These algorithms couple the key effects of macroscopic bubble rearrangement, thin film drainage, and membrane rupture. For each of the mechanisms, we construct consistent and accurate algorithms, and couple them together to work across the wide range of space and time scales that occur in foam dynamics. These algorithms include second order finite difference projection methods for computing incompressible fluid flow on the macroscale, second order finite element methods to solve thin film drainage equations in the lamellae and Plateau borders, multiphase Voronoi Implicit Interface Methods to track interconnected membrane boundaries and capture topological changes, and Lagrangian particle methods for conservative liquid redistribution during rearrangement and rupture. We derive a full set of numerical approximations that are coupled via interface jump conditions and flux boundary conditions, and show convergence for the individual mechanisms. We demonstrate our approach by computing a variety of foam dynamics, including coupled evolution of three-dimensional bubble clusters attached to an anchored membrane and collapse of a foam cluster.

97 MATHEMATICS AND COMPUTING↗

Numerical simulation of the geometrical-optics reduction of CE2 and comparisons to quasilinear dynamics

Zonal flows have been observed to appear spontaneously from turbulence in a number of physical settings. A complete theory for their behavior is still lacking. Recently, a number of studies have investigated the dynamics of zonal flows using quasilinear (QL) theories and the statistical framework of a second-order cumulant expansion (CE2). A geometrical-optics (GO) reduction of CE2, derived under an assumption of separation of scales between the fluctuations and the zonal flow, is studied here numerically. The reduced model, CE2-GO, has a similar phase-space mathematical structure to the traditional wave-kinetic equation, but that wave-kinetic equation has been shown to fail to preserve enstrophy conservation and to exhibit an ultraviolet catastrophe. CE2-GO, in contrast, preserves nonlinear conservation of both energy and enstrophy. We show here how to retain these conservation properties in a pseudospectral simulation of CE2-GO. We then present nonlinear simulations of CE2-GO and compare with direct simulations of quasilinear (QL) dynamics. Furthermore, we find that CE2-GO retains some similarities to QL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effective drift velocity from turbulent transport by vorticity

Here, we highlight the differing roles of vorticity and strain in the transport of coarse-grained scalars at length scales larger than ℓ by smaller scale (subscale) turbulence. We use the first term in a multiscale gradient expansion due to which exhibits excellent correlation with the exact subscale physics when the partitioning length ℓ is any scale smaller than that of the spectral peak. We show that unlike subscale strain, which acts as an anisotropic diffusion/antidiffusion tensor, subscale vorticity's contribution is solely a conservative advection of coarse-grained quantities by an eddy-induced nondivergent velocity, v * , that is proportional to the curl of vorticity. Therefore, material (Lagrangian) advection of coarse-grained quantities is accomplished not by the coarse-grained flow velocity, u¯ ℓ , but by the effective velocity, u¯ ℓ + v * , the physics of which may improve commonly used LES models.

42 ENGINEERING↗

A finite element model of conduction, convection, and phase change near a solid/melt interface

Detailed understanding of heat transfer and fluid flow is required for many aerospace thermal systems. These systems often include phase change and operate over a range of accelerations or effective gravitational fields. An approach to analyzing such systems is presented which requires the simultaneous solution of the conservation laws of energy, momentum, and mass, as well as an equation of state. The variable property form of the governing equations are developed in two-dimensional Cartesian coordinates for a Newtonian fluid. A numerical procedure for solving the governing equations is presented and implemented in a computer program. The Galerkin form of the finite element method is used to solve the spatial variation of the field variables, along with the implicit Crank-Nicolson time marching algorithm. Quadratic Langrangian elements are used for the internal energy and the two components of velocity. Linear Lagrangian elements are used for the pressure. The location of the solid/liquid interface as well as the temperatures are determined form the calculated internal energy and pressure. This approach is quite general in that it can describe heat transfer without phase change, phase change with a sharp interface, and phase change without an interface. Analytical results from this model are compared to those of other researchers studying transient conduction, convection, and phase change and are found to be in good agreement. The numerical procedure presented requires significant computer resources, but this is not unusual when compared to similar studies by other researchers. Several methods are suggested to reduce the computational times.

Viterna, Larry A.↗

Development of high-accuracy convection schemes for sequential solvers

An exploration is conducted of the applicability of such high resolution schemes as TVD to the resolving of sharp flow gradients using a sequential solution approach borrowed from pressure-based algorithms. It is shown that by extending these high-resolution shock-capturing schemes to a sequential solver that treats the equations as a collection of scalar conservation equations, the speed of signal propagation in the solution has to be coordinated by assigning the local convection speed as the characteristic speed for the entire system. A higher amount of dissipation is therefore needed to eliminate oscillations near discontinuities.

Thakur, Siddharth↗

Moist wind relationships

The equations describing the temporal and spatial behavior of the kinematic moisture and heat flux are described. In these nonlinear equations, the contribution by diabatic processes to the large-scale flux is composed of two parts. One part is associated with a Rayleigh damping term while the other arises from temporal and spatial changes in the pressure gradient term. It was found that the influence of diabatic processes on large-scale moisture fluxes depends greatly on the degree of balance between forcing and damping terms in the governing kinematic flux equations. The existence of a near balance requires a reduction in the large-scale horizontal geostrophic wind speed. Based on an examination of the moisture flux equations, it is argued that reductions in the large-scale horizontal wind speed observed within major cumulus cloud systems help conserve large-scale moisture fluxes. The deviation of the wind from geostrophic conditions is easily estimated for a near balanced state. This wind modification induces secondary vertical circulations that contribute to convergence, creating or supporting long-lived mesoscale flows. We believe this process to be a major supporter of the mesoscale circulations observed in severe storms and squall lines. In the tropics the wind modification has an antitriptic relationship. These diagnostic findings suggest possible modifications to the wind field in the application of a cumulus parameterization, and may be important in diabatic initialization of numerical weather prediction models.

Raymond, William H.↗

Computation of three-dimensional mixed convective boundary layer flow

The paper presents the numerical solution of heat and mass transfer during cross-flow (orthogonal) mixed convection. In this class of flow, a buoyancy-driven transport in the vertical direction and a forced convective flow in the horizontal direction results in a three-dimensional boundary layer structure adjacent to the plate. The rates of heat and mass transfer are determined by a combined influence of the two transport processes. The equations for the conservation of mass, momentum, energy, and species concentration were solved along with appropriate boundary conditions to determine the distributions of velocity components, temperature, and concentration across the thickness of the boundary layer at different locations on the plate. Results were expressed in dimensionless form using Reynolds number, Richardson number for heat transfer, Richardson number for mass transfer, Prandtl number, and Schmidt number as parameters. It was found that the transport is dominated by buoyancy at smaller vertical locations and at larger distances away from the forced convection leading edge. Effects of forced convection appeared to be very strong at smaller horizontal distances from the leading edge. The cross stream forced convection enhanced the rate of heat and mass transfer by a very significant amount.

Gadepalli, Prashandt↗

Aeroelastic Response and Protection of Space Shuttle External Tank Cable Trays

Sections of the Space Shuttle External Tank Liquid Oxygen (LO2) and Liquid Hydrogen (LH2) cable trays are shielded from potentially damaging airloads with foam Protuberance Aerodynamic Load (PAL) Ramps. Flight standard design LO2 and LH2 cable tray sections were tested with and without PAL Ramp models in the United States Air Force Arnold Engineering Development Center s (AEDC) 16T transonic wind tunnel to obtain experimental data on the aeroelastic stability and response characteristics of the trays and as part of the larger effort to determine whether the PAL ramps can be safely modified or removed. Computational Fluid Dynamic simulations of the full-stack shuttle launch configuration were used to investigate the flow characeristics around and under the cable trays without the protective PAL ramps and to define maximum crossflow Mach numbers and dynamic pressures experienced during launch. These crossflow conditions were used to establish wind tunnel test conditions which also included conservative margins. For all of the conditions and configurations tested, no aeroelastic instabilities or unacceptable dynamic response levels were encountered and no visible structural damage was experienced by any of the tested cable tray sections. Based upon this aeroelastic characterization test, three potentially acceptable alternatives are available for the LO2 cable tray PAL Ramps: Mini-Ramps, Tray Fences, or No Ramps. All configurations were tested to maximum conditions, except the LH2 trays at -15 deg. crossflow angle. This exception is the only caveat preventing the proposal of acceptable alternative configurations for the LH2 trays as well. Structural assessment of all tray loads and tray response measurements from launches following the Shuttle Return To Flight with the existing PAL Ramps will determine the acceptability of these PAL Ramp alternatives.

Edwards, John W.↗

Computational Fluid Dynamics Modeling of a Supersonic Nozzle and Integration into a Variable Cycle Engine Model

This paper covers the development of an integrated nonlinear dynamic simulation for a variable cycle turbofan engine and nozzle that can be integrated with an overall vehicle Aero-Propulso-Servo-Elastic (APSE) model. A previously developed variable cycle turbofan engine model is used for this study and is enhanced here to include variable guide vanes allowing for operation across the supersonic flight regime. The primary focus of this study is to improve the fidelity of the model's thrust response by replacing the simple choked flow equation convergent-divergent nozzle model with a MacCormack method based quasi-1D model. The dynamic response of the nozzle model using the MacCormack method is verified by comparing it against a model of the nozzle using the conservation element/solution element method. A methodology is also presented for the integration of the MacCormack nozzle model with the variable cycle engine.

Thrust↗

Computational Fluid Dynamics Modeling of a Supersonic Nozzle and Integration into a Variable Cycle Engine Model

This paper covers the development of an integrated nonlinear dynamic simulation for a variable cycle turbofan engine and nozzle that can be integrated with an overall vehicle Aero-Propulso-Servo-Elastic (APSE) model. A previously developed variable cycle turbofan engine model is used for this study and is enhanced here to include variable guide vanes allowing for operation across the supersonic flight regime. The primary focus of this study is to improve the fidelity of the model's thrust response by replacing the simple choked flow equation convergent-divergent nozzle model with a MacCormack method based quasi-1D model. The dynamic response of the nozzle model using the MacCormack method is verified by comparing it against a model of the nozzle using the conservation element/solution element method. A methodology is also presented for the integration of the MacCormack nozzle model with the variable cycle engine.

Control Theory↗

MAVEN Case Studies of Plasma Dynamics in Low‐Altitude Crustal Magnetic Field at Mars 1: Dayside Ion Spikes Associated With Radial Crustal Magnetic Fields

We report for the first time, simultaneous ion, electron, magnetic field vector and electric field wave measurements made possible by Mars Atmosphere and Volatile EvolutioN, during ion energy flux spikes in low‐altitude radial crustal magnetic fields on the Mars dayside. Observations show energetic electrons and ions (E > 25 eV) precipitating on magnetic field lines assumed as closed. Ions (E < 1.4 keV) display broad velocity distributions toward Mars, showing ions flowing from higher altitude possibly after magnetic reconnection or loss cone filling from pitch angle scattering effects. Precipitating ions (E < 1.4 keV) show nonadiabatic features depending on ion mass and energy and returning ions (E < 1.4 keV) show evidence of conserving the first adiabatic invariant in a mirror field. We observe magnetic field perturbations up to 60 nT, electric field wave amplitudes up to 38 mV/m, and brief periods of peaked electron spectra. At ∼175 km and at times Mars Atmosphere and Volatile EvolutioN is below the mirroring altitude of electrons, we observe mirroring and transverse heating of H+ ions alongside increased electric field wave amplitude fluctuations. It suggests field aligned potential drops result from different mirror altitudes of ions and electrons. Ions E > 1.4 keV (O+) occur as injected accelerated ion beams and ions heated after energization or deceleration. Energy dispersed kilo‐electron‐volt ions suggest a selection effect in radial magnetic fields for lower‐energy Marsward ions, compared to reflection of higher‐energy anti‐Sunward ions. Precipitating kilo‐electron‐volt ions show energy deposition rates of 3.6 ×10(exp -6) W/sq m and sputtering escape rates from precipitating O+ ions of 1.5 ×10(exp 5)/(sq cm.s) and 2.1 ×10(exp 6)/(sq cm.s) are calculated.

Plasma Dynamics↗

Comparing Commercial and Research Computational Fluid Dynamic Codes Using High-Order Workshop Benchmark Problems

The commercial computational fluid dynamic (CFD) code ANSYS Fluent and multiple research CFD codes (ez4d with uses the conservation element and solution element (CESE) method and codes that use the flux reconstruction (FR) method) were tested using three different benchmark problems from the International Workshop for High-Order CFD Methods. The benchmark problems included the transonic Ringleb flow, vortex transport by uniform flow, and laminar boundary layer on a flat plate. Simulation results from all three benchmark problems showed that the Fluent solutions had less error than the ez4d solutions for a given degree of freedom. As expected, both the Fluent and ez4d solutions had larger errors for a given degree of freedom than the simulations that used the FR method because both Fluent and ez4d utilized a second-order scheme whereas the FR codes utilized a fourth-order scheme.

CFD↗

Uncertainty Quantification on SAM Simulations of EBR-II Loss-of-Flow Tests

The System Analysis Module (SAM) is a modern system analysis tool being developed at Argonne National Laboratory for advanced non-LWR safety analysis. It aims to provide fast-running, whole-plant transient analyses capability with improved-fidelity for SFR, LFR, and MSR/FHR. In prior work, SAM has performed well in modeling the EBR-II benchmark based on the SHRT-17 and SHRT-45R protected and unprotected loss-of-flow experiments. However, underlying uncertainties in the modeling parameters need to be quantified to gain confidence in any best-estimate simulations. In this work, a Python coupling framework between Dakota and SAM has been developed to investigate these uncertainties, alongside code enhancements in SAM for user-flexibility in perturbing physical closure models and fluid equations of state. 25 significant input parameters covering boundary conditions, geometry, physical closure models, and equations of state were identified with their uncertainties defined from literature. Two SAM transient models for SHRT-17 and SHRT-45R were perturbed on these input parameters using Latin Hypercube Sampling and the uncertainties in predicted maximum fuel, cladding, and coolant temperatures were quantified as the output responses. Spearman rank correlation of input parameters to the predicted temperatures indicated that the most sensitive parameters include the initial core power, pipe wall drag, fuel & coolant conductivity, and peak channel flow area. Sensitivity of input parameters was found to differ between the two transients and core channel modeling strategies. The overall perturbations did not result in large uncertainties in the predicted temperatures, despite the efforts made to obtain accurate yet conservative input uncertainties.

Mui, Travis↗

A numerical study of latent thermal energy storage in a phase change material/carbon panel

To reduce the energetic dependence of building, it has become necessary to explore and develop new materials promoting energy conservation. Because of their high storage capacity, phase change materials (PCMs) are efficient to store thermal energy. In this paper, a 3D model was studied for simulation of energy storing cycles to predict the performances of PCM loaded panels. Carbon was used as supporting material for the PCM. The simulation was based on the enthalpy method using Ansys Fluent software. The panel was exposed to a daily heat flow including the effects of convection and radiation. The results show that the temperature decreased of approximately 2.5°C with a time shift about 2 hours. The steady state was reached after four cycles. Thus, after four cycles the PCM showed its effects on the temperature conditioning.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Loss Minimization With Optimal Power Dispatch in Multi-Frequency HVac Power Systems

Low-frequency high-voltage ac transmission scheme has recently been proposed as an alternative approach for bulk power transmission. This paper proposes a multi-period optimal power flow (OPF) for a multi-frequency HVac transmission system that interconnects both conventional 50/60-Hz and low-frequency grids using back-to-back converters with a centralized control scheme. The OPF objective is to minimize system losses by determining the optimal dispatch for generators, shunt capacitors, and converters. The OPF constraints include the operational constraints of all HVac grid and converter stations. The resulting mixed-integer nonlinear programing problem is solved using a proposed framework based on the predictor-corrector primal-dual interior-point method. The proposed OPF formulation and solution approach are verified using a multi-frequency HVac transmission system that is modified from the IEEE 57-bus system. Here, the results with the optimal dispatch from the proposed method during a simulated day show a significant loss reduction and an improved voltage regulation compared to those when an arbitrary dispatch is chosen.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗