Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Element Modelling”

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

Quantifying spectral albedo effects on bifacial photovoltaic module measurements and system model predictions

We provide a comprehensive analysis of the effect of spectral albedo on photovoltaic (PV) module measurements and system model predictions. We demonstrate how to account for albedo in indoor bifacial device measurements by adjusting the applied irradiance using the scaled rear irradiance method, exemplified on fabricated silicon heterojunction (SHJ) modules. System model performance is studied using a detailed 3D finite-element model, DUET, for fixed-tilt and horizontal single-axis tracked (SAT) arrays between 15 and 75°N. Spectral effects cause variations in measured SHJ module short-circuit current up to 2% and efficiency variation up to 0.3% abs. We further demonstrate that rear-side spectral mismatch factors (SMMs) resulting from including or omitting spectral albedo in PV system modeling vary between ±13%, while total (front+rear) SMMs vary up to 3%, depending on the deployment configuration and latitude. SAT array SMMs are weakly correlated with latitude, while fixed-tilt array SMMs increase with latitude, driven by an increasing proportion of ground-reflected light on the front-side of modules. Ground-reflections can constitute between 2% and 32% of total incident module irradiance, with notably high (>10%) contributions for fixed-tilt arrays at high latitude. Effects of spectral albedo are most significant for: (1) fixed-tilt deployments at high latitudes, (2) wide bandgap technologies such as perovskite and cadmium telluride cells, (3) albedos which vary steeply over the technology's absorption range, and (4) high albedo ground covers. Overall, we demonstrate that omitting spectral albedo effects can result in PV measurement and system-level modeling uncertainties on the order of several percent in these cases.

14 SOLAR ENERGY↗

Dissipation and Bathymetric Sensitivities in an Unstructured Mesh Global Tidal Model

Abstract The mechanisms and geographic distribution of global tidal dissipation in barotropic tidal models are examined using a high resolution unstructured mesh finite element model. Mesh resolution varies between 2 and 25 km and is especially focused on inner shelves and steep bathymetric gradients. Tidal response sensitivities to bathymetric changes are examined to put into context response sensitivities to frictional processes. We confirm that the Ronne Ice Shelf dramatically affects Atlantic tides but also find that bathymetry in the Hudson Bay system is a critical control. We follow a sequential frictional parameter optimization process and use TPXO9 data‐assimilated tidal elevations as a reference solution. From simulated velocities and depths, dissipation within the global model is estimated and allows us to pinpoint dissipation at high resolution. Boundary layer dissipation is extremely focused with 1.4% of the ocean accounting for 90% of the total. Internal tide friction is much more distributed with 16.7% of the ocean accounting for 90% of the total. Often highly regional dissipation can impact basin‐scale and even ocean wide tides. Optimized boundary layer friction parameters correlate very well with the physical characteristics of the locality with high friction factors associated with energetic tidal regions, deep ocean island chains, and ice covered areas. Global complex M 2 tide errors are 1.94 cm in deep waters. Total global boundary layer and internal tide dissipation are estimated, respectively, at 1.83 and 1.49 TW. This continues the trend in the literature toward attributing more dissipation to internal tides.

54 ENVIRONMENTAL SCIENCES↗

Broadband/Wideband Magnetoelectric Response

A broadband/wideband magnetoelectric (ME) composite offers new opportunities for sensing wide ranges of both DC and AC magnetic fields. The broadband/wideband behavior is characterized by flat ME response over a given AC frequency range and DC magnetic bias. The structure proposed in this study operates in the longitudinal-transversal (L-T) mode. In this paper, we provide information on (i) how to design broadband/wideband ME sensors and (ii) how to control the magnitude of ME response over a desired frequency and DC bias regime. A systematic study was conducted to identify the factors affecting the broadband/wideband behavior by developing experimental models and validating them against the predictions made through finite element modeling. A working prototype of the sensor with flat bands for both DC and AC magnetic field conditions was successfully obtained. These results are quite promising for practical applications such as current probe, low-frequency magnetic field sensing, and ME energy harvester.

Park, Chee-Sung↗

NASTRAN thermal analyzer: A general purpose finite element heat transfer computer program

