Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary element 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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

DeepONet-Assisted Optimization of Surface Topography for Transition Delay in a Mach 4.5 Boundary Layer

We use deep learning, an ensemble variational technique (EnVar), and direct numerical simulations(DNS) to design an optimal topography for a two-dimensional roughness element that delays the on-set of laminar-turbulent transition in a Mach 4.5 flat-plate boundary layer. Deep operator networks (DeepONets), which have the known ability to learn complex nonlinear operators within dynamical systems, are used for machine learning. For the baseline configuration of a smooth flat plate, the second-mode waves at the DNS inflow cause a quick nonlinear breakdown of the high-speed boundary layer within the computational domain. Results reported in the present study validate the ability of DeepONets to model the transition delay via a given topography of the roughness element. The computing cost to optimize the rough-ness element for minimal skin-friction drag is substantially lowered by the DeepONets-based reduced-order model. In comparison to the baseline method of EnVar optimization based on DNS alone, the DeepONets-based EnVar optimizer is able to delay transition past the outflow boundary of the computational domain while utilizing almost 5–6 times fewer DNS.

Machine Learning↗
DOE OSTI2023

Slip localization behavior at triple junctions in nickel-base superalloys

Incipient slip localization in the vicinity of hundreds of grain boundary triple junctions (TJs) in a lightly deformed nickel-base superalloy IN718 is studied using a combination of three-dimensional (3D) crystal plasticity finite element (CPFE) modeling, high resolution digital image correlation (HR-DIC) and 3D electron back-scatter diffraction tomography (3D EBSD). A 3D reconstruction method enables identification of thousands of TJs and correspondence of any observed slip bands with their originating TJ lines below the specimen surface. Here, we present a large-scale CPFE model of the experimental 3D microstructure composed of high-fidelity representation of the TJ lines and the boundaries and interiors of the parent grains and use it to calculate the local micromechanical response and slip activity of all TJs at the onset of macroscopic yielding. Statistical analysis of the calculated quantities reveal TJs develop larger stress concentration and grain-average re-orientation than grain interiors and grain boundaries, however no substantial differences in cumulative slip were found among these microstructural regions. We find that TJs with observed slip bands generate lower grain-average re-orientation, fewer active slip systems, and more localized slip on a single system than those without. The distinctions in the reorientation and slip activity are stronger in TJs that experience more intense slip.

36 MATERIALS SCIENCE↗
NASA NTRS2010

Numerical Computations of Hypersonic Boundary-Layer over Surface Irregularities

Surface irregularities such as protuberances inside a hypersonic boundary layer may lead to premature transition on the vehicle surface. Early transition in turn causes large localized surface heating that could damage the thermal protection system. Experimental measurements as well as numerical computations aimed at building a knowledge base for transition Reynolds numbers with respect to different protuberance sizes and locations have been actively pursued in recent years. This paper computationally investigates the unsteady wake development behind large isolated cylindrical roughness elements and the scaled wind-tunnel model of the trip used in a recent flight measurement during the reentry of space shuttle Discovery. An unstructured mesh, compressible flow solver based on the space-time conservation element, solution element (CESE) method is used to perform time-accurate Navier-Stokes calculations for the flow past a roughness element under several wind-tunnel conditions. For a cylindrical roughness element with a height to the boundary-layer thickness ratio from 0.8 to 2.5, the wake flow is characterized by a mushroom-shaped centerline streak and horse-shoe vortices. While time-accurate solutions converged to a steady-state for a ratio of 0.8, strong flow unsteadiness is present for a ratio of 1.3 and 2.5. Instability waves marked by distinct disturbance frequencies were found in the latter two cases. Both the centerline streak and the horse-shoe vortices become unstable downstream. The oscillatory vortices eventually reach an early breakdown stage for the largest roughness element. Spectral analyses in conjunction with the computed root mean square variations suggest that the source of the unsteadiness and instability waves in the wake region may be traced back to possible absolute instability in the front-side separation region.

Chang, Chau-Lyan↗
NASA NTRS1999

High-Resolution Genuinely Multidimensional Solution of Conservation Laws by the Space-Time Conservation Element and Solution Element Method

In this overview paper, we review the basic principles of the method of space-time conservation element and solution element for solving the conservation laws in one and two spatial dimensions. The present method is developed on the basis of local and global flux conservation in a space-time domain, in which space and time are treated in a unified manner. In contrast to the modern upwind schemes, the approach here does not use the Riemann solver and the reconstruction procedure as the building blocks. The drawbacks of the upwind approach, such as the difficulty of rationally extending the 1D scalar approach to systems of equations and particularly to multiple dimensions is here contrasted with the uniformity and ease of generalization of the Conservation Element and Solution Element (CE/SE) 1D scalar schemes to systems of equations and to multiple spatial dimensions. The assured compatibility with the simplest type of unstructured meshes, and the uniquely simple nonreflecting boundary conditions of the present method are also discussed. The present approach has yielded high-resolution shocks, rarefaction waves, acoustic waves, vortices, ZND detonation waves, and shock/acoustic waves/vortices interactions. Moreover, since no directional splitting is employed, numerical resolution of two-dimensional calculations is comparable to that of the one-dimensional calculations. Some sample applications displaying the strengths and broad applicability of the CE/SE method are reviewed.

