SEARCH · Engineering Papers
Results for “domain walls”
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.
Numerical study of large-eddy breakup and its effect on the drag characteristics of boundary layers
The break-up of a field of eddies by a flat-plate obstacle embedded in a boundary layer is studied using numerical solutions to the two-dimensional Navier-Stokes equations. The flow is taken to be incompressible and unsteady. The flow field is initiated from rest. A train of eddies of predetermined size and strength are swept into the computational domain upstream of the plate. The undisturbed velocity profile is given by the Blasius solution. The disturbance vorticity generated at the plate and wall, plus that introduced with the eddies, mix with the background vorticity and is transported throughout the entire flow. All quantities are scaled by the plate length, the unidsturbed free-stream velocity, and the fluid kinematic viscosity. The Reynolds number is 1000, the Blasius boundary layer thickness is 2.0, and the plate is positioned a distance of 1.0 above the wall. The computational domain is four units high and sixteen units long.
PagMYB128 regulates secondary cell wall formation by direct activation of cell wall biosynthetic genes during wood formation in poplar
The biosynthesis of cellulose, lignin, and hemicelluloses in plant secondary cell walls (SCWs) is regulated by a hierarchical transcriptional regulatory network. Here, this network features orthologous transcription factors shared between poplar and Arabidopsis, highlighting a foundational similarity in their genetic regulation. However, knowledge on the discrepant behavior of the transcriptional-level molecular regulatory mechanisms between poplar and Arabidopsis remains limited. In this study, we investigated the function of PagMYB128 during wood formation and found it had broader impacts on SCW formation compared to its Arabidopsis ortholog, AtMYB103. Transgenic poplar trees overexpressing PagMYB128 exhibited significantly enhanced xylem development, with fiber cells and vessels displaying thicker walls, and an increase in the levels of cellulose, lignin, and hemicelluloses in the wood. In contrast, plants with dominant repression of PagMYB128 demonstrated the opposite phenotypes. RNA sequencing and reverse transcription – quantitative polymerase chain reaction showed that PagMYB128 could activate SCW biosynthetic gene expression, and chromatin immunoprecipitation along with yeast one-hybrid, and effector–reporter assays showed this regulation was direct. Further analysis revealed that PagSND1 (SECONDARY WALL-ASSOCIATED NAC-DOMAIN PROTEIN1) directly regulates PagMYB128 but not cell wall metabolic genes, highlighting the pivotal role of PagMYB128 in the SND1-driven regulatory network for wood development, thereby creating a feedforward loop in SCW biosynthesis.
Shear-induced polydomain structures of nematic lyotropic chromonic liquid crystal disodium cromoglycate
Lyotropic chromonic liquid crystals (LCLCs) represent aqueous dispersions of organic disk-like molecules that form cylindrical aggregates. Despite the growing interest in these materials, their flow behavior is poorly understood. Here, we explore the effect of shear on dynamic structures of the nematic LCLC, formed by 14 wt% water dispersion of disodium cromoglycate (DSCG). We employ in situ polarizing optical microscopy (POM) and small-angle and wide-angle X-ray scattering (SAXS/WAXS) to obtain independent and complementary information on the director structures over a wide range of shear rates. The DSCG nematic shows a shear-thinning behavior with two shear-thinning regions (Region I at $\dot{\gamma}$ < 1 s -1 and Region III at $\dot{\gamma}$ > 10 s -1 ) separated by a pseudo-Newtonian Region II (1 s-1 < $\dot{\gamma}$ < 10 s -1 ). The material is of a tumbling type. In Region I, $\dot{\gamma}$ < 1 s -1 , the director realigns along the vorticity axis. An increase of [small gamma, Greek, dot above] above 1 s -1 triggers nucleation of disclination loops. The disclinations introduce patches of the director that deviates from the vorticity direction and form a polydomain texture. Extension of the domains along the flow and along the vorticity direction decreases with the increase of the shear rate to 10 s -1 . Above 10 s -1 , the domains begin to elongate along the flow. At $\dot{\gamma}$ > 100 s -1 , the texture evolves into periodic stripes in which the director is predominantly along the flow with left and right tilts. The period of stripes decreases with an increase of $\dot{\gamma}$. The shear-induced transformations are explained by the balance of the elastic and viscous energies. In particular, nucleation of disclinations is associated with an increase of the elastic energy at the walls separating nonsingular domains with different director tilts. The uncovered shear-induced structural effects would be of importance in the further development of LCLC applications.
Inertial Effects and Anisotropy for the Flow in a Domain of Close Packed Spheres with a Bounding Wall
Not provided.
Ancient Origin of Acetyltransferases Catalyzing O -acetylation of Plant Cell Wall Polysaccharides
Abstract Members of the domain of unknown function 231/trichome birefringence–like (TBL) family have been shown to be O-acetyltransferases catalyzing the acetylation of plant cell wall polysaccharides, including pectins, mannan, xyloglucan and xylan. However, little is known about the origin and evolution of plant cell wall polysaccharide acetyltransferases. Here, we investigated the biochemical functions of TBL homologs from Klebsormidium nitens, a representative of an early divergent class of charophyte green algae that are considered to be the closest living relatives of land plants, and Marchantia polymorpha, a liverwort that is an extant representative of an ancient lineage of land plants. The genomes of K. nitens and Marchantia polymorpha harbor two and six TBL homologs, respectively. Biochemical characterization of their recombinant proteins expressed in human embryonic kidney 293 cells demonstrated that the two K. nitens TBLs exhibited acetyltransferase activities acetylating the pectin homogalacturonan (HG) and hence were named KnPOAT1 and KnPOAT2. Among the six M. polymorpha TBLs, five (MpPOAT1 to 5) possessed acetyltransferase activities toward pectins and the remaining one (MpMOAT1) catalyzed 2-O- and 3-O-acetylation of mannan. While MpPOAT1,2 specifically acetylated HG, MpPOAT3,4,5 could acetylate both HG and rhamnogalacturonan-I. Consistent with the acetyltransferase activities of these TBLs, pectins isolated from K. nitens and both pectins and mannan from M. polymorpha were shown to be acetylated. These findings indicate that the TBL genes were recruited as cell wall polysaccharide O-acetyltransferases as early as in charophyte green algae with activities toward pectins and they underwent expansion and functional diversification to acetylate various cell wall polysaccharides during evolution of land plants.
Molecular-gas-dynamics simulations of turbulent Couette flow over a mean-free-path-scale permeable substrate
Here we report flow statistics and visualizations from molecular-gas-dynamics simulations using the direct simulation Monte Carlo (DSMC) method for turbulent Couette flow in a minimal domain where the lower wall is replaced by an idealized permeable fibrous substrate representative of thermal-protection-system materials for which the Knudsen number is O(10 -1 ). Comparisons are made with smooth-wall DSMC simulations and smooth-wall direct numerical simulations (DNS) of the Navier-Stokes equations for the same conditions. Roughness, permeability, and noncontinuum effects are assessed. In the range of Reynolds numbers considered herein, the scalings of the skin friction on the permeable substrate and of the mean flow within the substrate suggest that they are dominated by viscous effects. While the regenerative cycle characteristic of smooth-wall turbulence remains intact for all cases considered, we observe that the near-wall velocity fluctuations are modulated by the permeable substrate with a wavelength equal to the pore spacing. Additionally, the flow within the substrate shows significant rarefaction effects, resulting in an apparent permeability that is 13% larger than the intrinsic permeability. In contrast, the smooth-wall DSMC and DNS simulations exhibit remarkably good agreement for the statistics examined, despite the Knudsen number based on the viscous length scale being as large as O(10 -1 ). This latter result is at variance with classical estimates for the breakdown of the continuum assumption and calls for further investigations into the interaction of noncontinuum effects and turbulence.
Enhancement of second-order optical nonlinearities and nanoscale periodic domain patterning in ferroelectric boron-substituted aluminum nitride thin films
The discovery and development of CMOS-compatible, nonlinear optical materials is essential to produce integrated photonic devices with advanced functionalities. AlN is a strong candidate for on-chip device demonstration due to its intrinsic second-order optical nonlinearities, large bandgap, and well-established fabrication techniques. However, AlN is not easily phase matched for the largest coefficient d 33 ; the coefficients that could potentially be dispersion phase-matched, d 31 and d 15 , have weak nonlinearities. This work investigates ferroelectric Al 1-x B x N (x = 0 to 0.11) for viability as a large bandgap nonlinear optical material with unique suitability towards ultraviolet light generation using second harmonic generation. The linear and nonlinear optical properties are characterized accounting for material anisotropy. With increasing B concentration, a large enhancement from near negligible values to d 31 = 0.9 ± 0.1 pm/V and d 15 = 1.2 ± 0.1 pm/V is observed. This compares favorably to other large bandgap materials like β-Ba(BO 2 ) 2 , where the largest nonlinear coefficient is d 22 ~ 2.3 pm/V at 800 nm. This is accompanied by a change in the bandgap from 6.1 eV to 5.8 eV as B substitution goes from 0 to 11%. A periodically poled, quasi-phase-matched ferroelectric domain pattern with 400 nm domain size and a wall roughness of <16 nm is demonstrated.
RANS Simulation of Variable Density Turbulent Round Jets with Coflow using xRAGE Hydrodynamic Code and BHR Turbulence Models
This work for the fiscal year 2023 (FY23) is a continuation of previous efforts to evaluate the BHR turbulence models for their ability to accurately simulate variable density turbulent round jets with coflow. As before, RANS simulations are carried out using the xRAGE hydrodynamic code. The following are some of the previous findings. Israel showed that i) the symmetry boundary conditions for the BHR models in axisymmetric simulations were in error, ii) three grids of different resolutions did not lead to converging solutions, and iii) BHR 3.1 simulation exhibited instabilities and did not reach a steady state. Saenz and Rauenzahn derived and implemented into xRAGE the correct BHR boundary conditions at the symmetry axis. Cline conducted sensitivity studies with various parameters including the BHR models (versions 2, 3.1 and 4), gravity, material pressure, specific heat, initial turbulent kinetic energy and initial turbulent length scale, and found that the largest factor impacting on simulation results was the BHR model version, followed by the initial turbulent length scale. In addition, freeze boundary conditions at the exit and the side wall of the computational domain were used to remove anomalous flow behavior. Cline adjusted the inlet jet velocity, initial turbulent kinetic energy and initial turbulent length scale to obtain the best reasonable match with the experimental data of Charonko and Prestridge. It was found that the BHR 2 and 3.1 models performed in a similar manner, but the BHR 2 model produced much lower levels of density-specific-volume covariance and turbulent kinetic energy. The main focus for the FY23 is to study the effects of computational parameters related to the boundary conditions, mesh, domain size and timestep size. The reasoning behind this is that, unless simulation results can be shown to be reasonably independent from the aforementioned computational parameters, it would be difficult to attribute any discrepancies between simulation and experimental results to turbulence models. This important aspect has largely been overlooked in the previous years, and therefore will be studied comprehensively here. Additionally, effects of varying the initial turbulent length scale will be examined because it was previously identified as a major factor affecting the flow fields.
Studies of heat-source driven natural convection: A numerical investigation
Thermal convection driven by uniform volumetric energy sources was studied in a horizontal fluid layer bounded from above by a rigid, isothermal surface and from below by a rigid, zero heat-flux surface. The side walls of the fluid domain were assumed to be rigid and perfectly insulating. The computations were formally restricted to two-dimensional laminar convection but were carried out for a range of Rayleigh numbers which spans the regimes of laminar and turbulent flow. The results of the computations consists of streamline and isotherm patterns, horizontally averaged temperature distributions, and horizontally averaged Nusselt numbers at the upper surface. Flow and temperature fields do not exhibit a steady state, but horizontally averaged Nusselt numbers reach limiting, quasi-steady values for all Rayleigh numbers considered. Correlations of the Nusselt number in terms of the Rayleigh and Prandtl numbers were determined.
Cell separations and the demixing of aqueous two phase polymer solutions in microgravity
Partition in phase separated aqueous polymer solutions is a cell separation procedure thought to be adversely influenced by gravity. In preparation for performing cell partitioning experiments in space, and to provide general information concerning the demixing of immiscible liquids in low gravity, a series of phase separated aqueous polymer solutions have been flown on two shuttle flights. Fluorocarbon oil and water emulsions were also flown on the second flight. The aqueous polymer emulsions, which in one g demix largely by sedimentation and convection due to the density differences between the phases, demixed more slowly than on the ground and the final disposition of the phases was determined by the wetting of the container wall by the phases. The demixing behavior and kinetics were influenced by the phase volume ratio, physical properties of the systems and chamber wall interaction. The average domain size increased linearly with time as the systems demixed.
Lattice Boltzmann Method for Spacecraft Propellant Slosh Simulation
A scalable computational approach to the simulation of propellant tank sloshing dynamics in microgravity is presented. In this work, we use the lattice Boltzmann equation (LBE) to approximate the behavior of two-phase, single-component isothermal flows at very low Bond numbers. Through the use of a non-ideal gas equation of state and a modified multiple relaxation time (MRT) collision operator, the proposed method can simulate thermodynamically consistent phase transitions at temperatures and density ratios consistent with typical spacecraft cryogenic propellants, for example, liquid oxygen. Determination of the tank forces and moments relies upon the global momentum conservation of the fluid domain, and a parametric wall wetting model allows tuning of the free surface contact angle. Development of the interface is implicit and no interface tracking approach is required. Numerical examples illustrate the method's application to predicting bulk fluid motion including lateral propellant slosh in low-g conditions.
Lattice Boltzmann Method for Spacecraft Propellant Slosh Simulation
A scalable computational approach to the simulation of propellant tank sloshing dynamics in microgravity is presented. In this work, we use the lattice Boltzmann equation (LBE) to approximate the behavior of two-phase, single-component isothermal flows at very low Bond numbers. Through the use of a non-ideal gas equation of state and a modified multiple relaxation time (MRT) collision operator, the proposed method can simulate thermodynamically consistent phase transitions at temperatures and density ratios consistent with typical spacecraft cryogenic propellants, for example, liquid oxygen. Determination of the tank forces and moments is based upon a novel approach that relies on the global momentum conservation of the closed fluid domain, and a parametric wall wetting model allows tuning of the free surface contact angle. Development of the interface is implicit and no interface tracking approach is required. A numerical example illustrates the method's application to prediction of bulk fluid behavior during a spacecraft ullage settling maneuver.
Dynamic Smagorinsky Modeled Large-Eddy Simulations of Turbulence Using Tetrahedral Meshes
Eddy-resolving numerical computations of turbulent flows are emerging as viable alternatives to Reynolds Averaged Navier-Stokes (RANS) calculations for flows with an intrinsically steady mean state due to the advances in large-scale parallel computing. In these computations, medium to large turbulent eddies are resolved by the numerics while the smaller or subgrid scales are either modeled or taken care of by the inherent numerical dissipation. To advance the state of the art of unstructured-mesh turbulence simulation capabilities, large eddy simulations (LES) using the dynamic Smagorinsky model (DSM) on tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method. In contrast to what has been reported in the literature, the present implementation of dynamic models allows for active backscattering without any ad-hoc limiting of the eddy viscosity calculated from the subgrid-scale model. For the benchmark problems involving compressible isotropic turbulence decay as well as the shock/turbulent boundary layer interaction benchmark problems, no numerical instability associated with kinetic energy growth is observed and the volume percentage of the backscattering portion accounts for about 38-40% of the simulation domain. A slip-wall model in conjunction with the implemented DSM is used to simulate a relatively high Reynolds number Mach 2.85 turbulent boundary layer over a 30° ramp with several tetrahedral meshes and a wall-normal spacing of either Δγ& = 10 or Δγ& = 20. The computed mean wall pressure distribution, separation region size, mean velocity profiles, and Reynolds stress agree reasonably well with experimental data.
Modelling above-ground biomass stock over Norway using national forest inventory data with ArcticDEM and Sentinel-2 data
Boreal forests constitute a large portion of the global forest area, yet they are undersampled through field surveys, and only a few remotely sensed data sources provide structural information wall-to-wall throughout the boreal domain. ArcticDEM is a collection of high-resolution (2 m) space-borne stereogrammetric digital surface models (DSM) covering the entire land area north of 60° of latitude. The free-availability of ArcticDEM data offers new possibilities for aboveground biomass mapping (AGB) across boreal forests, and thus it is necessary to evaluate the potential for these data to map AGB over alternative open-data sources (i.e., Sentinel-2). This study was performed over the entire land area of Norway north of 60° of latitude, and the Norwegian national forest inventory (NFI) was used as a source of field data composed of accurately geolocated field plots (n=7710) systematically distributed across the study area. Separate random forest models were fitted using NFI data, and corresponding remotely sensed data consisting of either: i) a canopy height model (ArcticCHM) obtained by subtracting a high-quality digital terrain model (DTM) from the ArcticDEM DSM height values, ii) Sentinel-2 (S2), or iii) a combination of the two (ArcticCHM+S2). Furthermore, we assessed the effect of the forest- and terrain-specific factors on the models’ predictive accuracy. The best model (,i.e., ArcticCHM+S2) explained nearly 60% of the variance of the training set, which translated in the largest accuracy in terms of root mean square error (RMSE=41.4 t/ha). This result highlights the synergy between 3D and multispectral data in AGB modelling. Furthermore, this study showed that despite the importance of ArcticCHM variables, the S2 model performed slightly better than ArcticCHM model. This finding highlights some of the limitations of ArcticDEM, which, despite the unprecedented spatial resolution, is highly heterogeneous due to the blending of multiple acquisitions across different years and seasons. We found that both forest- and terrain-specific characteristics affected the uncertainty of the ArcticCHM+S2 model and concluded that the combined use of ArcticCHM and Sentinel-2 represents a viable solution for AGB mapping across boreal forests. The synergy between the two data sources allowed for a reduction of the saturation effects typical of multispectral data while ensuring the spatial consistency in the output predictions due to the removal of artifacts and data voids present in ArcticCHM data. While the main contribution of this study is to provide the first evidence of the best-case-scenario (i.e., availability of accurate terrain models) that ArcticDEM data can provide for large-scale AGB modelling, it remains critically important for other studies to investigate how ArcticDEM may be used in areas where no DTMs are available as is the case for large portions of the boreal zone.
Modelling of edge plasma dynamics with active wall boundary conditions
A self-consistent 2D model is presented for transport in boundary plasma and plasma-facing material walls. Plasma dynamics in the domain is represented by a 2D collisional plasma fluid model in the edge-plasma code UEDGE, and transport of hydrogen and heat in the wall is represented by a system of reaction–diffusion equations in the 1D wall code FACE. To account for variation of parameters along the wall, in the coupled model multiple instances of the FACE code run in parallel. Here, the coupled model provides a tool for investigating a range of dynamic plasma–material interactions phenomena in 2D. For demonstration of its capability, one application of particular interest is the role of active wall in tokamak strike point sweeping proposed for mitigation of divertor heat loads. In the present study, the coupled calculations are applied to investigation of the impact of heat and hydrogen transport in the material wall on the divertor plasma and target heat load during sweeping of the target strike point for parameters of a high-power tokamak.
Simulation of RCC Crack Growth Due to Carbon Oxidation in High-Temperature Gas Environments
The carbon wall oxidation technique coupled with a CFD technique was employed to study the flow in the expanding crack channel caused by the oxidation of the channel carbon walls. The recessing 3D surface morphing procedure was developed and tested in comparison with the arcjet experimental results. The multi-block structured adaptive meshing was used to model the computational domain changes due to the wall recession. Wall regression rates for a reinforced carbon-carbon (RCC) samples, that were tested in a high enthalpy arcjet environment, were computationally obtained and used to assess the channel expansion. The test geometry and flow conditions render the flow regime as the transitional to continuum, therefore Navier-Stokes gas dynamic approach with the temperature jump and velocity slip correction to the boundary conditions was used. The modeled mechanism for wall material loss was atomic oxygen reaction with bare carbon. The predicted channel growth was found to agree with arcjet observations. Local gas flow field results were found to affect the oxidation rate in a manner that cannot be predicted by previous mass loss correlations. The method holds promise for future modeling of materials gas-dynamic interactions for hypersonic flight.
Thermocryogenic buckling and stress analyses of a partially filled cryogenic tank subjected to cylindrical strip heating
Thermocryogenic buckling and stress analyses were conducted on a horizontally oriented cryogenic tank using the finite element method. The tank is a finite-length circular cylindrical shell with its two ends capped with hemispherical shells. The tank is subjected to cylindrical strip heating in the region above the liquid-cryogen fill level and to cryogenic cooling below the fill level (i.e., under thermocryogenic loading). The effects of cryogen fill level on the buckling temperature and thermocryogenic stress field were investigated in detail. Both the buckling temperature and stress magnitudes were relatively insensitive to the cryogen fill level. The buckling temperature, however, was quite sensitive to the radius-to-thickness ratio. A mechanical stress analysis of the tank also was conducted when the tank was under: (1) cryogen liquid pressure loading; (2) internal pressure loading; and (3) tank-wall inertia loading. Deformed shapes of the cryogenic tanks under different loading conditions were shown, and high-stress domains were mapped on the tank wall for the strain-gage installations. The accuracies of solutions from different finite element models were compared.