Transition Front Research Articles

KPP transition fronts in a one-dimensional two-patch habitat.

This paper is concerned with the existence of transition fronts for a one-dimensional two patch model with KPP reaction terms. Density and flux conditions are imposed at the interface between the two patches. We first construct a pair of suitable super- and sub solutions by making full use of information of the leading edges of two KPP fronts and gluing them through the interface conditions. Then, an entire solution obtained thanks to a limiting argument is shown to be a transition front moving from one patch to the other one. This propagating solution admits asymptotic past and future speeds, and it connects two different fronts, each associated with one of the two patches. The paper thus provides the first example of a transition front for a KPP-type two-patch model with interface conditions.

High-resolution Observations of H i in the IC 63 Reflection Nebula

Photodissociation regions (PDRs), where the (far-)ultraviolet light from hot young stars interact with the gas in surrounding molecular clouds, provide laboratories for understanding the nature and role of feedback by star formation on the interstellar medium. While the general nature of PDRs is well understood—at least under simplified conditions—the detailed dynamics and chemistry of these regions, including gas clumping, evolution over time, etc., can be very complex. We present interferometric observations of the 21 cm atomic hydrogen line, combined with [C ii] 158 μm observations, toward the nearby reflection nebula IC 63. We find a clumpy H i structure in the PDR, and a ring morphology for the H i emission at the tip of IC 63. We further unveil kinematic substructure, of the order of 1 km s−1, in the PDR layers and several legs that will disperse IC 63 in <0.5 Myr. We find that the dynamics in the PDR explain the observed clumpy H i distribution and lack of a well-defined H i/H2 transition front. However, it is currently not possible to conclude whether H i self-absorption and nonequilibrium chemistry also contribute to this clumpy morphology and missing H i/H2 transition front.

Modal Instabilities over Blunted Cones at Angle of Attack in Hypersonic Flow

The effects of nose radius and angle of attack on the linear modal amplification over blunt circular cones at hypersonic speeds are computationally investigated. The three-dimensional laminar flow solutions over a 1.5-m-long, 7°-half-angle cone with 5.080, 9.525, and 25.40 mm nosetip radii are computed for selected angle-of-attack values and freestream conditions that match the Mach 10 experiments at a freestream unit Reynolds number of 17.1 million per meter conducted within the Hypervelocity Wind Tunnel 9 at the Air Force Arnold Engineering Development Complex (AEDC). Results indicate that the linear amplification of stationary and traveling crossflow waves along inflection lines increases with the angle of attack and decreases with the nosetip radius. The trend in transition front with respect to increasing angle of attack and bluntness is found to be consistent with the predicted increase in the amplification factors for Mack’s second mode (MM) disturbances along the streamline trajectories for the small and medium bluntness cases. The increase in the MM amplification along the windward ray for higher angles of attack is shown to be the result of a progressively earlier entropy-layer swallowing. Computations also indicate that the transition amplification factor along the windward ray increases with the angle of attack and decreases with the nosetip radius. Furthermore, the transition -factors along the leeward ray are below 2 and, therefore, rather small to lead to transition.

Compressibility corrections to extend a smooth local correlation-based transition model to transonic flows

AbstractThis paper presents progress towards a transition modelling capability for use in the numerical solution of the Reynolds-averaged Navier-Stokes equations that provides accurate predictions for transonic flows and is thus suitable for use in the design of wings for aircraft flying at transonic speeds. To this end, compressibility corrections are developed and investigated to extend commonly used empirical correlations to transonic flight conditions while retaining their accuracy at low speeds. A compressibility correction for Tollmien-Schlichting instabilities is developed and applied to a smooth local correlation-based transition model and a stationary crossflow instability compressibility correction is included by adding a new crossflow source term function. Two- and three-dimensional transonic transition test cases demonstrate that the Tollmien-Schlichting compressibility correction produces substantially improved agreement with the experimental transition locations, particularly for higher Reynolds number applications where the effects of flow compressibility are expected to be more significant, such as the NASA CRM-NLF wing-body configuration, while the crossflow compressibility correction prevents an inaccurate, upstream transition front. The compressibility corrections and modifications do not significantly affect the numerical behaviour of the model, which provides an efficient alternative to non-local and higher-fidelity approaches, and can be applied to other transport-equation-based transition models with low-speed empirical correlations without affecting their predictive capability in the incompressible regime.

