SEARCH · Search NASA

Results for “boundary value 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,423 records · Page 79

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Lattice determination of 𝐼=0 and 2 𝜋⁢𝜋 scattering phase shifts with a physical pion mass

Phase shifts for 𝑠-wave 𝜋⁢𝜋 scattering in both the 𝐼 =0 and 𝐼 =2 channels are determined from a lattice QCD calculation performed on 741 gauge configurations obeying G-parity boundary conditions with a physical pion mass and lattice size of 32 3 ×64. These results support our recent 2021 study of direct 𝐶⁢𝑃 violation in 𝐾 → 𝜋⁢𝜋 decay, improving our earlier 2015 calculation. The phase shifts are determined for both stationary and moving 𝜋⁢𝜋 systems, at three (𝐼 = 0) and four (𝐼 = 2) different total momenta. We implement several 𝜋⁢𝜋 interpolating operators including a scalar bilinear “𝜎” operator and paired single-pion bilinear operators with the constituent pions carrying various relative momenta. Several techniques, including correlated fitting and a bootstrap determination of p-values have been used to refine the results and a comparison with the generalized eigenvalue problem method is given. A detailed systematic error analysis is performed which allows phase shift results to be presented at a fixed energy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1983

Spherical accretion onto quasars

The problem of time-steady accretion of gas irradiated by a self-consistently generated quasar-like continuum is investigated. A controversy in the literature is resolved by a careful distinction between the existence of solutions and their stability. Specific criteria are found for when thermal instabilities may disrupt single-phase steady flow. As a function of initial temperature (1,000,000-100,000,000 K), the boundaries of stable flow in the luminosity are determined relative to the Eddington (L/LE)-efficiency plane. The thermal and dynamical history of these steady flows is analyzed. The observational data that define quasars and type I Seyfert galaxies are matched to find values of efficiency and L/LE required for consistent spherically symmetric steady-flow models. The efficiency is constrained to be less than 1/100, and L/LE resides between several times 1/1000-1/100.

Krolik, J. H.↗

NASA NTRS1996

Kinetics of Diffusional Droplet Growth in a Liquid/Liquid Two-Phase System

We address the problem of diffusional interactions in a finite sized cluster of spherical particles for volume fractions, V(sub v) in the range 0-0.01. We determined the quasi-static monopole diffusion solution for n particles distributed at random in a continuous matrix. A global mass conservation condition is employed, obviating the need for any external boundary condition. The numerical results provide the instantaneous (snapshot) growth or shrinkage rate of each particle, precluding the need for extensive time-dependent computations. The close connection between these snapshot results and the coarsegrained kinetic constants are discussed. A square-root dependence of the deviations of the rate constants from their zero volume fraction value is found for the higher V(sub v) investigated. This behavior is consistent with predictions from diffusion Debye-Huckel screening theory. By contrast, a cube-root dependence, reported in earlier numerical studies, is found for the lower V(sub v) investigated. The roll-over region of the volume fraction where the two asymptotics merge depends on the number of particles, n, alone. A theoretical estimate for the roll-over point predicts that the corresponding V(sub v) varies as n(sup -2), in good agreement with the numerical results.

Glicksman, M. E.↗

NASA NTRS2004

On Nonlinear Combustion Instability in Liquid Propellant Rocket Motors

All liquid propellant rocket instability calculations in current use have limited value in the predictive sense and serve mainly as a correlating framework for the available data sets. The well-known n-t model first introduced by Crocco and Cheng in 1956 is still used as the primary analytical tool of this type. A multitude of attempts to establish practical analytical methods have achieved only limited success. These methods usually produce only stability boundary maps that are of little use in making critical design decisions in new motor development programs. Recent progress in understanding the mechanisms of combustion instability in solid propellant rockets"' provides a firm foundation for a new approach to prediction, diagnosis, and correction of the closely related problems in liquid motor instability. For predictive tools to be useful in the motor design process, they must have the capability to accurately determine: 1) time evolution of the pressure oscillations and limit amplitude, 2) critical triggering pulse amplitude, and 3) unsteady heat transfer rates at injector surfaces and chamber walls. The method described in this paper relates these critical motor characteristics directly to system design parameters. Inclusion of mechanisms such as wave steepening, vorticity production and transport, and unsteady detonation wave phenomena greatly enhance the representation of key features of motor chamber oscillatory behavior. The basic theoretical model is described and preliminary computations are compared to experimental data. A plan to develop the new predictive method into a comprehensive analysis tool is also described.

