Search NASA⌕ Search

SEARCH · Search NASA

Results for “EQUATIONS - ENERGY”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Radiation boundary condition and anisotropy correction for finite difference solutions of the Helmholtz equation

In this paper finite-difference solutions of the Helmholtz equation in an open domain are considered. By using a second-order central difference scheme and the Bayliss-Turkel radiation boundary condition, reasonably accurate solutions can be obtained when the number of grid points per acoustic wavelength used is large. However, when a smaller number of grid points per wavelength is used excessive reflections occur which tend to overwhelm the computed solutions. Excessive reflections are due to the incompability between the governing finite difference equation and the Bayliss-Turkel radiation boundary condition. The Bayliss-Turkel radiation boundary condition was developed from the asymptotic solution of the partial differential equation. To obtain compatibility, the radiation boundary condition should be constructed from the asymptotic solution of the finite difference equation instead. Examples are provided using the improved radiation boundary condition based on the asymptotic solution of the governing finite difference equation. The computed results are free of reflections even when only five grid points per wavelength are used. The improved radiation boundary condition has also been tested for problems with complex acoustic sources and sources embedded in a uniform mean flow. The present method of developing a radiation boundary condition is also applicable to higher order finite difference schemes. In all these cases no reflected waves could be detected. The use of finite difference approximation inevita bly introduces anisotropy into the governing field equation. The effect of anisotropy is to distort the directional distribution of the amplitude and phase of the computed solution. It can be quite large when the number of grid points per wavelength used in the computation is small. A way to correct this effect is proposed. The correction factor developed from the asymptotic solutions is source independent and, hence, can be determined once and for all. The effectiveness of the correction factor in providing improvements to the computed solution is demonstrated in this paper.

Tam, Christopher K. W.↗
NASA NTRS1995

The limitation and applicability of Musher-Sturman equation to two dimensional lower hybrid wave collapse

The existence of localized regions of intense lower hybrid waves in the auroral ionosphere recently observed by rocket and satellite experiments can be understood by the study of a non-linear two-timescale coupling process. In this Letter, we demonstrate that the leading non-linear term in the standard Musher-Sturman equation vanishes identically in strict two-dimensions (normal to the magnetic field). Instead, the new two-dimensional equation is characterized by a much weaker non-linear term which arises from the ponderomotive force perpendicular to the magnetic field, particularly that due to the ions. The old and new equations are compared by means of time-evolution calculations of wave fields. The results exhibit a remarkable difference in the evolution of the waves as governed by the two equations. Such dissimilar outcomes motivate our investigation of the limitation of Musher-Sturman equation in quasi-two-dimensions. Only within all these limits can Musher-Sturman equation adequately describe the collapse of lower hybrid waves.

Tam, Sunny W. Y.↗
NASA NTRS1995

On absorbing boundary conditions for linearized Euler equations by a perfectly matched layer

Recently, Berenger introduced a Perfectly Matched Layer (PML) technique for absorbing electromagnetic waves. In the present paper, a perfectly matched layer is proposed for absorbing out-going two-dimensional waves in a uniform mean flow, generated by linearized Euler equations. It is well known that the linearized Euler equations support acoustic waves, which travel with the speed of sound relative to the mean flow, and vorticity and entropy waves, which travel with the mean flow. The PML equations to be used at a region adjacent to the artificial boundary for absorbing these linear waves are defined. Plane waves solutions to the PML equations are developed and wave propagation and absorption properties are given. It is shown that the theoretical reflection coefficients at an interface between the Euler and PML domains are zero, independent of the angle of incidence and frequency of the waves. As such, the present study points out a possible alternative approach for absorbing out-going waves of the Euler equations with little or no reflection in computation. Numerical examples that demonstrate the validity of the proposed PML equations are also presented.

Hu, Fang Q.↗
NASA NTRS1994

Computation of Separated and Unsteady Flows with One- and Two-Equation Turbulence Models