Numerical Simulation of bio-heat transfer for cryoablation of regularly shaped tumours in liver tissue using multiprobes

The present study focuses upon pre-operative planning strategy for cryosurgical treatment of multiple regularly shaped tumours inside a three dimensional liver tissue. Numerical simulations provide an optimal framework to predict the number of cryo-probes, their placement, operation time and thermal necrosis to the tumour and surrounding healthy tissues. An efficient cryosurgery process requires keeping the tumour cell under lethal temperature which is between −40 °C to −50 °C. The freezing process of undesired tumour tissues involves phase transition from liquid phase to solid phase, the accurate capturing of transition front and size or location of ice balls generated in the process are the important factors of cryosurgery. In this study, fixed domain heat capacity method has been utilized to take into account the latent heat of phase change in bio-heat transfer equation. The ice balls generated with different number of probes haven been analysed. Numerical simulations have been carried out using standard Finite Element Method with COMSOL 5.5 and results obtained are validated with previous studies.

Effect of environmental disturbances on crossflow instability

Wind tunnel experiments on the receptivity of three-dimensional boundary layers were performed in a range of freestream turbulence intensities, Tu, from 0.01%—the lowest level ever achieved in this type of work—up to 0.41%. This work confirms that for Tu=0.01%, and presumably below this level, the transition process is dominated by stationary modes. These are receptive to surface roughness and generate Type-I and Type-II secondary instabilities that eventually cause the transition to turbulence. The saturation amplitude of these stationary waves is highly sensitive to the level of environmental disturbances; the former is here recorded to be the highest in the literature, with the latter being the lowest. Travelling modes are still present; however, their influence on the transition process is marginal. At matched surface roughness levels, when the level of environmental disturbance is enhanced to Tuge 0.33%, the travelling modes acquire more importance, strongly influencing the laminar/turbulent transition process, whilst the initial amplitude and growth of the stationary modes are hindered. For this level of Tu, is the interaction of steady and unsteady disturbances that produces highly amplified waves (Type-III), that quickly lead to nonlinear growth and anticipated turbulence. Finally, a simple rule of thumb is proposed, where the transition front was found to move forward by roughly 10% chord for an increase in one order of magnitude in the Tu levels.Graphical abstract

Spatial buffer-area assisted antiferroelectric-ferroelectric transition in NaNbO3

Lead-free antiferroelectric materials can be used in environmental-friendly energy applications and are receiving tremendous attention owing to their high energy-storage density, which is achieved through a phase transition between antiferroelectric and ferroelectric state. Here, by slowing down the antiferroelectric-ferroelectric transition via controllable beam irradiation in a transmission electron microscope, we demonstrate that such a transition in NaNbO3 undergoes by means of forming a buffer-area (about 40–150 nm wide) at the on-going transition front. The spatial buffer-area approximately maintains its shape during moving and are characterized to be the antiferroelectric phase but exhibits a gradual variation of lattice spacing along the [010] proceeding direction of phase transition. With help of such buffer-area, the overall stress or barrier between the two phases could be considerably reduced compared to the direct change between the two structures. Therefore, our findings provide new insights for understanding the dynamical energy-storage process in antiferroelectric materials.

Reaction-diffusion fronts in funnel-shaped domains

