Search NASA⌕ Search

SEARCH · Search NASA

Results for “MATHEMATICAL SOLUTIONS”

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 343 records · Page 19

The B-dot Earth Average Magnetic Field

The average Earth's magnetic field is solved with complex mathematical models based on mean square integral. Depending on the selection of the Earth magnetic model, the average Earth's magnetic field can have different solutions. This paper presents a simple technique that takes advantage of the damping effects of the b-dot controller and is not dependent of the Earth magnetic model; but it is dependent on the magnetic torquers of the satellite which is not taken into consideration in the known mathematical models. Also the solution of this new technique can be implemented so easily that the flight software can be updated during flight, and the control system can have current gains for the magnetic torquers. Finally, this technique is verified and validated using flight data from a satellite that it has been in orbit for three years.

Capo-Lugo, Pedro A.↗

Superposed hyperbolic kink and pulse solutions of coupled $φ$ 4 , NLS and mKdV equations

Here, in this paper, we obtain novel solutions of a coupled $φ$ 4 , a coupled nonlinear Schrödinger equation and a coupled modified Korteweg de Vries equation which can be re-expressed as a linear superposition of either the sum or the difference of two hyperbolic pulse solutions or the sum of either a two-kink or a kink and an antikink solution. These results demonstrate that the notion of superposed solutions extends to coupled nonlinear equations as well.

97 MATHEMATICS AND COMPUTING↗

Vortex methods for separated flows

The numerical solution of the Euler or Navier-Stokes equations by Lagrangian vortex methods is discussed. The mathematical background is presented and includes the relationship with traditional point-vortex studies, convergence to smooth solutions of the Euler equations, and the essential differences between two and three-dimensional cases. The difficulties in extending the method to viscous or compressible flows are explained. Two-dimensional flows around bluff bodies are emphasized. Robustness of the method and the assessment of accuracy, vortex-core profiles, time-marching schemes, numerical dissipation, and efficient programming are treated. Operation counts for unbounded and periodic flows are given, and two algorithms designed to speed up the calculations are described.

Spalart, Philippe R.↗

Large Deformation Behavior of Long Shallow Cylindrical Composite Panels

An exact solution is presented for the large deformation response of a simply supported orthotropic cylindrical panel subjected to a uniform line load along a cylinder generator. The cross section of the cylinder is circular and deformations up to the fully snapped through position are investigated. The orthotropic axes are parallel to the generator and circumferential directions. The governing equations are derived using laminated plate theory, nonlinear strain-displacement relations, and applying variational principles. The response is investigated for the case of a panel loaded exactly at midspan and for a panel with the load offset from midspan. The mathematical formulation is one dimensional in the circumferential coordinate. Solutions are obtained in closed-form. An experimental apparatus was designed to load the panels. Experimental results of displacement controlled tests performed on graphite-epoxy curved panels are compared with analytical predictions.

Carper, Douglas M.↗

Code Verification and Solution Verification framework in pin-resolved neutron transport code MPACT

Program verification in scientific computing encompasses the application of formal and mathematical techniques to a scientific computing code for its credibility, accuracy, and validity. Code Verification identifies bugs and performance issues in the software development stage. Solution Verification assesses the applicability of the code and the accuracy of the solution to problems of interest. Both activities utilize application cases and quantify the error against prescribed acceptance criteria. However, simply executing more application cases does not guarantee stronger or more comprehensive credibility. Here, we establish a verification framework that involves Code Verification and Solution Verification, both of which work together such that the overarching goal of “converge to the correct answer for the intended application” can be reasonably inferred. The application of such a verification framework is demonstrated using the pin-resolved neutron transport code MPACT, where standard unit tests and regression tests are covered, and where the Method of Exact Solutions and the Method of Manufactured Solutions are successfully used. Additionally, the applicability of Method of Manufactured Solutions is extended to the OECD/NEA C5G7 benchmark problems of practical material and geometric configurations. Solution Verification activities are demonstrated on a practical hierarchy of application models of increasing complexity ranging from 2D pin cell problems to 3D assembly problems. The convergence behavior and rate of convergence with respect to each individual variable are studied and provided. This framework can be adapted broadly to other fields involving scientific computing codes.

97 MATHEMATICS AND COMPUTING↗

Backscattering and Nonparaxiality Arrest Collapse of Damped Nonlinear Waves

