Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Approximate boundary”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Slip velocity boundary conditions for the lattice Boltzmann modeling of microchannel flows

Abstract Slip flows in ducts are important in numerous engineering applications, most notably in microchannel flows. Compared to the standard no‐slip Dirichlet condition, the case of slip formulates as a Robin‐type condition for the fluid tangential velocity. Such an increase in mathematical complexity is accompanied by a more challenging numerical transcription. The present work concerns with this topic, addressing the modeling of the slip velocity boundary condition in the lattice Boltzmann method (LBM) applied to steady slow viscous flows inside ducts of nontrivial shapes. As novelty, we extend the newly revised local second‐order boundary (LSOB) Dirichlet fluid flow method [ Philos. Trans. R. Soc. A 378, 20190404 (2020)] to implement the slip velocity condition within the two‐relaxation‐time (TRT) framework. The LSOB follows an in‐node philosophy where its operation principle seeks to explicitly reconstruct the unknown boundary populations in the form of a third‐order accurate Chapman–Enskog expansion, where the wall slip condition is built‐in as a normal Taylor‐type condition. The key point of this approach is that the required first‐ and second‐order momentum derivatives, rather than computed through nonlocal finite difference approximations, are locally determined through a simple local linear algebra procedure, whose formulation is particularly aided by the TRT symmetry argument. To express the obtained derivatives, two approaches are considered, called and , which operate with node and wall variables, respectively. These two formulations are developed to prescribe the physical slip condition over plane and curved walls, including the corners. Their consistency and accuracy characteristics are examined against alternative linkwise strategies to impose the wall slip velocity, such as the kinetic‐based diffusive bounce‐back scheme, the central linear interpolation slip scheme, and the multireflection slip scheme. The several slip schemes are tested over different 3D microchannel configurations, with walls not conforming with the LBM uniform mesh. Numerical tests confirm the advanced accuracy characteristics of the proposed LSOB slip boundary scheme, revealing the added challenge of the wall slip modeling, and that parabolic accuracy is a necessary requirement to reach second‐order accuracy within this problem class.

Silva, Goncalo↗
DOE OSTI2023

Stochastic Real-Time Second-Order Green’s Function Theory for Neutral Excitations in Molecules and Nanostructures

We present a real-time second-order Green’s function (GF) method for computing excited states in molecules and nanostructures, with a computational scaling of O(N e 3 ), where N e is the number of electrons. The cubic scaling is achieved by adopting the stochastic resolution of the identity to decouple the 4-index electron repulsion integrals. To improve the time propagation and the spectral resolution, we adopt the dynamic mode decomposition technique and assess the accuracy and efficiency of the combined approach for a chain of hydrogen dimer molecules of different lengths. We find that the stochastic implementation accurately reproduces the deterministic results for the electronic dynamics and excitation energies. Furthermore, we provide a detailed analysis of the statistical errors, bias, and long-time extrapolation. Overall, the approach offers an efficient route to investigate excited states in extended systems with open or closed boundary conditions.

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

On the symmetries of cosmological perturbations

The space of inflationary models is vast, containing wide varieties of mechanisms, symmetries, and spectra of particles. Consequently, the space of observational signatures is similarly complex. Hence, it is natural to look for boundaries of the space of models and their signatures. In this paper, we explore the possible symmetries associated with the primordial cosmological perturbations and their correlators in the asymptotic future. Assuming the observed homogeneity, isotropy and (approximate) scale invariance, we prove three main results. First, correlation functions of scalar metric fluctuations are uniquely characterized by soft theorems and are free from ambiguity under field redefinitions. Second, whatever the particle content and interactions, when the standard soft theorems apply, invariance under de Sitter boosts (linearly realized conformal invariance) is only possible if all connected correlators vanish identically, i.e. if the theory is free. Third, conformal invariance is the largest set of linearly realized (bosonic) symmetries of the correlators of any single scalar, irrespectively of any soft theorems or particle content.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Properties of Arctic liquid and mixed-phase clouds from shipborne Cloudnet observations during ACSE 2014

