Search NASA⌕ Search

SEARCH · Search NASA

Results for “Newton Method”

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 217 records · Page 12

Practical Aspects of Krylov Subspace Iterative Methods in CFD

Implementation issues associated with the application of Krylov subspace iterative methods, such as Newton-GMRES, are presented within the framework of practical computational fluid dynamic (CFD) applications. This paper categorizes, evaluates, and contrasts the major ingredients (function evaluations, matrix-vector products, and preconditioners) of Newton-GMRES Krylov subspace methods in terms of their effect on the local linear and global nonlinear convergence, memory requirements, and accuracy. The discussion focuses on Newton-GMRES in both a structured multi-zone incompressible Navier-Stokes solver and an unstructured mesh finite-volume Navier-Stokes solver. Approximate versus exact matrix-vector products, effective preconditioners, and other pertinent issues are addressed.

Pulliam, Thomas H.↗

The exact numerical calculation of propeller noise

This paper presents a method for computing the acoustic pressure signature of a propeller due to blade thickness and steady surface pressure loading including the effect of forward motion. The tip helical speed is restricted to subsonic. The formulation is based on the Ffowcs Williams-Hawkings exact result concerning the sound generation by moving bodies. The blade surface is described in a helical coordinate system by the blade angles, positions of the leading and trailing edges along the local pitch helix, and the camber and thickness distributions along the chord. The retarded distance of the moving acoustic source on an element of surface area of the propeller blade is obtained numerically by Newton's method. The surface integration over the propeller blades is carried out in the helical coordinate system using a double Gauss-Legendre formula. Good agreement was found between analytic and experimental results, both in the time and frequency domains. Some samples showing the effects of propeller geometry and tip Mach number on the acoustic signature are included.

Woan, C. J.↗

Global Asymptotic Behavior of Iterative Implicit Schemes

The global asymptotic nonlinear behavior of some standard iterative procedures in solving nonlinear systems of algebraic equations arising from four implicit linear multistep methods (LMMs) in discretizing three models of 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations (ODEs) is analyzed using the theory of dynamical systems. The iterative procedures include simple iteration and full and modified Newton iterations. The results are compared with standard Runge-Kutta explicit methods, a noniterative implicit procedure, and the Newton method of solving the steady part of the ODEs. Studies showed that aside from exhibiting spurious asymptotes, all of the four implicit LMMs can change the type and stability of the steady states of the differential equations (DEs). They also exhibit a drastic distortion but less shrinkage of the basin of attraction of the true solution than standard nonLMM explicit methods. The simple iteration procedure exhibits behavior which is similar to standard nonLMM explicit methods except that spurious steady-state numerical solutions cannot occur. The numerical basins of attraction of the noniterative implicit procedure mimic more closely the basins of attraction of the DEs and are more efficient than the three iterative implicit procedures for the four implicit LMMs. Contrary to popular belief, the initial data using the Newton method of solving the steady part of the DEs may not have to be close to the exact steady state for convergence. These results can be used as an explanation for possible causes and cures of slow convergence and nonconvergence of steady-state numerical solutions when using an implicit LMM time-dependent approach in computational fluid dynamics.

Yee, H. C.↗

Numerical implementation of the MOOSE subchannel module (SCM) algorithm

The subchannel module of moose (SCM), previously referred to as Pronghorn-SC, is a subchannel code designed to resolve single-phase flow fields and calculate the relevant flow variables, in nuclear reactor fuel-pin assemblies. The assemblies it models are: water-cooled, with bare fuel-pins in a square ducted, quadrilateral lattice and liquid-metal/water-cooled, with bare/wire-wrapped fuel pins in a hexagonal ducted, triangular lattice. Previous publications have presented the development, validation and verification of SCM. This work presents an comparative overview of the different solvers implemented within SCM. This includes the solver algorithm description, and a performance comparison. The two test cases chosen to demonstrate the solver performance were taken from the PSBT enthalpy mixing benchmark and the ORNL-19 pin benchmark.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Learning and discovering multiple solutions using physics-informed neural networks with random initialization and deep ensemble