We consider bistable reaction-diffusion equations in funnel-shaped domains of RN made up of straight parts and conical parts with positive opening angles. We study the large time dynamics of entire solutions emanating from a planar front in the straight part of such a domain and moving into the conical part. We show a dichotomy between blocking and spreading, by proving especially some new Liouville type results on stable solutions of semilinear elliptic equations in the whole space RN. We also show that any spreading solution is a transition front having a global mean speed, which is the unique speed of planar fronts, and that it converges at large time in the conical part of the domain to a well-formed front whose position is approximated by expanding spheres. Moreover, we provide sufficient conditions on the size R of the straight part of the domain and on the opening angle α of the conical part, under which the solution emanating from a planar front is blocked or spreads completely in the conical part. We finally show the openness of the set of parameters (R,α) for which the propagation is complete.

Computational study of boundary layer effects on stochastic rotor blade vortex shedding noise

This work entailed the use of lattice-Boltzmann simulations to investigate the effects of different sUAS rotor blade surface treatments, volumetric grid resolutions, and surface roughness heights on broadband noise prediction. Aerodynamic results showed trailing edge flow separation and a mispredicted laminar separation bubble for the laminar-to-turbulent transitional wall-function case and the fully-turbulent wall-function case as well as a transition front for the laminar-to-turbulent transitional wall-function case. It was seen that the current version of PowerFLOW (6-2021) cannot predict turbulent structures associated with turbulent wall-bounded flow unless a triggering mechanism (e.g., boundary layer trip) is employed, possibly explaining aerodynamic force and broadband noise underprediction. Tonal noise results showed thickness (i.e., geometry dependent) noise dominance at out-of-plane observers located above the rotor plane for the fundamental blade passage frequency and at all observer locations for the second harmonic of the blade passage frequency. Extrapolation techniques were also successfully used over broadband noise results to separate regions of underpredicted broadband noise caused by inadequate spatial resolution from other underpredicted regions caused by the misprediction of flow physics using the current computational strategy.

Unsteady disturbances in a swept wing boundary layer due to plasma forcing

This work investigates the response of a transitional boundary layer to spanwise-invariant dielectric barrier discharge plasma actuator (PA) forcing on a 45° swept wing at a chord Reynolds number of 2.17×106. Two important parameters of the PA operation are scrutinized, namely, the forcing frequency and the streamwise location of forcing. An array of passive discrete roughness elements is installed near the leading edge to promote and condition a set of critical stationary crossflow (CF) instability modes. Numerical solutions of the boundary layer equations and linear stability theory are used in combination with the experimental pressure distribution to provide predictions of critical stationary and traveling CF instabilities. The laminar–turbulent transition front is visualized and quantified by means of infrared thermography. Measurements of velocity fields are performed using hotwire anemometry scans at specific chordwise locations. The results demonstrate the inherent introduction of unsteady velocity disturbances by the plasma forcing. It is shown that, depending on actuator frequency and location, these disturbances can evolve into typical CF instabilities. Positive traveling low-frequency type III modes are generally amplified by PA in all tested cases, while the occurrence of negative traveling high-frequency type I secondary modes is favored when PA is operating at high frequency and at relatively downstream locations, with respect to the leading edge.

Instability and Transition of Hypersonic Boundary Layer on a Blunted Delta Wing

The instability and transition of hypersonic boundary layer on a blunted delta wing is performed by numerical simulation of experiments conducted in the quiet wind tunnel. Disturbances are calculated using linear stability theory to obtain the integrated disturbance growth rates and amplification factors from laminar steady base flows. Comparisons are made between the results of the stability analysis and the experimentally measured transition characteristics. The results show that the generalized inflections existing in velocity profiles reveal both the streamwise instability and crossflow instability of the boundary layer on the delta wing. In the parameter spaces of frequency and wave number, there is only one unstable region in which the first-mode and second-mode waves are integrated with each other at angle of attack, whereas the unstable regions of the first-mode and second-mode waves are separately distributed at angle of attack of . The numerically predicted factor distribution is in good agreement with the measured transition front, and the factor is about 3 when boundary-layer transition occurs. The low-frequency unstable waves with about 25–30 kHz play a key role in dominating the boundary-layer transition of the delta wing.

