Search NASA⌕ Search

SEARCH · Search NASA

Results for “encounter geometry”

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 253 records · Page 14

Probing Fundamental Mechanisms of Plastic Deformation with High Energy X-rays (Final Progress Report)

Overall, substantial progress was made on developing BCDI techniques in the direction of imaging defects in multi-crystal samples. Section 2.1 describes a new method to reconstruct the 3D atomic displacement field within a grain from multiple coherently resolved diffraction peaks. The key development is the ability to recover a best fit to the grain shape at the same time as the internal displacement (strain) variation with no prior knowledge of the microstructure. This is complemented by the development of a method, Sec. 2.2, to separate out sets of Laue reflections that belong to a single grain (and thus be able to index its orientation) in the commonly encountered situation in which the beam illuminates several grains in a given location. This is essential for enabling collection of coherent diffraction data (at APS 34-ID-C) from neighboring grains. All experiments depend on the quality of samples used. Although good progress has been made with thin film samples, we are working with a broader range of geometries. A new sample design uses a high-Z metal (e.g., Pt) deposited in a sub-micron trench in a low-Z metal (e.g., Al) so that the grain size is prevented from exceeding the current limit for BCDI. We also report in section 3.1 on successfully applying NF-HEDM to a fine-grain Mg sample that would normally be considered out of scope for the method. This represents progress towards being able to map out polycrystalline samples in 3D before performing BCDI on individual grains. We further report in section 3.2 on preliminary NF-HEDM characterization of a sample of UO2 that has sub-micron grains; the data imply the capability of mapping fine-grain polycrystals at 1-ID-E, again for the purposes of BCDI in the bulk. Lastly, our plan is to load grains in polycrystalline samples by inducing thermal stress. We also report the experimental measurements of the Ti- 6Al-4V alloy obtained at the European Synchrotron Radiation Facility (ESRF), which had been indexed and mapped via HEDM at APS in section 9. As discussed below, however, technique development pre-occupied the team throughout the project which meant that we did not make any significant progress towards answering the original hypotheses as stated in the proposal. In summary, we assert that we have made good progress in methods and samples and future work taking advantage of this progress is likely to obtain and analyze BCDI data from multiple adjacent grains in 3D polycrystals. Obtaining such data depends on re-gaining access to the Advanced Photon Source, which is anticipated in late summer or the Fall. In all cases, work has involved close collaboration with faculty, staff, postdocs and students from Brigham Young University (BYU), the Los Alamos National Laboratory and APS beamlines 1-ID and 34-ID.

36 MATERIALS SCIENCE↗

A Design and Fabrication Options Trade Study for Metallic Fuel without Internal Sodium Bonding

Nuclear fuels using alloys of uranium, or metallic fuels, have many beneficial properties. The classical metallic fuel design uses a loose fitting cylindrical “slug” of fuel placed inside stainless-steel cladding tubes where the gap is filled with sodium. This sodium bond is liquid at operating temperature and conducts heat from the slug to the cladding, especially in early life before fuel swells into contact with the cladding. Despite the benefits of sodium bonding, there is a desire to develop metallic fuel technologies without it chiefly to reduce chemical reaction hazards in spent fuel storage from sodium fast reactors operating on once-through fuel cycles. Elimination of the sodium bond may also help unlock potential benefits for fuel fabrication, reactor neutronics, and compatibility with other types of reactors. Creating a sodium-free metallic fuel revolves around the problem of manufacturing fuel slug geometries which are in close contact with the cladding at beginning of life to facilitate heat transport while alleviating fuel-cladding chemical interactions (FCCI) at this interface and providing enough free volume to accommodate fuel swelling. Accelerating development and qualification of this fuel system will require careful selection of design and manufacturing options. To this end, a design trade-off study was performed to evaluate candidate options. Several design and manufacturing options were assessed, weighted, scored, and ranked with respect to fabrication, normal reactor operation, off-normal scenarios, and back-end considerations. This effort was performed both for “baseline” needs, which represented a once-through fuel cycle at temperatures and burnups known to be viable for sodium-bonded metallic fuel, and for “enhanced” needs to represent opportunities for closed fuel cycles and/or more aggressive temperatures/burnups. The outcomes of this study prioritized a baseline technology using U-Zr alloy with additives to mitigate FCCI, produced in annular slug geometry by continuous casting, clad in austenitic stainless-steel alloy, and followed by a final step to swage the cladding down to close the gap. This study prioritized an enhanced fuel technology using U-Mo alloy, also produced by continuous casting into an annular geometry, followed by coating/plating with an FCCI barrier on the slug, again with a final step to swage the cladding diameter down using oxide dispersion strengthened steel. It was noted that development of the enhanced fuel technology would entail more risk, thus U-Zr alloy was put forth as a backup to U-Mo if challenges are encountered with FCCI barriers, and advanced ferritic/martensitic steels are put forth as a backup to oxide dispersion strengthened steels if swaging and welding are found unworkable.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Spacecraft Charging Issues for Launch Vehicles

