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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The Flux-Differencing Discontinuous Galerkin Method Applied to an Idealized Fully Compressible Nonhydrostatic Dry Atmosphere

Dynamical cores used to study the circulation of the atmosphere employ various numerical methods ranging from finite-volume, spectral element, global spectral, and hybrid methods. In this work, we explore the use of Flux-Differencing Discontinuous Galerkin (FDDG) methods to simulate a fully compressible dry atmosphere at various resolutions. We show that the method offers a judicious compromise between high-order accuracy and stability for large-eddy simulations and simulations of the atmospheric general circulation. In particular, filters, divergence damping, diffusion, hyperdiffusion, or sponge-layers are not required to ensure stability; only the numerical dissipation naturally afforded by FDDG is necessary. We apply the method to the simulation of dry convection in an atmospheric boundary layer and in a global atmospheric dynamical core in the standard benchmark of Held and Suarez

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1990

The effect of coatings and liners on heat transfer in a dry shaft-bush tribosystem

The temperatures due to frictional heating within a solid lubricated or coated journal bearing were analyzed by using a finite element method. A solid model of the shaft-bush tribocontact was generated with an eight-node, three-dimensional, first-order isoparametric heat-transfer element and the Patran solid modeler software. The Patmar (Patran-Marc) translator was used to help develop the Marc-based finite element program for the system; this software was used on the Cray X-MP supercomputer to perform a finite element analysis of the contact. The analysis was performed for various liner materials, for thin, hard, wear-resistant coated bearings, and for different geometries and thermal cooling boundary conditions. The analyses indicated that thermal conductivity of the liner or coating material is the most vital thermal parameter that controls the interface temperature. In addition to design variations, the proximity of the cooling source to the heat-flux-generating interface is critically important to the temperature control in the system.

Ghosh, Mihir K.↗
NASA NTRS2001

Investigation of Conjugate Heat Transfer in Turbine Blades and Vanes

We report on work carried out to develop a 3-D coupled Finite Volume/BEM-based temperature forward/flux back (TFFB) coupling algorithm to solve the conjugate heat transfer (CHT) which arises naturally in analysis of systems exposed to a convective environment. Here, heat conduction within a structure is coupled to heat transfer to the external fluid which is convecting heat into or out of the solid structure. There are two basic approaches to solving coupled fluid structural systems. The first is a direct coupling where the solution of the different fields is solved simultaneously in one large set of equations. The second approach is a loose coupling strategy where each set of field equations is solved to provide boundary conditions for the other. The equations are solved in turn until an iterated convergence criterion is met at the fluid-solid interface. The loose coupling strategy is particularly attractive when coupling auxiliary field equations to computational fluid dynamics codes. We adopt the latter method in which the BEM is used to solve heat conduction inside a structure which is exposed to a convective field which in turn is resolved by solving the NASA Glenn compressible Navier-Stokes finite volume code Glenn-HT. The BEM code features constant and bi-linear discontinuous elements and an ILU-preconditioned GMRES iterative solver for the resulting non-symmetric algebraic set arising in the conduction solution. Interface of flux and temperature is enforced at the solid/fluid interface, and a radial-basis function scheme is used to interpolated information between the CFD and BEM surface grids. Additionally, relaxation is implemented in passing the fluxes from the conduction solution to the fluid solution. Results from a simple test example are reported.

Kassab, A. J.↗
DOE OSTI2024

An Efficient Numerical Algorithm for Solving Coupled Time-Dependent Ginzburg-Landau Equation for Superconductivity and Elasticity

A decoupled finite element algorithm is developed for simulating the vortex dynamics on an elastic superconductor which couples the time-dependent Ginzburg- Landau equation with the complex-valued superconducting order parameter and the vector-valued magnetic potential, and the elasticity equation. We present an iterative algorithm for the decoupled system arising from the time and spatial discretization using a combination of preconditioner, algebraic multigrid method (AMG) and preconditioned conjugate gradient method (PCG). The iterative algorithm allows us to perform large-scale three-dimensional simulations of mesoscale pattern formation during superconducting phase transitions with arbitrary elastic boundary conditions. Here, the performance and efficiency of the algorithm are numerically verified by several benchmark problems, exhibiting up to two orders of magnitude improvement depending on the scale of discrete system compared to the exact solver.

