Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mesh Generation”

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,099 records · Page 61

Computational Fluid Dynamics Ventilation Study for the Human Powered Centrifuge at the International Space Station

The Human Powered Centrifuge (HPC) is a facility that is planned to be installed on board the International Space Station (ISS) to enable crew exercises under the artificial gravity conditions. The HPC equipment includes a "bicycle" for long-term exercises of a crewmember that provides power for rotation of HPC at a speed of 30 rpm. The crewmember exercising vigorously on the centrifuge generates the amount of carbon dioxide of about two times higher than a crewmember in ordinary conditions. The goal of the study is to analyze the airflow and carbon dioxide distribution within Pressurized Multipurpose Module (PMM) cabin when HPC is operating. A full unsteady formulation is used for airflow and CO2 transport CFD-based modeling with the so-called sliding mesh concept when the HPC equipment with the adjacent Bay 4 cabin volume is considered in the rotating reference frame while the rest of the cabin volume is considered in the stationary reference frame. The rotating part of the computational domain includes also a human body model. Localized effects of carbon dioxide dispersion are examined. Strong influence of the rotating HPC equipment on the CO2 distribution detected is discussed.

Son, Chang H.↗

Finite-time correlations boost large voltage angle fluctuations in electric power grids

Abstract Decarbonization in the energy sector has been accompanied by an increased penetration of new renewable energy sources in electric power systems. Such sources differ from traditional productions in that, first, they induce larger, undispatchable fluctuations in power generation and second, they lack inertia. Recent measurements have indeed reported long, non-Gaussian tails in the distribution of local voltage frequency data. Large frequency deviations may induce grid instabilities, leading in worst-case scenarios to cascading failures and large-scale blackouts. In this article, we investigate how correlated noise disturbances, characterized by the cumulants of their distribution, propagate through meshed, high-voltage power grids. For a single source of fluctuations, we show that long noise correlation times boost non-Gaussian voltage angle fluctuations so that they propagate similarly to Gaussian fluctuations over the entire network. However, they vanish faster, over short distances if the noise fluctuates rapidly. We furthermore demonstrate that a Berry–Esseen theorem leads to the vanishing of non-Gaussianities as the number of uncorrelated noise sources increases. Our predictions are corroborated by numerical simulations on realistic models of power grids.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells are created using polygon-clipping algorithms. The grid is stored in a binary-tree data structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded: A gradient-limited, linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The more robust of a series of viscous flux functions is used to provide the viscous fluxes at the cell interfaces. Adaptively-refined solutions of the Navier-Stokes equations using the Cartesian, cell-based approach are obtained and compared to theory, experiment and other accepted computational results for a series of low and moderate Reynolds number flows.

Coirier, William J.↗

Flexible Microsensor Array for the Monitoring and Control of Plant Growth System

Testing for plant experiments in space has begun to explore active nutrient delivery concepts in which water and nutrients are replenished on a continuous basis for long-term growth. The goal of this study is to develop a novel microsensor array to provide information on the dissolved oxygen environment in the plant root zone for the optimum control of hydroponics and solid substrate plant cultivation systems in the space environment. Miniaturized polarographic dissolved oxygen sensors have been designed and fabricated on a flexible Kapton (trademark) (polyimide) substrate. Two capabilities of the new microsensor array were explored. First, measurements of dissolved oxygen in the plant root zone in hydroponics and solid substrate culture systems were made. The microsensor array was fabricated on a flexible substrate, and then cut out into a mesh type to make a suspended array that could be placed either in a hydroponics system or in a solid substrate cultivation system to measure the oxygen environments. Second, the in situ self-diagnostic and self-calibration capability (two-point for oxygen) was adopted by dynamically controlling the microenvironment in close proximity to the microsensors. With a built-in generating electrode that surrounds the microsensor, two kinds of microenvironments (oxygen-saturated and oxygen-depleted phases) could be established by water electrolysis depending on the polarity of the generating electrode. The unique features of the new microsensor array (small size, multiple sensors, flexibility and self-diagnosis) can have exceptional benefits for the study and optimization of plant cultivation systems in both terrestrial and microgravity environments. The in situ self-diagnostic and self-calibration features of the microsensor array will also enable continuous verification of the operability during entire plant growth cycles. This concept of automated control of a novel chemical monitoring system will minimize crew time required for maintenance, as well as reduce volume, mass, and power consumption by eliminating bulky diagnosis systems including calibrant (fluid and gas) reservoir and flow system hardware.

Kim, Chang-Soo↗

Computational Non-Destructive Evaluation Improving Ultrasonic Interrogation of Complex Geometry Composite Parts

