Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element simulation”

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 289 records · Page 16

Microstructure and mesh sensitivities of mesoscale surrogate driving force measures for transgranular fatigue cracks in polycrystals

The number of cycles required to form and grow microstructurally small fatigue cracks in metals exhibits substantial variability, particularly for low applied strain amplitudes. This variability is commonly attributed to the heterogeneity of cyclic plastic deformation within the microstructure, and presents a challenge to minimum life design of fatigue resistant components. Our paper analyzes sources of variability that contribute to the driving force of transgranular fatigue cracks within nucleant grains. We also employ crystal plasticity finite element simulations that explicitly render the polycrystalline microstructure and Fatigue Indicator Parameters (FIPs) averaged over different volume sizes and shapes relative to the anticipated fatigue damage process zone. Volume averaging is necessary to both achieve description of a finite fatigue damage process zone and to regularize mesh dependence in simulations. Furthermore, results from constant amplitude remote applied straining are characterized in terms of the extreme value distributions of volume averaged FIPs. Grain averaged FIP values effectively mitigate mesh sensitivity, but they smear out variability within grains. Furthermore, volume averaging over bands that encompass critical transgranular slip planes appear to present the most attractive approach to mitigate mesh sensitivity while preserving variability within grains.

36 MATERIALS SCIENCE↗

Accumulation of collagen molecular unfolding is the mechanism of cyclic fatigue damage and failure in collagenous tissues

Overuse injuries to dense collagenous tissues are common, but their etiology is poorly understood. The predominant hypothesis that micro-damage accumulation exceeds the rate of biological repair is missing a mechanistic explanation. Here, we used collagen hybridizing peptides to measure collagen molecular damage during tendon cyclic fatigue loading and computational simulations to identify potential explanations for our findings. Our results revealed that triple-helical collagen denaturation accumulates with increasing cycles of fatigue loading, and damage is correlated with creep strain independent of the cyclic strain rate. Finite-element simulations demonstrated that biphasic fluid flow is a possible fascicle-level mechanism to explain the rate dependence of the number of cycles and time to failure. Molecular dynamics simulations demonstrated that triple-helical unfolding is rate dependent, revealing rate-dependent mechanisms at multiple length scales in the tissue. The accumulation of collagen molecular denaturation during cyclic loading provides a long-sought “micro-damage” mechanism for the development of overuse injuries.

59 BASIC BIOLOGICAL SCIENCES↗

The direct-current response of electrically conducting fractures excited by a grounded current source

Hydraulic fracture stimulation of low permeability reservoir rocks is an established and cross–cutting technology for enhancing hydrocarbon production in sedimentary formations and increasing heat exchange in crystalline geothermal systems. Whereas the primary measure of success is the ability to keep the newly generated fractures sufficiently open, long–term reservoir management requires a knowledge of the spatial extent, morphology, and distribution of the fractures — knowledge primarily informed by microseismic and ground deformation monitoring. To minimize the uncertainty associated with interpreting such data, we investigate through numerical simulation the usefulness of direct-current (DC) resistivity data for characterizing subsurface fractures with elevated electrical conductivity by considering a geophysical experiment consisting of a grounded current source deployed in a steel cased borehole. In doing so, the casing efficiently energizes the fractures with steady current. Finite element simulations of this experiment for a horizontal well intersecting a small set of vertical fractures indicate that the fractures manifest electrically in (at least) two ways: (1) a local perturbation in electric potential proximal to the fracture set, with limited farfield expression and (2) an overall reduction in the electric potential along the borehole casing due to enhanced current flow through the fractures into the surrounding formation. The change in casing potential results in a measurable effect that can be observed far from fractures themselves. Under these conditions, our results suggest that farfield, timelapse measurements of DC potentials can be interpreted by simple, linear inversion for a Coulomb charge distribution along the borehole path, including a local charge perturbation due to the fractures. As a result, this approach offers an inexpensive method for detecting and monitoring the time-evolution of electrically conducting fractures while ultimately providing an estimate of their effective conductivity — the latter providing an important measure independent of seismic methods on fracture shape, size, and hydraulic connectivity.

