SEARCH · Engineering Papers

Results for “singular equations”

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 433 records · Page 24

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Scaling of Shear Rheology of Concentrated Charged Colloidal Suspensions across Glass Transition

Electrostatic interparticle interactions are a key component in controlling and designing rheological characteristics of concentrated charged colloidal suspensions. Herein, we investigate electroviscous effects on shear rheology using highly charged silica particles. By fixing the volume fraction but varying the salinity, the system undergoes a glass transition as evidenced by the evolution of the yield stress and zero-shear viscosity. We show that the steady shear viscosities obey a critical scaling relation that scales the flow curves into a super- and a sub-critical branch with glass transition salinity serving as the bifurcation point; we also demonstrate an isoviscosity scaling that collapses all isoviscosity lines into a single master curve that exhibits no singularity. Based on each scaling relation, in conjunction with common modeling equations, the quantitative relationships between the shear viscosity, stress, and salinity are established. This study demonstrates a new framework to model the steady shear rheology of concentrated charged colloids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Two-way coupled particle-turbulence interaction: Effect of numerics and resolution on fluid and particle statistics

Euler-Lagrange point-particle simulation has emerged as a premier methodology for studying dispersed particle-laden flows. This method's popularity stems from its ability to resolve fine-scale fluid structures while also tracking individual particles at reduced cost using an appropriate particle acceleration model. However, the point-particle model has known convergence issues in that refinement of the fluid grid can lead to changes in the predicted statistics. The reasons for non-convergence are twofold: the point-particle two-way coupling force in the Navier-Stokes equations requires a numerical regularization and, without careful implementation, yields a singular force on the fluid with grid refinement. The second factor that yields grid-dependent statistics is that the point-particle force model in general depends on the undisturbed fluid velocity. When the undisturbed fluid velocity is not robustly modeled in a grid-insensitive way, the calculated force for both particles and fluid will be grid-dependent, contaminating their respective statistics. While the first issue regarding regularizing the point-particle source term has received attention in the literature, the consequences of robustly modeling the undisturbed velocity in the context of grid refinement of turbulence has received little attention. Here, we consider decaying homogeneous isotropic turbulence laden with particles at different Stokes numbers. For a given Stokes number, we systematically refine the grid and demonstrate that explicitly modeling the undisturbed fluid velocity yields relative grid insensitivity for the energy of the particle and fluid phases, as well as acceleration of the particles. We also demonstrate that an appropriately defined dissipation rate is also grid-insensitive when an undisturbed fluid velocity correction is used. In contrast, when the undisturbed fluid velocity is modeled using the conventional approach of interpolating the local fluid velocity to the particle location, we show this procedure yields divergent statistics with grid refinement. In particular, we show that higher-order interpolation of the fluid velocity in two-way coupled problems is worse than lower-order interpolation, in the absence of a correction procedure to estimate the undisturbed fluid velocity. We also examine velocity derivative statistics of the fluid phase and demonstrate that these statistics are not in general convergent even when the undisturbed fluid velocity is explicitly modeled. Collectively, the observations in this paper are used to present a philosophy on the types of questions which are answerable with point-particle methods.

42 ENGINEERING↗

NASA NTRS1974

An extension of the Laplace transform to Schwartz distributions

A characterization of the Laplace transform is developed which extends the transform to the Schwartz distributions. The class of distributions includes the impulse functions and other singular functions which occur as solutions to ordinary and partial differential equations. The standard theorems on analyticity, uniqueness, and invertibility of the transform are proved by using the characterization as the definition of the Laplace transform. The definition uses sequences of linear transformations on the space of distributions which extends the Laplace transform to another class of generalized functions, the Mikusinski operators. It is shown that the sequential definition of the transform is equivalent to Schwartz' extension of the ordinary Laplace transform to distributions but, in contrast to Schwartz' definition, does not use the distributional Fourier transform. Several theorems concerning the particular linear transformations used to define the Laplace transforms are proved. All the results proved in one dimension are extended to the n-dimensional case, but proofs are presented only for those situations that require methods different from their one-dimensional analogs.