Spacecraft charging is well known threat to successful long term spacecraft operations and instrument reliability in orbits that spend significant time in hot electron environments. In recent years, spacecraft charging has increasingly been recognized as a potentially significant engineering issue for launch vehicles used to deploy spacecraft using (a) low Earth orbit (LEO), high inclination flight trajectories that pass through the auroral zone, (b) geostationary transfer orbits that require exposures to the hot electron environments in the Earths outer radiation belts, and (c) LEO escape trajectories using multiple phasing orbits through the Earths radiation belts while raising apogee towards a final Earth escape geometry. Charging becomes an issue when significant areas of exposed insulating materials or ungrounded conductors are used in the launch vehicle design or the payload is designed for use in a benign charging region beyond the Earths magnetosphere but must survive passage through the strong charging regimes of the Earths radiation belts. This presentation will first outline the charging risks encountered on typical launch trajectories used to deploy spacecraft into Earth orbit and Earth escape trajectories. We then describe the process used by NASAs Launch Services Program to evaluate when surface and internal charging is a potential risk to a NASA mission. Finally, we describe the options for mitigating charging risks including modification of the launch vehicle andor payload design and controlling the risk through operational launch constraints to avoid significant charging environments.

Electrostatics↗

Spacecraft Charging Issues for Launch Vehicles

Spacecraft charging is well known threat to successful long term spacecraft operations and instrument reliability in orbits that spend significant time in hot electron environments. In recent years, spacecraft charging has increasingly been recognized as a potentially significant engineering issue for launch vehicles used to deploy spacecraft using (a) low Earth orbit (LEO), high inclination flight trajectories that pass through the auroral zone, (b) geostationary transfer orbits that require exposures to the hot electron environments in the Earths outer radiation belts, and (c) LEO escape trajectories using multiple phasing orbits through the Earths radiation belts while raising apogee towards a final Earth escape geometry. Charging becomes an issue when significant areas of exposed insulating materials or ungrounded conductors are used in the launch vehicle design or the payload is designed for use in a benign charging region beyond the Earths magnetosphere but must survive passage through the strong charging regimes of the Earths radiation belts. This presentation will first outline the charging risks encountered on typical launch trajectories used to deploy spacecraft into Earth orbit and Earth escape trajectories. We then describe the process used by NASAs Launch Services Program to evaluate when surface and internal charging is a potential risk to a NASA mission. Finally, we describe the options for mitigating charging risks including modification of the launch vehicle and/or payload design and controlling the risk through operational launch constraints to avoid significant charging environments

electrostatics↗

Sensitivity Analysis of Salt Storage Cavern Mechanical Integrity Test Parameters