Himansu, Ananda↗
NASA NTRS1985

On the penetration of a hot diapir through a strongly temperature-dependent viscosity medium

The ascent of a hot spherical body through a fluid with a strongly temperature-dependent viscosity has been studied using an axisymmetric finite element method. Numerical solutions range over Peclet numbers of 0.1 - 1000 from constant viscosity up to viscosity variations of 100,000. Both rigid and stress-free boundary conditions were applied at the surface of the sphere. The dependence of drag on viscosity variation was shown to have no dependence on the stress boundary condition except for a Stokes flow scaling factor. A Nusselt number parameterization based on the stress-free constant viscosity functional dependence on the Peclet number scaled by a parameter depending on the viscosity structure fits both stress-free and rigid boundary condition data above viscosity variations of 100. The temperature scale height was determined as a function of sphere radius. For the simple physical model studied in this paper pre-heating is required to reduce the ambient viscosity of the country rock to less than 10 to the 22nd sq cm/s in order for a 10 km diapir to penetrate a distance of several radii.

Daly, S. F.↗
NASA NTRS1987

Component mode synthesis and large deflection vibration of complex structures. Volume 2: Single-mode large deflection vibrations of beams and plates using finite element method

A finite element method is presented for the large amplitude vibrations of complex structures that can be modelled with beam and rectangular plate elements subjected to harmonic excitation. Both inplane deformation and inertia are considered in the formulation. Derivation of the harmonic force and nonlinear stiffness matrices for a beam and a rectangular plate element are presented. Solution procedures and convergence characteristics of the finite element method are described. Nonlinear response to uniform and concentrated harmonic loadings and improved nonlinear free vibration results are presented for beams and rectangular plates of various boundary conditions.

Mei, Chuh↗
NASA NTRS1987

Multiple-mode nonlinear free and forced vibrations of beams using finite element method

Effects of large deflection geometric nonlinearity and multiple modes for free and forced vibrations of beams are investigated. Inplane displacement and inplane inertia are included in the formulation. The finite element method is employed. A harmonic force matrix is derived for forced vibration analysis. The relations of amplitude, frequency ratio and forcing intensity for beams of different boundary conditions and various load conditions are presented.

Mei, Chuh↗
NASA NTRS1988

Analysis of crack propagation in roller bearings using the boundary integral equation method - A mixed-mode loading problem

Crack propagation in a rotating inner raceway of a high-speed roller bearing is analyzed using the boundary integral method. The model consists of an edge plate under plane strain condition upon which varying Hertzian stress fields are superimposed. A multidomain boundary integral equation using quadratic elements was written to determine the stress intensity factors KI and KII at the crack tip for various roller positions. The multidomain formulation allows the two faces of the crack to be modeled in two different subregions, making it possible to analyze crack closure when the roller is positioned on or close to the crack line. KI and KII stress intensity factors along any direction were computed. These calculations permit determination of crack growth direction along which the average KI times the alternating KI is maximum.

Ghosn, L. J.↗
NASA NTRS1989

A mesh regeneration method using quadrilateral and triangular elements for compressible flows

An adaptive remeshing method using both triangular and quadrilateral elements suitable for high-speed viscous flows is presented. For inviscid flows, the method generates completely unstructured meshes. For viscous flows, structured meshes are generated for boundary layers, and unstructured meshes are generated for inviscid flow regions. Examples of inviscid and viscous adaptations for high-speed flows are presented.

Vemaganti, G. R.↗
NASA NTRS1989

Compact finite volume methods for the diffusion equation

The paper describes an approach to treating initial-boundary-value problems by finite volume methods in which the parallel between differential and difference arguments is closely maintained. By using intrinsic geometrical properties of the volume elements, it is possible to describe discrete versions of the div, curl, and grad operators which lead, using summation-by-parts techniques, to familiar energy equations as well as the div curl = 0 and curl grad = 0 identities. For the diffusion equation, these operators describe compact schemes whose convergence is assured by the energy equations and which yield both the potential and the flux vector with second-order accuracy. A simplified potential form is especially useful for obtaining numerical results by multigrid and ADI methods.

Rose, Milton E.↗
NASA NTRS1991

A new finite element formulation for computational fluid dynamics. X - The compressible Euler and Navier-Stokes equations