Efficiency↗
NASA NTRS1991

Time-domain finite elements in optimal control with application to launch-vehicle guidance

A time-domain finite element method is developed for optimal control problems. The theory derived is general enough to handle a large class of problems including optimal control problems that are continuous in the states and controls, problems with discontinuities in the states and/or system equations, problems with control inequality constraints, problems with state inequality constraints, or problems involving any combination of the above. The theory is developed in such a way that no numerical quadrature is necessary regardless of the degree of nonlinearity in the equations. Also, the same shape functions may be employed for every problem because all strong boundary conditions are transformed into natural or weak boundary conditions. In addition, the resulting nonlinear algebraic equations are very sparse. Use of sparse matrix solvers allows for the rapid and accurate solution of very difficult optimization problems. The formulation is applied to launch-vehicle trajectory optimization problems, and results show that real-time optimal guidance is realizable with this method. Finally, a general problem solving environment is created for solving a large class of optimal control problems. The algorithm uses both FORTRAN and a symbolic computation program to solve problems with a minimum of user interaction. The use of symbolic computation eliminates the need for user-written subroutines which greatly reduces the setup time for solving problems.

Bless, Robert R.↗
NASA NTRS1991

Heat transfer in the turbulent boundary layer with a step change in surface roughness

Measurements of Stanton numbers, velocity profiles, temperature profiles, and turbulence intensity profiles are reported for turbulent flat plate boundary layer flows with a step change in surface roughness. The first 0.9 m length of the test surface is roughened with 1.27 mm diameter hemispheres spaced 2 base diameters apart in a staggered array. The remaining 1.5 m length is smooth. The experiments show that the step change from a rough to a smooth surface has a dramatic effect on the convective heat transfer. In many cases, the Stanton number drops below the smooth-wall correlation immediately downstream of the change in roughness. The Stanton number measurements are compared with predictions using the discrete element method with excellent results.

Taylor, Robert P.↗
NASA NTRS2013

Recent Development in the CESE Method for the Solution of the Navier-Stokes Equations Using Unstructured Triangular or Tetrahedral Meshes With High Aspect Ratio

In the multidimensional CESE development, triangles and tetrahedra turn out to be the most natural building blocks for 2D and 3D spatial meshes. As such the CESE method is compatible with the simplest unstructured meshes and thus can be easily applied to solve problems with complex geometries. However, because the method uses space-time staggered stencils, solution decoupling may become a real nuisance in applications involving unstructured meshes. In this paper we will describe a simple and general remedy which, according to numerical experiments, has removed any possibility of solution decoupling. Moreover, in a real-world viscous flow simulation near a solid wall, one often encounters a case where a boundary with high curvature or sharp corner is surrounded by triangular/tetrahedral meshes of extremely high aspect ratio (up to 106). For such an extreme case, the spatial projection of a space-time compounded conservation element constructed using the original CESE design may become highly concave and thus its centroid (referred to as a spatial solution point) may lie far outside of the spatial projection. It could even be embedded beyond a solid wall boundary and causes serious numerical difficulties. In this paper we will also present a new procedure for constructing conservation elements and solution elements which effectively overcomes the difficulties associated with the original design. Another difficulty issue which was addressed more recently is the wellknown fact that accuracy of gradient computations involving triangular/tetrahedral grids deteriorates rapidly as the aspect ratio of grid cells increases. The root cause of this difficulty was clearly identified and several remedies to overcome it were found through a rigorous mathematical analysis. However, because of the length of the current paper and the complexity of mathematics involved, this new work will be presented in another paper.

Chang, Sin-Chung↗
NASA NTRS1979

Multi-level adaptive finite element methods. 1: Variation problems