Finite element simulation was employed in modeling the ultrasound (UT) pressure pulse propagation through a coupled liquid-composite medium to reproduce experimental data. From the simulation point of view, the proposed approach is challenging when there is a large simulation domain. For example, it is shown that a sub-micron wavelength of an ultrasonic wave requires a mesh size of several microns and this in turn requires significant computational resources, as well as special care in modeling. Some of the simulation results are presented considering that such modeling should reproduce experimental data for a healthy and faulty composite structure with complex geometry. Many possible experimental setups are simulated to demonstrate the non-destructive testing technique. This setup includes the generation of pressure pulse propagating through the tested composite plate and possible scattering by discontinuities (area of different impedance) that may be present in the panel. This scattered pulse together with the baseline pressure pulse generates a signature on the probe element which can be used to locate the position of defects in the structures.

modeling↗

Generation of spiral bevel gears with conjugate tooth surfaces and tooth contact analysis

A new method for generation of spiral bevel gears is proposed. The main features of this method are as follows: (1) the gear tooth surfaces are conjugated and can transform rotation with zero transmission errors; (2) the tooth bearing contact is localized; (3) the center of the instantaneous contact ellipse moves in a plane that has a fixed orientation; (4) the contact normal performs in the process of meshing a parallel motion; (5) the motion of the contact ellipse provides improved conditions of lubrication; and (6) the gears can be manufactured by use of Gleason's equipment.

Litvin, Faydor L.↗

Generation Method Improves Spiral Bevel Gears

Conjugate tooth surfaces reduce kinematic errors. Method for generation of spiral bevel gears provides conjugated gear-tooth surfaces and improved bearing contact. Conjugated surfaces preferable because maintain constant gear ratios during tooth-contact cycles. Changing gear ratios of nonconjugated surfaces give rise to kinematical errors in transfer of motion from driving gears to driven gears. Errors major sources of noise in power transmissions. Computer program developed to simulate cutting and meshing processes for pinions and gears to minimize kinematical errors.

Litvin, Faydor L.↗

Understanding the impact of the gas diffusion layer structure on catalyst utilization in the PEM water electrolyzer

A multiphysics half-cell model of a polymer electrolyte membrane water electrolyzer (PEMWE) was developed to probe impacts of the detailed 3-dimensional pore structure of the gas diffusion layer (GDL) on performance characteristics. We show that pores in the titanium GDL mesh led to significant underutilization of the catalyst layer (CL), with only 45% of the catalyst effectively utilized. This contradicts the assumption of uniform electron flow across the CL, as shown in graphical abstract (a), as near-zero current was observed near GDL pore regions and the current distribution in CL was influenced by GDL structure, as shown in graphical abstract (b). Instead, oxygen generation was primarily concentrated under the solid titanium regions, diffusing out around the pore walls. High current density peaks were also noted at the GDL-catalyst contact, correlating with degradation hotspots that were directly observed in companion experiments, as shown in graphical abstract (c). Collectively, these findings point to the critical importance of the heterogeneous GDL porous architecture not only for PEMWE efficiency but also for uneven degradation of the CL.

COMSOL↗

Numerical techniques for large cosmological N-body simulations

Techniques for carrying out large N-body simulations of the gravitational evolution of clustering in the fundamental cube of an infinite periodic universe are described and compared. The accuracy of the forces derived from several commonly used particle mesh schemes is examined, showing how submesh resolution can be achieved by including short-range forces between particles by direct summation techniques. The time integration of the equations of motion is discussed, and the accuracy of the codes for various choices of 'time' variable and time step is tested by considering energy conservation as well as by direct analysis of particle trajectories. Methods for generating initial particle positions and velocities corresponding to a growing mode representation of a specified power spectrum of linear density fluctuations are described. The effects of force resolution are studied and different simulation schemes are compared. An algorithm is implemented for generating initial conditions by varying the number of particles, the initial amplitude of density fluctuations, and the initial peculiar velocity field.

Efstathiou, G.↗

Raster-Based Approach to Solar Pressure Modeling