The ability of one- and two-equation turbulence models to predict unsteady separated flows over airfoils is evaluated. An implicit, factorized, upwind-biased numerical scheme is used for the integration of the compressible, Reynolds averaged Navier-Stokes equations. The turbulent eddy viscosity is obtained from the computed mean flowfield by integration of the turbulent field equations. The two-equation turbulence models are discretized in space with an upwind-biased, second order accurate total variation diminishing scheme. One and two-equation turbulence models are first tested for a separated airfoil flow at fixed angle of incidence. The same models are then applied to compute the unsteady flowfields about airfoils undergoing oscillatory motion at low subsonic Mach numbers. Experimental cases where the flow has been tripped at the leading edge and where natural transition was allowed to occur naturally are considered. The more recently developed field-equation turbulence models capture the physics of unsteady separated flow significantly better than the standard kappa-epsilon and kappa-omega models. However, certain differences in the hysteresis effects are obtained. For an untripped high-Reynolds-number flow, it was found necessary to take into account the leading edge transitional flow region in order to capture the correct physical mechanism that leads to dynamic stall.

Ekaterinaris, John A.↗
NASA NTRS1994

Progress in Computational Aeroelasticity Using High Fidelity Flow and Structural Equations on Parallel Computers

Aeroelasticity which involves strong coupling of fluids, structures and controls is an important element in designing an aircraft. Computational aeroelasticity using low fidelity methods such as the linear aerodynamic flow equations coupled with the modal structural equations are well advanced. Though these low fidelity approaches are computationally less intensive, they are not adequate for the analysis of modern aircraft such as High Speed Civil Transport (HSCT) and Advanced Subsonic Transport (AST) which can experience complex flow/structure interactions. HSCT can experience vortex induced aeroelastic oscillations whereas AST can experience transonic buffet associated structural oscillations. Both aircraft may experience a dip in the flutter speed at the transonic regime. For accurate aeroelastic computations at these complex fluid/structure interaction situations, high fidelity equations such as the Navier-Stokes for fluids and the finite-elements for structures are needed. Computations using these high fidelity equations require large computational resources both in memory and speed. Current conventional supercomputers have reached their limitations both in memory and speed. As a result, parallel computers have evolved to overcome the limitations of conventional computers. This paper will address the transition that is taking place in computational aeroelasticity from conventional computers to parallel computers. The paper will address special techniques needed to take advantage of the architecture of new parallel computers. Results will be illustrated from computations made on iPSC/860 and IBM SP2 computer by using ENASERO code that directly couples the Euler/Navier-Stokes flow equations with high resolution finite-element structural equations.

Guruswamy, Guru P.↗
NASA NTRS1994

On the Connection Between One-and Two-Equation Models of Turbulence

A formalism will be presented that allows the transformation of two-equation eddy viscosity turbulence models into one-equation models. The transformation is based on an assumption that is widely accepted over a large range of boundary layer flows and that has been shown to actually improve predictions when incorporated into two-equation models of turbulence. Based on that assumption, a new one-equation turbulence model will be derived. The new model will be tested in great detail against a previously introduced one-equation model and against its parent two-equation model.

Menter, F. R.↗
NASA NTRS1994

Fluid/Structure Interaction Studies of Aircraft Using High Fidelity Equations on Parallel Computers

Abstract Aeroelasticity which involves strong coupling of fluids, structures and controls is an important element in designing an aircraft. Computational aeroelasticity using low fidelity methods such as the linear aerodynamic flow equations coupled with the modal structural equations are well advanced. Though these low fidelity approaches are computationally less intensive, they are not adequate for the analysis of modern aircraft such as High Speed Civil Transport (HSCT) and Advanced Subsonic Transport (AST) which can experience complex flow/structure interactions. HSCT can experience vortex induced aeroelastic oscillations whereas AST can experience transonic buffet associated structural oscillations. Both aircraft may experience a dip in the flutter speed at the transonic regime. For accurate aeroelastic computations at these complex fluid/structure interaction situations, high fidelity equations such as the Navier-Stokes for fluids and the finite-elements for structures are needed. Computations using these high fidelity equations require large computational resources both in memory and speed. Current conventional super computers have reached their limitations both in memory and speed. As a result, parallel computers have evolved to overcome the limitations of conventional computers. This paper will address the transition that is taking place in computational aeroelasticity from conventional computers to parallel computers. The paper will address special techniques needed to take advantage of the architecture of new parallel computers. Results will be illustrated from computations made on iPSC/860 and IBM SP2 computer by using ENSAERO code that directly couples the Euler/Navier-Stokes flow equations with high resolution finite-element structural equations.