This article presents Cloudnet retrievals of Arctic clouds from measurements conducted during a 3-month research expedition along the Siberian shelf during summer and autumn 2014. During autumn, we find a strong reduction in the occurrence of liquid clouds and an increase for both mixed-phase and ice clouds at low levels compared to summer. About 80 % of all liquid clouds observed during the research cruise show a liquid water path below the infrared black body limit of approximately 50 g m -2 . The majority of mixed-phase and ice clouds had an ice water path below 20 g m -2 . Cloud properties are analysed with respect to cloud-top temperature and boundary layer structure. Changes in these parameters have little effect on the geometric thickness of liquid clouds while mixed-phase clouds during warm-air advection events are generally thinner than when such events were absent. Cloud-top temperatures are very similar for all mixed-phase clouds. However, more cases of lower cloud-top temperature were observed in the absence of warm-air advection. Profiles of liquid and ice water content are normalized with respect to cloud base and height. For liquid water clouds, the liquid water content profile reveals a strong increase with height with a maximum within the upper quarter of the clouds followed by a sharp decrease towards cloud top. Liquid water content is lowest for clouds observed below an inversion during warm-air advection events. Most mixed-phase clouds show a liquid water content profile with a very similar shape to that of liquid clouds but with lower maximum values during events with warm air above the planetary boundary layer. The normalized ice water content profiles in mixed-phase clouds look different from those of liquid water content. They show a wider range in maximum values with the lowest ice water content for clouds below an inversion and the highest values for clouds above or extending through an inversion. The ice water content profile generally peaks at a height below the peak in the liquid water content profile – usually in the centre of the cloud, sometimes closer to cloud base, likely due to particle sublimation as the crystals fall through the cloud.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Computational Discovery of Intermolecular Singlet Fission Materials Using Many-Body Perturbation Theory

Intermolecular singlet fission (SF) is the conversion of a photogenerated singlet exciton into two triplet excitons residing on different molecules. SF has the potential to enhance the conversion efficiency of solar cells by harvesting two charge carriers from one high-energy photon, whose surplus energy would otherwise be lost to heat. The development of commercial SF-augmented modules is hindered by the limited selection of molecular crystals that exhibit intermolecular SF in the solid state. Computational exploration may accelerate the discovery of new SF materials. The GW approximation and Bethe–Salpeter equation (GW+BSE) within the framework of many-body perturbation theory is the current state-of-the-art method for calculating the excited-state properties of molecular crystals with periodic boundary conditions. In this Review, we discuss the usage of GW+BSE to assess candidate SF materials as well as its combination with low-cost physical or machine learned models in materials discovery workflows. We demonstrate three successful strategies for the discovery of new SF materials: (i) functionalization of known materials to tune their properties, (ii) finding potential polymorphs with improved crystal packing, and (iii) exploring new classes of materials. In addition, three new candidate SF materials are proposed here, which have not been published previously.

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

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Dynamical cluster approximation study of electron localization in the extended Hubbard model

We perform a detailed study of the phase transitions and mechanisms of electron localization in the extended Hubbard model using the dynamical cluster approximation on a 2×2 cluster. We explore the interplay of charge order and Mott physics. We find that a nearest-neighbor Coulomb interaction V causes “screening” effects close to the Mott phase transition, pushing the phase boundary to larger values of U. We also demonstrate the different effects of V on correlations in metallic and insulating regimes, and document the different correlation aspects of charge order and Mott states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

NSFnets (Navier-Stokes flow nets): Physics-informed neural networks for the incompressible Navier-Stokes equations

