SEARCH · Engineering Papers

Results for “finite difference”

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 1,009 records · Page 56

All sources NASA NTRS DOE OSTI

NASA NTRS1991

On the accuracy of shape sensitivity

The calculation of sensitivity of the response of a structure modeled by finite elements to shape variation is known to be subject to numerical difficulties. The accuracy of a given method is typically measured against the yard stick of finite-difference sensitivity calculation. The present paper demonstrates with a simple example that this approach may be flawed because of discretization errors associated with the finite element mesh. Seven methods for calculating sensitivity derivatives are compared for a two-material beam problem with a moving interface. It is found that as the mesh is refined, displacement sensitivity derivatives converge more slowly than the displacements. Six of the methods agree fairly well, but the adjoint variational surface method provides substantially different results. However, the difference is found to reflect convergence from another direction to the same answer rather than reduced accuracy. Additionally, it is observed that small derivatives are particularly prone to accuracy problems.

Haftka, R. T.↗

NASA NTRS1975

Pressure- and buoyancy-driven thermal convection in a rectangular enclosure

Results are presented for unsteady laminar thermal convection in compressible fluids at various reduced levels of gravity in a rectangular enclosure which is heated on one side and cooled on the opposite side. The results were obtained by solving numerically the equations of conservation for a viscous, compressible, heat-conducting, ideal gas in the presence of a gravitational body force. The formulation differs from the Boussinesq simplification in that the effects of variable density are completely retained. A conservative, explicit, time-dependent, finite-difference technique was used and good agreement was found for the limited cases where direct comparison with previous investigations was possible. The solutions show that the thermally induced motion is acoustic in nature at low levels of gravity and that the unsteady-state rate of heat transfer is thereby greatly enhanced relative to pure conduction. The nonlinear variable density profile skews the streamlines towards the cooler walls but is shown to have little effect on the steady-state isotherms.

Spradley, L. W.↗

NASA NTRS1989

Aerodynamic interaction between vortical wakes and lifting two-dimensional bodies

Unsteady rotor wake interactions with the empennage, tail boom, and other aerodynamic surfaces of a helicopter have a significant influence on its aerodynamic performance, the ride quality, and vibration. A numerical method for computing the aerodynamic interaction between an interacting vortex wake and the viscous flow about arbitrary two-dimensional bodies was developed to address this helicopter problem. The method solves for the flow field velocities on a body-fitted computational mesh using finite-difference techniques. The interacting vortex wake is represented by an array of discrete vortices which, in turn, are represented by a finite-core model. The evolution of the interacting vortex wake is calculated by Lagrangian techniques. The viscous flow field of the two-dimensional body is calculated on an Eulerian grid. The flow around circular and elliptic cylinders in the absence of an interacting vortex wake was calculated. These results compare very well with other numerical results and with results obtained from experiment and thereby demonstrate the accuracy of the viscous solution. The interaction of a rotor wake with the flow about a 4 to 1 elliptic cylinder at 45 degree incidence was calculated for a Reynolds number of 3000. The results demonstrate the significant variations in the lift and drag on the elliptic cylinder in the presence of the interacting rotor wake.

Stremel, Paul M.↗

NASA NTRS1992

The semi-discrete Galerkin finite element modelling of compressible viscous flow past an airfoil

A method is developed to solve the two-dimensional, steady, compressible, turbulent boundary-layer equations and is coupled to an existing Euler solver for attached transonic airfoil analysis problems. The boundary-layer formulation utilizes the semi-discrete Galerkin (SDG) method to model the spatial variable normal to the surface with linear finite elements and the time-like variable with finite differences. A Dorodnitsyn transformed system of equations is used to bound the infinite spatial domain thereby permitting the use of a uniform finite element grid which provides high resolution near the wall and automatically follows boundary-layer growth. The second-order accurate Crank-Nicholson scheme is applied along with a linearization method to take advantage of the parabolic nature of the boundary-layer equations and generate a non-iterative marching routine. The SDG code can be applied to any smoothly-connected airfoil shape without modification and can be coupled to any inviscid flow solver. In this analysis, a direct viscous-inviscid interaction is accomplished between the Euler and boundary-layer codes, through the application of a transpiration velocity boundary condition. Results are presented for compressible turbulent flow past NACA 0012 and RAE 2822 airfoils at various freestream Mach numbers, Reynolds numbers, and angles of attack. All results show good agreement with experiment, and the coupled code proved to be a computationally-efficient and accurate airfoil analysis tool.