Türbülans Yoğunluğunun Geçiş Yer Tahminine Etkisinin Bölgesel Korelasyon Geçiş Modeli ile İncelenmesi

Aerodynamic parameters are among the most important parameters in flight performance analyses of air/sea vehicles. Free stream turbulence intensity has significant importance on aerodynamic parameters obtained with computational fluid dynamic analyses and wind tunnel tests measurements. In this study, the effect of free stream turbulence intensity on the transition locations of spheroid geometry using with the widely accepted local correlation transition model. The computations are performed at 6.5 x 106 Reynolds number with 0 and 5 degree angle of attacks. The computations at 5 degree angle of attack are compared with an available experimental study. The effect of free stream turbulence intensity is investigated with axial force coefficient, normal force coefficient, and surface friction coefficient distributions. It is seen that axial and normal force coefficients increase with increasing free stream turbulence intensity. The differences in force coefficients obtained with the free stream intensity range used in the study shall create noteworthy effects in flight performance analyses. When surface friction coefficients are investigated, the transition model estimates the transition locations earlier while free stream turbulence intensity increases as expected. However, the transition front geometry is obtained significantly different with respect to the experimental results.

Entire solutions of combustion reaction-diffusion equations in exterior domains

In this paper, we consider entire solutions of combustion reaction-diffusion equations in exterior domains. By constructing appropriate super- and sub-solutions, we show that any given transition front of homogeneous combustion reaction-diffusion equation can emanate an entire solution in exterior domains. In addition, we prove that the entire solution is a complete propagation in the exterior domain with some suitable spatial geometrical conditions.

Estimation of Discretization Uncertainty Using the γ−Reθ Transition Model for Transitional Flows on 6:1 Spheroid

Abstract This paper aims to estimate the surface mesh size related discretization uncertainties using the γ−Reθ transition model combined with the shear stress transport (SST) k–ω turbulence model. For comparison, this work employs an available experimental study performed with a 6:1 prolate spheroid. The grid convergence index (GCI) study is performed for axial force, surface skin friction, and pressure coefficients with three levels of meshes. The transition model estimates the axial force coefficients (CX), approximately half of which are obtained using fully turbulent calculations with higher GCI values. The GCI values around the axial force coefficients for the level-2 mesh are less than 1% based on fully turbulent calculations. However, with the transition model, these values for the same mesh level increase to 10%. While the GCI values of surface pressure coefficients are very small based on both fully turbulent and transition model calculations, these coefficients show differences at the trailing part of the spheroid. Significant differences are also observed in the surface friction coefficients. While the model captures drastic changes in terms of transition in the surface friction coefficients at the suction side of the spheroid, such drastic change is not observed in fully turbulent calculations. On the other hand, there is no sign of any transition phenomenon at the pressure side, contrary to the observations of experimental measurements. The transition model is not able to estimate the transition front geometry correctly. The GCI values of the surface friction coefficients increase dramatically, up to 765% around the transition regions.

Modeling heat and mass transfer during artificial ground freezing considering the influence of water seepage

In the system of artificial ground freezing under water seepage, the dynamic liquid–ice phase equilibrium relationship considering the contribution of water seepage pressure was established based on the classical thermodynamics fundamentals. The freezing characteristic function was represented by combining the Young–Laplace law with van Genuchten model. A thermal-hydraulic model of saturated porous media considering dynamic liquid–ice phase equilibrium was established and validated thorough laboratory tests under various seepage velocities. The temperature field evolution laws of sand, silty and clay grounds were studied based on the thermal-hydraulic model. The temperature field development laws of primary side tunnel of subway cross passages in silty-clay layer were analyzed by applying the thermal-hydraulic model. The relationship between water seepage pressure and the decrease of water freezing point temperature was obtained according to the dynamic liquid–ice phase equilibrium. Compared with the simulation results of considering static liquid–ice phase equilibrium relationship, considering the dynamic liquid–ice phase equilibrium relationship changes the distribution of water-ice phase transition front and has significant influence on the temperature evolution in silty ground. However, that influence in sand ground and clay ground is insignificant because of the relatively lower water seepage pressure in sand ground and poor thermal conductivity in clay ground. Due to the cold energy accumulation bring by water seepage, the freezing rate in the downstream of freezing pipes is faster than that in the upstream of freezing pipes. This effects are exacerbated by lower decrease of water freezing point temperature caused by minor water seepage pressure in the downstream of freezing pipes. Four representative temperature measured points were selected to acquire the relationship between simulation and on-site temperatures. It was observed that the trends of simulation results reflect the tunnel freezing process well.

