# Crispy Filled

## Founded Year

2017## Stage

Shareholder Liquidity | Alive## About Crispy Filled

Crispy Filled is a prepackaged frozen food provider. The company offers Rosti, a line of vegetarian products made from potatoes and cheese, and are sold throughout the nation. Rosti products are gluten-free, low cholesterol, contain no trans fat, and are each 120 calories or less. The company was formerly known as Rosti and changed its name to Crispy Filled in June 2022. It was founded in 2017 and is based in Clackamas, Oregon.

## Missing: Crispy Filled's Product Demo & Case Studies

Promote your product offering to tech buyers.

Reach 1000s of buyers who use CB Insights to identify vendors, demo products, and make purchasing decisions.

## Missing: Crispy Filled's Product & Differentiators

Don’t let your products get skipped. Buyers use our vendor rankings to shortlist companies and drive requests for proposals (RFPs).

## Crispy Filled Patents

Crispy Filled has filed 17 patents.

The 3 most popular patent topics include:

- Graphical projections
- Light machine guns
- Lomanotidae

Application Date | Grant Date | Title | Related Topics | Status |
---|---|---|---|---|

1/16/2018 | 9/10/2019 | Graphical projections, Light machine guns, Cephalopod zootomy, Cuttlefish, Lomanotidae | Grant |

Application Date | 1/16/2018 |
---|---|

Grant Date | 9/10/2019 |

Title | |

Related Topics | Graphical projections, Light machine guns, Cephalopod zootomy, Cuttlefish, Lomanotidae |

Status | Grant |

## Latest Crispy Filled News

Apr 13, 2023