The critical nonlinear Schrodinger equation (NLS) models the propagation of intense laser light in Kerr media. This equation is derived from the more comprehensive nonlinear Helmholtz equation (NLH) by employing the paraxial approximation and neglecting the backscattered waves. It is known that if the input power of the laser beam (i.e., L(sub 2) norm of the initial solution) is sufficiently high, then the NLS model predicts that the beam will self-focus to a point (i.e.. collapse) at a finite propagation distance. Mathematically, this behavior corresponds to the formation of a singularity in the solution of the NLS. A key question which has been open for many years is whether the solution to the NLH, i.e., the 'parent' equation, may nonetheless exist and remain regular everywhere, in particular for those initial conditions (input powers) that lead to blowup in the NLS. In the current study, we address this question by introducing linear damping into both models and subsequently comparing the numerical solutions of the damped NLH (boundary-value problem) with the corresponding solutions of the damped NLS (initial-value problem). Linear damping is introduced in much the same way as done when analyzing the classical constant-coefficient Helmholtz equation using the limiting absorption principle. Numerically, we have found that it provides a very efficient tool for controlling the solutions of both the NLH and NHS. In particular, we have been able to identify initial conditions for which the NLS solution does become singular. whereas the NLH solution still remains regular everywhere. We believe that our finding of a larger domain of existence for the NLH than that for the NLS is accounted for by precisely those mechanisms, that have been neglected when deriving the NLS from the NLH, i.e., nonparaxiality and backscattering.

Fibich, G.↗

Solution of the wave equation for open surfaces involving a line integral over the edge

A simple mathematical model of a stationary source distribution for the supersonic-propeller noise-prediction formula of Farassat (1983) is developed to test the validity of the formula solutions. The conventional thickness source term is used in place of the Isom thickness formula; the relative importance of the line and surface integrals in the solutions is evaluated; and the numerical results are compared with those obtained with a conventional retarded-time solution in tables. Good agreement is obtained over elevation angles from 10 to 90 deg, and the line-integral contribution is found to be significant at all elevation angles and of the same order of magnitude as the surface-integral contribution at angles less than 30 deg. The amplitude-normalized directivity patterns for the four cases computed (x = 1.5 or 10; k = 5.0 or 50) are presented graphically.

Farassat, F.↗

Nonlinear rotordynamics analysis

Three analytic consequences of the nonlinear Jeffcott equations are examined. The primary application of these analyses is directed toward understanding the excessive vibrations recorded in the Liquid Oxygen (LOX) pump of the Space Shuttle Main Engine (SSME) during hot firing ground testing. The first task is to provide bounds on the coefficients of the equations which delimit the two cases of numerical solution as a circle or an annulus. The second task examines the mathematical generalization to multiple forcing functions, which includes the special problems of mass imbalance, side force, rubbing, and combination of these forces. Finally, stability and boundedness of the steady-state solutions is discussed and related to the corresponding linear problem.

Day, W. B.↗

Solution of Nonlinear Systems

The problem of solving systems of nonlinear equations has been relatively neglected in the mathematical literature, especially in the textbooks, in comparison to the corresponding linear problem. Moreover, treatments that have an appearance of generality fail to discuss the nature of the solutions and the possible pitfalls of the methods suggested. Probably it is unrealistic to expect that a unified and comprehensive treatment of the subject will evolve, owing to the great variety of situations possible, especially in the applied field where some requirement of human or mechanical efficiency is always present. Therefore we attempt here simply to pose the problem and to describe and partially appraise the methods of solution currently in favor.

Turner, L. R.↗

Steady hydromagnetic flows in open magnetic fields. I - A class of analytic solutions

In the case of an establishment of theoretical models of the hydromagnetic solar wind, the inclusion of the effects of the magnetic field in the solar wind makes it extremely dificult to solve the mathematical problem. This paper has the objective to present a set of particular analytic solutions. The general formulation of Tsinganos (1982) is used to identify a class of analytic solutions to the equations of steady hydromagnetic flows in spherical coordinates. Flow in an open magnetic field are studied, taking into account the problem in dimensionless form, the special case of radial flows with alpha = 0, general radial flows, illustrative examples for flows in which alpha is not equal to 0, a parametric study of nonradial flows in which alpha is not equal to zero, variations in the parameter nu, and variations in the initial speed eta.

Low, B. C.↗

Thermodynamically consistent algorithms for models of incompressible multiphase polymer solutions with a variable mobility

Here we present a general strategy for developing structure and property preserving numerical algorithms for thermodynamically consistent models of incompressible multiphase polymer solutions with a variable mobility. We first present a formalism to derive thermodynamically consistent, incompressible, multiphase polymer models. Then, we develop the general strategy, known as the supplementary variable method, to devise thermodynamically consistent numerical approximations to the models. We illustrate the numerical strategy using newly developed models of incompressible diblock copolymer solutions coupled with an electric and a magnetic field, respectively. Mesh refinement is conducted to verify convergence rates of the developed schemes. Some numerical examples are given to exhibit underlying dynamics absent from and driven by the external fields, respectively, highlighting differences between models with the variable and constant mobilities.

97 MATHEMATICS AND COMPUTING↗

Numerical Solution Of Navier-Stokes Equations

Three-dimensional, incompressible flow simulated on computer. Report discusses mathematical basis of IN3SD computer code and application of code to several test problems. Relies heavily on use of pseudocompressibility to solve numerically three-dimensional Navier-Stokes equations of viscous, incompressible flow. With further study and some modification, valuable computational tool for solution of engineering problems in fluid dynamics.

Rogers, Stuart E.↗

Vortex methods for separated flows