In the last 50 years there has been a tremendous progress in solving numerically the Navier-Stokes equations using finite differences, finite elements, spectral, and even meshless methods. Yet, in many real cases, we still cannot incorporate seamlessly (multi-fidelity) data into existing algorithms, and for industrial-complexity applications the mesh generation is time consuming and still an art. Moreover, solving ill-posed problems (e.g., lacking boundary conditions) or inverse problems is often prohibitively expensive and requires different formulations and new computer codes. Here, we employ physics-informed neural networks (PINNs), encoding the governing equations directly into the deep neural network via automatic differentiation, to overcome some of the aforementioned limitations for simulating incompressible laminar and turbulent flows. We develop the Navier-Stokes flow nets (NSFnets) by considering two different mathematical formulations of the Navier-Stokes equations: the velocity-pressure (VP) formulation and the vorticity-velocity (VV) formulation. Since this is a new approach, we first select some standard benchmark problems to assess the accuracy, convergence rate, computational cost and flexibility of NSFnets; analytical solutions and direct numerical simulation (DNS) databases provide proper initial and boundary conditions for the NSFnet simulations. The spatial and temporal coordinates are the inputs of the NSFnets, while the instantaneous velocity and pressure fields are the outputs for the VP-NSFnet, and the instantaneous velocity and vorticity fields are the outputs for the VV-NSFnet. This is unsupervised learning and, hence, no labeled data are required beyond boundary and initial conditions and the fluid properties. The residuals of the VP or VV governing equations, together with the initial and boundary conditions, are embedded into the loss function of the NSFnets. No data is provided for the pressure to the VP-NSFnet, which is a hidden state and is obtained via the incompressibility constraint without extra computational cost. Unlike the traditional numerical methods, NSFnets inherit the properties of neural networks (NNs), hence the total error is composed of the approximation, the optimization, and the generalization errors. Here, we empirically attempt to quantify these errors by varying the sampling (“residual”) points, the iterative solvers, and the size of the NN architecture. For the laminar flow solutions, we show that both the VP and the VV formulations are comparable in accuracy but their best performance corresponds to different NN architectures. The initial convergence rate is fast but the error eventually saturates to a plateau due to the dominance of the optimization error. For the turbulent channel flow, we show that NSFnets can sustain turbulence at , but due to expensive training we only consider part of the channel domain and enforce velocity boundary conditions on the subdomain boundaries provided by the DNS data base. We also perform a systematic study on the weights used in the loss function for balancing the data and physics components, and investigate a new way of computing the weights dynamically to accelerate training and enhance accuracy. In the last part, we demonstrate how NSFnets should be used in practice, namely for ill-posed problems with incomplete or noisy boundary conditions as well as for inverse problems. We obtain reasonably accurate solutions for such cases as well without the need to change the NSFnets and at the same computational cost as in the forward well-posed problems. As a result, we also present a simple example of transfer learning that will aid in accelerating the training of NSFnets for different parameter settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Taylor approximation variance reduction for approximation errors in PDE-constrained Bayesian inverse problems

In numerous applications, surrogate models are used as a replacement for accurate parameter-to-observable mappings when solving large-scale inverse problems governed by partial differential equations (PDEs). The surrogate model may be a computationally cheaper alternative to the accurate parameter-to-observable mappings and/or may ignore additional unknowns or sources of uncertainty. The Bayesian approximation error (BAE) approach provides a means to account for the induced uncertainties and approximation errors, i.e. the errors between the accurate parameter-to-observable mapping and the surrogate. The statistics of these errors are, however, in general unknown a priori, and are thus calculated using Monte Carlo sampling. Although the sampling is typically carried out offline, i.e. before considering the data, the process can still represent a computational bottleneck. In this work, we develop a scalable computational approach for reducing the costs associated with the sampling stage of the BAE approach. Specifically, we consider the Taylor expansion of the accurate and surrogate forward models with respect to the uncertain parameter fields either as a control variate for variance reduction or as a means to directly and efficiently approximate the mean and covariance of the approximation errors. We propose efficient methods for evaluating the expressions for the mean and covariance of the Taylor approximations based on linear(-ized) PDE solves. Furthermore, the proposed approach is independent of the dimension of the uncertain parameter, depending instead on the intrinsic dimension of the data, ensuring scalability to high-dimensional problems. The potential benefits of the proposed approach are demonstrated for two high-dimensional inverse problems governed by PDE examples, namely for the estimation of a distributed Robin boundary coefficient in a linear diffusion problem, and for a coefficient estimation problem governed by a nonlinear diffusion problem.

Bayesian approximation error↗
DOE OSTI2025