The program not only can render temperature distributions in solids subjected to various thermal boundary conditions, including effects of diffuse-gray thermal radiation, but is fully compatible in capacity and in the finite-element model representation with that of its structural counterpart in the NASTRAN system. The development history of the finite-element approach for determining temperatures is summarized. The scope of analysis capability, program structure, features, and limitations are given with the objective of providing NASTRAN users with an overall veiw of the NASTRAN thermal analyzer.

Lee, H.↗

The LDCM actuator for vibration suppression

A linear dc motor (LDCM) has been proposed as an actuator for the COFS I mast and the COFS program ground test Mini-Mast. The basic principles of operation of the LDCM as an actuator for vibration suppression in large flexible structures are reviewed. Because of force and stroke limitations, control loops are required to stabilize the actuator, which results in a non-standard actuator-plant configuration. A simulation model that includes LDCM actuator control loops and a finite element model of the Mast is described, with simulation results showing the excitation capability of the actuator.

Ide, Eric N.↗

A study of fiber volume fraction effects in notched unidirectional SCS-6/Ti-15V-3Cr-3Al-3Sn composite

Notched unidirectional SCS-6/Ti-15-3 composite of three different fiber volume fractions (vf = 0.15, 0.37, and 0.41) was investigated for various room temperature microstructural and material properties including: fatigue crack initiation, fatigue crack growth, and fracture toughness. While the matrix hardness is similar for all fiber volume fractions, the fiber/matrix interfacial shear strength and matrix residual stress increases with fiber volume fraction. The composite fatigue crack initiation stress is shown to be matrix controlled and occurs when the net maximum matrix stress approaches the endurance limit stress of the matrix. A model is presented which includes residual stresses and presents the composite initiation stress as a function of fiber volume fraction. This model predicts a maximum composite initiation stress at vf approximately 0.15 which agrees with the experimental data. The applied composite stress levels were increased as necessary for continued crack growth. The applied Delta(K) values at crack arrest increase with fiber volume fraction by an amount better approximated using an energy based formulation rather than when scaled linear with modulus. After crack arrest, the crack growth rate exponents for vf37 and vf41 were much lower and toughness much higher, when compared to the unreinforced matrix, because of the bridged region which parades with the propagating fatigue crack. However, the vf15 material exhibited a higher crack growth rate exponent and lower toughness than the unreinforced matrix because once the bridged fibers nearest the crack mouth broke, the stress redistribution broke all bridged fibers, leaving an unbridged crack. Degraded, unbridged behavior is modeled using the residual stress state in the matrix ahead of the crack tip. Plastic zone sizes were directly measured using a metallographic technique and allow prediction of an effective matrix stress intensity which agrees with the fiber pressure model if residual stresses are considered. The sophisticated macro/micro finite element models of the 0.15 and 0.37 fiber volume fractions presented show good agreement with experimental data and the fiber pressure model when an estimated effective fiber/matrix debond length is used.

Covey, Steven J.↗

Comparison of numerical simulation and experimental data for steam-in-place sterilization

A complex problem involving convective flow of a binary mixture containing a condensable vapor and noncondensable gas in a partially enclosed chamber was modelled and results compared to transient experimental values. The finite element model successfully predicted transport processes in dead-ended tubes with inside diameters of 0.4 to 1.0 cm. When buoyancy driven convective flow was dominant, temperature and mixture compositions agreed with experimental data. Data from 0.4 cm tubes indicate diffusion to be the primary air removal method in small diameter tubes and the diffusivity value in the model to be too large.

Young, Jack H.↗

Design and Stress Analysis of Low-Noise Adjusted Bearing Contact Spiral Bevel Gears

An integrated computerized approach for design and stress analysis of low-noise spiral bevel gear drives with adjusted bearing contact has been developed. The computation procedure is an iterative process, requiring four separate steps that provide: (a) a parabolic function of transmission errors that is able to reduce the effect of errors of alignment, and (b) reduction of the shift of bearing contact caused by misalignment. Application of finite element analysis permits the contact and bending stresses to be determined and investigate the formation of the bearing contact. The design of finite element models and boundary conditions is automated and does not require an intermediate CAD computer program. A commercially available finite element analysis computer program with contact capability was used to conduct the stress analysis. The theory developed is illustrated with numerical examples.