30 DIRECT ENERGY CONVERSION↗

Prediction of particle orientation in simple upsetting process of NdFeB magnets

The magnetic properties of NdFeB magnets are strongly affected by crystallographic texture which is highly associated with particle orientation. This study proposed a method for predicting the particle orientation in the simple upsetting process of NdFeB magnets. The method is based on finite element simulation with flow net analysis. The magnets in a cylindrical form were compressed by two flat dies in a chamber filled with argon at 750°C. Three forming speeds were taken into account in order to obtain flow stress curves used in simulations. The micrographs of the cross sections of the deformed magnets show that the particle deformation significantly increases with the compression. The phenomenon was also predicted by the proposed method. Both simulated and experimental results show that the inhomogeneity of the texture of the NdFeB magnets can be increased by the simple upsetting process. The predicted particle orientations were in a good agreement with those examined in the deformed magnets. The proposed method for predicting particle orientations can also be used in other forming processes of NdFeB magnets.

36 MATERIALS SCIENCE↗

An Evaluation of Different Plasticity and Failure Laws in Simulating Puncture in 7075-T651 Aluminum

In this work, various material models were studied for their ability to simulate puncture in a thin aluminum 7075-T651 plate due to low-velocity probe impact. Material models were generated by mixing and matching various work hardening laws with different failure criteria, and several hybrid material models were investigated. Finite element simulations of aluminum impact-response, based on each material model, were employed to predict the energy required for puncture and final plate tear-out geometry. Probes of different size and shape were used to impose various loading regimes, and numerical predictions were compared to experimental results from a previous study. It was found that no single combination of hardening and failure laws yielded universally accurate data, but that several material models could be used more reliably than others. Further, the importance of obtaining unique parameter-sets for work-hardening and failure criteria was illustrated.

36 MATERIALS SCIENCE↗

Finite element-integral simulation of static and flight fan noise radiation from the JT15D turbofan engine

An iterative finite element integral technique is used to predict the sound field radiated from the JT15D turbofan inlet. The sound field is divided into two regions: the sound field within and near the inlet which is computed using the finite element method and the radiation field beyond the inlet which is calculated using an integral solution technique. The velocity potential formulation of the acoustic wave equation was employed in the program. For some single mode JT15D data, the theory and experiment are in good agreement for the far field radiation pattern as well as suppressor attenuation. Also, the computer program is used to simulate flight effects that cannot be performed on a ground static test stand.

Baumeister, K. J.↗

Blockage effects in the chemotaxis of diffusiophoretic particles

Transport mechanisms at the micro- and nano-scale play an essential role in regulating intracellular organization. Recent work indicates that directed motion of constituents inside cells can emerge through diffusiophoretic transport, in which colloidal particles move under the influence of chemical gradients. Here, we examine how blockers—passive or actively consuming—reshape those gradients and thereby influence the motion of diffusiophoretic particles. By combining analytical solutions with finite element simulations, we first show that a single blocker can distort a background gradient enough to create or eliminate stagnation points, significantly modifying particle transport. We then introduce a second, explicitly sized blocker at one of these stagnation points and measure how its finite radius alters the diffusiophoretic velocity field for a test particle. Even moderate changes in the second blockers size can cause noticeable shifts in the substrate distribution, highlighting the importance of accounting for explicit particle radii under crowded or consumption-driven conditions. Our findings underscore that subtle geometric variations—such as the radii and positions of two or more blockers—can profoundly affect diffusiophoretic motion, providing a more complete picture of how blocking and crowding phenomena shape intracellular transport.