Sims, J. D.↗

DOE OSTI2017

Coulomb Problem for Z > Z{sub cr} in Doped Graphene

The dynamics of charge carriers in doped graphene, i.e., graphene with a gap in the energy spectrum depending on the substrate, in the presence of a Coulomb impurity with charge Z is considered within the effective two-dimensional Dirac equation. The wave functions of carriers with conserved angular momentum J = M + 1/2 are determined for a Coulomb potential modified at small distances. This case, just as any two-dimensional physical system, admits both integer and half-integer quantization of the orbital angular momentum in plane, M = 0, ±1, ±2, …. For J = 0, ±1/2, ±1, critical values of the effective charge Z{sub cr}(J, n) are calculated for which a level with angular momentum J and radial quantum numbers n = 0 and n = 1 reaches the upper boundary of the valence band. For Z < Z{sub cr} (J, n = 0), the energy of a level is presented as a function of charge Z for the lowest values of orbital angular momentum M, the level with J = 0 being the first to descend to the band edge. For Z>Z{sub cr} (J, n = 0), scattering phases are calculated as a function of hole energy for several values of supercriticality, as well as the positions ε{sub 0} and widths γ of quasistationary states as a function of supercriticality. The values of ε{sub 0}* and width γ* are pointed out for which quasidiscrete levels may show up as Breit–Wigner resonances in the scattering of holes by a supercritical impurity. Since the phases are real, the partial scattering matrix is unitary, so that the radial Dirac equation is consistent even for Z > Z{sub cr}. In this single-particle approximation, there is no spontaneous creation of electron–hole pairs, and the impurity charge cannot be screened by this mechanism.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1989

Automatic control of cryogenic wind tunnels

Inadequate Reynolds number similarity in testing of scaled models affects the quality of aerodynamic data from wind tunnels. This is due to scale effects of boundary-layer shock wave interaction which is likely to be severe at transonic speeds. The idea of operation of wind tunnels using test gas cooled to cryogenic temperatures has yielded a quantrum jump in the ability to realize full scale Reynolds number flow similarity in small transonic tunnels. In such tunnels, the basic flow control problem consists of obtaining and maintaining the desired test section flow parameters. Mach number, Reynolds number, and dynamic pressure are the three flow parameters that are usually required to be kept constant during the period of model aerodynamic data acquisition. The series of activity involved in modeling, control law development, mechanization of the control laws on a microcomputer, and the performance of a globally stable automatic control system for the 0.3-m Transonic Cryogenic Tunnel (TCT) are discussed. A lumped multi-variable nonlinear dynamic model of the cryogenic tunnel, generation of a set of linear control laws for small perturbation, and nonlinear control strategy for large set point changes including tunnel trajectory control are described. The details of mechanization of the control laws on a 16 bit microcomputer system, the software features, operator interface, the display and safety are discussed. The controller is shown to provide globally stable and reliable temperature control to + or - 0.2 K, pressure to + or - 0.07 psi and Mach number to + or - 0.002 of the set point value. This performance is obtained both during large set point commands as for a tunnel cooldown, and during aerodynamic data acquisition with intrusive activity like geometrical changes in the test section such as angle of attack changes, drag rake movements, wall adaptation and sidewall boundary-layer removal. Feasibility of the use of an automatic Reynolds number control mode with fixed Mach number control is demonstrated.

Balakrishna, S.↗

DOE OSTI2017

Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part I

In this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes—wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries—and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid–structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier–Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau–Rayleigh instability in water ripple free surface flow. In conclusion, these last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part II

Here, in this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes—wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries—and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid–structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier–Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau–Rayleigh instability in water ripple free surface flow. In these last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1991

The effects of suction on the nonlinear stability of the three-dimensional boundary layer above a rotating disc

There exist two types of stationary instability of the flow over a rotating disc corresponding to the upper branch, inviscid mode and the lower branch mode, which has a triple deck structure, of the neutral stability curve. A theoretical study of the linear problem and an account of the weakly nonlinear properties of the lower branch modes have been undertaken by Hall and MacKerrell respectively. Motivated by recent reports of experimental sightings of the lower branch mode and an examination of the role of suction on the linear stability properties of the flow here, the effects are studied of suction on the nonlinear disturbance described by MacKerrell. The additional analysis required in order to incorporate suction is relatively straightforward and enables the derivation of an amplitude equation which describes the evolution of the mode. For each value of the suction, a threshold value of the disturbance amplitude is obtained; modes of size greater than this threshold grow without limit as they develop away from the point of neutral stability.

