Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “velocity model evaluation”

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 703 records · Page 39

Experimental investigation of shock-cell noise reduction for dual-stream nozzles in simulated flight comprehensive data report. Volume 2: Laser velocimeter data, static pressures and shadowgraph photos

Parameters which contribute to supersonic jet shock noise were investigated for the purpose of determining means to reduce such noise generation to acceptable levels. Six dual-stream test nozzles with varying flow passage and plug closure designs were evaluated under simulated flight conditions in an anechoic chamber. All nozzles had combined convergent-divergent or convergent flow passages. Mean velocity and turbulence velocity measurements of 25 selected flow conditions were performed employing a laser Doppler velocimeter. Static pressure measurements were made to define the actual convergence-divergence condition. Test point definition, tabulation of aerodynamic test conditions, velocity histograms, and shadowgraph photographs are presented. Flow visualization through shadowgraph photography can contribute to the development of an analytical prediction model for shock noise from coannular plug nozzles.

Yamamoto, K.↗

On the effects of islands' geometry and size on inducing sea breeze circulation

Model scaling of the characteristics of the sea breeze circulation involved with circular and slab elongated-symmetric midlatitude islands was performed. Larger horizontal and vertical velocities were indicated over the small circular islands as compared to those over corresponding elongated islands. The circulation characteristics of both types of islands become similar when the half-width of the island is about the same as the distance of inland penetration of the sea breeze in a nonisland case. Evaluation of the results through a scale analysis is presented.

Mahrer, Y.↗

Comparison of Land Skin Temperature from a Land Model, Remote Sensing, and In-situ Measurement

Land skin temperature (Ts) is an important parameter in the energy exchange between the land surface and atmosphere. Here hourly Ts from the Community Land Model Version 4.0, MODIS satellite observations, and in-situ observations in 2003 were compared. Compared with the in-situ observations over four semi-arid stations, both MODIS and modeled Ts show negative biases, but MODIS shows an overall better performance. Global distribution of differences between MODIS and modeled Ts shows diurnal, seasonal, and spatial variations. Over sparsely vegetated areas, the model Ts is generally lower than the MODIS observed Ts during the daytime, while the situation is opposite at nighttime. The revision of roughness length for heat and the constraint of minimum friction velocity from Zeng et al. [2012] bring the modeled Ts closer to MODIS during the day, and have little effect on Ts at night. Five factors contributing to the Ts differences between the model and MODIS are identified, including the difficulty in properly accounting for cloud cover information at the appropriate temporal and spatial resolutions, and uncertainties in surface energy balance computation, atmospheric forcing data, surface emissivity, and MODIS Ts data. These findings have implications for the cross-evaluation of modeled and remotely sensed Ts, as well as the data assimilation of Ts observations into Earth system models.

Modeling↗

JT8D-100 turbofan engine, phase 1

The JT8D turbofan engine, widely used in short and medium range transport aircraft, contributes substantially to airport community noise. The jet noise is predominant in the JT8D engine and may be reduced in a modified engine, without loss of thrust, by increasing the airflow to reduce jet velocity. A configuration study evaluated the effects of fan airflow, fan pressure ratio, and bypass ratio on noise, thrust, and fuel comsumption. The cycle selected for the modified engine was based upon an increased diameter, single-stage fan and two additional core engine compressor stages, which replace the existing two-stage fan. Modifications were also made to the low pressure turbine to provide the increased torque required by the larger diameter fan. The resultant JT8D-100 engine models have the following characteristics at take-off thrust, compared to the current JT8D engine: Airflow and bypass ratio are increased, and fan pressure ratio and engine speed are reduced. The resultant engine is also longer, larger in diameter, and heavier than the JT8D base model, but these latter changes are compensated by the increased thrust and decreased fuel comsumption of the modified engine, thus providing the capability for maintaining the performance of the current JT8D-powered aircraft.

Source record↗

Development of a fire test facility for graphite fiber-reinforced composites