In this work we explore the capability of physics-informed neural networks (PINNs) to discover multiple solutions. Many real-world phenomena governed by nonlinear differential equations (DEs), such as fluid flow, exhibit multiple solutions under the same conditions, yet capturing this solution multiplicity remains a significant challenge. A key difficulty lies in providing appropriate initial conditions or guesses, as widely used time-marching schemes and Newton’s method are highly sensitive to these choices when solving complex computational problems. While machine learning models, particularly PINNs, have shown promise in solving DEs, their ability to capture multiple solutions remains underexplored. In this work, we propose a simple and practical approach using PINNs to learn and discover multiple solutions. We first demonstrate that PINNs, when combined with random initialization and deep ensemble method—originally developed for uncertainty quantification—can effectively uncover multiple solutions to nonlinear ordinary and partial DEs. Although training large ensembles of PINNs may appear computationally demanding, this can be done efficiently using vectorization techniques supported by modern deep learning frameworks, allowing many networks to be trained simultaneously. Our approach highlights the critical role of initialization in shaping solution diversity, addressing an often-overlooked aspect of machine learning for scientific computing. Furthermore, we propose utilizing PINN-generated solutions as initial conditions or initial guesses for conventional numerical solvers to enhance accuracy and efficiency in capturing multiple solutions. Extensive numerical experiments, including the Allen–Cahn equation and cavity flow, where our approach successfully identifies both stable and unstable solutions, validate the effectiveness of our method. These findings establish a general and efficient framework for addressing solution multiplicity in nonlinear DEs.

97 MATHEMATICS AND COMPUTING↗

Computer program for network synthesis by frequency response fit

Computer program synthesizes a passive network by minimizing the difference in desired and actual frequency response. The program solves for the critical points of the error function /weighted least squares fit between calculated and desired frequency response/ by the multivariable Newton-Raphson method with components constrained to an admissible region.

Green, S.↗

Analysis of an arched outer-race ball bearing considering centrifugal forces

A Newton-Raphson method of iteration was used in evaluating the radial and axial projection of the distance between the ball center and the outer raceway groove curvature center (V and W). Fatigue life evaluations were made. The similar analysis of a conventional bearing can be directly obtained from the arched bearing analysis by simply letting the amount of arching be zero (g = 0) and not considering equations related to the unloaded half of the outer race. The analysis was applied to a 150-mm angular contact ball bearing. Results for life, contact loads, and angles are shown for a conventional bearing (g = 0) and two arched bearings (g = 0.127 mm (0.005 in.), and 0.254 mm (0.010 in.)). The results indicate that an arched bearing is highly desirable for high speed applications. In particular, for a DN value of 3 million (20,000 rpm) and an applied axial load of 4448 N (1000 lb), an arched bearing shows an improvement in life of 306 percent over that of a conventional bearing. At 4.2 million DN (28,000 rpm), the corresponding improvement is 340 percent. It was also found for low speeds, the arched bearing does not offer the advantages that it does for high speed applications.

Hamrock, B. J.↗

Analysis of an arched outer-race ball bearing considering centrifugal forces.

Thrust-load analysis of a 150-mm angular contact ball bearing, taking into account centrifugal forces but omitting gyroscopics, elastohydrodynamics, and thermal effects. A Newton-Raphson method of iteration is used to evaluate the radial and axial projection of the distance between the ball center and the outer raceway groove curvature center. Fatigue life of the bearing is evaluated. Results for life, contact loads, and angles are given for a conventional bearing and two arched bearings.

Hamrock, B. J.↗

Thermochemical nonequilibrium in atomic hydrogen at elevated temperatures

A numerical study of the nonequilibrium flow of atomic hydrogen in a cascade arc was performed to obtain insight into the physics of the hydrogen cascade arc. A rigorous mathematical model of the flow problem was formulated, incorporating the important nonequilibrium transport phenomena and atomic processes which occur in atomic hydrogen. Realistic boundary conditions, including consideration of the wall electrostatic sheath phenomenon, were included in the model. The governing equations of the asymptotic region of the cascade arc were obtained by writing conservation of mass and energy equations for the electron subgas, an energy conservation equation for heavy particles and an equation of state. Finite-difference operators for variable grid spacing were applied to the governing equations and the resulting system of strongly coupled, stiff equations were solved numerically by the Newton-Raphson method.

Scott, R. K.↗

Extraction from flight data of longitudinal aerodynamic coefficients for F-8 aircraft with supercritical wing