Litvin, Faydor L.↗

Explicit Pore Pressure Material Model in Carbon-Cloth Phenolic

An explicit material model that uses predicted pressure in the pores of a carbon-cloth phenolic (CCP) composite has been developed. This model is intended to be used within a finite-element model to predict phenomena specific to CCP components of solid-fuel-rocket nozzles subjected to high operating temperatures and to mechanical stresses that can be great enough to cause structural failures. Phenomena that can be predicted with the help of this model include failures of specimens in restrained-thermal-growth (RTG) tests, pocketing erosion, and ply lifting

Gutierrez-Lemini, Danton↗

Comparison of Multiscale Method of Cells-Based Models for Predicting Elastic Properties of Filament Wound C/C-SiC

Three different multiscale models, based on the method of cells (generalized and high fidelity) micromechanics models were developed and used to predict the elastic properties of C/C-SiC composites. In particular, the following multiscale modeling strategies were employed: Concurrent multiscale modeling of all phases using the generalized method of cells, synergistic (two-way coupling in space) multiscale modeling with the generalized method of cells, and hierarchical (one-way coupling in space) multiscale modeling with the high fidelity generalized method of cells. The three models are validated against data from a hierarchical multiscale finite element model in the literature for a repeating unit cell of C/C-SiC. Furthermore, the multiscale models are used in conjunction with classical lamination theory to predict the stiffness of C/C-SiC plates manufactured via a wet filament winding and liquid silicon infiltration process recently developed by the German Aerospace Institute.

multiscale modeling↗

Magnetic Shield Design Modeling and Validation for SWOT Spacecraft Ka-Band Extended Interaction Klystron

Two Extended Interaction Klystrons (EIKs) containing strong permanent magnets were modeled magnetically in a representative spacecraft geometry using commercial finite element modeling techniques and were validated against measurements made at varying distances. Initial modeling results for the 63 A-m2 (amperes per square meter) dipole moment magnets showed that magnetic shields would be necessary in order to meet magnetic field requirements for the Surface Water and Ocean Topography (SWOT) spacecraft, which contains components that are susceptible to external DC (Direct Current) magnetic fields. JPL (Jet Propulsion Laboratory) and the EIK vendor proposed cold rolled steel and mu-metal as potential shield materials along with proposed thicknesses of 0.5 millimeters and 1.5 millimeters. Magnetic shields made from each of these materials were designed and modeled in software, taking high-field saturation into account. Prototype magnetic shields with these parameters were then built, measured with an existing EIK, and compared against modeling results. For single-axis field measurements along the dipole axis, modeling results were within 7 gauss of the measured values at 10 centimeters from the magnet, and converged to less than 1.5 gauss at distances greater than 14 centimeters from the magnet. Three-axis field measurements at locations of interest showed that model correlation improved to within 4 gauss at 11 centimeters and 2 gauss for distances ranging between 15 centimeters and 36 centimeters.

Gonzales, Edward↗

Aerodynamic and Structural Evaluation of an SMA Slat-Cove Filler Using Computational and Experimental Tools at Model Scale

Transport class aircraft produce a significant amount of airframe noise during approach and landing due to exposed geometric discontinuities that are hidden during cruise. The leading-edge slat is a primary contributor to this noise. In previous work, use of a slat-cove filler (SCF) has proven to reduce airframe noise by filling the cove aft of the slat, eliminating the circulation region within the cove. The goal of this work is to extend and improve upon past experimental and computational efforts on the evaluation of a scaled high-lift wing with a super elastic shape memory alloy (SMA) SCF. Recent turbulence measurements of the Texas A&M University 3ft-by-4ft wind tunnel allow for more accurate representation of the flow through the test section in computational fluid dynamics (CFD) analysis. The finite volume models used in CFD analysis are coupled to structural finite element models using a framework compatible with an SMA constitutive model and significant deformation, enabling fluid-structure interaction (FSI) analysis of the SCF. Both fully-deployed and retraction/deployment cases are considered. The displacement of the SCF on the experimental model is measured at various stages of retraction/deployment using a laser displacement sensor and digital image correlation system. Due to a lack of structural stiffness in the 3D-printed plastic slat during retraction and SCF stowage, a rigid steel slat is incorporated into the physical model and preliminary wind tunnel tests are conducted at multiple angles of attack in both retracted and deployed configurations.