A general numerical strategy for solving partial differential equations and other functional problems by cycling between coarser and finer levels of discretization is described. Optimal discretization schemes are provided together with very fast general solvers. It is described in terms of finite element discretizations of general nonlinear minimization problems. The basic processes (relaxation sweeps, fine-grid-to-coarse-grid transfers of residuals, coarse-to-fine interpolations of corrections) are directly and naturally determined by the objective functional and the sequence of approximation spaces. The natural processes, however, are not always optimal. Concrete examples are given and some new techniques are reviewed. Including the local truncation extrapolation and a multilevel procedure for inexpensively solving chains of many boundary value problems, such as those arising in the solution of time-dependent problems.

Brandt, A.↗
NASA NTRS1972

Two-dimensional finite-element temperature variance analysis

The finite element method is extended to thermal analysis by forming a variance analysis of temperature results so that the sensitivity of predicted temperatures to uncertainties in input variables is determined. The temperature fields within a finite number of elements are described in terms of the temperatures of vertices and the variational principle is used to minimize the integral equation describing thermal potential energy. A computer calculation yields the desired solution matrix of predicted temperatures and provides information about initial thermal parameters and their associated errors. Sample calculations show that all predicted temperatures are most effected by temperature values along fixed boundaries; more accurate specifications of these temperatures reduce errors in thermal calculations.

Heuser, J. S.↗
DOE OSTI2024

Comparison of Gas Phase Fragmentation Behaviors of Nuclear Fuel Cycle Ligands in Lanthanide and Americium Metal Ligand Nitrate Clusters

Introduction (120 words): Transport of metal ions across the aqueous-organic phase boundary is an essential step in a hydrometallurgical nuclear fuel reprocessing strategy. The study of transport agents for nuclear fuel elements is imperative to guide the design of ligands that boost the separation efficiency of the recovery process from fission products. However, limited studies have been made on the chemistry of these transport agents when complexing with transuranic elements in gas-phase where all surrounding factors are essentially excluded. This work investigates the reagent ligand complexations to transuranic and other metals and their dissociations in the gas phase. Comparisons are made between 4f and 5f elements and between ligands. Methods (120 words): (N,N-diisobutylcarbamoylmethyl)phenyloctylphosphine oxide (CMPO) and N,N,N',N'-tetraoctyldiglycolamide (TODGA) have been selected to complex with metal nitrates. The actinide americium and lanthanides neodymium, samarium, and europium were investigated as part of this work. The lanthanides were selected to act as size and electron configuration analogues of the minor actinides. Metal complexes with two ligands and two nitrates ([M(NO3)2(CMPO)2]+, for example) are studied in Bruker micrOTOF-Q II mass spectrometer equipped with collision-induced dissociation capability. The comparisons of the mass spectra are made in groups of homogenous ligands and mixed TODGA-CMPO ligands clusters. Comparisons are also made based on the complexed metals (Am and lanthanides). Preliminary data (300 words): Collision-induced dissociation mass spectrometry data are collected on two ligands complexed with metal nitrates where the two ligands are homogenous, with (CMPO)2 or (TODGA)2, or heterogeneous, with (TODGA)(CMPO). Several fragmentation patterns are observed among complexes with the CMPO ligand whereas the TODGA ligand commonly dissociates intact from the complex. Most of the metal complexes exhibit similar fragmentation patterns, but there are a few notable deviations in fragmentation patterns between the Am and Ln-bearing complexes. For the [M(CMPO)2(NO3)2]+ complexes, the initial loss of nitrate in the form of nitric acid is observed in all four complexes. However, [Am(CMPO)2(NO3)2]+ exhibits an additional fragmentation not found in the lanthanide complexes. Also, a significantly different ratio of the second nitric acid loss is found in the Am complex. These deviations may indicate the different interaction behaviors between actinides and lanthanides. The [M(TODGA)2(NO3)2]+ complexes exhibit the fragmentation as the loss of one TODGA ligand as an intact form and the loss of nitrate as nitric acid. The Am complex exhibits an additional fragmentation after losing the TOGDA ligand, which is not observed among the Ln complexes. The heterogeneous [M(TODGA)(CMPO)(NO3)2]+ complexes exhibit both similarities and differences between the Am and Ln complexes. For example, the heterogenous Am complex does not exhibit the loss of an intact TODGA ligand while all three Ln complexes do. This indicates that TODGA may bind more strongly to Am than Ln. Additionally, the intensity of the loss of CMPO ligand (as partially or whole) is found to be significantly larger than that of the loss of TODGA (as partially of whole) indicating that TODGA is bound to the metal significantly stronger than CMPO. Planned computational analysis will help understand the deviation in fragmentation behaviors between americium and lanthanide metal centers, or between TODGA and CMPO ligands. Novel aspect (20 words): Gas-phase actinide and lanthanide complex formation and fragmentation provide insight into the coordination environment differences of f-element metals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Comparison of Gas Phase Fragmentation Behaviors of Nuclear Fuel Cycle Ligands in Lanthanide and Americium Metal Ligand Nitrate Clusters

