Search NASASearch

SEARCH · Search NASA

Results for “Conformal mesh”

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 19 records

Balanced Central Schemes for the Shallow Water Equations on Unstructured Grids

We present a two-dimensional, well-balanced, central-upwind scheme for approximating solutions of the shallow water equations in the presence of a stationary bottom topography on triangular meshes. Our starting point is the recent central scheme of Kurganov and Petrova (KP) for approximating solutions of conservation laws on triangular meshes. In order to extend this scheme from systems of conservation laws to systems of balance laws one has to find an appropriate discretization of the source terms. We first show that for general triangulations there is no discretization of the source terms that corresponds to a well-balanced form of the KP scheme. We then derive a new variant of a central scheme that can be balanced on triangular meshes. We note in passing that it is straightforward to extend the KP scheme to general unstructured conformal meshes. This extension allows us to recover our previous well-balanced scheme on Cartesian grids. We conclude with several simulations, verifying the second-order accuracy of our scheme as well as its well-balanced properties.

Bryson, Steve

Finite area method for nonlinear conical flows

A fully conservative finite area method for the computation of steady inviscid flow about general conical bodies at incidence is described. The procedure utilizes the potential approximation and implements a body conforming mesh generator. The conical potential is assumed to have its best linear variation inside each mesh cell and a secondary interlocking cell system is used to establish the flux balance required to conserve mass. In the supersonic regions the scheme is desymmetrized by adding appropriate artificial viscosity in conservation form. The algorithm is nearly an order of a magnitude faster than present Euler methods and predicts known results accurately and qualitative features such as nodal point lift off correctly. Results are compared with those of other investigations.

Sritharan, S. S.

Finite area method for nonlinear supersonic conical flows

A fully conservative numerical method for the computation of steady inviscid supersonic flow about general conical bodies at incidence is described. The procedure utilizes the potential approximation and implements a body conforming mesh generator. The conical potential is assumed to have its best linear variation inside each mesh cell; a secondary interlocking cell system is used to establish the flux balance required to conserve mass. In the supersonic regions the scheme is symmetrized by adding artificial viscosity in conservation form. The algorithm is nearly an order of a magnitude faster than present Euler methods and predicts known results accurately and qualitative features such as nodal point lift off correctly. Results are compared with those of other investigators.

Sritharan, S. S.

Calculations Of Transonic Flow About A Wing

Report describes calculations of transonic airflows about wing in wind tunnel. Basic equations of flow used in study are Reynolds-averaged Navier-Stokes equations in strong conservation-law form. Equations of flow incorporated into finite-difference computer code called TNS (Transonic Navier-Stokes). Computational grid generated by solution of partial differential equations yielding smooth meshes conforming to surfaces of wing and wind tunnel.

Holst, Terry L.

Finite-difference algorithms for the time-domain Maxwell's equations - A numerical approach to RCS analysis

The applications of two CFD-based finite-difference methods to computational electromagnetics are investigated. In the first method, the time-domain Maxwell's equations are solved using the explicit Lax-Wendroff scheme and in the second method, the second-order wave equations satisfying the Maxwell's equations are solved using the implicit Crank-Nicolson scheme. The governing equations are transformed to a generalized curvilinear coordinate system and solved on a body-conforming mesh using the scattered-field formulation. The induced surface current and the bistatic radar cross section are computed and the results are validated for several two-dimensional test cases involving perfectly-conducting scatterers submerged in transverse-magnetic plane waves.

Vinh, Hoang

Mesh generation by conformal and quasiconformal mappings

It is pointed out that many recent advances in the finite-difference solution of elliptic equations have been limited to regions whose boundary contours coincide with coordinate lines of the Cartesian coordinate system. The reason for this is related to the fact that in the case of an arbitrary curvilinear coordinate system the original equation becomes much more complex. However, there is no added complexity if an orthogonal coordinate system is generated from a conformal mapping. In the present investigation, a finite difference method developed for the construction of conformal mappings has been generalized to construct quasi-conformal mappings. It is expected that the use of more sophisticated numerical algorithms could lead to improvements in both speed and accuracy. Quasi-conformal mappings have applications not only in the solution of elliptic equations but also in other areas such as orthogonal mesh generation on surfaces and the solution of certain fluid flow problems.

Mastin, C. W.

Finite-volume scheme for transonic potential flow about airfoils and bodies in an arbitrarily shaped channel

A conservative finite-volume difference scheme is developed for the potential equation to solve transonic flow about airfoils and bodies in an arbitrarily shaped channel. The scheme employs a mesh which is a nearly conformal O mesh about the airfoil and nearly orthogonal at the channel walls. The mesh extends to infinity upstream and downstream, where the mapping is singular. Special procedures are required to treat the singularities at infinity, including computation of the metrics near those points. Channels with exit areas different from inlet areas are solved; a body with a sting mount is an example of such a case.