Bassom, Andrew P.↗

NASA NTRS1975

Nonlinear vibration of cylindrical shells

The large amplitude vibration of a thin-walled cylindrical shell is analyzed. A perturbation method is used to solve the steady-state forced vibration problem. The simply-supported boundary conditions and the circumferential periodicity condition are satisfied. The resulting solution indicates that in addition to the fundamental mode, the response contains asymmetric modes as well as axisymmetric modes with the frequency twice that of the fundamental mode. Vibrations involving a single driven mode response are investigated. The results indicate that the nonlinearity is either softening or hardening depending on the mode. The vibrations involving both a driven mode and a companion mode are also investigated. An experimental investigation is conducted, and the results are in qualitative agreement with the theory. 'Nonstationary' response is detected at some frequencies for large amplitude response in which the amplitude drifts from one value to another. Various nonlinear phenomena are observed and comparisons with the theoretical results are made.

Chen, J. C.↗

NASA NTRS1999

Electrical Impact of SiC Structural Crystal Defects on High Electric Field Devices

Commercial epilayers are known to contain a variety of crystallographic imperfections. including micropipes, closed core screw dislocations. low-angle boundaries, basal plane dislocations, heteropolytypic inclusions, and non-ideal surface features like step bunching and pits. This paper reviews the limited present understanding of the operational impact of various crystal defects on SiC electrical devices. Aside from micropipes and triangular inclusions whose densities have been shrinking towards manageably small values in recent years, many of these defects appear to have little adverse operational and/or yield impact on SiC-based sensors, high-frequency RF, and signal conditioning electronics. However high-power switching devices used in power management and distribution circuits have historically (in silicon experience) demanded the highest material quality for prolonged safe operation, and are thus more susceptible to operational reliability problems that arise from electrical property nonuniformities likely to occur at extended crystal defects. A particular emphasis is placed on the impact of closed-core screw dislocations on high-power switching devices, because these difficult to observe defects are present in densities of thousands per cm,in commercial SiC epilayers. and their reduction to acceptable levels seems the most problematic at the present time.

Neudeck, Philip G.↗

DOE OSTI2016

On Some Criticality Benchmarks in Spherical Geometry

In a paper published in 2003, Sood et al. reported numerical results for a set of problems that are useful as benchmarks for criticality code verification. Here, we reconsider a subset of those problems, namely the one-group problems for homogeneous spheres with linearly anisotropic angular distribution of secondary neutrons, and we use an improved version of the P{sub N} method to obtain very accurate numerical results for them. The standard way of implementing the P{sub N} method for spheres has been shown to be numerically unstable by Aronson. Recently, in a joint work with C. E. Siewert of NCSU and J. R. Thomas, Jr. of VPI, we have developed a new version of the PN method for spheres that overcomes the difficulties pointed out by Aronson. In our implementation, we use singular-value decomposition (SVD) to solve, in a least squares sense, the system of equations obtained from the boundary condition for the coefficients of the approximation. In doing so, we can go to extremely high orders of the P{sub N} approximation (tens of thousands) without losing accuracy. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

Uniform Decay for Solutions of an Axially Moving Viscoelastic Beam

The paper deals with an axially moving viscoelastic structure modeled as an Euler–Bernoulli beam. The aim is to suppress the transversal displacement (transversal vibrations) that occur during the axial motion of the beam. It is assumed that the beam is moving with a constant axial speed and it is subject to a nonlinear force at the right boundary. We prove that when the axial speed of the beam is smaller than a critical value, the dissipation produced by the viscoelastic material is sufficient to suppress the transversal vibrations. It is shown that the rate of decay of the energy depends on the kernel which arise in the viscoelastic term. We consider a general kernel and notice that solutions cannot decay faster than the kernel.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1953

Method for Calculation of Laminar Heat Transfer in Air Flow Around Cylinders of Arbitrary Cross Section (including Large Temperature Differences and Transpiration Cooling)