Scholten, William↗

Magnetic Shield Design Modeling and Validation for SWOT Spacecraft Ka-Band Extended Interaction Klystron

Two Extended Interaction Klystrons (EIKs) containing strong permanent magnets were modeled magnetically in a representative spacecraft geometry using commercial finite element modeling techniques and were validated against measurements made at varying distances. Initial modeling results for the 63 A-m2 dipole moment magnets showed that magnetic shields would be necessary in order to meet magnetic field requirements for the Surface Water and Ocean Topography (SWOT) spacecraft, which contains components that are susceptible to external DC magnetic fields. JPL and the EIK vendor proposed cold rolled steel and mu-metal as potential shield materials along with proposed thicknesses of 0.5 mm and 1.5 mm. Magnetic shields made from each of these materials were designed and modeled in software, taking highfield saturation into account. Prototype magnetic shields with these parameters were then built, measured with an existing EIK, and compared against modeling results. For single-axis field measurements along the dipole axis, modeling results were within 7 gauss of the measured values at 10 cm from the magnet, and converged to less than 1.5 gauss at distances greater than 14 cm from the magnet. Three-axis field measurements at locations of interest showed that model correlation improved to within 4 gauss at 11 cm and 2 gauss for distances ranging between 15 cm and 36 cm.

Gonzales, Edward↗

Lattice Structured Lightweight Structural Materials

The development of lightweight structural materials is crucial for enhancing the performance and deployment feasibility of fission batteries. This study aims to produce lightweight structural materials whose strength-to-weight ratios exceed those of current widely used structural materials. To achieve this, advanced modeling and simulation tools were employed to design lattice structures with different lattice parameters and different lattice types. A process was successfully developed for transforming lattice-structured models into Multiphysics Object Oriented Simulation Environment (MOOSE) inputs. Finite element modeling (FEM) was used to simulate the uniaxial tensile testing of the lattice-structured parts to investigate the stress distribution at a given displacement. The modeling results showed that the lattice-structured sample displayed a lower Young’s modulus in comparison to the solid material; the increase in solid shell thickness and blend radius enhances the mechanical performance; and the effect of unit cell size on macro scale stress is minimal. Tensile testing was conducted on the solid and lattice-structured materials fabricated by laser powder bed fusion (LPBF) additive manufacturing. The experimental results agreed well with the model prediction. The approach of using modeling as a guiding tool for preliminary material design can significantly save time and cost for new material development.

lightweight material↗

Cohesive Modeling of Transverse Cracking in Laminates with a Single Layer of Elements per Ply

This study aims to bridge the gap between classical understanding of transverse cracking in cross-ply laminates and recent computational methods for the modeling of progressive laminate failure. Specifically, the study investigates under what conditions a finite element model with cohesive X-FEM cracks can reproduce the in situ effect for the ply strength. It is shown that it is possible to do so with a single element across the thickness of the ply, provided that the interface stiffness is properly selected. The optimal value for this interface stiffness is derived with an analytical shear lag model. It is also shown that, when the appropriate statistical variation of properties has been applied, models with a single element through the thickness of a ply can predict the density of transverse matrix cracks

VanDerMeer, Frans P.↗

Hierarchical Poly Tree Configurations for the Solution of Dynamically Refined Finte Element Models

This paper demonstrates how a multilevel substructuring technique, called the Hierarchical Poly Tree (HPT), can be used to integrate a localized mesh refinement into the original finite element model more efficiently. The optimal HPT configurations for solving isoparametrically square h-, p-, and hp-extensions on single and multiprocessor computers is derived. In addition, the reduced number of stiffness matrix elements that must be stored when employing this type of solution strategy is quantified. Moreover, the HPT inherently provides localize 'error-trapping' and a logical, efficient means with which to isolate physically anomalous and analytically singular behavior.

Gute, G. D.↗

A micromechanics model for bread dough