A space-time element method is presented for solving the compressible Euler and Navier-Stokes equations. The proposed formulation includes the variational equation, predictor multi-corrector algorithms and boundary conditions. The variational equation is based on the time-discontinuous Galerkin method, in which the physical entropy variables are employed. A least-squares operator and a discontinuity-capturing operator are added, resulting in a high-order accurate and unconditionally stable method. Implicit/explicit predictor multi-corrector algorithms, applicable to steady as well as unsteady problems, are presented; techniques are developed to enhance their efficiency. Implementation of boundary conditions is addressed; in particular, a technique is introduced to satisfy nonlinear essential boundary conditions, and a consistent method is presented to calculate boundary fluxes. Numerical results are presented to demonstrate the performance of the method.

Shakib, Farzin↗
NASA NTRS1992

Subsonic High-Lift Flight Research on the NASA Transport System Research Vehicle (TSRV)

Flight tests are being conducted on the Transport Systems Research Vehicle (B737-100 aircraft) at the NASA Langley Research Center as part of a multi-phased research program to obtain detailed flow characteristics on a multi-element, high-lift flap system. Recent test results obtained for a full-chord wing section including the slat, main-wing, and flap elements are presented. In addition, program status and future plans are discussed. Pressure distributions were obtained using pressure belts over the slat, main-wing, and flap elements at one spanwise station. Test conditions included a range of flap deflections, chord Reynolds numbers (10 to 21 million), and Mach numbers (0.16 to 0.40). In addition to the pressure distributions, Preston-tube measurements on the slat upper surface indicated relaminarization of turbulent flows which were introduced by the pressure belt on the slat leading-edge surface. Comparisons of the in-flight pressure measurements with predictions from a two-dimensional, viscous multi-element method modified by simple-sweep theory showed reasonable agreement; however, the overprediction of pressures on the flap elements shows a need for better modeling of confluent boundary layers and inclusion of three-dimensional effects.

Flight Testing↗
NASA NTRS1993

In-flight pressure distributions and skin-friction measurements on a subsonic transport high-lift wing section

Flight experiments are being conducted as part of a multiphased subsonic transport high-lift research program for correlation with wind-tunnel and computational results. The NASA Langley Transport Systems Research Vehicle (B737-100 aircraft) is used to obtain in-flight flow characteristics at full-scale Reynolds numbers to contribute to the understanding of 3-D high-lift, multi-element flows including attachment-line transition and relaminarization, confluent boundary-layer development, and flow separation characteristics. Flight test results of pressure distributions and skin friction measurements were obtained for a full-chord wing section including the slat, main-wing, and triple-slotted, Fowler flap elements. Test conditions included a range of flap deflections, chord Reynolds numbers (10 to 21 million), and Mach numbers (0.16 to 0.40). Pressure distributions were obtained at 144 chordwise locations of a wing section (53-percent wing span) using thin pressure belts over the slat, main-wing, and flap elements. Flow characteristics observed in the chordwise pressure distributions included leading-edge regions of high subsonic flows, leading-edge attachment-line locations, slat and main-wing cove-flow separation and reattachment, and trailing-edge flap separation. In addition to the pressure distributions, limited skin-friction measurements were made using Preston-tube probes. Preston-tube measurements on the slat upper surface suggested relaminarization of the turbulent flow introduced by the pressure belt on the slat leading-edge surface when the slat attachment line was laminar. Computational analysis of the in-flight pressure measurements using two-dimensional, viscous multielement methods modified with simple-sweep theory showed reasonable agreement. However, overprediction of the pressures on the flap elements suggests a need for better detailed measurements and improved modeling of confluent boundary layers as well as inclusion of three-dimensional viscous effects in the analysis.

Yip, Long P.↗
NASA NTRS1993

Analysis of microstrip dipoles and slots transversely coupled to a microstrip line using the FDTD method

Printed dipole elements and their complement, linear slots, are elementary radiators that have found use in low-profile antenna arrays. Low-profile antenna arrays, in addition to their small size and low weight characteristics, offer the potential advantage of low-cost, high-volume production with easy integration with active integrated circuit components. The design of such arrays requires that the radiation and impedance characteristics of the radiating elements be known. The FDTD (Finite-Difference Time-Domain) method is a general, straight-forward implementation of Maxwell's equations and offers a relatively simple way of analyzing both printed dipole and slot elements. Investigated in this work is the application of the FDTD method to the analysis of printed dipole and slot elements transversely coupled to an infinite transmission line in a multilayered configuration. Such dipole and slot elements may be used in dipole and slot series-fed-type linear arrays, where element offsets and interelement line lengths are used to obtain the desired amplitude distribution and beam direction, respectively. The design of such arrays is achieved using transmission line theory with equivalent circuit models for the radiating elements. In an equivalent circuit model, the dipole represents a shunt impedance to the transmission line, where the impedance is a function of dipole offset, length, and width. Similarly, the slot represents a series impedance to the transmission line. The FDTD method is applied to single dipole and slot elements transversely coupled to an infinite microstrip line using a fixed rectangular grid with Mur's second order absorbing boundary conditions. Frequency-dependent circuit and scattering parameters are obtained by saving desired time-domain quantities and using the Fourier transform. A Gaussian pulse excitation is applied to the microstrip transmission line, where the resulting reflected signal due to the presence of the radiating element is used to determine the equivalent element impedance.