Response of the Jovian Magnetosphere‐Ionosphere System to the Interplanetary Magnetic Field Discontinuity: A Simulation Study

AbstractThe response of the magnetosphere‐ionosphere (MI) system at Jupiter to the interplanetary magnetic field (IMF) tangential discontinuity is studied via global Magnetohydrodynamics (MHD) simulations. Our results show that the IMF discontinuity, which has a transition front perpendicular to the Sun‐Jupiter line, has strong effects on the Jovian MI system depending on local times and radial distances to Jupiter. The magnetospheric significant response to the discontinuity is delayed by ∼2–6 hr after the encounter of the discontinuity front with the bow shock, while the delay time for the ionospheric field‐aligned currents (FACs) is ∼3–8 hr. The outer magnetosphere, particularly at dayside, tends to respond slightly quicker than the inner regions. At the ionosphere, the response of downward FACs is found to occur earlier than that of the upward ones by a time ranging from less than 1 hr to several hours, which is also more prominent at dayside sectors. Using a new MI mapping method with Jupiter's JRM09 magnetic field model, we reproduce for the first time the ionospheric FACs consistent with the morphology of the observed ultraviolet (UV) main auroral emissions. This indicates a potential application of the simulated upward FACs on the morphology of the UV main auroral ovals.

Rheology of fibrous gels under compression

A number of biological tissues and synthetic gels consist of a fibrous network infused with liquid. There have been a few experimental studies of the rheological properties of such gels under applied compressive strain. Their results suggest that a plot of rheological moduli as a function of applied compressive strain has a long plateau flanked by a steeply increasing curve for large compressive strains and a slowly decreasing curve for small strains. In this paper we explain these trends in rheological properties using a chemo-elastic model characterized by a double-well strain energy function for the underlying fibrous network. The wells correspond to rarefied and densified phases of the fibrous network at low and high strains, respectively. These phases can co-exist across a movable transition front in the gel in order to accommodate overall applied compression. We find that the rheological properties of fibrous gels share similarities with a Kelvin–Voigt visco-elastic solid. The storage modulus has its origins in the elasticity of the fibrous network, while the loss modulus is determined by the dissipation caused by liquid flow through pores. The rheological properties can depend on the number of phase transition fronts present in a compressed sample. Our analysis may explain the dependence of storage and loss moduli of fibrin gels on the loading history. We also point to the need for combining rheological measurements on gels with a microstructural analysis that could shed light on various dissipation mechanisms.

Self-regulated non-reciprocal motions in single-material microstructures.

Living cilia stir, sweep and steer via swirling strokes of complex bending and twisting, paired with distinct reverse arcs1,2. Efforts to mimic such dynamics synthetically rely on multimaterial designs but face limits to programming arbitrary motions or diverse behaviours in one structure3-8. Here we show how diverse, complex, non-reciprocal, stroke-like trajectories emerge in a single-material system through self-regulation. When a micropost composed of photoresponsive liquid crystal elastomer with mesogens aligned oblique to the structure axis is exposed to a static light source, dynamic dances evolve as light initiates a travelling order-to-disorder transition front, transiently turning the structure into a complex evolving bimorph that twists and bends via multilevel opto-chemo-mechanical feedback. As captured by our theoretical model, the travelling front continuously reorients the molecular, geometric and illumination axes relative to each other, yielding pathways composed from series of twisting, bending, photophobic and phototropic motions. Guided by the model, herewe choreograph a wide range of trajectories by tailoring parameters, including illumination angle, light intensity, molecular anisotropy, microstructure geometry, temperature and irradiation intervals and duration. We further show how this opto-chemo-mechanical self-regulation serves as a foundation for creating self-organizing deformation patterns in closely spaced microstructure arrays via light-mediated interpost communication, as well as complex motions of jointed microstructures, with broad implications for autonomous multimodal actuators in areas such as soft robotics7,9,10, biomedical devices11,12 and energy transduction materials13, and for fundamental understanding of self-regulated systems14,15.