The solution of heat-transfer problems has become vital for many aeronautical applications. The shapes of objects to be cooled can often be approximated by cylinders of various cross sections with flow normal to the axis as, for instance heat transfer on gas-turbine blades and on air foils heated for deicing purposes. A laminar region always exists near the stagnation point of such objects. A method previously presented by E. R. G. Eckert permits the calculation of local heat transfer around the periphery of cylinders of arbitrary cross section in the laminar region for flow of a fluid with constant property values with an accuracy sufficient for engineering purposes. The method is based on exact solutions of the boundary-layer equations for incompressible wedge-type flow and on the postulate that at any point on the cylinder the boundary-layer growth is the same as that on a wedge with comparable flow conditions. This method is extended herein to take into account the influence of large temperature differences between the cylinder wall and the flow as well as the influence of transpiration cooling when the same medium as the outside flow is used as coolant.

Eckert, E R↗

NASA NTRSDate not supplied

Evaluating Subgrid-Scale Properties of Low Clouds over the Beaufort Sea in Arctic System Reanalysis using ARISE Airborne In Situ Observations

Arctic low clouds influence the evolution of the Arctic system through their effects on radiative fluxes, boundary layer mixing, stability, turbulence, humidity, and precipitation. Unfortunately, atmospheric models and retrospective analysis (reanalysis) products struggle to accurately simulate the occurrence and properties of low clouds in the Arctic. To better understand the magnitude and possible causes of this problem, we analyze the disagreements between the Arctic System Reanalysis version 2 (ASR) and data taken during the September 2014 Arctic Radiation-IceBridge Sea and Ice Experiment (ARISE) airborne campaign conducted over the Beaufort Sea. Previous analysis of grid scale-averaged meteorological variables revealed a lack of cloud water produced in ASR, related to warm/dry biases, a failure to produce sufficient cloud water in conditions that are observed to be favorable for large cloud water values, and sampling issues with the ARISE flight paths. However, ARISE collected meteorological observations at one second intervals, offering spatial data resolution comparable to the subgrid scale of ASR. Subgrid scale properties and processes must be parameterized in conventional atmospheric models and reanalyses such as ASR, and often rely on assumed distributions of meteorological variables about the grid box mean state. To enable an examination of the realism of ASR’s subgrid scale properties, we use high resolution data from ARISE determine the “subgrid scale” distributions of cloud water and thermodynamic variables in the real atmosphere. Preliminary results show that the subgrid scale distribution of thermodynamic variables about the grid scale mean resembles the triangular distribution assumed in many bulk microphysical parameterization schemes, but also with tails resembling Gaussian distributions. The main goal is to understand the partitioning of the total available atmospheric water into water vapor and cloud water. Models typically covert some percentage of total water in a grid box to cloud water based on the estimated amount of supersaturation from the assumed vapor distribution. The ARISE measurements can reveal the relationship between subgrid scale supersaturation and the partition between cloud water and water vapor.

Jason Brant Dodson↗

DOE OSTI2017

Variational Integration for Ideal Magnetohydrodynamics and Formation of Current Singularities

Coronal heating has been a long-standing conundrum in solar physics. Parker's conjecture that spontaneous current singularities lead to nanoflares that heat the corona has been controversial. In ideal magnetohydrodynamics (MHD), can genuine current singularities emerge from a smooth 3D line-tied magnetic field? To numerically resolve this issue, the schemes employed must preserve magnetic topology exactly to avoid artificial reconnection in the presence of (nearly) singular current densities. Structure-preserving numerical methods are favorable for mitigating numerical dissipation, and variational integration is a powerful machinery for deriving them. However, successful applications of variational integration to ideal MHD have been scarce. In this thesis, we develop variational integrators for ideal MHD in Lagrangian labeling by discretizing Newcomb's Lagrangian on a moving mesh using discretized exterior calculus. With the built-in frozen-in equation, the schemes are free of artificial reconnection, hence optimal for studying current singularity formation. Using this method, we first study a fundamental prototype problem in 2D, the Hahm-Kulsrud-Taylor (HKT) problem. It considers the effect of boundary perturbations on a 2D plasma magnetized by a sheared field, and its linear solution is singular. We find that with increasing resolution, the nonlinear solution converges to one with a current singularity. The same signature of current singularity is also identified in other 2D cases with more complex magnetic topologies, such as the coalescence instability of magnetic islands. We then extend the HKT problem to 3D line-tied geometry, which models the solar corona by anchoring the field lines in the boundaries. The effect of such geometry is crucial in the controversy over Parker's conjecture. The linear solution, which is singular in 2D, is found to be smooth. However, with finite amplitude, it can become pathological above a critical system length. The nonlinear solution turns out smooth for short systems. Nonetheless, the scaling of peak current density vs. system length suggests that the nonlinear solution may become singular at a finite length. With the results in hand, we cannot confirm or rule out this possibility conclusively, since we cannot obtain solutions with system lengths near the extrapolated critical value.