South, Jerry C., Jr.

Finite-volume scheme for transonic potential flow about airfoils and bodies in an arbitrarily-shaped channel

A conservative finite-volume difference scheme is developed for the potential equation to solve transonic flow about airfoils and bodies in an arbitrary channel. The scheme employs a mesh which is a nearly-conformal 'O' mesh about the airfoil and nearly orthogonal at the channel walls. The mesh extends to infinity upstream and downstream, where the mapping is singular. Special procedures are required to treat the singularities at infinity, including computation of the metrics near those points. Channels with exit areas different from inlet areas are solved; a body with a sting mount is an example of such a case.

South, J. C., Jr.

LAVA Voronoi Mesher for Wall-Modeled Large-Eddy Simulations

The unstructured Voronoi mesher currently being developed within the Launch, Ascent, and Vehicle Aerodynamics (LAVA) software framework at NASA Ames Research Center is described in detail. The discussions include, but are not limited to, the strategies used in the seeding and smoothing phases to ensure a high-quality mesh for Wall-Modeled Large-Eddy Simulations (WMLES), the methods behind the cell clipping algorithm responsible for conforming the mesh to a complex geometry, and the approach to create a global mesh from the distinct Voronoi cells. Applications and benefits of the Voronoi meshing approach are also presented. The automated meshing paradigm introduced was able to significantly reduce the time necessary to create a high-quality mesh around a complex geometry when compared against the current LAVA curvilinear overset meshing standard. For example, the task of generating a single mesh around the high-lift common research model requires the full dedication of a mesh generation expert for a period between one and two months. A family of meshes with different refinement levels up to six hundred million cells can be generated by a single engineer in a day or two. This technology has the potential for decreasing the turnaround time for conducting WMLES around complex geometries, as well as facilitating comprehensive mesh refinement studies.

Victor C B Sousa

First benchmark of the Unstructured Grid Adaptation Working Group

Unstructured grid adaptation is a technology that holds the potential to improve the automation and accuracy of computational fluid dynamics and other computational disciplines. Difficulty producing the highly anisotropic elements necessary for simulation on complex curved geometries that satisfies a resolution request has limited this technology's widespread adoption. The Unstructured Grid Adaptation Working Group is an open gathering of researchers working on adapting simplicial meshes to conform to a metric field. Current members span a wide range of institutions including academia, industry, and national laboratories. The purpose of this group is to create a common basis for understanding and improving mesh adaptation. We present our first major contribution: a common set of benchmark cases, including input meshes and analytic metric specifications, that are publicly available to be used for evaluating any mesh adaptation code. We also present the results of several existing codes on these benchmark cases, to illustrate their utility in identifying key challenges common to all codes and important differences between available codes. Future directions are defined to expand this benchmark to mature the technology necessary to impact practical simulation workflows.

Ibanez, Daniel

A New Approach to Model Delamination Growth in Fatigue Using the Virtual Crack Closure Technique Without Re-Meshing

A modeling approach is proposed to simulate delamination propagation in fatigue that combines the Virtual Crack Closure Technique (VCCT) with a progressive nodal release strategy. The progressive nodal release alleviates the artificial stress concentrations found when using the VCCT with instantaneous release to model 3D delamination without re-meshing. This enables crack shapes that do not conform to the underlying mesh to be readily simulated. The progressive release is formulated such that no artificial growth or healing occurs as a result of a change in loading conditions. Furthermore, expressions for calculating energy release rates at partially released nodes are proposed and applied to arbitrarily shaped delamination fronts. The shape of the delamination front is considered explicitly in the energy release rate calculations, and is determined by considering the normalized crack positions associated with the partial nodal release. Contrary to previous implementations, the local delamination propagation direction is not assumed to follow mesh lines, but instead is computed as part of the iterative procedure proposed to determine the maximum energy release along the crack front. The results obtained suggest the approach can accurately simulate Mode I, Mixed Mode I/II and Mode II fatigue delamination growth and be used to simulate growth of arbitrarily shaped cracks that do not conform to the underlying mesh.

N V De Carvalho

Numerical calculation of steady inviscid full potential compressible flow about wind turbine blades

An exact nonlinear mathematical model that accounts for three-dimensional cascade effects about the inner portions of the rotor blades and compressibility effects about the tip regions of the blades was derived. An artificially time dependent version was iteratively solved by a finite volume technique involving an artificial viscosity and a three-level consecutive mesh refinement. The exact boundary conditions were applied by generating a boundary conforming periodic computation mesh.

Dulikravich, D. S.