Song, Zehao [Northwestern Univ., Evanston, IL (Uni↗

Heat-affected zone liquation crack on resistance spot welded TWIP steels

In this study, the heat affected zone (HAZ) liquation crack and segregation behavior of the resistance spot welded twinning induced plasticity (TWIP) steel have been reported. Cracks appeared in the post-welded joints that originated at the partially melted zone (PMZ) and propagated from the PMZ through the heat affected zone (HAZ) to the base metal (BM). The crack length and crack opening widths were observed increasing with heat input; and the welding current was identified to be the most influencing parameter for crack formation. Cracks appeared at the PMZ when nugget diameter reached at 4.50 mm or above; and the liquation cracks were found to occur along two sides of the notch tip in the sheet direction rather than in the electrode direction. Cracks were backfilled with the liquid films which has lamellar structure and supposed to be the eutectic constituent. Co-segregation of alloy elements such as, C and Mn were detected on the liquid films by electron-probe microanalysis (EPMA) line scanning and element map which suggests that the liquid film was enrich of Mn and C. The eutectic constituent was identified by analyzing the calculated phase diagram along with thermal temperature history of finite element simulation. Preliminary experimental results showed that cracks have less/no significant effect on the static cross-tensile strength (CTS) and the tensile-shear strength (TSS). In addition, possible ways to avoid cracking were discussed. - Highlights: • The HAZ liquation crack during resistance spot welding of TWIP steel was examined. • Cracks were completely backfilled and healed with divorced eutectic secondary phase. • Co-segregation of C and Mn was detected in the cracked zone. • Heat input was the most influencing factor to initiate liquation crack. • Cracks have less/no significant effect on static tensile properties.

36 MATERIALS SCIENCE↗

Simulation of emergent compaction banding fronts caused by frictional boundaries

This study examines the role of boundary friction in promoting heterogeneous compaction in soft rock specimens loaded at high confining pressure outside the domain of compaction localisation. An elastoplastic constitutive model characterised by tunable hardening/softening behaviour is used to conduct the analyses. Finite-element simulations suggest that material instability is a non-necessary condition for the emergence of compaction fronts. Such fronts propagated as a result of a severe deviation in the local responses induced by frictional constraints. Overall, these findings suggest that boundary effects can bias the assessment of the extent of the compaction localisation domain. Experimental countermeasures and informed model calibration procedures are therefore necessary to minimise such bias and enable more accurate predictions of soft rock compaction.

42 ENGINEERING↗

Plastic flow anisotropy drives shear fracture

Abstract Fracture of initially crack-free bodies often occurs due to plastic instabilities known as shear bands. Previous computer simulations advanced a myriad of mechanisms to rationalize shear banding. However, they were restricted to planar geometries. We investigate the relevance of anisotropic plasticity by picking an axisymmetric tensile test rig, in which shear localization is rarely observed. The three-dimensional finite-element simulations of shear banding in this type of specimens are the first of their kind. The micromechanical modeling covers a range of competing mechanisms believed to be responsible for such failure. We show that anisotropic plasticity can effectively trigger shear bands thereby causing failure of ductile solids. Our results enable shear fracture to be rationalized in ductile rocks and mitigated against in designing advanced materials.

60 APPLIED LIFE SCIENCES↗

Measurement of a radial flow profile with eddy current flow meters and deep neural networks

Eddy current flow meters (ECFMs) measure flows of conductive fluids. Recent interest in ECFMs has increased due to applications in advanced nuclear reactors. ECFMs are well suited for such applications, as they can provide non-invasive measurements of flow in fluids that are often difficult to measure. Traditionally, ECFMs are operated using an alternating current at a single frequency, limiting ECFMs to measure average fluid velocities, blockages, or voids. Here, we expand the capabilities of ECFMs by measuring the fluid radial velocity profile of liquid mercury. To accomplish this, we made several ECFM sensitivity measurements at a range of frequencies. Different frequencies vary the electromagnetic skin depth of the device. By adjusting frequencies, we probed the fluid velocity at various radial locations and constructed a flow-velocity profile. The relationship between the ECFM measurements and velocity profile is nonlinear and requires solving an inverse problem. Using electromagnetic finite-element simulations to train a deep neural network (DNN), we created a model that provides a stable general relationship between the sensitivity measurements of an ECFM and the fluid velocity profile. Using ECFM measurements of liquid mercury, our DNN model calculates a flow profile that agrees well with computational fluid dynamics (CFD) simulations. This technique has potential to improve flow monitoring for optimization, safe operation of conductive fluid loops, and/or validating complex CFD models.

47 OTHER INSTRUMENTATION↗

Magnetic tweezers with high permeability electromagnets for fast actuation of magnetic beads

As a powerful and versatile scientific instrument, magnetic tweezers have been widely used in biophysical research areas, such as mechanical cell properties and single molecule manipulation. If one wants to steer bead position, the nonlinearity of magnetic properties and the strong position dependence of the magnetic field in most magnetic tweezers lead to quite a challenge in their control. In this article, we report multi-pole electromagnetic tweezers with high permeability cores yielding high force output, good maneuverability, and flexible design. For modeling, we adopted a piece-wise linear dependence of magnetization on field to characterize the magnetic beads. We implemented a bi-linear interpolation of magnetic field in the work space, based on a lookup table obtained from finite element simulation. The electronics and software were custom-made to achieve high performance. In addition, the effects of dimension and defect on structure of magnetic tips also were inspected. In a workspace with size of 0.1 × 0.1 mm{sup 2}, a force of up to 400 pN can be applied on a 2.8 μm superparamagnetic bead in any direction within the plane. Because the magnetic particle is always pulled towards a tip, the pulling forces from the pole tips have to be well balanced in order to achieve control of the particle’s position. Active video tracking based feedback control is implemented, which is able to work at a speed of up to 1 kHz, yielding good maneuverability of the magnetic beads.

36 MATERIALS SCIENCE↗

Thermal analysis of an indirectly heat pulsed non-volatile phase change material microwave switch

We show the finite element simulation of the melt/quench process in a phase change material (GeTe, germanium telluride) used for a radio frequency switch. The device is thermally activated by an independent NiCrSi (nickel chrome silicon) thin film heating element beneath a dielectric separating it electrically from the phase change layer. A comparison is made between the predicted and experimental minimum power to amorphize (MPA) for various thermal pulse powers and pulse time lengths. By including both the specific heat and latent heat of fusion for GeTe, we find that the MPA and the minimum power to crystallize follow the form of a hyperbola on the power time effect plot. We also find that the simulated time at which the entire center GeTe layer achieves melting accurately matches the MPA curve for pulse durations ranging from 75–1500 ns and pulse powers from 1.6–4 W.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Respiratory-Induced Haemodynamic Changes: A Contributing Factor to IVC Filter Penetration

PurposeThe purpose of the study is to evaluate the influence of respiratory-induced vena caval hemodynamic changes on filter migration/penetration.Materials and MethodsAfter placement of either a Gunther Tulip or Celect IVC filter, 101 consecutive patients scheduled for filter retrieval were prospectively enrolled in this study. Pre-retrieval CT scans were used to assess filter complications and to calculate cross-sectional area in three locations: at level of filter strut fixation, 3 cm above and 3 cm below. A 3D finite element simulation was constructed on these data and direct IVC pressure was recorded during filter retrieval. Cross-sectional areas and pressures of the vena cava were measured during neutral breathing and in Valsalva maneuver and identified filter complications were recorded. A statistical analysis of these variables was then performed.ResultsDuring Valsalva maneuvers, a 60 % decrease of the IVC cross-sectional area and a fivefold increase in the IVC pressure were identified (p < 0.001). There was a statistically significant difference in the reduction of the cross-sectional area at the filter strut level (p < 0.001) in patient with filter penetration. Difficulty in filter retrieval was higher in penetrated or tilted filters (p < 0.001; p = 0.005). 3D computational models showed significant IVC deformation around the filter during Valsalva maneuver.ConclusionCaval morphology and hemodynamics are clearly affected by Valsalva maneuvers. A physiological reduction of IVC cross-sectional area is associated with higher risk of filter penetration, despite short dwell times. Physiologic data should be used to improve future filter designs to remain safely implanted over longer dwell times.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Design, performance, and cost savings of using GF-PC additively manufactured tooling for stamping of HSS 590 sheet metal

Conventional sheet metal forming tooling in the automotive industry is made up of hardened steel and used for mass-production. Prototype tooling made of metal is durable, but it is only used for a small number of parts despite its high cost, contributing heavily to the vehicle development cost. Additive Manufacturing (AM) offers a low-cost alternative for fabricating tooling suitable for low-volume production. Here, this study investigates the feasibility of using AM polymer composite tooling for the stamping of 1.5mm thickness HSS 590 steel sheets through a two-pronged approach – experimental and numerical analysis. Sheet metal stamping experiments were conducted with glass fiber reinforced polycarbonate (GF-PC) AM polymer tooling and their performance was evaluated based on various metrics such as tool deformation and part accuracy. Finite element simulations of the stamping process were performed, which accurately captured the tool deformations seen in experiments with the use of anisotropic material models. The effective use of simulations in optimizing process parameters to achieve the desired final part geometry is also demonstrated. The GF-PC AM polymer composite tooling was found to be well suited for low volume production applications with potential cost savings.

36 MATERIALS SCIENCE↗

Exemplar for simulation challenges: Large-deformation micromechanics of Sylgard 184/glass microballoon syntactic foams.

Sylgard® 184/Glass Microballoon (GMB) potting material is currently used in many NW systems. Analysts need a macroscale constitutive model that can predict material behavior under complex loading and damage evolution. To address this need, ongoing modeling and experimental efforts have focused on study of damage evolution in these materials. Micromechanical finite element simulations that resolve individual GMB and matrix components promote discovery and better understanding of the material behavior. With these simulations, we can study the role of the GMB volume fraction, time-dependent damage, behavior under confined vs. unconfined compression, and the effects of partial damage. These simulations are challenging and push the boundaries of capability even with the high performance computing tools available at Sandia. We summarize the major challenges and the current state of this modeling effort, as an exemplar of micromechanical modeling needs that can motivate advances in future computing efforts.

36 MATERIALS SCIENCE↗

Effect of proton irradiation on AlGaN/GaN high electron mobility transistor off-state drain breakdown voltage

The effect of proton irradiation on the off-state drain breakdown voltage of AlGaN/GaN high electron mobility transistors (HEMTs) grown on Si substrates was studied by irradiating protons from the backside of the samples through via holes fabricated directly under the active area of the HEMTs. There was no degradation of drain current nor enhancement of off-state drain voltage breakdown voltage observed for HEMTs irradiated with 275 keV protons, for which the defects created by the proton irradiation were intentionally placed in the GaN buffer. HEMTs with defects positioned in the 2 dimensional electron gas channel region and AlGaN barrier using 330 keV protons not only showed degradation of both drain current and extrinsic transconductance but also exhibited an improvement of the off-state drain breakdown voltage. Finite-element simulations showed the enhancement of the latter were due to a reduction in electric field strength at the gate edges by introduction of charged defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

TEAMER: Sandia and CalWave Torque Spring Assessment

Power-take-off (PTO) systems for wave energy converters (WEC) require restoring forces for efficient hydrodynamic as well as mechanical to electric power transfer. Implementation of an effective spring mechanism that can provide a restoring torque on e.g. a rotational, winch type PTO shaft is challenging. SNL in collaboration with CalWave, will use numerical finite element simulations to study application of mechanical torsion springs as a highly efficient restoring force element for CalWave's rotational PTO. Material studies will include common steel as a baseline and potentially advanced materials such as fiber composite materials with different cross-sectional profiles. Test report including description and results of FEA (finite element analysis) models. Excel sheet summarizes which scripts were used for which tables and figures in the report. Part of TEAMER request for technical support (RFTS 1) award to CalWave and Sandia National Lab (SNL).

16 TIDAL AND WAVE POWER↗