High‐Speed Slot‐Die Coating with Donor‐Priority Rapid Aggregation Kinetics for Improved Morphology and Efficiency in Ecofriendly Organic Solar Cells

Abstract Solution‐processable organic solar cells (OSCs) represent a promising renewable photovoltaic technology with significant potential for eco‐compatible production. While high power conversion efficiencies (PCEs) have been achieved in OSCs, scaling this technology for high‐throughput manufacturing remains challenging. Key reason lies in the lack of efficient control strategies for the complex and long‐duration morphology evolution during high‐speed coating process with ecofriendly solvents. Here, a donor‐priority rapid aggregation process (DP‐RAP) scheme is proposed to solve this issue by adjusting the aggregation kinetics of donor and acceptor components. DP‐RAP enables blends with a nanoscale fiber network structure and favorable crystallinity, which contributes to balanced carrier transport and reduced recombination losses. As a result, the PCE is improved from 14.3% (reference) to 17.4% (DP‐RAP) for ultra‐high speed coated PM6:BTP‐eC9 devices in atmosphere, which is one of the highest values for non‐halogenated solvent‐processed solar cells at coating speeds of 500 mm s −1 . Moreover, the DP‐RAP based devices remain a stable PCE of approximately 17.4% across a broad range of coating speeds (20–500 mm s −1 ), illustrating its tolerance to the varied manufacturing conditions. This work highlights a promising avenue for the high‐speed, ecofriendly production of efficient OSCs, pushing the boundaries of practical manufacturing in renewable energy technologies.

Chemistry↗
DOE OSTI2026

Evaluating CO 2 mitigation strategies in SAF biorefineries: Techno-economic and life cycle analysis

The aviation sector requires scalable decarbonization strategies, and lignocellulosic sustainable aviation fuel (SAF) represents a promising pathway. This study comparatively evaluates the techno-economic analysis and life cycle assessment (LCA) of three CO 2 management strategies integrated within a U.S.-based gasification–Fischer–Tropsch SAF biorefinery: (i) catalytic hydrogenation of captured CO 2 to methanol, (ii) geological CO 2 sequestration, and (iii) mineralization to sodium bicarbonate (NaHCO 3 ). Techno-economic analysis indicates that methanol synthesis requires approximately 26% higher capital investment and 33% higher operating costs than mineralization. Although methanol co-production generates the highest gross revenue, NaHCO 3 production reduces the SAF minimum selling price by approximately 38% relative to both methanol synthesis and geological sequestration pathways, reflecting a more balanced cost allocation through mineral co-product valorization. Geological sequestration lowers operating costs by nearly 50% compared with methanol synthesis but remains highly dependent on carbon credit mechanisms. LCA reveals substantial divergence in climate performance. Relative to methanol synthesis, sequestration improves net greenhouse gas performance by approximately 163%, transitioning the system from net-positive to net-negative emissions. Mineralization further enhances carbon mitigation, achieving roughly 85% greater carbon reduction than sequestration and over sixfold improvement relative to methanol synthesis within the defined system boundary. Sensitivity analysis identified hydrogen price, co-product market value, and process emissions as dominant drivers. Under baseline assumptions, CO 2 mineralization is found to offer the most balanced pathway.

Carbon capture and storage↗
DOE OSTI2025

Theory and numerics of subspace approximation of eigenvalue problems

Large-scale eigenvalue problems arise in various fields of science and engineering and demand computationally efficient solutions. In this study, we investigate the subspace approximation for parametric linear eigenvalue problems, aiming to mitigate the computational burden associated with high-fidelity systems. Furthermore, we provide general error estimates under non-simple eigenvalue conditions, establishing some theoretical foundations for understanding the convergence behavior of subspace approximations. Numerical examples, including problems with one-dimensional to three-dimensional spatial domain and one-dimensional to two-dimensional parameter domain, are presented to demonstrate the efficacy of reduced basis method in handling parametric variations in boundary conditions and coefficient fields to achieve significant computational savings while maintaining high accuracy, making them promising tools for practical applications in large-scale eigenvalue computations.

Eigenvalue problems↗
DOE OSTI2022

