Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nuclear science and engineering - General studies of nuclear reactors”

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.

3,432 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Anomalous lattice expansion of bcc Ta in helium-loaded diamond anvil cells

Anomalous lattice expansion of bcc Ta is observed in helium-loaded diamond anvil cell (DAC) experiments near 6 GPa and above 400 K, while co-loaded W, Mo, and Au exhibit normal thermal expansion. The Ta volume expansion increases with temperature, reaching ~14% near 650 K before plateauing, and is largely retained after recovery to ambient conditions. Time-resolved measurements at 500 K reveal progressive expansion over several hours, consistent with gradual, thermally activated He incorporation. The largely retained expansion during subsequent cooling and holding at 400 K indicates that the incorporated He is predominantly kinetically trapped. The data reveal two distinct regimes: below ~550 K, modest and spatially uniform expansion consistent with interstitial He trapping at isolated, randomly distributed vacancies; above ~550 K, He-vacancy cluster growth producing dramatically larger expansion, peak broadening, and a core–shell microstructure revealed by micron-resolution X-ray diffraction mapping. These results demonstrate that He is not always a passive pressure medium in DAC experiments and reveal a previously unrecognized regime of He-mediated lattice modification in Ta, complementing conventional ion-implantation studies relevant to nuclear and fusion structural materials.

Diamond anvil cell↗
DOE OSTI2026

Impacts of Multidimensional Progenitor Perturbations on Core-collapse Supernova Explosions

Numerical studies of core-collapse supernovae have demonstrated the importance of nonradial motions in precollapse progenitors on the explosion outcome. We use the Chimera neutrino radiation hydrodynamics code running seven two-dimensional simulations of 15 M⊙ progenitors with different progenitor structures introduced by different one- and two-dimensional precollapse stellar evolution environments to examine the impacts of stellar structure and nonspherical motion in the precollapse progenitor on the development of explosions. We compare the explosion evolution of these models in terms of shock dynamics, diagnostic energy, neutrino heating, accretion, explosion geometry, nuclear abundances, and turbulent convection. We also analyze how stochastic variation impacts our simulations. Contrary to results reported in prior studies examining the impacts of multidimensional progenitors, we observe similar shock revival times and explosion development in our simulations despite differences in initial compositions and structures. We find no discernible impact from the accretion of nonradial perturbations from a multi-D progenitor onto the stalled shock in the revival and strength of explosion, as fully developed neutrino-driven convection behind the stalled shock is similar for all our models. For models with physically sourced noise in the iron core, a strong oscillation of the shock occurs after bounce and deflects infall laterally, and accelerates the saturation of the lateral turbulent kinetic energy. An examination of model stochasticity shows that any prior expected impacts on explosive outcome due to convection-related perturbations lie below the detectable threshold of numerical variation.

Chen, Chien-Hui [North Carolina State University]↗
NASA NTRS2026

Understanding and Controlling Hybrid Electric Gas Turbine Engine Transient Dynamics

The electrification of the gas turbine engine is known to increase the flexibility of aircraft architectures by enabling the generation of electrical power to distribute to other electrically based, thrust producing subsystems. It also has potential for direct performance benefits in the gas turbine engine itself, both at steady state and dynamically. Although the design focus of the gas turbine engine performance is primarily at steady state, it is often the instabilities occurring during transients that cause disequilibrium and constrain performance improvements. Instabilities arise due to the disequilibrium of the energy storage mechanisms within the traditional engine system and likewise for the electrified engine system. The primary energy storage mechanisms in the traditional system are the rotational inertia, gas path volumes, the thermal masses that make up the mechanical structure, and now, the electrical power system will provide additional contributions. Understanding the effect each of these energy storage mechanisms has on the others and controlling them appropriately allows for the suppression of state changes within the turbomachinery components to the degree that the components remain near steady state, thus reducing the disequilibrium within the system. The ability to tightly regulate the state changes of turbomachinery components, such as the compressor, minimizes the excursion of the compressor operating point from the operating line (operability), allowing for higher performing, more efficient compressor designs by decreasing the amount of stall margin needed for safe engine power level changes. Preliminary studies with the electrification of the turbine engine have shown that this is possible, and this concept can lead to design trades benefitting engine performance, weight, and volume. This paper explores in detail the energy storage mechanisms and control approaches for coordinating their state changes, and ultimately proposes that a higher performing, more efficient compressor design can result.