Title (20 word): Comparison of Gas Phase Fragmentation Behaviors of Nuclear Fuel Cycle Ligands in Lanthanide and Americium Metal Ligand Nitrate Clusters Introduction (120 words): Transport of metal ions across the aqueous-organic phase boundary is an essential step in a hydrometallurgical nuclear fuel reprocessing strategy. The study of transport agents for nuclear fuel elements is imperative to guide the design of ligands that boost the separation efficiency of the recovery process from fission products. However, limited studies have been made on the chemistry of these transport agents when complexing with transuranic elements in gas-phase where all surrounding factors are essentially excluded. This work investigates the reagent ligand complexations to transuranic and other metals and their dissociations in the gas phase. Comparisons are made between 4f and 5f elements and between ligands. Methods (120 words): (N,N-diisobutylcarbamoylmethyl)phenyloctylphosphine oxide (CMPO) and N,N,N',N'-tetraoctyldiglycolamide (TODGA) have been selected to complex with metal nitrates. The actinide americium and lanthanides neodymium, samarium, and europium were investigated as part of this work. The lanthanides were selected to act as size and electron configuration analogues of the minor actinides. Metal complexes with two ligands and two nitrates ([M(NO3)2(CMPO)2]+, for example) are studied in Bruker micrOTOF-Q II mass spectrometer equipped with collision-induced dissociation capability. The comparisons of the mass spectra are made in groups of homogenous ligands and mixed TODGA-CMPO ligands clusters. Comparisons are also made based on the complexed metals (Am and lanthanides). Preliminary data (300 words): Collision-induced dissociation mass spectrometry data are collected on two ligands complexed with metal nitrates where the two ligands are homogenous, with (CMPO)2 or (TODGA)2, or heterogeneous, with (TODGA)(CMPO). Several fragmentation patterns are observed among complexes with the CMPO ligand whereas the TODGA ligand commonly dissociates intact from the complex. Most of the metal complexes exhibit similar fragmentation patterns, but there are a few notable deviations in fragmentation patterns between the Am and Ln-bearing complexes. For the [M(CMPO)2(NO3)2]+ complexes, the initial loss of nitrate in the form of nitric acid is observed in all four complexes. However, [Am(CMPO)2(NO3)2]+ exhibits an additional fragmentation not found in the lanthanide complexes. Also, a significantly different ratio of the second nitric acid loss is found in the Am complex. These deviations may indicate the different interaction behaviors between actinides and lanthanides. The [M(TODGA)2(NO3)2]+ complexes exhibit the fragmentation as the loss of one TODGA ligand as an intact form and the loss of nitrate as nitric acid. The Am complex exhibits an additional fragmentation after losing the TOGDA ligand, which is not observed among the Ln complexes. The heterogeneous [M(TODGA)(CMPO)(NO3)2]+ complexes exhibit both similarities and differences between the Am and Ln complexes. For example, the heterogenous Am complex does not exhibit the loss of an intact TODGA ligand while all three Ln complexes do. This indicates that TODGA may bind more strongly to Am than Ln. Additionally, the intensity of the loss of CMPO ligand (as partially or whole) is found to be significantly larger than that of the loss of TODGA (as partially of whole) indicating that TODGA is bound to the metal significantly stronger than CMPO. Planned computational analysis will help understand the deviation in fragmentation behaviors between americium and lanthanide metal centers, or between TODGA and CMPO ligands. Novel aspect (20 words): Gas-phase actinide and lanthanide complex formation and fragmentation provide insight into the coordination environment differences of f-element metals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