Turbulent broadening of electron heat-flux width in electromagnetic gyrokinetic simulations of a helical scrape-off layer model

Here, we demonstrate that cross field transport in the scrape-off layer (SOL) can be moderately increased by electromagnetic effects in high-beta regimes, resulting in broadening of the electron heat-flux width on the endplates. This conclusion is taken from full- f electromagnetic gyrokinetic simulations of a helical SOL model that roughly approximates the SOL of the National Spherical Torus Experiment. The simulations have been performed with the Gkeyll code, which recently became the first code to demonstrate the capability to simulate electromagnetic gyrokinetic turbulence on open magnetic field lines with sheath boundary conditions. We scan the source rate and thus β, so that the normalized pressure gradient (the MHD ballooning parameter $α ∝ ∂β/∂r ∝ β/Lρ$) is scanned over an experimentally relevant range, α = 0.3 - 1.5. While there is little change in the pressure gradient scale length L p near the midplane as beta is increased, a 10% increase in cross field transport near the midplane results in an increase in the electron heat-flux width λ q and a 25% reduction of the peak electron heat flux to the endplates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Structural Aspects of Neutron Survival Probabilities

The neutron survival probability (and related quantities including probabilities of extinction and initiation) is a central element of the broader stochastic theory of neutron populations and finds application in fields including reactor start-up, analysis of reactor power bursts and criticality accidents, and safeguards. In a full neutron transport formulation, the equation governing the single-neutron survival probability is a backward or adjoint-like integro-partial differential equation with the added complexity of being highly nonlinear. Analogous formulations of this equation exist in the context of many approximate theories of neutron transport, with the point kinetics formulation having received significant theoretical attention since the 1940s. This work continues this tradition by providing a novel analysis of the single-neutron survival probability equation using the tools of boundary layer theory. The analysis reveals that the “fully dynamic” solution of the single-neutron survival probability equation—and some key probability distributions derived from it—may be cast as a singular perturbation around the underlying quasi-static single-neutron probability of initiation. In this perturbation solution, the expansion parameter is the ratio of the neutron generation time to a macroscopic time scale characterizing the overall system evolution; this interpretation illuminates some of the fundamental structural aspects of neutron survival phenomena.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Physics‐Informed Machine Learning Method for Large‐Scale Data Assimilation Problems

Abstract We develop a physics‐informed machine learning approach for large‐scale data assimilation and parameter estimation and apply it for estimating transmissivity and hydraulic head in the two‐dimensional steady‐state subsurface flow model of the Hanford Site given synthetic measurements of said variables. In our approach, we extend the physics‐informed conditional Karhunen‐Loéve expansion (PICKLE) method to modeling subsurface flow with unknown flux (Neumann) and varying head (time‐dependent Dirichlet) boundary conditions. We demonstrate that the PICKLE method is comparable in accuracy with the standard maximum a posteriori (MAP) method, but is significantly faster than MAP for large‐scale problems. Both methods use a mesh to discretize the computational domain. In MAP, the parameters and states are discretized on the mesh; therefore, the size of the MAP parameter estimation problem directly depends on the mesh size. In PICKLE, the mesh is used to evaluate the residuals of the governing equation, while the parameters and states are approximated by the truncated conditional Karhunen‐Loéve expansions with the number of parameters controlled by the smoothness of the parameter and state fields, and not by the mesh size. For a considered example, we demonstrate that the computational cost of PICKLE increases near linearly (as N 1.15 ) with the number of grid nodes N , while that of MAP increases much faster (as N 3.28 ). We also show that once trained for one set of Dirichlet boundary conditions (i.e., one river stage), the PICKLE method provides accurate estimates of the hydraulic head for any value of the Dirichlet boundary conditions (i.e., for any river stage).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Precise stellarator quasi-symmetry can be achieved with electromagnetic coils