Numerical calculation of steady inviscid full potential compressible flow about wind turbine blades

The air flow through a propeller-type wind turbine rotor is characterized by three-dimensional rotating cascade effects about the inner portions of the rotor blades and compressibility effects about the tip regions of the blades. In the case of large rotor diameter and/or increased rotor angular speed, the existence of small supersonic zones terminated by weak shocks is possible. An exact nonlinear mathematical model (called a steady Full Potential Equation - FPE) that accounts for the above phenomena has been rederived. An artificially time dependent version of FPE was iteratively solved by a finite volume technique involving an artificial viscosity and a three-level consecutive mesh refinement. The exact boundary conditions were applied by generating a boundary conforming periodic computation mesh.

Dulikravich, D. S.

A new conformal absorbing boundary condition for finite element meshes and parallelization of FEMATS

Some of the progress toward the development and parallelization of an improved version of the finite element code FEMATS is described. This is a finite element code for computing the scattering by arbitrarily shaped three dimensional surfaces composite scatterers. The following tasks were worked on during the report period: (1) new absorbing boundary conditions (ABC's) for truncating the finite element mesh; (2) mixed mesh termination schemes; (3) hierarchical elements and multigridding; (4) parallelization; and (5) various modeling enhancements (antenna feeds, anisotropy, and higher order GIBC).

Chatterjee, A.

A finite element solution of unsteady transonic flow problems for three-dimensional wings and bodies

A computerized finite element method has been developed to solve a form of the three-dimensional inviscid, irrotational unsteady transonic equation for pointed and rectangular wings and wing-body combinations. The equation is a 'parabolized' version of the full unsteady transonic equation which includes the effect of the advancing wave and the acceleration in the steady flow field thus including a thickness effect. Automatic mesh generation using conformal mapping has the advantage of allowing relatively arbitrary cross-sections to be closely approximated while taking advantage of the body's similarity to a more simple configuration. Plunge, pitch, and roll modes are computed

Kimble, K. R.

Aeroelastic analysis of wings using the Euler equations with a deforming mesh

Modifications to the CFL3D three-dimensional unsteady Euler/Navier-Stokes code for the aeroelastic analysis of wings are described. The modifications involve including a deforming mesh capability which can move the mesh to continuously conform to the instantaneous shape of the aeroelastically deforming wing, and including the structural equations of motion for their simultaneous time-integration with the governing flow equations. Calculations were performed using the Euler equations to verify the modifications to the code and as a first-step toward aeroelastic analysis using the Navier-Stokes equations. Results are presented for the NACA 0012 airfoil and a 45 deg sweptback wing to demonstrate applications of CFL3D for generalized force computations and aeroelastic analysis. Comparisons are made with published Euler results for the NACA 0012 airfoil and with experimental flutter data for the 45 deg sweptback wing to assess the accuracy of the present capability. These comparisons show good agreement and, thus, the CFL3D code may be used with confidence for aeroelastic analysis of wings. The paper describes the modifications that were made to the code and presents results and comparisons which assess the capability.

Robinson, Brian A.

Aeroelastic analysis of wings using the Euler equations with a deforming mesh

Modifications to the CFL3D three dimensional unsteady Euler/Navier-Stokes code for the aeroelastic analysis of wings are described. The modifications involve including a deforming mesh capability which can move the mesh to continuously conform to the instantaneous shape of the aeroelastically deforming wing, and including the structural equations of motion for their simultaneous time-integration with the governing flow equations. Calculations were performed using the Euler equations to verify the modifications to the code and as a first step toward aeroelastic analysis using the Navier-Stokes equations. Results are presented for the NACA 0012 airfoil and a 45 deg sweptback wing to demonstrate applications of CFL3D for generalized force computations and aeroelastic analysis. Comparisons are made with published Euler results for the NACA 0012 airfoil and with experimental flutter data for the 45 deg sweptback wing to assess the accuracy of the present capability. These comparisons show good agreement and, thus, the CFL3D code may be used with confidence for aeroelastic analysis of wings.

Robinson, Brian A.

Combining Progressive Nodal Release with the Virtual Crack Closure Technique to Model Fatigue Delamination Growth Without Re-Meshing

The present work summarizes an approach to model mixed-mode 3D fatigue crack growth using the Virtual Crack Closure Technique (VCCT) without requiring re-meshing. It is demonstrated that the proposed approach can be used to simulate crack shapes that do not conform to the underlying mesh. The proposed approach relies solely on Paris Law characterization data to model delamination growth. Growth is determined as a post-processing step at the end of each increment, and hence no convergence issues associated with the progressive nodal release are encountered. This approach can be readily applied using standard solid element formulations and is implemented via an interface user element in Abaqus/Standard.

Carvalho, N. V. De