The numerical solution of the Euler or Navier-Stokes equations by Lagrangian vortex methods is discussed. The mathematical background is presented in an elementary fashion and includes the relationship with traditional point-vortex studies, the convergence to smooth solutions of the Euler equations, and the essential differences between two- and three-dimensional cases. The difficulties in extending the method to viscous or compressible flows are explained. The overlap with the excellent review articles available is kept to a minimum and more emphasis is placed on the area of expertise, namely two-dimensional flows around bluff bodies. When solid walls are present, complete mathematical models are not available and a more heuristic attitude must be adopted. The imposition of inviscid and viscous boundary conditions without conformal mappings or image vortices and the creation of vorticity along solid walls are examined in detail. Methods for boundary-layer treatment and the question of the Kutta condition are discussed. Practical aspects and tips helpful in creating a method that really works are explained. The topics include the robustness of the method and the assessment of accuracy, vortex-core profiles, timemarching schemes, numerical dissipation, and efficient programming. Calculations of flows past streamlined or bluff bodies are used as examples when appropriate.

Spalart, Philippe R.↗

SANS quantification of bound water in water-soluble polymers across multiple concentration regimes

Contrast-variation small-angle neutron scattering (CV-SANS) is a widely used technique for quantifying hydration water in soft matter systems, but it is predominantly applied in the dilute regime or for systems with a well-defined structure factor. Here, CV-SANS was used to quantify the number of hydration water molecules associating with three water-soluble polymers with different critical solution temperatures and types of water–solute interactions in dilute, semidilute, and concentrated solution through the exploration of novel methods of data fitting and analysis. Multiple SANS fitting workflows with varying levels of model assumptions were evaluated and compared to give insight into SANS model selection. These fitting pathways ranged from general, model-free algorithms to more standard form and structure factor fitting. In addition, Monte Carlo bootstrapping was evaluated as a method to estimate parameter uncertainty through simulation of technical replicates. The most robust fitting workflow for dilute solutions was found to be form factor fitting without CV-SANS ( i.e. polymer in 100% D 2 O). For semidilute and concentrated solutions, while the model-free approach can be mathematically defined for CV-SANS data, the addition of a structure factor imposes physical constraints on the optimization problem, suggesting that the optimal fitting pathway should include appropriate form and structure factor models. The measured hydration numbers were consistent with the number of tightly bound water molecules associated with each monomer unit, and the concentration dependence of the hydration number was largely governed by the chemistry-specific interactions between water and polymer. Polymers with weaker water–polymer interactions ( i.e. those with fewer hydration water molecules) were found to have more bound water at higher concentrations than those with stronger water–polymer interactions due to the increase in the number of forced water–polymer contacts in the concentrated system. This SANS-based method to count hydration water molecules can be applied to polymers in any concentration regime, which will lead to improved understanding of water–polymer interactions and their impact on materials design.

36 MATERIALS SCIENCE↗

Convolutional Neural Networks for Problems in Transport Phenomena: A Theoretical Minimum

Convolutional neural network (CNN), a deep learning algorithm, has gained popularity in technological applications that rely on interpreting images (typically, an image is a 2D field of pixels). Transport phenomena is the science of studying different fields representing mass, momentum, or heat transfer. Some of the common fields are species concentration, fluid velocity, pressure, and temperature. Each of these fields can be expressed as an image(s). Consequently, CNNs can be leveraged to solve specific scientific problems in transport phenomena. Herein, we show that such problems can be grouped into three basic categories: (a) mapping a field to a descriptor (b) mapping a field to another field, and (c) mapping a descriptor to a field. After reviewing the representative transport phenomena literature for each of these categories, we illustrate the necessary steps for constructing appropriate CNN solutions using sessile liquid drops as an exemplar problem. If sufficient training data is available, CNNs can considerably speed up the solution of the corresponding problems. Finally, the present discussion is meant to be minimalistic such that readers can easily identify the transport phenomena problems where CNNs can be useful as well as construct and/or assess such solutions.

97 MATHEMATICS AND COMPUTING↗

The structure of synchronously rotating close binaries built on polytropic model nu = 3

Mathematical procedures are developed which can provide an exact solution for the structure of synchronously rotating close binaries built on the polytropic model with nu equal to 3. Results are discussed for separations of centers of 10.0, 4.0, 3.0, and 2.8 solar radii. For the present discussion, both components are assumed to have the solar mass and, in the zeroth approximation, the solar radius. The present results are compared with those obtained by Chandrasekhar (1933) using first-order perturbation theory.

Green, L. C.↗

On estimating gravity anomalies: A comparison of least squares collocation with least squares techniques

The least squares collocation algorithm for estimating gravity anomalies from geodetic data is shown to be an application of the well known regression equations which provide the mean and covariance of a random vector (gravity anomalies) given a realization of a correlated random vector (geodetic data). It is also shown that the collocation solution for gravity anomalies is equivalent to the conventional least-squares-Stokes' function solution when the conventional solution utilizes properly weighted zero a priori estimates. The mathematical and physical assumptions underlying the least squares collocation estimator are described, and its numerical properties are compared with the numerical properties of the conventional least squares estimator.

Argentiero, P.↗