Guruswamy, Guru↗
NASA NTRS1995

On the Dynamics of Some Discretizations of Convection-Diffusion Equations

Numerical discretizations of differential equations which model physical processes can possess dynamics quite different from that of the equations themselves. Recently the emphasis has been on the the dynamics of numerical discretizations for Ordinary Differential Equations (ODEs). For Partial Differential Equations (PDEs) using a method of lines approach the situation is more complex. First, the spatial discretisation may introduce dynamics not present in the original equations; second, the solution of the resulting system of ODEs is open to the modified dynamics of the ODE solver used. These two effects may interact in a complex manner. In this talk we present some results of our recent work on the dynamics of discretizations of convection-diffusion equations, including those produced using Total Variation Diminishing (TVD) schemes and adaptive grid techniques. A more general overview of the area may be found on our accompanying poster presentation.

Sweby, Peter K.↗
NASA NTRS2001

A Highly Accurate Technique for the Treatment of Flow Equations at the Polar Axis in Cylindrical Coordinates using Series Expansions

Numerical methods for solving the flow equations in cylindrical or spherical coordinates should be able to capture the behavior of the exact solution near the regions where the particular form of the governing equations is singular. In this work we focus on the treatment of these numerical singularities for finite-differences methods by reinterpreting the regularity conditions developed in the context of pseudo-spectral methods. A generally applicable numerical method for treating the singularities present at the polar axis, when nonaxisymmetric flows are solved in cylindrical, coordinates using highly accurate finite differences schemes (e.g., Pade schemes) on non-staggered grids, is presented. Governing equations for the flow at the polar axis are derived using series expansions near r=0. The only information needed to calculate the coefficients in these equations are the values of the flow variables and their radial derivatives at the previous iteration (or time) level. These derivatives, which are multi-valued at the polar axis, are calculated without dropping the accuracy of the numerical method using a mapping of the flow domain from (0,R)*(0,2pi) to (-R,R)*(0,pi), where R is the radius of the computational domain. This allows the radial derivatives to be evaluated using high-order differencing schemes (e.g., compact schemes) at points located on the polar axis. The proposed technique is illustrated by results from simulations of laminar-forced jets and turbulent compressible jets using large eddy simulation (LES) methods. In term of the general robustness of the numerical method and smoothness of the solution close to the polar axis, the present results compare very favorably to similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r=0, following the method proposed recently by Mohseni and Colonius (1). Extension of the method described here for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical or spherical coordinates with finite-differences schemes of various level of accuracy is immediate.

Constantinescu, George S.↗
NASA NTRS2002

Local Discontinuous Galerkin Methods for Partial Differential Equations with Higher Order Derivatives

In this paper we review the existing and develop new continuous Galerkin methods for solving time dependent partial differential equations with higher order derivatives in one and multiple space dimensions. We review local discontinuous Galerkin methods for convection diffusion equations involving second derivatives and for KdV type equations involving third derivatives. We then develop new local discontinuous Galerkin methods for the time dependent bi-harmonic type equations involving fourth derivatives, and partial differential equations involving fifth derivatives. For these new methods we present correct interface numerical fluxes and prove L(exp 2) stability for general nonlinear problems. Preliminary numerical examples are shown to illustrate these methods. Finally, we present new results on a post-processing technique, originally designed for methods with good negative-order error estimates, on the local discontinuous Galerkin methods applied to equations with higher derivatives. Numerical experiments show that this technique works as well for the new higher derivative cases, in effectively doubling the rate of convergence with negligible additional computational cost, for linear as well as some nonlinear problems, with a local uniform mesh.