Reynolds number influences in aeronautics

Reynolds number, a measure of the ratio of inertia to viscous forces, is a fundamental similarity parameter for fluid flows and therefore, would be expected to have a major influence in aerodynamics and aeronautics. Reynolds number influences are generally large, but monatomic, for attached laminar (continuum) flow; however, laminar flows are easily separated, inducing even stronger, non-monatomic, Reynolds number sensitivities. Probably the strongest Reynolds number influences occur in connection with transitional flow behavior. Transition can take place over a tremendous Reynolds number range, from the order of 20 x 10(exp 3) for 2-D free shear layers up to the order of 100 x 10(exp 6) for hypersonic boundary layers. This variability in transition behavior is especially important for complex configurations where various vehicle and flow field elements can undergo transition at various Reynolds numbers, causing often surprising changes in aerodynamics characteristics over wide ranges in Reynolds number. This is further compounded by the vast parameterization associated with transition, in that any parameter which influences mean viscous flow development (e.g., pressure gradient, flow curvature, wall temperature, Mach number, sweep, roughness, flow chemistry, shock interactions, etc.), and incident disturbance fields (acoustics, vorticity, particulates, temperature spottiness, even electro static discharges) can alter transition locations to first order. The usual method of dealing with the transition problem is to trip the flow in the generally lower Reynolds number wind tunnel to simulate the flight turbulent behavior. However, this is not wholly satisfactory as it results in incorrectly scaled viscous region thicknesses and cannot be utilized at all for applications such as turbine blades and helicopter rotors, nacelles, leading edge and nose regions, and High Altitude Long Endurance and hypersonic airbreathers where the transitional flow is an innately critical portion of the problem.

Bushnell, Dennis M.↗
NASA NTRS1995

Micromechanical models for textile composites

Numerical and analytical micromechanical models are presented to predict the thermoelastic behavior of a textile composite. In the numerical model, the unit-cell is discretized with finite elements, and periodic boundary conditions are imposed between opposite faces of the unit-cell. For a thin textile composite, stress gradients effects through the thickness are demonstrated. The consequent difference in the stiffness and strength behavior of thick and thin composites are discussed. The numerical model is implemented to predict 3-D thermo-elastic constants for a thick textile composite, and the plate thermo-mechanical properties for a thin textile composite. The numerical model is extended to compute the thermal residual microstresses due to processing to predict the composite failure envelopes. An analytical model - Selective Averaging Method (SAM) - is proposed, which is based on a judicious combination of stiffness and compliance averaging to estimate the 3-D elastic constants. Both the models are tested and verified for several examples by comparing the stiffness properties with elasticity solutions and available results.

Sankar, Bhavani V.↗
NASA NTRS2002

Techniques Optimized for Reducing Instabilities in Advanced Nickel-Base Superalloys for Turbine Blades

