Search NASA⌕ Search

SEARCH · Search NASA

Results for “Integral solution”

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 397 records · Page 22

Airframe Noise Prediction by Acoustic Analogy: Revisited

The present work follows a recent survey of airframe noise prediction methodologies. In that survey, Lighthill s acoustic analogy was identified as the most prominent analytical basis for current approaches to airframe noise research. Within this approach, a problem is typically modeled with the Ffowcs Williams and Hawkings (FW-H) equation, for which a geometry-independent solution is obtained by means of the use of the free-space Green function (FSGF). Nonetheless, the aeroacoustic literature would suggest some interest in the use of tailored or exact Green s function (EGF) for aerodynamic noise problems involving solid boundaries, in particular, for trailing edge (TE) noise. A study of possible applications of EGF for prediction of broadband noise from turbulent flow over an airfoil surface and the TE is, therefore, the primary topic of the present work. Typically, the applications of EGF in the literature have been limited to TE noise prediction at low Mach numbers assuming that the normal derivative of the pressure vanishes on the airfoil surface. To extend the application of EGF to higher Mach numbers, the uniqueness of the solution of the wave equation when either the Dirichlet or the Neumann boundary condition (BC) is specified on a deformable surface in motion. The solution of Lighthill s equation with either the Dirichlet or the Neumann BC is given for such a surface using EGFs. These solutions involve both surface and volume integrals just like the solution of FW-H equation using FSGF. Insight drawn from this analysis is evoked to discuss the potential application of EGF to broadband noise prediction. It appears that the use of a EGF offers distinct advantages for predicting TE noise of an airfoil when the normal pressure gradient vanishes on the airfoil surface. It is argued that such an approach may also apply to an airfoil in motion. However, for the prediction of broadband noise not directly associated with a trailing edge, the use of EGF does not appear to offer any advantages over the use of FSGF at the present stage of development. It is suggested here that the applications of EGF for airframe noise analysis be continued. As an example pertinent to airframe noise prediction, the Fast Scattering Code of NASA Langley is utilized to obtain the EGF numerically on the surface of a three dimensional wing with a flap and leading edge slat in uniform rectilinear motion. The interpretation and use of these numerical Green functions are then discussed.

Farassat, F.↗

An integration scheme to determine the dynamic response of a launch vehicle with several payloads

The coupled equations of motion of a system composed of a launch vehicle and multiple payloads are derived. In the process, it is shown how superfluous interface degrees of freedom on the booster side can be accommodated in the formulation. The discrete system equations are directly integrated, avoiding the solution of an expensive system eigenvalue problem. A modified Newmark-Chan-Beta numerical integration scheme is used to obtain the response. The unique form of the equations of motion allow for a quick and accurate solution both for the response and the internal loads. Although completely general, the presented technique is most effective when the number of interface degrees of freedom is relatively small compared to the overall number of degrees of freedom of the system. The technique is applied to the case of the S.T.S.-S.T.-OMS Kit system. Finally, an evaluation of the presented method is included.

Engels, R. C.↗

Enabling a Weight Efficient Power System for an Electrified Turbofan Through Gearbox Design and Control

Electrified Aircraft Propulsion (EAP) concepts could enable various benefits through a variety of use-cases. Benefits are sought in the form of reduced fuel burn and emissions. Many of these concepts involve the electrification of gas turbine engines. Electrification is accomplished through the integration of electric machines (EMs) with the engine shafts, providing the ability to inject and/or extract power as desired. Energy storage is also a common feature. An often-overlooked feature of the analysis is the means of integrating the EMs with the engine shafts. A trivial solution is to allow each shaft to have its own dedicated EM through independent geartrains. However, non-trivial mechanical integration solutions could provide benefits when considered in coordination with control logic to manage the operation of the propulsion system. Here, such a solution is considered and is shown to have the potential to reduce weight for a relevant conceptual electrified propulsion system. Weight saving benefits are demonstrated for various EAP use-cases. In particular, the mild hybrid application was shown to benefit from a power system weight reduction of 47% and an overall system weight reduction of 29%.

Hybrid Electric Propulsion↗

Enabling a Weight Efficient Power System for an Electrified Turbofan Through Gearbox Design and Control

Electrified Aircraft Propulsion (EAP) concepts could enable various benefits through a variety of use-cases. Benefits are sought in the form of reduced fuel burn and emissions. Many of these concepts involve the electrification of gas turbine engines. Electrification is accomplished through the integration of electric machines (EMs) with the engine shafts, providing the ability to inject and/or extract power as desired. Energy storage is also a common feature. An often-overlooked feature of the analysis is the means of integrating the EMs with the engine shafts. A trivial solution is to allow each shaft to have its own dedicated EM through independent geartrains. However, non-trivial mechanical integration solutions could provide benefits when considered in coordination with control logic to manage the operation of the propulsion system. Here, such a solution is considered and is shown to have the potential to reduce weight for a relevant conceptual electrified propulsion system. Weight saving benefits are demonstrated for various EAP use-cases. In particular, the mild hybrid application was shown to benefit from a power system weight reduction of 47% and an overall system weight reduction of 29%.