Yan, Jue↗
NASA NTRS2010

Forces Associated with Nonlinear Nonholonomic Constraint Equations

A concise method has been formulated for identifying a set of forces needed to constrain the behavior of a mechanical system, modeled as a set of particles and rigid bodies, when it is subject to motion constraints described by nonholonomic equations that are inherently nonlinear in velocity. An expression in vector form is obtained for each force; a direction is determined, together with the point of application. This result is a consequence of expressing constraint equations in terms of dot products of vectors rather than in the usual way, which is entirely in terms of scalars and matrices. The constraint forces in vector form are used together with two new analytical approaches for deriving equations governing motion of a system subject to such constraints. If constraint forces are of interest they can be brought into evidence in explicit dynamical equations by employing the well-known nonholonomic partial velocities associated with Kane's method; if they are not of interest, equations can be formed instead with the aid of vectors introduced here as nonholonomic partial accelerations. When the analyst requires only the latter, smaller set of equations, they can be formed directly; it is not necessary to expend the labor to form the former, larger set first and subsequently perform matrix multiplications.

Roithmayr, Carlos M.↗
NASA NTRS2018

Output Measurement Equations for Flexible Aircraft Flight Dynamics

A summary of output measurement equations for onboard sensors used in flight testing of flexible aircraft is presented. These equations include the effects of structural flexibility and are considerably more complex than the standard equations for rigid-body aircraft. The output equations discussed include accelerations from linear accelerometers, strains, angular rates, angular accelerations, Euler angles, true airspeed, and air flow angles. The output equations are derived in full form and then simplified in some cases. Linearized output equations, suitable for state-space or transfer function models, are also developed. Example flight test data from the X-56A subscale aeroelastic demonstrator is discussed, for reference.

Grauer, Jared A.↗
NASA NTRS2022

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗
NASA NTRSDate not supplied

Evaluation of the Accuracy of the Load Prediction Equations of Low-Load Balance Calibration Data

Accuracy and reliability of the load prediction equations of a low-load calibration data set of a force balance were investigated. First, independent load prediction equations were generated from the data of a full-load and low-load machine calibration. Then, the low- load data set was processed as a set of precision check loads for the load prediction equations that were obtained from the full-load and low-load calibrations. Finally, the load prediction equations of the two calibrations were applied to four manual check load data sets that were recorded between 2009 and 2022. The residuals of the predicted calibration and check loads were compared in all cases. No systematic improvement of the load prediction accuracy was observed when the load prediction equations of the low-load calibration were applied to the check load data sets. Therefore, it is recommended to apply a full-load calibration to a six-component balance even if the balance is not used across its entire design envelope. This approach has the advantage that the resulting load prediction equations are less likely to be applied outside of the calibration load ranges. In addition, it is expected that the numerical estimates of the gage sensitivities are more reliable assuming that the uncertainties of the applied loads and measured outputs of the calibration data are more or less constant across the load range of the balance.

strain-gage balance↗
NASA NTRSDate not supplied

Evaluation of the Accuracy of the Load Prediction Equations of Low-Load Balance Calibration Data

Accuracy and reliability of the load prediction equations of a low-load calibration data set of a force balance were investigated. First, independent load prediction equations were generated from the data of a full-load and low-load machine calibration. Then, the low-load data set was processed as a set of precision check loads for the load prediction equations that were obtained from the full-load and low-load calibrations. Finally, the load prediction equations of the two calibrations were applied to four manual check load data sets that were recorded between 2009 and 2022. The residuals of the predicted calibration and check loads were compared in all cases. No systematic improvement of the load prediction accuracy was observed when the load prediction equations of the low-load calibration were applied to the check load data sets. Therefore, it is recommended to apply a full-load calibration to a six-component balance even if the balance is not used across its entire design envelope. This approach has the advantage that the resulting load prediction equations are less likely to be applied outside of the calibration load ranges. In addition, it is expected that the numerical estimates of the bridge sensitivities are more reliable assuming that the uncertainties of the applied loads and measured outputs of the calibration data are more or less constant across the load range of the balance.