Meade, Andrew J., Jr.↗

NASA NTRS1996

A Magnetically Suspended Wheel for a Miniature Gyro Made Using Planar Fabrication Technologies

The technical feasibility of a magnetically suspended rotating wheel for miniature gyro applications was investigated under a NASA SBIR contract. A concept for a configuration for a system of compact, lightweight magnetic actuators capable of generating the necessary suspension forces and fabrication using millimachining planar fabrication technologies was developed. Both capacitive and electromagnetic position sensing concepts were developed for implementing a closed loop control system for supporting the wheel. A finite difference technique, implemented in a spreadsheet environment, for analyzing the force characteristics of the actuator was used and the results verified with Finite Element Analysis.

Dauwalter, Charles R.↗

NASA NTRS1996

Algebraic Multigrid by Smoothed Aggregation for Second and Fourth Order Elliptic Problems

Multigrid methods are very efficient iterative solvers for system of algebraic equations arising from finite element and finite difference discretization of elliptic boundary value problems. The main principle of multigrid methods is to complement the local exchange of information in point-wise iterative methods by a global one utilizing several related systems, called coarse levels, with a smaller number of variables. The coarse levels are often obtained as a hierarchy of discretizations with different characteristic meshsizes, but this requires that the discretization is controlled by the iterative method. To solve linear systems produced by existing finite element software, one needs to create an artificial hierarchy of coarse problems. The principal issue is then to obtain computational complexity and approximation properties similar to those for nested meshes, using only information in the matrix of the system and as little extra information as possible. Such algebraic multigrid method that uses the system matrix only was developed by Ruge. The prolongations were based on the matrix of the system by partial solution from given values at selected coarse points. The coarse grid points were selected so that each point would be interpolated to via so-called strong connections. Our approach is based on smoothed aggregation introduced recently by Vanek. First the set of nodes is decomposed into small mutually disjoint subsets. A tentative piecewise constant interpolation (in the discrete sense) is then defined on those subsets as piecewise constant for second order problems, and piecewise linear for fourth order problems. The prolongation operator is then obtained by smoothing the output of the tentative prolongation and coarse level operators are defined variationally.

Vanek, Petr↗

NASA NTRS2011

Ablator Modeling: Why Not Much Has Changed Over the Past 45 plus Years

In the late 1950s, the earliest models describing the thermal response of ablative materials were based on the heat of ablation concept, which is an empirical approach that was reasonable for the types of materials of interest at that time. In the early-mid 60s the models were expanded to include pyrolysis since organic resin composites became the TPS materials of interest. However, surface recession was still predominantly modeled via empirical correlation. The development of the 1-D CMA finite difference code in the mid-late 60s introduced the thermochemical ablation approach for gas/surface interactions. Since that time investigators have developed finite volume and finite element codes, in 1-D, 2-D and 3-D, but the basic modeling has not evolved significantly. Models describing internal gas pressure due to pyrolysis, particle impact erosion, in-depth radiant transport, etc., have been added to address specific problems, but the fundamental modeling has not evolved. The reasons for this stagnation, as viewed by the author, will be described.

Laub, B.↗

NASA NTRS2011

Ascent Heating Thermal Analysis on Spacecraft Adaptor Fairings

When the Crew Exploration Vehicle (CEV) is launched, the spacecraft adaptor (SA) fairings that cover the CEV service module (SM) are exposed to aero heating. Thermal analysis is performed to compute the fairing temperatures and to investigate whether the temperatures are within the material limits for nominal ascent aeroheating case. The ascent heating is analyzed by using computational fluid dynamics (CFD) and engineering codes at Marshall Space Flight Center. The aeroheating environment data used for this work is known as Thermal Environment 3 (TE3) heating data. One of the major concerns is with the SA fairings covering the CEV SM and the SM/crew launch vehicle (CLV) flange interface. The TE3 heating rate is a function of time, wall temperature, and the spatial locations. The implementation of the TE3 heating rate as boundary conditions in the thermal analysis becomes challenging. The ascent heating thermal analysis on SA fairings and SM/CLV flange interface are performed using two commercial software packages: Cullimore & Ring (C&R) Thermal Desktop (TD) 5.1 and MSC Patran 2007r1 b. TD is the pre-and post-processor for SINDA, which is a finite-difference-based solver. In TD, the geometry is built and meshed, the boundary conditions are defined, and then SINDA is used to compute temperatures. MSC Pthermal is a finite-element- based thermal solver. MSC Patran is the pre- and post-processor for Pthermal. Regarding the boundary conditions, the convection, contact resistance, and heat load can be imposed in different ways in both programs. These two software packages are used to build the thermal model for the same analysis to validate each other and show the differences in the modeling details.