This report examines the sensitivity of salt cavern Mechanical Integrity Tests (MIT) to uncertainties in key test parameters. MIT's are used by cavern operators to detect and quantify leak rates in access wells in underground salt storage caverns and involve a suite of measured and assumed parameters that have a direct impact on the sensitivity of the testing to detect actual leaks from the cavern storage system. Determining the sensitivity of the testing to these different parameters provides a basis for understanding the results from, and informing the design criteria for, this type of testing. Without fully understanding the sensitivity of the test to the testing parameters, it is possible that the test results may not accurately reflect the integrity of the cavern system; an actual leak may be missed, or an intact system may be interpreted as leaking. This report reviews the main parameters included in MITs and examines how selected changes in their values can impact test results. The deviations used in the sensitivity analyses were designed to be within the ranges believed to be similar to those which may be encountered during testing. The results show that small, plausible fluctuations in some of the parameters measured values can have a significant impact on the testing results. Of the parameters studied here, the sensitivity analyses showed the order of importance to be (from highest to lowest): nitrogen-oil interface depth measurement, well bore temperature, well head pressure, and finally the internal geometry of the testing interval.

42 ENGINEERING↗

Transonic airfoil and wing design using Navier-Stokes codes

An iterative design method has been implemented into 2D and 3D Navier-Stokes codes for the design of airfoils or wings with given target pressure distributions. The method begins with the analysis of an initial geometry, and obtains the analysis pressure distributions of that geometry. The differences between analysis pressures and target pressures are used to drive geometry changes through the use of a streamline curvature method. This paper describes the procedure that makes the iterative design method work for Navier-Stokes codes. Examples of 2D airfoil design, and 3D wing design are included. It is demonstrated that the method is highly effective for airfoil or wing design at flow conditions where no substantial separation occurs. Problems encountered in the airfoil design with shock induced flow separations are discussed.

Yu, N. J.↗

First Galileo image of asteroid 243 Ida

The second spacecraft encounter with an asteroid has yielded an unprecedentedly high resolution portrait of 243 Ida. On 28 Aug. 1993, Galileo obtained an extensive data set on this small member of the Koronis family. Most of the data recorded on the tape recorder will be returned to Earth in spring 1994. A five-frame mosaic of Ida was acquired with good illumination geometry a few minutes before closest approach; it has a resolution of 31 to 38 m/pixel amd was played back during Sept. 1993. Preliminary analyses of this single view of Ida are summarized.

Chapman, C. R.↗

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↗

Theoretical evaluation of a V/STOL fighter model utilizing the PAN AIR code

The PAN AIR computer code was investigated as a tool for predicting closely coupled aerodynamic and propulsive flowfields of arbitrary configurations. The NASA/Ames V/STOL fighter model, a configuration of complex geometry, was analyzed with the PAN AIR code. A successful solution for this configuration was obtained when the nozzle exit was treated as an impermeable surface and no wakes were included around the nozzle exit. When separated flow was simulated from the end of the nacelle, requiring the use of wake networks emanating from the nozzle exit, a number of problems were encountered. A circular body nacelle model was used to investigate various techniques for simulating the exhaust plume in PAN AIR. Several approaches were tested and eliminated because they could not correctly simulate the interference effects. Only one plume modeling technique gave good results. A PAN AIR computation that used a plume shape and inflow velocities obtained from the Navier-Stokes solution for the plume produced results for the effects of power that compared well with experimental data.

Howell, G. A.↗

Modeling of Wake-vortex Aircraft Encounters

There are more people passing through the world's airports today than at any other time in history. With this increase in civil transport, airports are becoming capacity limited. In order to increase capacity and thus meet the demands of the flying public, the number of runways and number of flights per runway must be increased. In response to the demand, the National Aeronautics and Space Administration (NASA), in conjunction with the Federal Aviation Administration (FAA), airport operators, and the airline industry are taking steps to increase airport capacity without jeopardizing safety. Increasing the production per runway increases the likelihood that an aircraft will encounter the trailing wake-vortex of another aircraft. The hazard of a wake-vortex encounter is that heavy load aircraft can produce high intensity wake turbulence, through the development of its wing-tip vortices. A smaller aircraft following in the wake of the heavy load aircraft will experience redistribution of its aerodynamic load. This creates a safety hazard for the smaller aircraft. Understanding this load redistribution is of great importance, particularly during landing and take-off. In this research wake-vortex effects on an encountering 10% scale model of the B737-100 aircraft are modeled using both strip theory and vortex-lattice modeling methods. The models are then compared to wind tunnel data that was taken in the 30ft x 60ft wind tunnel at NASA Langley Research Center (LaRC). Comparisons are made to determine if the models will have acceptable accuracy when parts of the geometry are removed, such as the horizontal stabilizer and the vertical tail. A sensitivity analysis was also performed to observe how accurately the models could match the experimental data if there was a 10% error in the circulation strength. It was determined that both models show accurate results when the wing, horizontal stabilizer, and vertical tail were a part of the geometry. When the horizontal stabilizer and vertical tail were removed there were difficulties modeling the sideforce coefficient and pitching moment. With the removal of only the vertical tail unacceptable errors occurred when modeling the sideforce coefficient and yawing moment. Lift could not be modeled with either the full geometry or the reduced geometry attempts.