The High-Speed Research (HSR) Airfoil Alloy program developed fourth-generation single-crystal superalloys with up to an 85 F increase in creep rupture capability over current production airfoil alloys. Recent results have been generated at the NASA Glenn Research Center on these fourth-generation alloys, but in coated form, for subsonic turbine blade applications under NASA's Ultra-Efficient Engine Technology (UEET) Program. One goal for UEET is to optimize the airfoil alloy/thermal barrier coating system for 3100 F turbine inlet temperatures. The state-of-the art turbine blade airfoil system consists of a superalloy single crystal that provides the basic mechanical performance of the airfoil. A thermal barrier coating is used to reduce the temperature of the base superalloy, and a bondcoat is deposited between the base material and the thermal barrier coating. The bondcoat improves the oxidation and corrosion resistance of the base superalloy and improves the spallation resistance of the thermal barrier coating. A commercial platinum aluminide bondcoat was applied to the HSR-developed alloys, and a diffusion zone developed as a result of interaction between the bondcoat and the superalloy. Optimized strength is obtained for superalloys when the refractory element content is high and the limits of microstructural stability are approached or exceeded slightly. For fourthgeneration alloys, instability leads to the formation of topologically close packed (TCP) phases, which form internally in the superalloy, and a secondary reaction zone (SRZ), which forms under the diffusion zone. There was a concern that excessive quantities of either TCP or SRZ might decrease the mechanical properties of the superalloy, with SRZ thought to be particularly detrimental and its formation unpredictable. Thus, an SRZreduction effort was initiated in the NASA UEET Program so that methods developed during the HSR project could be optimized further to reduce or eliminate the SRZ. An SRZ is a three-phase constituent composed of TCP and stringers of gamma phase in a matrix of gamma prime. An incoherent grain boundary separates the SRZ from the gammagamma prime microstructure of the superalloy. The SRZ is believed to form as a result of local chemistry changes in the superalloy due to the application of the diffusion aluminide bondcoat. Locally high surface stresses also appear to promote the formation of the SRZ. Thus, techniques that change the local alloy chemistry or reduce surface stresses have been examined for their effectiveness in reducing SRZ. These SRZ-reduction steps are performed on the test specimen or the turbine blade before the bondcoat is applied. Stressrelief heat treatments developed at NASA Glenn have been demonstrated to reduce significantly the amount of SRZ that develops during subsequent high-temperature exposures. Stress-relief heat treatments reduce surface stresses by recrystallizing a thin surface layer of the superalloy. However, in alloys with very high propensities to form SRZ, stress relief heat treatments alone do not eliminate SRZ entirely. Thus, techniques that modify the local chemistry under the bondcoat have been emphasized and optimized successfully at Glenn. One such technique is carburization, which changes the local chemistry by forming submicron carbides near the surface of the superalloy. Detailed characterizations have demonstrated that the depth and uniform distribution of these carbides are enhanced when a stress relief treatment and an appropriate surface preparation are employed in advance of the carburization treatment. Even in alloys that have the propensity to develop a continuous SRZ layer beneath the diffusion zone, the SRZ has been completely eliminated or reduced to low, manageable levels when this combination of techniques is utilized. Now that the techniques to mitigate SRZ have been established at Glenn, TCP phase formation is being emphasized in ongoing work under the UEET Program. The limitsf stability of the fourth-generation alloys with respect to TCP phase formation are currently being defined along with high-temperature creep rupture properties. In addition, a regression model is being developed at Glenn for the prediction of the presence of TCP phase in the microstructure and SRZ under the diffusion zone. The model is based on a design-of-experiments methodology with emphasis on the potential synergistic effects of alloying elements.

MacKay, Rebecca A.↗
DOE OSTI2025

Combining four-point bending and corrosion to map stress-dependent corrosion susceptibility of 316H stainless steel in FLiNaK

Understanding the effect of stress on the corrosion of steels in molten salts is important for material screening and safety assessment of molten salt reactors. We present a method that combines four-point bending with corrosion testing, enabling the mapping of corrosion susceptibility as a function of local stress from a single specimen. Guided by finite element analysis, a four-point bending setup was used to bend a 316H stainless-steel bar under controlled elastic stress during a 100-hour exposure in FLiNaK at 700 °C. Post-exposure characterization using scanning electron microscopy and energy-dispersive X-ray spectroscopy revealed a pronounced correlation between stress and corrosion. Regions under tensile stress exhibited significantly deeper attack depths and greater chromium depletion along grain boundaries, whereas compressive zones showed shallower corrosion cracks. The stress effect is attributed to stress concentration under tensile loading, which promotes chromium dissolution and crack propagation. This technique can be readily applied to other structural alloys, enabling quantitative measurement of corrosion susceptibility as a function of local stress.