Hybrid Electric Propulsion↗

Numerical solution for subcritical flows by a transonic integral equation method

The method presented makes use of a division of the region of integration into closed rectangular elements. The velocity is taken to be constant in each element. The integral equation is reduced to a matrix equation which can be solved by an appropriate iteration approach. The derivation and solution of the matrix equation are discussed and the matrix elements are considered. The described concepts were implemented for a nonlifting parabolic-arc airfoil.

Ogana, W.↗

Numerical methods of solving a system of multi-dimensional nonlinear equations of the diffusion type

The principles of conservation and stability of difference schemes achieved using the iteration control method were examined. For the schemes obtained of the predictor-corrector type, the conversion was proved for the control sequences of approximate solutions to the precise solutions in the Sobolev metrics. Algorithms were developed for reducing the differential problem to integral relationships, whose solution methods are known, were designed. The algorithms for the problem solution are classified depending on the non-linearity of the diffusion coefficients, and practical recommendations for their effective use are given.

Agapov, A. V.↗

Three-dimensional analysis of chevron-notched specimens by boundary integral method

The chevron-notched short bar and short rod specimens was analyzed by the boundary integral equations method. This method makes use of boundary surface elements in obtaining the solution. The boundary integral models were composed of linear triangular and rectangular surface segments. Results were obtained for two specimens with width to thickness ratios of 1.45 and 2.00 and for different crack length to width ratios ranging from 0.4 to 0.7. Crack opening displacement and stress intensity factors determined from displacement calculations along the crack front and compliance calculations were compared with experimental values and with finite element analysis.

Mendelson, A.↗

Preliminary Goddard geopotential using optical tracking data and a comparison with SAO models

A preliminary Goddard Space Flight Center (GSFC) geopotential and center of mass station coordinate solution was obtained from satellite orbital data using numerical integration theory. This geodetic solution is a prelude to a more general solution which will combine the 1971 International Satellite Geodesy Experiment (ISAGEX) laser data with the present data being employed. The present GSFC geopotential solution consists of the spherical harmonic coefficients through degree and order eight with higher order satellite resonant coefficients. The solution represents a first iteration result from 17 satellites with approximately 150 weekly orbital arcs containing some 40,000 optical observations. The GSFC preliminary result is compared with final results from the Smithsonian Astrophysical Observatory (SAO) solutions including the 1969 SAO Standard Earth II solution. One aspect of interest for the comparison is that SAO uses an analytic theory for the orbital solution whereas GSFC uses a numerical integration theory. The comparison of geopotential results shows that good agreement exists in general but that there are some areas of minor differences.

Lerch, F. J.↗

Parallel-in-Time Solution of Allen-Cahn Equations by Integrating Operator Learning into the Parareal Method

While recent advances in deep learning have shown promising efficiency gains in solving time-dependent partial differential equations (PDEs), matching the accuracy of conventional numerical solvers still remains a challenge. One strategy to improve the accuracy of deep learning-based solutions for time-dependent PDEs is to use the learned model as the coarse propagator in the Parareal method and a traditional numerical method as the fine solver. However, successful integration of deep learning into the Parareal method requires consistency between the coarse and fine solvers, particularly for PDEs exhibiting rapid changes such as sharp transitions. Here, to ensure this consistency, we propose using convolutional neural networks (CNNs) to learn the fully discrete time-stepping operator defined by the same numerical scheme employed as the fine solver. We demonstrate the effectiveness of the proposed method in solving the classical and mass-conservative Allen–Cahn (AC) equations. Through iterative updates in the Parareal algorithm, our approach achieves a significant computational speedup compared to traditional fine solvers while converging to high-accuracy solutions. Our results highlight that the proposed hybrid Parareal algorithm effectively accelerates simulations, particularly when implemented on multiple GPUs, and converges to the desired accuracy in only a few iterations. Another advantage of our method is that the CNN model is trained on trajectory-based data generated from random initial conditions, such that the trained model can be used to solve the AC equations with various initial conditions without retraining. This work demonstrates the potential of integrating neural network methods into parallel-in-time frameworks for efficient and accurate simulations of time-dependent PDEs.

97 MATHEMATICS AND COMPUTING↗

Electromagnetic inverse scattering

A three-dimensional electromagnetic inverse scattering identity, based on the physical optics approximation, is developed for the monostatic scattered far field cross section of perfect conductors. Uniqueness of this inverse identity is proven. This identity requires complete scattering information for all frequencies and aspect angles. A nonsingular integral equation is developed for the arbitrary case of incomplete frequence and/or aspect angle scattering information. A general closed-form solution to this integral equation is developed, which yields the shape of the scatterer from such incomplete information. A specific practical radar solution is presented. The resolution of this solution is developed, yielding short-pulse target resolution radar system parameter equations. The special cases of two- and one-dimensional inverse scattering and the special case of a priori knowledge of scatterer symmetry are treated in some detail. The merits of this solution over the conventional radar imaging technique are discussed.