Receptivity of crossflow instability to discrete roughness amplitude and location

The effect of discrete roughness elements on the development and breakdown of stationary crossflow instability on a swept wing is explored. Receptivity to various element heights and chordwise locations is explored using a combination of experimental and theoretical tools. Forcing configurations, determined based on linear stability predictions, are manufactured and applied on the wing in a low turbulence facility. Measurements are performed using infrared thermography, quantifying the transition front location, and planar particle image velocimetry, providing a reconstruction of stationary crossflow instabilities and their associated growth. Measurements are corroborated with simulations based on nonlinear parabolised stability equations. Results confirm the efficacy of discrete roughness elements in introducing and conditioning stationary crossflow instabilities. Primary instability amplitudes and resulting laminar-turbulent transition location are found to strongly depend on both roughness amplitude and chordwise location. The Reynolds number based on element height is found to satisfactorily approximate the initial forcing amplitude, revealing the importance of local velocity effects in non-zero-pressure gradient flows. Direct estimation of initial perturbation amplitudes from nonlinear simulations suggests the existence of pertinent flow mechanisms in the element vicinity, active in conditioning the onset of modal instabilities. Dedicated velocimetry planes, elucidate the development of a momentum deficit wake which rapidly decays downstream of the element followed by mild growth, representing the first experimental evidence of transient behaviour in swept wing boundary layers. The outcome of this work identifies a strong scalability of the transition dynamics to roughness amplitude and location, warranting the upscaling of roughness elements to more accessible, measurable and spatially resolved configurations in future experiments.

Amorphous-to-crystal transition in the layer-by-layer growth of bivalve shell prisms

Biomineralization integrates complex physical and chemical processes bio-controlled by the living organisms through ionic concentration regulation and organic molecules production. It allows tuning the structural, optical and mechanical properties of hard tissues during ambient-condition crystallisation, motivating a deeper understanding of the underlying processes. By combining state-of-the-art optical and X-ray microscopy methods, we investigated early-mineralized calcareous units from two bivalve species, Pinctada margaritifera and Pinna nobilis, revealing chemical and crystallographic structural insights. In these calcite units, we observed ring-like structural features correlated with a lack of calcite and an increase of amorphous calcium carbonate and proteins contents. The rings also correspond to a larger crystalline disorder and a larger strain level. Based on these observations, we propose a temporal biomineralization cycle, initiated by the production of an amorphous precursor layer, which further crystallizes with a transition front progressing radially from the unit centre, while the organics are expelled towards the prism edge. Simultaneously, along the shell thickness, the growth occurs following a layer-by-layer mode. These findings open biomimetic perspectives for the design of refined crystalline materials. Statement of significanceCalcareous biominerals are amongst the most present forms of biominerals. They exhibit astonishing structural, optical and mechanical properties while being formed at ambient synthesis conditions from ubiquitous ions, motivating the deep understanding of biomineralization. Here, we unveil the first formation steps involved in the biomineralization cycle of prismatic units of two bivalve species by applying a new multi-modal non-destructive characterization approach, sensitive to chemical and crystalline properties. The observations of structural features in mineralized units of different ages allowed the derivation of a temporal sequence for prism biomineralization, involving an amorphous precursor, a radial crystallisation front and a layer-by-layer sequence. Beyond these chemical and physical findings, the herein introduced multi-modal approach is highly relevant to other biominerals and bio-inspired studies.