Magnetic fields with quasi-symmetry are known to provide good confinement of charged particles and plasmas, but the extent to which quasi-symmetry can be achieved in practice has remained an open question. Recent work [M. Landreman and E. Paul, Phys. Rev. Lett. 128, 035001, 2022] reports the discovery of toroidal magnetic fields that are quasi-symmetric to orders-of-magnitude higher precision than previously known fields. We show that these fields can be accurately produced using electromagnetic coils of only moderate engineering complexity, that is, coils that have low curvature and that are sufficiently separated from each other. Our results demonstrate that these new quasi-symmetric fields are relevant for applications requiring the confinement of energetic charged particles for long time scales, such as nuclear fusion. The coils’ length plays an important role for how well the quasi-symmetric fields can be approximated. For the longest coil set considered and a mean field strength of 1 T, the departure from quasi-symmetry is of the order of Earth’s magnetic field. Additionally, we find that magnetic surfaces extend far outside the plasma boundary used by Landreman and Paul, providing confinement far from the core. Simulations confirm that the magnetic fields generated by the new coils confine particles with high kinetic energy substantially longer than previously known coil configurations. In particular, when scaled to a reactor, the best found configuration loses only 0.04% of energetic particles born at midradius when following guiding center trajectories for 200 ms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Enhanced thermoelectric efficiency in nanocrystalline bismuth telluride nanotubes

We report on the thermal and thermoelectric properties of individual nanocrystalline Bi 2 Te 3 nanotubes synthesized by the solution phase method using 3ω method and a microfabricated testbench. Measurements show that the nanotubes offer improved ZT compared to bulk Bi 2 Te 3 near room temperature due to an enhanced Seebeck coefficient and suppressed thermal conductivity. This improvement in ZT originates from the nanocrystalline nature and low dimensionality of the nanotubes. Domain boundary filtering of low-energy electrons provides an enhanced Seebeck coefficient. The scattering of phonons at the surface of the nanotube leads to suppressed thermal conductivity. These have been theoretically analyzed using the Boltzmann equation based on the relaxation time approximation and Landauer approach. This work clearly demonstrates the possibility of achieving enhancement in thermoelectric efficiency by combining nanocrystalline and low-dimensional systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Strain-Driven Evolution of Structural and Carrier Recombination Dynamics in Germanium

Strain engineering in group-IV semiconductors, in particular germanium (Ge), offers a promising route for advancing next-generation optoelectronic device performance. High minority carrier lifetime serves as a direct indicator of superior material quality and device efficiency. Several factors govern the minority carrier lifetime in both strained Ge and strained GeSn, such as strain level, epilayer thickness, defects, surface roughness, and Sn alloying. In this study, the impact of biaxial strain on effective minority carrier lifetimes of tensile-strained Ge (0% ..epsilon..-Ge/Al?As, 0.78% ..epsilon..-Ge/In0.115?Ga0.885?As, and 0.93% ..epsilon..-Ge/In0.16?Ga0.84?As) and compressively strained GeSn (-0.53% c-Ge0.97?Sn0.03/Al?As) epilayers, grown by molecular beam epitaxy, was experimentally investigated via the microwave-reflection photoconductive decay (..mu..-PCD) technique as a simultaneous function of thickness and strain. The structural and strain-relaxation properties were analyzed by high-resolution x-ray diffraction and Raman spectroscopy, while surface morphology analysis using atomic force microscopy provided insights into the surface roughness of strained Ge-based epilayers. Defect characterization by cross-section transmission electron microscopy revealed superior crystalline quality of ..epsilon..-Ge in one epilayer, while twin boundary formation was observed in the other. We investigated the carrier recombination dynamics in epitaxially strained Ge by extracting the bulk recombination lifetime and surface recombination velocity using the thickness-dependent minority carrier lifetime. High effective minority carrier lifetimes of approximately 60-340 ns were achieved at high strain levels with minimal dependence on the thickness of the strained Ge-based epilayers by means of the ..mu..-PCD technique. Minority carrier lifetimes in Ge-based epilayers exhibit a quadratic dependence on the biaxial strain, ranging from -0.88% (compressive) to +1.6% (tensile). The analysis of the carrier recombination dynamics in strained Ge-based heterostructures offers a pathway to designing group-IV based optoelectronic devices.

36 MATERIALS SCIENCE↗