Wang, Xiao Yen↗

NASA NTRS1976

Exploiting symmetries for efficient post-buckling analysis of composite plates

Two aspects of the post-buckling analysis of composite plates are considered in this paper. The first pertains to identifying the different types of symmetry exhibited by the pre- and post-buckling responses of rectangular and skewed composite plates. A procedure is presented for exploiting these symmetries in finite-element and finite-difference analyses. The procedure can be used with existing programs having multipoint constraint capability. The second aspect pertains to the post-buckling response of biaxially loaded composite plates. Numerical results from finite-element analysis are presented which show the range of parameters for which shear deformation and/or anisotropy are important in buckling and post-buckling analysis.

Noor, A. K.↗

NASA NTRS1988

Numerical simulation of axisymmetric turbulent flow in combustors and diffusors

A method for predicting turbulent flow in combustors and diffusers is developed. The Navier-Stokes equations, incorporating a turbulence kappa-epsilon model equation, were solved in a nonorthogonal curvilinear coordinate system. The solution applied the finite volume method to discretize the differential equations and utilized the SIMPLE algorithm iteratively to solve the differenced equations. A zonal grid method, wherein the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. However, treatment of the zonal boundaries required special handling. Boundary overlap and interpolating techniques were used and an adjustment of the flow variables was required to assure conservation of mass, momentum and energy fluxes. The numerical accuracy was assessed using different finite differencing methods, i.e., hybrid, quadratic upwind and skew upwind, to represent the convection terms. Flows in different geometries of combustors and diffusers were simulated and results compared with experimental data and good agreement was obtained.

Yung, Chain Nan↗

DOE OSTI2022

Angle- and polarization-resolved luminescence from suspended and hexagonal boron nitride encapsulated MoSe 2 monolayers

The polarized photoluminescence from atomically thin transition metal dichalcogenides is a frequently applied tool to scrutinize optical selection rules and valley physics, yet it is known to sensibly depend on a variety of internal and external material and sample properties. In this work, we apply combined angle- and polarization-resolved spectroscopy to explore the interplay of excitonic physics and phenomena arising from the commonly utilized encapsulation procedure on the optical properties of atomically thin M o S e 2 . We probe monolayers prepared in both suspended and encapsulated manners. We show that the hBN encapsulation significantly enhances the linear polarization of exciton photoluminescence emission at large emission angles. This degree of linear polarization of excitons can increase up to ∼ <#comment/> 17 % <#comment/> in the hBN encapsulated samples. As we confirm by finite-difference time-domain simulations, it can be directly connected to the optical anisotropy of the hBN layers. In comparison, the linear polarization at finite exciton momenta is significantly reduced in a suspended M o S e 2 monolayer, and becomes notable only in cryogenic conditions. This phenomenon strongly suggests that the effect is rooted in the k-dependent anisotropic exchange coupling inherent in 2D excitons. Our results have strong implications on further studies on valley contrasting selection rules and valley coherence phenomena using standard suspended and encapsulated samples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1986

Least squares finite element simulation of transonic flows

Finite difference approximation of transonic flow problems is a well-developed and largely successful approach. Nevertheless, there is still a real need to develop finite element methods for applications arising from fluid-structure interactions and problems with complicated boundaries. In this paper a least squares based finite element scheme is introduced. It is shown that, if suitably formulated, such an approach can lead to physically meaningful results. Bottlenecks that arise from such schemes are also discussed.

Chen, T. F.↗

DOE OSTI2022

Droplet formation simulation using mixed finite elements

Droplet formation happens in finite time due to the surface tension force. The linear stability analysis is useful to estimate the size of a droplet but fails to approximate the shape of the droplet. This is due to a highly nonlinear flow description near the point where the first pinch-off happens. A one-dimensional axisymmetric mathematical model was first developed by Eggers and Dupont [“Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] using asymptotic analysis. This asymptotic approach to the Navier–Stokes equations leads to a universal scaling explaining the self-similar nature of the solution. Numerical models for the one-dimensional model were developed using the finite difference [Eggers and Dupont, “Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] and finite element method [Ambravaneswaran et al., “Drop formation from a capillary tube: Comparison of one-dimensional and two-dimensional analyses and occurrence of satellite drops,” Phys. Fluids 14, 2606–2621 (2002)]. The focus of this study is to provide a robust computational model for one-dimensional axisymmetric droplet formation using the Portable, Extensible Toolkit for Scientific Computation. The code is verified using the Method of Manufactured Solutions and validated using previous experimental studies done by Zhang and Basaran [“An experimental study of dynamics of drop formation,” Phys. Fluids 7, 1184–1203 (1995)]. The present model is used for simulating pendant drops of water, glycerol, and paraffin wax, with an aspiration of extending the application to simulate more complex pinch-off phenomena.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1976