Abstract Non-Newtonian fluids have a viscosity that varies with applied stress. Elastoviscoplastic fluids, the elastic, viscous and plastic properties of which are interconnected in a non-trivial way, belong to this category. We have performed numerical simulations to investigate turbulence in elastoviscoplastic fluids at very high Reynolds-number values, as found in landslides and lava flows, focusing on the effect of plasticity. We find that the range of active scales in the energy spectrum reduces when increasing the fluid plasticity; when plastic effects dominate, a new scaling range emerges between the inertial range and the dissipative scales. An extended self-similarity analysis of the structure functions reveals that intermittency is present and grows with the fluid plasticity. The enhanced intermittency is caused by the non-Newtonian dissipation rate, which also exhibits an intermittent behaviour. These findings have relevance to catastrophic events in natural flows, such as landslides and lava flows, where the enhanced intermittency results in stronger extreme events, which are thus more destructive and difficult to predict. Main Many fluids in nature and industry exhibit a nonlinear relationship between shear stress and shear rate, which is referred to as non-Newtonian behaviour. Several non-Newtonian features can exist, and they are often present simultaneously. In this Article we focus on the so-called elastoviscoplastic (EVP) fluids, which are fluids with elastic, viscous and plastic properties. EVP materials combine solid-like behaviour and fluid-like response depending on the value of the applied stress: they behave like a solid when the applied stress is below a critical value known as the ‘yield stress’, and flow like a liquid otherwise 1 . The elastic nature of these materials is present in their solid as well as liquid states 2 . Such fluids are common in everyday life (examples include toothpaste, jam, cosmetics and mud), and turbulent flows of EVP fluids are found in many industrial processes, including sewage treatment, crude oil transportation, concrete pumping and mud drilling 3 , 4 , 5 , as well as in nature as landslides and lava flows 6 , 7 . A great deal of work has been done in the past to properly characterize the viscoelastic behaviour of a fluid in both laminar and turbulent flows 8 , 9 , 10 , 11 , 12 , 13 , while the effect of plasticity has been studied mainly in low-Reynolds-number laminar conditions 1 , 14 , 15 . Little is known about the plastic behaviour of an EVP fluid in turbulence. Rosti et al. 16 studied a turbulent channel flow of an EVP fluid, finding that the shape of the mean velocity profile controls the regions where the fluid is unyielded, forming plugs around the channel centreline that grow in size as the yield stress increases, similar to what is observed in a laminar condition. However, the presence of the plug region has an opposite effect on drag for laminar and turbulent flow configurations, resulting in drag reduction in the turbulent case and drag increase in the laminar one; the turbulent drag behaviour is due to the tendency of the turbulent flow to relaminarize, overall leading to a strongly nonlinear relation between yield stress and drag coefficient. Simulation results were then employed by Le Clainche et al. 17 , using high-order dynamic mode decomposition, to study the near-wall dynamics, comparing them to those in Newtonian and viscoelastic fluids. Their work revealed that both elasticity and plasticity have similar effects on the near-wall coherent structures, where the flow is characterized by long streaks disturbed for short periods by localized perturbations. A recent experimental study by Mitishita et al. 18 on a turbulent duct flow of Carbopol solution de facto verified the numerical results obtained by Rosti et al. 16 on the effect of plasticity on the mean flow profile and Reynolds stresses. Additionally, they observed an increase in the energy content at large scales and a decrease at small scales, when compared with a Newtonian fluid. Mitishita et al. reported a −7/2 scaling in the energy spectra at high wavenumbers during Carbopol flows compared to −5/3 scaling in the case of water flows. The newly observed scaling was attributed either to the decrease in the inertial effect in the presence of Carbopol solutions, which shrinks the inertial range of scales because the Reynolds numbers are much lower than in water flows, or to the elastic effects that become important at large wavenumbers where the fluid experiences high frequencies. Moreover, the shear-thinning effects that Carbopol solutions exhibit affected the anisotropy and the overall flow behaviour. The elastic and shear-thinning effects are rheological features of Carbopol solutions and cannot be eliminated experimentally. Homogeneous and isotropic turbulent flows have long been a focus of turbulence research for their simple theoretical analysis and the generality of their results. To this end, as has been extensively done in the past for viscoelastic flows, here we study tri-periodic homogeneous flow, where the celebrated K41 theory by Kolmogorov 19 can be directly applied to a classical Newtonian fluid. In this Article we study a homogeneous isotropic turbulent (HIT) flow of an EVP fluid at high Reynolds number, as shown in Fig. 1 . We aim to answer the following fundamental question: how does the Kolmogorov theory change when the fluid is EVP? We will mainly focus on its plastic behaviour and investigate how the yield stress affects the multiscale energy distribution and balance, and how the turbulent energy cascade is altered by the fluid’s plasticity. Our results show profound modifications of the classical picture predicted by the K41 theory for Newtonian fluids, with the emergence of a new scaling range, the dominance of the non-Newtonian flux and dissipation at small and intermediate scales, and enhanced intermittency of the flow. Fig. 1: Instantaneous colourmaps of the turbulent fluid dissipation ϵf in homogeneous isotropic turbulence of an EVP fluid at different Bingham numbers. Yielded regions are shown with a black–red–yellow colourscale, and unyielded regions with black–grey–white. Results To investigate the problem at hand, we performed massive three-dimensional direct numerical simulations (DNS) of HIT where we solve the flow equations fully coupled with the constitutive equation of the EVP fluid, within a tri-periodic domain of size L, using 1,024 grid points per side, as described in more detail in the Methods . The flow is controlled by four main parameters: the Reynolds number ReΛ, the Weissenberg number WiΛ, the viscosity ratio α and the Bingham number BiΛ, all based on the root-mean-square velocity fluctuations u′ and Taylor’s microscale Λ. We use the definitions \({{\rm{Re}}_{\varLambda }\equiv \rho {u}^{{\prime} }\varLambda /{\mu }_{{{{\rm{t}}}}},\,{\rm{Wi}}_{\varLambda }\equiv \lambda {\mu }_{{{{\rm{t}}}}}/\rho {\varLambda }_{0}^{2},\,\alpha ={\mu }_{{{{\rm{n}}}}}/{\mu }_{{{{\rm{t}}}}}}\) and \({\rm{Bi}}_{\varLambda }\equiv {\tau }_{\rm{y}}{\varLambda }_{0}/{\mu }_{{{{\rm{t}}}}}{u}_{0}^{{\prime} }\), where ρ is the fluid density, μt ≡ μf + μn is the total dynamic viscosity with μf being the fluid viscosity and μn the non-Newtonian one, λ is the relaxation time, τy is the yield stress, and subscript 0 denotes quantities from the BiΛ = 0 case. The Reynolds number describes the ratio of inertial to viscous forces, and we limit our analysis to high-Reynolds-number flows, achieving a Taylor microscale Reynolds number ReΛ ≈ 435 for the Newtonian flow, at which statistics of the flow have been found to be universal and exhibiting a proper scale separation, with an extensive inertial range of scales extended to almost two decades of wavenumbers. The Reynolds number explored here is the highest in DNS of HIT of non-Newtonian fluids. The Weissenberg number describes the ratio of elastic to viscous forces, and here we limit the analysis to WiΛ ≪ 1 (that is, WiΛ ≈ 10−3), to ensure that elastic effects are subdominant and all the observed changes are due to plasticity. We also fix a value of α = 0.1 to represent a dilute concentration of polymers, in accordance with previous works on the subject 16 , 20 . Thus, the key control parameter we vary is BiΛ, which describes the ratio of the yield stress to the viscous stress, and thus correlates with the prevalence of unyielded regions. Figure 2 depicts the turbulent kinetic energy spectra of the cases analysed. The BiΛ = 0 case is similar to the Newtonian case shown in Supplementary Fig. 1 , confirming that the effect of elasticity is subdominant and can be ignored. A clear E ~ κ−5/3 range is visible for more than one decade, showing that ReΛ is high enough to achieve scale separation, with the spectra exhibiting an inertial range of scales followed by a dissipative range. As BiΛ increases, the inertial range is limited to the large scales (small wavenumbers κ), with the energy increasing at large scales and decreasing at small scales. A clear deviation from Kolmogorov scaling becomes noticeable for BiΛ > 1, resulting in the emergence of a new apparent scaling of E ~ κ−2.3 that is shown more clearly by plotting compensated energy spectra (as shown in Supplementary Fig. 3 ). The difference in scaling between the experimental work (−7/2) 18 and the current study (−2.3) is mainly due to the higher values of Reynolds and Bingham numbers considered here. The abrupt change in the spectra with BiΛ is consistent with the bulk flow properties (ReΛ and the volume fraction of the unyielded regions Φ) shown in the inset of Fig. 2 : for the cases where BiΛ < 1, ReΛ remains relatively unaltered, with Φ always close to zero, whereas when BiΛ further increases, the microscale Reynolds number ReΛ and the volume Φ of the unyielded regions rapidly increase with a similar trend. Fig. 2: Turbulent kinetic energy spectra of EVP flows with various Bingham numbers. Different BIngham numbers are plotted in colours from dark to light; BiΛ = 0, 0.0025, 0.025, 0.25, 2.5, 12.5 and 25 are plotted in black, purple, dark blue, light blue, dark green, light green and orange, respectively. The expected Kolmogorov scaling for a Newtonian fluid is shown by a grey dashed line, and the grey dash-dotted line represents an apparent new non-Newtonian scaling E ~ κ−2.3 that emerges at large BiΛ. Inset: variation of the mean values of the microscale Reynolds number ReΛ (right axis, plotted as squares) and the volume fraction of the unyielded regions Φ (left axis, plotted as diamonds) as a function of BiΛ. Error bars report the s.d. of ReΛ in time, measured using 103 samples. Plastic effects start to appear for BiΛ ≳ 1, suggesting that Λ is the relevant length scale of the problem. To fully characterize the change in the energy spectra, we study the turbulent kinetic energy balance, which in wavenumber space can be expressed as $${{{{{\mathcal{F}}}}}_{{{{\rm{inj}}}}}(\kappa )+\varPi (\kappa )+{{{\mathcal{D}}}}(\kappa )+{{{\mathcal{N}}}}(\kappa )=\langle {\epsilon }_{{{{\rm{f}}}}}\rangle +\langle {\epsilon }_{{{{\rm{n}}}}}\rangle =\langle {\epsilon }_{{{{\rm{t}}}}}\rangle },$$ (1) where \({{{{\mathcal{F}}}}}_{{{{\rm{inj}}}}}\) is the turbulence production introduced by the external forcing (injected at the largest scale, κL ≡ 2π/L), and \({\varPi ,\,{{{\mathcal{D}}}}}\) and \({{{\mathcal{N}}}}\) are the nonlinear energy flux, the fluid dissipation and the non-Newtonian contribution, respectively. In addition to the classical bulk fluid dissipation rate ϵf, here we have a non-Newtonian dissipation ϵn, which is the rate of removal of turbulent kinetic energy from the flow due to the non-Newtonian extra stress tensor (Supplementary Sections I and II provide a derivation of this equation). Figure 3 shows the turbulent kinetic energy balance for a few representative values of BiΛ. When comparing with Supplementary Fig. 1b , the BiΛ = 0 case closely follows the classical Newtonian turbulent flow, where energy is carried by Π from the large to small scales before being dissipated by the fluid viscosity \({{{\mathcal{D}}}}\). The contribution of the nonlinear convective term Π, which appears as an almost horizontal plateau at relatively large scales, progressively decreases with BiΛ and shrinks towards larger scales, consistent with the reduction of the extension of the inertial range observed in Fig. 2 . The reduced energy flux with BiΛ is also accompanied by a decrease of the fluid dissipation \({{{\mathcal{D}}}}\), which is instead compensated by the increase of non-Newtonian contribution \({{{\mathcal{N}}}}\). At small scales (large κ), the relative importance of the non-Newtonian contribution increases with BiΛ, becoming comparable to the fluid dissipation for BiΛ ≈ 2.5 and eventually becoming the dominant term for BiΛ ≳ 12.5, corresponding to the emergence of the new scaling in the energy spectrum shown in Fig. 2 ; indeed, the non-Newtonian contribution can be interpreted as a combination of a pure energy flux (giving rise to the new scaling region) and a pure dissipative term, as recently suggested by Rosti and others 21 . Regarding the direction of energy flux, Supplementary Fig. 4 shows that we have a direct cascade of energy from large to small scales for all BiΛ (refs. 22 , 23 ). Fig. 3: Scale-by-scale energy balance for different BiΛ. Plotted are the energy flux of the nonlinear convective term Π (dashed lines), solvent dissipation \({{{\mathcal{D}}}}\) (dotted lines) and the non-Newtonian contribution \({{{\mathcal{N}}}}\) (solid lines) for BiΛ = 0 (black), BiΛ = 2.5 (dark green), BiΛ = 12.5 (light green) and BiΛ = 25 (orange). Each term is normalized by the total dissipation rate 〈ϵt〉. \({{{\mathcal{N}}}}\) grows at intermediate and small scales when BiΛ is increased, eventually becoming the dominant contribution. We extend the analysis done in the spectral domain by computing the longitudinal structure functions defined as Sp(r) = 〈(Δu(r))p〉, where p is the order of the structure function and Δu(r) = u(x + r) − u(x) is the difference in the fluid velocity across a length scale r, projected in the direction of r. According to K41, \({S}_{{{{\rm{p}}}}}(r) \sim {\left(\langle {\epsilon }_{{{{\rm{t}}}}}\rangle r\right)}^{p/3}\); however, when the structure functions are displayed as a function of r, as shown in Fig. 4a , they deviate from the K41 prediction as p increases. This phenomenon is thought to be due to the intermittency of the flow, that is, extreme events that are localized in space and time, and thereby break Kolmogorov’s hypothesis of self-similarity in the inertial range 24 . Intermittency results in the scaling exponent of r being a nonlinear concave function of p (instead of p/3) 25 . For the EVP fluid, two scaling regions appear at large BiΛ, with scaling consistent with those from the energy spectra, and with intermittency present in both scaling regions. The role of intermittency in the scaling exponents can be better appreciated when the structure functions are displayed in their extended self-similarity form, obtained by plotting one structure function against another 26 . In Fig. 4b , S4 and S6 are plotted against S2 for all Bingham numbers considered. We note a clear power-law scaling, which deviates from Kolmogorov’s prediction, even for the BiΛ = 0 case (shown in black). The departure from Kolmogorov’s prediction progressively grows as the plasticity of the fluid increases, suggesting that the flow becomes more intermittent due to its plasticity. This becomes more obvious when we plot Sn compensated by the intermittency correction at BiΛ = 0 against S2 (Supplementary Fig. 5 ). Also, intermittency appears to act equally in the two scaling regions present at large BiΛ. Fig. 4: Analysis of structure functions. a, Dependence of the longitudinal velocity structure functions S2 (pluses), S4 (squares) and S6 (hexagons) on separation distance r. Symbol colour denotes BiΛ and is the same as in Fig. 2 . Dashed lines show scalings predicted by K41, and dash-dotted lines show scalings predicted using the new non-Newtonian scaling E ~ κ−2.3. b, The two scalings collapse onto a single line when we plot the structure functions S4 and S6 in extended self-similarity form, that is, against S2. To easily see changes in gradient, we have normalized the structure functions by their values at r ≈ Λ. The dotted line shows a best fit through the data for BiΛ = 0, which deviates from the K41 prediction (dashed line) due to intermittency. Increasing BiΛ further increases this deviation. The inset in b reports the multifractal spectrum of the energy dissipation rate carried out by the fluid ϵf. Intermittency originates from the multifractal nature of the turbulent dissipation rate 24 . For Newtonian fluids, this can be quantified by the multifractal spectrum of the energy dissipation rate, ϵf (refs. 24 , 27 ), which we report in the inset of Fig. 4b . This graph demonstrates that F(α) is nearly identical for all BiΛ cases except for minor variations at small and large values of α. This implies that the fluid dissipation rate is not the cause of the enhanced intermittency observed in the extended self-similarity analysis. In the present flow, the turbulent kinetic energy is dissipated by two different terms, ϵf and ϵn, as seen in Fig. 1 . We thus investigate their respective behaviour by looking at their probability distribution functions (PDFs; Fig. 5 ). We name the non-Newtonian contribution ϵn a ‘dissipation’ because, on average, it removes energy from the flow, giving rise to the positive-skewed distributions in Fig. 5b ; however, unlike the fluid dissipation, it can take positive or negative values at particular locations in space and time. Figure 5a shows that the distribution of ϵf narrows as BiΛ increases 28 ; on the other hand, from Fig. 5b , we see that that the distribution of ϵn substantially broadens as BiΛ increases. Because the non-Newtonian dissipation becomes dominant for the largest BiΛ (as shown in Fig. 3 ), we can thus infer that the extreme values of ϵn are indeed the source of the enhanced intermittency observed from the structure function analysis in Fig. 4 . Fig. 5: PDFs of the fluid dissipation rate ϵf and non-Newtonian dissipation rate ϵn. a,b, PDF of the fluid dissipation rate ϵf (a) and of the non-Newtonian dissipation rate ϵn (b) averaged over time, and normalized by 〈ϵ0〉, the total dissipation of the BiΛ = 0 flow. As BiΛ increases, the PDF of ϵf narrows slightly, whereas the PDF of ϵn widens substantially.

### Crispy Filled Frequently Asked Questions (FAQ)

When was Crispy Filled founded?

Crispy Filled was founded in 2017.

Where is Crispy Filled's headquarters?

Crispy Filled's headquarters is located at 12420 SE Carpenter Drive, Clackamas.

What is Crispy Filled's latest funding round?

Crispy Filled's latest funding round is Shareholder Liquidity.

Who are the investors of Crispy Filled?

Investors of Crispy Filled include L2M Brand Accelerator .

## Discover the right solution for your team

The CB Insights tech market intelligence platform analyzes millions of data points on vendors, products, partnerships, and patents to help your team find their next technology solution.