Smith, Sonya T.↗

NRC Multiphysics Analysis Capability Deployment FY21: Part 3

This report details the progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Support for Advanced Non-Light Water Reactors.” The deliverables completed for this report are: Deliverable 1c: the capability to model gas mixtures was added to Pronghorn. A test problem mimicking the conditions achieved in a depressurized loss of forced cooling (DLOFC) event was solved with both RELAP-5 and Pronghorn. Pronghorn employed a finite vol ume method with the Kurganov-Tadmor discretization. The comparison between the mass fraction spatial profiles computed with RELAP-5 and Pronghorn clearly shows the presence of numerical artifacts (i.e., overly diffusive behavior at low Mach numbers). We confirmed that the problem disappears at higher Mach numbers. We recommend future work on the implementation of a low Mach finite volume formulation to better treat low Mach number problems. Deliverable 2a: we demonstrated two approaches to model the radiation/conduction/natural convection heat transfer across a stagnant gas for the PBMR-400 design using Pronghorn. The first approach is based on the net radiation method, which relies on the computation of view factors with the Multiphysics Object-Oriented Simulation Environment (MOOSE) ray tracing capability. The second method is a traditional thermal resistance approach. The test problems include both 2D and 3D geometries. In all cases, the results show very good agreement during a DLOFC transient. This confirms that the faster thermal resistance method produces solutions that are equivalent to the net radiation method for this geometry. Deliverable 3d: we demonstrated the use of the advection kernel for the delayed neutron precursor equation in Griffin with a 2D MSFR model. The results appear physical but further verification is recommended. We also recommend the addition of conjugate heat transfer to compute the temperatures and model the thermomechanic behavior of the reflectors and other structures. Significant memory and performance issues were encountered in the 3D axisymmetric model. Future work is recommended in this area. Task 8g: this task allows multidimensional MOOSE applications to be coupled to system codes (RELAP-7 and SAM). We implemented a faster multiphysics iteration coupling algorithm, which provides an overall 6× acceleration of the 3D-1D coupling of the core multidi- mensional fluid flow solver and the 1D primary and secondary loop model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Slat Noise Control Using a Slat Gap Filler

The leading edge slat of a high-lift system is one of the main noise contributors on many commercial aircraft during approach. This paper continues our previous studies on the gap filler for passive noise control on the 30P30N high-lift airfoil. An improved implementation of the gap filler was applied to minimize the effects of flow leakage encountered in the previous work, which resulted in spurious noise content in the far-field acoustic spectra. To evaluate the effect of passive flow control on the acoustics generated by the unsteady flow field, anechoic wind tunnel experiments are conducted on the two-dimensional, three-element high-lift airfoil with a gap filler mounted to the slat. The slat geometry modification associated with the gap filler alters the flow field in the cove region that dominates the generation of the acoustic field. A single angle of attack (훼푘=8◦) and three flow speeds corresponding to Reynolds numbersof푅푒푐=1.2푒6,1.5푒6, and1.71푒6are selected as the test conditions. Steady surface pressure measurements are conducted to assess the effect of the treatments on the overall lift. Acoustic array measurements are used to evaluate the influence of the gap filler on the radiated noise. Delay and Sum beamforming is applied to locate the noise sources on the model and to provide the integrated spectra. The gap filler is found to eliminate the narrowband peaks in the acoustic spectra and, also, to yield a 10 dB reduction in the broadband noise in comparison with the baseline case. Time-resolved Particle Image Velocimetry results show that the flow features are significantly altered with the presence of the gap filler, which leads to a more stable slat cove shear layer and, thus, to weaker pressure and velocity fluctuations.