A number of modifications were added to the Model 25 facility to make it specifically applicable for composite material screening tests. Most significant was the development of hardware for trapping fibers released during test and isolating them for quantitative measurement. Capability was added for increasing test section velocities and increasing the range of air/fuel ratios available from very rich to very lean. A provision was added for agitation of the test specimen and the combustion gases by a pulsating gas supply. A variety of specimen configurations was evaluated by an extensive series of tests on a graphite-epoxy reference composite which is representative of material currently used in advanced aircraft. A standard test technique and specimen configuration were established which was satisfactory for the reference composite. A test program was then performed on several alternate composite materials to verify the validity of the test method for relative evaluation of a variety of composite materials.

Alexander, J. G.↗

Method for the prediction of the installation aerodynamics of a propfan at subsonic speeds: User manual

This document is the user's guide for the method developed earlier for predicting the slipstream wing interaction at subsonic speeds. The analysis involves a subsonic panel code (HESS code) modified to handle the propeller onset flow. The propfan slipstream effects are superimposed on the normal flow boundary condition and are applied over the surface washed by the slipstream. The effects of the propeller slipstream are to increase the axial induced velocity, tangential velocity, and a total pressure rise in the wake of the propeller. Principles based on blade performance theory, momentum theory, and vortex theory were used to evaluate the slipstream effects. The code can be applied to any arbitrary three dimensional geometry, expressed in the form of HESS input format. The code can handle a propeller alone configuration or a propeller/nacelle/airframe configuration, operating up to high subcritical Mach numbers over a range of angles of attack. Inclusion of a viscous modelling is briefly outlined. Wind tunnel results/theory comparisons are included as examples for the application of the code to a generic supercritical wing/overwing Nacelle with a powered propfan. A sample input/output listing is provided.

Chandrasekaran, B.↗

A Variational Assimilation Method for Satellite and Conventional Data: a Revised Basic Model 2B

A variational objective analysis technique that modifies observations of temperature, height, and wind on the cyclone scale to satisfy the five 'primitive' model forecast equations is presented. This analysis method overcomes all of the problems that hindered previous versions, such as over-determination, time consistency, solution method, and constraint decoupling. A preliminary evaluation of the method shows that it converges rapidly, the divergent part of the wind is strongly coupled in the solution, fields of height and temperature are well-preserved, and derivative quantities such as vorticity and divergence are improved. Problem areas are systematic increases in the horizontal velocity components, and large magnitudes of the local tendencies of the horizontal velocity components. The preliminary evaluation makes note of these problems but detailed evaluations required to determine the origin of these problems await future research.

Achtemeier, Gary L.↗

Hanford Waste Treatment Plant LAB Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation

The Waste Treatment Plant laboratory (LAB) facility stack monitor locations were qualified using scale model stacks to mitigate the risk of identifying that sampling locations do not meet the qualification criteria on the full-scale stack. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale model and its sampling location were geometrically similar to the actual stack and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. Verification tests of the LAB stacks were performed at normal operating conditions. The minimum 1/6 DV value, along with the maximum 6 DV value from the scale model testing determines the range of conditions for which the full-scale stack may be operated and remain in compliance with the stack verification criterion. Based on these DV values, the corresponding stack flow rates for each of the LAB stacks are 625 to 47,237 scfm for LB-C2, 1,704 to 103,131 scfm for LB-S1, and 467 to 18,088 scfm for LB-S2. The remaining criteria for the stack verification to be considered valid involve the flow angle and velocity uniformity results. First, the flow angle at the full-scale stack must be =20°. Second, the velocity uniformity at the full-scale stack must be =20% coefficient of variance (COV). Finally, the velocity uniformity results for the actual and scale model stack tests must agree within 5% COV. These criteria were met through the full-scale stack tests at the LAB facility. Flow angle results were primarily less than 10 degrees, except for the LB-C2 Fan A results, which were an average of 17.6 degrees; all flow angle results were within the =20° criterion. The velocity uniformity results for each test condition averaged between 1.5 and 3.5% COV, which were all within the range of the target %COV values from the scale model tests. Based on these stack verification test results, the three LAB filtered exhaust stack sampling locations meet the qualification criteria provided in the ANSI/HPS N13.1-1999 standard for all fan operating configurations. This includes single-fan as well as dual-fan operations for LB-C2, each of the dual-fan operating conditions for LB-S1, and each single-fan operating condition for LB-S2. Further changes to the system configuration or operating conditions that are outside the bounds described in this report may require additional tests or analyses to determine compliance with the standard.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Hanford Waste Treatment Plant LAB Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation

The Waste Treatment Plant laboratory (LAB) facility stack monitor locations were qualified using scale model stacks to mitigate the risk of identifying that sampling locations do not meet the qualification criteria on the full-scale stack. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale model and its sampling location were geometrically similar to the actual stack, and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. Verification tests of the LAB stacks were performed at normal operating conditions. The minimum 1/6 DV value, along with the maximum 6 DV value from the scale model testing, determines the range of conditions for which the full-scale stack may be operated and remain in compliance with the stack verification criterion. A practical range for the full-scale stack qualification uses the average DV through 6 DV from the scale model tests to compute the corresponding flow rates. Table S1 lists the operating flow rates along with the average and maximum qualified stack flow rates for each of the LAB facility stacks. For each stack, the operating flow is below the maximum qualified stack flow, which means that the scale model test results are acceptable for stack qualification. The remaining criteria for the stack verification to be considered valid involve the flow angle and velocity uniformity results. First, the flow angle at the full-scale stack must be ≤20°. Second, the velocity uniformity at the full-scale stack must be ≤20% coefficient of variance (COV). Finally, the velocity uniformity results for the actual and scale model stack tests must agree within 5% COV. These criteria were met through the full-scale stack tests at the LAB facility. Flow angle results were primarily less than 10°, except for the LB-C2 Fan A results, which were an average of 13.7°; all flow angle results were within the ≤20° criterion. The velocity uniformity results for each test condition averaged between 1.5 and 4.1% COV, which were all within the range of the target percent coefficient of variation values from the scale model tests. Based on these stack verification test results, the three LAB filtered exhaust stack sampling locations meet the qualification criteria provided in the ANSI/HPS N13.1-1999 standard for all fan operating configurations. This includes single-fan as well as dual-fan operations for LB-C2, each of the dual-fan operating conditions for LB-S1, and each single-fan operating condition for LB-S2. Further changes to the system configuration or operating conditions that are outside the bounds described in this report may require additional tests or analyses to determine compliance with the standard.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Influence of bulk turbulence and entrance boundary layer thickness on the curved duct flow field

The influence of bulk turbulence and boundary layer thickness on the secondary flow development in a square, 90 degree turning duct was investigated. A three-dimensional laser velocimetry system was utilized to measure the mean and fluctuating components of velocity at six cross-planes in the duct. The results from this investigation, with entrance boundary layer thickness of 20 percent, were compared with the thin boundary layer results documented in NASA CR-174811. The axial velocity profiles, cross-flow velocities, and turbulence intensities were compared and evaluated with regard to the influence of bulk turbulence intensity and boundary layer thickness, and the influence was significant. The results of this investigation expand the 90 degree curved duct experimental data base to higher turbulence levels and thicker entrance boundary layers. The experimental results provide a challenging benchmark data base for computational fluid dynamics code development and validation. The variation of inlet bulk turbulence intensity provides additional information to aid in turbulence model evaluation.

Crawford, R. A.↗

Pollutant emissions from and within a model gas turbine combustor at elevated pressures and temperatures

Conventional and advanced gas turbine engines are coming under increased scrutiny regarding pollutant emissions. This, in turn, has created a need to obtain in-situ experimental data at practical conditions, as well as exhaust data, and to obtain the data in combustors that reflect modern designs. The in-situ data are needed to (1) assess the effects of design modifications on pollutant formation, and (2) develop a detailed data base on combustor performance for the development and verification of computer modeling. This paper reports on a novel high pressure, high temperature facility designed to acquire such data under controlled conditions and with access (optical and extractive) for in-situ measurements. To evaluate the utility of the facility, a model gas turbine combustor was selected which features practical hardware design, two rows of jets (primary and dilution) with four jets in each row, and advanced wall cooling techniques with laser drilled effusive holes. The dome is equipped with a flat-vaned swirler with vane angles of 60 degrees. Data are obtained at combustor pressures ranging from 2 to 10 atmospheres of pressure, levels of air preheat to 427 C, combustor reference velocities from 10.0 to 20.0 m/s, and an overall equivalence ratio of 0.3. Exit plane and in-situ measurements are presented for HC, O2, CO2, CO, and NO(x). The exit plane emissions of NO(x) correspond to levels reported from practical combustors and the in-situ data demonstrate the utility and potential for detailed flow field measurements.