316H↗
NASA NTRS1990

Evaluation of an adaptive unstructured remeshing technique for integrated fluid-thermal-structural analysis

An adaptive unstructured remeshing technique is evaluated for integrated fluid-thermal-structural analysis. The technique is combined with the finite element method to solve: (1) the Navier-Stokes equations for high-speed compressible flow; (2) the energy equation for the structural-thermal response; and (3) the quasi-static equilibrium equations for the structural response. The remeshing technique and the analysis solution procedure are described. The effectiveness of the approach is evaluated with two application studies. The flow analysis of Mach 8 shock-shock interference on a three-inch-diameter cylinder is used as the first application study to demonstrate the capability of the remeshing technique and to examine proper remeshing indicators for the inviscid and boundary layer regions. The applicability of the approach for the thermal and structural analyses of the structure is evaluated in the second application study of a 0.25-inch-diameter convectively cooled leading edge subjected to intense aerodynamic heating. Issues associated with remeshing indicators for thermal stress problems are identified.

Dechaumphai, Pramote↗
NASA NTRS1990

Evaluation of an adaptive unstructured remeshing technique for integrated fluid-thermal-structural analysis

An adaptive unstructured remeshing technique is evaluated for integrated fluid-thermal-structural analysis. The technique is combined with the finite element method to solve: (1) the Navier-Stokes equations for high-speed compressible flow, (2) the energy equation for the structural-thermal response, and (3) the quasi-static equilibrium equations for the structural response. The remeshing technique and the analysis solution procedure are described. The effectiveness of the approach is evaluated with two applications studies. The flow analysis of Mach 8 shock-shock interference on a three-inch-diameter cylinder is used as the first application study to demonstrate the capability of the remeshing technique and to examine proper remeshing indicators for the inviscid and boundary layer regions. The applicability of the approach for the thermal and structural analyses of the structure is evaluated in the second application study of a 0.25-inch-diameter convectively cooled leading edge subjected to intense aerodynamic heating. Issues associated with remeshing indicators for thermal stress problems are identified.

Dechaumphai, Pramote↗
DOE OSTI2023

Meshfree Multiphysics Damage Modeling of Li-ion Battery Materials

Through repeated charging and discharging cycles, the electrodes of a Li-ion battery experience significant swelling and contraction due to the movement of lithium, also known as intercalation. This research focuses specifically on the chemo-mechanical cracking in the cathode. Because cathode particles are comprised of many randomly-oriented grains, which have highly anisotropic material properties, the expansion and contraction is very non-uniform. As a result, stress concentrations tend to form between grains, which necessitates the modeling of crack propagation largely along grain boundaries and material interfaces. Chemo-mechanical damage models are generally simulated by one of two methods: the cohesive zone model (CZM) and the continuous damage model (CDM). The CZM is more accurate at capturing the sharp discontinuities of a crack but is very computationally expensive and intractable for large-scale models as a result. Conversely, the CDM is easily computed but not well-suited to easily allow for discontinuous field variables, which are inherent across a crack. This study aims to improve the CDM's ability to capture discontinuous cracks. Current versions of the CDM use the finite element method (FEM), which is one of the most widely used approaches for spatial discretization. This research investigates the use of the Reproducing Kernel Particle Method (RKPM), a meshfree method, for spatial discretization and aims to achieve a coupled chemo-mechanical crack propagation model that enhances accuracy while maintaining high computational efficiency. This model features a fully coupled, iterative electrochemistry solution, which informs the meshfree damage model of impending crack formation. To gain a better understanding of how the meshfree model compares to the current FEM model in capturing cathode crack propagation, further investigation is needed.

damage modeling↗