A three-dimensional, time-dependent model of Mobile Bay

A three-dimensional, time-variant mathematical model for momentum and mass transport in estuaries was developed and its solution implemented on a digital computer. The mathematical model is based on state and conservation equations applied to turbulent flow of a two-component, incompressible fluid having a free surface. Thus, bouyancy effects caused by density differences between the fresh and salt water, inertia from thare river and tidal currents, and differences in hydrostatic head are taken into account. The conservation equations, which are partial differential equations, are solved numerically by an explicit, one-step finite difference scheme and the solutions displayed numerically and graphically. To test the validity of the model, a specific estuary for which scaled model and experimental field data are available, Mobile Bay, was simulated. Comparisons of velocity, salinity and water level data show that the model is valid and a viable means of simulating the hydrodynamics and mass transport in non-idealized estuaries.

Pitts, F. H.↗

NASA NTRS1979

Documentation of the Fourth Order Band Model

A general circulation model is presented which uses quadratically conservative, fourth order horizontal space differences on an unstaggered grid and second order vertical space differences with a forward-backward or a smooth leap frog time scheme to solve the primitive equations of motion. The dynamic equations for motion, finite difference equations, a discussion of the structure and flow chart of the program code, a program listing, and three relevent papers are given.

Kalnay-Rivas, E.↗

NASA NTRS1980

A new approach in cascade flow analysis using the finite element method

A new approach in analyzing the potential flow past cascades and single airfoils using the finite element method is developed. In this analysis the circulation around the airfoil is not externally imposed but is directly computed in the numerical solution. Different finite element discretization patterns, orders of piecewise approximation, and grid sizes are used in the solution. The results obtained are compared with existing experimental measurements and exact solutions in cascades and single airfoils.

Baskharone, E.↗

NASA NTRS1999

Interaction of Porosity with an Advancing Solid/Liquid Interface: a Real-Time Investigation

Problems associated with formation of porosity during solidification continue to have a daily impact on the metal forming industry. Several past investigations have dealt with the nucleation and growth aspects of porosity. However, investigations related to the interaction of porosity with that of a solidification front has been limited mostly to organic analogues. In this paper we report on real time experimental observations of such interactions in metal alloys. Using a state of the art X-Ray Transmission Microscope (XTM) we have been able to observe and record the dynamics of the interaction. This includes distortion of the solid/liquid interface near a poro.sity, solute segr,egation patterns surrounding a porosity and the change in shape of the porosity during interaction with an advancing solid/liquid interface. Results will be presented for different Al alloys and growth conditions. The experimental data will be compared to theory using a recently developed 2D numerical model. The model employs a finite difference approach where the solid/liquid interface is defined through the points at which the interface intersects the grid lines. The transport variables are calculated at these points and the motion of the solidification front is determined by the magnitude of the transport variables. The model accounts for the interplay of the thermal and solutal field and the influence of capilarity to predict the shape of the solid/liquid interface with time in the vicinity of porosity. One can further calculate the perturbation of the solutal field by the presence of porosity in the melt.

Sen, S.↗

NASA NTRS2013

The Effect of Scale Dependent Discretization on the Progressive Failure of Composite Materials Using Multiscale Analyses

A multiscale modeling methodology, which incorporates a statistical distribution of fiber strengths into coupled micromechanics/ finite element analyses, is applied to unidirectional polymer matrix composites (PMCs) to analyze the effect of mesh discretization both at the micro- and macroscales on the predicted ultimate tensile (UTS) strength and failure behavior. The NASA code FEAMAC and the ABAQUS finite element solver were used to analyze the progressive failure of a PMC tensile specimen that initiates at the repeating unit cell (RUC) level. Three different finite element mesh densities were employed and each coupled with an appropriate RUC. Multiple simulations were performed in order to assess the effect of a statistical distribution of fiber strengths on the bulk composite failure and predicted strength. The coupled effects of both the micro- and macroscale discretizations were found to have a noticeable effect on the predicted UTS and computational efficiency of the simulations.

Ricks, Trenton M.↗