Drennan, S. A.↗

DHC-6 Twin Otter Tailplane Airfoil Section Testing in the Ohio State University 7x10 Wind Tunnel

Ice contaminated tailplane stall (ICTS) has been found to be responsible for 16 accidents with 139 fatalities over the last three decades, and is suspected to have played a role in other accidents and incidents. The need for fundamental research in this area has been recognized at three international conferences sponsored by the FAA since 1991. In order to conduct such research, a joint NASA/FAA Tailplane Icing Program was formed in 1994: the Ohio State University has played an important role in this effort. The program employs icing tunnel testing, dry wind tunnel testing, flight testing, and analysis using a six-degrees-of-freedom computer code tailored to this problem. A central goal is to quantify the effect of tailplane icing on aircraft stability and control to aid in the analysis of flight test procedures to identify aircraft susceptibility to ICTS. This report contains the results ot testing of a full scale 2D model of a tailplane section of NASA's Icing Research Aircraft, with and without ice shapes, in an Ohio State University 7 x 10 Low Speed wind tunnel in 1994. The results have been integrated into a comprehensive database of aerodynamic coefficients and stability and control derivatives that will permit detailed analysis of flight test results with the analytical computer program. The testing encompassed a full range of angles of attack and elevator deflections, as well as two velocities to evaluate Reynolds number effects. Lift, drag, pitching moment, and hinge moment coefficients were obtained. In addition. instrumentation for use during flight testing was verified to be effective, all components showing acceptable fidelity. Comparison of clean and iced airfoil results show the ice shapes causing a significant decrease in the magnitude of CLmax (from -1.3 to -0.64) and associated stall angle (from -18.6 deg to -8.2 deg). Furthermore, the ice shapes caused an increase in hinge moment coefficient of approximately 0.02, the change being markedly abrupt for one of the ice shapes. A noticeable effect of elevator deflection is that magnitude of the stall angle is decreased for negative (upward) elevator deflections. All these result are consistent with observed tailplane phenomena. and constitute an effective set of data for comprehensive analysis of ICTS

Hiltner, Dale↗

Simulation of ENSO Teleconnections to Precipitation Extremes over the United States in the High-Resolution Version of E3SM

Abstract We evaluate the simulated teleconnection of El Niño–Southern Oscillation (ENSO) to winter season precipitation extremes over the United States in a long (98 years) 1950 control high-resolution version (HR; 25-km nominal atmosphere model horizontal resolution) of the U.S. Department of Energy’s (DOE) Energy Exascale Earth System Model version 1 (E3SMv1). The model bias and spatial pattern of ENSO teleconnections to mean and extreme precipitation in HR overall are similar to the low-resolution model’s (LR; 110 km) historical simulation (four-member ensemble, 1925–59). However, over the southeastern United States (SE-U.S.), HR produces stronger El Niño–associated extremes, reducing LR’s model bias. Both LR and HR produce weaker than observed increase in storm track activity during El Niño events there, but HR improves the ENSO-associated variability of moisture transport over SE-U.S. During El Niño, stronger vertical velocities in HR produce stronger large-scale precipitation, causing larger latent heating of the troposphere that pulls in more moisture from the Gulf of Mexico into the SE-U.S. This positive feedback also contributes to the stronger mean and extreme precipitation response in HR. Over the Pacific Northwest, LR’s bias of stronger than observed La Niña associated extremes is amplified in HR. Both models simulate stronger than observed moisture transport from the Pacific Ocean into the region during La Niña years. The amplified HR bias there is due to stronger orographically driven vertical updrafts that create stronger large-scale precipitation, despite weaker La Niña–induced storm track activity. Significance Statement New high-resolution Earth system models (ESMs) solve mathematical equations of fluid flow at much smaller spatial scales than prevalent ESMs, and thus are prohibitively expensive to compute. However, they can be useful for simulating accurate details of regional climate extremes that are driven by naturally occurring climate oscillations like El Niño–Southern Oscillation (ENSO). Here, we evaluate the simulation of ENSO-driven precipitation extremes over the United States in the high-resolution version of the U.S. Department of Energy’s new Energy Exascale Earth System Model version 1. We find that the high-resolution model improves upon its low-resolution counterpart over the southeastern United States by producing a better transport of moisture into the region from the Gulf of Mexico during El Niño. Over the U.S. Pacific Northwest, the high-resolution model simulates the atmospheric flow in more detail over the complex mountainous terrain. However, it also brings in more moisture from the Pacific Ocean just like the low-resolution model. This causes it to produce precipitation extremes during La Niña years there that are stronger than that observed in the real world.