An algorithm has been developed to take advantage of the graphics processing hardware in modern computers to efficiently compute high-fidelity solar pressure forces and torques on spacecraft, taking into account the possibility of self-shading due to the articulation of spacecraft components such as solar arrays. The process is easily extended to compute other results that depend on three-dimensional attitude analysis, such as solar array power generation or free molecular flow drag. The impact of photons upon a spacecraft introduces small forces and moments. The magnitude and direction of the forces depend on the material properties of the spacecraft components being illuminated. The parts of the components being lit depends on the orientation of the craft with respect to the Sun, as well as the gimbal angles for any significant moving external parts (solar arrays, typically). Some components may shield others from the Sun. The purpose of this innovation is to enable high-fidelity computation of solar pressure and power generation effects of illuminated portions of spacecraft, taking self-shading from spacecraft attitude and movable components into account. The key idea in this innovation is to compute results dependent upon complicated geometry by using an image to break the problem into thousands or millions of sub-problems with simple geometry, and then the results from the simpler problems are combined to give high-fidelity results for the full geometry. This process is performed by constructing a 3D model of a spacecraft using an appropriate computer language (OpenGL), and running that model on a modern computer's 3D accelerated video processor. This quickly and accurately generates a view of the model (as shown on a computer screen) that takes rotation and articulation of spacecraft components into account. When this view is interpreted as the spacecraft as seen by the Sun, then only the portions of the craft visible in the view are illuminated. The view as shown on the computer screen is composed of up to millions of pixels. Each of those pixels is associated with a small illuminated area of the spacecraft. For each pixel, it is possible to compute its position, angle (surface normal) from the view direction, and the spacecraft material (and therefore, optical coefficients) associated with that area. With this information, the area associated with each pixel can be modeled as a simple flat plate for calculating solar pressure. The vector sum of these individual flat plate models is a high-fidelity approximation of the solar pressure forces and torques on the whole vehicle. In addition to using optical coefficients associated with each spacecraft material to calculate solar pressure, a power generation coefficient is added for computing solar array power generation from the sum of the illuminated areas. Similarly, other area-based calculations, such as free molecular flow drag, are also enabled. Because the model rendering is separated from other calculations, it is relatively easy to add a new model to explore a new vehicle or mission configuration. Adding a new model is performed by adding OpenGL code, but a future version might read a mesh file exported from a computer-aided design (CAD) system to enable very rapid turnaround for new designs

Wright, Theodore W. II↗

LS-DYNA User-Defined Internal Ballistic Modeling

LS-DYNA Explicit Finite Element Analysis software and CADPROG Internal Ballistic Analysis code are integrated as a single modeling tool to analyze pyro-mechanical devices. Coupling these two codes allows the ballistic calculations to be performed without modeling the propellant in an FEA mesh. At each LS-DYNA time step, CADPROG is run analytically using the kinematic data (stroke, velocity, chamber volume) fed back from the LS-DYNA simulation. The calculated pressure is then applied to a piston or any applicable surface mesh in LS-DYNA model interactively. All external loads such as friction, damping loads, assisting/resisting loads, locking/unlocking loads are modeled in LS-DYNA and coupled with pressure loads from CADPROG to create a fully defined equation of motion Both codes are written in Fortran 77. CADPROG is integrated in the LS-DYNA source code as a user-defined subroutine and recompiled to generate a custom executable.

pyrotechnics↗

Radsource Mr: Mixed Reality Planning Tool For Radioactive Recovery

The RadSource MR system leverages Meta Quest 3's advanced mixed reality capabilities to create a comprehensive spatial planning platform for end-of-life sealed radioactive source recovery operations. The application utilizes the Quest 3's high-resolution passthrough cameras and spatial mapping algorithms to generate accurate 3D environmental models. Core technical components include: (1) Real-time spatial measurement algorithms calculating distances, angles, slopes, and surface areas with sub-centimeter accuracy; (2) Virtual object placement system allowing users to position digital representations of recovery equipment (trailers, containment vessels, protective barriers) within the real environment; (3) Voice recording and annotation system for hands-free documentation in protective equipment; (4) 3D mesh capture and storage capabilities for post-operation analysis and regulatory documentation. (5) Procedure documentation is available for viewing in Mixed Reality, providing an innovative and convenient way to access the information during pre-visit and inspection activities. (6) Support for screen capture for the view for real world and virtual objects together to use it later for planning. The system integrates computer vision techniques for environmental understanding, spatial mathematics for precise measurements, and human-computer interaction principles optimized for hazardous environment operations. Data persistence allows teams to save and share planning sessions across multiple stakeholders while maintaining operational security requirements.