Yang Zhang↗

Reactive Hydrodynamics in Rotating Spherical and Cylindrical Geometry

In turbulent spray combustion among many complex interactions between local flow structures called turbulent eddies and droplets are those associated with rotation of droplets. In general, for a complete statistical description of turbulent sprays, consideration of at least four degrees of freedom respectively associated with translational, rotational, vibrational (pulsational), and internal motions of the droplet are needed. Clearly the interactions between all degrees of freedom of the droplets and those for the gaseous background field will be exceedingly complex. For example, one type of interaction between the translational and the rotational velocity of droplets results in droplet helicity, H(d) = w(d).v(d), the significance of which in turbulent spray combustion is yet to be recognized. The role of droplet rotation in turbulent spray combustion modeling and its impact on the evaporation of liquid fuel droplets was recently investigated. Also, the impact of rotation on combustion of solid particles such as is encountered in pulverized coal combustion has been emphasized. The problem of viscous flow around a rotating sphere discussed above also occurs in other areas of physical sciences such as astrophysics and geophysics. Consequently, the subject has been addressed in many classical as well as more recent investigations. According to these investigations, the rotation of a rigid sphere in an otherwise quiescent, unconfined environment results in the motion of the fluid towards the poles. The polar flows from the northern and southern hemispheres move along helical trajectories towards the equatorial plane. Eventually, the polar flows collide at the equatorial plane, thus producing a sheet of rotating fluid that is radially ejected outward on this plane. Therefore, a droplet induces a strained flow field as a result of its rotation. Since the spatial extent of equatorial jets could easily exceed many droplet diameters, interactions between neighboring droplets are enhanced as a result of their rotation. Also, the equatorial jet substantially alters the spherical geometry of the diffusion flame surface that surrounds a rotating droplet. The objective of the research is to gain more knowledge about the hydrodynamics within and around rotating spherical and cylindrical body of fluid, and the behavior of diffusion or premixed flame surfaces that could surround such symmetric body of rotating fluid.

Sohrab, Siavash H.↗

The quality control programme for ITk strip tracker module assembly

The assembly of the ATLAS Inner Tracker requires the construction of 19,000 silicon strip sensor detector modules in eight different geometries. Modules will be assembled and tested at 31 institutes on four continents from sensors, readout chips, and flexes. In order to adhere to the module specifications defined for sufficient tracking performance, a rigorous programme of quality control (QC) was established to cover components at every stage of assembly. Here, this contribution presents an overview of the QC programme for ITk strip tracker modules, issues encountered during the pre-production phase (5% of the production volume), and their solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Combustion instability analysis for liquid propellant rocket engines

The multi-dimensional numerical model has been developed to analyze the nonlinear combustion instabilities in liquid-fueled engines. The present pressure-based approach can handle the implicit pressure-velocity coupling in a non-iterative way. The additional scalar conservation equations for the chemical species, the energy, and the turbulent transport quantities can be handled by the same predictor-corrector sequences. This method is time-accurate and it can be applicable to the all-speed, transient, multi-phase, and reacting flows. Special emphasis is given to the acoustic/vaporization interaction which may act as the crucial rate-controlling mechanism in the liquid-fueled rocket engines. The subcritical vaporization is modeled to account for the effects of variable thermophysical properties, non-unitary Lewis number in the gas-film, the Stefan flow effect, and the effect of transient liquid heating. The test cases include the one-dimenisonal fast transient non-reacting and reacting flows, and the multi-dimensional combustion instabilities encountered in the liquid-fueled rocket thrust chamber. The present numerical model successfully demonstrated the capability to simulate the fast transient spray-combusting flows in terms of the limiting-cycle amplitude phenomena, correspondence between combustion and acoustics, and the steep-fronted wave and flame propagation. The investigated parameters include the spray initial conditions, air-fuel mixture ratios, and the engine geometry. Stable and unstable operating conditions are found for the liquid-fueled combustors. Under certain conditions, the limiting cycle behavior of the combusting flowfields is obtained. The numerical results indicate that the spray vaporization processes play an important role in releasing thermal energy and driving the combustion instability.