transient dynamics↗
DOE OSTI2025

Glancing Angle Deposition in Gas Sensing: Bridging Morphological Innovations and Sensor Performances

Glancing Angle Deposition (GLAD) has emerged as a versatile and powerful nanofabrication technique for developing next-generation gas sensors by enabling precise control over nanostructure geometry, porosity, and material composition. Through dynamic substrate tilting and rotation, GLAD facilitates the fabrication of highly porous, anisotropic nanostructures, such as aligned, tilted, zigzag, helical, and multilayered nanorods, with tunable surface area and diffusion pathways optimized for gas detection. This review provides a comprehensive synthesis of recent advances in GLAD-based gas sensor design, focusing on how structural engineering and material integration converge to enhance sensor performance. Key materials strategies include the construction of heterojunctions and core–shell architectures, controlled doping, and nanoparticle decoration using noble metals or metal oxides to amplify charge transfer, catalytic activity, and redox responsiveness. GLAD-fabricated nanostructures have been effectively deployed across multiple gas sensing modalities, including resistive, capacitive, piezoelectric, and optical platforms, where their high aspect ratios, tailored porosity, and defect-rich surfaces facilitate enhanced gas adsorption kinetics and efficient signal transduction. These devices exhibit high sensitivity and selectivity toward a range of analytes, including NO2, CO, H2S, and volatile organic compounds (VOCs), with detection limits often reaching the parts-per-billion level. Emerging innovations, such as photo-assisted sensing and integration with artificial intelligence for data analysis and pattern recognition, further extend the capabilities of GLAD-based systems for multifunctional, real-time, and adaptive sensing. Finally, current challenges and future research directions are discussed, emphasizing the promise of GLAD as a scalable platform for next-generation gas sensing technologies.

Chemistry↗
NASA NTRSDate not supplied

Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)

Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.

Bohdan O Wesely↗
DOE OSTI2026

Large Language Model for Validation, Optical Calibration, and Learning (VOCAL) Distributed Temperature Sensing Interface

Distributed temperature sensing (DTS) using fiber optic sensors (FOS) offers a promising method for temperature measurements in advanced reactors, such as sodium fast reactors and molten salt cooled reactors. To support the calibration and validation of DTS measurements, Argonne National Laboratory developed the Validation, Optical Calibration, and Learning (VOCAL) software package. This report describes the integration of a local large language model (LLM) with a retrieval-augmented generation (RAG) system into the VOCAL interface to serve as an interactive user assistant. The LLM framework enhances the VOCAL platform’s accessibility to users by explaining interface components, clarifying inputs and outputs, and answering user queries dynamically in real-time. The accuracy of the LLM assistant performance was evaluated with 20 queries regarding the interface and its parameters using experimental data from the Thermal Hydraulic Experimental Test Article (THETA) facility. Results demonstrate that the LLM achieved a 95% accuracy rate, with a BERTScore of 0.8816 and SBERT value of 0.7417. Furthermore, validation of the RAG system within the LLM framework showed optimal accuracy with k-values between 1 and 2 using the k-refinement convergence test. The prompt perturbation analysis demonstrated good initial consistency for the RAG system, exhibiting the highest accuracy under punctuation variations and the greatest sensitivity under query reordering. Notably, the model’s errors were limited to data retrieval failures rather than factual hallucinations, reinforcing its baseline reliability. The integration of LLM provides a highly accurate, userfriendly enhancement to the VOCAL platform without disrupting its core computational capabilities for FOS calibration and validation.

Hong, Evan↗
NASA NTRS2026

Multi-Agent Swarm State of the Art Report