54 ENVIRONMENTAL SCIENCES↗

EXPRES. I. HD 3651 as an Ideal RV Benchmark

The next generation of exoplanet-hunting spectrographs should deliver up to an order of magnitude improvement in radial velocity (RV) precision over the standard 1 m s{sup −1} state-of-the-art spectrographs. This advance is critical for enabling the detection of Earth-mass planets around Sun-like stars. New calibration techniques such as laser frequency combs and stabilized etalons ensure that the instrumental stability is well characterized. However, additional sources of error include stellar noise, undetected short-period planets, and telluric contamination. To understand and ultimately mitigate error sources, the contributing terms in the error budget must be isolated to the greatest extent possible. Here, we introduce a new high-cadence RV program, the Extreme Precision Spectrograph (EXPRES) 100 Earths Survey, which aims to identify rocky planets around bright, nearby G and K dwarfs. We also present a benchmark case: the 62 day orbit of a Saturn-mass planet orbiting the chromospherically quiet star, HD 3651. The combination of high eccentricity (0.6) and a moderately long orbital period ensures significant dynamical clearing of any inner planets. Our Keplerian model for this planetary orbit has a residual rms of 58 cm s{sup −1} over a ∼6 month time baseline. By eliminating significant contributors to the RV error budget, HD 3651 serves as a standard for evaluating the long-term precision of extreme precision RV programs.

79 ASTRONOMY AND ASTROPHYSICS↗

Mesh-Sequenced Realizations for Evaluation of Subgrid-Scale Models for Turbulent Combustion

This paper develops a new approach for analysis of subgrid closures for turbulent combustion as modeled using direct quadrature (finite-rate chemistry) techniques. The approach, termed multiresolution analysis through mesh-sequenced realizations (MRA-MSR), conducts simultaneous, constrained large-eddy simulations on a set of hierarchically coarsened meshes. Furthermore, the availability of underlying fine-mesh (subgrid) data corresponding to coarse-mesh locations allows a clearer assessment of the effects of unresolved fluctuations on apparent reactivity. A key to MRA-MSR is the correlation of eddy structures at coarser mesh levels, which is facilitated by the transfer of filtered fine-mesh velocity information. A seven-mesh MRA-MSR hierarchy using three resolution levels is applied to one of the Sydney bluff-body stabilized methane–hydrogen flames. Analysis of the simultaneously evolved data at different resolution levels reveals several interesting trends. First, at high Damköhler numbers, there is clear evidence of attenuation of apparent reactivity due to the effects of unresolved fluctuations. Secondly, single-point, single-time filtered density functions of a normalized subgrid Damköhler number show a characteristic beta probability density function (PDF) form and display evidence of scale similarity. Interrogation of the MRA-MSR database also shows that the recently-proposed least-squares minimization (LSM) turbulence–chemistry interaction model can account for the observed diminishment in reactivity at high Damköhler numbers but cannot reduce scatter significantly. A new form of the LSM model, which makes use of the normalized subgrid Damköhler number beta PDF distribution, performs slightly better than the original model, illustrating the potential of MRA-MSR both in assessing existing closure concepts and in developing new ones.

42 ENGINEERING↗

Evaluation of Model Microphysics within Precipitation Bands of Extratropical Cyclones