The mechanical behaviour of dough and gluten was studied in an effort to investigate whether bread dough can be treated as a two phase (starch and gluten) composite material. The dough and gluten show rate dependent behaviour under tension, compression and shear tests, and non-linear unloading-reloading curves under cyclic compression tests. There is evidence from cryo-Scanning Electron Microscopy (SEM) that damage in the form of debonding between starch and gluten occurs when the sample is stretched. A composite finite element model was developed using starch as filler and gluten as matrix. The interaction between the starch and gluten was modelled as cohesive contact. The finite element analysis predictions agree with trends seen in experimental test data on dough and gluten, further evidence that debonding of starch and gluten is a possible damage mechanism in dough.

36 MATERIALS SCIENCE↗

Hydrogeologic Modeling at the Sylvania Corning FUSRAP Site - 13419

A comparative groundwater hydrogeologic modeling analysis is presented herein to simulate potential contaminant migration pathways in a sole source aquifer in Nassau County, Long Island, New York. The source of contamination is related to historical operations at the Sylvania Corning Plant ('Site'), a 9.49- acre facility located at 70, 100 and 140 Cantiague Rock Road, Town of Oyster Bay in the westernmost portion of Hicksville, Long Island. The Site had historically been utilized as a nuclear materials manufacturing facility (e.g., cores, slug, and fuel elements) for reactors used in both research and electric power generation in early 1950's until late 1960's. The Site is contaminated with various volatile organic and inorganic compounds, as well as radionuclides. The major contaminants of concern at the Site are tetrachloroethene (PCE), trichloroethene (TCE), nickel, uranium, and thorium. These compounds are present in soil and groundwater underlying the Site and have migrated off-site. The Site is currently being investigated as part of the Formerly Utilized Sites Remedial Action Program (FUSRAP). The main objective of the current study is to simulate the complex hydrogeologic features in the region, such as numerous current and historic production well fields; large, localized recharge basins; and, multiple aquifers, and to assess potential contaminant migration pathways originating from the Site. For this purpose, the focus of attention was given to the underlying Magothy formation, which has been impacted by the contaminants of concern. This aquifer provides more than 90% of potable water supply in the region. Nassau and Suffolk Counties jointly developed a three-dimensional regional groundwater flow model to help understand the factors affecting groundwater flow regime in the region, to determine adequate water supply for public consumption, to investigate salt water intrusion in localized areas, to evaluate the impacts of regional pumping activity, and to better understand the contaminant transport and fate mechanisms through the underlying aquifers. This regional model, developed for the N.Y. State Department of Environmental Conservation (NYSDEC) by Camp Dresser and McKee (CDM), uses the finite element model DYNFLOW developed by CDM, Cambridge, Massachusetts. The coarseness of the regional model, however, could not adequately capture the hydrogeologic heterogeneity of the aquifer. Specifically, the regional model did not adequately capture the interbedded nature of the Magothy aquifer and, as such, simulated particles tended to track down-gradient from the Site in relatively straight lines while the movement of groundwater in such a heterogeneous aquifer is expected to proceed along a more tortuous path. This paper presents a qualitative comparison of site-specific groundwater flow modeling results with results obtained from the regional model. In order to assess the potential contaminant migration pathways, a particle tracking method was employed. Available site-specific and regional hydraulic conductivity data measured in-situ with respect to depth and location were incorporated into the T-PROG module in GMS model to define statistical variation to better represent the actual stratigraphy and layer heterogeneity. The groundwater flow characteristics in the Magothy aquifer were simulated with the stochastic hydraulic conductivity variation as opposed to constant values as employed in the regional model. Contaminant sources and their exact locations have been fully delineated at the Site during the Remedial Investigation (RI) phase of the project. Contaminant migration pathways originating from these source locations at the Site are qualitatively traced within the sole source aquifer utilizing particles introduced at source locations. Contaminant transport mechanism modeled in the current study is based on pure advection (i.e., plug flow) and mechanical dispersion while molecular diffusion effects are neglected due to relatively high groundwater velocities encountered in the aquifer. In addition, fate of contaminants is ignored hereby to simulate the worst-case scenario, which considers the contaminants of concern as tracer-like compounds for modeling purposes. The results of the modeling analysis are qualitatively compared with the County's regional model, and patterns of contaminant migration in the region are presented. (authors)

54 ENVIRONMENTAL SCIENCES↗