Kim, Y. M.↗

Operating condition and geometry effects on low-frequency afterburner combustion instability in a turbofan at altitude

Three afterburner configurations were tested in a low-bypass-ratio turbofan engine to determine the effect of various fuel distributions, inlet conditions, flameholder geometry, and fuel injection location on combustion instability. Tests were conducted at simulated flight conditions of Mach 0.75 and 1.3 at altitudes from 11,580 to 14,020 m (38,000 to 46,000 ft). In these tests combustion instability with frequency from 28 to 90 Hz and peak-to-peak pressure amplitude up to 46.5 percent of the afterburner inlet total pressure level was encountered. Combustion instability was suppressed in these tests by varying the fuel distribution in the afterburner.

Cullom, R. R.↗

Some mechanical design aspects of the European Robotic Arm

The European Robotic Arm (ERA) is a contribution to the Russian Segment of the International Space Station Alpha. It will start operating on the Russian Segment during the assembly phase. ERA is designed and produced by a large industrial consortium spread over Europe with Fokker Space & Systems as prime contractor. In this paper, we will describe some of the overall design aspects and focus on the development of several mechanisms within ERA. The operation of ERA during the approach of its end effector towards the grapple interface and the grapple operation is discussed, with a focus on mechanisms. This includes the geometry of the end effector leading edge, which is carefully designed to provide the correct and complete tactile information to a torque-force sensor (TFS). The data from this TFS are used to steer the arm such that forces and moments are kept below 20 N and 20 N.m respectively during the grappling operation. Two hardware models of the end effector are built. The problems encountered are described as well as their solutions. The joints in the wrists and the elbow initially used a harmonic drive lubricated by MoS2. During development testing, this combination showed an insufficient lifetime in air to survive the acceptance test program. The switch-over to a system comprising planetary gearboxes with grease lubrication is described. From these development efforts, conclusions are drawn and recommendations are given for the design of complex space mechanisms.

Lambooy, Peter J.↗

Derivation of Effective Properties Based on Porous Scale Simulations Using Filtering Techniques

This study presents a method for derivation of effective properties at the interface and in-depth of porous materials. The method defines a Representative Elementary Volume (REV) and applies filtering techniques to computer effective properties such as porosity and flow quantities, such as velocity and pressure. The script, developed to process the data was tested on the VTK type files that contain the mesh information and the flow solution. The method allows to choose between two types of filters, such as cellular and top-hat and define the size of the REV and number of samples along the domain. Extraction of the REV from the domain is performed to exact boundaries requested for the user. This is done using a triangulation technique and cutting through the cells to comply to the requested boundaries of the volume. The method can be applied to both structured and unstructured meshes. Filtering the material porosity and flow quantities involves integration of the numerical data. The algorithm provides three integration methods, such as Riemann sum, Monte Carlo and Quadrature rule to perform the integration. The Monte-Carlo technique permits the use of either uniform or linearly spaced distribution of points. The Quadrature rule is currently applicable to tetrahedral element types. The Monte Carlo and Quadrature rule methods require interpolation of the flow quantities at the sample points. For interpolation, two methods were tested and are readily available, Gaussian interpolation and re-sampling. It has been shown that re-sampling method has better consistency and acceptable accuracy in interpolation of the data. The algorithm was written in Python language and uses a number of modules. The major module besides numpy is PyVista. It is used to process the computational domain, clip the REV and interpolate the data. Quadrature rule integration was performed using a quadpy module. ParaView software was used externally to convert the flow solution to the VTK (or more specifically VTU) format. Integration of ParaView in the same environment with PyVista encountered problems and could not be implemented in this work. The developed algorithm is expected to be applicable to unstructured meshes and more complex porous structures as soon as the data can be passed in VTK type format. With the report is provided Python script for filtering the solution and a Matlab script for simple generation and processing of 2-D and 3-D porous channel geometries. The two scripts don't communicate.

Alexsander Zibitsker↗