Recent studies evaluating the bulk microphysical schemes (BMPs) within cloud resolving models (CRMs) have indicated large uncertainties and errors in the amount and size distributions of snow and cloud ice aloft. The snow prediction is sensitive to the snow densities, habits, and degree of riming within the BMPs. Improving these BMPs is a crucial step toward improving both weather forecasting and climate predictions. Several microphysical schemes in the Weather Research and Forecasting (WRF) model down to 1.33‐km grid spacing are evaluated using aircraft, radar, and ground in situ data from the Global Precipitation Mission Coldseason Precipitation Experiment (GCPEx) experiment, as well as a few years (15 winter storms) of surface measurements of riming, crystal habit, snow density, and radar measurements at Stony Brook, NY (SBNY on north shore of Long Island) during the 2009-2012 winter seasons. Surface microphysical measurements at SBNY were taken every 15 to 30 minutes using a stereo microscope and camera, and snow depth and snow density were also recorded. During these storms, a vertically-pointing Ku-band radar was used to observe the vertical evolution of reflectivity and Doppler vertical velocities. A Particle Size and Velocity (PARSIVEL) disdrometer was also used to measure the surface size distribution and fall speeds of snow at SBNY. For the 15 cases at SBNY, the WSM6, Morrison (MORR), Thompson (THOM2), and Stony Brook (SBU-YLIN) BMPs were validated. A non-spherical snow assumption (THOM2 and SBU-YLIN) simulated a more realistic distribution of reflectivity than spherical snow assumptions in the WSM6 and MORR schemes. The MORR, WSM6, and SBU-YLIN schemes are comparable to the observed velocity distribution in light and moderate riming periods. The THOM2 is approx 0.25 m/s too slow with its velocity distribution in these periods. In heavier riming, the vertical Doppler velocities in the WSM6, THOM2, and MORR schemes were approx 0.25 m/s too slow, while the SBU-YLIN was 0.25 to 0.5 m/s too fast. Overall, the BMPs simulate a size distribution close to the observed for D < 4 mm in the dendritic, plates, and mixed habit periods. The model BMPs underestimate the size distribution when large aggregates were observed. For D > 6 mm in the dendrites, side planes, and mixed habit periods, the BMPs are likely not simulating enough aggregation to create a larger size distribution, although the MORR (double moment) scheme seemed to perform best. These SBNY results will be compared with some results from GCPEx for a warm frontal snow band observed at 18 February 2012.

Colle, Brian A.↗

Transient solute drag and strain aging of dislocations

The transient drag force exerted by mobile solutes on a moving dislocation is computed using continuum theory. These mobile solutes form so-called Cottrell atmospheres around dislocations during static and dynamic strain aging. We evaluate the evolution of the drag force exerted by the atmosphere under two velocity time-histories: impulsive acceleration to a chosen velocity and a constant acceleration rate. A particular focus is on the conditions under which the stationary limit assumed by theories of dynamic strain aging is obeyed. According to our results, two conditions—one on the dislocation velocity and one on the acceleration rate—must be satisfied for the stationary limit to hold. Using the Orowan relation and a line tension model, we obtain estimates for the temperature, stress, strain rate, and dislocation density regimes where the stationary limit is valid, and compare these results with experiments for a few material systems.

36 MATERIALS SCIENCE↗

Inflation alternative via the gravitational field of a singularity

Abstract We explore the scenario that the observable Universe emerged from the vicinity of a negative mass ring singularity, and all content of the Universe travels at the same group velocity close to the speed of light on a geodesic trajectory along the axis of rotation of the singularity. In appropriate coordinate parametrization and evaluated on the trajectory, we find that the metric tensor in the vicinity of the trajectory exhibits a conformal scale factor a ( η ) with contraction and subsequent expansion properties that solve the horizon problem. We then introduce a static flow of gravitating radiation along the trajectory (perturbatively with respect to the mass scale of the singularity) to model a homogeneous radiation dominated Universe. Solving the Einstein field equations with a physically motivated ansatz of metric perturbation then reveals that the effective conformal scale factor indeed grows asymptotically with the same power law as expected in a conventional radiation dominated Universe.

Astronomy & Astrophysics↗