The Next-Generation Multi-Agent Swarm (NGS) Study conducted by NASA’s Ames Research Center for NASA’s Space Technology Mission Directorate (STMD) will develop a comprehensive understanding of emerging multi-agent swarm capabilities. The study aims to identify existing swarm capabilities and asses their potential for persistent lunar space situational awareness, surface monitoring, and distributed autonomy demonstrations. A key objective is to inform the design of a next-generation multi-agent swarm that can perform autonomous distributed remote sensing, position, navigation, and timing (PNT) services, automated deployment that leverages autonomy, edge computing, and interoperable networking to enable cooperative operations without the need for immediate human operation. This study will address specific shortfalls identified by STMD, including intelligent multi-agent constellations, autonomy, edge computation, position, navigation, and timing for small spacecraft, small spacecraft propulsion, and space situational awareness (1625, 1438, 1433, 1557, 1431, 1430, 1589). The NASA Ames Mission Design Center (MDC) will provide subject matter expertise to support systems engineering trades, while experts in autonomy and spacecraft swarms in NASA’s Intelligent Systems Division will lead the study and focus on identifying emerging next-generation swarm capabilities. The study objectives include: capturing the current state-of-the-art for multi-agent swarm capabilities, evaluating technologies and creating technology roadmaps, and developing at least one new technology demonstration mission concept. This initial NGS study report surveys the current state of the art in technology areas relevant for the next-generation multi-agent swarm design. Our primary focus is on surveying relevant deployed space systems1, supplemented with selective analysis of relevant proposed missions and technology developments that have yet to fly.

agent↗
DOE OSTI2026

Dressed-State Hamiltonian Engineering in a Strongly Interacting Solid-State Spin Ensemble

In quantum science applications, ranging from many-body physics to quantum metrology, dipolar interactions in spin ensembles are often controlled via Floquet engineering. However, this technique typically reduces the interaction strength between spins and effectively weakens the coupling to a target sensing field, limiting the metrological sensitivity. In this Letter, we develop and demonstrate an alternative method that directly tunes the native dipolar interaction in an ensemble of nitrogen-vacancy (NV) centers in diamond, thereby overcoming these limitations inherent to Floquet engineering. Our approach utilizes dressed-state qubit encoding under a bias magnetic field applied perpendicular to the crystal lattice orientation. This method leads to a 3.2× enhancement of the dimensionless coherence parameter JT 2 compared to state-of-the-art Floquet engineering and a 2.6× (8.3 dB) enhanced sensitivity in ac magnetometry. Furthermore, our results provide a powerful Hamiltonian engineering tool for future studies with NV ensembles and other interacting higher-spin (S > $\frac{1}{2}$) systems.

Quantum control↗
DOE OSTI2026

Software Validation Work With The ZPPR-15 Data

The analysis activities for fast reactors involve using many different pieces of software that are relied upon for their predictive capabilities. For this software to be considered reliable, documented proof that the predictions of the software are accurate is required. In this manuscript, the validation work that covers some of the Argonne software used in fast reactor design activities is discussed and displayed. This validation work includes neutron and gamma flux distributions, reaction rate distributions, and reactivity worth. In an ideal world, a reactor development program would have access to a comprehensive set of experimental facilities to help inform the design aspects of the reactor itself. While thermal-hydraulics experiments, and to a limited degree mechanical experiments, can be carried out today for validation needs, neutronics related experimental facilities are rather impractical because of the lack of experimental facilities. Given the desired time table for construction of new reactors, the reconstitution or creation of new neutronic experimental facilities is untenable and thus those reactor development programs must rely upon any available experimental measurements that are qualitatively similar to the design. While a methodology has been proposed to assess the similarity between the past experimental measurements and the reactor itself, that aspect is beyond the scope of this manuscript. In this manuscript, the focus is entirely placed on the analysis results for a series of experiments carried out at the ZPPR facility in Idaho in the mid-1980s. In this regard, this manuscript only shows the validation of the stated neutronics software for specific loadings of the ZPPR reactor. Because of the fuel form, its proposed enrichment, and the material content of the reactor core, the ZPPR-15 experiments were identified as potential validation data for the reactor. The ZPPR-15 experiments were intended as mockups of a 330 MWe Integral Fast Reactor program which was a follow on program to the Clinch River Breeder Reactor. In the ZPPR-15 series of experiments, measurements of the neutron spectrum, control rod worth, sodium void worth, foil reaction rate distributions, Doppler worth of heated samples, gamma dose, and axial expansion worth were all carried out and published. In many cases, these reactivity coefficients are good candidates to validate the reactivity coefficient calculation scheme used by the analysis software and included in the safety analysis activities of fast reactor development projects today. This manuscript discusses the modeling methodology and accuracy of the calculated experimental results using the LANL software MCNP and the ANL software package ARC (Argonne Reactor Codes). As will be shown, for many of the experimental measurements, the two software packages are found to be good predictive analysis tools for those experiments. In other cases, problems with the analysis methodology or underlying cross section data are exposed which indicates where predictive analysis is not as reliable. Finally, in some of the measurements the conclusion is reached that the experimental measurement cannot be reproduced with the analysis software as it is simply too difficult.