Tulintseff, A. N.↗
NASA NTRS1973

Elastic buckling analysis for composite stiffened panels and other structures subjected to biaxial inplane loads

An exact linear analysis method is presented for predicting buckling of structures with arbitrary uniform cross section. The structure is idealized as an assemblage of laminated plate-strip elements, curved and planar, and beam elements. Element edges normal to the longitudinal axes are assumed to be simply supported. Arbitrary boundary conditions may be specified on any external longitudinal edge of plate-strip elements. The structure or selected elements may be loaded in any desired combination of inplane transverse compression or tension side load and axial compression load. The analysis simultaneously considers all possible modes of instability and is applicable for the buckling of laminated composite structures. Numerical results correlate well with the results of previous analysis methods.

Viswanathan, A. V.↗
NASA NTRSDate not supplied

DeepONet-Assisted Optimization of Surface Topography for Transition Delay in A Mach 4.5 Boundary Layer

We use deep learning, an ensemble variationaltechnique (EnVar), and direct numerical simulations(DNS) to design an optimal topography for a two-dimensional roughness element that delays the on-set of laminar-turbulent transition in a Mach 4.5 flat-plate boundary layer. Deep operator networks (Deep-ONets), which have the known ability to learn com-plex nonlinear operators within dynamical systems,are used for machine learning. For the baseline config-uration of a smooth flat plate, the second-mode wavesat the DNS inflow cause a quick nonlinear breakdownof the high-speed boundary layer within the computa-tional domain. Results reported in the present studyvalidate the ability of DeepONets to model the tran-sition delay via a given topography of the roughnesselement. The computing cost to optimize the rough-ness element for minimal skin-friction drag is substan-tially lowered by the DeepONets-based reduced-ordermodel. In comparison to the baseline method of EnVaroptimization based on DNS alone, the DeepONets-based EnVar optimizer is able to delay transition pastthe outflow boundary of the computational domainwhile utilizing almost 5–6 times fewer DNS.

Machine Learning↗
NASA NTRS1994

High-lift flow-physics flight experiments on a subsonic civil transport aircraft (B737-100)

As part of the subsonic transport high-lift program, flight experiments are being conducted using NASA Langley's B737-100 to measure the flow characteristics of the multi-element high-lift system at full-scale high-Reynolds-number conditions. The instrumentation consists of hot-film anemometers to measure boundary-layer states, an infra-red camera to detect transition from laminar to turbulent flow, Preston tubes to measure wall shear stress, boundary-layer rakes to measure off-surface velocity profiles, and pressure orifices to measure surface pressure distributions. The initial phase of this research project was recently concluded with two flights on July 14. This phase consisted of a total of twenty flights over a period of about ten weeks. In the coming months the data obtained in this initial set of flight experiments will be analyzed and the results will be used to finalize the instrumentation layout for the next set of flight experiments scheduled for Winter and Spring of 1995. The main goal of these upcoming flights will be: (1) to measure more detailed surface pressure distributions across the wing for a range of flight conditions and flap settings; (2) to visualize the surface flows across the multi-element wing at high-lift conditions using fluorescent mini tufts; and (3) to measure in more detail the changes in boundary-layer state on the various flap elements as a result of changes in flight condition and flap deflection. These flight measured results are being correlated with experimental data measured in ground-based facilities as well as with computational data calculated with methods based on the Navier-Stokes equations or a reduced set of these equations. Also these results provide insight into the extent of laminar flow that exists on actual multi-element lifting surfaces at full-scale high-life conditions. Preliminary results indicate that depending on the deflection angle, the slat and flap elements have significant regions of laminar flow over a wide range of angles of attack. Boundary-layer transition mechanisms that were observed include attachment-line contamination on the slat and inflectional instability on the slat and fore flap. Also, the results agree fairly well with the predictions reported in a paper presented at last year's AIAA Fluid Dynamics Conference. The fact that extended regions of laminar flow are shown to exist on the various elements of the high-lift system raises the question what the effect is of loss of laminar flow as a result of insect contamiantion, rain or ice accumulation on high-life performance.

Vandam, Cornelis P.↗