Bojarski, N. N.↗

Effect of friction in wedging of elastic solids

In this paper the contact problem for an elastic wedge of arbitrary angle is considered. It is assumed that the external load is applied to the medium through a rigid wedge and the coefficient of friction between the loading wedge and the elastic solid is constant. The problem is reduced to a singular integral equation of the second kind with the contact pressure as the unknown function. An effective numerical solution of the integral equation is described and the results of three examples are presented. The comparison of these results with those obtained from the frictionless wedge problem indicates that generally friction has the tendency of reducing the peak values of the stress intensity factors calculated at the wedge apex and at the end points of the contact area.

Erdogan, F.↗

Inverse scattering for an exterior Dirichlet program

Scattering due to a metallic cylinder which is in the field of a wire carrying a periodic current is considered. The location and shape of the cylinder is obtained with a far field measurement in between the wire and the cylinder. The same analysis is applicable in acoustics in the situation that the cylinder is a soft wall body and the wire is a line source. The associated direct problem in this situation is an exterior Dirichlet problem for the Helmholtz equation in two dimensions. An improved low frequency estimate for the solution of this problem using integral equation methods is presented. The far field measurements are related to the solutions of boundary integral equations in the low frequency situation. These solutions are expressed in terms of mapping function which maps the exterior of the unknown curve onto the exterior of a unit disk. The coefficients of the Laurent expansion of the conformal transformations are related to the far field coefficients. The first far field coefficient leads to the calculation of the distance between the source and the cylinder.

Hariharan, S. I.↗

Effects of gravity-gradient torque on the rotational motion of a triaxial satellite in a precessing elliptic orbit.

A method of general perturbations, based on the use of Lie series to generate approximate canonical transformations, is applied to study the effects of gravity-gradient torque on the rotational motion of a triaxial, rigid satellite. The center of mass of the satellite is constrained to move in an elliptic orbit about an attracting point mass. The orbit, which has a constant inclination, is free to precess and spin. The method of general perturbations is used to obtain the Hamiltonian for the nonresonant secular and long-period rotational motion of the satellite. The differential equations derivable from the transformed Hamiltonian are integrable, and the solution for the long-term motion may be expressed in terms of Jacobian elliptic functions and elliptic integrals. Geometrical aspects of the long-term rotational motion are discussed, and a comparison of theoretical results with observations is made.

Cochran, J. E.↗

Long-term motion in a restricted problem of rotational motion

A method of general perturbations, based on the use of Lie series to generate approximate canonical transformations, is applied to study the long-term effects of gravity-gradient torque and orbital evolution on the rotational motion of a triaxial, rigid satellite. The center of mass of the satellite is constrained to move in an elliptic orbit about an attracting point mass. The orbit, which has a constant inclination, is constrained to precess and spin with constant rates. The method of general perturbations is used to obtain the Hamiltonian for the nonresonant secular and long-period rotational motion of the satellite to second order in n/omega sub 0, where n is the orbital mean motion of the center of mass and omega sub 0 is a reference value of the magnitude of the satellite's rotational angular velocity. The differential equations derivable from the transformed Hamiltonian are integrable, and the solution for the long-term motion may be expressed in terms of Jacobian elliptic functions and elliptic integrals. Geometrical aspects of the long-term rotational motion are discussed, and a comparison of theoretical results with observations is made.

Cochran, J. E.↗

A numerical study of transient heat and mass transfer in crystal growth

A numerical analysis of transient heat and solute transport across a rectangular cavity is performed. Five nonlinear partial differential equations which govern the conservation of mass, momentum, energy and solute concentration related to crystal growth in solution, are simultaneously integrated by a numerical method based on the SIMPLE algorithm. Numerical results showed that the flow, temperature and solute fields are dependent on thermal and solutal Grashoff number, Prandtl number, Schmidt number and aspect ratio. The average Nusselt and Sherwood numbers evaluated at the center of the cavity decrease markedly when the solutal buoyancy force acts in the opposite direction to the thermal buoyancy force. When the solutal and thermal buoyancy forces act in the same direction, however, Sherwood number increases significantly and yet Nusselt number decreases. Overall effects of convection on the crystal growth are seen to be an enhancement of growth rate as expected but with highly nonuniform spatial growth variations.

Han, Samuel Bang-Moo↗

A new payload integration method

This paper introduces a new payload integration method. The system equations of motion are derived in terms of interface-restrained booster and payload modes. These system equations are directly integrated, avoiding the solution of an expensive eigenvalue problem. A modified Newmark-Chan-Beta numerical scheme is used to perform this integration. Due to the unique form of the system equations it is possible to achieve significant savings. A similar savings is feasible for the computation of the internal loads. This new approach does not involve approximations. Although completely general, the presented technique is most effective when the number of interface degrees of freedom is relatively small. Both determinate and indeterminate interfaces are allowed.

Engels, R. C.↗