Aliberti, Gerardo↗
DOE OSTI2026

Equilibrium and non-equilibrium effects in high pressure phase transformations of carbon

The behavior of carbon in the range 1–100 GPa and 1–10 kK is central to problems in planetary interiors, inertial confinement fusion targets, and high-pressure synthesis of carbon-based materials, but experiments in this regime are difficult and often provide only indirect constraints on phase behavior. As a result, phase boundary loci, structure, and limits of metastability at high pressure remain uncertain. In this work, machine-learning enhanced atomistic simulations are used to address this knowledge gap. We determine the melt line up to 100 GPa, the graphite-diamond phase boundary up to the melt line, and analyze structure of the coexisting phases. We show that the coexisting liquid evolves smoothly with pressure without evidence for a first-order liquid–liquid transition. Orientation-resolved graphite melting simulations indicate that basal-plane interfaces develop a dewetting layer and undergo layer-by-layer melting, producing kinetic hysteresis and an apparent orientation dependence of the melt line. Non-equilibrium quenches from the melt are used to construct a kinetically limiting graphite–diamond phase boundary for rapid quenches from above the melt line, and show that graphite is metastable at pressures of up to ≈ 25 GPa. These results provide bounds on equilibrium and metastable behavior in carbon relevant for interpreting high-pressure experiments and for designing synthesis pathways to specific carbon microstructures.