Price, D. R.↗

NASA NTRS1977

Interacting turbulent boundary layer over a wavy wall

This paper is concerned with the two-dimensional supersonic flow of a thick turbulent boundary layer over a train of relatively small wave-like protuberances. The flow conditions and the geometry are such that there exists a strong interaction between the viscous and inviscid flow. The problem cannot be solved without inclusion of interaction effects due to the occurrence of the separation singularity in classical boundary layer methods. Here the interacting boundary layer equations are solved numerically using a time-like relaxation method with turbulence effects represented by the inclusion of the eddy viscosity model of Cebeci and Smith. Results are presented for flow over a train of up to six waves for Mach numbers of 2.5 and 3.5, Reynolds numbers of 10,000,000/m and 32,000,000/m, and wall temperature ratios of 0.4 and 0.8. Limited comparisons with independent experimental and analytical results are also given.

Polak, A.↗

NASA NTRS1982

Higher-Order Panel Method for Aerodynamic Flow Analysis

PANAIR uses high-order panel method to predict inviscid subsonic or supersonic flows about arbitrary configuration. Panel method solves linear partial differential equation numerically by approximating configuration surface with panels on which unknown "singularity strengths" are defined. PANAIR includes advanced software technology as well as advanced aerodynamic technology.

Erickson, L.↗

NASA NTRS1982

The motion of a drop in a rigid surface

The assumptions of a region of slip near a moving contact line (to remove the force singularity) and a constant contact angle are used to obtain the equation for the shape of a thin drop of liquid resting on a rigid plane. Three asymptotic expansions are matched together to obtain an expression for the rate at which the drop spreads. Some cases of sliding motion are also examined. Although the technique is presented here for thin drops only, it can also be applied to drops of arbitrary size.

Hocking, L. M.↗

NASA NTRS1984

Effect of boundary layers on solid walls in three-dimensional subsonic wind tunnels

A solution for the tunnel wall boundary layer effects for three-dimensional subsonic tunnels is presented. The model potentials are represented with simple singularities placed on the centerline of the tunnel and Laplace's equation in cylindrical coordinates is solved for either the conventional homogeneous slotted-wall boundary condition, the solid-wall viscous boundary condition, or a combination of them. The most pronounced wall boundary layer effect is on solid blockage for completely closed wind tunnels. Boundary layers on the wall reduce the blockage from the solid-wall, no-boundary-layer case in a manner similar to opening slots in a solid wall. Additionally, for solid-wall tunnel configurations, the streamline curvature interference factor is reduced by a significant amount, whereas the lift interference factor at the model station does not depend on the boundary layer parameter. For combination wall configurations, the slot effect of the horizontal walls dominates the viscous effect of the solid sidewalls.

Adcock, J. B.↗

NASA NTRS1988

Deconvolution results for wide field-of-view radiometer measurements of reflected solar radiation

The measurement of reflected solar radiation by spacecraft wide-field-of-view radiometers is considered analytically, and the method proposed by Smith (1987) for the solution of the discretized governing equations is described and demonstrated. The theoretical basis of the singular-value decomposition method (partitioning the albedo field into observable and unobservable components) is outlined; its application to Nimbus 7 ERB data is explained in detail; and numerical results are presented in graphs and maps and briefly characterized.

Smith, G. Louis↗

NASA NTRS2021

Thermal Inspection of a Composite Fuselage Section Using theMethod of Proper Generalized Decomposition

Proper Generalized Decomposition (PGD) is a reduced order modeling technique for the simulation of physical systems whose governing equations depend on boundary conditions, initial conditions, material properties, and geometric parameters. It uses separated representations of system covariates combined with an iterative approximation method known as successive enrichment in order to compute an accurate parameter-dependent approximation to the full governing equations. PGD can also be used as an alternative to the Singular Value Decomposition (SVD) of a matrix and therefore as an alternative to PCA thermography. In this paper PGD was used to analyze data derived from the inspection of a composite fuselage forward section using flash thermography, and the results were compared against the standard PCA approach.