The longitudinal aerodynamic derivatives of the F-8 aircraft with supercritical wing were obtained from flight data by a parameter-extraction algorithm at Mach numbers of 0.8, 0.9, and 0.98. A set of derivatives were obtained from which calculated aircraft responses were correlated almost identically with actual flight responses. In general, the trends of the extracted derivatives obtained by the algorithm agreed with those obtained by a Newton-Raphson method and with preliminary data from the Langley 8-foot transonic pressure tunnel. The wind-tunnel damping derivatives were, however, substantially higher than the converged damping derivatives possibly because of Reynolds number differences between flight and model tests.

Williams, J. L.↗

A modified Newton-Raphson analysis of flight measurements of the trailing vortices of a heavy jet transport

A modified Newton-Raphson method has been used to determine the parameters in equations describing the vortex flow to obtain a best match with flight measurements of the flow behind the C-5A airplane. The flight measurements were made using a specially instrumented T-33 airplane which passed as closely as possible to the centers of the trailing vortices at several distances behind the C-5A airplane. The flight measurements were transformed to flow velocity relative to an inertial frame of reference. The assumed form of the flow consisted of the superposition of two counterrotating, finite core vortices. The positions of the vortex centers, their total circulation, the effective eddy viscosity and measurement bias were the parameters adjusted. The assumed form of vortex flow fit well the measured velocities for the numerous sets of data, both flaps up and down for the C-5A airplane. The resulting values of total circulation, however, were about two-thirds that expected of a wing with an elliptical loading. A partial explanation of the less than expected circulation is a dip in the spanwise lift distribution at the airplane's center line. The distance between the trailing vortices at the smallest times encountered is somewhat less than that expected for an elliptical wing loading.

Taylor, L. W., Jr.↗

Evaluation of importance of lateral acceleration derivatives in extraction of lateral-directional derivatives at high angles of attack

A theoretical investigation was conducted to determine the importance of the lateral acceleration (beta) derivatives in the extraction of lateral-directional stability derivatives for swept wing airplanes at high angles of attack. Representative values of lateral acceleration derivatives in yaw and roll (Cn beta and Cl beta) were used in a computer program to generate representative flight motions at several angles of attack and altitudes. The computer-generated motions were then subjected to a parameter identification process based on a modified Newton-Raphson method. Two identification techniques were evaluated, one which included the beta derivatives and one which neglected them. The results of the study indicate that omission of the beta derivatives from mathematical models used in the derivative-extraction techniques can produce erroneous values for the lateral-directional stability derivatives particularly at high angles of attack, where the beta derivatives are large. The largest errors occur in the dynamic derivatives, but large errors may also occur in the static derivatives for cases in which the beta derivatives have large effects on the flight motions of the airplane. In addition, the resulting identified mathematical models provide poor motion prediction as well as erroneous predictions of dynamic modal characteristics. These results strongly indicate that the effects of beta derivatives should be considered in any attempt to extract lateral-directional aerodynamic parameters at high angles of attack.

Nguyen, L. T.↗

Preliminary design of composite wings for buckling, strength and displacement constraints

An unstiffened panel buckling constraint for balanced, symmetric laminated composites is included on the global design level in a mathematical programming structural optimization procedure for designing wing structures. Constraints are introduced by penalty functions, and Newton's method based on approximate second derivatives of the penalty terms is used as the search algorithm to obtain minimum-mass designs. Constraint approximations used during the optimization process contribute to the computational efficiency of the procedure. A criterion is developed that identifies the appropriate conservative form of the constraint approximations that are used with the optimization procedure. Minimum-mass design results are obtained for a multispar high-aspect-ratio wing subjected to material strength, minimum-gage, displacement, panel buckling and twist constraints. The material systems considered for the examples are all graphite-epoxy, graphite-epoxy with boron-epoxy spar caps, and all aluminum. The composite material designs are shown to have an advantage over the aluminum designs since they can often satisfy additional constraints with only small mass increases.

Starnes, J. H., Jr.↗

NEWSUMT: A FORTRAN program for inequality constrained function minimization, users guide

A computer program written in FORTRAN subroutine form for the solution of linear and nonlinear constrained and unconstrained function minimization problems is presented. The algorithm is the sequence of unconstrained minimizations using the Newton's method for unconstrained function minimizations. The use of NEWSUMT and the definition of all parameters are described.

Miura, H.↗