Lyu, Yanjun [Department of Materials Science and E↗
DOE OSTI2026

Single particle – MC-ICP-MS for isotopic analysis of uranium particles

Single particle – multi-collector – inductively coupled plasma – mass spectrometry (SP-MC-ICP-MS) was employed to measure a suspension of 1 µm U3O8 particles (∼1.3 pg total U/particle) to determine their individual isotopic compositions of 234U/238U, 235U/238U, and 236U/238U. The effects of different detector combinations for 235U and 238U, including secondary electron multipliers (SEM) and Faraday detectors (1011 and 1013 Ω amplifiers), were explored for accuracy and precision optimization on the observed 235U/238U. The minor isotopic ratios (i.e., 234U/238U and 236U/238U) were analyzed such that the 234U and 236U were monitored on SEM detectors and the 238U was monitored on a Faraday (1011 Ω) detector. Various integration times (5, 10, 25, and 50 ms) were investigated in all detector configurations to gain a better understanding of their impact on sensitivity, accuracy, and precision. For 235U/238U ratios of 1 µm U3O8 particles, a dual Faraday detector measurement with 1011 Ω was the optimal choice; measurement of 1021 particles yielded an average 235U/238U ratio of 0.00170 (14), a −1.8% relative difference (% RD) from the reference value. The minor isotopic compositions were determined to be 0.0000070 (14) and 0.0000758 (48) for the 234U/238U and 236U/238U, respectively. These measurements correspond to <8% and <1% RD from their reference value for the 234U/238U and 236U/238U, respectively. SP-MC-ICP-MS was also able to provide insight into measurement sensitivity. In these individual particles, merely 15 and 165 atto-grams (ag) of 234U and 236U were present (calculated). Initial limits of detection for SP-MC-ICP-MS were determined to be ∼1.0 ag (when measured via SEM detectors). This valuable approach is applicable to areas including nuclear forensics, nuclear safeguards, and geochemical analysis, which require high-precision measurements of uranium within micron-sized particles.

Manard, Benjamin [ORNL] (ORCID:0000000207400627)↗
DOE OSTI2026

Ecohydrology of urban environments

Urban watersheds include a gradient of highly managed to unmanaged ecosystems that vary in vegetation cover, built and natural drainage systems, and impervious cover. A distinct difference to undeveloped ecosystems is the dominance of material and energy import and export, the disruption of natural catenae, the covariation of soils, terrain and canopy conditions. In this chapter, we discuss the distinct characteristics of urban ecohydrology compared to non-urban systems. We also investigate emergent principles of unmanaged watersheds that underscore the promotion of ecohydrological resilience, and the need to abstract and adapt these principles to support the goal-oriented restoration of impaired urban ecosystems.

42 ENGINEERING↗
NASA NTRS1993

High temperature resistant jet fuels

The thermal stability of aviation fuels is evaluated according to the ASTM D 2341 method. This characteristic is linked to the presence of generators of instability, mainly diolephins and N and S derivates. The objectives of the present work are as follows: to obtain, in the laboratory, kerosene which is thermally stable; the isolation and identification of generators of instability; and the regeneration of thermally unstable kerosenes.

Luis M. Peloche↗
DOE OSTI2027

Low thermal budget dopant activation in shallow junctions of implanted Si via Tunable plasma enhanced annealing

Low thermal budget (LTB) annealing has become increasingly critical for shallow junctions as semiconductor devices are scaled down. Plasma Enhanced Annealing (PEA) has emerged as a promising LTB annealing process, however its mechanisms remain unclear. In this study, arsenic and boron implanted silicon wafers were subjected to helium or argon ion bombardment generated by annealing plasmas under controlled ion energies and doses in a home-built annealing reactor. The structural and electrical changes were characterized by four-point probe measurements, secondary ion mass spectrometry, Raman spectroscopy and spectroscopic ellipsometry. A pronounced reduction in sheet resistance, accompanied by a decrease in the damage layer thickness demonstrates the surface selective annealing capability of PEA process. Comparative studies using He and Ar plasmas and related ion bombardment reveal a dependence of activation efficacy on both ion species and the dopant implantation profile. Furthermore, a multi-step annealing mechanism is proposed, consisting of an initial (low dose) structural recovery stage followed by the generation of “extended” defects and concluded at large doses by an enhanced dopant activation. The proposed annealing kinetics are thought to be controlled by ion flux, dose and energy. Finally, these results highlight PEA as an effective, surface-selective, LTB approach for shallow-dopant activation and near surface annealing, and provide mechanistic insights into PEA processes.

dopant activation↗
DOE OSTI2026

The effect of dogbone sample size in the tensile testing of TATB-based plastic‑bonded explosive materials

Mechanical properties are of interest for many plastic‑bonded explosive (PBX) materials with tensile properties being of particular interest. Direct tensile measurements using dogbone-shaped samples are considered the gold standard, but they are fairly large, making testing more costly and less desirable from a safety perspective. We investigated whether the measured tensile strength depends on the dogbone specimen size, which to our knowledge, has not been reported in the literature for PBX materials. Understanding this should inform the feasibility of employing smaller samples and how sample size should be considered when comparing PBX dogbone values in the literature. The TATB-based PBX dogbone sample size was varied by (a) scaling all dimensions proportionally and (b) varying only the length of the samples. It was observed that the measured tensile peak stress (strength) was a function of the sample size, and was more dependent on the diameter (cross-sectional area) than the length of the samples. Since peak stress is calculated as peak force normalized to the diameter of the sample, one might not expect an explicit diameter dependence for the peak stress. Therefore, these results suggest there may be an additional strengthening effect as the sample diameter is increased.

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

A Novel Approach to Investigate Thermal Protection Systems Materials

The Koo Research Group (KRG) at The University of Texas at Austin (UT) and KAI has specialized in “Ablation Research” for more than fifteen years. Recently, the group has developed several incredibly unique capabilities that can advance “Thermal Protection Systems (TPS) Materials Research & Development” using an integrated experimental and numerical approach. The paper aims to introduce the methodology KRG has developed to solve this challenging problem. It will discuss how the KRG develops “Process-Properties-Performance” relationships of novel TPS materials in a systematical approach using (a) processing and fabrication, (b) thermal characterization of properties, (c) aerothermal testing, (d) microstructures characterization and analysis, and (e) numerical modeling. Progress and challenges of this research will also be discussed.

Engineering↗