Nondestructive Evaluation↗

NASA NTRS1995

Fully-Explicit and Self-Consistent Algebraic Reynolds Stress Models

A fully-explicit, self-consistent algebraic expression for the Reynolds stress, which is the exact solution to the Reynolds stress transport equation in the 'weak equilibrium' limit for two-dimensional mean flows for all linear and some quasi-linear pressure-strain models, is derived. Current explicit algebraic Reynolds stress models derived by employing the 'weak equilibrium' assumption treat the production-to-dissipation (P/epsilon) ratio implicitly, resulting in an effective viscosity that can be singular away from the equilibrium limit. In the present paper, the set of simultaneous algebraic Reynolds stress equations are solved in the full non-linear form and the eddy viscosity is found to be non-singular. Preliminary tests indicate that the model performs adequately, even for three dimensional mean flow cases. Due to the explicit and non-singular nature of the effective viscosity, this model should mitigate many of the difficulties encountered in computing complex turbulent flows with the algebraic Reynolds stress models.

Girimaji, Sharath S.↗

NASA NTRS1977

Mathematical and computational studies of equilibrium capillary free surfaces

The results of several independent studies are presented. The general question is considered of whether a wetting liquid always rises higher in a small capillary tube than in a larger one, when both are dipped vertically into an infinite reservoir. An analytical investigation is initiated to determine the qualitative behavior of the family of solutions of the equilibrium capillary free-surface equation that correspond to rotationally symmetric pendent liquid drops and the relationship of these solutions to the singular solution, which corresponds to an infinite spike of liquid extending downward to infinity. The block successive overrelaxation-Newton method and the generalized conjugate gradient method are investigated for solving the capillary equation on a uniform square mesh in a square domain, including the case for which the solution is unbounded at the corners. Capillary surfaces are calculated on the ellipse, on a circle with reentrant notches, and on other irregularly shaped domains using JASON, a general purpose program for solving nonlinear elliptic equations on a nonuniform quadrilaterial mesh. Analytical estimates for the nonexistence of solutions of the equilibrium capillary free-surface equation on the ellipse in zero gravity are evaluated.

Albright, N.↗

NASA NTRS1976

Sound propagation in choked ducts

The linearized equations describing the propagation of sound in variable area ducts containing flow are shown to be singular when the duct mean flow is sonic. The singularity is removed when previously ignored nonlinear terms are retained. The results of a numerical study, for the case of plane waves propagating in a one-dimensional converging-diverging duct, show that the sound field is adequately described by the linearized equations only when the axial mean flow Mach number at the duct throat M sub th 0.6. For M sub th 0.6, the numerical results showed that acoustic energy flux was not conserved. An attempt was made to extend the study to include the nonlinear behavior of the sound field. Meaningful results were not obtained due, primarily, to numerical difficulties.

Hersh, A. S.↗

NASA NTRS1985

Application of singular value decomposition to structural dynamics systems with constraints

Singular value decomposition is used to construct a coordinate transformation for a linear dynamic system subject to linear, homogeneous constraint equations. The method is compared with two commonly used methods, namely classical Gaussian elimination and Walton-Steeves approach. Although the classical method requires fewer numerical operations, the singular value decomposition method is more accurate and convenient in eliminating the dependent coordinates. Numerical examples are presented to demonstrate the application of the method.

Juang, J.-N.↗

NASA NTRS1993

Energy transfer and constrained simulations in isotropic turbulence