strain-gage balance↗
NASA NTRSDate not supplied

Construction of an Exact Pressure-Equilibrium Scheme for the Five-Equation Two-Phase Flow Model With Thermal Relaxation

Numerical simulation of compressible multiphase flows based on the four-equation (homogeneous relaxation) model is known to suffer from two fundamental difficulties with (a) wave propagation and (b) pressure equilibrium preservation. First, the mixture sound speed exhibits non-monotonic dependency with respect to the volume fraction, which leads to robustness issues in the resolution of shocks and acoustic wave propagation across two-phase regions. This difficulty can be mitigated by solving Allaire’s five-equation model augmented with infinitely fast phasic temperature equilibrium, from which solutions of the four-equation model can be recovered. However, when temperature is non-uniform, this augmented five-equation formulation still fails to preserve pressure equilibrium across material interfaces. In this work, we propose a fully conservative numerical scheme that exactly preserves pressure equilibrium at the discrete level for the augmented five-equation model, for arbitrary initial distributions of temperature and volume fraction. Combined with the monotonic sound speed property of the five-equation formulation, the proposed pressure-equilibrium preserving scheme significantly improves robustness in the presence of strong multiphase interactions, including shock–interface interactions and advection of material interfaces.

ESG↗
DOE OSTI2025

Constrained or unconstrained? Neural-network-based equation discovery from data

Throughout many fields, practitioners often rely on differential equations to model systems. Yet, for many applications, the theoretical derivation of such equations and/or the accurate resolution of their solutions may be intractable. Instead, recently developed methods, including those based on parameter estimation, operator subset selection, and neural networks, allow for the data-driven discovery of both ordinary and partial differential equations (PDEs), on a spectrum of interpretability. The success of these strategies is often contingent upon the correct identification of representative equations from noisy observations of state variables and, as importantly and intertwined with that, the mathematical strategies utilized to enforce those equations. Specifically, the latter has been commonly addressed via unconstrained optimization strategies. Representing the PDE as a neural network, we propose to discover the PDE (or the associated operator) by solving a constrained optimization problem and using an intermediate state representation similar to a physics-informed neural network (PINN). The objective function of this constrained optimization problem promotes matching the data, while the constraints require that the discovered PDE is satisfied at a number of spatial collocation points. We present a penalty method and a widely used trust-region barrier method to solve this constrained optimization problem, and we compare these methods on numerical examples. Our results on several example problems demonstrate that the latter constrained method outperforms the penalty method, particularly for higher noise levels or fewer collocation points. This work motivates further exploration into using sophisticated constrained optimization methods in scientific machine learning, as opposed to their commonly used, penalty-method or unconstrained counterparts. For both of these methods, we solve these discovered neural network PDEs with classical methods, such as finite difference methods, as opposed to PINNs-type methods relying on automatic differentiation. Here, we briefly highlight how simultaneously fitting the data while discovering the PDE improves the robustness to noise and other small, yet crucial, implementation details.

Data-driven discovery↗
DOE OSTI2021

Efficient quantum algorithm for dissipative nonlinear differential equations

Significance Nonlinear differential equations appear in many domains and are notoriously difficult to solve. Whereas previous quantum algorithms for general nonlinear differential equations have complexity exponential in the evolution time, we give the first quantum algorithm for dissipative nonlinear differential equations that is efficient provided the dissipation is sufficiently strong relative to nonlinear and forcing terms and the solution does not decay too rapidly. We also establish a lower bound showing that differential equations with sufficiently weak dissipation have worst-case complexity exponential in time, giving an almost tight classification of the quantum complexity of simulating nonlinear dynamics. Furthermore, numerical results for the Burgers equation suggest that our algorithm may potentially address complex nonlinear phenomena even in regimes with weaker dissipation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