Khadka, Rajiv [Idaho National Laboratory (INL), Id↗

Direct Simulation of Flow over A Wall-Mounted Hump

Simulating separated flows at high Reynolds numbers using Reynolds Averaged Navier-Stokes (RANS) modeled equations remains a challenge in aeronautics. The main hindrance to progress stems from the lack of extended de- tailed data pertinent to the root cause of failure of RANS models, as well as the little progress in RANS modeling innovations in the past several decades. The goal of the current effort is to generate data for separated flow at a Reynolds numbers where conventional models are challenged. We use Direct Numerical Simulations to model turbulent flow over the wall-mounted hump configuration to investigate the physics of flow separation and boundary layer recovery, as well as provide data relevant to the modeling community. A chord-based Reynolds number of Re(sub c) = 47,500 is considered with a turbulent inflow profile of Re(sub θ) = 1,400 (θ /c = 3%). We use FDL3DI, a code that solves the compressible Navier-Stokes equations using high-order compact-difference scheme and filter, with the standard re- cycling/rescaling method of generating turbulent boundary layers as inflow to the computational domain. Two differ- ent configurations of the upper-wall are analyzed for two sets of boundary conditions (slip and no-slip). The results are compared with the available higher Re(sub c) (= 936,000, Re(sub θ) = 7,200, θ /c = 0.77%) experiment for major flow features. The simulated lower Rec allows for DNS-like mesh resolutions, and adequately wide spans. The results from these simulations show earlier separation and delayed reattachment compared to Re(sub c) = 936,000, and significantly higher skin friction in the forebody of the hump. We also find that the upper-wall shape and boundary condition influence pressure distribution over the hump, whereas skin friction is only influenced by the boundary condition.

Wall-Mounted↗

LS-DYNA User-Defined Internal Ballistic Modeling

- LS-DYNA Explicit Finite Element Analysis software and CADPROG Internal Ballistic Analysis code are integrated as a single modeling tool to analyze pyro-mechanical devices. - Combining these two codes allows the ballistic calculations to be performed without modeling the propellant in an FEA mesh. At each LS-DYNA time step, CADPROG is analytically run using the kinematic data (stroke, velocity, chamber volume) fed back from the LS-DYNA simulation. Calculated pressure is then applied to a piston or any applicable surface mesh in LS-DYNA model interactively. -All external loads such as friction, damping loads, assisting/resisting loads, locking/unlocking loads are modeled in LS-DYNA and coupled with pressure loads from CADPROG to create a fully defined equation of motion -Both codes are written in Fortran 77. CADPROG is inserted in LS-DYNA source code as a user-defined subroutine and re-compiled to generate a custom executable.

pyrotechnics↗

Eight cm technology thruster development

A structural integrated ion thruster with 8-cm beam diameter (SIT-8) was developed for attitude control and stationkeeping of synchronous satellites. As optimized, the system demonstrates a thrust T=1.14 mlb (not corrected for beam V sub B = 1200 V (I sub sp = 2200 sec) total propellant utilization efficiency nu sub u = 59.8% (is approximately 72% without auxiliary pulse-igniter electrode), and electrical efficiency n sub E 61.9%. The thruster incorporates a wire-mesh anode and tantalum cover surfaces to control discharge chamber flake formation and employs an auxiliary pulse-igniter electrode for hollow-cathode ignition. When the SIT-8 is integrated with the compatible SIT-5 propellant tankage, the system envelope is 35 cm long by 13 cm flange bolt circle with a mass of 9.8 kg including 6.8 kg of mercury propellant. Two thrust vectoring systems which generate beam deflections in two orthogonal directions were also developed under the program and tested with the 8-cm thruster. One system vectors the beam over + or - 10 degrees by gimbaling of the entire thruster (not including tankage), while the other system vectors the beam over + or - 7 degrees by translating the accel electrode relative to the screen electrode.

Hyman, J., Jr.↗

Three dimensional spline-generated coordinate transformations for grids around wing-body configurations

A direct algebraic method was developed and applied to generate three dimensional grids around wing-body configurations. The method used is a generalized transfinite interpolation method which generates the desired coordinate transformation using geometric data only on the boundaries of the domain of interest. The geometric data that can be specified includes not only coordinates on the boundaries but also out-of-surface parametric derivatives that give a very precise control over the transformation in the vicinity of the surface. In addition to this, the method gives good control over the stretching of the mesh between different boundaries.

Eriksson, L. E.↗

Aeroacoustic Computations of a Generic Low Boom Concept in Landing Configuration: Part 1 - Aerodynamic Simulations

Computational results are presented for a 15%-scale, full-span, generic low-boom concept aircraft model simulated as installed in the NASA Langley Research Center 14- by 22-Foot Subsonic Tunnel and in free-air conditions. The simulations were conducted with the lattice Boltzmann solver PowerFLOW® to capture the time-accurate characteristics of the flow. The aerodynamic behavior of the aircraft was investigated with various combinations of flap and landing gear deployments. To replicate the flow speed used during tests of the model, a Mach number of approximately 0.20 was chosen while varying the angle of attack through a limited range. Analyses were performed on local flow quantities, global forces, as well as flow distributions between configurations. Aerodynamic quantities were shown to be sensitive to mesh resolution levels, especially with flap and landing gear systems retracted. When comparing wind tunnel to free-air simulations, the aerodynamic quantities sampled were found to be in good agreement, suggesting that tunnel blockage and model support effects did not adversely impact the aerodynamic behavior of the model during the wind tunnel test. The results of this study improve our understanding of the highly complex, vortical flow generated by configurations with swept, low aspect-ratio wings at speeds encountered during landing operations.

low boom↗

A Modular Conjugate Heat Transfer Optimization Framework for Thermal Management of Electric Aircraft

Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management, as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.

heat transfer↗