The defining characteristic of turbulent flows is their ability to dissipate energy, even in the limit of zero viscosity. The Euler equations, if constrained in such a way that the velocity derivatives remain bounded, conserve energy. But when they arise as the limit of the Navier-Stokes (NS) equations, when the Reynolds number goes to infinity, there is persuasive empirical evidence that the gradients become singular as just the right function of Re for the dissipation to remain non-zero and to approach a well defined limit. It is generally believed that this limiting value of the dissipation is a property of the Euler equations themselves, independent of the particular dissipative mechanism involved, and that it can be normalized with the large scale properties of the turbulent flow (e.g. the kinetic energy per unit volume u'(exp 2)/2, and the integral scale L) without reference to the Reynolds number or to other dissipative quantities. This is usually taken to imply that the low wave number end of the energy spectrum, far from the dissipative range, is also independent of the particular mechanism chosen to dispose of the energy transfer. In the following sections, we present some numerical experiments on the effect of substituting different dissipation models into the truncated Euler equations. We will see that the effect is mainly felt in the 'near dissipation' range of the energy spectrum, but that this range can be quite wide in some cases, contaminating a substantial range of wave numbers. In the process, we will develop a 'practical' approximation to the subgrid energy transfer in isotropic turbulence, and we will gain insight into the structure of the nonlinear interactions among turbulent scales of comparable size, and into the nature of energy backscatter. Some considerations on future research directions are offered at the end.

Jimenez, Javier↗

DOE OSTI2024

Inflation alternative via the gravitational field of a singularity

Abstract We explore the scenario that the observable Universe emerged from the vicinity of a negative mass ring singularity, and all content of the Universe travels at the same group velocity close to the speed of light on a geodesic trajectory along the axis of rotation of the singularity. In appropriate coordinate parametrization and evaluated on the trajectory, we find that the metric tensor in the vicinity of the trajectory exhibits a conformal scale factor a ( η ) with contraction and subsequent expansion properties that solve the horizon problem. We then introduce a static flow of gravitating radiation along the trajectory (perturbatively with respect to the mass scale of the singularity) to model a homogeneous radiation dominated Universe. Solving the Einstein field equations with a physically motivated ansatz of metric perturbation then reveals that the effective conformal scale factor indeed grows asymptotically with the same power law as expected in a conventional radiation dominated Universe.

Astronomy & Astrophysics↗

NASA NTRS1978

Multi-Level Adaptive Techniques (MLAT) for singular-perturbation problems

The multilevel (multigrid) adaptive technique, a general strategy of solving continuous problems by cycling between coarser and finer levels of discretization is described. It provides very fast general solvers, together with adaptive, nearly optimal discretization schemes. In the process, boundary layers are automatically either resolved or skipped, depending on a control function which expresses the computational goal. The global error decreases exponentially as a function of the overall computational work, in a uniform rate independent of the magnitude of the singular-perturbation terms. The key is high-order uniformly stable difference equations, and uniformly smoothing relaxation schemes.

Brandt, A.↗

NASA NTRS1983

Compact finite difference schemes for the Euler and Navier-Stokes equations

Implicit compact finite difference schemes for the Euler equations are described which furnish equivalent treatment of the conservation and nonconservation forms; a simple modification yields an entropy-producing scheme. An extension of the scheme also treats the compressible Navier-Stokes equations; when the viscosity and heat conduction coefficients are negligible only the boundary data appropriate to the Euler equation influence the solution to any significant extent, a result consistent with singular perturbation theory.

Rose, M. E.↗

NASA NTRS1983

Separation of time scales in aircraft trajectory optimization

Two methods for analyzing the time-scale properties of aircraft trajectory optimization problems are presented. Time-scale properties must be identified before solutions can be obtained by using singular perturbation methods. Both methods only require a knowledge of the state equations, the aircraft characteristics, and the bounds on the state and control variables. Although these methods give only rough estimates of time-scale separation, they do not require that an 'exact' optimal trajectory be known, as do the more rigorous methods, and they are an improvement on the ad hoc methods currently in use. The two methods are applied to an example problem for a high performance aircraft.

Ardema, M. D.↗