solar corona↗

DOE OSTI2023

Hydro-Code Implementation and Testing of a Kinetic Phase Transition Framework

In this report we describe the Kinetic Phase Transition (KPT) framework that has been worked out over the last 10 years (from around 2014) and the implementation of it into three different codes, the one-dimensional hydro- LASLO and the three-dimensional magneto-hydro- ALEGRA, Sandia codes, via subroutines in the LAMBDA Equations of State and constitutive models package, and Flag, an arbitrary Lagrangian-Eulerian multiphysics code developed within the Lagrangian Applications project (LAP) at LANL. We discuss the introduction of phase mass (and/or volume) fractions that are needed in a code for it to be ‘phase aware’, that is, not only the thermodynamic state is known in each point but also the mixture of the materials’ phases in that point. Further we point to the need of a full Equations of State for each phase in a material to achieve phase awareness and we review the equilibrium phase model, where a phase mixture is at its lowest Gibbs free energy state, to make this point clear. Contrasting the kinetic phase transition to this equilibrium model seamlessly introduce us to the KPT framework that is subsequently thoroughly discussed. While the determination of the total state and the states and mass fractions of phases in each point is a problem that can borrow many of its numerical details from Eulerian codes and mixture of materials (not phases), the update of mass fractions with time in a KPT framework needs a new set of considerations. General for any update model is that we need to prevent mass fractions from becoming unphysical (negative or their sum to be larger than one). We have solved this problem by implementing a subdivision of the hydro time step that prevents the phase from being fully present to not present at all in one subdivided time step by limiting the size of the subdivided time step. This scheme also corrects numerical problems from abrupt changes in parameter values, the so called Gibbs phenomena, that gives rise to slushing between phases in the KPT framework. Interspersed throughout the report are discussions on different thermodynamics considerations. EOS validity windows, limitations on the EOS phase space, are needed for the KPT framework and are discussed separately and exemplified. The KPT framework described in this report has been verified by code comparison, but validation is still an active area of research. There is room for improvement in the update model, both in the model for determination of rates and in how to prevent the mass fractions from becoming unphysical. In addition, the parameters in the KPT update model and the placement of the phase boundary in the EOS phase space, and interactions with other constitutive models, are closely related and interfering with each other. One possible way forward is to simultaneously develop KPT parameters, EOS, and constitutive models for each material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2014

Mean Flow Augmented Acoustics in Rocket Systems

Oscillatory motion in solid rocket motors and liquid engines has long been a subject of concern. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process and gas dynamics. The customary approach to modeling acoustic waves inside a rocket chamber is to apply the classical inhomogeneous wave equation to the combustion gas. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while the acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The converging section of a rocket nozzle, where gradients in pressure, density, and velocity become large, is a notable region where this approach is not applicable. The expulsion of unsteady energy through the nozzle of a rocket is identified as the predominate source of acoustic damping for most rocket systems. An accurate model of the acoustic behavior within this region where acoustic modes are influenced by the presence of a steady mean flow is required for reliable stability predictions. Recently, an approach to address nozzle damping with mean flow effects was implemented by French [1]. This new approach extends the work originated by Sigman and Zinn [2] by solving the acoustic velocity potential equation (AVPE) formulated by perturbing the Euler equations [3]. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, (del squared)(psi) − (lambda/c)(exp 2)(psi) − M(dot)[M(dot)(del)(del(psi))] − 2(lambda(M/c) + (M(dot)del(M))(dot)del(psi)−2(lambda)(psi)[M(dot)del(1/c)]=0 (1) with M as the Mach vector, c as the speed of sound, and lambda as the complex eigenvalue. French apply the finite volume method to solve the steady flow field within the combustion chamber and nozzle with inviscid walls. The complex eigenvalues and eigenvector are determined with the use of the ARPACK eigensolver. The present study employs the COMSOL Multphysics framework to solve the coupled eigenvalue problem using the finite element approach. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluated the gas flow inside of a solid rocket motor using St. Robert's law modeling solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are used by the coefficient form of the mathematics module to determine the eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. Pertinent results from these analyses are the complex valued eigenvalue and eigenvectors. Comparisons are made to the French results to evaluate the modeling approach. A comparison of the French results with that of the present analysis is displayed in figures 1 and 2, respectively. The graphic shows the first tangential eigenvector's real (a) and imaginary (b) values.

Fischbach, Sean R.↗